|
1
|
********************************************
|
|
2
|
* GLAST Science Analysis Software Log File *
|
|
3
|
* ---------------------------------------- *
|
|
4
|
* v2r9p0 (28/03/2011) *
|
|
5
|
********************************************
|
|
6
|
<sys> Linux kepler.cesr.fr 2.6.31.5 #2 SMP Tue Dec 21 15:34:16 CET 2010 x86_64 x86_64 x86_64 GNU/Linux
|
|
7
|
<sys> System runs Mosix
|
|
8
|
<sys> ST version : "09-32-05"
|
|
9
|
<sys> Add custom IRFs "P6_V11_C1_DIFFUSE,P6_V11_C2_DIFFUSE,P6_V11_DIFFUSE,P6_V11_F1_DIFFUSE,P6_V11_F2_DIFFUSE,P6_V11_I1_DIFFUSE,P6_V11_I2_DIFFUSE,P6_V3_B1_DIFFUSE,P6_V3_B2_DIFFUSE,P6_V3_C1_DIFFUSE,P6_V3_C2_DIFFUSE,P6_V3_F1_DIFFUSE,P6_V3_F2_DIFFUSE,P6_V3_I1_DIFFUSE,P6_V3_I2_DIFFUSE,P7_V6_MEV" to list of accessible IRFs
|
|
10
|
<par> Default IRF: "P6_V11_DIFFUSE"
|
|
11
|
<par> Set default IRF: "P7REP_SOURCE_V15"
|
|
12
|
<cmd> export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtselect rad="15.0" zmax="180.0" evtable="EVENTS" emax="20000.0" outfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/ft1.fits" tmax="0.0" evclsmax="10" evclass=INDEF tmin="0.0" convtype="-1" chatter="4" phasemax="1.0" emin="200.0" gui="no" clobber="yes" evclsmin="0" phasemin="0.0" ra="83.6331" debug="no" dec="22.0145" infile="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ft1.fits" mode="ql"
|
|
13
|
<===> This is gtselect version ScienceTools-09-32-05
|
|
14
|
<===> Applying filter string: ((EVENT_CLASS/4)%2 == 1) && angsep(RA,DEC,83.6331,22.0145) < 15 && 200 <= ENERGY && ENERGY <= 20000 && 0 <= ZENITH_ANGLE && ZENITH_ANGLE <= 100 && gtifilter()
|
|
15
|
<===> Done.
|
|
16
|
<fin> gtselect ok.
|
|
17
|
<cmd> export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtbin sctable="SC_DATA" nypix="200" binsz="0.1" coordsys="GAL" axisrot="0.0" evtable="EVENTS" emax="20000.0" dtime="0.0" snratio="1.0" mode="ql" evfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/ft1.fits" lcemin="30.0" decfield="DEC" ebinfile="NONE" lcemax="200000.0" enumbins="20" scfile="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ft2.fits" tbinalg="LIN" ebinalg="LOG" xref="184.557460101" algorithm="CCUBE" outfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/cntmap.fits" yref="-5.78434644902" tbinfile="NONE" gui="no" emin="200.0" nxpix="200" clobber="yes" efield="ENERGY" debug="no" proj="CAR" tstart="0.0" chatter="4" tstop="0.0" tfield="TIME" rafield="RA" denergy="100.0"
|
|
18
|
<===> This is gtbin version ScienceTools-09-32-05
|
|
19
|
<fin> gtbin ok.
|
|
20
|
<cmd> export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtexpcube2 irfs="P7REP_SOURCE_V15" thmin="0.0" emax="20000.0" nypix="360" binsz="0.5" table="EXPOSURE" enumbins="20" xref="0.0" thmax="180.0" nxpix="720" proj="CAR" coordsys="GAL" bincalc="EDGE" infile="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ltcube.fits" outfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/binned_expmap.fits" chatter="4" axisrot="0.0" emin="200.0" yref="0.0" clobber="yes" cmap="none" mode="ql" debug="no" ignorephi="no"
|
|
21
|
<===> This is gtexpcube2 version ScienceTools-09-32-05
|
|
22
|
<===> Computing binned exposure map....................!
|
|
23
|
<fin> gtexpcube2 ok.
|
|
24
|
<cmd> export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e -m4000 gtsrcmaps irfs="P7REP_SOURCE_V15" sctable="SC_DATA" psfcorr="yes" convol="yes" srcmdl="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/source_model.xml" emapbnds="yes" expcube="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ltcube.fits" resample="yes" debug="no" ptsrc="no" minbinsz="0.1" clobber="yes" bexpmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/binned_expmap.fits" outfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits" cmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/cntmap.fits" mode="ql" chatter="4" rfactor="2" gui="no" scfile="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ft2.fits"
|
|
25
|
<===> This is gtsrcmaps version ScienceTools-09-32-05
|
|
26
|
<===> IRF selection, P7REP_SOURCE_V15, and DSS keywords in /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/cntmap.fits[]
|
|
27
|
<===> are inconsistent. DSS keywords give
|
|
28
|
<===> Creating source named Extragal_diffuse
|
|
29
|
<===> Creating source named Galactic_diffuse
|
|
30
|
<===> Creating source named Crab
|
|
31
|
<===> Generating SourceMap for Extragal_diffuse
|
|
32
|
<===> resampling factor: 2
|
|
33
|
<===> mapsize: 960
|
|
34
|
<===> ....................!
|
|
35
|
<===> Generating SourceMap for Galactic_diffuseMapBase::readFitsFile: creating WcsMap2 object
|
|
36
|
<===>
|
|
37
|
<===> resampling factor: 2
|
|
38
|
<===> mapsize: 960
|
|
39
|
<===> ....................MapBased::deleteMap: /project-data/glast/skymaps/galprop/gll_iem_v05.fit
|
|
40
|
<===> SourceMap::makeDiffuseSource: called mapBaseObj->deleteMap()
|
|
41
|
<===> !
|
|
42
|
<===> Extragal_diffuse
|
|
43
|
<===> appending map for Extragal_diffuse
|
|
44
|
<===> Galactic_diffuse
|
|
45
|
<===> appending map for Galactic_diffuse
|
|
46
|
<fin> gtsrcmaps ok.
|
|
47
|
<cmd> cd /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtlike irfs="P7REP_SOURCE_V15" psfcorr="yes" expcube="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ltcube.fits" results="results.dat" evfile="none" expmap="none" statistic="BINNED" ftol="0.001" plot="no" tsmin="yes" gui="no" save="yes" specfile="counts_spectra.fits" sctable="SC_DATA" optimizer="NEWMINUIT" evtable="EVENTS" refit="no" sfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml" chatter="4" scfile="none" srcmdl="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/source_model.xml" toltype="ABS" clobber="yes" bexpmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/binned_expmap.fits" cmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits" mode="ql" debug="no" check_fit="yes"
|
|
48
|
<===> This is gtlike version ScienceTools-09-32-05
|
|
49
|
<===> IRF selection, P7REP_SOURCE_V15, and DSS keywords in /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits[]
|
|
50
|
<===> are inconsistent. DSS keywords give
|
|
51
|
<===> Creating source named Extragal_diffuse
|
|
52
|
<===> Creating source named Galactic_diffuse
|
|
53
|
<===> Creating source named Crab
|
|
54
|
<===> Generating SourceMap for Crab....................!
|
|
55
|
<===> 0 -26571.6 3755.8
|
|
56
|
<===> 1 -106401 99490.6
|
|
57
|
<===> 2 -1.35496e+06 1.36035e+06
|
|
58
|
<===> 3 -2.68207e+06 2.69076e+06
|
|
59
|
<===> 4 -69269.2 59681.2
|
|
60
|
<===> 5 -26469.1 5161.99
|
|
61
|
<===> 6 -109358 102207
|
|
62
|
<===> 7 -27402.8 3574.39
|
|
63
|
<===> 8 -27202.2 5309.07
|
|
64
|
<===> 9 -26854.7 3685.62
|
|
65
|
<===> 10 -26703.4 4413.85
|
|
66
|
<===> 11 -26531.8 4840.61
|
|
67
|
<===> 12 -26304.7 5089.48
|
|
68
|
<===> 13 -58500.3 47416.5
|
|
69
|
<===> 14 -27048.6 5496.64
|
|
70
|
<===> 15 -26269 5063.75
|
|
71
|
<===> 16 -26029.6 5159.31
|
|
72
|
<===> 17 -83004.7 74357.9
|
|
73
|
<===> 18 -26026.1 5056.14
|
|
74
|
<===> 19 -25995.8 4890.55
|
|
75
|
<===> 20 -25990.3 4788.24
|
|
76
|
<===> 21 -25989.8 4805.36
|
|
77
|
<===> 22 -25980.4 4839.41
|
|
78
|
<===> 23 -25980.1 4835.59
|
|
79
|
<===> 24 -25979.6 4836.36
|
|
80
|
<===> 25 -25979.6 4835.7
|
|
81
|
<===> 26 -25979.6 4831.82
|
|
82
|
<===> 27 -25979.6 4830.82
|
|
83
|
<===> 28 -25979.6 4830.82
|
|
84
|
<===> 29 -25979.6 4833.61
|
|
85
|
<===> 30 -25979.6 4828.03
|
|
86
|
<===> 31 -25984.9 4892.11
|
|
87
|
<===> 32 -25984.3 4776.12
|
|
88
|
<===> 33 -25979.6 4836.64
|
|
89
|
<===> 34 -25979.6 4825.06
|
|
90
|
<===> 35 -25979.6 4832.85
|
|
91
|
<===> 36 -25979.6 4828.79
|
|
92
|
<===> 37 -25979.8 4873.86
|
|
93
|
<===> 38 -25979.8 4789.64
|
|
94
|
<===> 39 -25979.6 4837.56
|
|
95
|
<===> 40 -25979.6 4824.12
|
|
96
|
<===> 41 -26001.6 5287.57
|
|
97
|
<===> 42 -26002.2 4400.21
|
|
98
|
<===> 43 -25979.8 4875.32
|
|
99
|
<===> 44 -25979.8 4786.58
|
|
100
|
<===> 45 -25979.6 4838.21
|
|
101
|
<===> 46 -25979.6 4823.43
|
|
102
|
<===> 47 -25979.6 4833.61
|
|
103
|
<===> 48 -25979.6 4828.03
|
|
104
|
<===> 49 -25979.6 4832.85
|
|
105
|
<===> 50 -25979.6 4828.79
|
|
106
|
<===> 51 -25979.6 4837.56
|
|
107
|
<===> 52 -25979.6 4824.12
|
|
108
|
<===> 53 -25979.6 4838.21
|
|
109
|
<===> 54 -25979.6 4823.43
|
|
110
|
<===> 55 -25979.6 4831.38
|
|
111
|
<===> 56 -25979.6 4830.26
|
|
112
|
<===> 57 -25979.6 4831.22
|
|
113
|
<===> 58 -25979.6 4830.41
|
|
114
|
<===> 59 -25979.6 4832.16
|
|
115
|
<===> 60 -25979.6 4829.47
|
|
116
|
<===> 61 -25979.6 4831.09
|
|
117
|
<===> 62 -25979.6 4830.55
|
|
118
|
<===> 63 -25979.6 4832.29
|
|
119
|
<===> 64 -25979.6 4829.34
|
|
120
|
<===> 65 -25979.6 4835.66
|
|
121
|
<===> 66 -25979.6 4840.36
|
|
122
|
<===> 67 -25979.6 4841.01
|
|
123
|
<===> 68 -25979.6 4839.6
|
|
124
|
<===> 69 -25979.6 4840.25
|
|
125
|
<===> 70 -25979.6 4844.96
|
|
126
|
<===>
|
|
127
|
<===> Minuit did successfully converge.
|
|
128
|
<===> # of function calls: 71
|
|
129
|
<===> minimum function Value: 25979.56816776
|
|
130
|
<===> minimum edm: 0.0007661844713016
|
|
131
|
<===> minimum internal state vector: LAVector parameters:
|
|
132
|
<===> -1.261489214475
|
|
133
|
<===> -0.4422987551761
|
|
134
|
<===> -1.496769208144
|
|
135
|
<===> -1.504080741126
|
|
136
|
<===>
|
|
137
|
<===> minimum internal covariance matrix: LASymMatrix parameters:
|
|
138
|
<===> 0.00026769123 -0.00019723322 -2.3849844e-05 4.841396e-07
|
|
139
|
<===> -0.00019723322 0.00027009987 2.0401873e-05 -1.2160875e-06
|
|
140
|
<===> -2.3849844e-05 2.0401873e-05 0.00018865593 -2.0456189e-05
|
|
141
|
<===> 4.841396e-07 -1.2160875e-06 -2.0456189e-05 2.9322636e-06
|
|
142
|
<===>
|
|
143
|
<===>
|
|
144
|
<===> # ext. || Name || type || Value || Error +/-
|
|
145
|
<===>
|
|
146
|
<===> 0 ||0_Integral || limited || 23.72764377996 ||1.7607803806
|
|
147
|
<===> 1 || 1_Index || limited || -2.140090702065 ||0.0525128279225
|
|
148
|
<===> 2 ||2_Normalization || limited || 1.369378051602 ||0.3591513710273
|
|
149
|
<===> 3 || 3_Value || limited || 1.112329672179 ||0.04036096008417
|
|
150
|
<===>
|
|
151
|
<===>
|
|
152
|
<===>
|
|
153
|
<===>
|
|
154
|
<===> 71 -25979.6 4830.82
|
|
155
|
<===> 72 -25979.6 4831.38
|
|
156
|
<===> 73 -25979.6 4830.26
|
|
157
|
<===> 74 -25979.6 4831.22
|
|
158
|
<===> 75 -25979.6 4830.41
|
|
159
|
<===> 76 -25979.6 4831.09
|
|
160
|
<===> 77 -25979.6 4830.55
|
|
161
|
<===> 78 -25979.6 4832.29
|
|
162
|
<===> 79 -25979.6 4829.34
|
|
163
|
<===> 80 -25979.6 4830.93
|
|
164
|
<===> 81 -25979.6 4830.7
|
|
165
|
<===> 82 -25979.6 4830.9
|
|
166
|
<===> 83 -25979.6 4830.73
|
|
167
|
<===> 84 -25979.6 4830.87
|
|
168
|
<===> 85 -25979.6 4830.76
|
|
169
|
<===> 86 -25979.6 4831.11
|
|
170
|
<===> 87 -25979.6 4830.52
|
|
171
|
<===> 88 -25979.6 4831.78
|
|
172
|
<===> 89 -25979.6 4831.64
|
|
173
|
<===> 90 -25979.6 4832.85
|
|
174
|
<===> 91 -25979.6 4831.49
|
|
175
|
<===> 92 -25979.6 4832.7
|
|
176
|
<===> 93 -25979.6 4832.56
|
|
177
|
<===>
|
|
178
|
<===> # of function calls: 94
|
|
179
|
<===> function Value: 25979.56816776
|
|
180
|
<===> expected distance to the Minimum (edm): 0.0007635878349732
|
|
181
|
<===> external parameters:
|
|
182
|
<===> # ext. || Name || type || Value || Error +/-
|
|
183
|
<===>
|
|
184
|
<===> 0 ||0_Integral || limited || 23.72764377996 ||1.752643597136
|
|
185
|
<===> 1 || 1_Index || limited || -2.140090702065 ||0.0522717745185
|
|
186
|
<===> 2 ||2_Normalization || limited || 1.369378051602 ||0.3582834903128
|
|
187
|
<===> 3 || 3_Value || limited || 1.112329672179 ||0.04026935761572
|
|
188
|
<===>
|
|
189
|
<===>
|
|
190
|
<===>
|
|
191
|
<===> covariance matrix:
|
|
192
|
<===> MnUserCovariance:
|
|
193
|
<===>
|
|
194
|
<===> 3.0719 -0.0669339 -0.0656889 0.00112112
|
|
195
|
<===> -0.0669339 0.00273246 0.00166842 -8.93026e-05
|
|
196
|
<===> -0.0656889 0.00166842 0.128371 -0.0125396
|
|
197
|
<===> 0.00112112 -8.93026e-05 -0.0125396 0.00162162
|
|
198
|
<===>
|
|
199
|
<===> MnUserCovariance Parameter correlations:
|
|
200
|
<===>
|
|
201
|
<===> 1 -0.730577 -0.104606 0.0158845
|
|
202
|
<===> -0.730577 1 0.0890832 -0.0424241
|
|
203
|
<===> -0.104606 0.0890832 1 -0.869112
|
|
204
|
<===> 0.0158845 -0.0424241 -0.869112 1
|
|
205
|
<===>
|
|
206
|
<===> global correlation coefficients :
|
|
207
|
<===> MnGlobalCorrelationCoeff:
|
|
208
|
<===>
|
|
209
|
<===> 0.738524
|
|
210
|
<===> 0.7318
|
|
211
|
<===> 0.874088
|
|
212
|
<===> 0.872858
|
|
213
|
<===>
|
|
214
|
<===>
|
|
215
|
<===> 94 -25979.6 4830.82
|
|
216
|
<===> .Computing TS values for each source (3 total)
|
|
217
|
<===> 95 -26730.6 4235.41
|
|
218
|
<===> 96 -26730.6 4235.41
|
|
219
|
<===> 97 -875108 878386
|
|
220
|
<===> 98 -3.23481e+06 3.2444e+06
|
|
221
|
<===> 99 -46628.8 34219.2
|
|
222
|
<===> 100 -26695.1 5029.12
|
|
223
|
<===> 101 -27357.7 7347.66
|
|
224
|
<===> 102 -26693.7 5020.25
|
|
225
|
<===> 103 -26689.8 4839.86
|
|
226
|
<===> 104 -26689.8 4827.9
|
|
227
|
<===> 105 -26689.7 4833.87
|
|
228
|
<===> 106 -26689.7 4835.33
|
|
229
|
<===> 107 -26689.7 4832.1
|
|
230
|
<===> 108 -26689.7 4832.1
|
|
231
|
<===> 109 -26689.9 4873.24
|
|
232
|
<===> 110 -26689.9 4792.14
|
|
233
|
<===> 111 -26689.7 4839.27
|
|
234
|
<===> 112 -26689.7 4824.96
|
|
235
|
<===> 113 -26708.8 5272.6
|
|
236
|
<===> 114 -26709.1 4413.84
|
|
237
|
<===> 115 -26689.9 4875.15
|
|
238
|
<===> 116 -26689.9 4789.27
|
|
239
|
<===> 117 -26689.7 4839.92
|
|
240
|
<===> 118 -26689.7 4824.28
|
|
241
|
<===> 119 -26689.7 4839.27
|
|
242
|
<===> 120 -26689.7 4824.96
|
|
243
|
<===> 121 -26689.7 4839.92
|
|
244
|
<===> 122 -26689.7 4824.28
|
|
245
|
<===> 123 -26689.7 4833.53
|
|
246
|
<===> 124 -26689.7 4830.67
|
|
247
|
<===> 125 -26689.7 4832.38
|
|
248
|
<===> 126 -26689.7 4831.81
|
|
249
|
<===> 127 -26689.7 4833.66
|
|
250
|
<===> 128 -26689.7 4830.53
|
|
251
|
<===> 129 -26689.8 4847.1
|
|
252
|
<===> 130 -26689.7 4832.1
|
|
253
|
<===> Generating SourceMap for Crab....................!
|
|
254
|
<===> ..!
|
|
255
|
<===>
|
|
256
|
<===> Photon fluxes are computed for the energy range 200 to 20000 MeV
|
|
257
|
<===>
|
|
258
|
<===> Crab:
|
|
259
|
<===> Integral: 23.7276 +/- 1.75268
|
|
260
|
<===> Index: -2.14009 +/- 0.0522729
|
|
261
|
<===> LowerLimit: 100
|
|
262
|
<===> UpperLimit: 500000
|
|
263
|
<===> TS value: 1420.34
|
|
264
|
<===> Flux: 1.07127e-06 +/- 5.77166e-08 photons/cm^2/s
|
|
265
|
<===>
|
|
266
|
<===> Extragal_diffuse:
|
|
267
|
<===> Normalization: 1.36938 +/- 0.358289
|
|
268
|
<===> Flux: 9.52653e-05 +/- 2.49164e-05 photons/cm^2/s
|
|
269
|
<===>
|
|
270
|
<===> Galactic_diffuse:
|
|
271
|
<===> Value: 1.11233 +/- 0.0402694
|
|
272
|
<===> MapBase::readFitsFile: creating WcsMap2 object
|
|
273
|
<===> Flux: 0.000332084 +/- 1.20219e-05 photons/cm^2/s
|
|
274
|
<===> WARNING: Fit may be bad in range [3169.79, 3990.52] (MeV)
|
|
275
|
<===>
|
|
276
|
<===>
|
|
277
|
<===> Total number of observed counts: 4832
|
|
278
|
<===> Total number of model events: 4830.82
|
|
279
|
<===> 131 -25979.56817 4830.81618
|
|
280
|
<===>
|
|
281
|
<===> -log(Likelihood): 25979.56817
|
|
282
|
<===>
|
|
283
|
<===> Writing fitted model to /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml
|
|
284
|
<===> Elapsed CPU time: 51.72
|
|
285
|
<fin> gtlike ok.
|
|
286
|
<log> gtlike refit requested! Restart.
|
|
287
|
<cmd> cd /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtlike irfs="P7REP_SOURCE_V15" psfcorr="yes" expcube="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ltcube.fits" results="results.dat" evfile="none" expmap="none" statistic="BINNED" ftol="0.001" plot="no" tsmin="yes" gui="no" save="yes" specfile="counts_spectra.fits" sctable="SC_DATA" optimizer="NEWMINUIT" evtable="EVENTS" refit="no" sfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml" chatter="4" scfile="none" srcmdl="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml" toltype="ABS" clobber="yes" bexpmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/binned_expmap.fits" cmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits" mode="ql" debug="no" check_fit="yes"
|
|
288
|
<===> This is gtlike version ScienceTools-09-32-05
|
|
289
|
<===> IRF selection, P7REP_SOURCE_V15, and DSS keywords in /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits[]
|
|
290
|
<===> are inconsistent. DSS keywords give
|
|
291
|
<===> Creating source named Crab
|
|
292
|
<===> Creating source named Extragal_diffuse
|
|
293
|
<===> Creating source named Galactic_diffuse
|
|
294
|
<===> Generating SourceMap for Crab....................!
|
|
295
|
<===> 0 -25979.6 4830.82
|
|
296
|
<===> 1 -27743.9 9901.31
|
|
297
|
<===> 2 -25984.2 5039.25
|
|
298
|
<===> 3 -25979.6 4841.11
|
|
299
|
<===> 4 -25979.6 4832.47
|
|
300
|
<===> 5 -26184.9 4713.47
|
|
301
|
<===> 6 -25980 4822.32
|
|
302
|
<===> 7 -25979.6 4831.95
|
|
303
|
<===> 8 -25979.6 4832.23
|
|
304
|
<===> 9 -25979.6 4831.54
|
|
305
|
<===> 10 -25979.6 4832.12
|
|
306
|
<===> 11 -25979.6 4832.12
|
|
307
|
<===> 12 -25979.6 4834.92
|
|
308
|
<===> 13 -25979.6 4829.33
|
|
309
|
<===> 14 -25985 4893.46
|
|
310
|
<===> 15 -25984.3 4777.38
|
|
311
|
<===> 16 -25979.6 4837.95
|
|
312
|
<===> 17 -25979.6 4826.36
|
|
313
|
<===> 18 -25979.6 4834.16
|
|
314
|
<===> 19 -25979.6 4830.09
|
|
315
|
<===> 20 -25979.8 4875.23
|
|
316
|
<===> 21 -25979.8 4790.9
|
|
317
|
<===> 22 -25979.6 4838.87
|
|
318
|
<===> 23 -25979.6 4825.42
|
|
319
|
<===> 24 -26001.6 5288.59
|
|
320
|
<===> 25 -26002.1 4401.75
|
|
321
|
<===> 26 -25979.8 4876.59
|
|
322
|
<===> 27 -25979.8 4787.91
|
|
323
|
<===> 28 -25979.6 4839.52
|
|
324
|
<===> 29 -25979.6 4824.73
|
|
325
|
<===> 30 -25979.6 4834.92
|
|
326
|
<===> 31 -25979.6 4829.33
|
|
327
|
<===> 32 -25979.6 4834.16
|
|
328
|
<===> 33 -25979.6 4830.09
|
|
329
|
<===> 34 -25979.6 4838.87
|
|
330
|
<===> 35 -25979.6 4825.42
|
|
331
|
<===> 36 -25979.6 4839.52
|
|
332
|
<===> 37 -25979.6 4824.73
|
|
333
|
<===> 38 -25979.6 4832.68
|
|
334
|
<===> 39 -25979.6 4831.56
|
|
335
|
<===> 40 -25979.6 4832.53
|
|
336
|
<===> 41 -25979.6 4831.72
|
|
337
|
<===> 42 -25979.6 4832.2
|
|
338
|
<===> 43 -25979.6 4832.04
|
|
339
|
<===> 44 -25979.6 4833.47
|
|
340
|
<===> 45 -25979.6 4830.78
|
|
341
|
<===> 46 -25979.6 4832.39
|
|
342
|
<===> 47 -25979.6 4831.85
|
|
343
|
<===> 48 -25979.6 4833.6
|
|
344
|
<===> 49 -25979.6 4830.64
|
|
345
|
<===> 50 -25979.6 4832.42
|
|
346
|
<===> 51 -25979.6 4831.83
|
|
347
|
<===> 52 -25979.6 4836.97
|
|
348
|
<===> 53 -25979.6 4841.67
|
|
349
|
<===> 54 -25979.6 4842.32
|
|
350
|
<===> 55 -25979.6 4840.91
|
|
351
|
<===> 56 -25979.6 4841.56
|
|
352
|
<===> 57 -25979.6 4846.27
|
|
353
|
<===>
|
|
354
|
<===> Minuit did successfully converge.
|
|
355
|
<===> # of function calls: 58
|
|
356
|
<===> minimum function Value: 25979.56744451
|
|
357
|
<===> minimum edm: 4.02477117386e-05
|
|
358
|
<===> minimum internal state vector: LAVector parameters:
|
|
359
|
<===> -1.261493909278
|
|
360
|
<===> -0.4420806567806
|
|
361
|
<===> -1.496976627901
|
|
362
|
<===> -1.504049894637
|
|
363
|
<===>
|
|
364
|
<===> minimum internal covariance matrix: LASymMatrix parameters:
|
|
365
|
<===> 0.00026769507 -0.00019719371 -2.389209e-05 4.797087e-07
|
|
366
|
<===> -0.00019719371 0.00026981508 2.0435687e-05 -1.2105976e-06
|
|
367
|
<===> -2.389209e-05 2.0435687e-05 0.00018948584 -2.0484263e-05
|
|
368
|
<===> 4.797087e-07 -1.2105976e-06 -2.0484263e-05 2.9282314e-06
|
|
369
|
<===>
|
|
370
|
<===>
|
|
371
|
<===> # ext. || Name || type || Value || Error +/-
|
|
372
|
<===>
|
|
373
|
<===> 0 ||0_Integral || limited || 23.72692923934 ||1.760767113393
|
|
374
|
<===> 1 || 1_Index || limited || -2.139105097084 ||0.05249055711983
|
|
375
|
<===> 2 ||2_Normalization || limited || 1.361718444533 ||0.3589337411366
|
|
376
|
<===> 3 || 3_Value || limited || 1.113358117175 ||0.04035182067449
|
|
377
|
<===>
|
|
378
|
<===>
|
|
379
|
<===>
|
|
380
|
<===>
|
|
381
|
<===> 58 -25979.6 4832.12
|
|
382
|
<===> 59 -25979.6 4832.68
|
|
383
|
<===> 60 -25979.6 4831.56
|
|
384
|
<===> 61 -25979.6 4832.2
|
|
385
|
<===> 62 -25979.6 4832.04
|
|
386
|
<===> 63 -25979.6 4832.39
|
|
387
|
<===> 64 -25979.6 4831.85
|
|
388
|
<===> 65 -25979.6 4832.42
|
|
389
|
<===> 66 -25979.6 4831.83
|
|
390
|
<===> 67 -25979.6 4832.23
|
|
391
|
<===> 68 -25979.6 4832.01
|
|
392
|
<===> 69 -25979.6 4832.14
|
|
393
|
<===> 70 -25979.6 4832.11
|
|
394
|
<===> 71 -25979.6 4832.18
|
|
395
|
<===> 72 -25979.6 4832.07
|
|
396
|
<===> 73 -25979.6 4832.18
|
|
397
|
<===> 74 -25979.6 4832.06
|
|
398
|
<===> 75 -25979.6 4832.76
|
|
399
|
<===> 76 -25979.6 4832.95
|
|
400
|
<===> 77 -25979.6 4832.98
|
|
401
|
<===> 78 -25979.6 4832.47
|
|
402
|
<===> 79 -25979.6 4832.5
|
|
403
|
<===> 80 -25979.6 4832.69
|
|
404
|
<===>
|
|
405
|
<===> # of function calls: 81
|
|
406
|
<===> function Value: 25979.56744451
|
|
407
|
<===> expected distance to the Minimum (edm): 3.983255451288e-05
|
|
408
|
<===> external parameters:
|
|
409
|
<===> # ext. || Name || type || Value || Error +/-
|
|
410
|
<===>
|
|
411
|
<===> 0 ||0_Integral || limited || 23.72692923934 ||1.751595617067
|
|
412
|
<===> 1 || 1_Index || limited || -2.139105097084 ||0.05221887255349
|
|
413
|
<===> 2 ||2_Normalization || limited || 1.361718444533 ||0.3582557203799
|
|
414
|
<===> 3 || 3_Value || limited || 1.113358117175 ||0.04028198719927
|
|
415
|
<===>
|
|
416
|
<===>
|
|
417
|
<===>
|
|
418
|
<===> covariance matrix:
|
|
419
|
<===> MnUserCovariance:
|
|
420
|
<===>
|
|
421
|
<===> 3.06822 -0.0668135 -0.0655504 0.00111039
|
|
422
|
<===> -0.0668135 0.00272693 0.0016649 -8.89968e-05
|
|
423
|
<===> -0.0655504 0.0016649 0.128351 -0.0125431
|
|
424
|
<===> 0.00111039 -8.89968e-05 -0.0125431 0.00162264
|
|
425
|
<===>
|
|
426
|
<===> MnUserCovariance Parameter correlations:
|
|
427
|
<===>
|
|
428
|
<===> 1 -0.730439 -0.104456 0.015737
|
|
429
|
<===> -0.730439 1 0.0889918 -0.0423084
|
|
430
|
<===> -0.104456 0.0889918 1 -0.869149
|
|
431
|
<===> 0.015737 -0.0423084 -0.869149 1
|
|
432
|
<===>
|
|
433
|
<===> global correlation coefficients :
|
|
434
|
<===> MnGlobalCorrelationCoeff:
|
|
435
|
<===>
|
|
436
|
<===> 0.738384
|
|
437
|
<===> 0.731659
|
|
438
|
<===> 0.87412
|
|
439
|
<===> 0.872894
|
|
440
|
<===>
|
|
441
|
<===>
|
|
442
|
<===> 81 -25979.6 4832.12
|
|
443
|
<===> .Computing TS values for each source (3 total)
|
|
444
|
<===> 82 -26730.5 4236.23
|
|
445
|
<===> 83 -26730.5 4236.23
|
|
446
|
<===> 84 -790565 793347
|
|
447
|
<===> 85 -3.33659e+06 3.34632e+06
|
|
448
|
<===> 86 -49467.9 37500.6
|
|
449
|
<===> 87 -26698 5089.92
|
|
450
|
<===> 88 -27506.2 7707.5
|
|
451
|
<===> 89 -26696.2 5078.52
|
|
452
|
<===> 90 -26689.9 4838.23
|
|
453
|
<===> 91 -26689.8 4833.59
|
|
454
|
<===> 92 -26689.7 4832.26
|
|
455
|
<===> 93 -26689.7 4832.26
|
|
456
|
<===> 94 -26690 4873.39
|
|
457
|
<===> 95 -26689.9 4792.31
|
|
458
|
<===> 96 -26689.7 4839.44
|
|
459
|
<===> 97 -26689.7 4825.12
|
|
460
|
<===> 98 -26708.8 5272.84
|
|
461
|
<===> 99 -26709.1 4413.94
|
|
462
|
<===> 100 -26689.9 4875.32
|
|
463
|
<===> 101 -26689.9 4789.43
|
|
464
|
<===> 102 -26689.7 4840.09
|
|
465
|
<===> 103 -26689.7 4824.44
|
|
466
|
<===> 104 -26689.7 4839.44
|
|
467
|
<===> 105 -26689.7 4825.12
|
|
468
|
<===> 106 -26689.7 4840.09
|
|
469
|
<===> 107 -26689.7 4824.44
|
|
470
|
<===> 108 -26689.7 4833.69
|
|
471
|
<===> 109 -26689.7 4830.83
|
|
472
|
<===> 110 -26689.7 4833.83
|
|
473
|
<===> 111 -26689.7 4830.7
|
|
474
|
<===> 112 -26689.8 4847.26
|
|
475
|
<===> 113 -26689.7 4832.26
|
|
476
|
<===> Generating SourceMap for Crab....................!
|
|
477
|
<===> ..!
|
|
478
|
<===>
|
|
479
|
<===> Photon fluxes are computed for the energy range 200 to 20000 MeV
|
|
480
|
<===>
|
|
481
|
<===> Crab:
|
|
482
|
<===> Integral: 23.7269 +/- 1.75163
|
|
483
|
<===> Index: -2.13911 +/- 0.05222
|
|
484
|
<===> LowerLimit: 100
|
|
485
|
<===> UpperLimit: 500000
|
|
486
|
<===> TS value: 1420.34
|
|
487
|
<===> Flux: 1.07194e-06 +/- 5.77353e-08 photons/cm^2/s
|
|
488
|
<===>
|
|
489
|
<===> Extragal_diffuse:
|
|
490
|
<===> Normalization: 1.36172 +/- 0.358261
|
|
491
|
<===> Flux: 9.47324e-05 +/- 2.49144e-05 photons/cm^2/s
|
|
492
|
<===>
|
|
493
|
<===> Galactic_diffuse:
|
|
494
|
<===> Value: 1.11336 +/- 0.040282
|
|
495
|
<===> MapBase::readFitsFile: creating WcsMap2 object
|
|
496
|
<===> Flux: 0.000332391 +/- 1.20257e-05 photons/cm^2/s
|
|
497
|
<===> WARNING: Fit may be bad in range [3169.79, 3990.52] (MeV)
|
|
498
|
<===>
|
|
499
|
<===>
|
|
500
|
<===> Total number of observed counts: 4832
|
|
501
|
<===> Total number of model events: 4832.12
|
|
502
|
<===> 114 -25979.56744 4832.121898
|
|
503
|
<===>
|
|
504
|
<===> -log(Likelihood): 25979.56744
|
|
505
|
<===>
|
|
506
|
<===> Writing fitted model to /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml
|
|
507
|
<===> Elapsed CPU time: 59.04
|
|
508
|
<fin> gtlike ok.
|
|
509
|
<cmd> export PFILES=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; export HOME=/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab; mosrun -e gtmodel irfs="P7REP_SOURCE_V15" expcube="/project-data/glast/flight/fssc/p7v6/lightcurve/weekly/bin_00000/source/ltcube.fits" resample="yes" outfile="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/modmap.fits" outtype="CMAP" chatter="4" srcmaps="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits" srcmdl="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/results_like.xml" convol="yes" gui="no" clobber="yes" bexpmap="/users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/binned_expmap.fits" mode="ql" debug="no" rfactor="2"
|
|
510
|
<===> This is gtmodel version ScienceTools-09-32-05
|
|
511
|
<===> IRF selection, P7REP_SOURCE_V15, and DSS keywords in /users-data/knodlseder/glast/analysis/gammalib/modelfit/week_00000_roi10_200_20000_20bin_extdiff_galdiff_Crab/srcmap.fits[]
|
|
512
|
<===> are inconsistent. DSS keywords give
|
|
513
|
<===> Creating source named Crab
|
|
514
|
<===> Creating source named Extragal_diffuse
|
|
515
|
<===> Creating source named Galactic_diffuse
|
|
516
|
<===> Generating SourceMap for Crab....................!
|
|
517
|
<fin> gtmodel ok.
|
|
518
|
<fin> SAS terminated properly.
|