fermi-lat / Likelihood

BSD 3-Clause "New" or "Revised" License
2 stars 1 forks source link

ERROR: Caught N10optimizers11OutOfBoundsE at the top level: Attempt to set the value outside of existing bounds #110

Open nyctophile-aayu opened 2 years ago

nyctophile-aayu commented 2 years ago

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