def calculateAvgMassEinsteinRadius(gz, qz): avgMass = 0.247 # avgMass = 0.20358470458734301 dL = cosmo.angular_diameter_distance(gz).to('m') dS = cosmo.angular_diameter_distance(qz).to('m') dLS = cosmo.angular_diameter_distance_z1z2(gz, qz).to('m') thetaE = 4 * const.G * u.Quantity( avgMass, 'solMass').to('kg') * dLS / dL / dS / const.c / const.c thetaEUnit = u.def_unit('theta_E', math.sqrt(thetaE.value) * u.rad) return thetaEUnit
def dLS(self): return cosmo.angular_diameter_distance_z1z2( self.__galaxy.redshift, self.__quasar.redshift).to('lyr')
def _scaleFactor(qR, gR): return (cosmo.angular_diameter_distance_z1z2(gR, qR) / cosmo.angular_diameter_distance(qR) / cosmo.angular_diameter_distance(gR)).to('1/m')