|
@ -402,7 +402,10 @@ class DensityProfileBEC2dModel(lmfit.Model): |
|
|
self.set_param_hint('sigmax_bec', min=0) |
|
|
self.set_param_hint('sigmax_bec', min=0) |
|
|
self.set_param_hint('sigmay_bec', min=0) |
|
|
self.set_param_hint('sigmay_bec', min=0) |
|
|
self.set_param_hint('sigma_th', min=0) |
|
|
self.set_param_hint('sigma_th', min=0) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
self.set_param_hint('atom_number_bec', expr=f'{self.prefix}amp_bec / 5 * 2 * 3.14159265359 * {self.prefix}sigmax_bec * {self.prefix}sigmay_bec') |
|
|
|
|
|
self.set_param_hint('atom_number_th', expr=f'{self.prefix}amp_th * 2 * 3.14159265359 * 1.20206 / 1.643 * {self.prefix}sigma_th * {self.prefix}sigma_th') |
|
|
|
|
|
self.set_param_hint('condensate_fraction', expr=f'{self.prefix}atom_number_bec / ({self.prefix}atom_number_bec + {self.prefix}atom_number_th)') |
|
|
|
|
|
|
|
|
def guess(self, data, x, y, pre_check=False, post_check=False, **kwargs): |
|
|
def guess(self, data, x, y, pre_check=False, post_check=False, **kwargs): |
|
|
"""Estimate and create initial model parameters for 2d bimodal fit, by doing a 1d bimodal fit along an integrated slice of the image |
|
|
"""Estimate and create initial model parameters for 2d bimodal fit, by doing a 1d bimodal fit along an integrated slice of the image |
|
|