Пример #1
0
 def m(self, theta, z):
     return (4.*self.S_s*self.rs_ang**2*nfw.gfunc(theta/self.rs_ang) + \
             self.mstar/pi/self.Mpc2deg**2) / self.S_cr(z)
Пример #2
0
 def alpha(self, theta, z):
     return (4.*self.S_s*theta/(theta/self.rs_ang)**2*nfw.gfunc(theta/self.rs_ang) + \
            self.mstar/pi*theta/(1./self.Mpc2deg)**2) / self.S_cr(z)
Пример #3
0
 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']
Пример #4
0
 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)
Пример #5
0
 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)
Пример #6
0
 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
Пример #7
0
 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']
Пример #8
0
 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)
Пример #9
0
 def Sigmabar(self, theta):
     return (4.*self.S_s*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 + \
             self.mstar/(theta/self.Mpc2deg)**2/pi) /1e12
Пример #10
0
 def Sigmabar(self, theta):
     return (4.*self.S_s*nfw.gfunc(theta/self.rs_ang)/(theta/self.rs_ang)**2 + \
             self.S_bulge*(sersic.M2d(theta, self.nser, self.reff_ang)/(theta/self.reff_ang)**2/pi)) /1e12
Пример #11
0
 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*sersic.M2d(theta, self.nser, self.reff_ang)) / self.S_cr(z)
Пример #12
0
 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)
Пример #13
0
 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']