void rf711_lagrangianmorph()
{
cHq3.setAttribute(
"NewPhysics",
true);
std::string
infilename = std::string(
gROOT->GetTutorialDir()) +
"/roofit/roofit/input_histos_rf_lagrangianmorph.root";
std::vector<std::string>
samplelist = {
"SM_NPsq0",
"cHq3_NPsq1",
"cHq3_NPsq2"};
file->Close();
std::unique_ptr<RooDataSet>
data{model.generate({
cHq3,
obsvar}, 1000000)};
hh_data->SetTitle(
"Morphing prediction");
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.05);
leg1->AddEntry(
"SM_NPsq0",
"SM",
"LP");
leg1->AddEntry(
"cHq3_NPsq1",
"c_{Hq^(3)}=1.0 at #Lambda^{-2}",
"LP");
leg1->AddEntry(
"cHq3_NPsq2",
"c_{Hq^(3)}=1.0 at #Lambda^{-4}",
"LP");
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.05);
leg2->AddEntry(
"morph_dh_cHq3=0.01",
"c_{Hq^{(3)}}=0.01",
"L");
leg2->AddEntry(
"morph_dh_cHq3=0.25",
"c_{Hq^{(3)}}=0.25",
"L");
leg2->AddEntry(
"morph_dh_cHq3=0.5",
"c_{Hq^{(3)}}=0.5",
"L");
gPad->SetLeftMargin(0.12);
gPad->SetRightMargin(0.18);
hh_data->GetYaxis()->SetTitle(
"c_{Hq^{(3)}}");
hh_data->GetYaxis()->SetRangeUser(0, 0.5);
hh_data->GetZaxis()->SetTitleOffset(1.8);
c1->SaveAs(
"rf711_lagrangianmorph.png");
}
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void data
R__EXTERN TStyle * gStyle
virtual bool add(const RooAbsArg &var, bool silent=false)
Add the specified argument to list.
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.
static void InvertPalette()
Invert the current color palette.
A ROOT file is an on-disk file, usually with extension .root, that stores objects in a file-system-li...
<div class="legacybox"><h2>Legacy Code</h2> TFolder is a legacy interface: there will be no bug fixes...
TH1 is the base class of all histogram classes in ROOT.
This class displays a legend box (TPaveText) containing several legend entries.
Mother of all ROOT objects.
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
void SetPalette(Int_t ncolors=kBird, Int_t *colors=nullptr, Float_t alpha=1.)
See TColor::SetPalette.
void SetNumberContours(Int_t number=20)
Set the default number of contour levels when drawing 2-d plots.
RooCmdArg YVar(const RooAbsRealLValue &var, const RooCmdArg &arg={})
RooCmdArg DataError(Int_t)
RooCmdArg Binning(const RooAbsBinning &binning)
RooCmdArg DrawOption(const char *opt)
RooCmdArg LineColor(TColorNumber color)
RooCmdArg XErrorSize(double width)
RooCmdArg MarkerSize(Size_t size)
RooCmdArg MarkerColor(TColorNumber 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
[#0] PROGRESS:InputArguments -- initializing physics inputs from file /github/home/ROOT-CI/build/tutorials/roofit/roofit/input_histos_rf_lagrangianmorph.root with object name(s) 'pTV'
[#0] PROGRESS:Caching -- creating cache from getCache function for 0x7fff484195e0
[#0] PROGRESS:Caching -- current storage has size 3
[#0] PROGRESS:ObjectHandling -- observable: pTV
[#0] PROGRESS:ObjectHandling -- binWidth: binWidth_pTV
[#1] INFO:DataHandling -- RooDataHist::adjustBinning(morph_dh_cHq3=0.01): fit range of variable pTV expanded to nearest bin boundaries: [10,600] --> [0,600]
[#1] INFO:InputArguments -- RooAbsData::plotOn(SM_NPsq0) INFO: dataset has non-integer weights, auto-selecting SumW2 errors instead of Poisson errors
[#1] INFO:InputArguments -- RooAbsData::plotOn(cHq3_NPsq1) INFO: dataset has non-integer weights, auto-selecting SumW2 errors instead of Poisson errors
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 24931.9 will supersede previous event count of 10852.3 for normalization of PDF projections
[#1] INFO:InputArguments -- RooAbsData::plotOn(cHq3_NPsq2) INFO: dataset has non-integer weights, auto-selecting SumW2 errors instead of Poisson errors
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 29789.2 will supersede previous event count of 24931.9 for normalization of PDF projections
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 18947.1 will supersede previous event count of 11104.6 for normalization of PDF projections
[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 30765.5 will supersede previous event count of 18947.1 for normalization of PDF projections
[#0] PROGRESS:Caching -- creating cache from getCache function for 0x5c2aab0
[#0] PROGRESS:Caching -- current storage has size 3
[#0] PROGRESS:ObjectHandling -- observable: pTV
[#0] PROGRESS:ObjectHandling -- binWidth: binWidth_pTV
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_SM_NPsq0_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_cHq3_NPsq1_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_cHq3_NPsq2_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)
[#0] PROGRESS:Caching -- creating cache from getCache function for 0x72f3580
[#0] PROGRESS:Caching -- current storage has size 3
[#0] PROGRESS:ObjectHandling -- observable: pTV
[#0] PROGRESS:ObjectHandling -- binWidth: binWidth_pTV
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_SM_NPsq0_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_cHq3_NPsq1_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(w_cHq3_NPsq2_morphfunc_Int[cHq3]) using numeric integrator RooIntegrator1D to calculate Int(cHq3)