def Runsrcprob(config):
config['space']['rad'] = config['srcprob']['rad']
Obs = Observation(config['out'], config, tag="srcprob")
_log("Summary")
Obs.printSum()
if config['srcprob']['FitsGeneration'] =='yes':
_log("Select data")
Obs.FirstCut()
if config['analysis']['ComputeDiffrsp'] == 'yes':
_log("Compute diffuse responses")
Obs.DiffResps()
_log("Run gtsrcprob")
Obs.SrcProb()
probfile=fits.open(Obs.Probfile)
srclist = open(config['srcprob']['srclist'],"r").readlines()
for src in srclist:
proba = probfile[1].data.field(string.split(src)[0])
energy = probfile[1].data.field("ENERGY")
ra = probfile[1].data.field("RA")
dec = probfile[1].data.field("DEC")
time = probfile[1].data.field("TIME")
indices = energy.argsort()
mes = Loggin.Message()
mes.info( "Results sorted by decreasing energy")
Print(indices,config,ra,dec,proba,energy,time)
print
mes.info( "Results sorted by decreasing probability")
indices = proba.argsort()
Print(indices,config,ra,dec,proba,energy,time)
def Analysis(folder, config, tag="", convtyp='-1', verbose = 1):
""" run an analysis"""
Obs = Observation(folder, config, convtyp, tag=tag)
if verbose:
utils._log('SUMMARY: ' + tag)
Obs.printSum()
FitRunner = FitMaker(Obs, config)##Class
if config['Spectrum']['FitsGeneration'] == 'yes':
FitRunner.GenerateFits() #Generates fits files
return FitRunner
def Analysis(folder, config, tag="", convtyp='-1', verbose=1):
""" run an analysis"""
Obs = Observation(folder, config, convtyp, tag=tag)
if verbose:
utils._log('SUMMARY: ' + tag)
Obs.printSum()
FitRunner = FitMaker(Obs, config) ##Class
if config['Spectrum']['FitsGeneration'] == 'yes':
FitRunner.GenerateFits() #Generates fits files
return FitRunner
def Analysis(folder, config, configgeneric=None, tag="", convtyp='-1', verbose = 1):
mes = Loggin.Message()
""" run an analysis"""
Obs = Observation(folder, config, tag=tag)
if verbose:
utils._log('SUMMARY: ' + tag)
Obs.printSum()
FitRunner = FitMaker(Obs, config)##Class
if config['Spectrum']['FitsGeneration'] == 'yes':
FitRunner.FirstSelection(configgeneric) #Generates fits files for the coarse selection
FitRunner.GenerateFits() #Generates fits files for the rest of the products
return FitRunner
def Analysis(folder, config, configgeneric=None, tag="", convtyp='-1', verbose = 1):
mes = Loggin.Message()
""" run an analysis"""
Obs = Observation(folder, config, tag=tag)
if verbose:
utils._log('SUMMARY: ' + tag)
Obs.printSum()
FitRunner = FitMaker(Obs, config)##Class
if config['Spectrum']['FitsGeneration'] == 'yes':
FitRunner.FirstSelection(configgeneric) #Generates fits files for the coarse selection
FitRunner.GenerateFits() #Generates fits files for the rest of the products
return FitRunner
def Analysis(folder, config, configgeneric=None, tag="", convtyp='-1', verbose = 1):
mes = Loggin.Message()
""" run an analysis"""
# If there are no xml files, create it and print a warning <--- This should be here?
#if len(glob.glob(config['file']['xml'].replace('.xml','*.xml')))==0:
if len(glob.glob(config['file']['xml']))==0: #.replace('.xml','*.xml')))==0:
mes.warning("Xml not found, creating one for the given config %s" %config['file']['xml'])
XmlMaker(config)
Obs = Observation(folder, config, tag=tag)
if verbose:
utils._log('SUMMARY: ' + tag)
Obs.printSum()
FitRunner = FitMaker(Obs, config)##Class
if config['Spectrum']['FitsGeneration'] == 'yes':
FitRunner.FirstSelection(configgeneric) #Generates fits files for the coarse selection
FitRunner.GenerateFits() #Generates fits files for the rest of the products
return FitRunner
def FindSrc(infile):
config = get_config(infile)
folder = config['out']
Obs = Observation(folder, config)
utils._log('SUMMARY: ')
Obs.printSum()
FitRunner = FitMaker(Obs, config)
if config["findsrc"]["FitsGeneration"]== "yes":
config['analysis']['likelihood'] = 'unbinned'
FitRunner.GenerateFits()
FitRunner._log('gtfindsrc', 'Optimize source position')
os.system("rm "+utils._dump_findsrcout(config))
Obs.FindSource()
try:
update_reg(config)
except:
pass
