automaticLightCurve.py 文件源码

def mergeGTIfiles(self):
        """
        Merge multiple GTI files when mergelongterm is True.
        Use gtselect.
        Assume the current workDir is longTerm/merged.
        """

        # Create list of GTI files
        if not self.daily:
            listname = self.workDir + '/' + self.src + '_gti.list'
        else:
            listname = self.workDir + '/' + self.src + '_daily_gti.list'
        filelist = open(listname, 'w')
        list = []
        if not self.daily:
            for file in glob.glob(self.workDir + '/../20????/' + self.src + '_gti.fits'):
                list.append(file)
        else:
            for file in glob.glob(self.workDir + '/../20????/' + self.src + '_daily_gti.fits'):
                list.append(file)
        # Sort the list of GTI files
        list = sorted(list)
        for item in list:
            filelist.write(item + '\n')
        filelist.close()

        fermi.filter['infile'] = '@' + listname
        if not self.daily:
            outfile = self.workDir + '/' + str(self.src) + '_gti.fits'
        else:
            outfile = self.workDir + '/' + str(self.src) + '_daily_gti.fits'
        fermi.filter['outfile'] = outfile

        # If outfile already exists, we re-create it
        if os.path.isfile(outfile):
            os.remove(outfile)

        fermi.filter['ra'] = self.ra
        fermi.filter['dec'] = self.dec
        fermi.filter['rad'] = self.roi
        fermi.filter['emin'] = self.emin
        fermi.filter['emax'] = self.emax
        fermi.filter['tmin'] = self.tstart
        fermi.filter['tmax'] = self.tstop
        fermi.filter['zmax'] = self.zmax
        fermi.filter['evclass'] = 128
        logging.info('Running gtselect')
        fermi.filter.run()