| In [2]: execfile("scripts/GRBREPORT.py") |
| Enter the GRBTRIGGERDATE (MET or YYYY-MM-DD HH:MM:SS [2008-08-03 18:31:23]243216766 Enter the TSTART with respect the trigget time (s) [-10] Enter the TSTOP with respect the trigget time (s) [60]300 Enter the RA of the GRB (deg) [281.833]119.88 Enter the DEC of the GRB (deg) [81.55]-56.6 Enter the R.O.I. (deg) [15] Enter the Minimum Energy (MeV) [10]100 Enter the Maximum Energy (MeV) [600000]100000 Enter the Number of Log bins [10] Enter the IRF [P6_V3_TRANSIENT]P6_V3_DIFFUSE Optional: FT1 file (or look in $INDIR):[] Optional: FT2 file (or look in $INDIR):[] Enter the Localization Error at 1-sigma (deg) [0.1] |
Note: In the example above, the "Optional: FT1 ... and FT1 ..." fields are left blank since the data was moved to the DATA/FITS/LAT subdirectory. These "Optional" fields provide a convenient method of creating a path to wherever you might prefer to place your data.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Updating the PAR file ... 243216766.0 243216766.0 080916009 243216766.0 243216766.0 080916009 Trigger according to GRB name : 243216778.0 GBM Directory not found DATA/FITS/GBM/GRB080916009 OUTPUT DIRECTORY: DATA/GRBOUT/080916009 <GTGRB.GRB.GRB instance at 0xb69f40cc> Searching in directory DATA/FITS/LAT File: DATA/FITS/LAT/ft1_243215766_243217766.fits TSTART=243215766.0 TSTOP=243217766.0 TTRIG.=243216766. found file DATA/FITS/LAT/ft1_243215766_243217766.fits Searching in directory DATA/FITS/LAT File: DATA/FITS/LAT/ft1_243215766_243217766.fits TSTART=243215766.0 TSTOP=243217766.0 TTRIG.=243216766.0 File: DATA/FITS/LAT/ft2_243215766_243217766.fits TSTART=243215766.0 TSTOP=243217766.6 TTRIG.=243216766.0 found file DATA/FITS/LAT/ft2_243215766_243217766.fits FT1 : DATA/FITS/LAT/ft1_243215766_243217766.fits FT2 : DATA/FITS/LAT/ft2_243215766_243217766.fits LAT DIR : DATA/FITS/LAT Creating root file... Open DATA/GRBOUT/080916009/080916009.root %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% In file: TMIN: 243215788.056 ( 2008-09-15T23:56:27.056210 ) TMAX: 243217765.82 ( 2008-09-16T00:29:24.819840 ) Emin: 10.5647 Emax: 926549.0 GRB Trigger time: 243216766.0 ( 2008-09-16T00:12:45 ) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Theta GRB = 48.8353293594 at time: -0.40000000596 GRBname: 080916009 T start: 243216756.0 T end: 243217066.0 (Duration: 310.0 ) Ra : 119.88 Dec : -56.6 THETA : 48.8353293594 Rad : 15.0 Emin : 100.0 Emax : 100000.0 IRFs: P6_V3_DIFFUSE
| plot angular separation? [y/N]y |
plot angular separation? [y/N]y
in DATA/FITS/LAT/ft2_243215766_243217766.fits
-978.011959434, 999.599999994, 1979
Time: 243216765.6 (-0.40000000596)
GRB Ra, dec= (119.9,-56.6)
SC Ra, dec= (75.6,-20.2)
Zenith Ra, dec= (62.0,12.1)
Lat in SAA? False
GRB THETA : 48.8353293594
FROM ZENITH : 83.6481497046
press enter to continue
Tips:
press enter to continue
| make select? [y/N]y ZMAX: [180]105 |
time -p fselect infile="DATA/FITS/LAT/ft1_243215766_243217766.fits"
outfile="DATA/GRBOUT/080916009/080916009_LAT_SEL.fits"
expr="EVENT_CLASS>2"
index="OBSOLETE"
histkw=yes
copyall=yes
keycopy=yes
clobber=yes
nrows=128 mode="ql"
real 4.79
user 0.01
sys 0.04
Selected : 8077 Events...
Selected : 817 Events...
time -p gtselect infile=DATA/GRBOUT/080916009/080916009_LAT_SEL.fits
outfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits ra=119.88 dec=-56.6
rad=15.0 tmin=243216756.0
tmax=243217066.0 emin=100.0
emax=100000.0
zmax=105.0
evclsmin=0
evclsmax=10
convtype=-1 phasemin=0.0
phasemax=1.0
evtable="EVENTS"
chatter=2
clobber=yes
debug=no gui=no mode="ql"
Done.
real 32.79
user 0.36
sys 0.29
Selected : 133 Events..
| Apply GTMKTTIME? [y/N]:y |
time -p gtmktime scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits sctable="SC_DATA" filter="IN_SAA!=T && LIVETIME>0 && (ANGSEP(RA_ZENITH,DEC_ZENITH,119.88,-56.6) + 15.0 < 105.0)" roicut=no evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits evtable="EVENTS" outfile="tmp.fits" apply_filter=yes overwrite=no header_obstimes=yes tstart=0.0 tstop=0.0 gtifile="default" chatter=2 clobber=yes debug=no gui=no mode="ql"
real 5.30
user 0.30
sys 0.08
CVS c0.png ipy_profile_glastgrb.py logfiles source_model.xml
DATA grblog.par ipy_profile_test.py matplotlibrc templateAnalysis.py
GRB080702B.py gtgrb_example ipy_user_conf.py logfiles opssim2_GRB090308937.py
GTGRB gtgrb_slac ipythonrc scripts
| Apply an energy dependent selection? [y/N]:N Make PHA2 and RSP [y/N]:N |
Note: If you plan to run rmfit, set the "Make PHA2 and RSP" option to "y"; gtgrb will make the pha2 files and save them in the correct format to be used in rmfit.
Now making a root tree with the selected events...
...opening DATA/GRBOUT/080916009/080916009_LAT_ROI.fits
Filename: DATA/GRBOUT/080916009/080916009_LAT_ROI.fits
No. Name Type Cards Dimensions Format
0 PRIMARY PrimaryHDU 31 () uint8
1 EVENTS BinTableHDU 176 133R x 22C [E, E, E, E, E, E, E, E, E, D, J, J, I, 3I, I, I, D, E, E, E, E, E]
2 GTI BinTableHDU 48 1R x 2C [D, D]
Find 133 events...
imported root!
243216769.645
243216770.588
243216771.315
243216771.381
243216771.459
243216771.616
243216772.174
243216772.265
243216772.269
243216772.296
| LIGHT CURVE BIN SIZE (seconds):1 |
===> 243216756.0 243217066.0 310.0
Make lightCurve for detector: LAT from 243216756.0 to 243217066.0
time -p gtbin evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT_lc.fits algorithm="LC" ebinalg="LOG" emin=30.0 emax=200000.0 ebinfile=NONE tbinalg="LIN" tstart=243216756.0 tstop=243217066.0 dtime=1.0 tbinfile=NONE coordsys="CEL" 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 v2r3p2
real 44.60
user 0.48
sys 0.08
Filename: DATA/GRBOUT/080916009/080916009_LAT_lc.fits
No. Name Type Cards Dimensions Format 0 PRIMARY PrimaryHDU 18 () uint8 1 RATE BinTableHDU 65 310R x 4C [D, D, J, E] 2 GTI BinTableHDU 41 1R x 2C [D, D]
NUMBER OF LAT LAT COUNTS: 133 Saved histogram in DATA/GRBOUT/080916009/080916009.root
plotting lightcurve from DATA/GRBOUT/080916009/080916009_LAT_lc.fits Info in <TCanvas::Print>: GIF file DATA/GRBOUT/080916009/080916009_LAT_lc.gif has been created SKY MAP Bin size (degrees):Info in <TCanvas::Print>: file ./c0.png has been created
Tip: If a blank window is displayed, click in it to refresh the screen.
| SKY MAP Bin size (degrees):1 |
ROI/binz: 15.0
Make Sky Map: 100.0 100000.0 10 LOG
binsz 1.0
nypix 15
nxpix 15
time -p gtbin evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT_map.fits algorithm="CMAP" ebinalg="LOG" emin=30.0 emax=200000.0 enumbins=0 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=243216756.0 tstop=243217066.0 dtime=1.0 tbinfile=NONE snratio=0.0 lcemin=0.0 lcemax=0.0 nxpix=15 nypix=15 binsz=1.0 coordsys="CEL" xref=119.88 yref=-56.6 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 v2r3p2
real 3.35
user 0.50
sys 0.05
| is ~/ds9 available? [y/N]N make GBM light curve? [y/N]N make LIKELIHOOD ANALYSIS? [y/N]y Include DIFFUSE COMPONENT? [y/N]y create XML MODEL? [y/N]y |
Creating xml file...
...adding a Point Source...
...adding the diffuse component...
... closing the xml file
time -p /nfs/farm/g/glast/u30/builds/rh9_gcc32opt/ScienceTools/ScienceTools-v9r15p4/Likelihood/v15r3p2/rh9_gcc32opt/gtltcube.exe evfile="DATA/GRBOUT/080916009/080916009_LAT_ROI.fits" evtable="EVENTS" scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits sctable="SC_DATA" outfile=DATA/GRBOUT/080916009/080916009_LAT_expCube.fits dcostheta=0.025 binsz=0.5 phibins=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 DATA/FITS/LAT/ft2_243215766_243217766.fits
.....................!
real 57.32
user 27.76
sys 7.10
time -p /nfs/farm/g/glast/u30/builds/rh9_gcc32opt/ScienceTools/ScienceTools-v9r15p4/Likelihood/v15r3p2/rh9_gcc32opt/gtexpmap.exe evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits evtable="EVENTS" scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits sctable="SC_DATA" expcube=DATA/GRBOUT/080916009/080916009_LAT_expCube.fits outfile=DATA/GRBOUT/080916009/080916009_LAT_expMap.fits irfs="P6_V3_DIFFUSE" edisp=no srcrad=25.0 nlong=100 nlat=100 nenergies=6 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.
Computing the ExposureMap using DATA/GRBOUT/080916009/080916009_LAT_expCube.fits
....................!
real 172.91
user 164.27
sys 0.90
| COMPUTE GTDIFFRESP ? [y/N]N Make likelihood fit? [y/N]y |
Note: Since the data was obtained from the Astro Server, the diffuse response does not have to be computed.
evtfile= DATA/GRBOUT/080916009/080916009_LAT_ROI.fits
scfile= DATA/FITS/LAT/ft2_243215766_243217766.fits
expmap= DATA/GRBOUT/080916009/080916009_LAT_expMap.fits
expcube= DATA/GRBOUT/080916009/080916009_LAT_expCube.fits
irf= P6_V3_DIFFUSE
xmlFile= source_model.xml
optimizer= NEWMINUIT
DATA/GRBOUT/080916009/080916009_LAT_spectrum.fits
Info in <Minuit2>: VariableMetricBuilder: matrix not pos.def, gdel > 0
Info in <Minuit2>: gdel = 4119.12
Info in <Minuit2>: negative or zero diagonal element in covariance matrix : i = 3
Info in <Minuit2>: added to diagonal of Error matrix a value : dg = 0.931739
Info in <Minuit2>: gdel = -53714.4
Minuit did successfully converge.
# of function calls: 95
minimum function Value: 682.900660468
minimum edm: 3.9381174e-06
minimum internal state vector: LAVector parameters:
-1.5708298
11.944436
-0.52779196
-181.7801
minimum internal covariance matrix: LASymMatrix parameters:
0.0081756925 -7.0653195e-05 -3.0830167e-06 1.1915959e-06
-7.0653195e-05 0.044468605 -0.0011844424 0.00048300839
-3.0830167e-06 -0.0011844424 0.0079085608 -0.003409511
1.1915959e-06 0.00048300839 -0.003409511 0.0051679466
# ext. || Name || type || Value || Error +/-
0 ||0_Normalization || limited || 1.0279574e-05 || 1.0216135
1 || 1_Value || limited || 2.0869557 || 0.60371172
2 ||2_Integral || limited || 248.18651 || 27.14511
3 || 3_Index || limited || -2.162119 || 0.092274
# of function calls: 23
function Value: 682.900660468
expected distance to the Minimum (edm): 3.928339e-06
external parameters:
# ext. || Name || type || Value || Error +/-
0 ||0_Normalization || limited || 1.0279574e-05 || 1.0216364
1 || 1_Value || limited || 2.0869557 || 0.60370325
2 ||2_Integral || limited || 248.18651 || 27.104966
3 || 3_Index || limited || -2.162119 || 0.09213705
covariance matrix:
MnUserCovariance:
1.14288e-06 3.63053e-07 1.66916e-06 -2.68405e-09
3.63053e-07 0.367171 -1.03585 0.00176055
1.66916e-06 -1.03585 735.645 -1.32834
-2.68405e-09 0.00176055 -1.32834 0.00849653
MnUserCovariance Parameter correlations:
1 0.000560448 5.75657e-05 -2.72376e-05
0.000560448 1 -0.0630271 0.0315205
5.75657e-05 -0.0630271 1 -0.531317
-2.72376e-05 0.0315205 -0.531317 1
global correlation coefficients :
MnGlobalCorrelationCoeff:
MnGlobalCorrelationCoeff:
0.000568149
0.0630723
0.533331
0.531321
log(like)= 682.900670123
Flux = 0.120057 (MeV/cm^2/s)
Flux = 0.000249 (ph/cm^2/s)
 |
 |
| make XSPEC ANALYSIS? [y/N]y compute LAT spectrum (XSPEC)? [y/N]y |
Make PHA1 FILE: 100.0 100000.0 10 LOG
time -p gtbin evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT.pha algorithm="PHA1" ebinalg="LOG" emin=100.0 emax=100000.0 enumbins=10 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=243216756.0 tstop=243217066.0 dtime=1.0 tbinfile=NONE snratio=0.0 lcemin=0.0 lcemax=0.0 nxpix=15 nypix=15 binsz=1.0 coordsys="CEL" xref=119.88 yref=-56.6 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 v2r3p2
real 1.20
user 0.49
sys 0.07
Generating LAT response matrix... 100.0 100000.0 100 LOG AT: 243216766.0
time -p gtrspgen respalg="GRB" specfile=DATA/GRBOUT/080916009/080916009_LAT.pha scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT.RSP irfs="P6_V3_DIFFUSE" sctable="SC_DATA" resptpl=DEFAULT chatter=2 clobber=yes debug=no gui=no mode="ql" time=243216756.0 thetacut=60.0 dcostheta=0.05 phinumbins=1 ebinalg="LOG" efield="ENERGY" emin=100.0 emax=100000.0 enumbins=100 denergy=50.0 ebinfile=NONE
This is gtrspgen version v3r1
real 5.53
user 3.29
sys 0.54
LAT
executing... DATA/GRBOUT/080916009/LAT_fit.xcm
***Warning: You do not have an up-to-date version of Xspec.init
in your ~/.xspec directory. It is recommended that you move
your old Xspec.init aside and start up XSPEC again. This will
place a new version Xspec.init in your directory, which you may
then modify using your old settings.
XSPEC version: 12.2.1
Build Date/Time: Tue Nov 22 14:34:40 2005
XSPEC12>@Make PHA1 FILE: 100.0 100000.0 10 LOG
time -p gtbin evfile=DATA/GRBOUT/080916009/080916009_LAT_ROI.fits scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT.pha algorithm="PHA1" ebinalg="LOG" emin=100.0 emax=100000.0 enumbins=10 denergy=0.0 ebinfile=NONE tbinalg="LIN" tstart=243216756.0 tstop=243217066.0 dtime=1.0 tbinfile=NONE snratio=0.0 lcemin=0.0 lcemax=0.0 nxpix=15 nypix=15 binsz=1.0 coordsys="CEL" xref=119.88 yref=-56.6 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 v2r3p2
real 1.20
user 0.49
sys 0.07
Generating LAT response matrix... 100.0 100000.0 100 LOG AT: 243216766.0
time -p gtrspgen respalg="GRB" specfile=DATA/GRBOUT/080916009/080916009_LAT.pha scfile=DATA/FITS/LAT/ft2_243215766_243217766.fits outfile=DATA/GRBOUT/080916009/080916009_LAT.RSP irfs="P6_V3_DIFFUSE" sctable="SC_DATA" resptpl=DEFAULT chatter=2 clobber=yes debug=no gui=no mode="ql" time=243216756.0 thetacut=60.0 dcostheta=0.05 phinumbins=1 ebinalg="LOG" efield="ENERGY" emin=100.0 emax=100000.0 enumbins=100 denergy=50.0 ebinfile=NONE
This is gtrspgen version v3r1
real 5.53
user 3.29
sys 0.54
LAT
executing... DATA/GRBOUT/080916009/LAT_fit.xcm
***Warning: You do not have an up-to-date version of Xspec.init
in your ~/.xspec directory. It is recommended that you move
your old Xspec.init aside and start up XSPEC again. This will
place a new version Xspec.init in your directory, which you may
then modify using your old settings.
XSPEC version: 12.2.1
Build Date/Time: Tue Nov 22 14:34:40 2005
| XSPEC12>@DATA/GRBOUT/080916009/LAT_fit.xcm |
Note: In this step, the "@" plus the path and .xcm file created by XSPEC (see path and file name highlighted in red above) is entered in order to generate the PGPLOT: data and folded model shown below.
!XSPEC12>data DATA/GRBOUT/080916009/080916009_LAT.pha;
1 spectrum in use
Source File: DATA/GRBOUT/080916009/080916009_LAT.pha
Net count rate (cts/s) for Spectrum:1 4.646e-01 +/- 4.097e-02
Assigned to Data Group 1 and Plot Group 1
Noticed Channels: 1-10
Telescope: GLAST Instrument: LAT Channel Type: PI
Exposure Time: 286.2 sec
Using Response (RMF) File DATA/GRBOUT/080916009/080916009_LAT.RSP
!XSPEC12>response DATA/GRBOUT/080916009/080916009_LAT.RSP;
Response successfully loaded.
!XSPEC12>backgrnd ;
reset background files
Syntax: backgrnd none
backgrnd [spectrum number] <filename{range}>
(many options: see help for full syntax)
!XSPEC12>setplot energy
!XSPEC12>model pow
/*
========================================================================
Model powerlaw<1> Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
1 1 powerlaw PhoIndex 1.00000 +/- 0.0
2 1 powerlaw norm 1.00000 +/- 0.0
________________________________________________________________________
Chi-Squared = 3.906198e+12 using 10 PHA bins.
Reduced chi-squared = 4.882748e+11 for 8 degrees of freedom
Null hypothesis probability = 0.000000e+00
Valid fit does not exist.
!XSPEC12>statistic cstat
C-statistic = 1.488200e+07 using 10 PHA bins.
Valid fit does not exist.
!XSPEC12>method leven
C-statistic = 1.488200e+07 using 10 PHA bins.
Valid fit does not exist.
!XSPEC12>fit 100
C-Statistic Lvl Par # 1 2
267.874 -1 1.0731 4.90616e-05
234.964 -1 1.12812 0.000123284
204.583 -1 1.18607 0.000300005
178.635 -1 1.24035 0.000688845
156.258 -1 1.2918 0.00147917
137.421 -1 1.33922 0.00292115
122.149 -1 1.38168 0.00522972
110.559 -1 1.41834 0.00845398
102.32 -1 1.44908 0.0124678
95.3362 -1 1.47452 0.0176494
88.6438 -1 1.4975 0.0246539
82.5159 -1 1.51958 0.0339866
76.9266 -1 1.54087 0.0462527
71.8304 -1 1.5614 0.0621646
67.1829 -1 1.58118 0.0825482
62.9435 -1 1.60023 0.108347
59.0746 -1 1.61857 0.140625
55.5421 -1 1.63623 0.180565
52.3147 -1 1.65322 0.22947
49.3642 -1 1.66957 0.288754
46.6649 -1 1.6853 0.359934
44.1935 -1 1.70043 0.444625
41.9291 -1 1.71499 0.544524
39.8526 -1 1.729 0.661396
37.947 -1 1.74248 0.79706
36.1967 -1 1.75544 0.953372
34.5879 -1 1.76792 1.13221
33.1078 -1 1.77993 1.33544
31.7452 -1 1.79149 1.56494
30.4897 -1 1.80263 1.82252
29.332 -1 1.81334 2.10996
28.2636 -1 1.82367 2.42895
27.2771 -1 1.83361 2.78113
26.3654 -1 1.8432 3.168
25.5223 -1 1.85243 3.59097
24.7421 -1 1.86133 4.05133
24.0197 -1 1.86992 4.55022
23.3502 -1 1.87819 5.08865
22.7296 -1 1.88617 5.66746
22.1538 -1 1.89387 6.28736
21.6193 -1 1.9013 6.94889
21.1229 -1 1.90847 7.65242
20.6616 -1 1.91539 8.39818
20.2327 -1 1.92208 9.18623
19.8338 -1 1.92853 10.0165
19.4625 -1 1.93476 10.8886
19.1168 -1 1.94078 11.8023
18.7947 -1 1.94659 12.757
18.4946 -1 1.95221 13.752
18.2147 -1 1.95764 14.7865
17.9537 -1 1.96289 15.8596
17.7102 -1 1.96797 16.9701
17.4828 -1 1.97288 18.117
17.2705 -1 1.97762 19.2991
17.0722 -1 1.98221 20.5149
16.8869 -1 1.98665 21.763
16.8606 -2 2.02387 32.4475
15.5804 -2 2.04538 44.5533
14.8866 -2 2.06319 57.6485
14.5192 -2 2.07761 70.7762
14.3252 -2 2.0891 83.1171
14.2214 -2 2.09813 94.1343
14.1647 -2 2.10516 103.584
14.1332 -2 2.11058 111.448
14.1154 -2 2.11475 117.846
14.1144 -3 2.12772 138.076
==============================
Variances and Principal Axes
1 2
4.37E-05 | -1.00 0.00
2.62E+04 | 0.00 -1.00
------------------------------
========================
Covariance Matrix
1 2
1.097e-02 1.693e+01
1.693e+01 2.622e+04
------------------------
========================================================================
Model powerlaw<1> Source No.: 1 Active/On
Model Model Component Parameter Unit Value
par comp
1 1 powerlaw PhoIndex 2.12772 +/- 0.104759
2 1 powerlaw norm 138.076 +/- 161.923
________________________________________________________________________
C-statistic = 14.11 using 10 PHA bins.
!XSPEC12>tclout model
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set mod [string trim $xspec_tclout]
!XSPEC12>tclout energies
!XSPEC12>lindex $xspec_tclout 0
!XSPEC12>set emin [lindex $xspec_tclout 0]
!XSPEC12>lindex $xspec_tclout end
!XSPEC12>set emax [lindex $xspec_tclout end]
!XSPEC12>open DATA/GRBOUT/080916009/LAT_fit.dat w
!XSPEC12>set fileid [open DATA/GRBOUT/080916009/LAT_fit.dat w]
!XSPEC12>puts $fileid "===== fit result"
!XSPEC12>tclout model
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set mod [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$mod"
!XSPEC12>puts $fileid "----- parameters"
!XSPEC12>tcloutr modcomp
!XSPEC12>source /afs/slac/g/glast/applications/astroTools/headas/i686-pc-linux-gnu-libc2.2.4/../spectral/scripts/tcloutr.tcl
!XSPEC12>proc tcloutr args {
global xspec_tclout
eval "tclout $args"
return $xspec_tclout
}
!XSPEC12>tcloutr modcomp
!XSPEC12>set comps [tcloutr modcomp]
!XSPEC12>tcloutr modpar
!XSPEC12>set pars [tcloutr modpar]
!XSPEC12>for {set i 1} {$i <= $pars} {incr i} {
tclout pinfo $i
set pn [string trim $xspec_tclout]
if {$comps>1} then {
puts $fileid "$i$pn"
} else {
puts $fileid "$pn"
}
tclout param $i
set par [string trim $xspec_tclout]
puts $fileid "$par"
set num [llength $xspec_tclout]
if {$num > 2} then {
tclout error $i
set err [string trim $xspec_tclout]
puts $fileid "$err"
} else {
puts $fileid "fixed"
}
}
!XSPEC12>set i 1
!XSPEC12>tclout pinfo $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set pn [string trim $xspec_tclout]
!XSPEC12>if {$comps>1} then {
puts $fileid "$i$pn"
} else {
puts $fileid "$pn"
}
!XSPEC12>puts $fileid "$pn"
!XSPEC12>tclout param $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set par [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$par"
!XSPEC12>llength $xspec_tclout
!XSPEC12>set num [llength $xspec_tclout]
!XSPEC12>if {$num > 2} then {
tclout error $i
set err [string trim $xspec_tclout]
puts $fileid "$err"
} else {
puts $fileid "fixed"
}
!XSPEC12>tclout error $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set err [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$err"
!XSPEC12>incr i
!XSPEC12>tclout pinfo $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set pn [string trim $xspec_tclout]
!XSPEC12>if {$comps>1} then {
puts $fileid "$i$pn"
} else {
puts $fileid "$pn"
}
!XSPEC12>puts $fileid "$pn"
!XSPEC12>tclout param $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set par [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$par"
!XSPEC12>llength $xspec_tclout
!XSPEC12>set num [llength $xspec_tclout]
!XSPEC12>if {$num > 2} then {
tclout error $i
set err [string trim $xspec_tclout]
puts $fileid "$err"
} else {
puts $fileid "fixed"
}
!XSPEC12>tclout error $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set err [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$err"
!XSPEC12>incr i
!XSPEC12>puts $fileid "----- covariance"
!XSPEC12>tclout covariance
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set cov [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$cov"
!XSPEC12>puts $fileid "----- simpars"
!XSPEC12>tclout simpars
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set god [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$god"
!XSPEC12>puts $fileid "----- statistics"
!XSPEC12>tclout stat
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set stat [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$stat"
!XSPEC12>puts $fileid "----- dof"
!XSPEC12>tclout dof
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set dof [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$dof"
!XSPEC12>puts $fileid "-----"
!XSPEC12>puts $fileid "^^^^^"
!XSPEC12>puts $fileid "===== spectrum"
!XSPEC12>tcloutr datasets
!XSPEC12>set nospec [tcloutr datasets]
!XSPEC12>for {set i 1} {$i <= $nospec} {incr i} {
puts $fileid "----- spectrum $i"
tclout plot ld x $i
set e [string trim $xspec_tclout]
puts $fileid "$e"
tclout plot ld xerr $i
set ee [string trim $xspec_tclout]
puts $fileid "$ee"
tclout plot ld y $i
set y [string trim $xspec_tclout]
puts $fileid "$y"
tclout plot ld yerr $i
set ye [string trim $xspec_tclout]
puts $fileid "$ye"
plot ld ratio $i
}
!XSPEC12>set i 1
!XSPEC12>puts $fileid "----- spectrum $i"
!XSPEC12>tclout plot ld x $i
!XSPEC12>plot ld
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set e [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$e"
!XSPEC12>tclout plot ld xerr $i
!XSPEC12>plot ld
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set ee [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$ee"
!XSPEC12>tclout plot ld y $i
!XSPEC12>plot ld
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set y [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$y"
!XSPEC12>tclout plot ld yerr $i
!XSPEC12>plot ld
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set ye [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$ye"
!XSPEC12>plot ld ratio $i
!XSPEC12>incr i
!XSPEC12>puts $fileid "^^^^^"
!XSPEC12>puts $fileid "===== flux"
!XSPEC12>flux 50. 300.
No overlap between matrix range ( 100000, 1e+08 )
and the requested range ( 50, 300 )
!XSPEC12>tcloutr datasets
!XSPEC12>set nospec [tcloutr datasets]
!XSPEC12>for {set i 1} {$i <= $nospec} {incr i} {
tclout flux $i
set god [string trim $xspec_tclout]
puts $fileid "$god"
}
!XSPEC12>set i 1
!XSPEC12>tclout flux $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set god [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$god"
!XSPEC12>incr i
!XSPEC12>flux 100000. 1000000000.
Upper range bound 1e+09 reset by matrix bound to 1e+08
Model Flux 0.00028129 photons (2.3355e-07 ergs/cm^2/s) range (1.0000e+05 - 1.0000e+08 keV)
!XSPEC12>tcloutr datasets
!XSPEC12>set nospec [tcloutr datasets]
!XSPEC12>for {set i 1} {$i <= $nospec} {incr i} {
tclout flux $i
set god [string trim $xspec_tclout]
puts $fileid "$god"
}
!XSPEC12>set i 1
!XSPEC12>tclout flux $i
!XSPEC12>string trim $xspec_tclout
!XSPEC12>set god [string trim $xspec_tclout]
!XSPEC12>puts $fileid "$god"
!XSPEC12>incr i
!XSPEC12>puts $fileid "^^^^^"
!XSPEC12>for {set i 1} {$i <= $nospec} {incr i} {
puts $fileid "0 0 0 0 0 0"
}
!XSPEC12>set i 1
!XSPEC12>puts $fileid "0 0 0 0 0 0"
!XSPEC12>incr i
!XSPEC12>puts $fileid "^^^^^"
!XSPEC12>close $fileid
!XSPEC12>cpd /xw
!XSPEC12>plot ld ratio
XSPEC12>
| XSPEC12>exit Do you really want to exit? (y) y XSPEC: quit compute GBM spectrum (XSPEC)? [y/N]N compute GBM+LAT spectrum (XSPEC)? [y/N]N Compute Upper Limit? [y/N]y Compute the background? [y/N]N |
Creating xml file...
...adding a Point Source...
... closing the xml file
evtfile= DATA/GRBOUT/080916009/080916009_LAT_ROI.fits
scfile= DATA/FITS/LAT/ft2_243215766_243217766.fits
expmap= nonetime -p fselect infile="DATA/GRBOUT/080916009/080916009_LAT_SEL.fits" outfile="DATA/GRBOUT/080916009/080916009_LAT_SEL.fits" expr="EVENT_CLASS>2" index="OBSOLETE" histkw=yes copyall=yes keycopy=yes clobber=yes nrows=128 mode="ql"
real 0.07
user 0.01
sys 0.02
Selected : 817 Events...
Selected : 817 Events...
| Compute upperlimit fixing the background? [y/N]N Compute the UL without background? [y/N]y |
Creating xml file...
...adding a Point Source...
... closing the xml file
evtfile= DATA/GRBOUT/080916009/080916009_LAT_ROI.fits
scfile= DATA/FITS/LAT/080916009_ft2_243215766_243217766.fits
expmap= none
expcube= none
irf= P6_V3_DIFFUSE
xmlFile= source_model.xml
optimizer= NEWMINUIT
DATA/GRBOUT/080916009/080916009_LAT_spectrum.fits
Minuit did successfully converge.
# of function calls: 53
minimum function Value: 732.390741987
minimum edm: 8.3954439e-07
minimum internal state vector: LAVector parameters:
-0.46486434
308.30895
minimum internal covariance matrix: LASymMatrix parameters:
0.0079344502 -0.00317588
-0.00317588 0.0045721411
# ext. || Name || type || Value || Error +/-
0 ||0_Integral || limited || 275.84924 || 28.132389
1 || 1_Index || limited || -2.1610412 || 0.086772616
# of function calls: 10
function Value: 732.390741987
expected distance to the Minimum (edm): 8.4013037e-07
external parameters:
# ext. || Name || type || Value || Error +/-
0 ||0_Integral || limited || 275.84924 || 28.093772
1 || 1_Index || limited || -2.1610412 || 0.086653247
covariance matrix:
MnUserCovariance:
790.302 -1.28036
-1.28036 0.00751449
MnUserCovariance Parameter correlations:
1 -0.525396
-0.525396 1
global correlation coefficients :
MnGlobalCorrelationCoeff:
0.525396
0.525396
log(like)= 732.390747553
Minuit did successfully converge.
# of function calls: 26
minimum function Value: 732.390741219
minimum edm: 7.1795853e-08
minimum internal state vector: LAVector parameters:
-0.46478068
0.43287407
minimum internal covariance matrix: LASymMatrix parameters:
0.007917531 -0.0031593897
-0.0031593897 0.0045614793
# ext. || Name || type || Value || Error +/-
0 ||0_Integral || limited || 275.88663 || 28.103597
1 || 1_Index || limited || -2.161037 || 0.086671368
# of function calls: 10
function Value: 732.390741219
expected distance to the Minimum (edm): 7.1779474e-08
external parameters:
# ext. || Name || type || Value || Error +/-
0 ||0_Integral || limited || 275.88663 || 28.097056
1 || 1_Index || limited || -2.161037 || 0.086651198
covariance matrix:
MnUserCovariance:
790.487 -1.28049
-1.28049 0.00751414
MnUserCovariance Parameter correlations:
1 -0.5254
-0.5254 1
global correlation coefficients :
MnGlobalCorrelationCoeff:
0.5254
0.5254
Warning: divide by zero encountered in double_scalars
TS = -1464.78148767
computing UL for GRB
0 275.92401561 -1.22074436604e-06 0.000276967237302
1 276.60151561 0.000329595973767 0.000277648060094
2 277.27901561 0.00123738017862 0.000278328891069
3 277.95651561 0.00271863386479 0.000279009726482
4 278.63401561 0.00476989198307 0.000279690566353
5 279.31151561 0.00738771659246 0.000280371410677
6 279.98901561 0.010568696605 0.000281052259453
7 280.66651561 0.0143094475186 0.000281733112676
8 281.34401561 0.018606611155 0.000282413970345
9 282.02151561 0.0234568553977 0.000283094832455
10 282.69901561 0.028856873936 0.000283775699005
11 283.37651561 0.0348033860106 0.00028445656999
12 284.05401561 0.0412931361619 0.000285137445409
13 284.73151561 0.0483228939826 0.000285818325258
14 285.40901561 0.0558894538734 0.000286499209534
15 286.08651561 0.063989634799 0.000287180098234
16 286.76401561 0.0726202800508 0.000287860991355
17 287.44151561 0.0817782570093 0.000288541888894
18 288.11901561 0.0914604569115 0.000289222790849
19 288.79651561 0.101663794619 0.000289903697215
20 289.47401561 0.112385208392 0.000290584607991
21 290.15151561 0.123621659662 0.000291265523173
22 290.82901561 0.135370132808 0.000291946442758
23 291.50651561 0.147627634941 0.000292627366743
24 292.18401561 0.160391195683 0.000293308295125
25 292.86151561 0.173657866952 0.000293989227902
26 293.53901561 0.187424722751 0.00029467016507
27 294.21651561 0.201688858957 0.000295351106626
28 294.89401561 0.216447393116 0.000296032052567
29 295.57151561 0.231697464233 0.00029671300289
30 296.24901561 0.247436232575 0.000297393957593
31 296.92651561 0.263660879467 0.000298074916671
32 297.60401561 0.280368607093 0.000298755880124
33 298.28151561 0.297556638306 0.000299436847946
34 298.95901561 0.315222216427 0.000300117820135
35 299.63651561 0.333362605059 0.000300798796689
36 300.31401561 0.351975087896 0.000301479777605
37 300.99151561 0.371056968537 0.000302160762878
38 301.66901561 0.390605570299 0.000302841752507
39 302.34651561 0.410618236037 0.000303522746488
40 303.02401561 0.431092327961 0.000304203744819
41 303.70151561 0.452025227458 0.000304884747496
42 304.37901561 0.473414334918 0.000305565754517
43 305.05651561 0.495257069556 0.000306246765878
44 305.73401561 0.517550869241 0.000306927781576
45 306.41151561 0.540293190324 0.000307608801609
46 307.08901561 0.56348150747 0.000308289825974
47 307.76651561 0.587113313493 0.000308970854667
48 308.44401561 0.611186119188 0.000309651887686
49 309.12151561 0.635697453168 0.000310332925028
50 309.79901561 0.660644861706 0.00031101396669
51 310.47651561 0.68602590857 0.000311695012668
52 311.15401561 0.711838174871 0.00031237606296
53 311.83151561 0.738079258902 0.000313057117564
54 312.50901561 0.764746775986 0.000313738176475
55 313.18651561 0.791838358323 0.000314419239691
56 313.86401561 0.819351654839 0.00031510030721
57 314.54151561 0.847284331036 0.000315781379027
58 315.21901561 0.875634068847 0.000316462455142
59 315.89651561 0.904398566483 0.000317143535549
60 316.57401561 0.933575538298 0.000317824620247
61 317.25151561 0.963162714638 0.000318505709232
62 317.92901561 0.993157841704 0.000319186802503
63 318.60651561 1.02355868141 0.000319867900055
64 319.28401561 1.05436301125 0.000320549001885
65 319.96151561 1.08556862414 0.000321230107992
66 320.63901561 1.11717332832 0.000321911218372
67 321.31651561 1.1491749472 0.000322592333022
68 321.99401561 1.1815713192 0.000323273451939
69 322.67151561 1.21436029767 0.000323954575121
70 323.34901561 1.24753975072 0.000324635702564
71 324.02651561 1.28110756113 0.000325316834266
72 324.70401561 1.31506162618 0.000325997970224
73 325.38151561 1.34939985754 0.000326679110435
74 326.05901561 1.38412018117 0.000327360254896
0.000326300048905 325.004237987
Closing root file...
| bash-2.05b$ ls DATA/GRBOUT/080916009/ |
080916009.root 080916009_LAT_ROI.fits 080916009_LAT_SEL.fits 080916009_LAT_lc.fits 080916009_LAT_spectrum.fits 080916009_LAT.RSP 080916009_LAT_ROI.root 080916009_LAT_expCube.fits 080916009_LAT_lc.gif LAT_fit.dat 080916009_LAT.pha 080916009_LAT_ROI_archive_0.fits 080916009_LAT_expMap.fits 080916009_LAT_map.fits LAT_fit.xcm
bash-2.05b$
| Owned by: Nicola Omodei |
| Last updated by: Chuck Patterson 11/30/2009 |
