{{lastupdated_at}} by {{lastupdated_by}}

{{>toc}}

h1. Validation of Pass 7 IRFs

h2. GammaLib-00-06-00

The initial validation of the Pass 7 IRFs has been done using GammaLib-00-06-00. Here the results for several datasets.

h3. week_00000_roi05_200_20000_10bin_extmod_galpropv02_Crab_ltcube

| *Parameter* | *gtlike* | *gammalib (with correction)* | *gammalib (w/o correction)* |
| Counts (for same flux) | 814.987 | 842.551 | 827.99 |
| Flux | 24.23 +/- 1.50 | 19.23 +/- 1.21 | - |
| Spectral index | 2.125 +/- 0.048 | 2.039 +/- 0.051 | - |
| Extragalactic diffuse | 2.286 +/- 0.462 | 2.896 +/- 0.442 | - |
| Galactic diffuse | 0.862 +/- 0.065 | 0.863 +/- 0.064 | - |

h4. gtlike results

<pre>
Crab:
Integral: 24.2369 +/- 1.49572
Index: -2.12467 +/- 0.0476363
LowerLimit: 100
UpperLimit: 500000
TS value: 1300.37
Flux: 2.41076e-06 +/- 1.50551e-07 photons/cm^2/s
Total number of model events: 814.987

Extragal_diffuse:
Normalization: 2.28565 +/- 0.462091
Flux: 0.000472695 +/- 9.55313e-05 photons/cm^2/s

Galactic_diffuse:
Value: 0.861504 +/- 0.0652024
Flux: 0.00041505 +/- 3.14108e-05 photons/cm^2/s

WARNING: Fit may be bad in range [100, 1258.93] (MeV)
WARNING: Fit may be bad in range [1995.26, 5011.87] (MeV)

Total number of observed counts: 2211
Total number of model events: 2211

-log(Likelihood): 10794.50987
</pre>

h4. GammaLib Results with efficiency correction

<pre>
=== GOptimizerLM ===
 Optimized function value ..: 10814.4
 Absolute precision ........: 1e-06
 Optimization status .......: converged
 Number of parameters ......: 13
 Number of free parameters .: 4
 Number of iterations ......: 9
 Lambda ....................: 1e-12
=== GObservations ===
 Number of observations ....: 1
 Number of predicted events : 2211
=== GModels ===
 Number of models ..........: 3
 Number of parameters ......: 13
=== GModelDiffuseSource ===
 Name ......................: Extragal_diffuse
 Instruments ...............: all
 Model type ................: "ConstantValue" * "FileFunction" * "Constant"
 Number of parameters ......: 3
 Number of spatial par's ...: 1
  Value ....................: 1 [0,10]  (fixed,scale=1,gradient)
 Number of spectral par's ..: 1
  Normalization ............: 2.89602 +/- 0.441842 [0,1000]  (free,scale=1,gradient)
 Number of temporal par's ..: 1
  Constant .................: 1 (relative value) (fixed,scale=1,gradient)
=== GModelDiffuseSource ===
 Name ......................: Galactic_diffuse
 Instruments ...............: all
 Model type ................: "MapCubeFunction" * "ConstantValue" * "Constant"
 Number of parameters ......: 3
 Number of spatial par's ...: 1
  Normalization ............: 1 [0.001,1000]  (fixed,scale=1,gradient)
 Number of spectral par's ..: 1
  Value ....................: 0.862952 +/- 0.0637931 [0,1000]  (free,scale=1,gradient)
 Number of temporal par's ..: 1
  Constant .................: 1 (relative value) (fixed,scale=1,gradient)
=== GModelPointSource ===
 Name ......................: Crab
 Instruments ...............: all
 Model type ................: "SkyDirFunction" * "PowerLaw2" * "Constant"
 Number of parameters ......: 7
 Number of spatial par's ...: 2
  RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
  DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
 Number of spectral par's ..: 4
  Integral .................: 1.92253e-06 +/- 1.20962e-07 [1e-14,0.0001] ph/cm2/s (free,scale=1e-07,gradient)
  Index ....................: -2.03863 +/- 0.0509223 [-5,5]  (free,scale=1,gradient)
  LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1)
  UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1)
 Number of temporal par's ..: 1
  Constant .................: 1 (relative value) (fixed,scale=1,gradient)
</pre>


h3. week_00000_P7_roi05_100_200_01bin_extmod_galprop_crab

Analysis performed with fixed spectral index over a narrow energy range of 100-200 MeV.
| *Parameter* | *gtlike* | *gammalib (with correction)* |
| -logL | 2395.656 | 2396.97 |
| Predicted counts | 658.017 | 658 |
| Flux | 26.10 +/- 3.10 | 15.58 +/- 1.95 |
| Spectral index | -2.1 | -2.1 |
| Extragalactic diffuse | 3.47 +/- 1.35 | 3.95 +/- 1.29 |
| Galactic diffuse | 0.68 +/- 0.36 | 0.71 +/- 0.40 |

h4. gtlike results

<pre>
Crab:
Integral: 26.1036 +/- 3.10488
Index: -2.1
LowerLimit: 100
UpperLimit: 500000
TS value: 108.858
Flux: 1.39271e-06 +/- 1.65655e-07 photons/cm^2/s

Extragal_diffuse:
Normalization: 3.47027 +/- 1.34533
Flux: 0.000427812 +/- 0.000165849 photons/cm^2/s

Galactic_diffuse:
Value: 0.679862 +/- 0.362457
Flux: 0.000134393 +/- 7.16491e-05 photons/cm^2/s

WARNING: Fit may be bad in range [100, 200] (MeV)

Total number of observed counts: 658
Total number of model events: 658.017

-log(Likelihood): 2395.656269
</pre>

h4. GammaLib Results with efficiency correction

<pre>
2012-07-11T21:55:18: === GOptimizerLM ===
2012-07-11T21:55:18:  Optimized function value ..: 2396.97
2012-07-11T21:55:18:  Absolute precision ........: 1e-06
2012-07-11T21:55:18:  Optimization status .......: converged
2012-07-11T21:55:18:  Number of parameters ......: 13
2012-07-11T21:55:18:  Number of free parameters .: 3
2012-07-11T21:55:18:  Number of iterations ......: 6
2012-07-11T21:55:18:  Lambda ....................: 1e-09
2012-07-11T21:55:18:  Maximum log likelihood ....: -2396.97
2012-07-11T21:55:18:  Npred .....................: 658
2012-07-11T21:55:18: === GModels ===
2012-07-11T21:55:18:  Number of models ..........: 3
2012-07-11T21:55:18:  Number of parameters ......: 13
2012-07-11T21:55:18: === GModelDiffuseSource ===
2012-07-11T21:55:18:  Name ......................: Extragal_diffuse
2012-07-11T21:55:18:  Instruments ...............: all
2012-07-11T21:55:18:  Model type ................: "ConstantValue" * "FileFunction" * "Constant"
2012-07-11T21:55:18:  Number of parameters ......: 3
2012-07-11T21:55:18:  Number of spatial par's ...: 1
2012-07-11T21:55:18:   Value ....................: 1 [0,10]  (fixed,scale=1,gradient)
2012-07-11T21:55:18:  Number of spectral par's ..: 1
2012-07-11T21:55:18:   Normalization ............: 3.95168 +/- 1.29363 [0,1000]  (free,scale=1,gradient)
2012-07-11T21:55:18:  Number of temporal par's ..: 1
2012-07-11T21:55:18:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
2012-07-11T21:55:18: === GModelDiffuseSource ===
2012-07-11T21:55:18:  Name ......................: Galactic_diffuse
2012-07-11T21:55:18:  Instruments ...............: all
2012-07-11T21:55:18:  Model type ................: "MapCubeFunction" * "ConstantValue" * "Constant"
2012-07-11T21:55:18:  Number of parameters ......: 3
2012-07-11T21:55:18:  Number of spatial par's ...: 1
2012-07-11T21:55:18:   Normalization ............: 1 [0.001,1000]  (fixed,scale=1,gradient)
2012-07-11T21:55:18:  Number of spectral par's ..: 1
2012-07-11T21:55:18:   Value ....................: 0.71251 +/- 0.398651 [0,1000]  (free,scale=1,gradient)
2012-07-11T21:55:18:  Number of temporal par's ..: 1
2012-07-11T21:55:18:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
2012-07-11T21:55:18: === GModelPointSource ===
2012-07-11T21:55:18:  Name ......................: Crab
2012-07-11T21:55:18:  Instruments ...............: all
2012-07-11T21:55:18:  Model type ................: "SkyDirFunction" * "PowerLaw2" * "Constant"
2012-07-11T21:55:18:  Number of parameters ......: 7
2012-07-11T21:55:18:  Number of spatial par's ...: 2
2012-07-11T21:55:18:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
2012-07-11T21:55:18:   DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
2012-07-11T21:55:18:  Number of spectral par's ..: 4
2012-07-11T21:55:18:   Integral .................: 1.55815e-06 +/- 1.94647e-07 [1e-14,0.0001] ph/cm2/s (free,scale=1e-07,gradient)
2012-07-11T21:55:18:   Index ....................: -2.1 [-5,5]  (fixed,scale=1,gradient)
2012-07-11T21:55:18:   LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1)
2012-07-11T21:55:18:   UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1)
2012-07-11T21:55:18:  Number of temporal par's ..: 1
2012-07-11T21:55:18:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
</pre>



h3. week_00000_P7_roi05_1000_2000_01bin_extmod_galprop_crab

Analysis performed with fixed spectral index over a narrow energy range from 1-2 GeV.
| *Parameter* | *gtlike* | *gammalib (with correction)* |
| -logL | 874.856 | 879.393 |
| Predicted counts | 213.042 | 213 |
| Flux | 27.72 +/- 3.30 | 25.58 +/- 3.21 |
| Spectral index | -2.1 | -2.1 |
| Extragalactic diffuse | 0.00019418 +/- 0.0433937 | 0 +/- 5.5129 |
| Galactic diffuse | 0.95 +/- 0.09 | 0.85 +/- 0.25 |

h4. gtlike results

<pre>
Crab:
Integral: 27.7185 +/- 3.29829
Index: -2.1
LowerLimit: 100
UpperLimit: 500000
TS value: 335.184
Flux: 1.17471e-07 +/- 1.39781e-08 photons/cm^2/s

Extragal_diffuse:
Normalization: 0.00019418 +/- 0.0433937
Flux: 7.9062e-10 +/- 1.7668e-07 photons/cm^2/s

Galactic_diffuse:
Value: 0.952101 +/- 0.092758
Flux: 2.87516e-05 +/- 2.8011e-06 photons/cm^2/s

WARNING: Fit may be bad in range [1000, 2000] (MeV)

Total number of observed counts: 213
Total number of model events: 213.042

-log(Likelihood): 874.8585611
</pre>


h4. GammaLib Results with efficiency correction

<pre>
2012-07-11T22:11:46: === GOptimizerLM ===
2012-07-11T22:11:46:  Optimized function value ..: 879.393
2012-07-11T22:11:46:  Absolute precision ........: 1e-06
2012-07-11T22:11:46:  Optimization status .......: converged
2012-07-11T22:11:46:  Number of parameters ......: 13
2012-07-11T22:11:46:  Number of free parameters .: 3
2012-07-11T22:11:46:  Number of iterations ......: 23
2012-07-11T22:11:46:  Lambda ....................: 1
2012-07-11T22:11:46:  Maximum log likelihood ....: -879.393
2012-07-11T22:11:46:  Npred .....................: 213
2012-07-11T22:11:46: === GModels ===
2012-07-11T22:11:46:  Number of models ..........: 3
2012-07-11T22:11:46:  Number of parameters ......: 13
2012-07-11T22:11:46: === GModelDiffuseSource ===
2012-07-11T22:11:46:  Name ......................: Extragal_diffuse
2012-07-11T22:11:46:  Instruments ...............: all
2012-07-11T22:11:46:  Model type ................: "ConstantValue" * "FileFunction" * "Constant"
2012-07-11T22:11:46:  Number of parameters ......: 3
2012-07-11T22:11:46:  Number of spatial par's ...: 1
2012-07-11T22:11:46:   Value ....................: 1 [0,10]  (fixed,scale=1,gradient)
2012-07-11T22:11:46:  Number of spectral par's ..: 1
2012-07-11T22:11:46:   Normalization ............: 0 +/- 5.5129 [0,1000]  (free,scale=1,gradient)
2012-07-11T22:11:46:  Number of temporal par's ..: 1
2012-07-11T22:11:46:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
2012-07-11T22:11:46: === GModelDiffuseSource ===
2012-07-11T22:11:46:  Name ......................: Galactic_diffuse
2012-07-11T22:11:46:  Instruments ...............: all
2012-07-11T22:11:46:  Model type ................: "MapCubeFunction" * "ConstantValue" * "Constant"
2012-07-11T22:11:46:  Number of parameters ......: 3
2012-07-11T22:11:46:  Number of spatial par's ...: 1
2012-07-11T22:11:46:   Normalization ............: 1 [0.001,1000]  (fixed,scale=1,gradient)
2012-07-11T22:11:46:  Number of spectral par's ..: 1
2012-07-11T22:11:46:   Value ....................: 0.84716 +/- 0.248145 [0,1000]  (free,scale=1,gradient)
2012-07-11T22:11:46:  Number of temporal par's ..: 1
2012-07-11T22:11:46:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
2012-07-11T22:11:46: === GModelPointSource ===
2012-07-11T22:11:46:  Name ......................: Crab
2012-07-11T22:11:46:  Instruments ...............: all
2012-07-11T22:11:46:  Model type ................: "SkyDirFunction" * "PowerLaw2" * "Constant"
2012-07-11T22:11:46:  Number of parameters ......: 7
2012-07-11T22:11:46:  Number of spatial par's ...: 2
2012-07-11T22:11:46:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
2012-07-11T22:11:46:   DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
2012-07-11T22:11:46:  Number of spectral par's ..: 4
2012-07-11T22:11:46:   Integral .................: 2.55776e-06 +/- 3.20673e-07 [1e-14,0.0001] ph/cm2/s (free,scale=1e-07,gradient)
2012-07-11T22:11:46:   Index ....................: -2.1 [-5,5]  (fixed,scale=1,gradient)
2012-07-11T22:11:46:   LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1)
2012-07-11T22:11:46:   UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1)
2012-07-11T22:11:46:  Number of temporal par's ..: 1
2012-07-11T22:11:46:   Constant .................: 1 (relative value) (fixed,scale=1,gradient)
</pre>