コード例 #1
0
ファイル: gtfunction.py プロジェクト: zblz/enrico
 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"
コード例 #2
0
ファイル: gtfunction.py プロジェクト: jlenain/enrico
 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"
コード例 #3
0
ファイル: gtfunction.py プロジェクト: dcfidalgo/enrico
    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"
コード例 #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"
コード例 #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"
コード例 #6
0
ファイル: gtfunction.py プロジェクト: mireianievas/enrico
    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"
コード例 #7
0
ファイル: gtdiffrsp_mp.py プロジェクト: rsnemmen/gtapps_mp
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
コード例 #8
0
ファイル: gtdiffrsp_mp.py プロジェクト: giacomov/gtburst
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