From 47a0e1a2c7bfe62e3a9a142a702e0a1e71fa88ef Mon Sep 17 00:00:00 2001 From: Jürgen Knödlseder Date: Wed, 8 Jan 2014 17:44:16 +0100 Subject: [PATCH] Add test files --- inst/lat/test/data/p7repv15/Crab.xml | 32 ++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt | 495 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/crab_results_ctools.xml | 33 +++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/crab_results_sciencetools.xml | 33 +++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_clean_back_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_clean_front_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_clean_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_source_back_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_source_front_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/iso_source_v05.txt | 40 ++++++++++++++++++++++++++++++++++++++++ inst/lat/test/data/p7repv15/observation_fermi.xml | 8 ++++++++ 11 files changed, 841 insertions(+), 0 deletions(-) create mode 100644 inst/lat/test/data/p7repv15/Crab.xml create mode 100644 inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt create mode 100644 inst/lat/test/data/p7repv15/crab_results_ctools.xml create mode 100644 inst/lat/test/data/p7repv15/crab_results_sciencetools.xml create mode 100644 inst/lat/test/data/p7repv15/iso_clean_back_v05.txt create mode 100644 inst/lat/test/data/p7repv15/iso_clean_front_v05.txt create mode 100644 inst/lat/test/data/p7repv15/iso_clean_v05.txt create mode 100644 inst/lat/test/data/p7repv15/iso_source_back_v05.txt create mode 100644 inst/lat/test/data/p7repv15/iso_source_front_v05.txt create mode 100644 inst/lat/test/data/p7repv15/iso_source_v05.txt create mode 100644 inst/lat/test/data/p7repv15/observation_fermi.xml diff --git a/inst/lat/test/data/p7repv15/Crab.xml b/inst/lat/test/data/p7repv15/Crab.xml new file mode 100644 index 0000000..7037973 --- /dev/null +++ b/inst/lat/test/data/p7repv15/Crab.xml @@ -0,0 +1,32 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt b/inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt new file mode 100644 index 0000000..defea23 --- /dev/null +++ b/inst/lat/test/data/p7repv15/compare_gtlike_ctlike.txt @@ -0,0 +1,495 @@ +LAT example analysis of Crab region +First analysis with Scince tools 09-32-05, then analysis of the same files with gtlike and ctlike: + +ra = 83.6331 +dec = 22.0145 + +rad = 7.5 +10x10 deg, 0.1 degperbin +zmax = 100 +emin = 100 MeV +emax = 10000 MeV + + +tmin = 239557446.6 MET in s +tmax= 240162199.6 MET in s +nebins=20 + + + + +gtlike results: + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +ctlike: + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +logfiles: + +------------------ running: pybinlike ------------------ +Parameters: + irfs = P7REP_SOURCE_V15 + expcube = ltcube_239557446.6_240162199.6_scaled.fits + srcmdl = Crab.xml + cmap = srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits + bexpmap = expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits + optimizer = Minuit + specfile = binned_specfile_py_Crab_ct.fits + sfile = binned_fitted_py_Crab_ct.xml + results = binned_results_py_Crab_ct.dat + scratchdir = none +{'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'} +{'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'} + +Parameters: + 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 + + Setting up binned likelihood + Initial model parameters + +Source maps: srcmap_Crab_ct_CEL_0.1dpb_AIT_100-10000MeV_P7REP_SOURCE_V15.fits +Exposure cube: ltcube_239557446.6_240162199.6_scaled.fits +Exposure map: expcube_239557446.6_240162199.6_100-10000_P7REP_SOURCE_V15.fits +IRFs: P7REP_SOURCE_V15 +Source model file: Crab.xml +Optimizer: Minuit +Tolerance = 0.001 +Tolerance type = 1 +Crab + Spectrum: PowerLaw2 +0 Integral: 1.000e+00 0.000e+00 1.000e-05 1.000e+03 ( 1.000e-07) +1 Index: -2.000e+00 0.000e+00 -5.000e+00 -1.000e-01 ( 1.000e+00) +2 LowerLimit: 1.000e+02 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed +3 UpperLimit: 5.000e+05 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed + +gll_iem_v05 + Spectrum: PowerLaw +4 Prefactor: 1.000e+00 0.000e+00 1.000e-05 1.000e+02 ( 1.000e+00) +5 Index: 0.000e+00 0.000e+00 -1.000e+00 1.000e+00 ( 1.000e+00) +6 Scale: 1.000e+02 0.000e+00 5.000e+01 2.000e+03 ( 1.000e+00) fixed + +iso_source_v05 + Spectrum: FileFunction +7 Normalization: 1.000e+00 0.000e+00 0.000e+00 1.000e+01 ( 1.000e+00) + + Performing likelihood fit + + ********** + ** 1 **SET PRINT 2.000 + ********** + ********** + ** 2 **SET NOWARN + ********** + + PARAMETER DEFINITIONS: + NO. NAME VALUE STEP SIZE LIMITS + 1 'Integral ' 1.0000 1.0000 .10000E-04 1000.0 + 2 'Index ' -2.0000 1.0000 -5.0000 -.10000 + 3 'Prefactor ' 1.0000 1.0000 .10000E-04 100.00 + 4 'Index ' .0000 1.0000 -1.0000 1.0000 + 5 'Normalizat' 1.0000 1.0000 .0000 10.000 + ********** + ** 3 **SET ERR .5000 + ********** + ********** + ** 4 **SET GRAD 1.000 + ********** + ********** + ** 5 **MINIMIZE 800.0 2.000 + ********** + + FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4. + START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. .20E-02 + + FCN= 11570.88 FROM MIGRAD STATUS=INITIATE 1 CALLS 2 TOTAL + EDM= unknown STRATEGY= 1 NO ERROR MATRIX + + EXT PARAMETER CURRENT GUESS STEP FIRST + NO. NAME VALUE ERROR SIZE DERIVATIVE + 1 Integral 1.00000 1.0000 .43760E-01 -5385.8 + 2 Index -2.0000 1.0000 .43485 -428.26 + 3 Prefactor 1.00000 1.0000 .14092 -6446.5 + 4 Index .0000 1.0000 1.5688 -733.82 + 5 Normalizat 1.00000 1.0000 .46267 -348.37 + ERR DEF= .500 + NO ERROR MATRIX + + FCN= 11498.56 FROM MIGRAD STATUS=PROGRESS 8 CALLS 9 TOTAL + EDM= .24E+07 STRATEGY= 1 NO ERROR MATRIX + + EXT PARAMETER CURRENT GUESS STEP FIRST + NO. NAME VALUE ERROR SIZE DERIVATIVE + 1 Integral 1.0532 1.0000 .16626E-02 -3483.6 + 2 Index -1.9689 1.0000 .13055E-01 -293.04 + 3 Prefactor 1.2158 1.0000 .20638E-01 1980.5 + 4 Index .28706 1.0000 .29116 961.06 + 5 Normalizat 1.0364 1.0000 .12022E-01 -107.65 + ERR DEF= .500 + + FCN= 10808.51 FROM MIGRAD STATUS=PROGRESS 44 CALLS 45 TOTAL + EDM= .44E-04 STRATEGY=1 ERROR MATRIX UNCERTAINTY= 11.8% + + EXT PARAMETER APPROXIMATE STEP FIRST + NO. NAME VALUE ERROR SIZE DERIVATIVE + 1 Integral 25.587 1.4987 -.40989E-03 -.74472 + 2 Index -2.1354 .46094E-01 .66442E-03 -.13219 + 3 Prefactor .72871 .12597 -.25291E-02 -1.8487 + 4 Index .86828E-01 .62329E-01 .10321E-01 -.28497 + 5 Normalizat 2.8758 .62827 .17529E-01 -.11416 + ERR DEF= .500 + + MIGRAD MINIMIZATION HAS CONVERGED. + + MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. + START COVARIANCE MATRIX CALCULATION. + EIGENVALUES OF SECOND-DERIVATIVE MATRIX: + -.6215E-01 .2198E+00 .7051E+00 .1404E+01 .2733E+01 + + FCN= 10808.51 FROM MIGRAD STATUS=CONVERGED 93 CALLS 94 TOTAL + EDM= .13E-03 STRATEGY= 1 ERR MATRIX NOT POS-DEF + + EXT PARAMETER APPROXIMATE STEP FIRST + NO. NAME VALUE ERROR SIZE DERIVATIVE + 1 Integral 25.587 1.7100 .15489E-01 -.74473 + 2 Index -2.1354 .50728E-01 .43485E-01 -.13214 + 3 Prefactor .72871 .45763 .85444E-02 -1.8487 + 4 Index .86828E-01 .25580 .51124E-01 -.28492 + 5 Normalizat 2.8758 1.2558 .13663 -.11415 + ERR DEF= .500 + + EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 5 ERR DEF= .500 + .292E+01 -.620E-01 -.430E+00 .246E+00 .898E+00 + -.620E-01 .257E-02 .113E-01 -.658E-02 -.242E-01 + -.430E+00 .113E-01 .210E+00 -.118E+00 -.561E+00 + .246E+00 -.658E-02 -.118E+00 .669E-01 .312E+00 + .898E+00 -.242E-01 -.561E+00 .312E+00 .162E+01 + ERR MATRIX NOT POS-DEF + + PARAMETER CORRELATION COEFFICIENTS + NO. GLOBAL 1 2 3 4 5 + 1 .79750 1.000 -.715 -.549 .556 .413 + 2 .73270 -.715 1.000 .487 -.501 -.375 + 3 .99776 -.549 .487 1.000 -.995 -.963 + 4 .99652 .556 -.501 -.995 1.000 .946 + 5 .97863 .413 -.375 -.963 .946 1.000 + ERR MATRIX NOT POS-DEF +Final values: + Integral = 25.587 + Index = -2.13544 + Prefactor = 0.728708 + Index = 0.0868276 + Normalizat = 2.87582 + ********** + ** 6 **HESSE + ********** + START COVARIANCE MATRIX CALCULATION. + EIGENVALUES OF SECOND-DERIVATIVE MATRIX: + .2895E-01 .2394E+00 .6905E+00 .1336E+01 .2705E+01 + COVARIANCE MATRIX CALCULATED SUCCESSFULLY + + FCN= 10808.51 FROM HESSE STATUS=OK 31 CALLS 125 TOTAL + EDM= .68E-04 STRATEGY= 1 ERROR MATRIX ACCURATE + + EXT PARAMETER INTERNAL INTERNAL + NO. NAME VALUE ERROR STEP SIZE VALUE + 1 Integral 25.587 1.5107 .68454E-04 -1.2495 + 2 Index -2.1354 .45941E-01 .14981E-03 .17002 + 3 Prefactor .72871 .14405 .37751E-04 -1.3999 + 4 Index .86828E-01 .71562E-01 .22668E-03 .86937E-01 + 5 Normalizat 2.8758 .68294 .60570E-03 -.43878 + ERR DEF= .500 + + EXTERNAL ERROR MATRIX. NDIM= 100 NPAR= 5 ERR DEF= .500 + .228E+01 -.426E-01 .644E-02 .399E-02 -.344E+00 + -.426E-01 .211E-02 -.447E-03 -.911E-04 .922E-02 + .644E-02 -.447E-03 .208E-01 -.958E-02 -.824E-01 + .399E-02 -.911E-04 -.958E-02 .513E-02 .329E-01 + -.344E+00 .922E-02 -.824E-01 .329E-01 .470E+00 + + PARAMETER CORRELATION COEFFICIENTS + NO. GLOBAL 1 2 3 4 5 + 1 .70952 1.000 -.613 .030 .037 -.332 + 2 .63226 -.613 1.000 -.068 -.028 .293 + 3 .97562 .030 -.068 1.000 -.929 -.834 + 4 .95052 .037 -.028 -.929 1.000 .669 + 5 .91975 -.332 .293 -.834 .669 1.000 +Minuit fit quality: 3 estimated distance: 6.76066e-05 +Minuit parameter uncertainties: + 1 1.5107 + 2 0.0459439 + 3 0.144061 + 4 0.0716234 + 5 0.685555 +Generating SourceMap for Crab....................! + + EXTERNAL ERROR MATRIX. NDIM= 5 NPAR= 5 ERR DEF= .500 + .228E+01 -.426E-01 .644E-02 .399E-02 -.344E+00 + -.426E-01 .211E-02 -.447E-03 -.911E-04 .922E-02 + .644E-02 -.447E-03 .208E-01 -.958E-02 -.824E-01 + .399E-02 -.911E-04 -.958E-02 .513E-02 .329E-01 + -.344E+00 .922E-02 -.824E-01 .329E-01 .470E+00 + Fitted model parameters + +Crab + Spectrum: PowerLaw2 +0 Integral: 2.559e+01 1.511e+00 1.000e-05 1.000e+03 ( 1.000e-07) +1 Index: -2.135e+00 4.594e-02 -5.000e+00 -1.000e-01 ( 1.000e+00) +2 LowerLimit: 1.000e+02 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed +3 UpperLimit: 5.000e+05 0.000e+00 1.000e+01 1.000e+06 ( 1.000e+00) fixed + +gll_iem_v05 + Spectrum: PowerLaw +4 Prefactor: 7.287e-01 1.441e-01 1.000e-05 1.000e+02 ( 1.000e+00) +5 Index: 8.683e-02 7.162e-02 -1.000e+00 1.000e+00 ( 1.000e+00) +6 Scale: 1.000e+02 0.000e+00 5.000e+01 2.000e+03 ( 1.000e+00) fixed + +iso_source_v05 + Spectrum: FileFunction +7 Normalization: 2.876e+00 6.856e-01 0.000e+00 1.000e+01 ( 1.000e+00) + +{'Crab': {'DEC': [22.0145], + 'Energy flux 1000': [0.0008451339118765881, 0.00012748598777983171], + 'Flux 1000': [1.8710656060226001e-07, 1.5607151735282336e-08], + 'Index': [-2.1354435334969293, 0.04594388482845476], + 'Integral': [25.587001794287488, 1.5107021246857777], + 'LowerLimit': [100.0], + 'Npred': 883.6820271808168, + 'RA': [83.63309999999998], + 'TS value': 1432.7716170280473, + 'UpperLimit': [500000.0]}, + 'gll_iem_v05': {'Energy flux 1000': [0.13014639519669133, + 0.016019766902211626], + 'Flux 1000': [4.553183162897382e-05, 3.6174781586347625e-06], + 'Index': [0.08682762462159843, 0.07162344769818912], + 'Npred': 916.193626511584, + 'Prefactor': [0.7287081943702254, 0.1440608280265844], + 'Scale': [100.0], + 'TS value': 636.2220754477894}, + 'iso_source_v05': {'Energy flux 1000': [0.09452113626537338, + 0.022532509845750396], + 'Flux 1000': [1.7700945209265112e-05, + 4.217913754868659e-06], + 'Normalization': [2.875818635678156, 0.6855552052125737], + 'Npred': 424.77274035060253, + 'TS value': 53.92443256788101}, + 'logLike': 10808.505288286657} + +Writing used config file to: ./usedconf/used_Crab_select_mktime_cmap_ccube_ltcube_scaleltcube_expcube2_srcmap_pybinlike.conf + +---------------------------------------------- +Total elapsed time was: 0.0427083333333 hours + + + + + + +2014-01-07T15:30:35: ******************************************************************************** +2014-01-07T15:30:35: * ctlike * +2014-01-07T15:30:35: * ---------------------------------------------------------------------------- * +2014-01-07T15:30:35: * Version: 00-05-00 * +2014-01-07T15:30:35: ******************************************************************************** +2014-01-07T15:35:17: +============+ +2014-01-07T15:35:17: | Parameters | +2014-01-07T15:35:17: +============+ +2014-01-07T15:35:17: infile ....................: observation_fermi.xml +2014-01-07T15:35:17: stat ......................: POISSON +2014-01-07T15:35:17: refit .....................: no +2014-01-07T15:35:17: caldb .....................: ${GAMMALIB}/share/caldb/cta +2014-01-07T15:35:17: irf .......................: kb_A_50h_v3 +2014-01-07T15:35:17: srcmdl ....................: Crab_pathadjusted.xml +2014-01-07T15:35:17: outmdl ....................: crab_results_ctools.xml +2014-01-07T15:35:17: chatter ...................: 2 +2014-01-07T15:35:17: clobber ...................: yes +2014-01-07T15:35:17: debug .....................: no +2014-01-07T15:35:17: mode ......................: ql +2014-01-07T15:35:17: +2014-01-07T15:35:17: +=============+ +2014-01-07T15:35:17: | Observation | +2014-01-07T15:35:17: +=============+ +2014-01-07T15:35:17: === GObservations === +2014-01-07T15:35:17: Number of observations ....: 1 +2014-01-07T15:35:17: Number of predicted events : 0 +2014-01-07T15:35:17: === GLATObservation === +2014-01-07T15:35:17: Name ......................: Crab +2014-01-07T15:35:17: Identifier ................: 00002 +2014-01-07T15:35:17: Instrument ................: LAT +2014-01-07T15:35:17: Statistics ................: Poisson +2014-01-07T15:35:17: Ontime ....................: 0 +2014-01-07T15:35:17: Livetime ..................: 0 +2014-01-07T15:35:17: === GLATResponse === +2014-01-07T15:35:17: Calibration database ......: /Users/schulz/software/gammalib/inst/lat/caldb/ +2014-01-07T15:35:17: Response name .............: P7REPSOURCE_V15 +2014-01-07T15:35:17: Section(s) ................: front & back +2014-01-07T15:35:17: === GLATLtCube === +2014-01-07T15:35:17: +2014-01-07T15:35:17: +=================================+ +2014-01-07T15:35:17: | Maximum likelihood optimisation | +2014-01-07T15:35:17: +=================================+ +2014-01-07T15:35:20: Initial iteration: func=11537.5, Lambda=0.001 +2014-01-07T15:35:22: Iteration 1: func=11259.1, Lambda=0.0001, delta=278.423, max(grad)=0 [-1] +2014-01-07T15:35:24: Iteration 2: func=11070.8, Lambda=1e-05, delta=188.277, max(grad)=0 [-1] +2014-01-07T15:35:26: Iteration 3: func=10915.6, Lambda=1e-06, delta=155.181, max(grad)=0 [-1] +2014-01-07T15:35:28: Iteration 4: func=10838.8, Lambda=1e-07, delta=76.8212, max(grad)=0 [-1] +2014-01-07T15:35:30: Iteration 5: func=10824.7, Lambda=1e-08, delta=14.1277, max(grad)=0 [-1] +2014-01-07T15:35:32: Iteration 6: func=10824.3, Lambda=1e-09, delta=0.369295, max(grad)=0.326075 [11] +2014-01-07T15:35:33: Iteration 7: func=10824.3, Lambda=1e-10, delta=0.00068754, max(grad)=0.000177934 [5] +2014-01-07T15:35:35: Iteration 8: func=10824.3, Lambda=1e-11, delta=6.65009e-06, max(grad)=0.00969042 [11] +2014-01-07T15:35:37: Iteration 9: func=10824.3, Lambda=1e-12, delta=9.97497e-08, max(grad)=0.00010425 [5] +2014-01-07T15:35:39: +2014-01-07T15:35:39: +=========================================+ +2014-01-07T15:35:39: | Maximum likelihood optimization results | +2014-01-07T15:35:39: +=========================================+ +2014-01-07T15:35:39: === GOptimizerLM === +2014-01-07T15:35:39: Optimized function value ..: 10824.3 +2014-01-07T15:35:39: Absolute precision ........: 1e-06 +2014-01-07T15:35:39: Optimization status .......: converged +2014-01-07T15:35:39: Number of parameters ......: 15 +2014-01-07T15:35:39: Number of free parameters .: 5 +2014-01-07T15:35:39: Number of iterations ......: 9 +2014-01-07T15:35:39: Lambda ....................: 1e-12 +2014-01-07T15:35:39: Maximum log likelihood ....: -10824.3 +2014-01-07T15:35:39: Observed events (Nobs) ...: 2225 +2014-01-07T15:35:39: Predicted events (Npred) ..: 2225 (Nobs - Npred = 2.88554e-07) +2014-01-07T15:35:39: === GModels === +2014-01-07T15:35:39: Number of models ..........: 3 +2014-01-07T15:35:39: Number of parameters ......: 15 +2014-01-07T15:35:39: === GModelSky === +2014-01-07T15:35:39: Name ......................: iso_source_v05 +2014-01-07T15:35:39: Instruments ...............: all +2014-01-07T15:35:39: Instrument scale factors ..: unity +2014-01-07T15:35:39: Observation identifiers ...: all +2014-01-07T15:35:39: Model type ................: DiffuseSource +2014-01-07T15:35:39: Model components ..........: "ConstantValue" * "FileFunction" * "Constant" +2014-01-07T15:35:39: Number of parameters ......: 3 +2014-01-07T15:35:39: Number of spatial par's ...: 1 +2014-01-07T15:35:39: Value ....................: 1 [0,10] (fixed,scale=1,gradient) +2014-01-07T15:35:39: Number of spectral par's ..: 1 +2014-01-07T15:35:39: Normalization ............: 3.48618 +/- 0.694252 [0,10] (free,scale=1,gradient) +2014-01-07T15:35:39: Number of temporal par's ..: 1 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient) +2014-01-07T15:35:39: === GModelSky === +2014-01-07T15:35:39: Name ......................: gll_iem_v05 +2014-01-07T15:35:39: Instruments ...............: all +2014-01-07T15:35:39: Instrument scale factors ..: unity +2014-01-07T15:35:39: Observation identifiers ...: all +2014-01-07T15:35:39: Model type ................: DiffuseSource +2014-01-07T15:35:39: Model components ..........: "MapCubeFunction" * "PowerLaw" * "Constant" +2014-01-07T15:35:39: Number of parameters ......: 5 +2014-01-07T15:35:39: Number of spatial par's ...: 1 +2014-01-07T15:35:39: Normalization ............: 1 [0.001,1000] (fixed,scale=1,gradient) +2014-01-07T15:35:39: Number of spectral par's ..: 3 +2014-01-07T15:35:39: Prefactor ................: 0.774972 +/- 0.156895 [1e-05,100] ph/cm2/s/MeV (free,scale=1,gradient) +2014-01-07T15:35:39: Index ....................: 0.0695083 +/- 0.0745843 [-1,1] (free,scale=1,gradient) +2014-01-07T15:35:39: PivotEnergy ..............: 100 [50,2000] MeV (fixed,scale=1,gradient) +2014-01-07T15:35:39: Number of temporal par's ..: 1 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient) +2014-01-07T15:35:39: === GModelSky === +2014-01-07T15:35:39: Name ......................: Crab +2014-01-07T15:35:39: Instruments ...............: all +2014-01-07T15:35:39: Instrument scale factors ..: unity +2014-01-07T15:35:39: Observation identifiers ...: all +2014-01-07T15:35:39: Model type ................: PointSource +2014-01-07T15:35:39: Model components ..........: "SkyDirFunction" * "PowerLaw2" * "Constant" +2014-01-07T15:35:39: Number of parameters ......: 7 +2014-01-07T15:35:39: Number of spatial par's ...: 2 +2014-01-07T15:35:39: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) +2014-01-07T15:35:39: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) +2014-01-07T15:35:39: Number of spectral par's ..: 4 +2014-01-07T15:35:39: Integral .................: 1.99229e-06 +/- 1.19129e-07 [1e-12,0.0001] ph/cm2/s (free,scale=1e-07,gradient) +2014-01-07T15:35:39: Index ....................: -2.05072 +/- 0.0491567 [-5,-0.1] (free,scale=1,gradient) +2014-01-07T15:35:39: LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1) +2014-01-07T15:35:39: UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1) +2014-01-07T15:35:39: Number of temporal par's ..: 1 +2014-01-07T15:35:39: Constant .................: 1 (relative value) (fixed,scale=1,gradient) +2014-01-07T15:35:39: +2014-01-07T15:35:39: +==============+ +2014-01-07T15:35:39: | Save results | +2014-01-07T15:35:39: +==============+ +2014-01-07T15:35:39: +2014-01-07T15:35:39: Application "ctlike" terminated after 304 wall clock seconds, consuming 276.453 seconds of CPU time. diff --git a/inst/lat/test/data/p7repv15/crab_results_ctools.xml b/inst/lat/test/data/p7repv15/crab_results_ctools.xml new file mode 100644 index 0000000..f7dfc5e --- /dev/null +++ b/inst/lat/test/data/p7repv15/crab_results_ctools.xml @@ -0,0 +1,33 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inst/lat/test/data/p7repv15/crab_results_sciencetools.xml b/inst/lat/test/data/p7repv15/crab_results_sciencetools.xml new file mode 100644 index 0000000..1609ee8 --- /dev/null +++ b/inst/lat/test/data/p7repv15/crab_results_sciencetools.xml @@ -0,0 +1,33 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inst/lat/test/data/p7repv15/iso_clean_back_v05.txt b/inst/lat/test/data/p7repv15/iso_clean_back_v05.txt new file mode 100644 index 0000000..8a047da --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_clean_back_v05.txt @@ -0,0 +1,40 @@ +56.2289 6.14984e-07 1.95438e-09 +71.1114 3.72408e-07 7.31677e-10 +89.9329 2.24515e-07 4.22726e-10 +113.736 1.34567e-07 2.53027e-10 +143.839 8.01971e-08 1.55795e-10 +181.91 4.73688e-08 9.4308e-11 +230.058 2.79075e-08 6.39124e-11 +290.949 1.58439e-08 3.84029e-11 +367.956 8.68704e-09 2.69651e-11 +465.346 4.71094e-09 1.58821e-11 +588.512 2.52936e-09 1.17147e-11 +744.277 1.30046e-09 6.68091e-12 +941.27 6.26596e-10 4.68236e-12 +1190.4 3.42749e-10 2.75409e-12 +1505.47 2.13221e-10 2.184e-12 +1903.94 1.17209e-10 1.29849e-12 +2407.87 5.80198e-11 8.7009e-13 +3045.17 3.34818e-11 5.30604e-13 +3851.16 1.99078e-11 3.772e-13 +4870.47 1.2684e-11 2.47964e-13 +6159.57 8.35673e-12 1.79189e-13 +7789.86 5.23671e-12 1.15874e-13 +9851.66 3.14674e-12 8.02699e-14 +12459.2 1.91127e-12 5.07799e-14 +15756.8 1.152e-12 3.61796e-14 +19927.3 6.66738e-13 2.20685e-14 +25201.6 3.88757e-13 1.59768e-14 +31871.8 2.08827e-13 9.27482e-15 +40307.6 1.0111e-13 6.10698e-15 +50976 5.74735e-14 3.70713e-15 +64468.2 3.95881e-14 3.2122e-15 +81531.5 1.93587e-14 1.78986e-15 +103111 6.71432e-15 9.9232e-16 +130402 3.46736e-15 5.54121e-16 +164916 2.488e-15 5.08807e-16 +208566 1.4619e-15 3.17424e-16 +263768 7.35647e-16 1.99379e-16 +333582 4.44108e-16 1.21822e-16 +421873 3.10646e-16 7.81281e-17 +533533 2.10663e-16 8.48059e-17 diff --git a/inst/lat/test/data/p7repv15/iso_clean_front_v05.txt b/inst/lat/test/data/p7repv15/iso_clean_front_v05.txt new file mode 100644 index 0000000..55c9f93 --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_clean_front_v05.txt @@ -0,0 +1,40 @@ +56.2289 7.00767e-07 2.05152e-09 +71.1114 3.5362e-07 6.5247e-10 +89.9329 1.7497e-07 3.83742e-10 +113.736 9.56748e-08 2.15936e-10 +143.839 5.74716e-08 1.40301e-10 +181.91 3.37344e-08 8.49085e-11 +230.058 1.99782e-08 5.81583e-11 +290.949 1.1179e-08 3.45801e-11 +367.956 5.92402e-09 2.3831e-11 +465.346 3.16501e-09 1.38598e-11 +588.512 1.69805e-09 1.02603e-11 +744.277 8.70565e-10 5.84559e-12 +941.27 4.20393e-10 4.12681e-12 +1190.4 2.25173e-10 2.3816e-12 +1505.47 1.31679e-10 1.78915e-12 +1903.94 7.59982e-11 1.0866e-12 +2407.87 4.40828e-11 7.50919e-13 +3045.17 2.75195e-11 4.84198e-13 +3851.16 1.67597e-11 3.26318e-13 +4870.47 1.11104e-11 2.19511e-13 +6159.57 7.6758e-12 1.55582e-13 +7789.86 4.94059e-12 1.02034e-13 +9851.66 2.993e-12 6.89941e-14 +12459.2 1.83156e-12 4.37814e-14 +15756.8 1.10818e-12 3.10795e-14 +19927.3 6.40836e-13 1.89177e-14 +25201.6 3.65477e-13 1.3356e-14 +31871.8 2.10468e-13 8.16437e-15 +40307.6 1.23355e-13 6.0373e-15 +50976 6.7553e-14 3.60071e-15 +64468.2 3.46979e-14 2.57714e-15 +81531.5 1.82312e-14 1.49538e-15 +103111 9.66072e-15 1.13249e-15 +130402 5.06368e-15 6.67285e-16 +164916 2.7834e-15 5.52066e-16 +208566 1.09541e-15 2.8482e-16 +263768 3.10259e-16 1.6174e-16 +333582 1.11669e-16 1.83793e-16 +421873 4.61306e-17 2.16381e-16 +533533 1.84268e-17 1.43993e-16 diff --git a/inst/lat/test/data/p7repv15/iso_clean_v05.txt b/inst/lat/test/data/p7repv15/iso_clean_v05.txt new file mode 100644 index 0000000..2912950 --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_clean_v05.txt @@ -0,0 +1,40 @@ +56.2289 6.5739e-07 1.4138e-09 +71.1114 3.65413e-07 4.89816e-10 +89.9329 2.00736e-07 2.85946e-10 +113.736 1.152e-07 1.66444e-10 +143.839 6.89144e-08 1.0502e-10 +181.91 4.05412e-08 6.35093e-11 +230.058 2.38891e-08 4.32175e-11 +290.949 1.34485e-08 2.58042e-11 +367.956 7.24769e-09 1.79489e-11 +465.346 3.89691e-09 1.0497e-11 +588.512 2.08799e-09 7.75267e-12 +744.277 1.06997e-09 4.41272e-12 +941.27 5.14818e-10 3.09991e-12 +1190.4 2.78497e-10 1.80526e-12 +1505.47 1.6826e-10 1.39626e-12 +1903.94 9.4569e-11 8.37462e-13 +2407.87 5.03616e-11 5.69928e-13 +3045.17 3.01447e-11 3.56621e-13 +3851.16 1.81404e-11 2.47008e-13 +4870.47 1.17935e-11 1.64283e-13 +6159.57 7.96723e-12 1.17544e-13 +7789.86 5.06304e-12 7.65707e-14 +9851.66 3.0549e-12 5.23379e-14 +12459.2 1.86351e-12 3.31707e-14 +15756.8 1.1267e-12 2.35973e-14 +19927.3 6.51709e-13 1.43751e-14 +25201.6 3.74743e-13 1.02486e-14 +31871.8 2.09642e-13 6.12858e-15 +40307.6 1.1352e-13 4.31756e-15 +50976 6.30806e-14 2.59093e-15 +64468.2 3.63955e-14 1.99971e-15 +81531.5 1.87688e-14 1.14987e-15 +103111 8.60458e-15 7.83817e-16 +130402 4.43249e-15 4.46745e-16 +164916 2.55162e-15 3.65312e-16 +208566 1.25221e-15 2.02086e-16 +263768 5.21725e-16 1.27884e-16 +333582 2.74784e-16 6.98538e-17 +421873 1.75191e-16 4.53195e-17 +533533 1.07328e-16 4.48274e-17 diff --git a/inst/lat/test/data/p7repv15/iso_source_back_v05.txt b/inst/lat/test/data/p7repv15/iso_source_back_v05.txt new file mode 100644 index 0000000..c3202ff --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_source_back_v05.txt @@ -0,0 +1,40 @@ +56.2289 6.83755e-07 1.67255e-09 +71.1114 4.36476e-07 6.59398e-10 +89.9329 2.79292e-07 3.98411e-10 +113.736 1.68791e-07 2.40891e-10 +143.839 9.57607e-08 1.45197e-10 +181.91 5.67609e-08 8.8243e-11 +230.058 3.49773e-08 6.15644e-11 +290.949 2.04359e-08 3.79016e-11 +367.956 1.1528e-08 2.70169e-11 +465.346 6.18807e-09 1.58656e-11 +588.512 3.15459e-09 1.14482e-11 +744.277 1.58637e-09 6.42133e-12 +941.27 7.60779e-10 4.53329e-12 +1190.4 4.17295e-10 2.6617e-12 +1505.47 2.61309e-10 2.0739e-12 +1903.94 1.5083e-10 1.2661e-12 +2407.87 8.21785e-11 8.53775e-13 +3045.17 5.05611e-11 5.46006e-13 +3851.16 3.11552e-11 3.79761e-13 +4870.47 2.06313e-11 2.57045e-13 +6159.57 1.37211e-11 1.82255e-13 +7789.86 9.06803e-12 1.22851e-13 +9851.66 6.05466e-12 8.96354e-14 +12459.2 3.76973e-12 5.80444e-14 +15756.8 2.2024e-12 4.0894e-14 +19927.3 1.25479e-12 2.46933e-14 +25201.6 6.84864e-13 1.70474e-14 +31871.8 3.90811e-13 1.01917e-14 +40307.6 2.45148e-13 7.48291e-15 +50976 1.6491e-13 5.17197e-15 +64468.2 1.21258e-13 4.10662e-15 +81531.5 8.29367e-14 2.89173e-15 +103111 5.42085e-14 2.16979e-15 +130402 3.40396e-14 1.42783e-15 +164916 2.03405e-14 1.04345e-15 +208566 1.22391e-14 6.57557e-16 +263768 7.24007e-15 4.6094e-16 +333582 4.60703e-15 2.96301e-16 +421873 3.14245e-15 1.8946e-16 +533533 2.11738e-15 2.04747e-16 diff --git a/inst/lat/test/data/p7repv15/iso_source_front_v05.txt b/inst/lat/test/data/p7repv15/iso_source_front_v05.txt new file mode 100644 index 0000000..a07d432 --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_source_front_v05.txt @@ -0,0 +1,40 @@ +56.2289 7.91481e-07 2.07496e-09 +71.1114 4.00085e-07 6.59868e-10 +89.9329 1.97799e-07 3.81463e-10 +113.736 1.09154e-07 2.15209e-10 +143.839 6.66767e-08 1.37176e-10 +181.91 4.00923e-08 8.41428e-11 +230.058 2.46076e-08 5.73707e-11 +290.949 1.40487e-08 3.4465e-11 +367.956 7.48291e-09 2.32467e-11 +465.346 4.06794e-09 1.36255e-11 +588.512 2.24474e-09 1.00711e-11 +744.277 1.16918e-09 5.77999e-12 +941.27 5.66706e-10 4.00777e-12 +1190.4 3.09224e-10 2.3318e-12 +1505.47 1.85977e-10 1.72559e-12 +1903.94 1.11821e-10 1.0756e-12 +2407.87 6.85752e-11 7.34621e-13 +3045.17 4.55161e-11 4.95739e-13 +3851.16 2.95205e-11 3.33077e-13 +4870.47 2.06427e-11 2.34089e-13 +6159.57 1.4654e-11 1.67248e-13 +7789.86 9.67737e-12 1.12273e-13 +9851.66 6.00056e-12 7.76956e-14 +12459.2 3.68184e-12 4.96075e-14 +15756.8 2.21869e-12 3.57453e-14 +19927.3 1.27902e-12 2.1877e-14 +25201.6 7.0362e-13 1.53921e-14 +31871.8 3.93642e-13 9.15151e-15 +40307.6 2.39977e-13 6.99649e-15 +50976 1.41487e-13 4.40287e-15 +64468.2 8.10192e-14 3.24656e-15 +81531.5 4.69263e-14 2.01564e-15 +103111 2.749e-14 1.53595e-15 +130402 1.54988e-14 9.37165e-16 +164916 8.44393e-15 6.88848e-16 +208566 4.60262e-15 4.05032e-16 +263768 2.4421e-15 2.8329e-16 +333582 1.45421e-15 1.72515e-16 +421873 9.65276e-16 1.11303e-16 +533533 6.28432e-16 1.16579e-16 diff --git a/inst/lat/test/data/p7repv15/iso_source_v05.txt b/inst/lat/test/data/p7repv15/iso_source_v05.txt new file mode 100644 index 0000000..948be1a --- /dev/null +++ b/inst/lat/test/data/p7repv15/iso_source_v05.txt @@ -0,0 +1,40 @@ +56.2289 7.30895e-07 1.30583e-09 +71.1114 4.24266e-07 4.7121e-10 +89.9329 2.44466e-07 2.79758e-10 +113.736 1.41763e-07 1.64058e-10 +143.839 8.21758e-08 1.00672e-10 +181.91 4.87323e-08 6.11292e-11 +230.058 2.98962e-08 4.2208e-11 +290.949 1.72277e-08 2.5629e-11 +367.956 9.4382e-09 1.77896e-11 +465.346 5.07763e-09 1.04103e-11 +588.512 2.67238e-09 7.59147e-12 +744.277 1.36168e-09 4.29333e-12 +941.27 6.55161e-10 3.0057e-12 +1190.4 3.57951e-10 1.75506e-12 +1505.47 2.19657e-10 1.33387e-12 +1903.94 1.29329e-10 8.22123e-13 +2407.87 7.46693e-11 5.57991e-13 +3045.17 4.77307e-11 3.66569e-13 +3851.16 3.02226e-11 2.50417e-13 +4870.47 2.06202e-11 1.73058e-13 +6159.57 1.4236e-11 1.23309e-13 +7789.86 9.4042e-12 8.29352e-14 +9851.66 6.01972e-12 5.87246e-14 +12459.2 3.71736e-12 3.77277e-14 +15756.8 2.21144e-12 2.69308e-14 +19927.3 1.26828e-12 1.63845e-14 +25201.6 6.9493e-13 1.14254e-14 +31871.8 3.92037e-13 6.80312e-15 +40307.6 2.41409e-13 5.08872e-15 +50976 1.52356e-13 3.35563e-15 +64468.2 9.94611e-14 2.57062e-15 +81531.5 6.33069e-14 1.71297e-15 +103111 3.98623e-14 1.29936e-15 +130402 2.41898e-14 8.3516e-16 +164916 1.41232e-14 6.15216e-16 +208566 8.27701e-15 3.80584e-16 +263768 4.75459e-15 2.66526e-16 +333582 2.97763e-15 1.69153e-16 +421873 2.02425e-15 1.08582e-16 +533533 1.35632e-15 1.16762e-16 diff --git a/inst/lat/test/data/p7repv15/observation_fermi.xml b/inst/lat/test/data/p7repv15/observation_fermi.xml new file mode 100644 index 0000000..35ca9b5 --- /dev/null +++ b/inst/lat/test/data/p7repv15/observation_fermi.xml @@ -0,0 +1,8 @@ + + + + + + + + -- libgit2 0.21.4