Processing CheckIfMinosChangesParrErr.C... RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University All rights reserved, please read http://roofit.sourceforge.net/license.txt M O D E L P A R A M E T E R S: ---------------------------------------------- massSigMeanTruth 60 massSigStdDevTruth 4 massBkgLamTruth -0.025 massRangeMin 0 massRangeMax 150 vggSigMeanTruth 120 vggSigStdDevTruth 5 vggBkgMeanTruth 100 vggBkgStdDevTruth 10 vggRangeMin 60 vggRangeMax 140 vexexSigLamTruth -0.0125 vexexBkgLamTruth -0.01 vexexRangeMin 0 vexexRangeMax 350 ---------------------------------------------- ModelBuilder::bldMassSig() ModelBuilder::bldMassbkg ModelBuilder::bldMass() ModelBuilder::bldVggSig ModelBuilder::bldVggBkg ModelBuilder::bldVgg ModelBuilder::bld_bkg_sig_full RooAddPdf::pdfFull[ nSigFit * pdfSigProd + nBkgFit * pdfBkgProd ] = 0.122684 RooRealVar::nSigFit = 20000 L(-INF - +INF) RooRealVar::nBkgFit = 80000 L(-INF - +INF) [#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization [#1] INFO:Minization -- The following expressions will be evaluated in cache-and-track mode: (pdfMassSig,pdfMassBkg) ********** ** 1 **SET PRINT 1 ********** ********** ** 2 **SET NOGRAD ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 massBkgLamFit -2.50000e-02 1.00000e+00 no limits 2 massSigMeanFit 6.00000e+01 1.00000e+00 no limits 3 massSigStdDevFit 4.00000e+00 1.00000e+00 no limits 4 nBkgFit 8.00000e+04 1.00000e+00 no limits 5 nSigFit 2.00000e+04 1.00000e+00 no limits ********** ** 3 **SET ERR 0.5 ********** ********** ** 4 **SET PRINT 1 ********** ********** ** 5 **SET STR 1 ********** NOW USING STRATEGY 1: TRY TO BALANCE SPEED AGAINST RELIABILITY ********** ** 6 **MIGRAD 2500 1 ********** FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4. START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03 FCN=-599824 FROM MIGRAD STATUS=INITIATE 24 CALLS 25 TOTAL EDM= unknown STRATEGY= 1 NO ERROR MATRIX EXT PARAMETER CURRENT GUESS STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 massBkgLamFit -2.50000e-02 1.00000e+00 1.00000e+00 9.54293e+03 2 massSigMeanFit 6.00000e+01 1.00000e+00 1.00000e+00 -1.55558e+01 3 massSigStdDevFit 4.00000e+00 1.00000e+00 1.00000e+00 -1.49824e+01 4 nBkgFit 8.00000e+04 1.00000e+00 1.00000e+00 -4.57909e-03 5 nSigFit 2.00000e+04 1.00000e+00 1.00000e+00 -6.48367e-03 ERR DEF= 0.5 MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. COVARIANCE MATRIX CALCULATED SUCCESSFULLY FCN=-599826 FROM MIGRAD STATUS=CONVERGED 117 CALLS 118 TOTAL EDM=0.000136139 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 massBkgLamFit -2.50277e-02 1.13798e-04 5.95889e-05 -2.74732e+01 2 massSigMeanFit 6.00220e+01 3.95307e-02 2.10741e-02 -4.09282e-02 3 massSigStdDevFit 4.02236e+00 3.59633e-02 1.79501e-02 2.98251e-02 4 nBkgFit 8.03107e+04 3.05358e+02 1.56861e+02 -1.19252e-05 5 nSigFit 2.01853e+04 1.82057e+02 8.96868e+01 4.69423e-05 ERR DEF= 0.5 EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5 1.295e-08 -7.622e-08 -5.783e-07 3.556e-03 -3.555e-03 -7.622e-08 1.563e-03 -9.354e-05 3.791e-01 -3.798e-01 -5.783e-07 -9.354e-05 1.293e-03 -2.131e+00 2.135e+00 3.556e-03 3.791e-01 -2.131e+00 9.324e+04 -1.293e+04 -3.555e-03 -3.798e-01 2.135e+00 -1.293e+04 3.314e+04 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 1 2 3 4 5 1 0.20347 1.000 -0.017 -0.141 0.102 -0.172 2 0.08138 -0.017 1.000 -0.066 0.031 -0.053 3 0.36007 -0.141 -0.066 1.000 -0.194 0.326 4 0.26951 0.102 0.031 -0.194 1.000 -0.233 5 0.38729 -0.172 -0.053 0.326 -0.233 1.000 ********** ** 7 **SET ERR 0.5 ********** ********** ** 8 **SET PRINT 1 ********** ********** ** 9 **HESSE 2500 ********** COVARIANCE MATRIX CALCULATED SUCCESSFULLY FCN=-599826 FROM HESSE STATUS=OK 31 CALLS 149 TOTAL EDM=0.000136639 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 massBkgLamFit -2.50277e-02 1.13829e-04 2.38356e-06 -2.50277e-02 2 massSigMeanFit 6.00220e+01 3.95367e-02 4.21482e-03 6.00220e+01 3 massSigStdDevFit 4.02236e+00 3.60038e-02 7.18005e-04 4.02236e+00 4 nBkgFit 8.03107e+04 3.05499e+02 6.27443e+00 8.03107e+04 5 nSigFit 2.01853e+04 1.82233e+02 3.58747e+00 2.01853e+04 ERR DEF= 0.5 EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5 1.296e-08 -7.557e-08 -5.842e-07 3.583e-03 -3.583e-03 -7.557e-08 1.563e-03 -9.733e-05 3.818e-01 -3.819e-01 -5.842e-07 -9.733e-05 1.296e-03 -2.153e+00 2.153e+00 3.583e-03 3.818e-01 -2.153e+00 9.333e+04 -1.302e+04 -3.583e-03 -3.819e-01 2.153e+00 -1.302e+04 3.321e+04 PARAMETER CORRELATION COEFFICIENTS NO. GLOBAL 1 2 3 4 5 1 0.20474 1.000 -0.017 -0.143 0.103 -0.173 2 0.08321 -0.017 1.000 -0.068 0.032 -0.053 3 0.36275 -0.143 -0.068 1.000 -0.196 0.328 4 0.27109 0.103 0.032 -0.196 1.000 -0.234 5 0.38940 -0.173 -0.053 0.328 -0.234 1.000 [#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization [#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization [#1] INFO:Minization -- The following expressions have been identified as constant and will be precalculated and cached: (pdfMassSig,pdfMassBkg) [#1] INFO:Fitting -- RooAbsPdf::fitTo(pdfMassSum) Calculating sum-of-weights-squared correction matrix for covariance matrix [#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization [#1] INFO:InputArguments -- Printing Yields RooArgList:: = (nSigFit,nBkgFit) [#1] INFO:InputArguments -- yield in pdf: nSigFit 20183.7 [#1] INFO:InputArguments -- yield in pdf: nBkgFit 80312 [#1] INFO:Eval -- Checking Likelihood normalization: [#1] INFO:Eval -- Yield of specie Sum of Row in Matrix Norm [#1] INFO:Eval -- 20183.7 20183.6 1 [#1] INFO:Eval -- 80312 80311.9 1 [#1] INFO:Eval -- Calculating sWeight RooDataSet::pdfFullData[mass,vgg,nSigFit_sw,L_nSigFit,nBkgFit_sw,L_nBkgFit] = 100496 entries RooExtendPdf::pdfExtVggSig[ pdf=pdfVggSig n=nSigFit ] = 0.000335463 ----- FIT WITH MINOS AND SUMW2ERROR ENABLED ----- [#0] WARNING:InputArguments -- RooAbsPdf::fitTo(pdfExtVggSig) WARNING: sum-of-weights correction does not apply to MINOS errors [#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization ********** ** 1 **SET PRINT 0 ********** ********** ** 2 **SET NOGRAD ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 nSigFit 2.01837e+04 1.70445e+02 no limits 2 vggSigMeanFit 1.20000e+02 1.00000e+00 no limits 3 vggSigStdDevFit 5.00000e+00 1.00000e+00 no limits ********** ** 3 **SET ERR 0.5 ********** ********** ** 4 **SET PRINT 0 ********** ********** ** 5 **SET STR 1 ********** ********** ** 6 **MIGRAD 1500 1 ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN=-119104 FROM MIGRAD STATUS=CONVERGED 44 CALLS 45 TOTAL EDM=4.0998e-05 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 nSigFit 2.01837e+04 1.42069e+02 3.38368e+01 9.24192e-07 2 vggSigMeanFit 1.20016e+02 3.46261e-02 8.24569e-03 -4.01338e-03 3 vggSigStdDevFit 4.91826e+00 2.45256e-02 5.83378e-03 2.60970e-01 ERR DEF= 0.5 ********** ** 7 **SET ERR 0.5 ********** ********** ** 8 **SET PRINT 0 ********** ********** ** 9 **HESSE 1500 ********** FCN=-119104 FROM HESSE STATUS=OK 16 CALLS 61 TOTAL EDM=4.11166e-05 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 nSigFit 2.01837e+04 1.42070e+02 1.35347e+00 2.01837e+04 2 vggSigMeanFit 1.20016e+02 3.46261e-02 1.64914e-03 1.20016e+02 3 vggSigStdDevFit 4.91826e+00 2.45256e-02 1.16676e-03 4.91826e+00 ERR DEF= 0.5 [#1] INFO:Fitting -- RooAbsPdf::fitTo(pdfExtVggSig) Calculating sum-of-weights-squared correction matrix for covariance matrix ********** ** 10 **SET ERR 0.5 ********** ********** ** 11 **SET PRINT 0 ********** ********** ** 12 **HESSE 1500 ********** FCN=-94226 FROM HESSE STATUS=OK 24 CALLS 85 TOTAL EDM=53813.2 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 nSigFit 2.01837e+04 1.18416e+02 2.76108e+01 2.01837e+04 2 vggSigMeanFit 1.20016e+02 3.37407e-02 3.36424e-02 1.20016e+02 3 vggSigStdDevFit 4.91826e+00 7.96455e-03 2.38018e-02 4.91826e+00 ERR DEF= 0.5 ********** ** 13 **MINOS 1500 1 ********** FCN=-119048 FROM MINOS STATUS=NEW MINIMU 99 CALLS 184 TOTAL EDM= unknown STRATEGY= 1 NO ERROR MATRIX EXT PARAMETER PARABOLIC MINOS ERRORS NO. NAME VALUE ERROR NEGATIVE POSITIVE 1 nSigFit 2.17232e+04 1.18416e+02 2 vggSigMeanFit 1.20016e+02 3.46222e-02 3 vggSigStdDevFit 4.91811e+00 2.45095e-02 ERR DEF= 0.5 NEW MINIMUM FOUND. GO BACK TO MINIMIZATION STEP. ================================================= V V V VVVVVVV VVVVV VVV V MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN=-119104 FROM MIGRAD STATUS=CONVERGED 46 CALLS 230 TOTAL EDM=1.12805e-07 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 nSigFit 2.01837e+04 1.42069e+02 3.37821e+01 -1.66898e-06 2 vggSigMeanFit 1.20016e+02 3.46250e-02 8.25161e-03 -5.44488e-06 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 5.84433e-03 9.69967e-03 ERR DEF= 0.5 FCN=-119104 FROM MINOS STATUS=SUCCESSFUL 18 CALLS 248 TOTAL EDM=1.12805e-07 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER PARABOLIC MINOS ERRORS NO. NAME VALUE ERROR NEGATIVE POSITIVE 1 nSigFit 2.01837e+04 1.42069e+02 -1.41702e+02 1.42436e+02 2 vggSigMeanFit 1.20016e+02 3.46250e-02 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 ERR DEF= 0.5 ********** ** 14 **MINOS 1500 2 ********** FCN=-119104 FROM MINOS STATUS=SUCCESSFUL 22 CALLS 270 TOTAL EDM=1.12805e-07 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER PARABOLIC MINOS ERRORS NO. NAME VALUE ERROR NEGATIVE POSITIVE 1 nSigFit 2.01837e+04 1.42069e+02 -1.41702e+02 1.42436e+02 2 vggSigMeanFit 1.20016e+02 3.46250e-02 -3.46251e-02 3.46256e-02 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 ERR DEF= 0.5 ********** ** 15 **MINOS 1500 3 ********** FCN=-119104 FROM MINOS STATUS=SUCCESSFUL 22 CALLS 292 TOTAL EDM=1.12805e-07 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER PARABOLIC MINOS ERRORS NO. NAME VALUE ERROR NEGATIVE POSITIVE 1 nSigFit 2.01837e+04 1.42069e+02 -1.41702e+02 1.42436e+02 2 vggSigMeanFit 1.20016e+02 3.46250e-02 -3.46251e-02 3.46256e-02 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 -2.44265e-02 2.46207e-02 ERR DEF= 0.5 [#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization ----- FIT WITH MINOS DISABLED AND SUMW2ERROR ENABLED ----- [#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization ********** ** 16 **SET PRINT 0 ********** ********** ** 17 **SET NOGRAD ********** PARAMETER DEFINITIONS: NO. NAME VALUE STEP SIZE LIMITS 1 nSigFit 2.01837e+04 1.18416e+02 no limits 2 vggSigMeanFit 1.20016e+02 3.37407e-02 no limits 3 vggSigStdDevFit 4.91826e+00 7.96455e-03 no limits ********** ** 18 **SET ERR 0.5 ********** ********** ** 19 **SET PRINT 0 ********** ********** ** 20 **SET STR 1 ********** ********** ** 21 **MIGRAD 1500 1 ********** MIGRAD MINIMIZATION HAS CONVERGED. MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX. FCN=-119104 FROM MIGRAD STATUS=CONVERGED 36 CALLS 37 TOTAL EDM=6.49285e-08 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 nSigFit 2.01837e+04 1.42069e+02 3.38368e+01 6.49565e-07 2 vggSigMeanFit 1.20016e+02 3.46250e-02 8.25186e-03 2.32165e-04 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 5.84478e-03 9.67866e-03 ERR DEF= 0.5 ********** ** 22 **SET ERR 0.5 ********** ********** ** 23 **SET PRINT 0 ********** ********** ** 24 **HESSE 1500 ********** FCN=-119104 FROM HESSE STATUS=OK 16 CALLS 53 TOTAL EDM=6.51364e-08 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 nSigFit 2.01837e+04 1.42070e+02 1.35347e+00 2.01837e+04 2 vggSigMeanFit 1.20016e+02 3.46250e-02 1.65037e-03 1.20016e+02 3 vggSigStdDevFit 4.91811e+00 2.45237e-02 2.33791e-04 4.91811e+00 ERR DEF= 0.5 [#1] INFO:Fitting -- RooAbsPdf::fitTo(pdfExtVggSig) Calculating sum-of-weights-squared correction matrix for covariance matrix ********** ** 25 **SET ERR 0.5 ********** ********** ** 26 **SET PRINT 0 ********** ********** ** 27 **HESSE 1500 ********** FCN=-94221.3 FROM HESSE STATUS=OK 24 CALLS 77 TOTAL EDM=53718.2 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER INTERNAL INTERNAL NO. NAME VALUE ERROR STEP SIZE VALUE 1 nSigFit 2.01837e+04 1.18416e+02 2.76108e+01 2.01837e+04 2 vggSigMeanFit 1.20016e+02 3.38118e-02 3.36676e-02 1.20016e+02 3 vggSigStdDevFit 4.91811e+00 7.98124e-03 4.76934e-03 4.91811e+00 ERR DEF= 0.5 [#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization ----- COMPARE THE FIT RESULTS ----- Have look at the error of vggSigStdDevFit in both cases WITH MINOS AND SUMW2ERROR ENABLED: RooFitResult: minimized FCN value: -94226, estimated distance to minimum: 53813.2 covariance matrix quality: Unknown, matrix was externally provided Status : MINIMIZE=0 HESSE=0 HESSE=0 MINOS=0 Floating Parameter FinalValue +/- Error -------------------- -------------------------- nSigFit 2.0184e+04 +/- 1.18e+02 vggSigMeanFit 1.2002e+02 +/- 3.37e-02 vggSigStdDevFit 4.9183e+00 +/- 7.96e-03 Covariance for nSigFit: 29052.3 Sqrt(cov): 170.447 Covariance for vggSigMeanFit: 0.00171763 Sqrt(cov): 0.0414443 Covariance for vggSigStdDevFit: 0.00779129 Sqrt(cov): 0.0882683 => Parameter error and Sqrt of covariance matrix element do NOT agree WITH MINOS DISABLED AND SUMW2ERROR ENABLED: RooFitResult: minimized FCN value: -94221.3, estimated distance to minimum: 53718.2 covariance matrix quality: Unknown, matrix was externally provided Status : MINIMIZE=0 HESSE=0 HESSE=0 Floating Parameter FinalValue +/- Error -------------------- -------------------------- nSigFit 2.0184e+04 +/- 1.70e+02 vggSigMeanFit 1.2002e+02 +/- 4.14e-02 vggSigStdDevFit 4.9181e+00 +/- 8.83e-02 Covariance for nSigFit: 29052.3 Sqrt(cov): 170.447 Covariance for vggSigMeanFit: 0.00171748 Sqrt(cov): 0.0414425 Covariance for vggSigStdDevFit: 0.00778953 Sqrt(cov): 0.0882583 => Paramter error and Sqrt of covariance matrix element do agree