def sigma_crit(self, z, frame='comoving'):
assert frame in ('comoving', 'physical','proper')
A = self._c**2 / (4*np.pi*self._G)
if frame == 'comoving':
A = A / (1+z)**2
return A / (self.angular_diameter_distance(z)*self.beta(z_s=z_s))
so:
[x] self.beta is not defined
[x] the (1+z)^2 in comoving frame should be the lens redshift
[x] z_s (when calling self.beta) does not exist
[x] If we want self.beta to replicate BaseLensing.beta then the Da multiplying it should be the lens Da, not source
currently:
so:
self.betais not defined(1+z)^2in comoving frame should be the lens redshiftz_s(when callingself.beta) does not existself.betato replicateBaseLensing.betathen the Da multiplying it should be the lens Da, not source