Hi all ROOT users,
I’m trying to fit (migrad) an histogram with the sum of 3 TF1. Despite the fitting function (red) is quite superimposed to the histogram, the fitting procedure fails, returning some NaN.
I tried to use SetParameter to improve the fit, but doesn’t seem to work.
I attach the whole terminal output and the “fitted” histogram (red = fit function, other colors = functions that have to be added to create the red one).
How can I make this working?
Thank you for your help!!
[code]Warning in TCanvas::Constructor: Deleting canvas with same name: c1
Bin Width: 1 Entries: 3.42975e+06
** 9 **SET PRINT 2
** 10 **SET NOGRAD
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 Baseline Mean [p.e] 1.20000e+01 3.60000e+00 no limits
2 Light Yield Mean [p.e./keV] 5.24000e-01 1.57200e-01 no limits
3 Baseline Variance [p.e]^2 1.00000e+00 3.00000e-01 no limits
4 Rel. SPE Variance 1.40000e-01 1.40000e-02 no limits
** 11 **FIX 4
FCN= unknown FROM FIX STATUS=RESET 0 CALLS 0 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS PHYSICAL LIMITS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 Baseline Mean [p.e] 1.20000e+01 3.60000e+00
2 Light Yield Mean [p.e./keV] 5.24000e-01 1.57200e-01
3 Baseline Variance [p.e]^2 1.00000e+00 3.00000e-01
4 Rel. SPE Variance 1.40000e-01 fixed
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
5 Rel. LY Variance 1.50000e-02 1.50000e-03 no limits
** 12 **FIX 5
FCN= unknown FROM FIX STATUS=RESET 0 CALLS 0 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS PHYSICAL LIMITS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 Baseline Mean [p.e] 1.20000e+01 3.60000e+00
2 Light Yield Mean [p.e./keV] 5.24000e-01 1.57200e-01
3 Baseline Variance [p.e]^2 1.00000e+00 3.00000e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
6 Bin Width 1.00000e+00 1.00000e-01 no limits
** 13 **FIX 6
FCN= unknown FROM FIX STATUS=RESET 0 CALLS 0 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS PHYSICAL LIMITS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 Baseline Mean [p.e] 1.20000e+01 3.60000e+00
2 Light Yield Mean [p.e./keV] 5.24000e-01 1.57200e-01
3 Baseline Variance [p.e]^2 1.00000e+00 3.00000e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
7 Constant 4.60000e-01 1.38000e-01 no limits
8 Single C14 Amplitude [Hz] 1.60000e+05 4.80000e+04 no limits
** 14 **SET ERR 1
** 15 **SET PRINT 2
** 16 **MIGRAD 2120 0.0001
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-07
FCN=106388 FROM MIGRAD STATUS=INITIATE 20 CALLS 21 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.20000e+01 3.60000e+00 3.60000e+00 -3.08228e+04
2 Light Yield Mean [p.e./keV] 5.24000e-01 1.57200e-01 1.57200e-01 -2.42928e+06
3 Baseline Variance [p.e]^2 1.00000e+00 3.00000e-01 3.00000e-01 -8.63031e+02
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.60000e-01 1.38000e-01 1.38000e-01 -3.64799e+03
8 Single C14 Amplitude [Hz] 1.60000e+05 4.80000e+04 4.80000e+04 -3.07527e+00
NO ERROR MATRIX
FCN=20412.3 FROM MIGRAD STATUS=PROGRESS 36 CALLS 37 TOTAL
EDM=261815 STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29972e+01 3.60000e+00 9.97179e-01 -1.16015e+04
2 Light Yield Mean [p.e./keV] 5.40249e-01 1.57200e-01 1.62490e-02 -1.11159e+06
3 Baseline Variance [p.e]^2 3.86676e+01 3.00000e-01 3.76676e+01 -3.89530e+02
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 5.30936e-01 1.38000e-01 7.09356e-02 -1.55992e+03
8 Single C14 Amplitude [Hz] 1.65241e+05 4.80000e+04 5.24070e+03 -9.83041e-01
FCN=368.99 FROM MIGRAD STATUS=PROGRESS 189 CALLS 190 TOTAL
EDM=9.59632 STRATEGY= 1 ERROR MATRIX UNCERTAINTY 14.4 per cent
EXT PARAMETER APPROXIMATE STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29567e+01 1.90111e-01 -1.51242e-02 -2.78408e+01
2 Light Yield Mean [p.e./keV] 5.92509e-01 2.19481e-03 2.52779e-04 -4.03298e+02
3 Baseline Variance [p.e]^2 1.82319e+00 4.41293e+00 -5.02781e-01 6.93787e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.99534e-01 1.72265e-02 -1.66956e-03 1.48224e+02
8 Single C14 Amplitude [Hz] 1.59902e+05 9.33374e+02 7.22902e+01 -5.92434e-03
FCN=366.689 FROM MIGRAD STATUS=PROGRESS 399 CALLS 400 TOTAL
EDM=4.28587 STRATEGY= 1 ERROR MATRIX UNCERTAINTY 0.0 per cent
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29021e+01 1.89396e-01 -2.83986e-05 -2.24139e+01
2 Light Yield Mean [p.e./keV] 5.93419e-01 2.16925e-03 4.71998e-07 -3.18705e+02
3 Baseline Variance [p.e]^2 1.41480e-02 4.35976e+00 -9.37873e-04 5.62172e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.93480e-01 1.72388e-02 -3.17645e-06 1.19818e+02
8 Single C14 Amplitude [Hz] 1.60163e+05 9.30380e+02 1.35898e-01 -4.77790e-03
MIGRAD FAILS TO FIND IMPROVEMENT
START COVARIANCE MATRIX CALCULATION.
EIGENVALUES OF SECOND-DERIVATIVE MATRIX:
3.7904e-03 1.8560e-02 5.6530e-01 9.9252e-01 3.4198e+00
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=366.679 FROM HESSE STATUS=OK 39 CALLS 643 TOTAL
EDM=4.2802 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29019e+01 1.83969e-01 1.50980e-04 -2.23887e+01
2 Light Yield Mean [p.e./keV] 5.93423e-01 2.14191e-03 2.13117e-06 -3.18347e+02
3 Baseline Variance [p.e]^2 5.71288e-03 4.35994e+00 6.65101e-04 5.61540e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.93452e-01 1.70353e-02 1.46003e-04 1.19683e+02
8 Single C14 Amplitude [Hz] 1.60164e+05 9.15524e+02 1.09119e+00 -4.77252e-03
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 7 8
1 0.99615 1.000 -0.765 0.420 0.366 -0.984
2 0.99434 -0.765 1.000 -0.891 -0.375 0.699
3 0.98716 0.420 -0.891 1.000 0.272 -0.357
7 0.40256 0.366 -0.375 0.272 1.000 -0.350
8 0.99186 -0.984 0.699 -0.357 -0.350 1.000
FCN=366.678 FROM MIGRAD STATUS=PROGRESS 663 CALLS 664 TOTAL
EDM=4.27912 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29019e+01 1.83969e-01 -2.84649e-05 -2.23859e+01
2 Light Yield Mean [p.e./keV] 5.93423e-01 2.14191e-03 4.72391e-07 -3.18307e+02
3 Baseline Variance [p.e]^2 4.77494e-03 4.35994e+00 -9.37941e-04 5.61470e-01
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.93449e-01 1.70353e-02 -3.17616e-06 1.19668e+02
8 Single C14 Amplitude [Hz] 1.60164e+05 9.15524e+02 1.36224e-01 -4.77193e-03
MIGRAD FAILS TO FIND IMPROVEMENT
MIGRAD TERMINATED WITHOUT CONVERGENCE.
FCN=366.677 FROM MIGRAD STATUS=FAILED 699 CALLS 700 TOTAL
EDM=-nan STRATEGY= 1 ERROR MATRIX UNCERTAINTY 100.0 per cent
EXT PARAMETER APPROXIMATE STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 Baseline Mean [p.e] 1.29018e+01 nan -nan -2.23831e+01
2 Light Yield Mean [p.e./keV] 5.93424e-01 nan -nan -3.18267e+02
3 Baseline Variance [p.e]^2 3.83864e-03 nan -nan -inf
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.93445e-01 nan -nan 1.19653e+02
8 Single C14 Amplitude [Hz] 1.60164e+05 nan -nan -4.77133e-03
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=1
-nan -nan -nan -nan -nan
-nan -nan -nan -nan -nan
-nan -nan -nan -nan -nan
-nan -nan -nan -nan -nan
-nan -nan -nan -nan -nan
ERR MATRIX APPROXIMATE
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 7 8
1 nan -nan -nan -nan -nan -nan
2 nan -nan -nan -nan -nan -nan
3 nan -nan -nan -nan -nan -nan
7 nan -nan -nan -nan -nan -nan
8 -nan -nan -nan -nan -nan -nan
ERR MATRIX APPROXIMATE
Info in : Finished to run MIGRAD - status 4
Warning in : Abnormal termination of minimization.
FCN=366.677 FROM MIGRAD STATUS=FAILED 699 CALLS 700 TOTAL
EDM=-nan STRATEGY= 1 ERROR MATRIX UNCERTAINTY 100.0 per cent
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 Baseline Mean [p.e] 1.29018e+01 nan
2 Light Yield Mean [p.e./keV] 5.93424e-01 nan
3 Baseline Variance [p.e]^2 3.83864e-03 nan
4 Rel. SPE Variance 1.40000e-01 fixed
5 Rel. LY Variance 1.50000e-02 fixed
6 Bin Width 1.00000e+00 fixed
7 Constant 4.93445e-01 nan
8 Single C14 Amplitude [Hz] 1.60164e+05 nan
[/code]
test_fit.pdf (8.07 KB)