Logo ROOT  
Reference Guide
rf712_lagrangianmorphfit.C File Reference

Detailed Description

View in nbviewer Open in SWAN
Performing a simple fit with RooLagrangianMorphFunc.

a morphing function is setup as a function of three variables and a fit is performed on a pseudo-dataset.

[#0] PROGRESS:InputArguments -- initializing physics inputs from file /home/sftnight/build/workspace/root-makedoc-master/rootspi/rdoc/src/master.build/tutorials/roofit/input_histos_rf_lagrangianmorph.root with object name(s) 'pTV'
[#0] PROGRESS:Caching -- creating cache from getCache function for 0x7ffe254e0430
[#0] PROGRESS:Caching -- current storage has size 10
[#0] PROGRESS:ObjectHandling -- observable: pTV
[#0] PROGRESS:ObjectHandling -- binWidth: binWidth_pTV
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(pseudo_dh): fit range of variable pTV expanded to nearest bin boundaries: [10,600] --> [0,600]
[#0] PROGRESS:Caching -- creating cache from getCache function for 0x55dda54fc980
[#0] PROGRESS:Caching -- current storage has size 10
[#0] PROGRESS:ObjectHandling -- observable: pTV
[#0] PROGRESS:ObjectHandling -- binWidth: binWidth_pTV
**********
** 1 **SET PRINT 1
**********
**********
** 2 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 cHDD 0.00000e+00 1.00000e-01 -1.00000e+01 1.00000e+01
2 cHl3 0.00000e+00 1.00000e-01 -1.00000e+01 1.00000e+01
3 cHq3 0.00000e+00 1.00000e-01 -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 1500 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=39609.1 FROM MIGRAD STATUS=INITIATE 14 CALLS 15 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 cHDD 0.00000e+00 1.00000e-01 1.00002e-02 -2.72888e-01
2 cHl3 0.00000e+00 1.00000e-01 1.00002e-02 -1.83697e+01
3 cHq3 0.00000e+00 1.00000e-01 1.00002e-02 -3.57539e+03
ERR DEF= 0.5
RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (39615.1) to force MIGRAD to back out of this region. Error log follows.
Parameter values: cHDD=0.223088 cHl3=2.02101 cHq3=0.00780325
RooNLLVar::nll_wrap_pdf_pseudo_dh[ parameters=(binWidth_pTV,cHDD,cHl3,cHq3,kSM,nNP0,nNP1,nNP2,nNP3,nNP4) ]
function value is NAN @ parameters=(binWidth_pTV = 0.05,cHDD = 0.223088 +/- 0.1,cHl3 = 2.02101 +/- 0.1,cHq3 = 0.00780325 +/- 0.1,kSM = 1,nNP0 = 1,nNP1 = 1,nNP2 = 1,nNP3 = 1,nNP4 = 1)
RooWrapperPdf::wrap_pdf[ inputFunction=morphfunc ]
p.d.f value is less than zero (-0.041358), trying to recover @ inputFunction=morphfunc=-0.0413585
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0413585
p.d.f value is less than zero (-0.043016), trying to recover @ inputFunction=morphfunc=-0.043016
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.043016
RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (39615.1) to force MIGRAD to back out of this region. Error log follows.
Parameter values: cHDD=0.56869 cHl3=4.96014 cHq3=0.0199009
RooNLLVar::nll_wrap_pdf_pseudo_dh[ parameters=(binWidth_pTV,cHDD,cHl3,cHq3,kSM,nNP0,nNP1,nNP2,nNP3,nNP4) ]
function value is NAN @ parameters=(binWidth_pTV = 0.05,cHDD = 0.56869 +/- 0.1,cHl3 = 4.96014 +/- 0.1,cHq3 = 0.0199009 +/- 0.1,kSM = 1,nNP0 = 1,nNP1 = 1,nNP2 = 1,nNP3 = 1,nNP4 = 1)
RooWrapperPdf::wrap_pdf[ inputFunction=morphfunc ]
p.d.f value is less than zero (-0.074541), trying to recover @ inputFunction=morphfunc=-0.0745407
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0745407
p.d.f value is less than zero (-1.173323), trying to recover @ inputFunction=morphfunc=-1.17332
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-1.17332
p.d.f value is less than zero (-1.063175), trying to recover @ inputFunction=morphfunc=-1.06317
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-1.06317
p.d.f value is less than zero (-1.234351), trying to recover @ inputFunction=morphfunc=-1.23435
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-1.23435
p.d.f value is less than zero (-1.674381), trying to recover @ inputFunction=morphfunc=-1.67438
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-1.67438
p.d.f value is less than zero (-0.075079), trying to recover @ inputFunction=morphfunc=-0.0750793
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0750793
... (remaining 10 messages suppressed)
RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (39615.1) to force MIGRAD to back out of this region. Error log follows.
Parameter values: cHDD=0.318765 cHl3=2.86725 cHq3=0.0111509
RooNLLVar::nll_wrap_pdf_pseudo_dh[ parameters=(binWidth_pTV,cHDD,cHl3,cHq3,kSM,nNP0,nNP1,nNP2,nNP3,nNP4) ]
function value is NAN @ parameters=(binWidth_pTV = 0.05,cHDD = 0.318765 +/- 0.1,cHl3 = 2.86725 +/- 0.1,cHq3 = 0.0111509 +/- 0.1,kSM = 1,nNP0 = 1,nNP1 = 1,nNP2 = 1,nNP3 = 1,nNP4 = 1)
RooWrapperPdf::wrap_pdf[ inputFunction=morphfunc ]
p.d.f value is less than zero (-0.070498), trying to recover @ inputFunction=morphfunc=-0.0704979
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0704979
p.d.f value is less than zero (-0.063999), trying to recover @ inputFunction=morphfunc=-0.0639991
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0639991
p.d.f value is less than zero (-0.140607), trying to recover @ inputFunction=morphfunc=-0.140607
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.140607
RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (39615.1) to force MIGRAD to back out of this region. Error log follows.
Parameter values: cHDD=0.193718 cHl3=1.7579 cHq3=0.00677582
RooNLLVar::nll_wrap_pdf_pseudo_dh[ parameters=(binWidth_pTV,cHDD,cHl3,cHq3,kSM,nNP0,nNP1,nNP2,nNP3,nNP4) ]
function value is NAN @ parameters=(binWidth_pTV = 0.05,cHDD = 0.193718 +/- 0.1,cHl3 = 1.7579 +/- 0.1,cHq3 = 0.00677582 +/- 0.1,kSM = 1,nNP0 = 1,nNP1 = 1,nNP2 = 1,nNP3 = 1,nNP4 = 1)
RooWrapperPdf::wrap_pdf[ inputFunction=morphfunc ]
p.d.f value is less than zero (-0.029110), trying to recover @ inputFunction=morphfunc=-0.0291104
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0291104
p.d.f value is less than zero (-0.011532), trying to recover @ inputFunction=morphfunc=-0.0115319
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0115319
RooAbsMinimizerFcn: Minimized function has error status.
Returning maximum FCN so far (39615.1) to force MIGRAD to back out of this region. Error log follows.
Parameter values: cHDD=0.820784 cHl3=1.67112 cHq3=0.0140305
RooNLLVar::nll_wrap_pdf_pseudo_dh[ parameters=(binWidth_pTV,cHDD,cHl3,cHq3,kSM,nNP0,nNP1,nNP2,nNP3,nNP4) ]
function value is NAN @ parameters=(binWidth_pTV = 0.05,cHDD = 0.820784 +/- 0.1,cHl3 = 1.67112 +/- 0.1,cHq3 = 0.0140305 +/- 0.1,kSM = 1,nNP0 = 1,nNP1 = 1,nNP2 = 1,nNP3 = 1,nNP4 = 1)
RooWrapperPdf::wrap_pdf[ inputFunction=morphfunc ]
p.d.f value is less than zero (-0.010720), trying to recover @ inputFunction=morphfunc=-0.0107196
getLogVal() top-level p.d.f evaluates to NaN @ inputFunction=morphfunc=-0.0107196
MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=39604.9 FROM MIGRAD STATUS=CONVERGED 97 CALLS 98 TOTAL
EDM=4.90989e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 cHDD 1.98078e-01 1.34076e+00 1.76827e-02 -1.24714e-02
2 cHl3 9.98258e-01 2.82061e-01 3.84123e-03 -2.35191e-01
3 cHq3 1.00089e-02 4.07058e-03 5.39370e-05 6.27101e+00
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 3 ERR DEF=0.5
1.809e+00 -5.353e-02 -1.409e-03
-5.353e-02 7.958e-02 1.262e-04
-1.409e-03 1.262e-04 1.657e-05
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3
1 0.28129 1.000 -0.141 -0.257
2 0.16034 -0.141 1.000 0.110
3 0.26790 -0.257 0.110 1.000
**********
** 7 **SET ERR 0.5
**********
**********
** 8 **SET PRINT 1
**********
**********
** 9 **HESSE 1500
**********
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=39604.9 FROM HESSE STATUS=OK 16 CALLS 114 TOTAL
EDM=4.93538e-05 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 cHDD 1.98078e-01 1.32824e+00 7.07309e-04 1.98091e-02
2 cHl3 9.98258e-01 2.81948e-01 1.53649e-04 9.99923e-02
3 cHq3 1.00089e-02 4.04044e-03 2.15748e-06 1.00089e-03
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 3 ERR DEF=0.5
1.775e+00 -4.527e-02 -1.221e-03
-4.527e-02 7.952e-02 1.208e-04
-1.221e-03 1.208e-04 1.633e-05
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3
1 0.24668 1.000 -0.121 -0.227
2 0.14509 -0.121 1.000 0.106
3 0.24026 -0.227 0.106 1.000
[#1] INFO:Fitting -- RooAbsPdf::fitTo(wrap_pdf) Calculating sum-of-weights-squared correction matrix for covariance matrix
**********
** 10 **SET ERR 0.5
**********
**********
** 11 **SET PRINT 1
**********
**********
** 12 **HESSE 1500
**********
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=40071 FROM HESSE STATUS=OK 22 CALLS 136 TOTAL
EDM=4.93731 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 cHDD 1.98078e-01 1.36717e+00 1.44291e-02 1.98091e-02
2 cHl3 9.98258e-01 1.81304e-01 3.13445e-03 9.99923e-02
3 cHq3 1.00089e-02 4.47315e-03 4.40126e-05 1.00089e-03
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 3 ERR DEF=0.5
1.881e+00 6.681e-02 -2.510e-03
6.681e-02 3.287e-02 -3.029e-06
-2.510e-03 -3.029e-06 2.001e-05
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3
1 0.48867 1.000 0.269 -0.409
2 0.29280 0.269 1.000 -0.004
3 0.42378 -0.409 -0.004 1.000
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(morphing_postfit_dh): fit range of variable pTV expanded to nearest bin boundaries: [0,600] --> [0,600]
[#1] INFO:InputArguments -- RooAbsData::plotOn(pseudo_dh) INFO: dataset has non-integer weights, auto-selecting SumW2 errors instead of Poisson errors
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 7389.24 will supercede previous event count of 7396.07 for normalization of PDF projections
#include <RooDataHist.h>
#include <RooFitResult.h>
#include <RooPlot.h>
#include <RooRealVar.h>
#include <TAxis.h>
#include <TCanvas.h>
#include <TH2.h>
#include <TStyle.h>
using namespace RooFit;
{
// C r e a t e v a r i a b l e s f o r
// m o r p h i n g f u n c t i o n
// ---------------------------------------------
std::string observablename = "pTV";
RooRealVar obsvar(observablename.c_str(), "observable of pTV", 10, 600);
RooRealVar kSM("kSM", "sm modifier", 1.0);
RooRealVar cHq3("cHq3", "EFT modifier", -10.0, 10.0);
cHq3.setAttribute("NewPhysics", true);
RooRealVar cHl3("cHl3", "EFT modifier", -10.0, 10.0);
cHl3.setAttribute("NewPhysics", true);
RooRealVar cHDD("cHDD", "EFT modifier", -10.0, 10.0);
cHDD.setAttribute("NewPhysics", true);
// I n p u t s n e e d e d f o r c o n f i g
// ---------------------------------------------
std::string infilename = std::string(gROOT->GetTutorialDir()) + "/roofit/input_histos_rf_lagrangianmorph.root";
std::vector<std::string> samplelist = {"SM_NPsq0", "cHq3_NPsq1", "cHq3_NPsq2", "cHl3_NPsq1",
"cHl3_NPsq2", "cHDD_NPsq1", "cHDD_NPsq2", "cHl3_cHDD_NPsq2",
"cHq3_cHDD_NPsq2", "cHl3_cHq3_NPsq2"};
// S e t u p C o n f i g
// ---------------------------------------------
config.fileName = infilename;
config.observableName = observablename;
config.folderNames = samplelist;
config.couplings.add(cHq3);
config.couplings.add(cHl3);
config.couplings.add(cHDD);
config.couplings.add(kSM);
// C r e a t e m o r p h i n g f u n c t i o n
// ---------------------------------------------
RooLagrangianMorphFunc morphfunc("morphfunc", "morphed dist. of pTV", config);
// C r e a t e p s e u d o d a t a h i s t o g r a m
// f o r f i t
// ---------------------------------------------
morphfunc.setParameter("cHq3", 0.01);
morphfunc.setParameter("cHl3", 1.0);
morphfunc.setParameter("cHDD", 0.2);
auto pseudo_hist = morphfunc.createTH1("pseudo_hist");
auto pseudo_dh = new RooDataHist("pseudo_dh", "pseudo_dh", RooArgList(obsvar), pseudo_hist);
// reset parameters to zeros before fit
morphfunc.setParameter("cHq3", 0.0);
morphfunc.setParameter("cHl3", 0.0);
morphfunc.setParameter("cHDD", 0.0);
// error set used as initial step size
cHq3.setError(0.1);
cHl3.setError(0.1);
cHDD.setError(0.1);
// W r a p p d f o n m o r p h f u n c a n d
// f i t t o d a t a h i s t o g r a m
// ---------------------------------------------
// wrapper pdf to normalise morphing function to a morphing pdf
RooWrapperPdf model("wrap_pdf", "wrap_pdf", morphfunc);
auto fitres = model.fitTo(*pseudo_dh, SumW2Error(true), Optimize(false), Save());
auto hcorr = fitres->correlationHist();
// E x t r a c t p o s t f i t d i s t r i b u t i o n
// a n d p l o t w i t h i n i t i a l
// h i s t o g r a m
// ---------------------------------------------
auto postfit_hist = morphfunc.createTH1("morphing_postfit_hist");
RooDataHist postfit_dh("morphing_postfit_dh", "morphing_postfit_dh", RooArgList(obsvar), postfit_hist);
auto frame0 = obsvar.frame(Title("Fitted histogram of p_{T}^{V}"));
postfit_dh.plotOn(frame0, Name("postfit_dist"), DrawOption("C"), LineColor(kBlue), DataError(RooAbsData::None),
pseudo_dh->plotOn(frame0, Name("input"));
// D r a w p l o t s o n c a n v a s
// ---------------------------------------------
TCanvas *c1 = new TCanvas("fig3", "fig3", 800, 400);
c1->Divide(2, 1);
c1->cd(1);
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.05);
model.paramOn(frame0, Layout(0.50, 0.75, 0.9), Parameters(config.couplings));
frame0->GetXaxis()->SetTitle("p_{T}^{V}");
frame0->Draw();
c1->cd(2);
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.15);
hcorr->SetMarkerSize(3.);
hcorr->SetTitle("correlation matrix");
hcorr->GetYaxis()->SetTitleOffset(1.4);
hcorr->GetYaxis()->SetLabelSize(0.1);
hcorr->GetXaxis()->SetLabelSize(0.1);
hcorr->GetYaxis()->SetBinLabel(1, "c_{HDD}");
hcorr->GetYaxis()->SetBinLabel(2, "c_{Hl^{(3)}}");
hcorr->GetYaxis()->SetBinLabel(3, "c_{Hq^{(3)}}");
hcorr->GetXaxis()->SetBinLabel(3, "c_{HDD}");
hcorr->GetXaxis()->SetBinLabel(2, "c_{Hl^{(3)}}");
hcorr->GetXaxis()->SetBinLabel(1, "c_{Hq^{(3)}}");
hcorr->Draw("colz text");
c1->SaveAs("rf712_lagrangianmorphfit.png");
}
@ kBlue
Definition: Rtypes.h:66
#define gROOT
Definition: TROOT.h:404
R__EXTERN TStyle * gStyle
Definition: TStyle.h:414
#define gPad
Definition: TVirtualPad.h:288
virtual bool add(const RooAbsArg &var, bool silent=false)
Add the specified argument to list.
RooArgList is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgList.h:22
The RooDataHist is a container class to hold N-dimensional binned data.
Definition: RooDataHist.h:39
Class RooLagrangianMorphing is a implementation of the method of Effective Lagrangian Morphing,...
RooRealVar represents a variable that can be changed from the outside.
Definition: RooRealVar.h:40
The RooWrapperPdf is a class that can be used to convert a function into a PDF.
Definition: RooWrapperPdf.h:24
The Canvas class.
Definition: TCanvas.h:23
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:1589
void SetPaintTextFormat(const char *format="g")
Definition: TStyle.h:370
RooCmdArg Parameters(const RooArgSet &params)
RooCmdArg Layout(double xmin, double xmax=0.99, double ymin=0.95)
RooCmdArg Save(bool flag=true)
RooCmdArg SumW2Error(bool flag)
RooCmdArg Optimize(Int_t flag=2)
RooCmdArg DataError(Int_t)
RooCmdArg DrawOption(const char *opt)
RooCmdArg XErrorSize(double width)
RooCmdArg LineColor(Color_t color)
return c1
Definition: legend1.C:41
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition: Common.h:18
const char * Name
Definition: TXMLSetup.cxx:67
const char * Title
Definition: TXMLSetup.cxx:68
std::vector< std::string > folderNames
Date
January 2022
Author
Rahul Balasubramanian

Definition in file rf712_lagrangianmorphfit.C.