def runAnalysisStep(self,bin=0,tmin=0,tmax=0,delete=True):
basename = self.commonConf['base'] + "_bin" + str(bin)
os.symlink(self.commonConf['base'] + ".list", basename+".list")
os.symlink(self.commonConf['base'] + "_SC.fits", basename+"_SC.fits")
analysisConfig = {"ra" : self.analysisConf["ra"],
"dec" : self.analysisConf["dec"],
"rad" : self.analysisConf["rad"],
"tmin" : tmin,
"tmax" : tmax,
"emin" : self.analysisConf["emin"],
"emax" : self.analysisConf["emax"],
"zmax" : self.analysisConf["zmax"],
"binsize" : self.analysisConf["binsize"],
"convtype" : self.analysisConf["convtype"]}
commonConfig = {"base" : basename,
"eventclass" : self.commonConf["eventclass"],
"binned" : False,
"irfs" : self.commonConf["irfs"],
"verbosity" : self.commonConf["verbosity"],
"multicore" : 0}
qA_bin = qA.quickAnalysis(basename,False,analysisConfig,commonConfig)
qA_bin.analysisConf['tmin'] = tmin
qA_bin.analysisConf['tmax'] = tmax
qA_bin.runSelect(True)
qA_bin.runGTI(True)
qA_bin.runLTCube(True)
qA_bin.runExpMap(True)
def runAnalysisStepMP(bininfo):
bin = bininfo[0]
tmin = bininfo[1]
commonConf = bininfo[2]
analysisConf = bininfo[3]
curveConf = bininfo[4]
tmax = tmin + float(curveConf['tstep'])
print bin,tmin,tmax
dir = "quickCurve_"+ str(curveConf['tstep']) + "_bin" + str(bin)
if not os.path.isdir(dir):
os.mkdir(dir)
analysisConfig = {"ra" : analysisConf["ra"],
"dec" : analysisConf["dec"],
"rad" : analysisConf["rad"],
"tmin" : tmin,
"tmax" : tmax,
"emin" : analysisConf["emin"],
"emax" : analysisConf["emax"],
"zmax" : analysisConf["zmax"],
"binsize" : analysisConf["binsize"],
"convtype" : analysisConf["convtype"],
"filter" : analysisConf["filter"],
"roicut" : analysisConf["roicut"],
"ltzmax" : analysisConf["ltzmax"]}
commonConfig = {"base" : commonConf["base"],
"eventclass" : commonConf["eventclass"],
"binned" : False,
"irfs" : commonConf["irfs"],
"verbosity" : commonConf["verbosity"],
"multicore" : 0}
qA_bin = qA.quickAnalysis(commonConf['base'],False,analysisConfig,commonConfig)
qA_bin.analysisConf['tmin'] = tmin
qA_bin.analysisConf['tmax'] = tmax
qA_bin.runSelect(True, False,
outfile = dir + "/" + commonConf['base'] + "_filtered.fits")
qA_bin.runGTI(True,
evfile=dir + "/" + commonConf['base']+'_filtered.fits',
outfile=dir + "/" + commonConf['base']+'_filtered_gti.fits')
qA_bin.runLTCube(True,
evfile = dir + "/" + commonConf['base']+'_filtered_gti.fits',
outfile = dir + "/" + commonConf['base']+'_ltcube.fits')
qA_bin.runExpMap(True,
evfile = dir + "/" + commonConf['base']+'_filtered_gti.fits',
expcube = dir + "/" + commonConf['base']+'_ltcube.fits',
outfile = dir + "/" + commonConf['base']+'_expMap.fits')
