Example #1
0
 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"
Example #2
0
 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"
Example #3
0
    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"
Example #4
0
    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"
Example #5
0
    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"
Example #6
0
    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"
Example #7
0
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
Example #8
0
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