def get_source_kappa(self):
self.sources['kappa'] = (self.halo_norm*gnfw.fast_Sigma(self.sources['r']/self.Mpc2deg, self.rs, self.gammadm) + \
self.S_bulge*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 (self.halo_norm*gnfw.fast_Sigma(theta/self.Mpc2deg, self.rs, self.gammadm) +\
self.S_bulge*deVaucouleurs.Sigma(theta/self.reff_ang, 1.))/1e12
def get_source_gammat(self):
self.sources['gammat'] = (self.halo_norm*(gnfw.fast_M2d(self.sources['r'], self.rs_ang, self.gammadm)/pi/(self.sources['r']/self.Mpc2deg)**2 - gnfw.fast_Sigma(self.sources['r']/self.Mpc2deg, self.rs, self.gammadm)) + \
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 kappa(self, theta, s_cr):
return (self.halo_norm*gnfw.fast_Sigma(theta/self.Mpc2deg, self.rs, self.gammadm) + \
self.S_bulge*deVaucouleurs.Sigma(theta/self.reff_ang, 1.)) /s_cr