catalog_generator.py 文件源码

def __init__(self, ap_paramList, ra1,dec1, ra2, dec2, background_density, z):
        '''
        @param ap_paramList[seed]: Seed for random number generator
        @param ra1:  Left right ascension
        @param dec1: Bottom declination
        @param ra2:  Right right ascension
        @param dec2: Top declination
        @param background_density: galaxy background density in galaxies/square degree
        @param z: Redshift of galaxy cluster
        '''

        self.ra1 = ra1
        self.dec1 = dec1
        self.ra2 = ra2
        self.dec2 = dec2
        self.background_density = background_density
        self.z = z

        self.__H0 = 70
        self.__Omega_m = 0.3

        self.__R0 = 2
        self.__Rs = 0.15 / 0.7
        self.__Rcore = 0.1 / 0.7
        self.__norm = 1/quad(nfw,0,self.__R0,args=(1,self.__Rs,self.__Rcore))[0]

        super(CatalogGenerator, self).__init__(ap_paramList)