Python Sigmaの例

コード例 #1

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

 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

コード例 #2

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

 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)

コード例 #3

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

 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)

コード例 #4

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

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

コード例 #5

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

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

コード例 #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_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']

コード例 #8

ファイル: 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']

コード例 #9

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

 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

コード例 #10

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

 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