2021-05-28T13:20:32: ******************************************************************************** 2021-05-28T13:20:32: * ctlike * 2021-05-28T13:20:32: * ---------------------------------------------------------------------------- * 2021-05-28T13:20:32: * Version: 1.7.3 * 2021-05-28T13:20:32: ******************************************************************************** 2021-05-28T13:20:32: +============+ 2021-05-28T13:20:32: | Parameters | 2021-05-28T13:20:32: +============+ 2021-05-28T13:20:32: inobs .....................: ./obs.xml 2021-05-28T13:20:32: inmodel ...................: ./model_TS.xml 2021-05-28T13:20:32: expcube ...................: [not queried] 2021-05-28T13:20:32: psfcube ...................: [not queried] 2021-05-28T13:20:32: edispcube .................: [not queried] 2021-05-28T13:20:32: bkgcube ...................: [not queried] 2021-05-28T13:20:32: caldb .....................: prod3b-v2 2021-05-28T13:20:32: irf .......................: South_z20_0.5h 2021-05-28T13:20:32: edisp .....................: no 2021-05-28T13:20:32: outmodel ..................: ./outmodel1.xml 2021-05-28T13:20:32: outcovmat .................: NONE 2021-05-28T13:20:32: statistic .................: DEFAULT 2021-05-28T13:20:32: refit .....................: no 2021-05-28T13:20:32: like_accuracy .............: 0.005 2021-05-28T13:20:32: max_iter ..................: 50 2021-05-28T13:20:32: fix_spat_for_ts ...........: no 2021-05-28T13:20:32: nthreads ..................: 0 2021-05-28T13:20:32: chatter ...................: 3 2021-05-28T13:20:32: clobber ...................: yes 2021-05-28T13:20:32: debug .....................: no 2021-05-28T13:20:32: mode ......................: ql 2021-05-28T13:20:32: logfile ...................: ./logfile1.log 2021-05-28T13:20:32: 2021-05-28T13:20:32: WARNING: Energy dispersion will *NOT* be considered for the computation. To consider 2021-05-28T13:20:32: energy dispersion, please set the "edisp" parameter to "yes". Be aware that 2021-05-28T13:20:32: using energy dispersion will considerably slow down the computations. 2021-05-28T13:20:32: 2021-05-28T13:20:32: +====================+ 2021-05-28T13:20:32: | Input observations | 2021-05-28T13:20:32: +====================+ 2021-05-28T13:20:32: === GObservations === 2021-05-28T13:20:32: Number of observations ....: 2 2021-05-28T13:20:32: Number of models ..........: 2 2021-05-28T13:20:32: Number of observed events .: 253794 2021-05-28T13:20:32: Number of predicted events : 0 2021-05-28T13:20:32: === GCTAObservation === 2021-05-28T13:20:32: Name ......................: src2330 2021-05-28T13:20:32: Identifier ................: 01 2021-05-28T13:20:32: Instrument ................: CTA 2021-05-28T13:20:32: Event file ................: ./event1.fits 2021-05-28T13:20:32: Event type ................: EventList 2021-05-28T13:20:32: Statistic .................: cstat 2021-05-28T13:20:32: Ontime ....................: 3420 s 2021-05-28T13:20:32: Livetime ..................: 3351.6 s 2021-05-28T13:20:32: Deadtime correction .......: 0.98 2021-05-28T13:20:32: User energy range .........: undefined 2021-05-28T13:20:32: === GCTAPointing === 2021-05-28T13:20:32: Pointing direction ........: (RA,Dec)=(186.151,-4.82077) 2021-05-28T13:20:32: === GCTAResponseIrf === 2021-05-28T13:20:32: Caldb mission .............: cta 2021-05-28T13:20:32: Caldb instrument ..........: prod3b-v2 2021-05-28T13:20:32: Response name .............: South_z20_0.5h 2021-05-28T13:20:32: Energy dispersion .........: Not used 2021-05-28T13:20:32: Safe energy range .........: undefined 2021-05-28T13:20:32: === GCTAEventList === 2021-05-28T13:20:32: Number of events ..........: 126897 (disposed in "./event1.fits") 2021-05-28T13:20:32: Time interval .............: 51544.5 - 51544.5395833333 days 2021-05-28T13:20:32: Energy interval ...........: 0.1 - 100 TeV 2021-05-28T13:20:32: Region of interest ........: RA=186.151, DEC=-4.82077 [0,0] Radius=2.16 deg 2021-05-28T13:20:32: === GSkyRegions === 2021-05-28T13:20:32: Number of regions .........: 0 2021-05-28T13:20:32: === GCTAObservation === 2021-05-28T13:20:32: Name ......................: src2330 2021-05-28T13:20:32: Identifier ................: 02 2021-05-28T13:20:32: Instrument ................: CTA 2021-05-28T13:20:32: Event file ................: ./event2.fits 2021-05-28T13:20:32: Event type ................: EventList 2021-05-28T13:20:32: Statistic .................: cstat 2021-05-28T13:20:32: Ontime ....................: 3420 s 2021-05-28T13:20:32: Livetime ..................: 3351.6 s 2021-05-28T13:20:32: Deadtime correction .......: 0.98 2021-05-28T13:20:32: User energy range .........: undefined 2021-05-28T13:20:32: === GCTAPointing === 2021-05-28T13:20:32: Pointing direction ........: (RA,Dec)=(186.151,-7.61578) 2021-05-28T13:20:32: === GCTAResponseIrf === 2021-05-28T13:20:32: Caldb mission .............: cta 2021-05-28T13:20:32: Caldb instrument ..........: prod3b-v2 2021-05-28T13:20:32: Response name .............: South_z20_0.5h 2021-05-28T13:20:32: Energy dispersion .........: Not used 2021-05-28T13:20:32: Safe energy range .........: undefined 2021-05-28T13:20:32: === GCTAEventList === 2021-05-28T13:20:32: Number of events ..........: 126897 (disposed in "./event2.fits") 2021-05-28T13:20:32: Time interval .............: 51544.5 - 51544.5395833333 days 2021-05-28T13:20:32: Energy interval ...........: 0.1 - 100 TeV 2021-05-28T13:20:32: Region of interest ........: RA=186.151, DEC=-7.61578 [0,0] Radius=2.16 deg 2021-05-28T13:20:32: === GSkyRegions === 2021-05-28T13:20:32: Number of regions .........: 0 2021-05-28T13:20:32: 2021-05-28T13:20:32: +=================================+ 2021-05-28T13:20:32: | Maximum likelihood optimisation | 2021-05-28T13:20:32: +=================================+ 2021-05-28T13:20:34: >Iteration 0: -logL=1482911.076, Lambda=1.0e-03 2021-05-28T13:20:34: Parameter "Prefactor" drives optimization step (step=0.00167941) 2021-05-28T13:20:35: >Iteration 1: -logL=1282289.201, Lambda=1.0e-03, delta=200621.876, step=1.7e-03, max(|grad|)=-2965.509270 [Prefactor:2] 2021-05-28T13:20:35: Parameter "Prefactor" does not drive optimization step anymore. 2021-05-28T13:20:35: Parameter "Index" drives optimization step (step=0.000757938) 2021-05-28T13:20:37: >Iteration 2: -logL=1282289.186, Lambda=1.0e-04, delta=0.015, step=7.6e-04, max(|grad|)=-126312.015803 [Prefactor:2] 2021-05-28T13:20:37: Parameter "Index" does not drive optimization step anymore. 2021-05-28T13:20:37: Parameter "Prefactor" drives optimization step (step=2.21465e-09) 2021-05-28T13:20:37: Parameter "Index" hits maximum: 5.04412 > 5 (1) 2021-05-28T13:20:38: Iteration 3: -logL=1282289.227, Lambda=1.0e-05, delta=-0.041, step=2.2e-09, max(|grad|)=270030.654272 [Prefactor:2] (stalled) 2021-05-28T13:20:38: Parameter "Prefactor" does not drive optimization step anymore. 2021-05-28T13:20:38: Parameter "Prefactor" hits minimum: -0.000308903 < 1e-07 (1) 2021-05-28T13:20:38: Parameter "Index" hits maximum: 1420.66 > 5 (2) 2021-05-28T13:20:39: >Iteration 4: -logL=1282288.364, Lambda=1.0e-04, delta=0.863, step=1.0e+00, max(|grad|)=276896.798511 [Prefactor:2] 2021-05-28T13:20:39: Parameter "Prefactor" hits minimum: -0.000312891 < 1e-07 (2) 2021-05-28T13:20:39: Parameter "Index" hits maximum: 1439.88 > 5 (3) 2021-05-28T13:20:41: Iteration 5: -logL=1282288.365, Lambda=1.0e-05, delta=-0.001, step=1.0e+00, max(|grad|)=276901.443904 [Prefactor:2] (stalled) 2021-05-28T13:20:41: Parameter "Prefactor" hits minimum: -0.000309536 < 1e-07 (3) 2021-05-28T13:20:41: Parameter "Index" hits maximum 5 more than 3 times. Fix parameter at maximum for now. 2021-05-28T13:20:42: >Iteration 6: -logL=1282288.365, Lambda=1.0e-04, delta=0.000, step=1.0e+00, max(|grad|)=276900.075884 [Prefactor:2] 2021-05-28T13:20:42: Free parameter "Index" after convergence was reached with frozen parameter. 2021-05-28T13:20:42: Parameter "Prefactor" hits minimum 1e-07 more than 3 times. Fix parameter at minimum for now. 2021-05-28T13:20:43: >Iteration 7: -logL=1282288.361, Lambda=1.0e-03, delta=0.003, step=1.0e+00, max(|grad|)=-120.467992 [Index:7] 2021-05-28T13:20:43: Free parameter "Prefactor" after convergence was reached with frozen parameter. 2021-05-28T13:20:43: Free parameter "Index" after convergence was reached with frozen parameter. 2021-05-28T13:20:43: 2021-05-28T13:20:43: +==================+ 2021-05-28T13:20:43: | Curvature matrix | 2021-05-28T13:20:43: +==================+ 2021-05-28T13:20:43: === GMatrixSparse === 2021-05-28T13:20:43: Number of rows ............: 10 2021-05-28T13:20:43: Number of columns .........: 10 2021-05-28T13:20:43: Number of nonzero elements : 9 2021-05-28T13:20:43: Number of allocated cells .: 521 2021-05-28T13:20:43: Memory block size .........: 512 2021-05-28T13:20:43: Sparse matrix fill ........: 0.09 2021-05-28T13:20:43: Pending element ...........: 0 2021-05-28T13:20:43: Fill stack size ...........: 0 (none) 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0.000532064741279702, 0, 0, 0.555318882851058, -1.12310105792725, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0.555318882851058, 0, 0, 255292.212934206, -382215.650721677, 0, 0 2021-05-28T13:20:43: 0, 0, 0, -1.12310105792725, 0, 0, -382215.650721677, 778455.658442613, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:43: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 2021-05-28T13:20:45: 2021-05-28T13:20:45: +====================================+ 2021-05-28T13:20:45: | Maximum likelihood re-optimisation | 2021-05-28T13:20:45: +====================================+ 2021-05-28T13:20:45: >Iteration 0: -logL=1282288.334, Lambda=1.0e-03 2021-05-28T13:20:46: >Iteration 1: -logL=1282288.324, Lambda=1.0e-03, delta=0.009, step=1.0e+00, max(|grad|)=0.023781 [Index:1] 2021-05-28T13:20:47: Iteration 2: -logL=1282288.324, Lambda=1.0e-04, delta=-0.000, step=1.0e+00, max(|grad|)=-0.000075 [Index:1] (stalled) 2021-05-28T13:20:48: Iteration 3: -logL=1282288.324, Lambda=1.0e-03, delta=-0.000, step=1.0e+00, max(|grad|)=0.000000 [Prefactor:0] (stalled) 2021-05-28T13:20:49: Iteration 4: -logL=1282288.324, Lambda=1.0e-02, delta=-0.000, step=1.0e+00, max(|grad|)=0.000000 [Index:1] (stalled) 2021-05-28T13:20:49: >Iteration 5: -logL=1282288.324, Lambda=1.0e-01, delta=0.000, step=1.0e+00, max(|grad|)=0.000000 [Prefactor:0] 2021-05-28T13:20:49: 2021-05-28T13:20:49: +============================================+ 2021-05-28T13:20:49: | Maximum likelihood re-optimisation results | 2021-05-28T13:20:49: +============================================+ 2021-05-28T13:20:49: === GOptimizerLM === 2021-05-28T13:20:49: Optimized function value ..: 1282288.324 2021-05-28T13:20:49: Absolute precision ........: 0.005 2021-05-28T13:20:49: Acceptable value decrease .: 2 2021-05-28T13:20:49: Optimization status .......: converged 2021-05-28T13:20:49: Number of parameters ......: 4 2021-05-28T13:20:49: Number of free parameters .: 2 2021-05-28T13:20:49: Number of iterations ......: 5 2021-05-28T13:20:49: Lambda ....................: 0.01 2021-05-28T13:20:49: 2021-05-28T13:20:49: +=========================================+ 2021-05-28T13:20:49: | Maximum likelihood optimisation results | 2021-05-28T13:20:49: +=========================================+ 2021-05-28T13:20:49: === GOptimizerLM === 2021-05-28T13:20:49: Optimized function value ..: 1282288.361 2021-05-28T13:20:49: Absolute precision ........: 0.005 2021-05-28T13:20:49: Acceptable value decrease .: 2 2021-05-28T13:20:49: Optimization status .......: converged 2021-05-28T13:20:49: Number of parameters ......: 10 2021-05-28T13:20:49: Number of free parameters .: 4 2021-05-28T13:20:49: Number of iterations ......: 7 2021-05-28T13:20:49: Lambda ....................: 0.0001 2021-05-28T13:20:49: Maximum log likelihood ....: -1282288.361 2021-05-28T13:20:49: Observed events (Nobs) ...: 253794.000 2021-05-28T13:20:49: Predicted events (Npred) ..: 253854.941 (Nobs - Npred = -60.9413234166277) 2021-05-28T13:20:49: === GModels === 2021-05-28T13:20:49: Number of models ..........: 2 2021-05-28T13:20:49: Number of parameters ......: 10 2021-05-28T13:20:49: === GModelSky === 2021-05-28T13:20:49: Name ......................: src2330 2021-05-28T13:20:49: Instruments ...............: all 2021-05-28T13:20:49: Test Statistic ............: -0.0737513825297356 2021-05-28T13:20:49: Observation identifiers ...: all 2021-05-28T13:20:49: Model type ................: PointSource 2021-05-28T13:20:49: Model components ..........: "PointSource" * "PowerLaw" * "Constant" 2021-05-28T13:20:49: Number of parameters ......: 6 2021-05-28T13:20:49: Number of spatial par's ...: 2 2021-05-28T13:20:49: RA .......................: 184.5492 [-360,360] deg (fixed,scale=1) 2021-05-28T13:20:49: DEC ......................: -5.8137 [-90,90] deg (fixed,scale=1) 2021-05-28T13:20:49: Number of spectral par's ..: 3 2021-05-28T13:20:49: Prefactor ................: 1e-23 +/- 7.01670377521098e-21 [1e-23,1e-13] ph/cm2/s/MeV (free,scale=1e-16,gradient) 2021-05-28T13:20:49: Index ....................: -5 +/- 342.953063165799 [-5,-0] (free,scale=-1,gradient) 2021-05-28T13:20:49: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) 2021-05-28T13:20:49: Number of temporal par's ..: 1 2021-05-28T13:20:49: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2021-05-28T13:20:49: Number of scale par's .....: 0 2021-05-28T13:20:49: === GCTAModelIrfBackground === 2021-05-28T13:20:49: Name ......................: Background 2021-05-28T13:20:49: Instruments ...............: CTA 2021-05-28T13:20:49: Observation identifiers ...: all 2021-05-28T13:20:49: Model type ................: "PowerLaw" * "Constant" 2021-05-28T13:20:49: Number of parameters ......: 4 2021-05-28T13:20:49: Number of spectral par's ..: 3 2021-05-28T13:20:49: Prefactor ................: 0.99706027335371 +/- 0.00384601729826523 [0.001,1000] ph/cm2/s/MeV (free,scale=1,gradient) 2021-05-28T13:20:49: Index ....................: -0.00280988323133382 +/- 0.00220313067859279 [-5,5] (free,scale=1,gradient) 2021-05-28T13:20:49: PivotEnergy ..............: 1000000 [10000,1000000000] MeV (fixed,scale=1000000,gradient) 2021-05-28T13:20:49: Number of temporal par's ..: 1 2021-05-28T13:20:49: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) 2021-05-28T13:20:49: 2021-05-28T13:20:49: +==============+ 2021-05-28T13:20:49: | Save results | 2021-05-28T13:20:49: +==============+ 2021-05-28T13:20:49: Model definition file .....: ./outmodel1.xml 2021-05-28T13:20:49: Covariance matrix file ....: NONE 2021-05-28T13:20:49: 2021-05-28T13:20:49: Application "ctlike" terminated after 17 wall clock seconds, consuming 24.3724 seconds of CPU time.