def DiffResps(self): """run gtdiffresp""" diffResps['evfile'] = self.eventfile diffResps['scfile'] = self.ft2 diffResps['srcmdl'] = self.xmlfile diffResps['irfs'] = self.irfs diffResps['convert'] = "no" diffResps['evclsmin'] = self.Configuration['analysis']['evclass'] diffResps['clobber'] = self.clobber diffResps.run() print "\ndone"
def DiffResps(self): """run gtdiffresp""" diffResps['evfile']=self.eventfile diffResps['scfile']=self.ft2 diffResps['srcmdl']=self.xmlfile diffResps['irfs']=self.irfs diffResps['convert']="no" diffResps['evclsmin']=self.Configuration['analysis']['evclass'] diffResps['clobber'] = self.clobber diffResps.run() print "\ndone"
def DiffResps(self): """run gtdiffresp""" diffResps['evfile']=self.eventfile diffResps['scfile']=self.ft2 diffResps['srcmdl']=self.xmlfile if self.Configuration['event']['irfs'] != 'CALDB': diffResps['evtype']= self.Configuration['event']['evtype'] diffResps['evclass']=self.Configuration['event']['evclass'] else : diffResps['evtype']= 'INDEF' diffResps['evclass']= 'INDEF' diffResps['irfs']=self.irfs diffResps['convert']="no" diffResps['clobber'] = self.clobber diffResps.run() print "\ndone"
def DiffResps(self): """run gtdiffresp""" diffResps['evfile'] = self.eventfile diffResps['scfile'] = self.ft2 diffResps['srcmdl'] = self.xmlfile if self.Configuration['event']['irfs'] != 'CALDB': diffResps['evtype'] = self.Configuration['event']['evtype'] diffResps['evclass'] = self.Configuration['event']['evclass'] else: diffResps['evtype'] = 'INDEF' diffResps['evclass'] = 'INDEF' diffResps['irfs'] = self.irfs diffResps['convert'] = "no" diffResps['clobber'] = self.clobber diffResps.run() print "\ndone"
def DiffResps(self): """run gtdiffresp""" if (self.clobber == "no" and os.path.isfile(self.diffrspflag)): #print("File exists and clobber is False") return (0) diffResps['evfile'] = self.mktimefile diffResps['scfile'] = self.ft2 diffResps['srcmdl'] = self.xmlfile diffResps['evtype'] = 'INDEF' diffResps['evclass'] = 'INDEF' diffResps['irfs'] = self.irfs diffResps['convert'] = "no" diffResps['clobber'] = self.clobber diffResps.run() with open(self.diffrspflag, "w") as diffrspflag: diffrspflag.write("") print "\ndone"
def DiffResps(self): """run gtdiffresp""" if (self.clobber=="no" and os.path.isfile(self.diffrspflag)): #print("File exists and clobber is False") return(0) diffResps['evfile']=self.mktimefile diffResps['scfile']=self.ft2 diffResps['srcmdl']=self.xmlfile diffResps['evtype']= 'INDEF' diffResps['evclass']= 'INDEF' diffResps['irfs']=self.irfs diffResps['convert']="no" diffResps['clobber'] = self.clobber diffResps.run() with open(self.diffrspflag,"w") as diffrspflag: diffrspflag.write("") print "\ndone"
def diffrsp(times): '''This is the atomic function that actually runs in the seperate threads. It takes a list as input where the first element is tmin, second is tmax, third is the spacecraft file, fourth is the event file, fifth is the source model and sixth is the IRF. It first uses gtselect with wide open cuts to divide up the event file then it runs gtdiffrsp on that event file. The temporary event file is deleted automatically. The function returns the name of the created event file which can be combined with other files and/or deleted later.''' print "Starting calculation on interval {} to {}".format( times[0], times[1]) osfilehandle, outfilename = tempfile.mkstemp(suffix=".fits") filter['rad'] = "INDEF" filter['evclass'] = 0 filter['evclsmin'] = 0 filter['evclsmax'] = 10 filter['infile'] = times[3] filter['outfile'] = outfilename filter['ra'] = "INDEF" filter['dec'] = "INDEF" filter['tmin'] = times[0] filter['tmax'] = times[1] filter['emin'] = 0 filter['emax'] = 400000 filter['zmax'] = 180 filter['convtype'] = -1 filter['chatter'] = 0 filter.run(print_command=False) diffResps['evfile'] = outfilename diffResps['scfile'] = times[2] diffResps['srcmdl'] = times[4] diffResps['irfs'] = times[5] diffResps['chatter'] = 0 diffResps.run(print_command=False) print "Completed calculation of interval {} to {}".format( times[0], times[1]) return outfilename
def diffrsp(times): '''This is the atomic function that actually runs in the seperate threads. It takes a list as input where the first element is tmin, second is tmax, third is the spacecraft file, fourth is the event file, fifth is the source model and sixth is the IRF. It first uses gtselect with wide open cuts to divide up the event file then it runs gtdiffrsp on that event file. The temporary event file is deleted automatically. The function returns the name of the created event file which can be combined with other files and/or deleted later.''' print "Starting calculation on interval {} to {}".format(times[0],times[1]) osfilehandle,outfilename = tempfile.mkstemp(suffix=".fits") filter['rad'] = "INDEF" filter['evclass'] = 0 filter['evclsmin'] = 0 filter['evclsmax'] = 10 filter['infile'] = times[3] filter['outfile'] = outfilename filter['ra'] = "INDEF" filter['dec'] = "INDEF" filter['tmin'] = times[0] filter['tmax'] = times[1] filter['emin'] = 0 filter['emax'] = 400000 filter['zmax'] = 180 filter['convtype'] = -1 filter['chatter'] = 0 filter.run(print_command=False) diffResps['evfile'] = outfilename diffResps['scfile'] = times[2] diffResps['srcmdl'] = times[4] diffResps['irfs'] = times[5] diffResps['chatter'] = 0 diffResps.run(print_command=False) print "Completed calculation of interval {} to {}".format(times[0],times[1]) return outfilename