-
Notifications
You must be signed in to change notification settings - Fork 1
Description
Hi guys,
I am receiving the same error, no matter if I am using MINUIT or NEWMINUIT. Is there anyone who can help me with this? I am inserting the values used and the error in calculating the likelihood. I have tried extending the time limits, reducing the energy limits. But the same issue arises. I've also tried CALDB and P8R3_SOURCE_V3 response functions as well. Also, tried reducing the search radius to 10.
SOURCE: CTA 102
Equatorial coordinates (degrees) | (338.152,11.7308)
Time range (MET) | (504144004,505872005)
Time range (Gregorian) | (2016-12-23 00:00:00,2017-01-12 00:00:00)
Energy range (MeV) | (100,300000)
Search radius (degrees) | 30
=======================================================
gtlike refit=yes plot=yes sfile=CTA102output_model.xml
Statistic to use (BINNED|UNBINNED) [UNBINNED]
Spacecraft file[spacecraft.fits]
Event file[CTA_102_region_filtered_gti.fits]
Unbinned exposure map[CTA_102_unbin_expmap.fits]
Exposure hypercube file[CTA_102_ltcube.fits]
Source model file[CTA102input_model.xml]
Response functions to use[CALDB]
Optimizer (DRMNFB|NEWMINUIT|MINUIT|DRMNGB|LBFGS) [NEWMINUIT] MINUIT
** 1 **SET PRINT .000
** 2 **SET NOWARN
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 'Prefactor ' 10.000 1.0000 .10000E-02 1000.0
2 'Index ' -2.0000 1.0000 -5.0000 -1.0000
3 'Prefactor ' 1.0000 1.0000 .0000 10.000
4 'Normalizat' 1.0000 1.0000 .10000E-01 10.000
** 3 **SET ERR .5000
** 4 **SET GRAD 1.000
** 5 **MINIMIZE 700.0 2.000
MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
FCN= 171995.0 FROM MIGRAD STATUS=CONVERGED 100 CALLS 101 TOTAL
EDM= .55E-03 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Prefactor 103.59 1.3555 .10063E-01 -1.2435
2 Index -1.9834 .70021E-02 .24293E-01 -4.6675
3 Prefactor 1.0025 2.5273 .46267 -.26492E-03
4 Normalizat 4.0377 .62158E-01 .11255 -.89310E-02
ERR DEF= .500
Final values:
Prefactor = 103.589
Index = -1.98343
Prefactor = 1.00254
Normalizat = 4.03771
** 6 **HESSE
FCN= 171995.0 FROM HESSE STATUS=OK 25 CALLS 126 TOTAL
EDM= .60E-03 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 Prefactor 103.59 1.4143 .54785E-02 -22979.
2 Index -1.9834 .72898E-02 .75327E-02 31594.
3 Prefactor 1.0025 2.5273 .92534E-01 -.92645
4 Normalizat 4.0377 .62506E-01 .49742E-03 4699.6
ERR DEF= .500
Minuit fit quality: 3 estimated distance: 0.000599406
Minuit parameter uncertainties:
1 1.41434
2 0.00728983
3 3.00338
4 0.0625079
Computing TS values for each source (3 total)
.Caught N10optimizers11OutOfBoundsE at the top level: Attempt to set the value outside of existing bounds.
==================================
gtlike refit=yes plot=yes sfile=CTA102output_model.xml
Statistic to use (BINNED|UNBINNED) [UNBINNED]
Spacecraft file[none] spacecraft.fits
Event file[none] CTA_102_region_filtered_gti.fits
Unbinned exposure map[none] CTA_102_unbin_expmap.fits
Exposure hypercube file[none] CTA_102_ltcube.fits
Source model file[] CTA102input_model.xml
Response functions to use[CALDB]
Optimizer (DRMNFB|NEWMINUIT|MINUIT|DRMNGB|LBFGS) [MINUIT] NEWMINUIT
Minuit did successfully converge.
of function calls: 42
minimum function Value: 171994.9872681
minimum edm: 4.75373262634e-06
minimum internal state vector: LAVector parameters:
-2.226172565269
3558.891291822
-0.9271294544361
47.31875575351
minimum internal covariance matrix: LASymMatrix parameters:
4.4614549e-05 -3.1660721e-05 -1.6288658e-06 -2.7433362e-05
-3.1660721e-05 3.6857014e-05 5.3607321e-06 1.5163935e-05
-1.6288658e-06 5.3607321e-06 0.42813035 -4.8388454e-06
-2.7433362e-05 1.5163935e-05 -4.8388454e-06 0.00035378235
ext. || Name || type || Value || Error +/-
0 ||0_Prefactor || limited || 103.5915549928 ||1.439244864446
1 || 1_Index || limited || -1.983442899588 ||0.007393908851233
2 ||2_Prefactor || limited || 1.000497345649 ||1.339317644674
3 ||3_Normalization || limited || 4.037793057884 ||0.06517445024886
of function calls: 93
function Value: 171994.9872681
expected distance to the Minimum (edm): 4.295608770915e-06
external parameters:
ext. || Name || type || Value || Error +/-
0 ||0_Prefactor || limited || 103.5915549928 ||1.410577529742
1 || 1_Index || limited || -1.983442899588 ||0.007265448214781
2 ||2_Prefactor || limited || 1.000497345649 ||5.205923884213
3 ||3_Normalization || limited || 4.037793057884 ||0.06257240777728
covariance matrix status: 3
covariance matrix:
MnUserCovariance:
1.98974 -0.00793135 2.03474e-06 -0.0227087
-0.00793135 5.27871e-05 1.95012e-09 6.37048e-05
2.03474e-06 1.95012e-09 18.0079 -1.52107e-06
-0.0227087 6.37048e-05 -1.52107e-06 0.00391552
MnUserCovariance Parameter correlations:
1 -0.7739 3.39922e-07 -0.257276
-0.7739 1 6.32507e-08 0.140124
3.39922e-07 6.32507e-08 1 -5.72827e-06
-0.257276 0.140124 -5.72827e-06 1
global correlation coefficients :
MnGlobalCorrelationCoeff:
0.788363
0.776304
5.84737e-06
0.273613
Computing TS values for each source (3 total)
.Caught N10optimizers11OutOfBoundsE at the top level: Attempt to set the value outside of existing bounds.
Please help asap.
Thanks in advance.
-Ayush