iLat_3C454.3: Examples

01:

the following Ft1 files are availabe
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243760583_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243766312_v002.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243772042_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243777771_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243781394_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243787695_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243793786_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243799818_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243805812_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243811794_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243817770_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243823740_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243829529_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft1/gll_ph_r0243834802_v000.fit
Requested tstart, tstop: 243756000.0 243842400.0
GTI range: 243760585.99 243840783.084
the following Ft2 files are availabe
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243760583_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243766312_v002.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243772042_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243777771_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243781394_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243787695_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243793786_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243799818_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243805812_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243811794_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243817770_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243823740_v000.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243829529_v001.fit
root://glast-rdr.slac.stanford.edu//glast/Data/Flight/Level1/LPA/prod/1.65/ft2/gll_pt_r0243834802_v000.fit

ft1file list and ft2file list are consistent
Fractional coverage of FT1 files: 0.773315453599
time -p fmerge infiles="@ft1merge_file_list" outfile="FT1_243756000_243842400.fits" columns=" " mextname=" " copyprime=yes lastkey=" " history=yes clobber=yes mode="ql"
real 3.21
user 0.76
sys 1.49
time -p fmerge infiles="@ft1merge_file_list" outfile="ft1merge_gti.fits" columns=" " mextname=" " copyprime=yes lastkey=" " history=yes clobber=yes mode="ql"
real 0.15
user 0.01
sys 0.01
Overwrite existing file 'FT1_243756000_243842400.fits'.
time -p fchecksum infile="FT1_243756000_243842400.fits" update=yes datasum=yes ok=yes mode="ql"
Checksum keywords updated successfully.
real 0.36
user 0.03
sys 0.05
time -p gtselect infile=FT1_243756000_243842400.fits outfile=tmp.fits ra=180.0 dec=0.0 rad=180.0 tmin=243756000.0 tmax=243842400.0 emin=30.0 emax=500000.0 evclsmin=0 evclsmax=10 convtype=-1 zmax=105.0 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql"
Done.
real 2.45
user 0.56
sys 0.20
time -p fcopy infile="gll_pt_r0243760583_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0000.fits" mode="ql"
real 0.02
user 0.00
sys 0.01
time -p fcopy infile="gll_pt_r0243766312_v002.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0001.fits" mode="ql"
real 0.03
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243772042_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0002.fits" mode="ql"
real 0.03
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243777771_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0003.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243781394_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0004.fits" mode="ql"
real 0.03
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243787695_v000.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0005.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243793786_v000.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0006.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243799818_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0007.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243805812_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0008.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243811794_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0009.fits" mode="ql"
real 0.02
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243817770_v000.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0010.fits" mode="ql"
real 0.04
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243823740_v000.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0011.fits" mode="ql"
real 0.08
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243829529_v001.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0012.fits" mode="ql"
real 0.03
user 0.00
sys 0.00
time -p fcopy infile="gll_pt_r0243834802_v000.fit[SC_DATA][LIVETIME>0]" outfile="ft2_filtered_0013.fits" mode="ql"
real 0.03
user 0.00
sys 0.00
time -p fmerge infiles="@ft2_file_list" outfile="ft2.fits" columns=" " mextname=" " copyprime=yes lastkey=" " history=yes clobber=yes mode="ql"
real 0.15
user 0.03
sys 0.01
Overwrite existing file 'ft2.fits'.
time -p gtmktime scfile=ft2.fits sctable="SC_DATA" filter="IN_SAA!=T && LIVETIME>0" roicut="yes" evfile=tmp.fits evtable="EVENTS" outfile="ft1.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 chatter=2 clobber=yes debug=no gui=no mode="ql"
DSS keywords required for ROI-based zenith angle cut
are not present in the FT1 file.
real 2.28
user 0.31
sys 0.15

01a

ft1.fits ft1list.list ft2.fits ft2list.list

02

rm: cannot lstat `ft1_diffuse.fits': No such file or directory
fcopy infile=ft1.fits"[EVENTS][CTBCLASSLEVEL>2]" outfile=ft1_diffuse.fits
Output File is: ft1_diffuse.fits

03

working ft2 file: ft2.fits

['START', 'STOP', 'SC_POSITION', 'LAT_GEO', 'LON_GEO', 'RAD_GEO', 'RA_ZENITH', 'DEC_ZENITH', 'B_MCILWAIN', 'L_MCILWAIN', 'GEOMAG_LAT', 'IN_SAA', 'RA_SCZ', 'DEC_SCZ', 'RA_SCX', 'DEC_SCX', 'LAT_MODE', 'LIVETIME', 'QSJ_1', 'QSJ_2', 'QSJ_3', 'QSJ_4']

03a

03b

ft1.fits ft1_diffuse.fits ft1list.list ft2.fits ft2list.list

03c

03d

03e

04

time -p gtselect infile=ft1_diffuse.fits outfile=region_0000_evt.fits ra=343.6566 dec=16.1494 rad=15.0 tmin=243760585.99 tmax=243840783.084 emin=100.0 emax=300000.0 evclsmin=0 evclsmax=10 convtype=-1 zmax=105.0 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql"
Done.
real 0.90
user 0.21
sys 0.03
time -p gtmktime scfile=ft2.fits sctable="SC_DATA" filter="angsep(RA_ZENITH,DEC_ZENITH,343.657, 16.149)+ 15.000<105.000" roicut="yes" evfile=region_0000_evt.fits evtable="EVENTS" outfile="foo.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 chatter=2 clobber=yes debug=no gui=no mode="ql"
real 0.76
user 0.16
sys 0.03
time -p gtselect infile=ft1_diffuse.fits outfile=region_0000_pgw_evt.fits ra=343.6566 dec=16.1494 rad=37.5 tmin=243760585.99 tmax=243840783.084 emin=100.0 emax=300000.0 evclsmin=0 evclsmax=10 convtype=-1 zmax=105.0 phasemin=0.0 phasemax=1.0 evtable="EVENTS" chatter=2 clobber=yes debug=no gui=no mode="ql"
Done.
real 0.89
user 0.21
sys 0.04
time -p gtbin evfile=region_0000_pgw_evt.fits scfile=NONE outfile=region_0000_pgw_map.fits algorithm="CMAP" ebinalg="LOG" ebinfile=NONE tbinalg="LIN" tbinfile=NONE nxpix=50 nypix=50 binsz=0.5 coordsys="CEL" xref=343.6566 yref=16.1494 axisrot=0.0 rafield="RA" decfield="DEC" proj="AIT" evtable="EVENTS" sctable="SC_DATA" efield="ENERGY" tfield="TIME" chatter=2 clobber=yes debug=no gui=no mode="ql"
This is gtbin version v2r0p5
real 1.13
user 0.22
sys 0.04
Counts Map Plot saved in: Counts_Map_343_16.png
time -p pgwave2D input_file="region_0000_pgw_map.fits" bgk_choise="n" input_bgk_file="" circ_square="s" N_iterations=1 SN_ratio="0" N_scale=1 scala="3" otpix="10" n_sigma=7.0 median_box=9.0 r_threshold=0.5 kappa=3.0 min_pix=3.0 m_num=10 border_size=4 s_box=5 fitsio_choice="n" recursive_choice="n" verbose_level=0

1. Initialization and dynamic memory allocation. Execution in progress...

2. External Image Loading... done.

3. *** SCAN (1) of (1) ***
		A. Gaussian Filter... Scan (1) of (1)
	3.1. Scale (1) of (1) with MH scale-factor set to ** 3 channels **
done.
		B. Median Filter... done.
		C. WAVELET TRANSFORM MAP building (WT of the input image)...done.
		D. THRESHOLD and OVER-THRESHOLD MAP building... done.
		E. Candidate-sources search... done.

	-> Number of identified candidate-sources = 1


	3.2Selection of sources found in two consecutive scales


4. Removal of double detections and border-sources.

****************************************************************************
5. DETECTED SOURCES: 1

warning: could not find bounding box

 Detected sources DS9 regions written on file: 	region_0000_pgw_map.reg

 Detected sources list written on file: 	region_0000_pgw_map.list

6. Clean.

7. Images ready, esecution ended.

real 2.13
user 1.44
sys 0.05
smoothed counts map + sources saved in region_0000_pgw_map_cmap.png
['PGW_0001']
Overwrite existing file 'region_0000_pgwsrc.fits'.
XML Model saved in: region_0000.xml

05

Counts_Map_343_16.png	   OTM_scan1_scale1.fits  foo.fits  ft1list.list  ft2list.list	region_0000.xml region_0000_pgw_evt.fits	region_0000_pgw_map.list  region_0000_pgwsrc.fits
Detected_scan1_scale1.reg  WT_scan1_scale1.fits   ft1.fits  ft2.fits	  iLat.history	region_0000_evt.fits  region_0000_pgw_map.fits	region_0000_pgw_map.reg

06

#[ID]   [X]	   [Y]     [RA]        [DEC]     [POS_ERR]     [SNR]     [K-signf]  [Counts]  [SigC]  [BKG]  [SigBkg]
   1    25.0     24.0     343.3967      15.8992    0.675       7.782       16.381       81      10       3       7

07

<?xml version="1.0" ?>
<source_library title="source library">
<source name="PGW_0001" type="PointSource">
<spectrum type="PowerLaw2">
<parameter error="0.00" free="1" max="1000" min="1e-05" name="Integral" scale="1e-06" value="2.000"/>
<parameter error="0.00" free="1" max="0" min="-5" name="Index" scale="1" value="-2.200"/>
<parameter free="0" max="5e5" min="20" name="LowerLimit" scale="1" value="100.000000"/>
<parameter free="0" max="3e5" min="20" name="UpperLimit" scale="1" value="300000.000000"/>
</spectrum>
<spatialModel type="SkyDirFunction">
<parameter free="0" max="360" min="-360" name="RA" scale="1" value="343.3967"/>
<parameter free="0" max="90" min="-90" name="DEC" scale="1" value="15.8992"/>
</spatialModel>
</source>
<source name="Galactic Diffuse" type="DiffuseSource">
<spectrum type="ConstantValue">
<parameter error="0." free="1" max="10.0" min="0.0" name="Value" scale="1.0" value="1.0"/>
</spectrum>
<spatialModel file="$(EXTFILESSYS)/galdiffuse/GP_gamma.fits" type="MapCubeFunction">
<parameter free="0" max="1000.0" min="0.001" name="Normalization" scale="1.0" value="1.0"/>
</spatialModel>
</source>
<source name="Extragalactic Diffuse" type="DiffuseSource">
<spectrum type="PowerLaw">
<parameter error="0.289864" free="1" max="100.0" min="1e-05" name="Prefactor" scale="1e-07" value="1.65213"/>
<parameter free="0" max="-1.0" min="-3.5" name="Index" scale="1.0" value="-2.1"/>
<parameter free="0" max="200.0" min="50.0" name="Scale" scale="1.0" value="100.0"/>
</spectrum>
<spatialModel type="ConstantValue">
<parameter free="0" max="10.0" min="0.0" name="Value" scale="1.0" value="1.0"/>
</spatialModel>
</source>
</source_library>

08

time -p gtltcube evfile="region_0000_evt.fits" evtable="EVENTS" scfile=ft2.fits sctable="SC_DATA" outfile=expCube_0000.fits dcostheta=0.025 binsize=1.0 tmin=0.0 tmax=0.0 file_version="1" zmax=180.0 chatter=2 clobber=yes debug=no gui=no mode="ql"
Working on file ft2.fits
.....................!
real 18.41
user 16.85
sys 0.33
Livetime Cube saved in: expCube_0000.fits
expCube file created

09

Please wait this task will take several minutes .....
time -p gtdiffrsp evfile=region_0000_evt.fits evtable="EVENTS" scfile=ft2.fits sctable="SC_DATA" srcmdl=diffuse_0000.xml irfs="P6_V1_DIFFUSE" edisp=no chatter=2 clobber=no debug=no gui=no mode="ql"
adding source Extragalactic Diffuse
adding source Galactic Diffuse
Working on...
region_0000_evt.fits.....................!
real 117.83
user 114.71
sys 0.33

10

time -p gtexpmap evfile=region_0000_evt.fits evtable="EVENTS" scfile=ft2.fits sctable="SC_DATA" expcube=expCube_0000.fits outfile=expMap_0000.fits irfs="P6_V1_DIFFUSE" edisp=no srcrad=25.0 nlong=100 nlat=100 nenergies=10 submap=no nlongmin=0 nlongmax=0 nlatmin=0 nlatmax=0 chatter=2 clobber=yes debug=no gui=no mode="ql"
The exposure maps generated by this tool are meant
to be used for *unbinned* likelihood analysis only.
Do not use them for binned analyses.

11

Counts_Map_343_16.png	   WT_scan1_scale1.fits  expMap_0000.fits  ft1_diffuse.fits  ft2list.list     region_0000_evt.fits	region_0000_pgw_map.list
Detected_scan1_scale1.reg  diffuse_0000.xml	 foo.fits	   ft1list.list      old	      region_0000_pgw_evt.fits	region_0000_pgw_map.reg
OTM_scan1_scale1.fits	   expCube_0000.fits	 ft1.fits	   ft2.fits	     region_0000.xml  region_0000_pgw_map.fits	region_0000_pgwsrc.fits

11a

12

Extragalactic Diffuse

  Spectrum: PowerLaw

  0      Prefactor:  1.763e+00  7.186e-01  1.000e-05  1.000e+02 ( 1.000e-07)

  1          Index: -2.100e+00  0.000e+00 -3.500e+00 -1.000e+00 ( 1.000e+00) fixed

  2          Scale:  1.000e+02  0.000e+00  5.000e+01  2.000e+02 ( 1.000e+00) fixed

Galactic Diffuse

  Spectrum: ConstantValue

  3          Value:  6.678e-01  4.196e-01  0.000e+00  1.000e+01 ( 1.000e+00)
PGW_0001

  Spectrum: PowerLaw2

  4       Integral:  2.742e+00  3.541e-01  1.000e-05  1.000e+03 ( 1.000e-06)

  5          Index: -2.396e+00  1.246e-01 -5.000e+00  0.000e+00 ( 1.000e+00)

  6     LowerLimit:  1.000e+02  0.000e+00  1.000e+01  5.000e+05 ( 1.000e+00) fixed

  7     UpperLimit:  3.000e+05  0.000e+00  1.000e+01  5.000e+05 ( 1.000e+00) fixed

#SRC   Tstart Tstop emin  emax TS Npred flux fluerr index indexerr

PGW_0001 243760616.6 243840775.6 100.0 300000.0 354.917541034 103.067280805 2.7419872173e-06 3.54060171614e-07 -2.39615292904 0.124612134072


  Unbinned likelihood results saved in: like_result_0000.xml

  ['ObsCounts', 'Extragalactic Diffuse', 'Galactic Diffuse', 'PGW_0001']

  Like counts spectra saved in :  like_specplot_0000.png

13

3C454-like-results.dat	   WT_scan1_scale1.fits  foo.fits	   ft2.fits		   like_spectrum_0000.fits  region_0000_pgw_evt.fits  region_0000_pgwsrc.fits
Counts_Map_343_16.png	   diffuse_0000.xml	 ft1.fits	   ft2list.list		   old			    region_0000_pgw_map.fits
Detected_scan1_scale1.reg  expCube_0000.fits	 ft1_diffuse.fits  like_result_0000.xml    region_0000.xml	    region_0000_pgw_map.list
OTM_scan1_scale1.fits	   expMap_0000.fits	 ft1list.list	   like_specplot_0000.png  region_0000_evt.fits     region_0000_pgw_map.reg

14

PGW_0001 243760616.6 243840775.6 100.0 300000.0 354.917541034 103.067280805 2.7419872173e-06 3.54060171614e-07 -2.39615292904 0.124612134072   False

[PGW-src-ID  Tstart     Tstop       Emin  Emax     TS(test stat)   Npred          Flux              Flux error         index         index error    Upper Limit]

Flux is the "Integral" parameter for PowerLaw2 model scaled for the "scale" value in the xml model file

  Npred is the number of photons predicted by the spectral model

  TS is the Test Statistic value

15

setdcdir /Data/Flight/Level1/LPA
setwd ./
getdcft 243756000.0  243842400.0 ft1.fits ft2.fits
$ls
setft2 ft2.fits
setft1 ft1.fits 
setevtcl 2
setft1 ft1_diffuse.fits
setirfs P6_V1_DIFFUSE
setene 100 300000 
setroi 343.6566  16.1494 15 25
$ls
xmlmodel 0.5
ltcube 
diffuse
expmap
$ ls
setulike region_0000_evt.fits expCube_0000.fits expMap_0000.fits region_0000.xml
ulike 3C454-like-results.dat
$ls
$cat 3C454-like-results.dat
quit