GammaLib

Here the reference gtlike result for a 5 degree ROI around the Crab position, including survey data until 1 June 2016:

<===> Crab:
<===> Integral: 29.24 +/- 0.0678804
<===> Index: -2.2446 +/- 0.0018335
<===> LowerLimit: 100
<===> UpperLimit: 500000
<===> TS value: 731874
<===> Flux: 5.52719e-06 +/- 1.59659e-08 photons/cm^2/s
<===> 
<===> Extragal_diffuse:
<===> Normalization: 2.66561 +/- 0.0232332
<===> Flux: 0.000784169 +/- 6.82459e-06 photons/cm^2/s
<===> 
<===> Galactic_diffuse:
<===> Value: 0.914973 +/- 0.00273087
<===> MapBase::readFitsFile: creating WcsMap2 object
<===> Flux: 0.000586389 +/- 1.75032e-06 photons/cm^2/s
<===> WARNING: Fit may be bad in range [60, 463.345] (MeV)
<===> WARNING: Fit may be bad in range [549.393, 5964.79] (MeV)
<===> WARNING: Fit may be bad in range [13979.5, 16575.7] (MeV)
<===> WARNING: Fit may be bad in range [19654, 27631.8] (MeV)
<===> WARNING: Fit may be bad in range [32763.4, 38848] (MeV)
<===> WARNING: Fit may be bad in range [253012, 300000] (MeV)
<===> 
<===> 
<===> Total number of observed counts: 1.24021e+06
<===> Total number of model events: 1.24021e+06
<===> 512  1249847.395  1240212.035
<===> 
<===> -log(Likelihood): -1249847.395
<===> 
<===> Writing fitted model to /users-data/knodlseder/glast/analysis/gammalib/proc/pass8_crab_roi05_60_300000_50bin_extmod_galmod_crab/results_like.xml
<===> Elapsed CPU time: 418.15

Here the corresponding ctlike run:

2016-06-21T15:07:42: +=========================================+
2016-06-21T15:07:42: | Maximum likelihood optimisation results |
2016-06-21T15:07:42: +=========================================+
2016-06-21T15:07:42: === GOptimizerLM ===
2016-06-21T15:07:42:  Optimized function value ..: -1249868.760
2016-06-21T15:07:42:  Absolute precision ........: 0.005
2016-06-21T15:07:42:  Acceptable value decrease .: 2
2016-06-21T15:07:42:  Optimization status .......: converged
2016-06-21T15:07:42:  Number of parameters ......: 13
2016-06-21T15:07:42:  Number of free parameters .: 4
2016-06-21T15:07:42:  Number of iterations ......: 9
2016-06-21T15:07:42:  Lambda ....................: 1e-12
2016-06-21T15:07:42:  Maximum log likelihood ....: 1249868.760
2016-06-21T15:07:42:  Observed events  (Nobs) ...: 1240212.000
2016-06-21T15:07:42:  Predicted events (Npred) ..: 1240211.999 (Nobs - Npred = 0.0014539)
2016-06-21T15:07:42: === GModels ===
2016-06-21T15:07:42:  Number of models ..........: 3
2016-06-21T15:07:42:  Number of parameters ......: 13
2016-06-21T15:07:42: === GModelSky ===
2016-06-21T15:07:42:  Name ......................: Extragal_diffuse
2016-06-21T15:07:42:  Instruments ...............: all
2016-06-21T15:07:42:  Instrument scale factors ..: unity
2016-06-21T15:07:42:  Observation identifiers ...: all
2016-06-21T15:07:42:  Model type ................: DiffuseSource
2016-06-21T15:07:42:  Model components ..........: "ConstantValue" * "FileFunction" * "Constant" 
2016-06-21T15:07:42:  Number of parameters ......: 3
2016-06-21T15:07:42:  Number of spatial par's ...: 1
2016-06-21T15:07:42:   Value ....................: 1 [0,10]  (fixed,scale=1,gradient)
2016-06-21T15:07:42:  Number of spectral par's ..: 1
2016-06-21T15:07:42:   Normalization ............: 2.71162 +/- 0.0227147 [0,1000]  (free,scale=1,gradient)
2016-06-21T15:07:42:  Number of temporal par's ..: 1
2016-06-21T15:07:42:   Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
2016-06-21T15:07:42: === GModelSky ===
2016-06-21T15:07:42:  Name ......................: Galactic_diffuse
2016-06-21T15:07:42:  Instruments ...............: all
2016-06-21T15:07:42:  Instrument scale factors ..: unity
2016-06-21T15:07:42:  Observation identifiers ...: all
2016-06-21T15:07:42:  Model type ................: DiffuseSource
2016-06-21T15:07:42:  Model components ..........: "MapCubeFunction" * "ConstantValue" * "Constant" 
2016-06-21T15:07:42:  Number of parameters ......: 3
2016-06-21T15:07:42:  Number of spatial par's ...: 1
2016-06-21T15:07:42:   Normalization ............: 1 [0.001,1000]  (fixed,scale=1,gradient)
2016-06-21T15:07:42:  Number of spectral par's ..: 1
2016-06-21T15:07:42:   Value ....................: 0.907689 +/- 0.00268146 [0,1000]  (free,scale=1,gradient)
2016-06-21T15:07:42:  Number of temporal par's ..: 1
2016-06-21T15:07:42:   Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
2016-06-21T15:07:42: === GModelSky ===
2016-06-21T15:07:42:  Name ......................: Crab
2016-06-21T15:07:42:  Instruments ...............: all
2016-06-21T15:07:42:  Instrument scale factors ..: unity
2016-06-21T15:07:42:  Observation identifiers ...: all
2016-06-21T15:07:42:  Model type ................: PointSource
2016-06-21T15:07:42:  Model components ..........: "SkyDirFunction" * "PowerLaw2" * "Constant" 
2016-06-21T15:07:42:  Number of parameters ......: 7
2016-06-21T15:07:42:  Number of spatial par's ...: 2
2016-06-21T15:07:42:   RA .......................: 83.6331 [-360,360] deg (fixed,scale=1)
2016-06-21T15:07:42:   DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1)
2016-06-21T15:07:42:  Number of spectral par's ..: 4
2016-06-21T15:07:42:   Integral .................: 2.88621e-06 +/- 6.62962e-09 [1e-14,0.0001] ph/cm2/s (free,scale=1e-07,gradient)
2016-06-21T15:07:42:   Index ....................: -2.23854 +/- 0.00175264 [-5,5]  (free,scale=1,gradient)
2016-06-21T15:07:42:   LowerLimit ...............: 100 [10,1e+06] MeV (fixed,scale=1)
2016-06-21T15:07:42:   UpperLimit ...............: 500000 [10,1e+06] MeV (fixed,scale=1)
2016-06-21T15:07:42:  Number of temporal par's ..: 1
2016-06-21T15:07:42:   Normalization ............: 1 (relative value) (fixed,scale=1,gradient)
2016-06-21T15:07:42: 
2016-06-21T15:07:42: +==============+
2016-06-21T15:07:42: | Save results |
2016-06-21T15:07:42: +==============+
2016-06-21T15:07:42: Save fitted parameters into file "crab_results.xml".
2016-06-21T15:07:42: 
2016-06-21T15:07:42: Application "ctlike" terminated after 23 wall clock seconds, consuming 12.1102 seconds of CPU time.

This looks pretty close, ctlike is however considerably faster than gtlike.