2020-01-15T14:18:00: ******************************************************************************** 2020-01-15T14:18:00: * ctulimit * 2020-01-15T14:18:00: * ---------------------------------------------------------------------------- * 2020-01-15T14:18:00: * Version: 1.6.3 * 2020-01-15T14:18:00: ******************************************************************************** 2020-01-15T14:18:00: +============+ 2020-01-15T14:18:00: | Parameters | 2020-01-15T14:18:00: +============+ 2020-01-15T14:18:00: inobs .....................: cube.fits 2020-01-15T14:18:00: inmodel ...................: binmodel.xml 2020-01-15T14:18:00: srcname ...................: hello 2020-01-15T14:18:00: parname ...................: 2020-01-15T14:18:00: expcube ...................: expcube.fits 2020-01-15T14:18:00: psfcube ...................: psfcube.fits 2020-01-15T14:18:00: edispcube .................: [not queried] 2020-01-15T14:18:00: bkgcube ...................: bkgcube.fits 2020-01-15T14:18:00: caldb .....................: [not queried] 2020-01-15T14:18:00: irf .......................: [not queried] 2020-01-15T14:18:00: edisp .....................: no 2020-01-15T14:18:00: confidence ................: 0.95 2020-01-15T14:18:00: sigma_min .................: 0.0 2020-01-15T14:18:00: sigma_max .................: 10.0 2020-01-15T14:18:00: eref ......................: 1.0 2020-01-15T14:18:00: emin ......................: 1.0 2020-01-15T14:18:00: emax ......................: 100.0 2020-01-15T14:18:00: tol .......................: 1e-5 2020-01-15T14:18:00: statistic .................: DEFAULT 2020-01-15T14:18:00: like_accuracy .............: 0.005 2020-01-15T14:18:00: max_iter ..................: 50 2020-01-15T14:18:00: chatter ...................: 4 2020-01-15T14:18:00: clobber ...................: yes 2020-01-15T14:18:00: debug .....................: yes 2020-01-15T14:18:00: mode ......................: ql 2020-01-15T14:18:00: logfile ...................: ctulimit.log 2020-01-15T14:18:00: 2020-01-15T14:18:00: WARNING: Energy dispersion will *NOT* be considered for the computation. To consider 2020-01-15T14:18:00: energy dispersion, please set the "edisp" parameter to "yes". Be aware that 2020-01-15T14:18:00: using energy dispersion will considerably slow down the computations. 2020-01-15T14:18:00: 2020-01-15T14:18:00: +===================+ 2020-01-15T14:18:00: | Input observation | 2020-01-15T14:18:00: +===================+ 2020-01-15T14:18:00: === GObservations === 2020-01-15T14:18:00: Number of observations ....: 1 2020-01-15T14:18:00: Number of models ..........: 2 2020-01-15T14:18:00: Number of observed events .: 61269 2020-01-15T14:18:00: Number of predicted events : 0 2020-01-15T14:18:00: === GCTAObservation === 2020-01-15T14:18:00: Name ......................: 2020-01-15T14:18:00: Identifier ................: 000001 2020-01-15T14:18:00: Instrument ................: CTA 2020-01-15T14:18:00: Event file ................: cube.fits 2020-01-15T14:18:00: Event type ................: CountsCube 2020-01-15T14:18:00: Statistic .................: cstat 2020-01-15T14:18:00: Ontime ....................: 18000 s 2020-01-15T14:18:00: Livetime ..................: 17640 s 2020-01-15T14:18:00: Deadtime correction .......: 0.98 2020-01-15T14:18:00: User energy range .........: undefined 2020-01-15T14:18:00: === GCTAPointing === 2020-01-15T14:18:00: Pointing direction ........: (RA,Dec)=(30.19,-33.49) 2020-01-15T14:18:00: === GCTAResponseCube === 2020-01-15T14:18:00: Energy dispersion .........: Not used 2020-01-15T14:18:00: === GCTACubeExposure === 2020-01-15T14:18:00: Filename ..................: expcube.fits 2020-01-15T14:18:00: Livetime ..................: 17640 sec 2020-01-15T14:18:00: === GEnergies === 2020-01-15T14:18:00: Number of energies ........: 36 2020-01-15T14:18:00: === GGti === 2020-01-15T14:18:00: Number of intervals .......: 1 2020-01-15T14:18:00: Ontime ....................: 18000 sec 2020-01-15T14:18:00: Elapsed time ..............: 18000 sec 2020-01-15T14:18:00: MJD range .................: 51544.5 - 51544.7083333333 s (TT) 2020-01-15T14:18:00: UTC range .................: 2000-01-01T11:58:56 - 2000-01-01T16:58:56 s (TT) 2020-01-15T14:18:00: Reference MJD .............: 51544.5 2020-01-15T14:18:00: === GSkyMap === 2020-01-15T14:18:00: Number of pixels ..........: 40000 2020-01-15T14:18:00: X axis dimension ..........: 200 2020-01-15T14:18:00: Y axis dimension ..........: 200 2020-01-15T14:18:00: Number of maps ............: 36 2020-01-15T14:18:00: Shape of maps .............: (36) 2020-01-15T14:18:00: === GCTACubePsf === 2020-01-15T14:18:00: Filename ..................: psfcube.fits 2020-01-15T14:18:00: === GEnergies === 2020-01-15T14:18:00: Number of energies ........: 36 2020-01-15T14:18:00: Number of delta bins ......: 200 2020-01-15T14:18:00: Delta range ...............: 1.875e-06 - 0.298501875 deg 2020-01-15T14:18:00: === GSkyMap === 2020-01-15T14:18:00: Number of pixels ..........: 121 2020-01-15T14:18:00: X axis dimension ..........: 11 2020-01-15T14:18:00: Y axis dimension ..........: 11 2020-01-15T14:18:00: Number of maps ............: 7200 2020-01-15T14:18:00: Shape of maps .............: (200, 36) 2020-01-15T14:18:00: === GCTACubeBackground === 2020-01-15T14:18:00: Filename ..................: bkgcube.fits 2020-01-15T14:18:00: === GEbounds === 2020-01-15T14:18:00: Number of intervals .......: 35 2020-01-15T14:18:00: Energy range ..............: 100 GeV - 199 TeV 2020-01-15T14:18:00: === GSkyMap === 2020-01-15T14:18:00: Number of pixels ..........: 40000 2020-01-15T14:18:00: X axis dimension ..........: 200 2020-01-15T14:18:00: Y axis dimension ..........: 200 2020-01-15T14:18:00: Number of maps ............: 35 2020-01-15T14:18:00: Shape of maps .............: (35) 2020-01-15T14:18:00: === GCTAResponseCache === 2020-01-15T14:18:00: Number of cached values ...: 0 2020-01-15T14:18:00: === GCTAEventCube === 2020-01-15T14:18:00: Number of events ..........: 61269 2020-01-15T14:18:00: Number of elements ........: 1400000 2020-01-15T14:18:00: Number of pixels ..........: 40000 2020-01-15T14:18:00: Number of energy bins .....: 35 2020-01-15T14:18:00: Time interval .............: 51544.5 - 51544.7083333333 days 2020-01-15T14:18:00: Pointing ..................: (RA,Dec)=(30.19,-33.49) 2020-01-15T14:18:00: === GEbounds === 2020-01-15T14:18:00: Number of intervals .......: 35 2020-01-15T14:18:00: Energy range ..............: 100 GeV - 199 TeV 2020-01-15T14:18:00: === GSkyMap === 2020-01-15T14:18:00: Number of pixels ..........: 40000 2020-01-15T14:18:00: X axis dimension ..........: 200 2020-01-15T14:18:00: Y axis dimension ..........: 200 2020-01-15T14:18:00: Number of maps ............: 35 2020-01-15T14:18:00: Shape of maps .............: (35) 2020-01-15T14:18:00: === GSkyRegions === 2020-01-15T14:18:00: Number of regions .........: 0 2020-01-15T14:18:00: === GModels === 2020-01-15T14:18:00: Number of models ..........: 2 2020-01-15T14:18:00: Number of parameters ......: 8 2020-01-15T14:18:00: === GModelSky === 2020-01-15T14:18:00: Name ......................: hello 2020-01-15T14:18:00: Instruments ...............: all 2020-01-15T14:18:00: Observation identifiers ...: all 2020-01-15T14:18:00: Model type ................: PointSource 2020-01-15T14:18:00: Model components ..........: "PointSource" * "FileFunction" * "Constant" 2020-01-15T14:18:00: Number of parameters ......: 4 2020-01-15T14:18:00: Number of spatial par's ...: 2 2020-01-15T14:18:00: RA .......................: 30.19 [-360,360] deg (fixed,scale=1) 2020-01-15T14:18:00: DEC ......................: -33.49 [-90,90] deg (fixed,scale=1) 2020-01-15T14:18:00: Number of spectral par's ..: 1 2020-01-15T14:18:00: Normalization ............: 1e-26 +/- 0 [1e-31,1e-16] (free,scale=1e-26,gradient) 2020-01-15T14:18:00: Number of temporal par's ..: 1 2020-01-15T14:18:00: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2020-01-15T14:18:00: Number of scale par's .....: 0 2020-01-15T14:18:00: === GCTAModelCubeBackground === 2020-01-15T14:18:00: Name ......................: BackgroundModel 2020-01-15T14:18:00: Instruments ...............: CTA, HESS, MAGIC, VERITAS 2020-01-15T14:18:00: Observation identifiers ...: all 2020-01-15T14:18:00: Model type ................: "PowerLaw" * "Constant" 2020-01-15T14:18:00: Number of parameters ......: 4 2020-01-15T14:18:00: Number of spectral par's ..: 3 2020-01-15T14:18:00: Prefactor ................: 1 +/- 0 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) 2020-01-15T14:18:00: Index ....................: 0 +/- 0 [-5,5] (free,scale=1,gradient) 2020-01-15T14:18:00: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) 2020-01-15T14:18:00: Number of temporal par's ..: 1 2020-01-15T14:18:00: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2020-01-15T14:18:00: 2020-01-15T14:18:00: +=============================+ 2020-01-15T14:18:00: | Compute best-fit likelihood | 2020-01-15T14:18:00: +=============================+ 2020-01-15T14:18:20: === GOptimizerLM === 2020-01-15T14:18:20: Optimized function value ..: 168654.670 2020-01-15T14:18:20: Absolute precision ........: 0.005 2020-01-15T14:18:20: Acceptable value decrease .: 2 2020-01-15T14:18:20: Optimization status .......: converged 2020-01-15T14:18:20: Number of parameters ......: 8 2020-01-15T14:18:20: Number of free parameters .: 3 2020-01-15T14:18:20: Number of iterations ......: 2 2020-01-15T14:18:20: Lambda ....................: 1e-05 2020-01-15T14:18:20: Maximum log likelihood ....: 168654.670 2020-01-15T14:18:20: === GModels === 2020-01-15T14:18:20: Number of models ..........: 2 2020-01-15T14:18:20: Number of parameters ......: 8 2020-01-15T14:18:20: === GModelSky === 2020-01-15T14:18:20: Name ......................: hello 2020-01-15T14:18:20: Instruments ...............: all 2020-01-15T14:18:20: Observation identifiers ...: all 2020-01-15T14:18:20: Model type ................: PointSource 2020-01-15T14:18:20: Model components ..........: "PointSource" * "FileFunction" * "Constant" 2020-01-15T14:18:20: Number of parameters ......: 4 2020-01-15T14:18:20: Number of spatial par's ...: 2 2020-01-15T14:18:20: RA .......................: 30.19 [-360,360] deg (fixed,scale=1) 2020-01-15T14:18:20: DEC ......................: -33.49 [-90,90] deg (fixed,scale=1) 2020-01-15T14:18:20: Number of spectral par's ..: 1 2020-01-15T14:18:20: Normalization ............: 5.26522319449013e-26 +/- 1.97622288842422e-24 [1e-31,1e-16] (free,scale=1e-26,gradient) 2020-01-15T14:18:20: Number of temporal par's ..: 1 2020-01-15T14:18:20: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2020-01-15T14:18:20: Number of scale par's .....: 0 2020-01-15T14:18:20: === GCTAModelCubeBackground === 2020-01-15T14:18:20: Name ......................: BackgroundModel 2020-01-15T14:18:20: Instruments ...............: CTA, HESS, MAGIC, VERITAS 2020-01-15T14:18:20: Observation identifiers ...: all 2020-01-15T14:18:20: Model type ................: "PowerLaw" * "Constant" 2020-01-15T14:18:20: Number of parameters ......: 4 2020-01-15T14:18:20: Number of spectral par's ..: 3 2020-01-15T14:18:20: Prefactor ................: 1.00229246180979 +/- 0.00996393841169144 [0.01,100] ph/cm2/s/MeV (free,scale=1,gradient) 2020-01-15T14:18:20: Index ....................: 0.00566383767416368 +/- 0.00569982550335703 [-5,5] (free,scale=1,gradient) 2020-01-15T14:18:20: PivotEnergy ..............: 1000000 MeV (fixed,scale=1000000,gradient) 2020-01-15T14:18:20: Number of temporal par's ..: 1 2020-01-15T14:18:20: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2020-01-15T14:18:20: 2020-01-15T14:18:20: +=========================+ 2020-01-15T14:18:20: | Find parameter brackets | 2020-01-15T14:18:20: +=========================+ 2020-01-15T14:18:33: Parameter range ...........: [5.26522319449013, 1981.48811161871] logL(1981.48811161871)=168667.21243489 2020-01-15T14:18:33: 2020-01-15T14:18:33: +=====================+ 2020-01-15T14:18:33: | Compute upper limit | 2020-01-15T14:18:33: +=====================+ 2020-01-15T14:18:33: Model name ................: hello 2020-01-15T14:18:33: Parameter name ............: Normalization 2020-01-15T14:18:33: Confidence level ..........: 95% 2020-01-15T14:18:33: Log-likelihood difference .: 1.92073 2020-01-15T14:18:33: Initial parameter range ...: [5.26522319449013, 1981.48811161871] 2020-01-15T14:18:43: Iteration 0 ...............: [5.26522319449013, 1981.48811161871] dlogL(993.376667406599)=2.86417350469856 2020-01-15T14:18:52: Iteration 1 ...............: [5.26522319449013, 993.376667406599] dlogL(499.320945300545)=-0.242350231972523 2020-01-15T14:19:01: Iteration 2 ...............: [499.320945300545, 993.376667406599] dlogL(746.348806353572)=1.21084060784779 2020-01-15T14:19:11: Iteration 3 ...............: [499.320945300545, 746.348806353572] dlogL(622.834875827058)=0.454152261081617 2020-01-15T14:19:20: Iteration 4 ...............: [499.320945300545, 622.834875827058] dlogL(561.077910563801)=0.097538116708165 2020-01-15T14:19:29: Iteration 5 ...............: [499.320945300545, 561.077910563801] dlogL(530.199427932173)=-0.0746197887929156 2020-01-15T14:19:38: Iteration 6 ...............: [530.199427932173, 561.077910563801] dlogL(545.638669247987)=0.0109220421290956 2020-01-15T14:19:48: Iteration 7 ...............: [530.199427932173, 545.638669247987] dlogL(537.91904859008)=-0.0319851393578574 2020-01-15T14:19:57: Iteration 8 ...............: [537.91904859008, 545.638669247987] dlogL(541.778858919034)=-0.0105653471255209 2020-01-15T14:20:07: Iteration 9 ...............: [541.778858919034, 545.638669247987] dlogL(543.70876408351)=0.000169933628058061 2020-01-15T14:20:26: Iteration 10 ..............: [541.778858919034, 543.70876408351] dlogL(542.743811501272)=-0.00519980280660093 2020-01-15T14:20:35: Iteration 11 ..............: [542.743811501272, 543.70876408351] dlogL(543.226287792391)=-0.00251546790241264 2020-01-15T14:20:44: Iteration 12 ..............: [543.226287792391, 543.70876408351] dlogL(543.467525937951)=-0.00117292554932646 2020-01-15T14:20:58: Iteration 13 ..............: [543.467525937951, 543.70876408351] dlogL(543.588145010731)=-0.000501520902616903 2020-01-15T14:21:07: Iteration 14 ..............: [543.588145010731, 543.70876408351] dlogL(543.648454547121)=-0.000165808596648276 2020-01-15T14:21:16: Iteration 15 ..............: [543.648454547121, 543.70876408351] dlogL(543.678609315315)=2.05388641916215e-06 2020-01-15T14:21:16: Final parameter ...........: 5.43679e-24 2020-01-15T14:21:16: 2020-01-15T14:21:16: +=====================+ 2020-01-15T14:21:16: | Upper limit results | 2020-01-15T14:21:16: +=====================+ 2020-01-15T14:21:16: Differential flux limit ...: 1.18288889756006e-20 ph/cm2/s/MeV at 1 TeV 2020-01-15T14:21:16: Integral flux limit .......: 1.61843265414486e-14 ph/cm2/s within [1-100] TeV 2020-01-15T14:21:16: Energy flux limit .........: 7.48232134622067e-14 erg/cm2/s within [1-100] TeV 2020-01-15T14:21:35: 2020-01-15T14:21:35: Application "ctulimit" terminated after 215 wall clock seconds, consuming 196.519 seconds of CPU time.