def Sigmabar(self, theta):
return (4.*self.S_s*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 + \
self.S_bulge*(deVaucouleurs.fast_M2d(theta/self.reff_ang)/(theta/self.reff_ang)**2/pi)) /1e12
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 alpha(self, theta, z):
return (4.*self.S_s*theta/(theta/self.rs_ang)**2*nfw.gfunc(theta/self.rs_ang) + \
self.S_bulge/pi*theta/(theta/self.reff_ang)**2*deVaucouleurs.fast_M2d(theta/self.reff_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.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 gammat(self, theta, s_cr):
return (self.halo_norm*(gnfw.fast_M2d(theta, self.rs_ang, self.gammadm)/pi/(theta/self.Mpc2deg)**2 - gnfw.fast_Sigma(theta/self.Mpc2deg, self.rs, self.gammadm)) + \
self.S_bulge*(deVaucouleurs.fast_M2d(theta/self.reff_ang)/(theta/self.reff_ang)**2/pi - \
deVaucouleurs.Sigma(theta/self.reff_ang, 1.))) / s_cr
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 Sigmabar(self, theta):
return (self.halo_norm*gnfw.fast_M2d(theta, self.rs_ang, self.gammadm)/(theta/self.Mpc2deg)**2/pi + \
self.S_bulge*(deVaucouleurs.fast_M2d(theta/self.reff_ang)/(theta/self.reff_ang)**2/pi)) /1e12
def alpha(self, theta, s_cr):
return (self.halo_norm/theta*gnfw.fast_M2d(theta/self.rs_ang) + \
self.S_bulge/pi*theta/(theta/self.reff_ang)**2*deVaucouleurs.fast_M2d(theta/self.reff_ang)) \
/ s_cr