Python fast_M2dの例

コード例 #1

ファイル: wl_lens_models.py プロジェクト: EiffL/bayesian_hierarchical_wl

 def alpha(self, theta, z):
     return (self.halo_norm/theta*gnfw.fast_M2d(theta/self.rs_ang) + \
            self.mstar/pi*theta/(1./self.Mpc2deg)**2) / self.S_cr(z)

コード例 #2

ファイル: wl_lens_models.py プロジェクト: EiffL/bayesian_hierarchical_wl

 def Sigmabar(self, theta):
     return (self.halo_norm*gnfw.fast_M2d(theta, self.rs_ang, self.gamma)/(theta/self.Mpc2deg)**2/pi + \
             self.mstar/(theta/self.Mpc2deg)**2/pi) /1e12

コード例 #3

ファイル: wl_lens_models.py プロジェクト: EiffL/bayesian_hierarchical_wl

 def gammat(self, theta, z):
     return (self.halo_norm*(gnfw.fast_M2d(theta, self.rs_ang, self.gamma)/pi/(theta/self.Mpc2deg)**2 - gnfw.fast_Sigma(theta/self.Mpc2deg, self.rs, self.gamma)) + \
             self.mstar/(theta/self.Mpc2deg)**2/pi) / self.S_cr(z)

コード例 #4

ファイル: wl_lens_models.py プロジェクト: EiffL/bayesian_hierarchical_wl

 def m(self, theta, z):
     return (self.halo_norm*gnfw.fast_M2d(theta, self.rs_ang, self.gamma) + self.mstar)/ (pi * self.S_cr(z) / self.Mpc2deg**2) 

コード例 #5

ファイル: wl_lens_models.py プロジェクト: EiffL/bayesian_hierarchical_wl

 def get_source_gammat(self):
     self.sources['gammat'] = (self.halo_norm*(gnfw.fast_M2d(self.sources['r'], self.rs_ang, self.gamma)/pi/(self.sources['r']/self.Mpc2deg)**2 - gnfw.fast_Sigma(self.sources['r']/self.Mpc2deg, self.rs, self.gamma)) + \
             self.mstar/(self.sources['r']/self.Mpc2deg)**2/pi) / self.sources['s_cr']

コード例 #6

ファイル: wl_lens_models.py プロジェクト: astrosonnen/bayesian_hierarchical_wl

 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

コード例 #7

ファイル: wl_lens_models.py プロジェクト: astrosonnen/bayesian_hierarchical_wl

 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']

コード例 #8

ファイル: wl_lens_models.py プロジェクト: astrosonnen/bayesian_hierarchical_wl

 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

コード例 #9

ファイル: wl_lens_models.py プロジェクト: astrosonnen/bayesian_hierarchical_wl

 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