Domenico Tiziani / gammalib

Commit 47a0e1a2c7bfe62e3a9a142a702e0a1e71fa88ef

Authored by Jürgen Knödlseder
1 parent ac040aff

Add test files

Inline Side-by-side
Showing 11 changed files with 841 additions and 0 deletions
inst/lat/test/data/p7repv15/Crab.xml 0 → 100644
  View file @47a0e1a
  1 +<source_library title="source library">
  2 +<source name="iso_source_v05" type="DiffuseSource">
  3 + <spectrum file="$(PACKAGE_SOURCE)/inst/lat/test/data/p7repv15/iso_source_v05.txt" type="FileFunction">
  4 + <parameter free="1" max="10" min="0" name="Normalization" scale="1" value="1.0" />
  5 + </spectrum>
  6 + <spatialModel type="ConstantValue">
  7 + <parameter free="0" max="10" min="0" name="Value" scale="1" value="1" />
  8 + </spatialModel>
  9 + </source>
  10 +<source name="gll_iem_v05" type="DiffuseSource">
  11 + <spectrum type="PowerLaw">
  12 + <parameter free="1" max="100" min="1e-05" name="Prefactor" scale="1" value="1." />
  13 + <parameter free="1" max="1" min="-1" name="Index" scale="1" value="0." />
  14 + <parameter free="0" max="2000" min="50" name="Scale" scale="1" value="100" />
  15 + </spectrum>
  16 + <spatialModel type="MapCubeFunction" file="$(PACKAGE_SOURCE)/test/data/test_cube.fits">
  17 + <parameter free="0" max="1000" min="0.001" name="Normalization" scale="1" value="1" />
  18 + </spatialModel>
  19 + </source>
  20 + <source type="PointSource" name="Crab">
  21 + <spectrum type="PowerLaw2">
  22 + <parameter scale="1e-07" name="Integral" min="1e-05" max="1000.0" value="1.0" free="1"/>
  23 + <parameter scale="1.0" name="Index" min="-5.0" max="-0.1" value="-2.0" free="1"/>
  24 + <parameter scale="1.0" name="LowerLimit" min="10.0" max="1000000.0" value="100.0" free="0"/>
  25 + <parameter scale="1.0" name="UpperLimit" min="10.0" max="1000000.0" value="500000.0" free="0"/>
  26 + </spectrum>
  27 + <spatialModel type="SkyDirFunction">
  28 + <parameter scale="1" name="RA" min="-360" max="360" value="83.6331" free="0"/>
  29 + <parameter scale="1" name="DEC" min="-90" max="90" value="22.0145" free="0"/>
  30 + </spatialModel>
  31 + </source>
  32 +</source_library>
inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt 0 → 100644
  View file @47a0e1a
  1 +LAT example analysis of Crab region
  2 +First analysis with Scince tools 09-32-05, then analysis of the same files with gtlike and ctlike:
  3 +
  4 +ra = 83.6331
  5 +dec = 22.0145
  6 +
  7 +rad = 7.5
  8 +10x10 deg, 0.1 degperbin
  9 +zmax = 100
  10 +emin = 100 MeV
  11 +emax = 10000 MeV
  12 +
  13 +
  14 +tmin = 239557446.6 MET in s
  15 +tmax= 240162199.6 MET in s
  16 +nebins=20
  17 +
  18 +
  19 +
  20 +
  21 +gtlike results:
  22 +
  23 +<?xml version="1.0" standalone="no"?>
  24 +<source_library title="source library">
  25 + <source name="Crab" type="PointSource">
  26 + <spectrum type="PowerLaw2">
  27 + <parameter error="1.510702125" free="1" max="1000" min="1e-05" name="Integral" scale="1e-07" value="25.58700179" />
  28 + <parameter error="0.04594388483" free="1" max="-0.1" min="-5" name="Index" scale="1" value="-2.135443533" />
  29 + <parameter free="0" max="1000000" min="10" name="LowerLimit" scale="1" value="100" />
  30 + <parameter free="0" max="1000000" min="10" name="UpperLimit" scale="1" value="500000" />
  31 + </spectrum>
  32 + <spatialModel type="SkyDirFunction">
  33 + <parameter free="0" max="360" min="-360" name="RA" scale="1" value="83.6331" />
  34 + <parameter free="0" max="90" min="-90" name="DEC" scale="1" value="22.0145" />
  35 + </spatialModel>
  36 + </source>
  37 + <source name="gll_iem_v05" type="DiffuseSource">
  38 + <spectrum type="PowerLaw">
  39 + <parameter error="0.144060828" free="1" max="100" min="1e-05" name="Prefactor" scale="1" value="0.7287081944" />
  40 + <parameter error="0.0716234477" free="1" max="1" min="-1" name="Index" scale="1" value="0.08682762462" />
  41 + <parameter free="0" max="2000" min="50" name="Scale" scale="1" value="100" />
  42 + </spectrum>
  43 + <spatialModel file="/lustre/fs8/group/hess/user/aschulz/lat_data/fssc_files/diffuse/p7rep/gll_iem_v05.fit" type="MapCubeFunction">
  44 + <parameter free="0" max="1000" min="0.001" name="Normalization" scale="1" value="1" />
  45 + </spatialModel>
  46 + </source>
  47 + <source name="iso_source_v05" type="DiffuseSource">
  48 + <spectrum file="/lustre/fs8/group/hess/user/aschulz/lat_data/fssc_files/diffuse/p7rep/iso_source_v05.txt" type="FileFunction">
  49 + <parameter error="0.6855552052" free="1" max="10" min="0" name="Normalization" scale="1" value="2.875818636" />
  50 + </spectrum>
  51 + <spatialModel type="ConstantValue">
  52 + <parameter free="0" max="10" min="0" name="Value" scale="1" value="1" />
  53 + </spatialModel>
  54 + </source>
  55 +</source_library>
  56 +
  57 +
  58 +
  59 +
  60 +ctlike:
  61 +
  62 +<?xml version="1.0" encoding="UTF-8" standalone="no"?>
  63 +<source_library title="source library">
  64 + <source name="iso_source_v05" type="DiffuseSource">
  65 + <spectrum type="FileFunction" file="/Users/schulz/software/gammalib/inst/lat/test/data/p7repv15/iso_source_v05.txt">
  66 + <parameter name="Normalization" value="3.48618" error="0.694252" scale="1" min="0" max="10" free="1" />
  67 + </spectrum>
  68 + <spatialModel type="ConstantValue">
  69 + <parameter name="Value" value="1" scale="1" min="0" max="10" free="0" />
  70 + </spatialModel>
  71 + </source>
  72 + <source name="gll_iem_v05" type="DiffuseSource">
  73 + <spectrum type="PowerLaw">
  74 + <parameter name="Prefactor" value="0.774972" error="0.156895" scale="1" min="1e-05" max="100" free="1" />
  75 + <parameter name="Index" value="0.0695083" error="0.0745843" scale="1" min="-1" max="1" free="1" />
  76 + <parameter name="Scale" value="100" scale="1" min="50" max="2000" free="0" />
  77 + </spectrum>
  78 + <spatialModel type="MapCubeFunction" file="/Users/schulz/software/gammalib/inst/lat/test/data/p7repv15/gll_iem_v05.fit">
  79 + <parameter name="Normalization" value="1" scale="1" min="0.001" max="1000" free="0" />
  80 + </spatialModel>
  81 + </source>
  82 + <source name="Crab" type="PointSource">
  83 + <spectrum type="PowerLaw2">
  84 + <parameter name="Integral" value="19.9229" error="1.19129" scale="1e-07" min="1e-05" max="1000" free="1" />
  85 + <parameter name="Index" value="-2.05072" error="0.0491567" scale="1" min="-5" max="-0.1" free="1" />
  86 + <parameter name="LowerLimit" value="100" scale="1" min="10" max="1e+06" free="0" />
  87 + <parameter name="UpperLimit" value="500000" scale="1" min="10" max="1e+06" free="0" />
  88 + </spectrum>
  89 + <spatialModel type="SkyDirFunction">
  90 + <parameter name="RA" value="83.6331" scale="1" min="-360" max="360" free="0" />
  91 + <parameter name="DEC" value="22.0145" scale="1" min="-90" max="90" free="0" />
  92 + </spatialModel>
  93 + </source>
  94 +</source_library>
  95 +
  96 +
  97 +
  98 +
  99 +
  100 +
  101 +
  102 +
  103 +logfiles:
  104 +
  105 +------------------ running: pybinlike ------------------
  106 +Parameters:
  107 + irfs = P7REP_SOURCE_V15
  108 + expcube = ltcube_239557446.6_240162199.6_scaled.fits
  109 + srcmdl = Crab.xml
  110 + cmap = srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits
  111 + bexpmap = expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits
  112 + optimizer = Minuit
  113 + specfile = binned_specfile_py_Crab_ct.fits
  114 + sfile = binned_fitted_py_Crab_ct.xml
  115 + results = binned_results_py_Crab_ct.dat
  116 + scratchdir = none
  117 +{'irfs': 'P7REP_SOURCE_V15', 'expcube': 'ltcube_239557446.6_240162199.6_scaled.fits', 'optimizer': 'Minuit', 'srcmdl': 'Crab.xml', 'sfile': 'binned_fitted_py_Crab_ct.xml', 'results': 'binned_results_py_Crab_ct.dat', 'bexpmap': 'expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits', 'scratchdir': 'none', 'cmap': 'srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits', 'specfile': 'binned_specfile_py_Crab_ct.fits'}
  118 +{'irfs': 'P7REP_SOURCE_V15', 'expcube': 'ltcube_239557446.6_240162199.6_scaled.fits', 'optimizer': 'Minuit', 'srcmdl': 'Crab.xml', 'sfile': 'binned_fitted_py_Crab_ct.xml', 'results': 'binned_results_py_Crab_ct.dat', 'bexpmap': 'expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits', 'scratchdir': 'none', 'cmap': 'srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits', 'specfile': 'binned_specfile_py_Crab_ct.fits'}
  119 +
  120 +Parameters:
  121 + irfs = P7REP_SOURCE_V15
  122 + expcube = ltcube_239557446.6_240162199.6_scaled.fits
  123 + optimizer = Minuit
  124 + srcmdl = Crab.xml
  125 + sfile = binned_fitted_py_Crab_ct.xml
  126 + results = binned_results_py_Crab_ct.dat
  127 + bexpmap = expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits
  128 + scratchdir = none
  129 + cmap = srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits
  130 + specfile = binned_specfile_py_Crab_ct.fits
  131 +
  132 + Setting up binned likelihood
  133 + Initial model parameters
  134 +
  135 +Source maps: srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits
  136 +Exposure cube: ltcube_239557446.6_240162199.6_scaled.fits
  137 +Exposure map: expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits
  138 +IRFs: P7REP_SOURCE_V15
  139 +Source model file: Crab.xml
  140 +Optimizer: Minuit
  141 +Tolerance = 0.001
  142 +Tolerance type = 1
  143 +Crab
  144 + Spectrum: PowerLaw2
  145 +0 Integral: 1.000e+00 0.000e+00 1.000e-05 1.000e+03 ( 1.000e-07)
  146 +1 Index: -2.000e+00 0.000e+00 -5.000e+00 -1.000e-01 ( 1.000e+00)
  147 +2 LowerLimit: 1.000e+02 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed
  148 +3 UpperLimit: 5.000e+05 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed
  149 +
  150 +gll_iem_v05
  151 + Spectrum: PowerLaw
  152 +4 Prefactor: 1.000e+00 0.000e+00 1.000e-05 1.000e+02 ( 1.000e+00)
  153 +5 Index: 0.000e+00 0.000e+00 -1.000e+00 1.000e+00 ( 1.000e+00)
  154 +6 Scale: 1.000e+02 0.000e+00 5.000e+01 2.000e+03 ( 1.000e+00) fixed
  155 +
  156 +iso_source_v05
  157 + Spectrum: FileFunction
  158 +7 Normalization: 1.000e+00 0.000e+00 0.000e+00 1.000e+01 ( 1.000e+00)
  159 +
  160 + Performing likelihood fit
  161 +
  162 + **********
  163 + ** 1 **SET PRINT 2.000
  164 + **********
  165 + **********
  166 + ** 2 **SET NOWARN
  167 + **********
  168 +
  169 + PARAMETER DEFINITIONS:
  170 + NO. NAME VALUE STEP SIZE LIMITS
  171 + 1 'Integral ' 1.0000 1.0000 .10000E-04 1000.0
  172 + 2 'Index ' -2.0000 1.0000 -5.0000 -.10000
  173 + 3 'Prefactor ' 1.0000 1.0000 .10000E-04 100.00
  174 + 4 'Index ' .0000 1.0000 -1.0000 1.0000
  175 + 5 'Normalizat' 1.0000 1.0000 .0000 10.000
  176 + **********
  177 + ** 3 **SET ERR .5000
  178 + **********
  179 + **********
  180 + ** 4 **SET GRAD 1.000
  181 + **********
  182 + **********
  183 + ** 5 **MINIMIZE 800.0 2.000
  184 + **********
  185 +
  186 + FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
  187 + START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. .20E-02
  188 +
  189 + FCN= 11570.88 FROM MIGRAD STATUS=INITIATE 1 CALLS 2 TOTAL
  190 + EDM= unknown STRATEGY= 1 NO ERROR MATRIX
  191 +
  192 + EXT PARAMETER CURRENT GUESS STEP FIRST
  193 + NO. NAME VALUE ERROR SIZE DERIVATIVE
  194 + 1 Integral 1.00000 1.0000 .43760E-01 -5385.8
  195 + 2 Index -2.0000 1.0000 .43485 -428.26
  196 + 3 Prefactor 1.00000 1.0000 .14092 -6446.5
  197 + 4 Index .0000 1.0000 1.5688 -733.82
  198 + 5 Normalizat 1.00000 1.0000 .46267 -348.37
  199 + ERR DEF= .500
  200 + NO ERROR MATRIX
  201 +
  202 + FCN= 11498.56 FROM MIGRAD STATUS=PROGRESS 8 CALLS 9 TOTAL
  203 + EDM= .24E+07 STRATEGY= 1 NO ERROR MATRIX
  204 +
  205 + EXT PARAMETER CURRENT GUESS STEP FIRST
  206 + NO. NAME VALUE ERROR SIZE DERIVATIVE
  207 + 1 Integral 1.0532 1.0000 .16626E-02 -3483.6
  208 + 2 Index -1.9689 1.0000 .13055E-01 -293.04
  209 + 3 Prefactor 1.2158 1.0000 .20638E-01 1980.5
  210 + 4 Index .28706 1.0000 .29116 961.06
  211 + 5 Normalizat 1.0364 1.0000 .12022E-01 -107.65
  212 + ERR DEF= .500
  213 +
  214 + FCN= 10808.51 FROM MIGRAD STATUS=PROGRESS 44 CALLS 45 TOTAL
  215 + EDM= .44E-04 STRATEGY=1 ERROR MATRIX UNCERTAINTY= 11.8%
  216 +
  217 + EXT PARAMETER APPROXIMATE STEP FIRST
  218 + NO. NAME VALUE ERROR SIZE DERIVATIVE
  219 + 1 Integral 25.587 1.4987 -.40989E-03 -.74472
  220 + 2 Index -2.1354 .46094E-01 .66442E-03 -.13219
  221 + 3 Prefactor .72871 .12597 -.25291E-02 -1.8487
  222 + 4 Index .86828E-01 .62329E-01 .10321E-01 -.28497
  223 + 5 Normalizat 2.8758 .62827 .17529E-01 -.11416
  224 + ERR DEF= .500
  225 +
  226 + MIGRAD MINIMIZATION HAS CONVERGED.
  227 +
  228 + MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
  229 + START COVARIANCE MATRIX CALCULATION.
  230 + EIGENVALUES OF SECOND-DERIVATIVE MATRIX:
  231 + -.6215E-01 .2198E+00 .7051E+00 .1404E+01 .2733E+01
  232 +
  233 + FCN= 10808.51 FROM MIGRAD STATUS=CONVERGED 93 CALLS 94 TOTAL
  234 + EDM= .13E-03 STRATEGY= 1 ERR MATRIX NOT POS-DEF
  235 +
  236 + EXT PARAMETER APPROXIMATE STEP FIRST
  237 + NO. NAME VALUE ERROR SIZE DERIVATIVE
  238 + 1 Integral 25.587 1.7100 .15489E-01 -.74473
  239 + 2 Index -2.1354 .50728E-01 .43485E-01 -.13214
  240 + 3 Prefactor .72871 .45763 .85444E-02 -1.8487
  241 + 4 Index .86828E-01 .25580 .51124E-01 -.28492
  242 + 5 Normalizat 2.8758 1.2558 .13663 -.11415
  243 + ERR DEF= .500
  244 +
  245 + EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 5 ERR DEF= .500
  246 + .292E+01 -.620E-01 -.430E+00 .246E+00 .898E+00
  247 + -.620E-01 .257E-02 .113E-01 -.658E-02 -.242E-01
  248 + -.430E+00 .113E-01 .210E+00 -.118E+00 -.561E+00
  249 + .246E+00 -.658E-02 -.118E+00 .669E-01 .312E+00
  250 + .898E+00 -.242E-01 -.561E+00 .312E+00 .162E+01
  251 + ERR MATRIX NOT POS-DEF
  252 +
  253 + PARAMETER CORRELATION COEFFICIENTS
  254 + NO. GLOBAL 1 2 3 4 5
  255 + 1 .79750 1.000 -.715 -.549 .556 .413
  256 + 2 .73270 -.715 1.000 .487 -.501 -.375
  257 + 3 .99776 -.549 .487 1.000 -.995 -.963
  258 + 4 .99652 .556 -.501 -.995 1.000 .946
  259 + 5 .97863 .413 -.375 -.963 .946 1.000
  260 + ERR MATRIX NOT POS-DEF
  261 +Final values:
  262 + Integral = 25.587
  263 + Index = -2.13544
  264 + Prefactor = 0.728708
  265 + Index = 0.0868276
  266 + Normalizat = 2.87582
  267 + **********
  268 + ** 6 **HESSE
  269 + **********
  270 + START COVARIANCE MATRIX CALCULATION.
  271 + EIGENVALUES OF SECOND-DERIVATIVE MATRIX:
  272 + .2895E-01 .2394E+00 .6905E+00 .1336E+01 .2705E+01
  273 + COVARIANCE MATRIX CALCULATED SUCCESSFULLY
  274 +
  275 + FCN= 10808.51 FROM HESSE STATUS=OK 31 CALLS 125 TOTAL
  276 + EDM= .68E-04 STRATEGY= 1 ERROR MATRIX ACCURATE
  277 +
  278 + EXT PARAMETER INTERNAL INTERNAL
  279 + NO. NAME VALUE ERROR STEP SIZE VALUE
  280 + 1 Integral 25.587 1.5107 .68454E-04 -1.2495
  281 + 2 Index -2.1354 .45941E-01 .14981E-03 .17002
  282 + 3 Prefactor .72871 .14405 .37751E-04 -1.3999
  283 + 4 Index .86828E-01 .71562E-01 .22668E-03 .86937E-01
  284 + 5 Normalizat 2.8758 .68294 .60570E-03 -.43878
  285 + ERR DEF= .500
  286 +
  287 + EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 5 ERR DEF= .500
  288 + .228E+01 -.426E-01 .644E-02 .399E-02 -.344E+00
  289 + -.426E-01 .211E-02 -.447E-03 -.911E-04 .922E-02
  290 + .644E-02 -.447E-03 .208E-01 -.958E-02 -.824E-01
  291 + .399E-02 -.911E-04 -.958E-02 .513E-02 .329E-01
  292 + -.344E+00 .922E-02 -.824E-01 .329E-01 .470E+00
  293 +
  294 + PARAMETER CORRELATION COEFFICIENTS
  295 + NO. GLOBAL 1 2 3 4 5
  296 + 1 .70952 1.000 -.613 .030 .037 -.332
  297 + 2 .63226 -.613 1.000 -.068 -.028 .293
  298 + 3 .97562 .030 -.068 1.000 -.929 -.834
  299 + 4 .95052 .037 -.028 -.929 1.000 .669
  300 + 5 .91975 -.332 .293 -.834 .669 1.000
  301 +Minuit fit quality: 3 estimated distance: 6.76066e-05
  302 +Minuit parameter uncertainties:
  303 + 1 1.5107
  304 + 2 0.0459439
  305 + 3 0.144061
  306 + 4 0.0716234
  307 + 5 0.685555
  308 +Generating SourceMap for Crab....................!
  309 +
  310 + EXTERNAL ERROR MATRIX. NDIM= 5 NPAR= 5 ERR DEF= .500
  311 + .228E+01 -.426E-01 .644E-02 .399E-02 -.344E+00
  312 + -.426E-01 .211E-02 -.447E-03 -.911E-04 .922E-02
  313 + .644E-02 -.447E-03 .208E-01 -.958E-02 -.824E-01
  314 + .399E-02 -.911E-04 -.958E-02 .513E-02 .329E-01
  315 + -.344E+00 .922E-02 -.824E-01 .329E-01 .470E+00
  316 + Fitted model parameters
  317 +
  318 +Crab
  319 + Spectrum: PowerLaw2
  320 +0 Integral: 2.559e+01 1.511e+00 1.000e-05 1.000e+03 ( 1.000e-07)
  321 +1 Index: -2.135e+00 4.594e-02 -5.000e+00 -1.000e-01 ( 1.000e+00)
  322 +2 LowerLimit: 1.000e+02 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed
  323 +3 UpperLimit: 5.000e+05 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed
  324 +
  325 +gll_iem_v05
  326 + Spectrum: PowerLaw
  327 +4 Prefactor: 7.287e-01 1.441e-01 1.000e-05 1.000e+02 ( 1.000e+00)
  328 +5 Index: 8.683e-02 7.162e-02 -1.000e+00 1.000e+00 ( 1.000e+00)
  329 +6 Scale: 1.000e+02 0.000e+00 5.000e+01 2.000e+03 ( 1.000e+00) fixed
  330 +
  331 +iso_source_v05
  332 + Spectrum: FileFunction
  333 +7 Normalization: 2.876e+00 6.856e-01 0.000e+00 1.000e+01 ( 1.000e+00)
  334 +
  335 +{'Crab': {'DEC': [22.0145],
  336 + 'Energy flux 1000': [0.0008451339118765881, 0.00012748598777983171],
  337 + 'Flux 1000': [1.8710656060226001e-07, 1.5607151735282336e-08],
  338 + 'Index': [-2.1354435334969293, 0.04594388482845476],
  339 + 'Integral': [25.587001794287488, 1.5107021246857777],
  340 + 'LowerLimit': [100.0],
  341 + 'Npred': 883.6820271808168,
  342 + 'RA': [83.63309999999998],
  343 + 'TS value': 1432.7716170280473,
  344 + 'UpperLimit': [500000.0]},
  345 + 'gll_iem_v05': {'Energy flux 1000': [0.13014639519669133,
  346 + 0.016019766902211626],
  347 + 'Flux 1000': [4.553183162897382e-05, 3.6174781586347625e-06],
  348 + 'Index': [0.08682762462159843, 0.07162344769818912],
  349 + 'Npred': 916.193626511584,
  350 + 'Prefactor': [0.7287081943702254, 0.1440608280265844],
  351 + 'Scale': [100.0],
  352 + 'TS value': 636.2220754477894},
  353 + 'iso_source_v05': {'Energy flux 1000': [0.09452113626537338,
  354 + 0.022532509845750396],
  355 + 'Flux 1000': [1.7700945209265112e-05,
  356 + 4.217913754868659e-06],
  357 + 'Normalization': [2.875818635678156, 0.6855552052125737],
  358 + 'Npred': 424.77274035060253,
  359 + 'TS value': 53.92443256788101},
  360 + 'logLike': 10808.505288286657}
  361 +
  362 +Writing used config file to: ./usedconf/used_Crab_select_mktime_cmap_ccube_ltcube_scaleltcube_expcube2_srcmap_pybinlike.conf
  363 +
  364 +----------------------------------------------
  365 +Total elapsed time was: 0.0427083333333 hours
  366 +
  367 +
  368 +
  369 +
  370 +
  371 +
  372 +2014-01-07T15:30:35: ********************************************************************************
  373 +2014-01-07T15:30:35: * ctlike *
  374 +2014-01-07T15:30:35: * ---------------------------------------------------------------------------- *
  375 +2014-01-07T15:30:35: * Version: 00-05-00 *
  376 +2014-01-07T15:30:35: ********************************************************************************
  377 +2014-01-07T15:35:17: +============+
  378 +2014-01-07T15:35:17: | Parameters |
  379 +2014-01-07T15:35:17: +============+
  380 +2014-01-07T15:35:17: infile ....................: observation_fermi.xml
  381 +2014-01-07T15:35:17: stat ......................: POISSON
  382 +2014-01-07T15:35:17: refit .....................: no
  383 +2014-01-07T15:35:17: caldb .....................: ${GAMMALIB}/share/caldb/cta
  384 +2014-01-07T15:35:17: irf .......................: kb_A_50h_v3
  385 +2014-01-07T15:35:17: srcmdl ....................: Crab_pathadjusted.xml
  386 +2014-01-07T15:35:17: outmdl ....................: crab_results_ctools.xml
  387 +2014-01-07T15:35:17: chatter ...................: 2
  388 +2014-01-07T15:35:17: clobber ...................: yes
  389 +2014-01-07T15:35:17: debug .....................: no
  390 +2014-01-07T15:35:17: mode ......................: ql
  391 +2014-01-07T15:35:17:
  392 +2014-01-07T15:35:17: +=============+
  393 +2014-01-07T15:35:17: | Observation |
  394 +2014-01-07T15:35:17: +=============+
  395 +2014-01-07T15:35:17: === GObservations ===
  396 +2014-01-07T15:35:17: Number of observations ....: 1
  397 +2014-01-07T15:35:17: Number of predicted events : 0
  398 +2014-01-07T15:35:17: === GLATObservation ===
  399 +2014-01-07T15:35:17: Name ......................: Crab
  400 +2014-01-07T15:35:17: Identifier ................: 00002
  401 +2014-01-07T15:35:17: Instrument ................: LAT
  402 +2014-01-07T15:35:17: Statistics ................: Poisson
  403 +2014-01-07T15:35:17: Ontime ....................: 0
  404 +2014-01-07T15:35:17: Livetime ..................: 0
  405 +2014-01-07T15:35:17: === GLATResponse ===
  406 +2014-01-07T15:35:17: Calibration database ......: /Users/schulz/software/gammalib/inst/lat/caldb/
  407 +2014-01-07T15:35:17: Response name .............: P7REPSOURCE_V15
  408 +2014-01-07T15:35:17: Section(s) ................: front & back
  409 +2014-01-07T15:35:17: === GLATLtCube ===
  410 +2014-01-07T15:35:17:
  411 +2014-01-07T15:35:17: +=================================+
  412 +2014-01-07T15:35:17: | Maximum likelihood optimisation |
  413 +2014-01-07T15:35:17: +=================================+
  414 +2014-01-07T15:35:20: Initial iteration: func=11537.5, Lambda=0.001
  415 +2014-01-07T15:35:22: Iteration 1: func=11259.1, Lambda=0.0001, delta=278.423, max(grad)=0 [-1]
  416 +2014-01-07T15:35:24: Iteration 2: func=11070.8, Lambda=1e-05, delta=188.277, max(grad)=0 [-1]
  417 +2014-01-07T15:35:26: Iteration 3: func=10915.6, Lambda=1e-06, delta=155.181, max(grad)=0 [-1]
  418 +2014-01-07T15:35:28: Iteration 4: func=10838.8, Lambda=1e-07, delta=76.8212, max(grad)=0 [-1]
  419 +2014-01-07T15:35:30: Iteration 5: func=10824.7, Lambda=1e-08, delta=14.1277, max(grad)=0 [-1]
  420 +2014-01-07T15:35:32: Iteration 6: func=10824.3, Lambda=1e-09, delta=0.369295, max(grad)=0.326075 [11]
  421 +2014-01-07T15:35:33: Iteration 7: func=10824.3, Lambda=1e-10, delta=0.00068754, max(grad)=0.000177934 [5]
  422 +2014-01-07T15:35:35: Iteration 8: func=10824.3, Lambda=1e-11, delta=6.65009e-06, max(grad)=0.00969042 [11]
  423 +2014-01-07T15:35:37: Iteration 9: func=10824.3, Lambda=1e-12, delta=9.97497e-08, max(grad)=0.00010425 [5]
  424 +2014-01-07T15:35:39:
  425 +2014-01-07T15:35:39: +=========================================+
  426 +2014-01-07T15:35:39: | Maximum likelihood optimization results |
  427 +2014-01-07T15:35:39: +=========================================+
  428 +2014-01-07T15:35:39: === GOptimizerLM ===
  429 +2014-01-07T15:35:39: Optimized function value ..: 10824.3
  430 +2014-01-07T15:35:39: Absolute precision ........: 1e-06
  431 +2014-01-07T15:35:39: Optimization status .......: converged
  432 +2014-01-07T15:35:39: Number of parameters ......: 15
  433 +2014-01-07T15:35:39: Number of free parameters .: 5
  434 +2014-01-07T15:35:39: Number of iterations ......: 9
  435 +2014-01-07T15:35:39: Lambda ....................: 1e-12
  436 +2014-01-07T15:35:39: Maximum log likelihood ....: -10824.3
  437 +2014-01-07T15:35:39: Observed events (Nobs) ...: 2225
  438 +2014-01-07T15:35:39: Predicted events (Npred) ..: 2225 (Nobs - Npred = 2.88554e-07)
  439 +2014-01-07T15:35:39: === GModels ===
  440 +2014-01-07T15:35:39: Number of models ..........: 3
  441 +2014-01-07T15:35:39: Number of parameters ......: 15
  442 +2014-01-07T15:35:39: === GModelSky ===
  443 +2014-01-07T15:35:39: Name ......................: iso_source_v05
  444 +2014-01-07T15:35:39: Instruments ...............: all
  445 +2014-01-07T15:35:39: Instrument scale factors ..: unity
  446 +2014-01-07T15:35:39: Observation identifiers ...: all
  447 +2014-01-07T15:35:39: Model type ................: DiffuseSource
  448 +2014-01-07T15:35:39: Model components ..........: "ConstantValue" * "FileFunction" * "Constant"
  449 +2014-01-07T15:35:39: Number of parameters ......: 3
  450 +2014-01-07T15:35:39: Number of spatial par's ...: 1
  451 +2014-01-07T15:35:39: Value ....................: 1 [0,10] (fixed,scale=1,gradient)
  452 +2014-01-07T15:35:39: Number of spectral par's ..: 1
  453 +2014-01-07T15:35:39: Normalization ............: 3.48618 +/- 0.694252 [0,10] (free,scale=1,gradient)
  454 +2014-01-07T15:35:39: Number of temporal par's ..: 1
  455 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  456 +2014-01-07T15:35:39: === GModelSky ===
  457 +2014-01-07T15:35:39: Name ......................: gll_iem_v05
  458 +2014-01-07T15:35:39: Instruments ...............: all
  459 +2014-01-07T15:35:39: Instrument scale factors ..: unity
  460 +2014-01-07T15:35:39: Observation identifiers ...: all
  461 +2014-01-07T15:35:39: Model type ................: DiffuseSource
  462 +2014-01-07T15:35:39: Model components ..........: "MapCubeFunction" * "PowerLaw" * "Constant"
  463 +2014-01-07T15:35:39: Number of parameters ......: 5
  464 +2014-01-07T15:35:39: Number of spatial par's ...: 1
  465 +2014-01-07T15:35:39: Normalization ............: 1 [0.001,1000] (fixed,scale=1,gradient)
  466 +2014-01-07T15:35:39: Number of spectral par's ..: 3
  467 +2014-01-07T15:35:39: Prefactor ................: 0.774972 +/- 0.156895 [1e-05,100] ph/cm2/s/MeV (free,scale=1,gradient)
  468 +2014-01-07T15:35:39: Index ....................: 0.0695083 +/- 0.0745843 [-1,1] (free,scale=1,gradient)
  469 +2014-01-07T15:35:39: PivotEnergy ..............: 100 [50,2000] MeV (fixed,scale=1,gradient)
  470 +2014-01-07T15:35:39: Number of temporal par's ..: 1
  471 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  472 +2014-01-07T15:35:39: === GModelSky ===
  473 +2014-01-07T15:35:39: Name ......................: Crab
  474 +2014-01-07T15:35:39: Instruments ...............: all
  475 +2014-01-07T15:35:39: Instrument scale factors ..: unity
  476 +2014-01-07T15:35:39: Observation identifiers ...: all
  477 +2014-01-07T15:35:39: Model type ................: PointSource
  478 +2014-01-07T15:35:39: Model components ..........: "SkyDirFunction" * "PowerLaw2" * "Constant"
  479 +2014-01-07T15:35:39: Number of parameters ......: 7
  480 +2014-01-07T15:35:39: Number of spatial par's ...: 2
  481 +2014-01-07T15:35:39: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
  482 +2014-01-07T15:35:39: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
  483 +2014-01-07T15:35:39: Number of spectral par's ..: 4
  484 +2014-01-07T15:35:39: Integral .................: 1.99229e-06 +/- 1.19129e-07 [1e-12,0.0001] ph/cm2/s (free,scale=1e-07,gradient)
  485 +2014-01-07T15:35:39: Index ....................: -2.05072 +/- 0.0491567 [-5,-0.1] (free,scale=1,gradient)
  486 +2014-01-07T15:35:39: LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1)
  487 +2014-01-07T15:35:39: UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1)
  488 +2014-01-07T15:35:39: Number of temporal par's ..: 1
  489 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient)
  490 +2014-01-07T15:35:39:
  491 +2014-01-07T15:35:39: +==============+
  492 +2014-01-07T15:35:39: | Save results |
  493 +2014-01-07T15:35:39: +==============+
  494 +2014-01-07T15:35:39:
  495 +2014-01-07T15:35:39: Application "ctlike" terminated after 304 wall clock seconds, consuming 276.453 seconds of CPU time.
inst/lat/test/data/p7repv15/crab_results_ctools.xml 0 → 100644
  View file @47a0e1a
  1 +<?xml version="1.0" encoding="UTF-8" standalone="no"?>
  2 +<source_library title="source library">
  3 + <source name="iso_source_v05" type="DiffuseSource">
  4 + <spectrum type="FileFunction" file="/Users/schulz/software/gammalib/inst/lat/test/data/p7repv15/iso_source_v05.txt">
  5 + <parameter name="Normalization" value="3.48618" error="0.694252" scale="1" min="0" max="10" free="1" />
  6 + </spectrum>
  7 + <spatialModel type="ConstantValue">
  8 + <parameter name="Value" value="1" scale="1" min="0" max="10" free="0" />
  9 + </spatialModel>
  10 + </source>
  11 + <source name="gll_iem_v05" type="DiffuseSource">
  12 + <spectrum type="PowerLaw">
  13 + <parameter name="Prefactor" value="0.774972" error="0.156895" scale="1" min="1e-05" max="100" free="1" />
  14 + <parameter name="Index" value="0.0695083" error="0.0745843" scale="1" min="-1" max="1" free="1" />
  15 + <parameter name="Scale" value="100" scale="1" min="50" max="2000" free="0" />
  16 + </spectrum>
  17 + <spatialModel type="MapCubeFunction" file="/Users/schulz/software/gammalib/inst/lat/test/data/p7repv15/gll_iem_v05.fit">
  18 + <parameter name="Normalization" value="1" scale="1" min="0.001" max="1000" free="0" />
  19 + </spatialModel>
  20 + </source>
  21 + <source name="Crab" type="PointSource">
  22 + <spectrum type="PowerLaw2">
  23 + <parameter name="Integral" value="19.9229" error="1.19129" scale="1e-07" min="1e-05" max="1000" free="1" />
  24 + <parameter name="Index" value="-2.05072" error="0.0491567" scale="1" min="-5" max="-0.1" free="1" />
  25 + <parameter name="LowerLimit" value="100" scale="1" min="10" max="1e+06" free="0" />
  26 + <parameter name="UpperLimit" value="500000" scale="1" min="10" max="1e+06" free="0" />
  27 + </spectrum>
  28 + <spatialModel type="SkyDirFunction">
  29 + <parameter name="RA" value="83.6331" scale="1" min="-360" max="360" free="0" />
  30 + <parameter name="DEC" value="22.0145" scale="1" min="-90" max="90" free="0" />
  31 + </spatialModel>
  32 + </source>
  33 +</source_library>
inst/lat/test/data/p7repv15/crab_results_sciencetools.xml 0 → 100644
  View file @47a0e1a
  1 +<?xml version="1.0" standalone="no"?>
  2 +<source_library title="source library">
  3 + <source name="Crab" type="PointSource">
  4 + <spectrum type="PowerLaw2">
  5 + <parameter error="1.510702125" free="1" max="1000" min="1e-05" name="Integral" scale="1e-07" value="25.58700179" />
  6 + <parameter error="0.04594388483" free="1" max="-0.1" min="-5" name="Index" scale="1" value="-2.135443533" />
  7 + <parameter free="0" max="1000000" min="10" name="LowerLimit" scale="1" value="100" />
  8 + <parameter free="0" max="1000000" min="10" name="UpperLimit" scale="1" value="500000" />
  9 + </spectrum>
  10 + <spatialModel type="SkyDirFunction">
  11 + <parameter free="0" max="360" min="-360" name="RA" scale="1" value="83.6331" />
  12 + <parameter free="0" max="90" min="-90" name="DEC" scale="1" value="22.0145" />
  13 + </spatialModel>
  14 + </source>
  15 + <source name="gll_iem_v05" type="DiffuseSource">
  16 + <spectrum type="PowerLaw">
  17 + <parameter error="0.144060828" free="1" max="100" min="1e-05" name="Prefactor" scale="1" value="0.7287081944" />
  18 + <parameter error="0.0716234477" free="1" max="1" min="-1" name="Index" scale="1" value="0.08682762462" />
  19 + <parameter free="0" max="2000" min="50" name="Scale" scale="1" value="100" />
  20 + </spectrum>
  21 + <spatialModel file="/lustre/fs8/group/hess/user/aschulz/lat_data/fssc_files/diffuse/p7rep/gll_iem_v05.fit" type="MapCubeFunction">
  22 + <parameter free="0" max="1000" min="0.001" name="Normalization" scale="1" value="1" />
  23 + </spatialModel>
  24 + </source>
  25 + <source name="iso_source_v05" type="DiffuseSource">
  26 + <spectrum file="/lustre/fs8/group/hess/user/aschulz/lat_data/fssc_files/diffuse/p7rep/iso_source_v05.txt" type="FileFunction">
  27 + <parameter error="0.6855552052" free="1" max="10" min="0" name="Normalization" scale="1" value="2.875818636" />
  28 + </spectrum>
  29 + <spatialModel type="ConstantValue">
  30 + <parameter free="0" max="10" min="0" name="Value" scale="1" value="1" />
  31 + </spatialModel>
  32 + </source>
  33 +</source_library>
inst/lat/test/data/p7repv15/iso_clean_back_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 6.14984e-07 1.95438e-09
  2 +71.1114 3.72408e-07 7.31677e-10
  3 +89.9329 2.24515e-07 4.22726e-10
  4 +113.736 1.34567e-07 2.53027e-10
  5 +143.839 8.01971e-08 1.55795e-10
  6 +181.91 4.73688e-08 9.4308e-11
  7 +230.058 2.79075e-08 6.39124e-11
  8 +290.949 1.58439e-08 3.84029e-11
  9 +367.956 8.68704e-09 2.69651e-11
  10 +465.346 4.71094e-09 1.58821e-11
  11 +588.512 2.52936e-09 1.17147e-11
  12 +744.277 1.30046e-09 6.68091e-12
  13 +941.27 6.26596e-10 4.68236e-12
  14 +1190.4 3.42749e-10 2.75409e-12
  15 +1505.47 2.13221e-10 2.184e-12
  16 +1903.94 1.17209e-10 1.29849e-12
  17 +2407.87 5.80198e-11 8.7009e-13
  18 +3045.17 3.34818e-11 5.30604e-13
  19 +3851.16 1.99078e-11 3.772e-13
  20 +4870.47 1.2684e-11 2.47964e-13
  21 +6159.57 8.35673e-12 1.79189e-13
  22 +7789.86 5.23671e-12 1.15874e-13
  23 +9851.66 3.14674e-12 8.02699e-14
  24 +12459.2 1.91127e-12 5.07799e-14
  25 +15756.8 1.152e-12 3.61796e-14
  26 +19927.3 6.66738e-13 2.20685e-14
  27 +25201.6 3.88757e-13 1.59768e-14
  28 +31871.8 2.08827e-13 9.27482e-15
  29 +40307.6 1.0111e-13 6.10698e-15
  30 +50976 5.74735e-14 3.70713e-15
  31 +64468.2 3.95881e-14 3.2122e-15
  32 +81531.5 1.93587e-14 1.78986e-15
  33 +103111 6.71432e-15 9.9232e-16
  34 +130402 3.46736e-15 5.54121e-16
  35 +164916 2.488e-15 5.08807e-16
  36 +208566 1.4619e-15 3.17424e-16
  37 +263768 7.35647e-16 1.99379e-16
  38 +333582 4.44108e-16 1.21822e-16
  39 +421873 3.10646e-16 7.81281e-17
  40 +533533 2.10663e-16 8.48059e-17
inst/lat/test/data/p7repv15/iso_clean_front_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 7.00767e-07 2.05152e-09
  2 +71.1114 3.5362e-07 6.5247e-10
  3 +89.9329 1.7497e-07 3.83742e-10
  4 +113.736 9.56748e-08 2.15936e-10
  5 +143.839 5.74716e-08 1.40301e-10
  6 +181.91 3.37344e-08 8.49085e-11
  7 +230.058 1.99782e-08 5.81583e-11
  8 +290.949 1.1179e-08 3.45801e-11
  9 +367.956 5.92402e-09 2.3831e-11
  10 +465.346 3.16501e-09 1.38598e-11
  11 +588.512 1.69805e-09 1.02603e-11
  12 +744.277 8.70565e-10 5.84559e-12
  13 +941.27 4.20393e-10 4.12681e-12
  14 +1190.4 2.25173e-10 2.3816e-12
  15 +1505.47 1.31679e-10 1.78915e-12
  16 +1903.94 7.59982e-11 1.0866e-12
  17 +2407.87 4.40828e-11 7.50919e-13
  18 +3045.17 2.75195e-11 4.84198e-13
  19 +3851.16 1.67597e-11 3.26318e-13
  20 +4870.47 1.11104e-11 2.19511e-13
  21 +6159.57 7.6758e-12 1.55582e-13
  22 +7789.86 4.94059e-12 1.02034e-13
  23 +9851.66 2.993e-12 6.89941e-14
  24 +12459.2 1.83156e-12 4.37814e-14
  25 +15756.8 1.10818e-12 3.10795e-14
  26 +19927.3 6.40836e-13 1.89177e-14
  27 +25201.6 3.65477e-13 1.3356e-14
  28 +31871.8 2.10468e-13 8.16437e-15
  29 +40307.6 1.23355e-13 6.0373e-15
  30 +50976 6.7553e-14 3.60071e-15
  31 +64468.2 3.46979e-14 2.57714e-15
  32 +81531.5 1.82312e-14 1.49538e-15
  33 +103111 9.66072e-15 1.13249e-15
  34 +130402 5.06368e-15 6.67285e-16
  35 +164916 2.7834e-15 5.52066e-16
  36 +208566 1.09541e-15 2.8482e-16
  37 +263768 3.10259e-16 1.6174e-16
  38 +333582 1.11669e-16 1.83793e-16
  39 +421873 4.61306e-17 2.16381e-16
  40 +533533 1.84268e-17 1.43993e-16
inst/lat/test/data/p7repv15/iso_clean_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 6.5739e-07 1.4138e-09
  2 +71.1114 3.65413e-07 4.89816e-10
  3 +89.9329 2.00736e-07 2.85946e-10
  4 +113.736 1.152e-07 1.66444e-10
  5 +143.839 6.89144e-08 1.0502e-10
  6 +181.91 4.05412e-08 6.35093e-11
  7 +230.058 2.38891e-08 4.32175e-11
  8 +290.949 1.34485e-08 2.58042e-11
  9 +367.956 7.24769e-09 1.79489e-11
  10 +465.346 3.89691e-09 1.0497e-11
  11 +588.512 2.08799e-09 7.75267e-12
  12 +744.277 1.06997e-09 4.41272e-12
  13 +941.27 5.14818e-10 3.09991e-12
  14 +1190.4 2.78497e-10 1.80526e-12
  15 +1505.47 1.6826e-10 1.39626e-12
  16 +1903.94 9.4569e-11 8.37462e-13
  17 +2407.87 5.03616e-11 5.69928e-13
  18 +3045.17 3.01447e-11 3.56621e-13
  19 +3851.16 1.81404e-11 2.47008e-13
  20 +4870.47 1.17935e-11 1.64283e-13
  21 +6159.57 7.96723e-12 1.17544e-13
  22 +7789.86 5.06304e-12 7.65707e-14
  23 +9851.66 3.0549e-12 5.23379e-14
  24 +12459.2 1.86351e-12 3.31707e-14
  25 +15756.8 1.1267e-12 2.35973e-14
  26 +19927.3 6.51709e-13 1.43751e-14
  27 +25201.6 3.74743e-13 1.02486e-14
  28 +31871.8 2.09642e-13 6.12858e-15
  29 +40307.6 1.1352e-13 4.31756e-15
  30 +50976 6.30806e-14 2.59093e-15
  31 +64468.2 3.63955e-14 1.99971e-15
  32 +81531.5 1.87688e-14 1.14987e-15
  33 +103111 8.60458e-15 7.83817e-16
  34 +130402 4.43249e-15 4.46745e-16
  35 +164916 2.55162e-15 3.65312e-16
  36 +208566 1.25221e-15 2.02086e-16
  37 +263768 5.21725e-16 1.27884e-16
  38 +333582 2.74784e-16 6.98538e-17
  39 +421873 1.75191e-16 4.53195e-17
  40 +533533 1.07328e-16 4.48274e-17
inst/lat/test/data/p7repv15/iso_source_back_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 6.83755e-07 1.67255e-09
  2 +71.1114 4.36476e-07 6.59398e-10
  3 +89.9329 2.79292e-07 3.98411e-10
  4 +113.736 1.68791e-07 2.40891e-10
  5 +143.839 9.57607e-08 1.45197e-10
  6 +181.91 5.67609e-08 8.8243e-11
  7 +230.058 3.49773e-08 6.15644e-11
  8 +290.949 2.04359e-08 3.79016e-11
  9 +367.956 1.1528e-08 2.70169e-11
  10 +465.346 6.18807e-09 1.58656e-11
  11 +588.512 3.15459e-09 1.14482e-11
  12 +744.277 1.58637e-09 6.42133e-12
  13 +941.27 7.60779e-10 4.53329e-12
  14 +1190.4 4.17295e-10 2.6617e-12
  15 +1505.47 2.61309e-10 2.0739e-12
  16 +1903.94 1.5083e-10 1.2661e-12
  17 +2407.87 8.21785e-11 8.53775e-13
  18 +3045.17 5.05611e-11 5.46006e-13
  19 +3851.16 3.11552e-11 3.79761e-13
  20 +4870.47 2.06313e-11 2.57045e-13
  21 +6159.57 1.37211e-11 1.82255e-13
  22 +7789.86 9.06803e-12 1.22851e-13
  23 +9851.66 6.05466e-12 8.96354e-14
  24 +12459.2 3.76973e-12 5.80444e-14
  25 +15756.8 2.2024e-12 4.0894e-14
  26 +19927.3 1.25479e-12 2.46933e-14
  27 +25201.6 6.84864e-13 1.70474e-14
  28 +31871.8 3.90811e-13 1.01917e-14
  29 +40307.6 2.45148e-13 7.48291e-15
  30 +50976 1.6491e-13 5.17197e-15
  31 +64468.2 1.21258e-13 4.10662e-15
  32 +81531.5 8.29367e-14 2.89173e-15
  33 +103111 5.42085e-14 2.16979e-15
  34 +130402 3.40396e-14 1.42783e-15
  35 +164916 2.03405e-14 1.04345e-15
  36 +208566 1.22391e-14 6.57557e-16
  37 +263768 7.24007e-15 4.6094e-16
  38 +333582 4.60703e-15 2.96301e-16
  39 +421873 3.14245e-15 1.8946e-16
  40 +533533 2.11738e-15 2.04747e-16
inst/lat/test/data/p7repv15/iso_source_front_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 7.91481e-07 2.07496e-09
  2 +71.1114 4.00085e-07 6.59868e-10
  3 +89.9329 1.97799e-07 3.81463e-10
  4 +113.736 1.09154e-07 2.15209e-10
  5 +143.839 6.66767e-08 1.37176e-10
  6 +181.91 4.00923e-08 8.41428e-11
  7 +230.058 2.46076e-08 5.73707e-11
  8 +290.949 1.40487e-08 3.4465e-11
  9 +367.956 7.48291e-09 2.32467e-11
  10 +465.346 4.06794e-09 1.36255e-11
  11 +588.512 2.24474e-09 1.00711e-11
  12 +744.277 1.16918e-09 5.77999e-12
  13 +941.27 5.66706e-10 4.00777e-12
  14 +1190.4 3.09224e-10 2.3318e-12
  15 +1505.47 1.85977e-10 1.72559e-12
  16 +1903.94 1.11821e-10 1.0756e-12
  17 +2407.87 6.85752e-11 7.34621e-13
  18 +3045.17 4.55161e-11 4.95739e-13
  19 +3851.16 2.95205e-11 3.33077e-13
  20 +4870.47 2.06427e-11 2.34089e-13
  21 +6159.57 1.4654e-11 1.67248e-13
  22 +7789.86 9.67737e-12 1.12273e-13
  23 +9851.66 6.00056e-12 7.76956e-14
  24 +12459.2 3.68184e-12 4.96075e-14
  25 +15756.8 2.21869e-12 3.57453e-14
  26 +19927.3 1.27902e-12 2.1877e-14
  27 +25201.6 7.0362e-13 1.53921e-14
  28 +31871.8 3.93642e-13 9.15151e-15
  29 +40307.6 2.39977e-13 6.99649e-15
  30 +50976 1.41487e-13 4.40287e-15
  31 +64468.2 8.10192e-14 3.24656e-15
  32 +81531.5 4.69263e-14 2.01564e-15
  33 +103111 2.749e-14 1.53595e-15
  34 +130402 1.54988e-14 9.37165e-16
  35 +164916 8.44393e-15 6.88848e-16
  36 +208566 4.60262e-15 4.05032e-16
  37 +263768 2.4421e-15 2.8329e-16
  38 +333582 1.45421e-15 1.72515e-16
  39 +421873 9.65276e-16 1.11303e-16
  40 +533533 6.28432e-16 1.16579e-16
inst/lat/test/data/p7repv15/iso_source_v05.txt 0 → 100644
  View file @47a0e1a
  1 +56.2289 7.30895e-07 1.30583e-09
  2 +71.1114 4.24266e-07 4.7121e-10
  3 +89.9329 2.44466e-07 2.79758e-10
  4 +113.736 1.41763e-07 1.64058e-10
  5 +143.839 8.21758e-08 1.00672e-10
  6 +181.91 4.87323e-08 6.11292e-11
  7 +230.058 2.98962e-08 4.2208e-11
  8 +290.949 1.72277e-08 2.5629e-11
  9 +367.956 9.4382e-09 1.77896e-11
  10 +465.346 5.07763e-09 1.04103e-11
  11 +588.512 2.67238e-09 7.59147e-12
  12 +744.277 1.36168e-09 4.29333e-12
  13 +941.27 6.55161e-10 3.0057e-12
  14 +1190.4 3.57951e-10 1.75506e-12
  15 +1505.47 2.19657e-10 1.33387e-12
  16 +1903.94 1.29329e-10 8.22123e-13
  17 +2407.87 7.46693e-11 5.57991e-13
  18 +3045.17 4.77307e-11 3.66569e-13
  19 +3851.16 3.02226e-11 2.50417e-13
  20 +4870.47 2.06202e-11 1.73058e-13
  21 +6159.57 1.4236e-11 1.23309e-13
  22 +7789.86 9.4042e-12 8.29352e-14
  23 +9851.66 6.01972e-12 5.87246e-14
  24 +12459.2 3.71736e-12 3.77277e-14
  25 +15756.8 2.21144e-12 2.69308e-14
  26 +19927.3 1.26828e-12 1.63845e-14
  27 +25201.6 6.9493e-13 1.14254e-14
  28 +31871.8 3.92037e-13 6.80312e-15
  29 +40307.6 2.41409e-13 5.08872e-15
  30 +50976 1.52356e-13 3.35563e-15
  31 +64468.2 9.94611e-14 2.57062e-15
  32 +81531.5 6.33069e-14 1.71297e-15
  33 +103111 3.98623e-14 1.29936e-15
  34 +130402 2.41898e-14 8.3516e-16
  35 +164916 1.41232e-14 6.15216e-16
  36 +208566 8.27701e-15 3.80584e-16
  37 +263768 4.75459e-15 2.66526e-16
  38 +333582 2.97763e-15 1.69153e-16
  39 +421873 2.02425e-15 1.08582e-16
  40 +533533 1.35632e-15 1.16762e-16
inst/lat/test/data/p7repv15/observation_fermi.xml 0 → 100644
  View file @47a0e1a
  1 +<observation_list title="observation list">
  2 + <observation name="Crab" id="00002" instrument="LAT">
  3 + <parameter name="CountsMap" file="srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits"/>
  4 + <parameter name="ExposureMap" file="expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits"/>
  5 + <parameter name="LiveTimeCube" file="ltcube_239557446.6_240162199.6_scaled.fits"/>
  6 + <parameter name="IRF" value="P7REP_SOURCE_V15"/>
  7 + </observation>
  8 +</observation_list>