Commit 42474a6f authored by Yan Zhaojun's avatar Yan Zhaojun
Browse files

debug

parent a5945932

csst_mci_sim/csst_mci_sim.py

+11 −11
Original line number Diff line number Diff line
@@ -1874,7 +1874,7 @@ class MCIsimulator(): 
        nlayccd = 0    

        

        

        self.log.info('load star catlog successfully 0000000')

        

        ############################################    

        

        ra  = self.information['ra_pnt0']
@@ -1893,12 +1893,11 @@ class MCIsimulator(): 
        

        self.log.info('dir_path:')

        self.log.info(self.information['dir_path'])

        self.log.info(time_jd)

        

        

        # EarthShine from straylight

        sl = StrayLight(self.information['dir_path'],jtime=time_jd, sat=np.array([x_sat, y_sat, z_sat]), radec=np.array([(ra*u.degree).value, (dec*u.degree).value]))

        

        self.log.info('load star catlog 3333333333333333333333333')

              

        

        earth_e = sl.caculateEarthShineFilter(filter='r')
@@ -1911,7 +1910,7 @@ class MCIsimulator(): 
        earthshine_wave0, earthshine_flux0 = ill2flux(earth_e+star_e, self.information['dir_path'])

        

        

        self.log.info('load star catlog 444444444')

        

        # sample as ifs wavelength

        wave_mci = np.linspace(2500, 11000, 8501) #np.arange(2500, 11000, 1)

        f1 = interp1d(wave0, zodi0)
@@ -1932,7 +1931,7 @@ class MCIsimulator(): 
        

        ########################################################################################

        

        self.log.info('load star catlog 5555555555555')

        

        #########################################################################################

        

        self.cal_sky_noise()
@@ -1957,7 +1956,7 @@ class MCIsimulator(): 
        

        if self.TianceEffect:      

            

            self.log.info('load star catlog 666666666666666')

            

            ra_list = self.star['ra_gaia'].tolist()

            dec_list = self.star['dec_gaia'].tolist()

            pmra_list = self.star['pmra_gaia'].tolist()
@@ -2145,7 +2144,7 @@ class MCIsimulator(): 
            ##cal_SED_photons(self, flux_arr):

            intscales,PSF_eff_weight=self.cal_SED_photons(star_flux) 

            

            self.log.info('load star catlog 99999999')

            

            

            gx=dict()

            gy=dict()
@@ -2229,7 +2228,7 @@ class MCIsimulator(): 
                cx0,cy0=centroid(stamp_img.array)

                

                

                self.log.info('galsim.PhotonArray.makeFromImage')

               

                

                if self.appFatt :
@@ -2268,7 +2267,6 @@ class MCIsimulator(): 
                photons.addTo(final_image[ch]) 

                

                

                self.log.info('load star catlog 11111111111')

                       

                ###############################

                ###### debug code
@@ -2967,6 +2965,8 @@ class MCIsimulator(): 


        print('total los nummber on CCD is: ',nlayccd )

        

        self.log.info('total los nummber on CCD is = %i' % (nlayccd))

    

        #################################################################  

        tab=Table()            

        tab['gal_ra']=np.asarray(gal_ra)