def xpselect(file_path, **kwargs):
"""Application for data subselection.
We want to (loosely) model this on
http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/help/gtselect.txt
"""
assert(file_path.endswith('.fits'))
event_select = xEventSelect(file_path, **kwargs)
outfile = event_select.get('outfile')
if os.path.exists(outfile) and not event_select.get('clobber'):
logger.info('Output file %s already exists.' % outfile)
logger.info('Remove the file or set "clobber = True" to overwite it.')
else:
event_select.select()
return outfile
file_map='casa_cmap.fits'
evt_file_path='casa.fits'
file_selected_path='casa_sel.fits'
file_selected_cube_path='casa_sel_mcube.fits'
outfile=None
for i,a in enumerate(sys.argv):
if '-rad' in a: rad = float(sys.argv[i+1])
elif '-dec' in a: dec = float(sys.argv[i+1])
elif '-ra' in a: ra = float(sys.argv[i+1])
elif '-o' in a: outfile = sys.argv[i+1]
pass
fig=xBinnedMap(file_map).plot(show=False,subplot=(1,2,1))#,figure=fig)
fig.show_circles(ra,dec,rad,lw=1)
rad*=60
evtSelect=xEventSelect(evt_file_path, ra=ra, dec=dec, rad=rad, outfile=file_selected_path,
emax=None, emin=None, mc=False, mcsrcid=[], phasemax=None, phasemin=None, phimax=None, phimin=None, tmax=None, tmin=None)
evtSelect.select()
evtBin=xEventBinningMCUBE(file_selected_path, ebins=1, outfile=file_selected_cube_path, evfile=file_selected_path,
emax=10.0,nypix=256, ebinfile=None, phasebins=50,ebinalg='LIN',
xref=None,phibins=75,nxpix=256,tbins=100,proj='TAN',tstart=None,tbinalg='LIN',algorithm='MCUBE',mc=False,binsz=2.5,yref=None,tbinfile=None,emin=1.0,tstop=None)
evtBin.bin_()
binModulation = xBinnedModulationCube(file_selected_cube_path)
binModulation.plot(show=False,xsubplot=1)
for fit in binModulation.fit_results:
print fit
angle = fit.phase
angle_err = fit.phase_error
visibility = fit.visibility
elif '-dec' in a: dec = float(sys.argv[i + 1])
elif '-ra' in a: ra = float(sys.argv[i + 1])
elif '-o' in a: outfile = sys.argv[i + 1]
pass
fig = xBinnedMap(file_map).plot(show=False, subplot=(1, 2, 1)) #,figure=fig)
fig.show_circles(ra, dec, rad, lw=1)
rad *= 60
evtSelect = xEventSelect(evt_file_path,
ra=ra,
dec=dec,
rad=rad,
outfile=file_selected_path,
emax=None,
emin=None,
mc=False,
mcsrcid=[],
phasemax=None,
phasemin=None,
phimax=None,
phimin=None,
tmax=None,
tmin=None)
evtSelect.select()
evtBin = xEventBinningMCUBE(file_selected_path,
ebins=1,
outfile=file_selected_cube_path,
evfile=file_selected_path,
emax=10.0,
nypix=256,