From 90aa4e17b647a123c6b1677a2287c2be6a1cdf7a Mon Sep 17 00:00:00 2001 From: "zhaojunyan@shao.ac.cn" Date: Sat, 13 Apr 2024 22:17:30 +0800 Subject: [PATCH] test --- csst_mci_sim/csst_mci_sim.py | 174 ++++++++++++++++++----------------- 1 file changed, 89 insertions(+), 85 deletions(-) diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py index b13f015..5cd11b4 100644 --- a/csst_mci_sim/csst_mci_sim.py +++ b/csst_mci_sim/csst_mci_sim.py @@ -1252,78 +1252,7 @@ class MCIsimulator(): """ self.configfile = configfile self.section = 'TEST' #####simulation section; - #load instrument model, these values are also stored in the FITS header - self.information = MCIinstrumentModel.MCIinformation() - self.information['G_filters']=["F275W", "F280N","NUV", "WU", "CBU", "F343N", "u", "F373N", "F395N", "F336W"] - self.information['R_filters']=["F487N", "F502N", "CBV", "r", "F656N", "F658N", "F467M", "F555W", "F606W", "F673N"] - self.information['I_filters']=["z", "y", "F815N", "CBI", "F925N", "F960M", "F968N", "F845M" ,"F850LP" ,"F814W"] - - #### load telescope efficiency data - print(self.information['dir_path']) - - self.tel_eff=np.load(self.information['dir_path']+'MCI_inputData/tel_eff/tel_eff.npy',allow_pickle=True).item() - - #### load MCI filter data - self.filterP=np.load(self.information['dir_path']+'MCI_inputData/MCI_filters/mci_filterPWTC.npy',allow_pickle=True).item() - - ########################################################################### - - now=datetime.utcnow() - #data_time=now.strftime("%Y-%m-%d-%H-%M-%S") - result_day=now.strftime("%Y-%m-%d") - #self.result_path='../MCI_simData_'+result_day - - # if os.getcwd() =='/home/zjyan/work/MCI/simulator': - # self.result_path='/data/mci_simData/'+result_day ### CSST1 - # else: - # self.result_path='/home/yan/MCI/simData/'+result_day ### my PC - if self.information['dir_path'] == '/nfsdata/share/simulation-unittest/mci_sim/': - self.result_path = self.information['dir_path'] + \ - 'mci_sim_result/'+result_day - else: - - home_path = os.environ['HOME'] - - if home_path == '/home/yan': - - self.result_path = '../MCI_simData/'+self.source+ss+result_day - else: - self.result_path = '/data/mcisimdata/'+result_day - - - - if os.path.isdir(self.result_path)==False: - os.mkdir(self.result_path) - os.mkdir(self.result_path+'/cali_Data') - os.mkdir(self.result_path+'/log_Data') - os.mkdir(self.result_path+'/ori_Cali') - os.mkdir(self.result_path+'/ori_Sky') - os.mkdir(self.result_path+'/sky_Data') - os.mkdir(self.result_path+'/PSF_Data') - - #update settings with defaults - self.information.update(dict(quadrant=int(0), - ccdx=int(0), - ccdy=int(0), - ccdxgap=1.643, - ccdygap=8.116, - xsize=2000, - ysize=2000, - fullwellcapacity=90000, - pixel_size=0.05, - dark=0.001, - exptime=300.0, - readouttime=4., - rdose=8.0e9, - ghostCutoff=22.0, - ghostRatio=5.e-5, - coveringfraction=0.1, #CR: - cosmicraylengths =FOLDER+'MCI_inputData/data/cdf_cr_length.dat', - cosmicraydistance=FOLDER+'MCI_inputData/data/cdf_cr_total.dat', - parallelTrapfile =FOLDER+'MCI_inputData/data/cdm_euclid_parallel.dat', - serialTrapfile =FOLDER+'MCI_inputData/data/cdm_euclid_serial.dat', - mode='same')) ############################################################################### def readConfigs(self,simnumber,source): @@ -1339,21 +1268,7 @@ class MCIsimulator(): #self.config.readfp(open(self.configfile)) self.config.read_file(open(self.configfile)) - now=datetime.utcnow() - - data_time=now.strftime("%Y-%m-%d-%H-%M-%S") - - print('simnumber', simnumber) - self.log = lg.setUpLogger(self.result_path+'/log_Data/MCIsim_'+source+'_'+data_time+'_Num_'+str(simnumber)+'.log') - - #print('logger.handlers=', self.log.handlers) - - self.log.info('-------STARTING A NEW SIMULATION------------') - - self.log.info('The exposure order is %i ' % simnumber ) - - #self.log.info(self.information) ################################################################################################3 @@ -1549,6 +1464,95 @@ class MCIsimulator(): self._createEmpty() + now=datetime.utcnow() + + data_time=now.strftime("%Y-%m-%d-%H-%M-%S") + + print('simnumber', simnumber) + + self.log = lg.setUpLogger(self.result_path+'/log_Data/MCIsim_'+source+'_'+data_time+'_Num_'+str(simnumber)+'.log') + + #print('logger.handlers=', self.log.handlers) + + self.log.info('-------STARTING A NEW SIMULATION------------') + + self.log.info('The exposure order is %i ' % simnumber ) + + #load instrument model, these values are also stored in the FITS header + self.information = MCIinstrumentModel.MCIinformation() + self.information['G_filters']=["F275W", "F280N","NUV", "WU", "CBU", "F343N", "u", "F373N", "F395N", "F336W"] + self.information['R_filters']=["F487N", "F502N", "CBV", "r", "F656N", "F658N", "F467M", "F555W", "F606W", "F673N"] + self.information['I_filters']=["z", "y", "F815N", "CBI", "F925N", "F960M", "F968N", "F845M" ,"F850LP" ,"F814W"] + + #### load telescope efficiency data + print(self.information['dir_path']) + + self.tel_eff=np.load(self.information['dir_path']+'MCI_inputData/tel_eff/tel_eff.npy',allow_pickle=True).item() + + #### load MCI filter data + self.filterP=np.load(self.information['dir_path']+'MCI_inputData/MCI_filters/mci_filterPWTC.npy',allow_pickle=True).item() + + ########################################################################### + + now=datetime.utcnow() + #data_time=now.strftime("%Y-%m-%d-%H-%M-%S") + result_day=now.strftime("%Y-%m-%d") + #self.result_path='../MCI_simData_'+result_day + + # if os.getcwd() =='/home/zjyan/work/MCI/simulator': + # self.result_path='/data/mci_simData/'+result_day ### CSST1 + # else: + # self.result_path='/home/yan/MCI/simData/'+result_day ### my PC + + if self.information['dir_path'] == '/nfsdata/share/simulation-unittest/mci_sim/': + self.result_path = self.information['dir_path'] + \ + 'mci_sim_result/'+result_day + else: + + home_path = os.environ['HOME'] + + if home_path == '/home/yan': + + self.result_path = '../MCI_simData/'+self.source+ss+result_day + else: + self.result_path = '/data/mcisimdata/'+result_day + + + + if os.path.isdir(self.result_path)==False: + os.mkdir(self.result_path) + os.mkdir(self.result_path+'/cali_Data') + os.mkdir(self.result_path+'/log_Data') + os.mkdir(self.result_path+'/ori_Cali') + os.mkdir(self.result_path+'/ori_Sky') + os.mkdir(self.result_path+'/sky_Data') + os.mkdir(self.result_path+'/PSF_Data') + + #update settings with defaults + self.information.update(dict(quadrant=int(0), + ccdx=int(0), + ccdy=int(0), + ccdxgap=1.643, + ccdygap=8.116, + xsize=2000, + ysize=2000, + fullwellcapacity=90000, + pixel_size=0.05, + dark=0.001, + exptime=300.0, + readouttime=4., + rdose=8.0e9, + ghostCutoff=22.0, + ghostRatio=5.e-5, + coveringfraction=0.1, #CR: + cosmicraylengths =FOLDER+'MCI_inputData/data/cdf_cr_length.dat', + cosmicraydistance=FOLDER+'MCI_inputData/data/cdf_cr_total.dat', + parallelTrapfile =FOLDER+'MCI_inputData/data/cdm_euclid_parallel.dat', + serialTrapfile =FOLDER+'MCI_inputData/data/cdm_euclid_serial.dat', + mode='same')) + + + self.log.info('Read in the configuration files and created an empty array') return -- GitLab