Logo ROOT  

Reference Guide
minuit2GausFit.C File Reference
Tutorials » Fit Tutorials

Detailed Description

View in nbviewer Open in SWAN Perform fits with different configurations using Minuit2

Do Fit 1
Do Fit 2
Minuit2Minimizer: Minimize with max-calls 1345 convergence for edm < 0.01 strategy 1
MnSeedGenerator: for initial parameters FCN = 130.4086890832
MnSeedGenerator: Initial state: - FCN = 130.4086890832 Edm = 91.9034 NCalls = 13
VariableMetric: start iterating until Edm is < 2e-05
VariableMetric: Initial state - FCN = 130.4086890832 Edm = 91.9034 NCalls = 13
VariableMetric: Iteration # 0 - FCN = 130.4086890832 Edm = 91.9034 NCalls = 13
VariableMetric: Iteration # 1 - FCN = 69.84374040416 Edm = 1.31429 NCalls = 21
VariableMetric: Iteration # 2 - FCN = 65.47146392405 Edm = 0.251272 NCalls = 29
VariableMetric: Iteration # 3 - FCN = 65.16360304303 Edm = 0.00342912 NCalls = 37
VariableMetric: Iteration # 4 - FCN = 65.15866988104 Edm = 6.85944e-05 NCalls = 45
VariableMetric: Iteration # 5 - FCN = 65.15859769243 Edm = 2.0024e-09 NCalls = 53
VariableMetric: After Hessian - FCN = 65.15859769243 Edm = 1.93774e-09 NCalls = 69
VariableMetric: Iteration # 6 - FCN = 65.15859769243 Edm = 1.93774e-09 NCalls = 69
Number of iterations 7
----------> Iteration 0
FVAL = 130.408689083 Edm = 91.9034258956 Nfcn = 13
Error matrix change = 1
Parameters : p0 = 45.7135 p1 = 0.0276 p2 = 1.03818
----------> Iteration 1
FVAL = 69.8437404042 Edm = 1.314294384 Nfcn = 21
Error matrix change = 0.556547
Parameters : p0 = 39.098 p1 = 0.017263 p2 = 0.97602
----------> Iteration 2
FVAL = 65.4714639241 Edm = 0.251272226796 Nfcn = 29
Error matrix change = 0.506801
Parameters : p0 = 36.0678 p1 = 0.00678987 p2 = 1.02566
----------> Iteration 3
FVAL = 65.163603043 Edm = 0.00342911860885 Nfcn = 37
Error matrix change = 0.25723
Parameters : p0 = 36.3576 p1 = 0.0108748 p2 = 1.03384
----------> Iteration 4
FVAL = 65.158669881 Edm = 6.85944158694e-05 Nfcn = 45
Error matrix change = 0.140183
Parameters : p0 = 36.3021 p1 = 0.0129287 p2 = 1.03428
----------> Iteration 5
FVAL = 65.1585976924 Edm = 2.00240394288e-09 Nfcn = 53
Error matrix change = 0.0876016
Parameters : p0 = 36.3132 p1 = 0.013082 p2 = 1.03413
----------> Iteration 6
FVAL = 65.1585976924 Edm = 1.93773754036e-09 Nfcn = 69
Error matrix change = 0
Parameters : p0 = 36.3132 p1 = 0.013082 p2 = 1.03413
Minuit2Minimizer : Valid minimum - status = 0
FVAL = 65.1585976924265395
Edm = 1.93773754035702457e-09
Nfcn = 69
Constant = 36.3132 +/- 1.52625
Mean = 0.013082 +/- 0.0347499
Sigma = 1.03413 +/- 0.0288039 (limited)
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 34.7869
MnSeedGenerator: for initial parameters FCN = 66.17168175294
MnSeedGenerator: Initial state: - FCN = 66.17168175294 Edm = 0.000115215 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.17168175294 Edm = 0.000115215 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.17168175294 Edm = 0.000115215 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.17156294352 Edm = 1.20757e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.17156294352 Edm = 1.20757e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 34.78692902863 ||
1 || Mean || free || 0.01282787988291 ||0.0348666316134
2 || Sigma || limited || 1.052525488272 ||0.02287829646426
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 34.7967
MnSeedGenerator: for initial parameters FCN = 66.15854209966
MnSeedGenerator: Initial state: - FCN = 66.15854209966 Edm = 2.4499e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15854209966 Edm = 2.4499e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15854209966 Edm = 2.4499e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15851755606 Edm = 2.34582e-09 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15851755606 Edm = 2.34582e-09 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 34.79672794636 ||
1 || Mean || free || 0.01282531657197 ||0.03487161339835
2 || Sigma || limited || 1.052412336455 ||0.02289211306185
------ end Minos scan for low interval for parameter Constant
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 37.8394
MnSeedGenerator: for initial parameters FCN = 66.14680616321
MnSeedGenerator: Initial state: - FCN = 66.14680616321 Edm = 0.000113488 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.14680616321 Edm = 0.000113488 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.14680616321 Edm = 0.000113488 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.14669579129 Edm = 9.98321e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.14669579129 Edm = 9.98321e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 37.839426645 ||
1 || Mean || free || 0.013138315697 ||0.03461465133039
2 || Sigma || limited || 1.016337090554 ||0.02202480214272
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 37.8486
MnSeedGenerator: for initial parameters FCN = 66.15854840966
MnSeedGenerator: Initial state: - FCN = 66.15854840966 Edm = 2.54315e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15854840966 Edm = 2.54315e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15854840966 Edm = 2.54315e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15852305107 Edm = 1.76672e-09 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15852305107 Edm = 1.76672e-09 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 37.84859116424 ||
1 || Mean || free || 0.01314314345261 ||0.03462106703996
2 || Sigma || limited || 1.016226092139 ||0.02198842496187
------ end Minos scan for up interval for parameter Constant
Minos: Lower error for parameter Constant : -1.51651
Minos: Upper error for parameter Constant : 1.53547
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to -0.0216679
MnSeedGenerator: for initial parameters FCN = 66.15867475548
MnSeedGenerator: Initial state: - FCN = 66.15867475548 Edm = 0.00108156 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15867475548 Edm = 0.00108156 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15867475548 Edm = 0.00108156 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15759334219 Edm = 2.34127e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.15759334219 Edm = 2.34127e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 36.27570815353 ||1.52283205643
1 || Mean || fixed || -0.02166791598046 ||
2 || Sigma || limited || 1.034173643509 ||0.02883404011912
------ end Minos scan for low interval for parameter Mean
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to 0.0478319
MnSeedGenerator: for initial parameters FCN = 66.15900963262
MnSeedGenerator: Initial state: - FCN = 66.15900963262 Edm = 0.0010711 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15900963262 Edm = 0.0010711 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15900963262 Edm = 0.0010711 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15794043745 Edm = 1.78958e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.15794043745 Edm = 1.78958e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 36.29027277098 ||1.528570300301
1 || Mean || fixed || 0.04783186596055 ||
2 || Sigma || limited || 1.033627406915 ||0.02876872699949
------ end Minos scan for up interval for parameter Mean
Minos: Lower error for parameter Mean : -0.0347674
Minos: Upper error for parameter Mean : 0.0347613
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.00533
MnSeedGenerator: for initial parameters FCN = 66.17114745115
MnSeedGenerator: Initial state: - FCN = 66.17114745115 Edm = 0.000116035 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.17114745115 Edm = 0.000116035 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.17114745115 Edm = 0.000116035 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.17102737617 Edm = 1.73265e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.17102737617 Edm = 1.73265e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.26272487446 ||1.208458195983
1 || Mean || free || 0.01314134487118 ||0.03476774211277
2 || Sigma || fixed || 1.005325497388 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.0055
MnSeedGenerator: for initial parameters FCN = 66.15855782889
MnSeedGenerator: Initial state: - FCN = 66.15855782889 Edm = 2.60516e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15855782889 Edm = 2.60516e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15855782889 Edm = 2.60516e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15853134832 Edm = 1.47818e-08 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15853134832 Edm = 1.47818e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.25655819797 ||1.219155880895
1 || Mean || free || 0.01313469447119 ||0.03474536599194
2 || Sigma || fixed || 1.005502857103 ||
------ end Minos scan for low interval for parameter Sigma
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.06293
MnSeedGenerator: for initial parameters FCN = 66.14463016447
MnSeedGenerator: Initial state: - FCN = 66.14463016447 Edm = 9.38204e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.14463016447 Edm = 9.38204e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.14463016447 Edm = 9.38204e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.14453941055 Edm = 1.37821e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.14453941055 Edm = 1.37821e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 35.34685842948 ||1.167433401845
1 || Mean || free || 0.01250932915253 ||0.03478840494679
2 || Sigma || fixed || 1.062933208465 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.06314
MnSeedGenerator: for initial parameters FCN = 66.15851550893
MnSeedGenerator: Initial state: - FCN = 66.15851550893 Edm = 3.22677e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15851550893 Edm = 3.22677e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15851550893 Edm = 3.22677e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15848385023 Edm = 2.14124e-08 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15848385023 Edm = 2.14124e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 35.34022746778 ||1.156471787584
1 || Mean || free || 0.01251042186767 ||0.03477993700921
2 || Sigma || fixed || 1.063137834893 ||
------ end Minos scan for up interval for parameter Sigma
Minos: Lower error for parameter Sigma : -0.0286274
Minos: Upper error for parameter Sigma : 0.0290102
Covariance Matrix:
Constant Mean Sigma
Constant 2.3294 0.0002427 -0.027611
Mean 0.0002427 0.0012076 -9.3319e-06
Sigma -0.027611 -9.3319e-06 0.00082969
Correlation Matrix:
Constant Mean Sigma
Constant 1 0.0045761 -0.62806
Mean 0.0045761 1 -0.0093231
Sigma -0.62806 -0.0093231 1
****************************************
Minimizer is Minuit2
Chi2 = 65.1586
NDf = 56
Edm = 1.93774e-09
NCalls = 69
Constant = 36.3132 +/- 1.52625 -1.51651 +1.53547 (Minos)
Mean = 0.013082 +/- 0.0347499 -0.0347674 +0.0347613 (Minos)
Sigma = 1.03413 +/- 0.0288039 -0.0286274 +0.0290102 (Minos) (limited)
Do Fit 3
****************************************
Minimizer is Minuit2
Chi2 = 65.1586
NDf = 56
Edm = 5.63237e-09
NCalls = 68
Constant = 36.3274 +/- 1.52734 -1.51763 +1.53654 (Minos)
Mean = 0.0130805 +/- 0.0347499 -0.0347659 +0.0347628 (Minos)
Sigma = 1.03373 +/- 0.0288151 -0.0286384 +0.0290217 (Minos) (limited)
Do Fit 4
Minuit2Minimizer: Minimize with max-calls 1345 convergence for edm < 0.01 strategy 1
MnSeedGenerator: for initial parameters FCN = 59.37893582261
MnSeedGenerator: Initial state: - FCN = 59.37893582261 Edm = 23.8051 NCalls = 13
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 59.37893582261 Edm = 23.8051 NCalls = 13
VariableMetric: Iteration # 0 - FCN = 59.37893582261 Edm = 23.8051 NCalls = 13
VariableMetric: Iteration # 1 - FCN = 44.94069265714 Edm = 0.773988 NCalls = 21
VariableMetric: Iteration # 2 - FCN = 43.43606117719 Edm = 0.0453628 NCalls = 29
VariableMetric: Iteration # 3 - FCN = 43.39363464345 Edm = 0.000175467 NCalls = 37
VariableMetric: Iteration # 4 - FCN = 43.39345546819 Edm = 9.74238e-08 NCalls = 44
VariableMetric: After Hessian - FCN = 43.39345546819 Edm = 9.98059e-08 NCalls = 60
VariableMetric: Iteration # 5 - FCN = 43.39345546819 Edm = 9.98059e-08 NCalls = 60
Number of iterations 6
----------> Iteration 0
FVAL = 59.3789358226 Edm = 23.8051155917 Nfcn = 13
Error matrix change = 1
Parameters : p0 = 45.7135 p1 = 0.0276 p2 = 1.03818
----------> Iteration 1
FVAL = 44.9406926571 Edm = 0.773987920349 Nfcn = 21
Error matrix change = 0.580853
Parameters : p0 = 40.6144 p1 = 0.0276007 p2 = 1.00093
----------> Iteration 2
FVAL = 43.4360611772 Edm = 0.0453627798765 Nfcn = 29
Error matrix change = 0.40441
Parameters : p0 = 38.2485 p1 = 0.0276002 p2 = 1.03476
----------> Iteration 3
FVAL = 43.3936346434 Edm = 0.000175467155603 Nfcn = 37
Error matrix change = 0.237745
Parameters : p0 = 38.4512 p1 = 0.0276009 p2 = 1.0378
----------> Iteration 4
FVAL = 43.3934554682 Edm = 9.7423803633e-08 Nfcn = 44
Error matrix change = 0.128267
Parameters : p0 = 38.427 p1 = 0.027601 p2 = 1.03819
----------> Iteration 5
FVAL = 43.3934554682 Edm = 9.98059014815e-08 Nfcn = 60
Error matrix change = 0
Parameters : p0 = 38.427 p1 = 0.027601 p2 = 1.03819
Minuit2Minimizer : Valid minimum - status = 0
FVAL = 43.3934554681902327
Edm = 9.98059014815227973e-08
Nfcn = 60
Constant = 38.427 +/- 1.48837
Mean = 0.027601 +/- 0.032831
Sigma = 1.03819 +/- 0.0232194 (limited)
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 36.9387
MnSeedGenerator: for initial parameters FCN = 43.90873925582
MnSeedGenerator: Initial state: - FCN = 43.90873925582 Edm = 0.000132178 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.90873925582 Edm = 0.000132178 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.90873925582 Edm = 0.000132178 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.90860183103 Edm = 2.09782e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.90860183103 Edm = 2.09782e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 36.93866931645 ||
1 || Mean || free || 0.02760135158809 ||0.0328310221178
2 || Sigma || limited || 1.051945513458 ||0.01945999630157
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 36.9607
MnSeedGenerator: for initial parameters FCN = 43.89326997499
MnSeedGenerator: Initial state: - FCN = 43.89326997499 Edm = 5.0917e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89326997499 Edm = 5.0917e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89326997499 Edm = 5.0917e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89321904908 Edm = 6.05599e-13 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89321904908 Edm = 6.05599e-13 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 36.96071319015 ||
1 || Mean || free || 0.02760128717335 ||0.03283102271151
2 || Sigma || limited || 1.051749143987 ||0.01945893847383
------ end Minos scan for low interval for parameter Constant
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 39.9154
MnSeedGenerator: for initial parameters FCN = 43.87938545703
MnSeedGenerator: Initial state: - FCN = 43.87938545703 Edm = 0.000138937 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.87938545703 Edm = 0.000138937 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.87938545703 Edm = 0.000138937 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.87925218963 Edm = 2.295e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.87925218963 Edm = 2.295e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 39.91541762539 ||
1 || Mean || free || 0.0276008466576 ||0.03283099763553
2 || Sigma || limited || 1.025134476821 ||0.01847912948869
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 39.937
MnSeedGenerator: for initial parameters FCN = 43.89329471521
MnSeedGenerator: Initial state: - FCN = 43.89329471521 Edm = 5.69203e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89329471521 Edm = 5.69203e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89329471521 Edm = 5.69203e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89323778415 Edm = 7.78066e-13 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89323778415 Edm = 7.78066e-13 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 39.9370187331 ||
1 || Mean || free || 0.0276007763412 ||0.03283099824302
2 || Sigma || limited || 1.024937319007 ||0.01847821131437
------ end Minos scan for up interval for parameter Constant
Minos: Lower error for parameter Constant : -1.46667
Minos: Upper error for parameter Constant : 1.51031
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to -0.00523004
MnSeedGenerator: for initial parameters FCN = 43.89345820803
MnSeedGenerator: Initial state: - FCN = 43.89345820803 Edm = 0.000261584 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89345820803 Edm = 0.000261584 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89345820803 Edm = 0.000261584 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89319680439 Edm = 6.02092e-11 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.89319680439 Edm = 6.02092e-11 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 38.40723377277 ||1.488068262532
1 || Mean || fixed || -0.005230041413473 ||
2 || Sigma || limited || 1.038717822636 ||0.02322032420128
------ end Minos scan for low interval for parameter Mean
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to 0.060432
MnSeedGenerator: for initial parameters FCN = 43.89345269077
MnSeedGenerator: Initial state: - FCN = 43.89345269077 Edm = 0.000256919 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89345269077 Edm = 0.000256919 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89345269077 Edm = 0.000256919 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.8931959431 Edm = 5.58583e-11 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.8931959431 Edm = 5.58583e-11 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 38.40721102864 ||1.488070519217
1 || Mean || fixed || 0.06043198080746 ||
2 || Sigma || limited || 1.038718498907 ||0.02322057576887
------ end Minos scan for up interval for parameter Mean
Minos: Lower error for parameter Mean : -0.0328395
Minos: Upper error for parameter Mean : 0.0328395
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.01497
MnSeedGenerator: for initial parameters FCN = 43.91322394171
MnSeedGenerator: Initial state: - FCN = 43.91322394171 Edm = 0.000116281 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.91322394171 Edm = 0.000116281 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.91322394171 Edm = 0.000116281 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.91310259725 Edm = 2.22312e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.91310259725 Edm = 2.22312e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 39.30500437194 ||1.242641962681
1 || Mean || free || 0.02760073943742 ||0.03283088908748
2 || Sigma || fixed || 1.014970026362 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.01541
MnSeedGenerator: for initial parameters FCN = 43.8931845953
MnSeedGenerator: Initial state: - FCN = 43.8931845953 Edm = 8.60941e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.8931845953 Edm = 8.60941e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.8931845953 Edm = 8.60941e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89309848173 Edm = 1.01991e-11 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89309848173 Edm = 1.01991e-11 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 39.28869836486 ||1.242675512086
1 || Mean || free || 0.02760072551055 ||0.03283088944825
2 || Sigma || fixed || 1.015413201704 ||
------ end Minos scan for low interval for parameter Sigma
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.06141
MnSeedGenerator: for initial parameters FCN = 43.87506262924
MnSeedGenerator: Initial state: - FCN = 43.87506262924 Edm = 0.000127983 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.87506262924 Edm = 0.000127983 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.87506262924 Edm = 0.000127983 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.87494052583 Edm = 2.67712e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.87494052583 Edm = 2.67712e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.58705454386 ||1.188298763707
1 || Mean || free || 0.02760170741526 ||0.03283107916862
2 || Sigma || fixed || 1.061408778615 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.06185
MnSeedGenerator: for initial parameters FCN = 43.89321154572
MnSeedGenerator: Initial state: - FCN = 43.89321154572 Edm = 9.64242e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89321154572 Edm = 9.64242e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89321154572 Edm = 9.64242e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89311509665 Edm = 6.67633e-12 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89311509665 Edm = 6.67633e-12 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.57055265832 ||1.18834606973
1 || Mean || free || 0.02760170074746 ||0.03283107928904
2 || Sigma || fixed || 1.061851004318 ||
------ end Minos scan for up interval for parameter Sigma
Minos: Lower error for parameter Sigma : -0.0227841
Minos: Upper error for parameter Sigma : 0.0236699
Covariance Matrix:
Constant Mean Sigma
Constant 2.2153 -3.264e-06 -0.019956
Mean -3.264e-06 0.0010779 8.9196e-08
Sigma -0.019956 8.9196e-08 0.00053915
Correlation Matrix:
Constant Mean Sigma
Constant 1 -6.6797e-05 -0.57743
Mean -6.6797e-05 1 0.00011701
Sigma -0.57743 0.00011701 1
****************************************
Minimizer is Minuit2
MinFCN = 43.3935
Chi2 = 74.4576
NDf = 97
Edm = 9.98059e-08
NCalls = 60
Constant = 38.427 +/- 1.48837 -1.46667 +1.51031 (Minos)
Mean = 0.027601 +/- 0.032831 -0.0328395 +0.0328395 (Minos)
Sigma = 1.03819 +/- 0.0232194 -0.0227841 +0.0236699 (Minos) (limited)
Do Fit 1
Do Fit 2
Minuit2Minimizer: Minimize with max-calls 1345 convergence for edm < 0.01 strategy 1
MnSeedGenerator: for initial parameters FCN = 130.4086890832
MnSeedGenerator: Initial state: - FCN = 130.4086890832 Edm = 91.9034 NCalls = 13
Fumili: start iterating until Edm is < 1e-06
Fumili: Initial state - FCN = 130.4086890832 Edm = 91.9034 NCalls = 13
Fumili: Iteration # 2 - FCN = 79.16955375104 Edm = 11.2614 NCalls = 14
Fumili: Iteration # 3 - FCN = 65.15899551326 Edm = 0.000440916 NCalls = 15
Fumili: Iteration # 4 - FCN = 65.15859944286 Edm = 1.63332e-06 NCalls = 16
Fumili: Iteration # 5 - FCN = 65.15859769852 Edm = 7.48888e-09 NCalls = 17
Fumili: After Hessian - FCN = 65.15859769852 Edm = 8.03809e-09 NCalls = 45
Number of iterations 6
----------> Iteration 0
FVAL = 130.408689083 Edm = 91.9034257196 Nfcn = 13
Error matrix change = 1
Parameters : p0 = 45.7135 p1 = 0.0276 p2 = 1.03818
----------> Iteration 1
FVAL = 79.169553751 Edm = 11.2613962604 Nfcn = 14
Error matrix change = 0.727775
Parameters : p0 = 36.1708 p1 = 0.0126892 p2 = 0.949066
----------> Iteration 2
FVAL = 65.1589955133 Edm = 0.000440916194411 Nfcn = 15
Error matrix change = 0.429675
Parameters : p0 = 36.2965 p1 = 0.0129513 p2 = 1.03469
----------> Iteration 3
FVAL = 65.1585994429 Edm = 1.63332385536e-06 Nfcn = 16
Error matrix change = 0.215296
Parameters : p0 = 36.3119 p1 = 0.0130821 p2 = 1.03417
----------> Iteration 4
FVAL = 65.1585976985 Edm = 7.48887590173e-09 Nfcn = 17
Error matrix change = 0.107683
Parameters : p0 = 36.3131 p1 = 0.0130818 p2 = 1.03413
----------> Iteration 5
FVAL = 65.1585976985 Edm = 8.03809330501e-09 Nfcn = 45
Error matrix change = 0
Parameters : p0 = 36.3131 p1 = 0.0130818 p2 = 1.03413
Minuit2Minimizer : Valid minimum - status = 0
FVAL = 65.1585976985221009
Edm = 8.03809330501340492e-09
Nfcn = 45
Constant = 36.3131 +/- 1.52625
Mean = 0.0130818 +/- 0.0347499
Sigma = 1.03413 +/- 0.0288039 (limited)
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 34.7868
MnSeedGenerator: for initial parameters FCN = 66.17179656804
MnSeedGenerator: Initial state: - FCN = 66.17179656804 Edm = 0.000114641 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.17179656804 Edm = 0.000114641 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.17179656804 Edm = 0.000114641 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.17167834955 Edm = 1.20177e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.17167834955 Edm = 1.20177e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 34.78684262717 ||
1 || Mean || free || 0.0128278444398 ||0.03486670438982
2 || Sigma || limited || 1.052526566199 ||0.02287829819503
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 34.7967
MnSeedGenerator: for initial parameters FCN = 66.15854261731
MnSeedGenerator: Initial state: - FCN = 66.15854261731 Edm = 2.49664e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15854261731 Edm = 2.49664e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15854261731 Edm = 2.49664e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15851760566 Edm = 2.38471e-09 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15851760566 Edm = 2.38471e-09 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 34.79672790901 ||
1 || Mean || free || 0.01282533083682 ||0.03487153809973
2 || Sigma || limited || 1.0524123382 ||0.02289210304275
------ end Minos scan for low interval for parameter Constant
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 37.8393
MnSeedGenerator: for initial parameters FCN = 66.14669137862
MnSeedGenerator: Initial state: - FCN = 66.14669137862 Edm = 0.00011292 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.14669137862 Edm = 0.00011292 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.14669137862 Edm = 0.00011292 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.1465815596 Edm = 9.93495e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.1465815596 Edm = 9.93495e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 37.83933786873 ||
1 || Mean || free || 0.01313833187114 ||0.03461457277685
2 || Sigma || limited || 1.01633809602 ||0.02202484016116
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 37.8486
MnSeedGenerator: for initial parameters FCN = 66.15854884808
MnSeedGenerator: Initial state: - FCN = 66.15854884808 Edm = 2.58951e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15854884808 Edm = 2.58951e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15854884808 Edm = 2.58951e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15852302745 Edm = 1.80676e-09 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15852302745 Edm = 1.80676e-09 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 37.84859114596 ||
1 || Mean || free || 0.01314312871966 ||0.03462088442543
2 || Sigma || limited || 1.01622608852 ||0.02198839735727
------ end Minos scan for up interval for parameter Constant
Minos: Lower error for parameter Constant : -1.51642
Minos: Upper error for parameter Constant : 1.53556
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to -0.0216681
MnSeedGenerator: for initial parameters FCN = 66.15868704027
MnSeedGenerator: Initial state: - FCN = 66.15868704027 Edm = 0.00108064 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15868704027 Edm = 0.00108064 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15868704027 Edm = 0.00108064 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15760652641 Edm = 2.3681e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.15760652641 Edm = 2.3681e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 36.27570694486 ||1.522818712605
1 || Mean || fixed || -0.02166814548753 ||
2 || Sigma || limited || 1.034173669298 ||0.0288344808871
------ end Minos scan for low interval for parameter Mean
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to 0.0478317
MnSeedGenerator: for initial parameters FCN = 66.15899735635
MnSeedGenerator: Initial state: - FCN = 66.15899735635 Edm = 0.00107023 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15899735635 Edm = 0.00107023 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15899735635 Edm = 0.00107023 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15792904156 Edm = 1.8096e-08 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.15792904156 Edm = 1.8096e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 36.29027381276 ||1.528579941891
1 || Mean || fixed || 0.04783166767322 ||
2 || Sigma || limited || 1.033627389593 ||0.02876838027371
------ end Minos scan for up interval for parameter Mean
Minos: Lower error for parameter Mean : -0.0347671
Minos: Upper error for parameter Mean : 0.0347615
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.00533
MnSeedGenerator: for initial parameters FCN = 66.17103425966
MnSeedGenerator: Initial state: - FCN = 66.17103425966 Edm = 0.00011552 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.17103425966 Edm = 0.00011552 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.17103425966 Edm = 0.00011552 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.1709147234 Edm = 1.72099e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.1709147234 Edm = 1.72099e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.26267108435 ||1.208409247736
1 || Mean || free || 0.01314135307985 ||0.03476823300593
2 || Sigma || fixed || 1.005327091355 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.0055
MnSeedGenerator: for initial parameters FCN = 66.15855728432
MnSeedGenerator: Initial state: - FCN = 66.15855728432 Edm = 2.56055e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15855728432 Edm = 2.56055e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15855728432 Edm = 2.56055e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15853125533 Edm = 1.46111e-08 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15853125533 Edm = 1.46111e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.25655770252 ||1.219155680359
1 || Mean || free || 0.01313471695752 ||0.03474567944057
2 || Sigma || fixed || 1.005502858425 ||
------ end Minos scan for low interval for parameter Sigma
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.06293
MnSeedGenerator: for initial parameters FCN = 66.14474254323
MnSeedGenerator: Initial state: - FCN = 66.14474254323 Edm = 9.3372e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.14474254323 Edm = 9.3372e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.14474254323 Edm = 9.3372e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.14465221649 Edm = 1.36711e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 66.14465221649 Edm = 1.36711e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 35.34680406136 ||1.167488843672
1 || Mean || free || 0.01250928549283 ||0.03478793021982
2 || Sigma || fixed || 1.062934869512 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.06314
MnSeedGenerator: for initial parameters FCN = 66.15851512063
MnSeedGenerator: Initial state: - FCN = 66.15851512063 Edm = 3.17346e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 1e-05
VariableMetric: Initial state - FCN = 66.15851512063 Edm = 3.17346e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 66.15851512063 Edm = 3.17346e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 66.15848398823 Edm = 2.12152e-08 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 66.15848398823 Edm = 2.12152e-08 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 35.34022778435 ||1.156471795901
1 || Mean || free || 0.01251040017838 ||0.03477934397344
2 || Sigma || fixed || 1.063137836914 ||
------ end Minos scan for up interval for parameter Sigma
Minos: Lower error for parameter Sigma : -0.0286291
Minos: Upper error for parameter Sigma : 0.0290085
Covariance Matrix:
Constant Mean Sigma
Constant 2.3294 0.00024274 -0.027611
Mean 0.00024274 0.0012076 -9.3326e-06
Sigma -0.027611 -9.3326e-06 0.00082969
Correlation Matrix:
Constant Mean Sigma
Constant 1 0.0045768 -0.62806
Mean 0.0045768 1 -0.0093238
Sigma -0.62806 -0.0093238 1
****************************************
Minimizer is Minuit2 / Fumili
Chi2 = 65.1586
NDf = 56
Edm = 8.03809e-09
NCalls = 45
Constant = 36.3131 +/- 1.52625 -1.51642 +1.53556 (Minos)
Mean = 0.0130818 +/- 0.0347499 -0.0347671 +0.0347615 (Minos)
Sigma = 1.03413 +/- 0.0288039 -0.0286291 +0.0290085 (Minos) (limited)
Do Fit 3
****************************************
Minimizer is Minuit2 / Fumili
Chi2 = 65.1586
NDf = 56
Edm = 1.52369e-08
NCalls = 45
Constant = 36.3272 +/- 1.52734 -1.51745 +1.53671 (Minos)
Mean = 0.0130818 +/- 0.0347499 -0.0347671 +0.0347615 (Minos)
Sigma = 1.03373 +/- 0.0288151 -0.0286415 +0.0290186 (Minos) (limited)
Do Fit 4
Minuit2Minimizer: Minimize with max-calls 1345 convergence for edm < 0.01 strategy 1
MnSeedGenerator: for initial parameters FCN = 59.37893582261
MnSeedGenerator: Initial state: - FCN = 59.37893582261 Edm = 23.8051 NCalls = 13
Fumili: start iterating until Edm is < 5e-07
Fumili: Initial state - FCN = 59.37893582261 Edm = 23.8051 NCalls = 13
Fumili: Iteration # 2 - FCN = 52.23535907444 Edm = 15.7973 NCalls = 14
Fumili: Iteration # 3 - FCN = 44.59114050271 Edm = 0.678155 NCalls = 15
Fumili: Iteration # 4 - FCN = 43.49291484693 Edm = 0.0638834 NCalls = 16
Fumili: Iteration # 5 - FCN = 43.40537058229 Edm = 0.00696015 NCalls = 17
Fumili: Iteration # 6 - FCN = 43.39560346139 Edm = 0.00120472 NCalls = 18
Fumili: Iteration # 7 - FCN = 43.39386966545 Edm = 0.000231498 NCalls = 19
Fumili: Iteration # 8 - FCN = 43.39353563705 Edm = 4.48954e-05 NCalls = 20
Fumili: Iteration # 9 - FCN = 43.39347091481 Edm = 8.70149e-06 NCalls = 21
Fumili: Iteration # 10 - FCN = 43.39345837804 Edm = 1.68518e-06 NCalls = 22
Fumili: Iteration # 11 - FCN = 43.39345595089 Edm = 3.26218e-07 NCalls = 23
Fumili: After Hessian - FCN = 43.39345595089 Edm = 5.82496e-07 NCalls = 51
Fumili: Iteration # 13 - FCN = 43.39345536839 Edm = 1.95619e-15 NCalls = 52
Fumili: After Hessian - FCN = 43.39345536839 Edm = 1.96796e-15 NCalls = 84
Number of iterations 14
----------> Iteration 0
FVAL = 59.3789358226 Edm = 23.8051156246 Nfcn = 13
Error matrix change = 1
Parameters : p0 = 45.7135 p1 = 0.0276 p2 = 1.03818
----------> Iteration 1
FVAL = 52.2353590744 Edm = 15.797316836 Nfcn = 14
Error matrix change = 0.946152
Parameters : p0 = 37.0455 p1 = 0.0276012 p2 = 0.975225
----------> Iteration 2
FVAL = 44.5911405027 Edm = 0.678155218069 Nfcn = 15
Error matrix change = 0.516813
Parameters : p0 = 37.6923 p1 = 0.0135425 p2 = 1.01848
----------> Iteration 3
FVAL = 43.4929148469 Edm = 0.0638834337408 Nfcn = 16
Error matrix change = 0.269485
Parameters : p0 = 38.0009 p1 = 0.0171462 p2 = 1.0397
----------> Iteration 4
FVAL = 43.4053705823 Edm = 0.00696014762982 Nfcn = 17
Error matrix change = 0.139355
Parameters : p0 = 38.2218 p1 = 0.0263423 p2 = 1.03895
----------> Iteration 5
FVAL = 43.3956034614 Edm = 0.00120472221679 Nfcn = 18
Error matrix change = 0.072946
Parameters : p0 = 38.3316 p1 = 0.0271707 p2 = 1.03891
----------> Iteration 6
FVAL = 43.3938696655 Edm = 0.000231498054964 Nfcn = 19
Error matrix change = 0.0380601
Parameters : p0 = 38.3845 p1 = 0.0275055 p2 = 1.03848
----------> Iteration 7
FVAL = 43.3935356371 Edm = 4.48954181602e-05 Nfcn = 20
Error matrix change = 0.0197548
Parameters : p0 = 38.4079 p1 = 0.0275531 p2 = 1.03834
----------> Iteration 8
FVAL = 43.3934709148 Edm = 8.70149166209e-06 Nfcn = 21
Error matrix change = 0.010197
Parameters : p0 = 38.4183 p1 = 0.0275828 p2 = 1.03826
----------> Iteration 9
FVAL = 43.393458378 Edm = 1.68517866991e-06 Nfcn = 22
Error matrix change = 0.00523948
Parameters : p0 = 38.4229 p1 = 0.0275925 p2 = 1.03822
----------> Iteration 10
FVAL = 43.3934559509 Edm = 3.26217613616e-07 Nfcn = 23
Error matrix change = 0.00268169
Parameters : p0 = 38.4249 p1 = 0.0275974 p2 = 1.03821
----------> Iteration 11
FVAL = 43.3934559509 Edm = 5.82496472019e-07 Nfcn = 51
Error matrix change = 0
Parameters : p0 = 38.4249 p1 = 0.0275974 p2 = 1.03821
----------> Iteration 12
FVAL = 43.3934553684 Edm = 1.95619311042e-15 Nfcn = 52
Error matrix change = 0.37759
Parameters : p0 = 38.4265 p1 = 0.027601 p2 = 1.0382
----------> Iteration 13
FVAL = 43.3934553684 Edm = 1.96795731974e-15 Nfcn = 84
Error matrix change = 0
Parameters : p0 = 38.4265 p1 = 0.027601 p2 = 1.0382
Minuit2Minimizer : Valid minimum - status = 0
FVAL = 43.3934553683858653
Edm = 1.96795731974180762e-15
Nfcn = 84
Constant = 38.4265 +/- 1.48835
Mean = 0.027601 +/- 0.0328313
Sigma = 1.0382 +/- 0.0232198 (limited)
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 36.9381
MnSeedGenerator: for initial parameters FCN = 43.90913752769
MnSeedGenerator: Initial state: - FCN = 43.90913752769 Edm = 0.000128912 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.90913752769 Edm = 0.000128912 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.90913752769 Edm = 0.000128912 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.90900349914 Edm = 2.04554e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.90900349914 Edm = 2.04554e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 36.93809827799 ||
1 || Mean || free || 0.02760135116835 ||0.03283131266545
2 || Sigma || limited || 1.051951097523 ||0.01946025022673
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 36.9607
MnSeedGenerator: for initial parameters FCN = 43.89327284463
MnSeedGenerator: Initial state: - FCN = 43.89327284463 Edm = 5.38546e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89327284463 Edm = 5.38546e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89327284463 Edm = 5.38546e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89321898134 Edm = 5.72157e-13 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89321898134 Edm = 5.72157e-13 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 36.96071328799 ||
1 || Mean || free || 0.02760128835915 ||0.03283131323332
2 || Sigma || limited || 1.051749140492 ||0.01945900753087
------ end Minos scan for low interval for parameter Constant
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #0 : Constant using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 0
MnFunctionCross: Run Migrad fixing parameters : Constant #0 to 39.9148
MnSeedGenerator: for initial parameters FCN = 43.87898700326
MnSeedGenerator: Initial state: - FCN = 43.87898700326 Edm = 0.000135588 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.87898700326 Edm = 0.000135588 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.87898700326 Edm = 0.000135588 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.87885694851 Edm = 2.2403e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.87885694851 Edm = 2.2403e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 39.91480282851 ||
1 || Mean || free || 0.02760084623686 ||0.03283128780438
2 || Sigma || limited || 1.02513959399 ||0.01847936065682
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 39.937
MnSeedGenerator: for initial parameters FCN = 43.89329728076
MnSeedGenerator: Initial state: - FCN = 43.89329728076 Edm = 5.99101e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89329728076 Edm = 5.99101e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89329728076 Edm = 5.99101e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89323736006 Edm = 7.27813e-13 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89323736006 Edm = 7.27813e-13 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || fixed || 39.93701808256 ||
1 || Mean || free || 0.02760077769483 ||0.0328312883912
2 || Sigma || limited || 1.024937322128 ||0.01847827295471
------ end Minos scan for up interval for parameter Constant
Minos: Lower error for parameter Constant : -1.46608
Minos: Upper error for parameter Constant : 1.5109
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to -0.00523033
MnSeedGenerator: for initial parameters FCN = 43.89345811834
MnSeedGenerator: Initial state: - FCN = 43.89345811834 Edm = 0.000252607 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89345811834 Edm = 0.000252607 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89345811834 Edm = 0.000252607 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89320568772 Edm = 5.81412e-11 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.89320568772 Edm = 5.81412e-11 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 38.40723387762 ||1.488055620279
1 || Mean || fixed || -0.005230333357774 ||
2 || Sigma || limited || 1.038717834184 ||0.02322048266436
------ end Minos scan for low interval for parameter Mean
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #1 : Mean using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 1
MnFunctionCross: Run Migrad fixing parameters : Mean #1 to 0.0604323
MnSeedGenerator: for initial parameters FCN = 43.89345258537
MnSeedGenerator: Initial state: - FCN = 43.89345258537 Edm = 0.000248021 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89345258537 Edm = 0.000248021 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89345258537 Edm = 0.000248021 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89320473128 Edm = 5.37527e-11 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.89320473128 Edm = 5.37527e-11 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 38.40721113403 ||1.488058006041
1 || Mean || fixed || 0.06043226962547 ||
2 || Sigma || limited || 1.038718510468 ||0.02322073426213
------ end Minos scan for up interval for parameter Mean
Minos: Lower error for parameter Mean : -0.0328395
Minos: Upper error for parameter Mean : 0.0328395
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS LOWER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of lower Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.01498
MnSeedGenerator: for initial parameters FCN = 43.91282758618
MnSeedGenerator: Initial state: - FCN = 43.91282758618 Edm = 0.000119483 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.91282758618 Edm = 0.000119483 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.91282758618 Edm = 0.000119483 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.91270290067 Edm = 2.28497e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.91270290067 Edm = 2.28497e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 39.30465313231 ||1.242626811845
1 || Mean || free || 0.02760073968046 ||0.0328311825415
2 || Sigma || fixed || 1.014978794059 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.01541
MnSeedGenerator: for initial parameters FCN = 43.89318034584
MnSeedGenerator: Initial state: - FCN = 43.89318034584 Edm = 8.24299e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89318034584 Edm = 8.24299e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89318034584 Edm = 8.24299e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89309789618 Edm = 9.64739e-12 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89309789618 Edm = 9.64739e-12 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 39.28869808509 ||1.242669954814
1 || Mean || free || 0.02760072050436 ||0.03283118304179
2 || Sigma || fixed || 1.01541321481 ||
------ end Minos scan for low interval for parameter Sigma
******************************************************************************************************
Minuit2Minimizer::GetMinosError - Run MINOS UPPER error for parameter #2 : Sigma using max-calls 1345, tolerance 0.01
--------- MnMinos -------
Determination of upper Minos error for parameter 2
MnFunctionCross: Run Migrad fixing parameters : Sigma #2 to 1.06142
MnSeedGenerator: for initial parameters FCN = 43.87545635866
MnSeedGenerator: Initial state: - FCN = 43.87545635866 Edm = 0.000131341 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.87545635866 Edm = 0.000131341 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.87545635866 Edm = 0.000131341 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.87533105126 Edm = 2.74651e-07 NCalls = 14
MnFunctionCross: Result after Migrad - FCN = 43.87533105126 Edm = 2.74651e-07 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.58672607686 ||1.18828439594
1 || Mean || free || 0.02760170767428 ||0.03283136776382
2 || Sigma || fixed || 1.061418364269 ||
MnFunctionCross: Run Migrad again (#2) : - parameter 0 fixed to 1.06185
MnSeedGenerator: for initial parameters FCN = 43.89320778165
MnSeedGenerator: Initial state: - FCN = 43.89320778165 Edm = 9.26238e-05 NCalls = 9
VariableMetric: start iterating until Edm is < 5e-06
VariableMetric: Initial state - FCN = 43.89320778165 Edm = 9.26238e-05 NCalls = 9
VariableMetric: Iteration # 0 - FCN = 43.89320778165 Edm = 9.26238e-05 NCalls = 9
VariableMetric: Iteration # 1 - FCN = 43.89311513272 Edm = 6.53023e-12 NCalls = 14
MnFunctionCross: Result after 2nd Migrad - FCN = 43.89311513272 Edm = 6.53023e-12 NCalls = 14
# ext. || Name || type || Value || Error +/-
0 || Constant || free || 37.57055285115 ||1.18834088042
1 || Mean || free || 0.02760169702021 ||0.03283136795789
2 || Sigma || fixed || 1.061851005188 ||
------ end Minos scan for up interval for parameter Sigma
Minos: Lower error for parameter Sigma : -0.0227933
Minos: Upper error for parameter Sigma : 0.0236607
Covariance Matrix:
Constant Mean Sigma
Constant 2.2152 -3.2649e-06 -0.019956
Mean -3.2649e-06 0.0010779 8.922e-08
Sigma -0.019956 8.922e-08 0.00053917
Correlation Matrix:
Constant Mean Sigma
Constant 1 -6.6816e-05 -0.57743
Mean -6.6816e-05 1 0.00011703
Sigma -0.57743 0.00011703 1
****************************************
Minimizer is Minuit2 / Fumili
MinFCN = 43.3935
Chi2 = 74.4574
NDf = 97
Edm = 1.96796e-15
NCalls = 84
Constant = 38.4265 +/- 1.48835 -1.46608 +1.5109 (Minos)
Mean = 0.027601 +/- 0.0328313 -0.0328395 +0.0328395 (Minos)
Sigma = 1.0382 +/- 0.0232198 -0.0227933 +0.0236607 (Minos) (limited)
#include "TH1.h"
#include "TF1.h"
#include "TCanvas.h"
#include "TRandom3.h"
#include "TVirtualFitter.h"
#include "TPaveLabel.h"
#include "TStyle.h"
#include <iostream>
#include <string>
void testGausFit( std::string type = "Minuit2", int n = 1000) {
gRandom = new TRandom3();
TVirtualFitter::SetDefaultFitter(type.c_str() );
std::string name;
name = "h1_" + type;
TH1D * h1 = new TH1D(name.c_str(),"Chi2 Fit",100, -5, 5. );
name = "h2_" + type;
TH1D * h2 = new TH1D(name.c_str(),"Chi2 Fit with Minos Error",100, -5, 5. );
name = "h3_" + type;
TH1D * h3 = new TH1D(name.c_str(),"Chi2 Fit with Integral and Minos",100, -5, 5. );
name = "h4_" + type;
TH1D * h4 = new TH1D(name.c_str(),"Likelihood Fit with Minos Error",100, -5, 5. );
gStyle->SetOptStat(1111111);
gStyle->SetOptFit(1111111);
for (int i = 0; i < n; ++i) {
double x = gRandom->Gaus(0,1);
h1->Fill( x );
h2->Fill( x );
h3->Fill( x );
h4->Fill( x );
}
std::string cname = type + "Canvas" ;
std::string ctitle = type + " Gaussian Fit" ;
TCanvas *c1 = new TCanvas(cname.c_str(),cname.c_str(),10,10,900,900);
c1->Divide(2,2);
c1->cd(1);
cout << "\nDo Fit 1\n";
h1->Fit("gaus","Q");
h1->Draw();
c1->cd(2);
cout << "\nDo Fit 2\n";
h2->Fit("gaus","VE");
h2->Draw();
c1->cd(3);
cout << "\nDo Fit 3\n";
h3->Fit("gaus","IE");
h3->Draw();
c1->cd(4);
cout << "\nDo Fit 4\n";
h4->Fit("gaus","VLE");
h4->Draw();
}
void minuit2GausFit() {
int n = 1000;
testGausFit("Minuit2",n);
testGausFit("Fumili2",n);
}
name
char name[80]
Definition: TGX11.cxx:109
type
int type
Definition: TGX11.cxx:120
gRandom
R__EXTERN TRandom * gRandom
Definition: TRandom.h:62
gStyle
R__EXTERN TStyle * gStyle
Definition: TStyle.h:410
TCanvas
The Canvas class.
Definition: TCanvas.h:27
TH1D
1-D histogram with a double per channel (see TH1 documentation)}
Definition: TH1.h:614
TH1::Fit
virtual TFitResultPtr Fit(const char *formula, Option_t *option="", Option_t *goption="", Double_t xmin=0, Double_t xmax=0)
Fit histogram with function fname.
Definition: TH1.cxx:3808
TH1::Fill
virtual Int_t Fill(Double_t x)
Increment bin with abscissa X by 1.
Definition: TH1.cxx:3275
TH1::Draw
virtual void Draw(Option_t *option="")
Draw this histogram with options.
Definition: TH1.cxx:2998
TRandom3
Random number generator class based on M.
Definition: TRandom3.h:27
TRandom::Gaus
virtual Double_t Gaus(Double_t mean=0, Double_t sigma=1)
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigm...
Definition: TRandom.cxx:263
TStyle::SetOptStat
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
Definition: TStyle.cxx:1590
TStyle::SetOptFit
void SetOptFit(Int_t fit=1)
The type of information about fit parameters printed in the histogram statistics box can be selected ...
Definition: TStyle.cxx:1542
TVirtualFitter::SetDefaultFitter
static void SetDefaultFitter(const char *name="")
static: set name of default fitter
Definition: TVirtualFitter.cxx:245
c1
return c1
Definition: legend1.C:41
x
Double_t x[n]
Definition: legend1.C:17
n
const Int_t n
Definition: legend1.C:16
h1
TH1F * h1
Definition: legend1.C:5
Author
Lorenzo Moneta

Definition in file minuit2GausFit.C.