34 std::string observablename =
"pTV";
35 RooRealVar obsvar(observablename.c_str(),
"observable of pTV", 10, 600);
37 RooRealVar cHq3(
"cHq3",
"EFT modifier", -10.0, 10.0);
38 cHq3.setAttribute(
"NewPhysics",
true);
39 RooRealVar cHl3(
"cHl3",
"EFT modifier", -10.0, 10.0);
40 cHl3.setAttribute(
"NewPhysics",
true);
41 RooRealVar cHDD(
"cHDD",
"EFT modifier", -10.0, 10.0);
42 cHDD.setAttribute(
"NewPhysics",
true);
46 std::string infilename = std::string(
gROOT->GetTutorialDir()) +
"/roofit/input_histos_rf_lagrangianmorph.root";
47 std::vector<std::string> samplelist = {
"SM_NPsq0",
"cHq3_NPsq1",
"cHq3_NPsq2",
"cHl3_NPsq1",
48 "cHl3_NPsq2",
"cHDD_NPsq1",
"cHDD_NPsq2",
"cHl3_cHDD_NPsq2",
49 "cHq3_cHDD_NPsq2",
"cHl3_cHq3_NPsq2"};
69 morphfunc.setParameter(
"cHq3", 0.01);
70 morphfunc.setParameter(
"cHl3", 1.0);
71 morphfunc.setParameter(
"cHDD", 0.2);
73 auto pseudo_hist = morphfunc.createTH1(
"pseudo_hist");
77 morphfunc.setParameter(
"cHq3", 0.0);
78 morphfunc.setParameter(
"cHl3", 0.0);
79 morphfunc.setParameter(
"cHDD", 0.0);
92 auto hcorr = fitres->correlationHist();
98 auto postfit_hist = morphfunc.createTH1(
"morphing_postfit_hist");
99 RooDataHist postfit_dh(
"morphing_postfit_dh",
"morphing_postfit_dh",
RooArgList(obsvar), postfit_hist);
101 auto frame0 = obsvar.frame(
Title(
"Fitted histogram of p_{T}^{V}"));
104 pseudo_dh->plotOn(frame0,
Name(
"input"));
112 gPad->SetLeftMargin(0.15);
113 gPad->SetRightMargin(0.05);
116 frame0->GetXaxis()->SetTitle(
"p_{T}^{V}");
120 gPad->SetLeftMargin(0.15);
121 gPad->SetRightMargin(0.15);
124 hcorr->SetMarkerSize(3.);
125 hcorr->SetTitle(
"correlation matrix");
126 hcorr->GetYaxis()->SetTitleOffset(1.4);
127 hcorr->GetYaxis()->SetLabelSize(0.1);
128 hcorr->GetXaxis()->SetLabelSize(0.1);
129 hcorr->GetYaxis()->SetBinLabel(1,
"c_{HDD}");
130 hcorr->GetYaxis()->SetBinLabel(2,
"c_{Hl^{(3)}}");
131 hcorr->GetYaxis()->SetBinLabel(3,
"c_{Hq^{(3)}}");
132 hcorr->GetXaxis()->SetBinLabel(3,
"c_{HDD}");
133 hcorr->GetXaxis()->SetBinLabel(2,
"c_{Hl^{(3)}}");
134 hcorr->GetXaxis()->SetBinLabel(1,
"c_{Hq^{(3)}}");
135 hcorr->Draw(
"colz text");
136 c1->SaveAs(
"rf712_lagrangianmorphfit.png");
R__EXTERN TStyle * gStyle
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.
Container class to hold N-dimensional binned data.
Class RooLagrangianMorphing is a implementation of the method of Effective Lagrangian Morphing,...
Variable that can be changed from the outside.
The RooWrapperPdf is a class that can be used to convert a function into a PDF.
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
void SetPaintTextFormat(const char *format="g")
RooCmdArg Parameters(const RooArgSet ¶ms)
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 PrintLevel(Int_t code)
RooCmdArg DrawOption(const char *opt)
RooCmdArg XErrorSize(double width)
RooCmdArg LineColor(Color_t color)
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
std::string observableName
std::vector< std::string > folderNames