Logo ROOT  

Reference Guide
 
Loading...
Searching...
No Matches
rf506_msgservice.C File Reference
Tutorials » RooFit Tutorials

Detailed Description

View in nbviewer Open in SWAN Organisation and simultaneous fits: tuning and customizing the RooFit message logging facility

␛[1mRooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby␛[0m
Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University
All rights reserved, please read http://roofit.sourceforge.net/license.txt
[#0] WARNING:InputArguments -- The parameter 's' with range [-10, 10] of the RooGaussian 'g' exceeds the safe range of (0, inf). Advise to limit its range.
Active Message streams
[0] MinLevel = PROGRESS Topic = Generation Minimization Plotting Fitting Integration LinkStateMgmt Eval Caching Optimization ObjectHandling InputArguments Tracing Contents DataHandling NumericIntegration FastEvaluations
[1] MinLevel = INFO Topic = Minimization Plotting Fitting Eval Caching ObjectHandling InputArguments DataHandling NumericIntegration
[2] MinLevel = INFO Topic = HistFactory
Active Message streams
[0] MinLevel = PROGRESS Topic = Generation Minimization Plotting Fitting Integration LinkStateMgmt Eval Caching Optimization ObjectHandling InputArguments Tracing Contents DataHandling NumericIntegration FastEvaluations
[1] MinLevel = INFO Topic = Minimization Plotting Fitting Eval Caching ObjectHandling InputArguments DataHandling NumericIntegration
[2] MinLevel = INFO Topic = HistFactory
[#1] INFO:Integration -- RooRealIntegral::ctor(g_Int[x]) Constructing integral of function g over observables(x) with normalization () with range identifier <none>
[#1] INFO:Integration -- g: Observable x is suitable for analytical integration (if supported by p.d.f)
[#1] INFO:Integration -- g: Function integrated observables (x) internally with code 1
[#1] INFO:Integration -- g: Observables (x) are analytically integrated with code 1
RooRealIntegral::g_Int[x][ Int gd[Ana](x) ] = 2.50663
Active Message streams
[0] MinLevel = PROGRESS Topic = Generation Minimization Plotting Fitting Integration LinkStateMgmt Eval Caching Optimization ObjectHandling InputArguments Tracing Contents DataHandling NumericIntegration FastEvaluations
[1] MinLevel = INFO Topic = Minimization Plotting Fitting Integration Eval Caching ObjectHandling InputArguments DataHandling NumericIntegration
[2] MinLevel = INFO Topic = HistFactory
[#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization
[#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (g,p)
[#0] WARNING:Minimization -- RooAbsMinimizerFcn::synchronize: WARNING: no initial error estimate available for f: using 0.1
[#0] WARNING:Minimization -- RooAbsMinimizerFcn::synchronize: WARNING: no initial error estimate available for m: using 2
[#0] WARNING:Minimization -- RooAbsMinimizerFcn::synchronize: WARNING: no initial error estimate available for p0: using 0.005
[#0] WARNING:Minimization -- RooAbsMinimizerFcn::synchronize: WARNING: no initial error estimate available for s: using 2
**********
** 1 **SET PRINT 1
**********
**********
** 2 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 f 5.00000e-01 1.00000e-01 0.00000e+00 1.00000e+00
2 m 0.00000e+00 2.00000e+00 -1.00000e+01 1.00000e+01
3 p0 1.00000e-02 5.00000e-03 0.00000e+00 1.00000e+00
4 s 1.00000e+00 2.00000e+00 -1.00000e+01 1.00000e+01
**********
** 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 2000 1
**********
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
prevFCN = 31.75882326 START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03
f=0.501, p0=0.01,
prevFCN = 31.76758673 f=0.499,
prevFCN = 31.75008506 f=0.5, m=0.02014,
prevFCN = 31.72980141 m=-0.02014,
prevFCN = 31.78749743 m=0.004235,
prevFCN = 31.75274724 m=-0.004235,
prevFCN = 31.76488389 m=0, p0=0.01005,
prevFCN = 31.75954549 p0=0.009948,
prevFCN = 31.75810351 p0=0.01013,
prevFCN = 31.76062402 p0=0.009872,
prevFCN = 31.75703788 p0=0.01, s=1.02,
prevFCN = 31.67915486 s=0.9799,
prevFCN = 31.84068856 s=1.002,
prevFCN = 31.75075447 s=0.998,
prevFCN = 31.76691405 m=0.004235, s=1,
prevFCN = 31.75274724 m=0.02118,
prevFCN = 31.72829189 m=0.06353,
prevFCN = 31.66614321 m=0.1906,
prevFCN = 31.47244052 m=0.5715,
prevFCN = 30.87227362 m=1.707,
prevFCN = 29.88047134 m=2.064,
prevFCN = 29.97604268 f=0.5008, m=1.707,
prevFCN = 29.88564172 f=0.4992,
prevFCN = 29.87531802 f=0.5, m=1.711,
prevFCN = 29.8803659 m=1.703,
prevFCN = 29.88060433 m=1.71,
prevFCN = 29.88039214 m=1.704,
prevFCN = 29.88056503 m=1.707, p0=0.01013,
prevFCN = 29.88274199 p0=0.009875,
prevFCN = 29.87821833 p0=0.01, s=1.002,
prevFCN = 29.87987611 s=0.9984,
prevFCN = 29.88107011 s=1.003,
prevFCN = 29.87926952 s=0.9967,
prevFCN = 29.88168765 FCN=29.8805 FROM MIGRAD STATUS=INITIATE 33 CALLS 34 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 f 5.00000e-01 1.00000e-01 0.00000e+00 3.38641e+00
2 m 1.70692e+00 2.00000e+00 1.71532e-01 -2.90086e-01
3 p0 1.00000e-02 5.00000e-03 0.00000e+00 1.80182e+00
4 s 1.00000e+00 2.00000e+00 0.00000e+00 -3.63523e+00
ERR DEF= 0.5
f=0.2774, m=1.724, p0=0.000388, s=1.276,
prevFCN = 28.97288883 f=0.2912, m=1.723, p0=0.0006269, s=1.257,
prevFCN = 28.97893369 f=0.278, m=1.724, p0=0.000388, s=1.276,
prevFCN = 28.97325909 f=0.2768,
prevFCN = 28.97252778 f=0.2774, m=1.727,
prevFCN = 28.97321652 m=1.722,
prevFCN = 28.97256819 m=1.724, p0=0.0004103,
prevFCN = 28.97321664 p0=0.0003662,
prevFCN = 28.97257031 p0=0.000388, s=1.279,
prevFCN = 28.97432524 s=1.273,
prevFCN = 28.97145596 s=1.282,
prevFCN = 28.97575564 s=1.269,
prevFCN = 28.97003607 f=0.2452, m=1.663, p0=5.498e-07, s=0.9876,
prevFCN = 28.86089971 f=0.2274, m=1.628, p0=0.0001029, s=0.8218,
prevFCN = 28.87647974 f=0.241, m=1.655, p0=3.204e-06, s=0.9491,
prevFCN = 28.85685477 f=0.2418,
prevFCN = 28.85692699 f=0.2403,
prevFCN = 28.85679737 f=0.241, m=1.659,
prevFCN = 28.85436009 m=1.651,
prevFCN = 28.8593705 m=1.655, p0=1.215e-06,
prevFCN = 28.85682558 p0=6.138e-06,
prevFCN = 28.85689784 p0=3.204e-06, s=0.9556,
prevFCN = 28.85804685 s=0.9426,
prevFCN = 28.85571545 s=0.9524,
prevFCN = 28.85746237 s=0.9458,
prevFCN = 28.85626117 f=0.2171, m=1.753, p0=0.0002379, s=0.6877,
prevFCN = 28.6750677 f=0.1311, m=2.143, p0=0.004887, s=-0.3611, RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (31.8407) to force MIGRAD to back out of this region. Error log follows.
Parameter values: f=0.131146 m=2.14309 p0=0.00488743 s=-0.361096
RooNLLVar::nll_model_modelData[ paramSet=(f,m,p0,s) ]
function value is NAN @ paramSet=(f = 0.131146,m = 2.14309,p0 = 0.00488743,s = -0.361096)
RooAddPdf::model[ f * g + [%] * p ]
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 7.09176e-64/-0.905133,p = 0.980378/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 6.44518e-55/-0.905133,p = 1.03835/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.703736/-0.905133,p = 1.01195/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.847231/-0.905133,p = 1.01149/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 4.52496e-09/-0.905133,p = 0.999534/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 1.58333e-80/-0.905133,p = 0.976642/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 8.36671e-221/-0.905133,p = 0.95429/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 6.49026e-20/-0.905133,p = 0.993885/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 3.6989e-92/-0.905133,p = 0.97426/20), !coefficients=(f = 0.131146)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.868573/-0.905133,p = 1.00954/20), !coefficients=(f = 0.131146)
RooGaussian::g[ x=x mean=m sigma=s ]
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-4.01476, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=7.84746, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=2.44579, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=2.35102, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-0.0953319, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-4.77917, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-9.35254, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-1.25126, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=-5.2665, mean=m=2.14309, sigma=s=-0.361096
p.d.f normalization integral is zero or negative: -0.905133 @ x=x=1.9514, mean=m=2.14309, sigma=s=-0.361096
prevFCN = 122.3539561 f=0.2285, m=1.706, p0=7.868e-05, s=0.8134,
prevFCN = 28.78647388 f=0.1125, m=2.24, p0=0.006973, s=-0.6231, RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (31.8407) to force MIGRAD to back out of this region. Error log follows.
Parameter values: f=0.112499 m=2.24008 p0=0.00697314 s=-0.623113
RooNLLVar::nll_model_modelData[ paramSet=(f,m,p0,s) ]
function value is NAN @ paramSet=(f = 0.112499,m = 2.24008,p0 = 0.00697314,s = -0.623113)
RooAddPdf::model[ f * g + [%] * p ]
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 1.31738e-22/-1.56191,p = 0.972004/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 2.60031e-18/-1.56191,p = 1.05472/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.946964/-1.56191,p = 1.01705/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.984276/-1.56191,p = 1.01639/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.000890568/-1.56191,p = 0.999335/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 2.78564e-28/-1.56191,p = 0.966674/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 6.92557e-76/-1.56191,p = 0.934783/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 1.5234e-07/-1.56191,p = 0.991275/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 3.06144e-32/-1.56191,p = 0.963276/20), !coefficients=(f = 0.112499)
getLogVal() top-level p.d.f evaluates to NaN @ !refCoefNorm=(), !pdfs=(g = 0.898243/-1.56191,p = 1.01361/20), !coefficients=(f = 0.112499)
RooGaussian::g[ x=x mean=m sigma=s ]
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-4.01476, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=7.84746, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=2.44579, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=2.35102, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-0.0953319, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-4.77917, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-9.35254, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-1.25126, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=-5.2665, mean=m=2.24008, sigma=s=-0.623113
p.d.f normalization integral is zero or negative: -1.561912 @ x=x=1.9514, mean=m=2.24008, sigma=s=-0.623113
prevFCN = 188.0319041 f=0.1613, m=1.997, p0=0.00245, s=0.03234,
prevFCN = 28.83468394 f=0.1957, m=1.844, p0=0.0007883, s=0.4446,
prevFCN = 28.29265339 f=0.1624, m=1.992, p0=0.002383, s=0.04546,
prevFCN = 28.59402569 f=0.1875, m=1.879, p0=0.001091, s=0.3493,
prevFCN = 28.21262283 f=0.1848, m=1.891, p0=0.001205, s=0.317,
prevFCN = 28.22913773 f=0.1881, m=1.879, p0=0.001091, s=0.3493,
prevFCN = 28.21167987 f=0.1869,
prevFCN = 28.21358338 f=0.1875, m=1.883,
prevFCN = 28.193551 m=1.876,
prevFCN = 28.23178915 m=1.881,
prevFCN = 28.2053519 m=1.878,
prevFCN = 28.21990745 m=1.879, p0=0.001046,
prevFCN = 28.21195074 p0=0.001136,
prevFCN = 28.21330938 p0=0.001091, s=0.3526,
prevFCN = 28.20611524 s=0.346,
prevFCN = 28.21961933 s=0.3499,
prevFCN = 28.2115002 s=0.3488,
prevFCN = 28.21375915 f=0.1722, m=1.965, p0=0.001851, s=0.1732,
prevFCN = 28.7434305 f=0.1862, m=1.887, p0=0.001148, s=0.3343,
prevFCN = 28.20584234 f=0.1867,
prevFCN = 28.20499847 f=0.1856,
prevFCN = 28.20669988 f=0.1862, m=1.888,
prevFCN = 28.19816797 m=1.885,
prevFCN = 28.21353051 m=1.887, p0=0.001103,
prevFCN = 28.20517217 p0=0.001193,
prevFCN = 28.20652621 p0=0.001148, s=0.3349,
prevFCN = 28.20444403 s=0.3338,
prevFCN = 28.20725615 f=0.1818, m=1.915, p0=0.001339, s=0.308,
prevFCN = 28.10982159 f=0.1645, m=2.028, p0=0.002253, s=0.2025,
prevFCN = 27.87419438 f=0.1563, m=2.083, p0=0.00279, s=0.1506,
prevFCN = 28.12016835 f=0.1689, m=1.998, p0=0.001992, s=0.23,
prevFCN = 27.88982196 f=0.165, m=2.028, p0=0.002253, s=0.2025,
prevFCN = 27.87267018 f=0.164,
prevFCN = 27.87573266 f=0.1645, m=2.029,
prevFCN = 27.86358098 m=2.027,
prevFCN = 27.88482321 m=2.028, p0=0.002191,
prevFCN = 27.87325244 p0=0.002317,
prevFCN = 27.87515036 p0=0.002253, s=0.203,
prevFCN = 27.86943774 s=0.202,
prevFCN = 27.8789932 s=0.2028,
prevFCN = 27.87149504 s=0.2022,
prevFCN = 27.87690726 f=0.1678, m=2.025, p0=0.002114, s=0.2145,
prevFCN = 27.78334859 f=0.1813, m=2.015, p0=0.001599, s=0.2624,
prevFCN = 27.57262788 f=0.1863, m=2.011, p0=0.00143, s=0.2798,
prevFCN = 27.54143343 f=0.1904, m=2.008, p0=0.001303, s=0.2937,
prevFCN = 27.52888695 f=0.193, m=2.006, p0=0.001222, s=0.3027,
prevFCN = 27.52564081 f=0.1936,
prevFCN = 27.52428232 f=0.1925,
prevFCN = 27.52701397 f=0.193, m=2.008,
prevFCN = 27.51916389 m=2.005,
prevFCN = 27.53214121 m=2.007,
prevFCN = 27.52075178 m=2.006,
prevFCN = 27.5305432 m=2.006, p0=0.001178,
prevFCN = 27.52494514 p0=0.001268,
prevFCN = 27.52634985 p0=0.001222, s=0.303,
prevFCN = 27.52549414 s=0.3024,
prevFCN = 27.52579186 s=0.3032,
prevFCN = 27.52538731 s=0.3022,
prevFCN = 27.52590768 f=0.1972, m=2.021, p0=0.001181, s=0.2996,
prevFCN = 27.43543975 f=0.2143, m=2.078, p0=0.001023, s=0.2869,
prevFCN = 27.09091162 f=0.2593, m=2.222, p0=0.0006772, s=0.2553,
prevFCN = 26.49313792 f=0.3114, m=2.377, p0=0.0003818, s=0.2209,
prevFCN = 26.77270369 f=0.2701, m=2.254, p0=0.0006078, s=0.248,
prevFCN = 26.44774367 f=0.2707,
prevFCN = 26.44785928 f=0.2695,
prevFCN = 26.44764262 f=0.2701, m=2.255,
prevFCN = 26.44771606 m=2.254,
prevFCN = 26.44779912 m=2.255,
prevFCN = 26.44772189 m=2.254,
prevFCN = 26.4477783 m=2.254, p0=0.000577,
prevFCN = 26.44724245 p0=0.0006393,
prevFCN = 26.4482581 p0=0.0006078, s=0.2485,
prevFCN = 26.44920896 s=0.2475,
prevFCN = 26.44629083 f=0.259, m=2.32, p0=9.987e-05, s=0.1948,
prevFCN = 26.48265893 f=0.2655, m=2.282, p0=0.0003445, s=0.2259,
prevFCN = 26.40622954 f=0.266,
prevFCN = 26.40626101 f=0.2649,
prevFCN = 26.40621084 f=0.2655, m=2.282,
prevFCN = 26.40686301 m=2.281,
prevFCN = 26.40561063 m=2.282, p0=0.0003218,
prevFCN = 26.40586112 p0=0.000368,
prevFCN = 26.40661063 p0=0.0003445, s=0.2264,
prevFCN = 26.4070037 s=0.2254,
prevFCN = 26.40547561 f=0.2619, m=2.269, p0=1.164e-05, s=0.2128,
prevFCN = 26.37813788 f=0.2603, m=2.263, p0=9.758e-06, s=0.2071,
prevFCN = 26.37648547 f=0.2607, m=2.265, p0=2.459e-06, s=0.2084,
prevFCN = 26.37624118 f=0.2613,
prevFCN = 26.3761512 f=0.2602,
prevFCN = 26.37634401 f=0.2607, m=2.265,
prevFCN = 26.37638659 m=2.264,
prevFCN = 26.37611051 m=2.265, p0=4.814e-06,
prevFCN = 26.37627918 p0=8.88e-07,
prevFCN = 26.37621583 p0=2.459e-06, s=0.2089,
prevFCN = 26.37618008 s=0.208,
prevFCN = 26.37632149 f=0.266, m=2.258, p0=2.905e-06, s=0.2111,
prevFCN = 26.37502847 f=0.2652, m=2.259, p0=2.835e-06, s=0.2107,
prevFCN = 26.37498935 f=0.2657,
prevFCN = 26.37499434 f=0.2646,
prevFCN = 26.37499719 f=0.2652, m=2.26,
prevFCN = 26.37500579 m=2.259,
prevFCN = 26.37498568 m=2.259, p0=5.355e-06,
prevFCN = 26.37503005 p0=1.11e-06,
prevFCN = 26.37496148 p0=2.835e-06, s=0.211,
prevFCN = 26.37501113 s=0.2104,
prevFCN = 26.37497995 MIGRAD MINIMIZATION HAS CONVERGED.
FCN=26.375 FROM MIGRAD STATUS=CONVERGED 166 CALLS 167 TOTAL
EDM=5.2123e-05 STRATEGY= 1 ERROR MATRIX UNCERTAINTY 4.2 per cent
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 f 2.65182e-01 1.69358e-01 1.01521e-02 -1.15686e-03
2 m 2.25938e+00 1.48943e-01 -5.27364e-04 1.95576e-01
3 p0 2.83539e-06 2.52439e-02 2.31288e-04 2.72050e-02
4 s 2.10717e-01 1.00428e-01 2.26902e-04 4.71766e-01
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 4 ERR DEF=0.5
3.021e-02 2.135e-03 9.360e-06 -1.319e-03
2.135e-03 2.219e-02 -1.105e-06 -1.134e-03
9.360e-06 -1.105e-06 2.887e-07 -8.247e-06
-1.319e-03 -1.134e-03 -8.247e-06 1.009e-02
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4
1 0.14158 1.000 0.082 0.100 -0.076
2 0.11285 0.082 1.000 -0.014 -0.076
3 0.17975 0.100 -0.014 1.000 -0.153
4 0.17988 -0.076 -0.076 -0.153 1.000
s=0.2107, **********
** 7 **SET ERR 0.5
**********
**********
** 8 **SET PRINT 1
**********
**********
** 9 **HESSE 2000
**********
prevFCN = 26.37498935 f=0.2657,
prevFCN = 26.37499434 f=0.2646,
prevFCN = 26.37499719 f=0.2652, m=2.26,
prevFCN = 26.37500579 m=2.259,
prevFCN = 26.37498568 m=2.259, p0=5.355e-06,
prevFCN = 26.37503005 p0=1.11e-06,
prevFCN = 26.37496148 p0=2.835e-06, s=0.211,
prevFCN = 26.37501113 s=0.2104,
prevFCN = 26.37497995 f=0.2653, s=0.2107,
prevFCN = 26.37498932 f=0.2651,
prevFCN = 26.37498989 f=0.2652, m=2.259,
prevFCN = 26.37499161 m=2.259,
prevFCN = 26.37498759 m=2.259, p0=3.276e-06,
prevFCN = 26.37499646 p0=2.427e-06,
prevFCN = 26.37498275 p0=2.835e-06, s=0.2108,
prevFCN = 26.37499272 s=0.2107,
prevFCN = 26.37498647 f=0.2657, m=2.26, s=0.2107,
prevFCN = 26.37501147 m=2.259, p0=5.355e-06,
prevFCN = 26.37503504 f=0.2652, m=2.26,
prevFCN = 26.37504649 f=0.2657, m=2.259, p0=2.835e-06, s=0.211,
prevFCN = 26.37501574 f=0.2652, m=2.26,
prevFCN = 26.37502857 m=2.259, p0=5.355e-06,
prevFCN = 26.37505184 COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=26.375 FROM HESSE STATUS=OK 23 CALLS 190 TOTAL
EDM=0.000117538 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 f 2.65182e-01 1.49160e-01 1.23217e-03 -4.88878e-01
2 m 2.25938e+00 1.40794e-01 5.14100e-05 2.27906e-01
3 p0 2.83539e-06 3.06275e-02 1.26031e-03 -1.56743e+00
4 s 2.10717e-01 9.42272e-02 3.30485e-05 2.10732e-02
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 4 ERR DEF=0.5
2.315e-02 -1.213e-03 -3.313e-08 5.042e-04
-1.213e-03 1.982e-02 8.012e-09 -1.075e-03
-3.313e-08 8.012e-09 3.510e-07 -3.512e-09
5.042e-04 -1.075e-03 -3.512e-09 8.879e-03
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4
1 0.06441 1.000 -0.057 -0.000 0.035
2 0.09725 -0.057 1.000 0.000 -0.081
3 0.00038 -0.000 0.000 1.000 -0.000
4 0.08659 0.035 -0.081 -0.000 1.000
p0=2.835e-06, s=0.2107, [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization
[#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization
[#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (g,p)
**********
** 10 **SET PRINT 1
**********
**********
** 11 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 f 2.65182e-01 1.49160e-01 0.00000e+00 1.00000e+00
2 m 2.25938e+00 1.40794e-01 -1.00000e+01 1.00000e+01
3 p0 2.83539e-06 3.06275e-02 0.00000e+00 1.00000e+00
MINUIT WARNING IN PARAMETR
============== VARIABLE3 BROUGHT BACK INSIDE LIMITS.
4 s 2.10717e-01 9.42272e-02 -1.00000e+01 1.00000e+01
**********
** 12 **SET ERR 0.5
**********
**********
** 13 **SET PRINT 1
**********
**********
** 14 **SET STR 1
**********
NOW USING STRATEGY 1: TRY TO BALANCE SPEED AGAINST RELIABILITY
**********
** 15 **MIGRAD 2000 1
**********
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
prevFCN = 26.37498935 START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03
f=0.2667,
prevFCN = 26.37503756 f=0.2636,
prevFCN = 26.3750457 f=0.2657,
prevFCN = 26.37499433 f=0.2646,
prevFCN = 26.37499718 f=0.2652, m=2.261,
prevFCN = 26.37506807 m=2.258,
prevFCN = 26.37501158 m=2.26,
prevFCN = 26.37500583 m=2.259,
prevFCN = 26.37498568 m=2.259, p0=6.546e-06,
prevFCN = 26.3750493 p0=6.547e-07,
prevFCN = 26.37495412 p0=5.353e-06,
prevFCN = 26.37503002 p0=1.111e-06,
prevFCN = 26.37496149 p0=2.835e-06, s=0.2117,
prevFCN = 26.37508388 s=0.2098,
prevFCN = 26.37499557 s=0.2111,
prevFCN = 26.37501161 s=0.2104,
prevFCN = 26.3749799 FCN=26.375 FROM MIGRAD STATUS=INITIATE 16 CALLS 17 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 f 2.65182e-01 1.49160e-01 3.51738e-01 -1.15686e-03
2 m 2.25938e+00 1.40794e-01 1.44537e-02 1.95575e-01
3 p0 2.83539e-06 3.06275e-02 1.74945e-01 2.72050e-02
4 s 2.10717e-01 9.42272e-02 9.42495e-03 4.71750e-01
ERR DEF= 0.5
f=0.2652, m=2.259, p0=1.388e-14, s=0.2103,
prevFCN = 26.37493073 f=0.2658,
prevFCN = 26.37493708 f=0.2647,
prevFCN = 26.37493718 f=0.2652, m=2.259,
prevFCN = 26.37493591 m=2.258,
prevFCN = 26.3749384 m=2.259, p0=3.967e-07,
prevFCN = 26.37493713 p0=3.964e-07,
prevFCN = 26.37493713 p0=1.388e-14, s=0.2106,
prevFCN = 26.37493624 s=0.21,
prevFCN = 26.37493816 MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
s=0.2103,
prevFCN = 26.37493073 f=0.2658,
prevFCN = 26.37493708 f=0.2647,
prevFCN = 26.37493718 f=0.2652, m=2.259,
prevFCN = 26.37493591 m=2.258,
prevFCN = 26.3749384 m=2.259, p0=3.967e-07,
prevFCN = 26.37493713 p0=3.964e-07,
prevFCN = 26.37493713 p0=1.388e-14, s=0.2106,
prevFCN = 26.37493624 s=0.21,
prevFCN = 26.37493816 f=0.2654, s=0.2103,
prevFCN = 26.37493097 f=0.2651,
prevFCN = 26.37493099 f=0.2652, m=2.259,
prevFCN = 26.37493073 m=2.259,
prevFCN = 26.37493123 m=2.259, p0=1.589e-08,
prevFCN = 26.37493098 p0=1.583e-08,
prevFCN = 26.37493098 p0=1.388e-14, s=0.2104,
prevFCN = 26.3749308 s=0.2102,
prevFCN = 26.37493117 f=0.2658, m=2.259, s=0.2103,
prevFCN = 26.37494295 m=2.259, p0=3.967e-07,
prevFCN = 26.37494349 f=0.2652, m=2.259,
prevFCN = 26.37494232 f=0.2658, m=2.259, p0=1.388e-14, s=0.2106,
prevFCN = 26.3749422 f=0.2652, m=2.259,
prevFCN = 26.37494248 m=2.259, p0=3.967e-07,
prevFCN = 26.37494265 COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=26.3749 FROM MIGRAD STATUS=CONVERGED 48 CALLS 49 TOTAL
EDM=1.75981e-07 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 f 2.65243e-01 1.49171e-01 1.23134e-03 -4.01360e-05
2 m 2.25899e+00 1.40574e-01 5.14960e-05 -2.42389e-02
3 p0 1.38778e-14 3.06306e-02 1.25947e-03** at limit **
4 s 2.10301e-01 9.37700e-02 3.36127e-05 -2.80291e-02
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 4 ERR DEF=0.5
2.315e-02 -1.209e-03 -3.656e-13 5.077e-04
-1.209e-03 1.976e-02 8.789e-14 -1.103e-03
-3.656e-13 8.789e-14 1.721e-15 -3.915e-14
5.077e-04 -1.103e-03 -3.915e-14 8.793e-03
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4
1 0.06443 1.000 -0.056 -0.000 0.036
2 0.09933 -0.056 1.000 0.000 -0.084
3 0.00006 -0.000 0.000 1.000 -0.000
4 0.08918 0.036 -0.084 -0.000 1.000
p0=1.388e-14, s=0.2103, **********
** 16 **SET ERR 0.5
**********
**********
** 17 **SET PRINT 1
**********
**********
** 18 **HESSE 2000
**********
prevFCN = 26.37493073 f=0.2654,
prevFCN = 26.37493097 f=0.2651,
prevFCN = 26.37493099 f=0.2652, m=2.259,
prevFCN = 26.37493073 m=2.259,
prevFCN = 26.37493123 m=2.259, p0=1.589e-08,
prevFCN = 26.37493098 p0=1.583e-08,
prevFCN = 26.37493098 p0=1.388e-14, s=0.2104,
prevFCN = 26.3749308 s=0.2102,
prevFCN = 26.37493117 f=0.2653, s=0.2103,
prevFCN = 26.37493073 f=0.2652,
prevFCN = 26.37493074 f=0.2652, m=2.259,
prevFCN = 26.37493069 m=2.259,
prevFCN = 26.37493079 m=2.259, p0=6.405e-10,
prevFCN = 26.37493074 p0=6.286e-10,
prevFCN = 26.37493074 p0=1.388e-14, s=0.2103,
prevFCN = 26.3749307 s=0.2103,
prevFCN = 26.37493077 f=0.2654, m=2.259, s=0.2103,
prevFCN = 26.37493101 m=2.259, p0=1.589e-08,
prevFCN = 26.37493123 f=0.2652, m=2.259,
prevFCN = 26.37493099 f=0.2654, m=2.259, p0=1.388e-14, s=0.2104,
prevFCN = 26.37493103 f=0.2652, m=2.259,
prevFCN = 26.37493085 m=2.259, p0=1.589e-08,
prevFCN = 26.37493105 COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=26.3749 FROM HESSE STATUS=OK 23 CALLS 72 TOTAL
EDM=1.75675e-07 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 f 2.65243e-01 1.49170e-01 2.46269e-04 -4.88741e-01
2 m 2.25899e+00 1.40589e-01 1.02992e-05 2.27865e-01
3 p0 1.38778e-14 3.06306e-02 2.51894e-04 -1.57080e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
4 s 2.10301e-01 9.37822e-02 6.72254e-06 2.10316e-02
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 4 ERR DEF=0.5
2.315e-02 -1.205e-03 -7.334e-14 5.122e-04
-1.205e-03 1.977e-02 1.745e-14 -1.122e-03
-7.334e-14 1.745e-14 1.721e-15 -7.911e-15
5.122e-04 -1.122e-03 -7.911e-15 8.795e-03
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4
1 0.06440 1.000 -0.056 -0.000 0.036
2 0.10040 -0.056 1.000 0.000 -0.085
3 0.00001 -0.000 0.000 1.000 -0.000
4 0.09060 0.036 -0.085 -0.000 1.000
p0=1.388e-14, s=0.2103, [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization
All Message streams
[0] MinLevel = PROGRESS Topic = Generation Minimization Plotting Fitting Integration LinkStateMgmt Eval Caching Optimization ObjectHandling InputArguments Tracing Contents DataHandling NumericIntegration FastEvaluations
[1] MinLevel = INFO Topic = Minimization Plotting Fitting Eval Caching ObjectHandling InputArguments DataHandling NumericIntegration
[2] MinLevel = INFO Topic = HistFactory
[3] MinLevel = DEBUG Topic = LinkStateMgmt
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server x(0x7ffe2a4b4f60) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server m(0x7ffe2a4b6b60) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server s(0x7ffe2a4b6778) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server x(0x7ffe2a4b4f60) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server m(0x7ffe2a4b6b60) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server s(0x7ffe2a4b6778) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server x(0x558d7c50c0d0) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server m(0x558d7c50c5c0) for value
[#3] DEBUG:LinkStateMgmt -- RooAbsArg::addServer(0x558d7d5f2c60,g): adding server s(0x558d7d5b81b0) for value
RooGaussian::g[ x=x mean=m sigma=s ] = 8.804e-26
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "RooPolynomial.h"
#include "RooAddPdf.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "RooPlot.h"
#include "RooMsgService.h"
using namespace RooFit;
void rf506_msgservice()
{
// C r e a t e p d f
// --------------------
// Construct gauss(x,m,s)
RooRealVar x("x", "x", -10, 10);
RooRealVar m("m", "m", 0, -10, 10);
RooRealVar s("s", "s", 1, -10, 10);
RooGaussian gauss("g", "g", x, m, s);
// Construct poly(x,p0)
RooRealVar p0("p0", "p0", 0.01, 0., 1.);
RooPolynomial poly("p", "p", x, p0);
// Construct model = f*gauss(x) + (1-f)*poly(x)
RooRealVar f("f", "f", 0.5, 0., 1.);
RooAddPdf model("model", "model", RooArgSet(gauss, poly), f);
RooDataSet *data = model.generate(x, 10);
// P r i n t c o n f i g u r a t i o n o f m e s s a g e s e r v i c e
// ---------------------------------------------------------------------------
// Print streams configuration
RooMsgService::instance().Print();
cout << endl;
// A d d i n g I n t e g r a t i o n t o p i c t o e x i s t i n g I N F O s t r e a m
// -----------------------------------------------------------------------------------------------
// Print streams configuration
RooMsgService::instance().Print();
cout << endl;
// Add Integration topic to existing INFO stream
RooMsgService::instance().getStream(1).addTopic(Integration);
// Construct integral over gauss to demonstrate new message stream
RooAbsReal *igauss = gauss.createIntegral(x);
igauss->Print();
// Print streams configuration in verbose, which also shows inactive streams
cout << endl;
RooMsgService::instance().Print();
cout << endl;
// Remove stream
RooMsgService::instance().getStream(1).removeTopic(Integration);
// E x a m p l e s o f p d f v a l u e t r a c i n g s t r e a m
// -----------------------------------------------------------------------
// Show DEBUG level message on function tracing, trace RooGaussian only
RooMsgService::instance().addStream(DEBUG, Topic(Tracing), ClassName("RooGaussian"));
// Perform a fit to generate some tracing messages
model.fitTo(*data, Verbose(kTRUE));
// Reset message service to default stream configuration
RooMsgService::instance().reset();
// Show DEBUG level message on function tracing on all objects, redirect output to file
RooMsgService::instance().addStream(DEBUG, Topic(Tracing), OutputFile("rf506_debug.log"));
// Perform a fit to generate some tracing messages
model.fitTo(*data, Verbose(kTRUE));
// Reset message service to default stream configuration
RooMsgService::instance().reset();
// E x a m p l e o f a n o t h e r d e b u g g i n g s t r e a m
// ---------------------------------------------------------------------
// Show DEBUG level messages on client/server link state management
RooMsgService::instance().addStream(DEBUG, Topic(LinkStateMgmt));
RooMsgService::instance().Print("v");
// Clone composite pdf g to trigger some link state management activity
RooAbsArg *gprime = gauss.cloneTree();
gprime->Print();
// Reset message service to default stream configuration
RooMsgService::instance().reset();
}
DEBUG
#define DEBUG
Definition Polynomial.cxx:40
f
#define f(i)
Definition RSha256.hxx:104
kTRUE
const Bool_t kTRUE
Definition RtypesCore.h:100
RooAbsArg
RooAbsArg is the common abstract base class for objects that represent a value and a "shape" in RooFi...
Definition RooAbsArg.h:69
RooAbsArg::cloneTree
virtual RooAbsArg * cloneTree(const char *newname=0) const
Clone tree expression of objects.
Definition RooAbsArg.cxx:2223
RooAbsArg::Print
virtual void Print(Option_t *options=0) const
Print the object to the defaultPrintStream().
Definition RooAbsArg.h:337
RooAbsReal
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
Definition RooAbsReal.h:64
RooAbsReal::createIntegral
RooAbsReal * createIntegral(const RooArgSet &iset, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none(), const RooCmdArg &arg6=RooCmdArg::none(), const RooCmdArg &arg7=RooCmdArg::none(), const RooCmdArg &arg8=RooCmdArg::none()) const
Create an object that represents the integral of the function over one or more observables listed in ...
Definition RooAbsReal.cxx:554
RooAddPdf
RooAddPdf is an efficient implementation of a sum of PDFs of the form.
Definition RooAddPdf.h:32
RooArgSet
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:35
RooDataSet
RooDataSet is a container class to hold unbinned data.
Definition RooDataSet.h:36
RooGaussian
Plain Gaussian p.d.f.
Definition RooGaussian.h:24
RooMsgService::instance
static RooMsgService & instance()
Return reference to singleton instance.
Definition RooMsgService.cxx:366
RooMsgService::getStream
StreamConfig & getStream(Int_t id)
Definition RooMsgService.h:152
RooMsgService::addStream
Int_t addStream(RooFit::MsgLevel level, const RooCmdArg &arg1=RooCmdArg(), const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg())
Add a message logging stream for message with given RooFit::MsgLevel or higher.
Definition RooMsgService.cxx:198
RooMsgService::Print
void Print(Option_t *options=0) const
Print configuration of message service.
Definition RooMsgService.cxx:554
RooPolynomial
RooPolynomial implements a polynomial p.d.f of the form.
Definition RooPolynomial.h:28
RooRealVar
RooRealVar represents a variable that can be changed from the outside.
Definition RooRealVar.h:39
x
Double_t x[n]
Definition legend1.C:17
RooFit
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition Common.h:18
rf506_msgservice
Definition rf506_msgservice.py:1
RooMsgService::StreamConfig::addTopic
void addTopic(RooFit::MsgTopic newTopic)
Definition RooMsgService.h:114
RooMsgService::StreamConfig::removeTopic
void removeTopic(RooFit::MsgTopic oldTopic)
Definition RooMsgService.h:118
m
auto * m
Definition textangle.C:8
Date
July 2008
Author
Wouter Verkerke

Definition in file rf506_msgservice.C.