def gammat(self, theta, z):
return (2.*self.S_s*(2.*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 - nfw.Ffunc(theta/self.rs_ang)) + \
self.mstar/(theta/self.Mpc2deg)**2/pi) / self.S_cr(z)
def kappa(self, theta, z):
return 2.*self.S_s*nfw.Ffunc(theta/self.rs_ang) / self.S_cr(z)
def get_source_gammat(self):
self.sources['gammat'] = (2.*self.S_s*(2.*nfw.gfunc(self.sources['r']/self.rs_ang)/(self.sources['r']/self.rs_ang)**2 + \
- nfw.Ffunc(self.sources['r']/self.rs_ang)) + \
self.mstar/(self.sources['r']/self.Mpc2deg)**2/pi) / self.sources['s_cr']
def get_source_kappa(self):
self.sources['kappa'] = 2.*self.S_s*nfw.Ffunc(self.sources['r']/self.rs_ang) / self.sources['s_cr']
def kappa(self, theta, z):
return (2.*self.S_s*nfw.Ffunc(theta/self.rs_ang) + self.S_bulge*deVaucouleurs.Sigma(theta/self.reff_ang, 1.)) /self.S_cr(z)
def Sigma(self, theta): # surface mass density in M_Sun/pc^2
return (2.*self.S_s*nfw.Ffunc(theta/self.rs_ang) + self.S_bulge*deVaucouleurs.Sigma(theta/self.reff_ang, 1.))/1e12
def get_source_kappa(self):
self.sources['kappa'] = (2.*self.S_s*nfw.Ffunc(self.sources['r']/self.rs_ang) + \
self.S_bulge*deVaucouleurs.Sigma(self.sources['r']/self.reff_ang, 1.)) / self.sources['s_cr']
def gammat(self, theta, z):
return (2.*self.S_s*(2.*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 - nfw.Ffunc(theta/self.rs_ang)) + \
self.S_bulge*(deVaucouleurs.fast_M2d(theta/self.reff_ang)/(theta/self.reff_ang)**2/pi - \
deVaucouleurs.Sigma(theta/self.reff_ang, 1.))) / self.S_cr(z)
def Sigma(self, theta): # surface mass density in M_Sun/pc^2
return (2.*self.S_s*nfw.Ffunc(theta/self.rs_ang))/1e12
def get_source_gammat(self):
self.sources['gammat'] = (2.*self.S_s*(2.*nfw.gfunc(self.sources['r']/self.rs_ang)/(self.sources['r']/self.rs_ang)**2 + \
- nfw.Ffunc(self.sources['r']/self.rs_ang)) + \
self.S_bulge*(deVaucouleurs.fast_M2d(self.sources['r']/self.reff_ang)/(self.sources['r']/self.reff_ang)**2/pi - \
deVaucouleurs.Sigma(self.sources['r']/self.reff_ang, 1.))) / self.sources['s_cr']
def Sigma(self, theta): # surface mass density in M_Sun/pc^2
return (2. * self.S_s * nfw.Ffunc(theta / self.rs_ang) + self.S_bulge *
sersic.Sigma(theta / self.reff_ang, self.nser, 1.)) / 1e12
def kappa(self, theta, z):
return (
2. * self.S_s * nfw.Ffunc(theta / self.rs_ang) + self.S_bulge *
sersic.Sigma(theta / self.reff_ang, self.nser, 1.)) / self.S_cr(z)
def gammat(self, theta, z):
return (2.*self.S_s*(2.*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 - nfw.Ffunc(theta/self.rs_ang)) + \
self.S_bulge*(sersic.M2d(theta, self.nser, self.reff_ang)/(theta/self.reff_ang)**2/pi - \
sersic.Sigma(theta/self.reff_ang, self.nser, 1.))) / self.S_cr(z)
def get_source_kappa(self):
self.sources['kappa'] = (2.*self.S_s*nfw.Ffunc(self.sources['r']/self.rs_ang) + \
self.S_bulge*sersic.Sigma(self.sources['r']/self.reff_ang, self.nser, 1.)) / self.sources['s_cr']
def get_source_gammat(self):
self.sources['gammat'] = (2.*self.S_s*(2.*nfw.gfunc(self.sources['r']/self.rs_ang)/(self.sources['r']/self.rs_ang)**2 + \
- nfw.Ffunc(self.sources['r']/self.rs_ang)) + \
self.S_bulge*(sersic.M2d(self.sources['r'], self.nser, self.reff_ang)/(self.sources['r']/self.reff_ang)**2/pi - \
sersic.Sigma(self.sources['r']/self.reff_ang, self.nser, 1.))) / self.sources['s_cr']
