#include <FitWithRooFit.h>

Public Member Functions
RooRealVar *	a0 ()

RooRealVar *	a1 ()

RooRealVar *	a2 ()

RooRealVar *	a3 ()

RooRealVar *	a4 ()

RooRealVar *	a5 ()

RooRealVar *	a6 ()

RooRealVar *	alpha ()

RooAbsPdf *	buildBackgroundModel (RooRealVar *x, const TString &backgroundType)
	Build the model for the specified background type. More...

RooAbsPdf *	buildModel (RooRealVar *x, const TString &signalType, const TString &backgroundType)
	Build the model to fit. More...

RooAbsPdf *	buildSignalModel (RooRealVar *x, const TString &signalType)
	Build the model for the specified signal type. More...

RooRealVar *	expCoeffa0 ()

RooRealVar *	expCoeffa1 ()

RooRealVar *	expCoeffa2 ()

RooRealVar *	fGCB ()

void	fit (TH1 *histo, const TString signalType, const TString backgroundType, const double &xMin=0., const double &xMax=0., bool sumW2Error=false)

	FitWithRooFit ()

RooRealVar *	fsig ()

RooRealVar *	gamma ()

RooRealVar *	gaussFrac ()

RooRealVar *	gaussFrac2 ()

rooPair	importTH1 (TH1 *histo, const double &inputXmin, const double &inputXmax)

void	initA0 (const double &value, const double &min, const double &max, const TString &name="a0", const TString &title="a0")

void	initA1 (const double &value, const double &min, const double &max, const TString &name="a1", const TString &title="a1")

void	initA2 (const double &value, const double &min, const double &max, const TString &name="a2", const TString &title="a2")

void	initA3 (const double &value, const double &min, const double &max, const TString &name="a3", const TString &title="a3")

void	initA4 (const double &value, const double &min, const double &max, const TString &name="a4", const TString &title="a4")

void	initA5 (const double &value, const double &min, const double &max, const TString &name="a5", const TString &title="a5")

void	initA6 (const double &value, const double &min, const double &max, const TString &name="a6", const TString &title="a6")

void	initAlpha (const double &value, const double &min, const double &max, const TString &name="alpha", const TString &title="alpha")

void	initExpCoeffA0 (const double &value, const double &min, const double &max, const TString &name="expCoeffa0", const TString &title="expCoeffa0")

void	initExpCoeffA1 (const double &value, const double &min, const double &max, const TString &name="expCoeffa1", const TString &title="expCoeffa1")

void	initExpCoeffA2 (const double &value, const double &min, const double &max, const TString &name="expCoeffa2", const TString &title="expCoeffa2")

void	initFGCB (const double &value, const double &min, const double &max, const TString &name="fGCB", const TString &title="fGCB")

void	initFsig (const double &value, const double &min, const double &max, const TString &name="fsig", const TString &title="signal fraction")

void	initGamma (const double &value, const double &min, const double &max, const TString &name="gamma", const TString &title="gamma")

void	initGaussFrac (const double &value, const double &min, const double &max, const TString &name="GaussFrac", const TString &title="GaussFrac")

void	initGaussFrac2 (const double &value, const double &min, const double &max, const TString &name="GaussFrac2", const TString &title="GaussFrac2")

void	initMean (const double &value, const double &min, const double &max, const TString &name="mean", const TString &title="mean")

void	initMean2 (const double &value, const double &min, const double &max, const TString &name="mean2", const TString &title="mean2")

void	initMean3 (const double &value, const double &min, const double &max, const TString &name="mean3", const TString &title="mean3")

void	initN (const double &value, const double &min, const double &max, const TString &name="n", const TString &title="n")

void	initSigma (const double &value, const double &min, const double &max, const TString &name="sigma", const TString &title="sigma")

void	initSigma2 (const double &value, const double &min, const double &max, const TString &name="sigma2", const TString &title="sigma2")

void	initSigma3 (const double &value, const double &min, const double &max, const TString &name="sigma3", const TString &title="sigma3")

RooRealVar *	mean ()

RooRealVar *	mean2 ()

RooRealVar *	mean3 ()

RooRealVar *	n ()

void	reinitializeParameters ()

RooRealVar *	sigma ()

RooRealVar *	sigma2 ()

RooRealVar *	sigma3 ()

Public Attributes
bool	useChi2_

Protected Attributes
RooRealVar *	a0_

RooRealVar *	a1_

RooRealVar *	a2_

RooRealVar *	a3_

RooRealVar *	a4_

RooRealVar *	a5_

RooRealVar *	a6_

RooRealVar *	alpha_

RooRealVar *	expCoeffa0_

RooRealVar *	expCoeffa1_

RooRealVar *	expCoeffa2_

RooRealVar *	fGCB_

RooRealVar *	fsig_

RooRealVar *	gamma_

RooRealVar *	gaussFrac2_

RooRealVar *	gaussFrac_

double	initVal_a0

double	initVal_a1

double	initVal_a2

double	initVal_a3

double	initVal_a4

double	initVal_a5

double	initVal_a6

double	initVal_alpha

double	initVal_expCoeffa0

double	initVal_expCoeffa1

double	initVal_expCoeffa2

double	initVal_fGCB

double	initVal_fsig

double	initVal_gamma

double	initVal_gaussFrac

double	initVal_gaussFrac2

double	initVal_mean

double	initVal_mean2

double	initVal_mean3

double	initVal_n

double	initVal_sigma

double	initVal_sigma2

double	initVal_sigma3

RooRealVar *	mean2_

RooRealVar *	mean3_

RooRealVar *	mean_

RooRealVar *	n_

RooRealVar *	sigma2_

RooRealVar *	sigma3_

RooRealVar *	sigma_

Detailed Description

Definition at line 51 of file FitWithRooFit.h.

Constructor & Destructor Documentation

◆ FitWithRooFit()

FitWithRooFit::FitWithRooFit ( )

inline

Definition at line 53 of file FitWithRooFit.h.

       : useChi2_(false),
         mean_(nullptr),
         mean2_(nullptr),
         mean3_(nullptr),
         sigma_(nullptr),
         sigma2_(nullptr),
         sigma3_(nullptr),
         gamma_(nullptr),
         gaussFrac_(nullptr),
         gaussFrac2_(nullptr),
         expCoeffa0_(nullptr),
         expCoeffa1_(nullptr),
         expCoeffa2_(nullptr),
         fsig_(nullptr),
         a0_(nullptr),
         a1_(nullptr),
         a2_(nullptr),
         a3_(nullptr),
         a4_(nullptr),
         a5_(nullptr),
         a6_(nullptr),
         alpha_(nullptr),
         n_(nullptr),
         fGCB_(nullptr) {}

Member Function Documentation

◆ a0()

RooRealVar* FitWithRooFit::a0 ( )

inline

Definition at line 236 of file FitWithRooFit.h.

References a0_.

236 { return a0_; }

FitWithRooFit::a0_

RooRealVar * a0_

Definition: FitWithRooFit.h:270

◆ a1()

RooRealVar* FitWithRooFit::a1 ( )

inline

Definition at line 237 of file FitWithRooFit.h.

References a1_.

237 { return a1_; }

FitWithRooFit::a1_

RooRealVar * a1_

Definition: FitWithRooFit.h:271

◆ a2()

RooRealVar* FitWithRooFit::a2 ( )

inline

Definition at line 238 of file FitWithRooFit.h.

References a2_.

238 { return a2_; }

FitWithRooFit::a2_

RooRealVar * a2_

Definition: FitWithRooFit.h:272

◆ a3()

RooRealVar* FitWithRooFit::a3 ( )

inline

Definition at line 239 of file FitWithRooFit.h.

References a3_.

239 { return a3_; }

FitWithRooFit::a3_

RooRealVar * a3_

Definition: FitWithRooFit.h:273

◆ a4()

RooRealVar* FitWithRooFit::a4 ( )

inline

Definition at line 240 of file FitWithRooFit.h.

References a4_.

240 { return a4_; }

FitWithRooFit::a4_

RooRealVar * a4_

Definition: FitWithRooFit.h:274

◆ a5()

RooRealVar* FitWithRooFit::a5 ( )

inline

Definition at line 241 of file FitWithRooFit.h.

References a5_.

241 { return a5_; }

FitWithRooFit::a5_

RooRealVar * a5_

Definition: FitWithRooFit.h:275

◆ a6()

RooRealVar* FitWithRooFit::a6 ( )

inline

Definition at line 242 of file FitWithRooFit.h.

References a6_.

242 { return a6_; }

FitWithRooFit::a6_

RooRealVar * a6_

Definition: FitWithRooFit.h:276

◆ alpha()

RooRealVar* FitWithRooFit::alpha ( )

inline

Definition at line 243 of file FitWithRooFit.h.

References alpha_.

243 { return alpha_; }

FitWithRooFit::alpha_

RooRealVar * alpha_

Definition: FitWithRooFit.h:277

◆ buildBackgroundModel()

RooAbsPdf * FitWithRooFit::buildBackgroundModel	(	RooRealVar *	x,
		const TString &	backgroundType
	)

Build the model for the specified background type.

Definition at line 489 of file FitWithRooFit.cc.

References a0_, a1_, a2_, a3_, a4_, a5_, a6_, beamvalidation::exit(), expCoeffa0_, expCoeffa1_, expCoeffa2_, fsig_, and x.

Referenced by buildModel().

                                                                                            {
   RooAbsPdf* background = nullptr;
   if (backgroundType == "exponential") {
     // Add an exponential for the background
     if ((expCoeffa1_ == nullptr) || (fsig_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize expCoeffa1 and fsig";
       exit(1);
     }
     background = new RooExponential("exponential", "exponential", *x, *expCoeffa1_);
   }
 
   if (backgroundType == "exponentialpol") {
     // Add an exponential for the background
     if ((expCoeffa0_ == nullptr) || (expCoeffa1_ == nullptr) || (expCoeffa2_ == nullptr) || (fsig_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize expCoeff and fsig";
       exit(1);
     }
     background = new RooGenericPdf("exponential",
                                    "exponential",
                                    "TMath::Exp(@1+@2*@0+@3*@0*@0)",
                                    RooArgList(*x, *expCoeffa0_, *expCoeffa1_, *expCoeffa2_));
   }
 
   else if (backgroundType == "chebychev0") {
     // Add a linear background
     if (a0_ == nullptr) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize a0";
       exit(1);
     }
     background = new RooChebychev("chebychev0", "chebychev0", *x, *a0_);
   } else if (backgroundType == "chebychev1") {
     // Add a 2nd order chebychev polynomial background
     if ((a0_ == nullptr) || (a1_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize a0 and a1";
       exit(1);
     }
     background = new RooChebychev("chebychev1", "chebychev1", *x, RooArgList(*a0_, *a1_));
   } else if (backgroundType == "chebychev3") {
     // Add a 3rd order chebychev polynomial background
     if ((a0_ == nullptr) || (a1_ == nullptr) || (a2_ == nullptr) || (a3_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize a0, a1, a2 and a3";
       exit(1);
     }
     background = new RooChebychev("3rdOrderPol", "3rdOrderPol", *x, RooArgList(*a0_, *a1_, *a2_, *a3_));
   }
 
   else if (backgroundType == "chebychev6") {
     // Add a 6th order chebychev polynomial background
     if ((a0_ == nullptr) || (a1_ == nullptr) || (a2_ == nullptr) || (a3_ == nullptr) || (a4_ == nullptr) ||
         (a5_ == nullptr) || (a6_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize a0, a1, a2, a3, "
              "a4, a5 and a6";
       exit(1);
     }
     background =
         new RooChebychev("6thOrderPol", "6thOrderPol", *x, RooArgList(*a0_, *a1_, *a2_, *a3_, *a4_, *a5_, *a6_));
   }
 
   return background;
 }

◆ buildModel()

RooAbsPdf * FitWithRooFit::buildModel	(	RooRealVar *	x,
		const TString &	signalType,
		const TString &	backgroundType
	)

Build the model to fit.

Definition at line 557 of file FitWithRooFit.cc.

References buildBackgroundModel(), buildSignalModel(), beamvalidation::exit(), fsig_, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, and x.

Referenced by fit().

                                                                                                             {
   RooAbsPdf* model = nullptr;
 
   RooAbsPdf* signal = buildSignalModel(x, signalType);
   RooAbsPdf* background = buildBackgroundModel(x, backgroundType);
 
   if ((signal != nullptr) && (background != nullptr)) {
     // Combine signal and background pdfs
     edm::LogPrint("FitWithRooFit") << "Building model with signal and backgound";
     model = new RooAddPdf("model", "model", RooArgList(*signal, *background), *fsig_);
   } else if (signal != nullptr) {
     edm::LogPrint("FitWithRooFit") << "Building model with signal";
     model = signal;
   } else if (background != nullptr) {
     edm::LogPrint("FitWithRooFit") << "Building model with backgound";
     model = background;
   } else {
     edm::LogWarning("FitWithRooFit") << "Nothing to fit";
     exit(0);
   }
   return model;
 }

◆ buildSignalModel()

RooAbsPdf * FitWithRooFit::buildSignalModel	(	RooRealVar *	x,
		const TString &	signalType
	)

Build the model for the specified signal type.

Definition at line 339 of file FitWithRooFit.cc.

References alpha_, beamvalidation::exit(), fGCB_, gamma_, gaussFrac2_, gaussFrac_, mean2_, mean3_, mean_, n_, sigma2_, sigma3_, sigma_, and x.

Referenced by buildModel().

                                                                                    {
   RooAbsPdf* signal = nullptr;
   if (signalType == "gaussian") {
     // Fit a Gaussian p.d.f to the data
     if ((mean_ == nullptr) || (sigma_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean and sigma";
       exit(1);
     }
     signal = new RooGaussian("gauss", "gauss", *x, *mean_, *sigma_);
   } else if (signalType == "doubleGaussian") {
     // Fit with double gaussian
     if ((mean_ == nullptr) || (sigma_ == nullptr) || (sigma2_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, sigma and sigma2";
       exit(1);
     }
     RooGaussModel* gaussModel = new RooGaussModel("gaussModel", "gaussModel", *x, *mean_, *sigma_);
     RooGaussModel* gaussModel2 = new RooGaussModel("gaussModel2", "gaussModel2", *x, *mean_, *sigma2_);
     signal = new RooAddModel("doubleGaussian", "double gaussian", RooArgList(*gaussModel, *gaussModel2), *gaussFrac_);
   } else if (signalType == "tripleGaussian") {
     // Fit with triple gaussian
     if ((mean_ == nullptr) || (mean2_ == nullptr) || (mean3_ == nullptr) || (sigma_ == nullptr) ||
         (sigma2_ == nullptr) || (sigma3_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, mean2, "
              "mean3, sigma, sigma2, sigma3";
       exit(1);
     }
     RooGaussModel* gaussModel = new RooGaussModel("gaussModel", "gaussModel", *x, *mean_, *sigma_);
     RooGaussModel* gaussModel2 = new RooGaussModel("gaussModel2", "gaussModel2", *x, *mean2_, *sigma2_);
     RooGaussModel* gaussModel3 = new RooGaussModel("gaussModel3", "gaussModel3", *x, *mean3_, *sigma3_);
     signal = new RooAddModel("tripleGaussian",
                              "triple gaussian",
                              RooArgList(*gaussModel, *gaussModel2, *gaussModel3),
                              RooArgList(*gaussFrac_, *gaussFrac2_));
   } else if (signalType == "breitWigner") {
     // Fit a Breit-Wigner
     if ((mean_ == nullptr) || (gamma_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean and gamma";
       exit(1);
     }
     signal = new RooBreitWigner("breiWign", "breitWign", *x, *mean_, *gamma_);
   } else if (signalType == "relBreitWigner") {
     // Fit a relativistic Breit-Wigner
     if ((mean_ == nullptr) || (gamma_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean and gamma";
       exit(1);
     }
     signal = new RooGenericPdf("Relativistic Breit-Wigner",
                                "RBW",
                                "@0/(pow(@0*@0 - @1*@1,2) + @2*@2*@0*@0*@0*@0/(@1*@1))",
                                RooArgList(*x, *mean_, *gamma_));
   } else if (signalType == "voigtian") {
     // Fit a Voigtian
     if ((mean_ == nullptr) || (sigma_ == nullptr) || (gamma_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, sigma and gamma";
       exit(1);
     }
     signal = new RooVoigtian("voigt", "voigt", *x, *mean_, *gamma_, *sigma_);
   } else if (signalType == "crystalBall") {
     // Fit a CrystalBall
     if ((mean_ == nullptr) || (sigma_ == nullptr) || (alpha_ == nullptr) || (n_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, sigma, "
              "alpha and n";
       exit(1);
     }
     signal = new RooCBShape("crystalBall", "crystalBall", *x, *mean_, *sigma_, *alpha_, *n_);
   } else if (signalType == "breitWignerTimesCB") {
     // Fit a Breit Wigner convoluted with a CrystalBall
     if ((mean_ == nullptr) || (mean2_ == nullptr) || (sigma_ == nullptr) || (gamma_ == nullptr) ||
         (alpha_ == nullptr) || (n_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, mean2, "
              "sigma, gamma, alpha and n";
       exit(1);
     }
     RooAbsPdf* bw = new RooBreitWigner("breiWigner", "breitWigner", *x, *mean_, *gamma_);
     RooAbsPdf* cb = new RooCBShape("crystalBall", "crystalBall", *x, *mean2_, *sigma_, *alpha_, *n_);
     signal = new RooFFTConvPdf("breitWignerTimesCB", "breitWignerTimesCB", *x, *bw, *cb);
   } else if (signalType == "relBreitWignerTimesCB") {
     // Fit a relativistic Breit Wigner convoluted with a CrystalBall
     if ((mean_ == nullptr) || (mean2_ == nullptr) || (sigma_ == nullptr) || (gamma_ == nullptr) ||
         (alpha_ == nullptr) || (n_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, mean2, "
              "sigma, gamma, alpha and n";
       exit(1);
     }
     RooGenericPdf* bw = new RooGenericPdf("Relativistic Breit-Wigner",
                                           "RBW",
                                           "@0/(pow(@0*@0 - @1*@1,2) + @2*@2*@0*@0*@0*@0/(@1*@1))",
                                           RooArgList(*x, *mean_, *gamma_));
     RooAbsPdf* cb = new RooCBShape("crystalBall", "crystalBall", *x, *mean2_, *sigma_, *alpha_, *n_);
     signal = new RooFFTConvPdf("relBreitWignerTimesCB", "relBreitWignerTimesCB", *x, *bw, *cb);
   } else if (signalType == "gaussianPlusCrystalBall") {
     // Fit a Gaussian + CrystalBall with the same mean
     if ((mean_ == nullptr) || (sigma_ == nullptr) || (alpha_ == nullptr) || (n_ == nullptr) || (sigma2_ == nullptr) ||
         (fGCB_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, sigma, "
              "sigma2, alpha, n and fGCB";
       exit(1);
     }
     RooAbsPdf* tempCB = new RooCBShape("crystalBall", "crystalBall", *x, *mean_, *sigma_, *alpha_, *n_);
     RooAbsPdf* tempGaussian = new RooGaussian("gauss", "gauss", *x, *mean_, *sigma2_);
 
     signal =
         new RooAddPdf("gaussianPlusCrystalBall", "gaussianPlusCrystalBall", RooArgList(*tempCB, *tempGaussian), *fGCB_);
   } else if (signalType == "voigtianPlusCrystalBall") {
     // Fit a Voigtian + CrystalBall with the same mean
     if ((mean_ == nullptr) || (sigma_ == nullptr) || (gamma_ == nullptr) || (alpha_ == nullptr) || (n_ == nullptr) ||
         (sigma2_ == nullptr) || (fGCB_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, gamma, "
              "sigma, sigma2, alpha, n and fGCB";
       exit(1);
     }
     RooAbsPdf* tempVoigt = new RooVoigtian("voigt", "voigt", *x, *mean_, *gamma_, *sigma_);
     RooAbsPdf* tempCB = new RooCBShape("crystalBall", "crystalBall", *x, *mean_, *sigma2_, *alpha_, *n_);
 
     signal =
         new RooAddPdf("voigtianPlusCrystalBall", "voigtianPlusCrystalBall", RooArgList(*tempCB, *tempVoigt), *fGCB_);
   } else if (signalType == "breitWignerPlusCrystalBall") {
     // Fit a Breit-Wigner + CrystalBall with the same mean
     if ((mean_ == nullptr) || (gamma_ == nullptr) || (alpha_ == nullptr) || (n_ == nullptr) || (sigma2_ == nullptr) ||
         (fGCB_ == nullptr)) {
       edm::LogError("FitWithRooFit")
           << "Error: one or more parameters are not initialized. Please be sure to initialize mean, gamma, "
              "sigma, alpha, n and fGCB";
       exit(1);
     }
     RooAbsPdf* tempBW = new RooBreitWigner("breitWign", "breitWign", *x, *mean_, *gamma_);
     RooAbsPdf* tempCB = new RooCBShape("crystalBall", "crystalBall", *x, *mean_, *sigma2_, *alpha_, *n_);
 
     signal =
         new RooAddPdf("breitWignerPlusCrystalBall", "breitWignerPlusCrystalBall", RooArgList(*tempCB, *tempBW), *fGCB_);
   }
 
   else if (signalType != "") {
     edm::LogError("FitWithRooFit") << "Unknown signal function: " << signalType << ". Signal will not be in the model";
   }
   return signal;
 }

◆ expCoeffa0()

RooRealVar* FitWithRooFit::expCoeffa0 ( )

inline

Definition at line 232 of file FitWithRooFit.h.

References expCoeffa0_.

232 { return expCoeffa0_; }

FitWithRooFit::expCoeffa0_

RooRealVar * expCoeffa0_

Definition: FitWithRooFit.h:266

◆ expCoeffa1()

RooRealVar* FitWithRooFit::expCoeffa1 ( )

inline

Definition at line 233 of file FitWithRooFit.h.

References expCoeffa1_.

233 { return expCoeffa1_; }

FitWithRooFit::expCoeffa1_

RooRealVar * expCoeffa1_

Definition: FitWithRooFit.h:267

◆ expCoeffa2()

RooRealVar* FitWithRooFit::expCoeffa2 ( )

inline

Definition at line 234 of file FitWithRooFit.h.

References expCoeffa2_.

234 { return expCoeffa2_; }

FitWithRooFit::expCoeffa2_

RooRealVar * expCoeffa2_

Definition: FitWithRooFit.h:268

◆ fGCB()

RooRealVar* FitWithRooFit::fGCB ( )

inline

Definition at line 245 of file FitWithRooFit.h.

References fGCB_.

245 { return fGCB_; }

FitWithRooFit::fGCB_

RooRealVar * fGCB_

Definition: FitWithRooFit.h:279

◆ fit()

void FitWithRooFit::fit	(	TH1 *	histo,
		const TString	signalType,
		const TString	backgroundType,
		const double &	xMin = `0.`,
		const double &	xMax = `0.`,
		bool	sumW2Error = `false`
	)

Definition at line 19 of file FitWithRooFit.cc.

References buildModel(), nano_mu_local_reco_cff::chi2, cuy::dh, timingPdfMaker::histo, importTH1(), DataMixerDataOnData_cff::Label, visualization-live-secondInstance_cfg::m, ReggeGribovPartonMC_EposLHC_2760GeV_PbPb_cfi::model, reinitializeParameters(), useChi2_, x, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

Referenced by trackingPlots.Iteration::modules(), and ZMassBinFit_OldTool().

                                          {
   reinitializeParameters();
 
   rooPair imported = importTH1(histo, xMin, xMax);
   RooRealVar x(imported.first);
   RooDataHist* dh = imported.second;
 
   // Make plot of binned dataset showing Poisson error bars (RooFit default)
   RooPlot* frame = x.frame(RooFit::Title("di-muon mass M(#mu^{+}#mu^{-}) [GeV]"));
   frame->SetName(TString(histo->GetName()) + "_frame");
   dh->plotOn(frame);
 
   // Build the composite model
   RooAbsPdf* model = buildModel(&x, signalType, backgroundType);
 
   std::unique_ptr<RooAbsReal> chi2{model->createChi2(*dh, RooFit::DataError(RooAbsData::SumW2))};
 
   // Fit the composite model
   // -----------------------
   // Fit with likelihood
   if (!useChi2_) {
     if (sumW2Error)
       model->fitTo(*dh, RooFit::Save(), RooFit::SumW2Error(kTRUE));
     else
       model->fitTo(*dh);
   }
   // Fit with chi^2
   else {
     edm::LogWarning("FitWithRooFit") << "FITTING WITH CHI^2";
     RooMinimizer m(*chi2);
     m.migrad();
     m.hesse();
     // RooFitResult* r_chi2_wgt = m.save();
   }
   model->plotOn(frame, RooFit::LineColor(kRed));
   model->plotOn(frame, RooFit::Components(backgroundType), RooFit::LineStyle(kDashed));
   model->paramOn(frame,
                  RooFit::Label("fit result"),
                  RooFit::Layout(0.65, 0.90, 0.90),
                  RooFit::Format("NEU", RooFit::AutoPrecision(2)));
 
   // TODO: fix next lines to get the prob(chi2) (ndof should be dynamically set according to the choosen pdf)
   // double chi2 = xframe.chiSquare("model","data",ndof);
   // double ndoff = xframeGetNbinsX();
   // double chi2prob = TMath::Prob(chi2,ndoff);
 
   // P l o t   a n d   f i t   a   R o o D a t a H i s t   w i t h   i n t e r n a l   e r r o r s
   // ---------------------------------------------------------------------------------------------
 
   // If histogram has custom error (i.e. its contents is does not originate from a Poisson process
   // but e.g. is a sum of weighted events) you can data with symmetric 'sum-of-weights' error instead
   // (same error bars as shown by ROOT)
   RooPlot* frame2 = x.frame(RooFit::Title("di-muon mass M(#mu^{+}#mu^{-}) [GeV]"));
   dh->plotOn(frame2, RooFit::DataError(RooAbsData::SumW2));
   model->plotOn(frame2, RooFit::LineColor(kRed));
   model->plotOn(frame2, RooFit::Components(backgroundType), RooFit::LineStyle(kDashed));
   model->paramOn(frame2,
                  RooFit::Label("fit result"),
                  RooFit::Layout(0.65, 0.90, 0.90),
                  RooFit::Format("NEU", RooFit::AutoPrecision(2)));
 
   // Please note that error bars shown (Poisson or SumW2) are for visualization only, the are NOT used
   // in a maximum likelihood fit
   //
   // A (binned) ML fit will ALWAYS assume the Poisson error interpretation of data (the mathematical definition
   // of likelihood does not take any external definition of errors). Data with non-unit weights can only be correctly
   // fitted with a chi^2 fit (see rf602_chi2fit.C)
 
   // Draw all frames on a canvas
   // if( canvas == 0 ) {
   //   canvas = new TCanvas("rf102_dataimport","rf102_dataimport",800,800);
   //   canvas->cd(1);
   // }
   // canvas->Divide(2,1);
   // canvas->cd(1) ; frame->Draw();
   // canvas->cd(2) ; frame2->Draw();
 
   frame2->Draw();
 }

◆ fsig()

RooRealVar* FitWithRooFit::fsig ( )

inline

Definition at line 235 of file FitWithRooFit.h.

References fsig_.

235 { return fsig_; }

FitWithRooFit::fsig_

RooRealVar * fsig_

Definition: FitWithRooFit.h:269

◆ gamma()

RooRealVar* FitWithRooFit::gamma ( )

inline

Definition at line 229 of file FitWithRooFit.h.

References gamma_.

Referenced by ZMassBinFit_OldTool().

229 { return gamma_; }

FitWithRooFit::gamma_

RooRealVar * gamma_

Definition: FitWithRooFit.h:263

◆ gaussFrac()

RooRealVar* FitWithRooFit::gaussFrac ( )

inline

Definition at line 230 of file FitWithRooFit.h.

References gaussFrac_.

230 { return gaussFrac_; }

FitWithRooFit::gaussFrac_

RooRealVar * gaussFrac_

Definition: FitWithRooFit.h:264

◆ gaussFrac2()

RooRealVar* FitWithRooFit::gaussFrac2 ( )

inline

Definition at line 231 of file FitWithRooFit.h.

References gaussFrac2_.

231 { return gaussFrac2_; }

FitWithRooFit::gaussFrac2_

RooRealVar * gaussFrac2_

Definition: FitWithRooFit.h:265

◆ importTH1()

rooPair FitWithRooFit::importTH1	(	TH1 *	histo,
		const double &	inputXmin,
		const double &	inputXmax
	)

Definition at line 5 of file FitWithRooFit.cc.

References timingPdfMaker::histo, x, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

Referenced by fit().

                                                                                              {
   double xMin = inputXmin;
   double xMax = inputXmax;
   if ((xMin == xMax) && xMin == 0) {
     xMin = histo->GetXaxis()->GetXmin();
     xMax = histo->GetXaxis()->GetXmax();
   }
   // Declare observable x
   RooRealVar x("InvMass", "di-muon mass M(#mu^{+}#mu^{-}) [GeV]", xMin, xMax);
   // Create a binned dataset that imports contents of TH1 and associates its contents to observable 'x'
   return (std::make_pair(x, new RooDataHist("dh", "dh", x, RooFit::Import(*histo))));
 }

◆ initA0()

void FitWithRooFit::initA0	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a0"`,
		const TString &	title = `"a0"`
	)

Definition at line 209 of file FitWithRooFit.cc.

References a0_, initVal_a0, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a0_ != nullptr)
     delete a0_;
   a0_ = new RooRealVar(name, title, value, min, max);
   initVal_a0 = value;
 }

◆ initA1()

void FitWithRooFit::initA1	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a1"`,
		const TString &	title = `"a1"`
	)

Definition at line 217 of file FitWithRooFit.cc.

References a1_, initVal_a1, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a1_ != nullptr)
     delete a1_;
   a1_ = new RooRealVar(name, title, value, min, max);
   initVal_a1 = value;
 }

◆ initA2()

void FitWithRooFit::initA2	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a2"`,
		const TString &	title = `"a2"`
	)

Definition at line 225 of file FitWithRooFit.cc.

References a2_, initVal_a2, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a2_ != nullptr)
     delete a2_;
   a2_ = new RooRealVar(name, title, value, min, max);
   initVal_a2 = value;
 }

◆ initA3()

void FitWithRooFit::initA3	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a3"`,
		const TString &	title = `"a3"`
	)

Definition at line 233 of file FitWithRooFit.cc.

References a3_, initVal_a3, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a3_ != nullptr)
     delete a3_;
   a3_ = new RooRealVar(name, title, value, min, max);
   initVal_a3 = value;
 }

◆ initA4()

void FitWithRooFit::initA4	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a4"`,
		const TString &	title = `"a4"`
	)

Definition at line 241 of file FitWithRooFit.cc.

References a4_, initVal_a4, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a4_ != nullptr)
     delete a4_;
   a4_ = new RooRealVar(name, title, value, min, max);
   initVal_a4 = value;
 }

◆ initA5()

void FitWithRooFit::initA5	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a5"`,
		const TString &	title = `"a5"`
	)

Definition at line 249 of file FitWithRooFit.cc.

References a5_, initVal_a5, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a5_ != nullptr)
     delete a5_;
   a5_ = new RooRealVar(name, title, value, min, max);
   initVal_a5 = value;
 }

◆ initA6()

void FitWithRooFit::initA6	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"a6"`,
		const TString &	title = `"a6"`
	)

Definition at line 257 of file FitWithRooFit.cc.

References a6_, initVal_a6, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (a6_ != nullptr)
     delete a6_;
   a6_ = new RooRealVar(name, title, value, min, max);
   initVal_a6 = value;
 }

◆ initAlpha()

void FitWithRooFit::initAlpha	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"alpha"`,
		const TString &	title = `"alpha"`
	)

Definition at line 265 of file FitWithRooFit.cc.

References alpha_, initVal_alpha, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (alpha_ != nullptr)
     delete alpha_;
   alpha_ = new RooRealVar(name, title, value, min, max);
   initVal_alpha = value;
 }

◆ initExpCoeffA0()

void FitWithRooFit::initExpCoeffA0	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"expCoeffa0"`,
		const TString &	title = `"expCoeffa0"`
	)

Definition at line 177 of file FitWithRooFit.cc.

References expCoeffa0_, initVal_expCoeffa0, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (expCoeffa0_ != nullptr)
     delete expCoeffa0_;
   expCoeffa0_ = new RooRealVar(name, title, value, min, max);
   initVal_expCoeffa0 = value;
 }

◆ initExpCoeffA1()

void FitWithRooFit::initExpCoeffA1	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"expCoeffa1"`,
		const TString &	title = `"expCoeffa1"`
	)

Definition at line 185 of file FitWithRooFit.cc.

References expCoeffa1_, initVal_expCoeffa1, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (expCoeffa1_ != nullptr)
     delete expCoeffa1_;
   expCoeffa1_ = new RooRealVar(name, title, value, min, max);
   initVal_expCoeffa1 = value;
 }

◆ initExpCoeffA2()

void FitWithRooFit::initExpCoeffA2	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"expCoeffa2"`,
		const TString &	title = `"expCoeffa2"`
	)

Definition at line 193 of file FitWithRooFit.cc.

References expCoeffa2_, initVal_expCoeffa2, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (expCoeffa2_ != nullptr)
     delete expCoeffa2_;
   expCoeffa2_ = new RooRealVar(name, title, value, min, max);
   initVal_expCoeffa2 = value;
 }

◆ initFGCB()

void FitWithRooFit::initFGCB	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"fGCB"`,
		const TString &	title = `"fGCB"`
	)

Definition at line 281 of file FitWithRooFit.cc.

References fGCB_, initVal_fGCB, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (fGCB_ != nullptr)
     delete fGCB_;
   fGCB_ = new RooRealVar(name, title, value, min, max);
   initVal_fGCB = value;
 }

◆ initFsig()

void FitWithRooFit::initFsig	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"fsig"`,
		const TString &	title = `"signal fraction"`
	)

Definition at line 201 of file FitWithRooFit.cc.

References fsig_, initVal_fsig, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (fsig_ != nullptr)
     delete fsig_;
   fsig_ = new RooRealVar(name, title, value, min, max);
   initVal_fsig = value;
 }

◆ initGamma()

void FitWithRooFit::initGamma	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"gamma"`,
		const TString &	title = `"gamma"`
	)

Definition at line 153 of file FitWithRooFit.cc.

References gamma_, initVal_gamma, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (gamma_ != nullptr)
     delete gamma_;
   gamma_ = new RooRealVar(name, title, value, min, max);
   initVal_gamma = value;
 }

◆ initGaussFrac()

void FitWithRooFit::initGaussFrac	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"GaussFrac"`,
		const TString &	title = `"GaussFrac"`
	)

Definition at line 161 of file FitWithRooFit.cc.

References gaussFrac_, initVal_gaussFrac, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

                                                                                                           {
   if (gaussFrac_ != nullptr)
     delete gaussFrac_;
   gaussFrac_ = new RooRealVar(name, title, value, min, max);
   initVal_gaussFrac = value;
 }

◆ initGaussFrac2()

void FitWithRooFit::initGaussFrac2	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"GaussFrac2"`,
		const TString &	title = `"GaussFrac2"`
	)

Definition at line 169 of file FitWithRooFit.cc.

References gaussFrac2_, initVal_gaussFrac2, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

                                                                                                           {
   if (gaussFrac2_ != nullptr)
     delete gaussFrac2_;
   gaussFrac2_ = new RooRealVar(name, title, value, min, max);
   initVal_gaussFrac2 = value;
 }

◆ initMean()

void FitWithRooFit::initMean	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"mean"`,
		const TString &	title = `"mean"`
	)

Definition at line 105 of file FitWithRooFit.cc.

References initVal_mean, SiStripPI::max, mean_, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (mean_ != nullptr)
     delete mean_;
   mean_ = new RooRealVar(name, title, value, min, max);
   initVal_mean = value;
 }

◆ initMean2()

void FitWithRooFit::initMean2	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"mean2"`,
		const TString &	title = `"mean2"`
	)

Definition at line 113 of file FitWithRooFit.cc.

References initVal_mean2, SiStripPI::max, mean2_, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (mean2_ != nullptr)
     delete mean2_;
   mean2_ = new RooRealVar(name, title, value, min, max);
   initVal_mean2 = value;
 }

◆ initMean3()

void FitWithRooFit::initMean3	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"mean3"`,
		const TString &	title = `"mean3"`
	)

Definition at line 121 of file FitWithRooFit.cc.

References initVal_mean3, SiStripPI::max, mean3_, SiStripPI::min, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

                                                                                                           {
   if (mean3_ != nullptr)
     delete mean3_;
   mean3_ = new RooRealVar(name, title, value, min, max);
   initVal_mean3 = value;
 }

◆ initN()

void FitWithRooFit::initN	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"n"`,
		const TString &	title = `"n"`
	)

Definition at line 273 of file FitWithRooFit.cc.

References initVal_n, SiStripPI::max, SiStripPI::min, n_, Skims_PA_cff::name, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (n_ != nullptr)
     delete n_;
   n_ = new RooRealVar(name, title, value, min, max);
   initVal_n = value;
 }

◆ initSigma()

void FitWithRooFit::initSigma	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"sigma"`,
		const TString &	title = `"sigma"`
	)

Definition at line 129 of file FitWithRooFit.cc.

References initVal_sigma, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, sigma_, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (sigma_ != nullptr)
     delete sigma_;
   sigma_ = new RooRealVar(name, title, value, min, max);
   initVal_sigma = value;
 }

◆ initSigma2()

void FitWithRooFit::initSigma2	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"sigma2"`,
		const TString &	title = `"sigma2"`
	)

Definition at line 137 of file FitWithRooFit.cc.

References initVal_sigma2, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, sigma2_, runGCPTkAlMap::title, and relativeConstraints::value.

Referenced by ZMassBinFit_OldTool().

                                                                                                           {
   if (sigma2_ != nullptr)
     delete sigma2_;
   sigma2_ = new RooRealVar(name, title, value, min, max);
   initVal_sigma2 = value;
 }

◆ initSigma3()

void FitWithRooFit::initSigma3	(	const double &	value,
		const double &	min,
		const double &	max,
		const TString &	name = `"sigma3"`,
		const TString &	title = `"sigma3"`
	)

Definition at line 145 of file FitWithRooFit.cc.

References initVal_sigma3, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, sigma3_, runGCPTkAlMap::title, and relativeConstraints::value.

                                                                                                           {
   if (sigma3_ != nullptr)
     delete sigma3_;
   sigma3_ = new RooRealVar(name, title, value, min, max);
   initVal_sigma3 = value;
 }

◆ mean()

RooRealVar* FitWithRooFit::mean ( )

inline

Definition at line 223 of file FitWithRooFit.h.

References mean_.

Referenced by ZMassBinFit_OldTool().

223 { return mean_; }

FitWithRooFit::mean_

RooRealVar * mean_

Definition: FitWithRooFit.h:257

◆ mean2()

RooRealVar* FitWithRooFit::mean2 ( )

inline

Definition at line 224 of file FitWithRooFit.h.

References mean2_.

Referenced by ZMassBinFit_OldTool().

224 { return mean2_; }

FitWithRooFit::mean2_

RooRealVar * mean2_

Definition: FitWithRooFit.h:258

◆ mean3()

RooRealVar* FitWithRooFit::mean3 ( )

inline

Definition at line 225 of file FitWithRooFit.h.

References mean3_.

225 { return mean3_; }

FitWithRooFit::mean3_

RooRealVar * mean3_

Definition: FitWithRooFit.h:259

◆ n()

RooRealVar* FitWithRooFit::n ( )

inline

Definition at line 244 of file FitWithRooFit.h.

References n_.

244 { return n_; }

FitWithRooFit::n_

RooRealVar * n_

Definition: FitWithRooFit.h:278

◆ reinitializeParameters()

void FitWithRooFit::reinitializeParameters ( )

Definition at line 289 of file FitWithRooFit.cc.

References a0_, a1_, a2_, a3_, a4_, a5_, a6_, alpha_, expCoeffa0_, expCoeffa1_, expCoeffa2_, fGCB_, fsig_, gamma_, gaussFrac2_, gaussFrac_, initVal_a0, initVal_a1, initVal_a2, initVal_a3, initVal_a4, initVal_a5, initVal_a6, initVal_alpha, initVal_expCoeffa0, initVal_expCoeffa1, initVal_expCoeffa2, initVal_fGCB, initVal_fsig, initVal_gamma, initVal_gaussFrac, initVal_gaussFrac2, initVal_mean, initVal_mean2, initVal_mean3, initVal_n, initVal_sigma, initVal_sigma2, initVal_sigma3, mean2_, mean3_, mean_, n_, sigma2_, sigma3_, and sigma_.

Referenced by fit().

                                            {
   if (mean_ != nullptr)
     mean_->setVal(initVal_mean);
   if (mean2_ != nullptr)
     mean2_->setVal(initVal_mean2);
   if (mean3_ != nullptr)
     mean3_->setVal(initVal_mean3);
   if (sigma_ != nullptr)
     sigma_->setVal(initVal_sigma);
   if (sigma2_ != nullptr)
     sigma2_->setVal(initVal_sigma2);
   if (sigma3_ != nullptr)
     sigma3_->setVal(initVal_sigma3);
   if (gamma_ != nullptr)
     gamma_->setVal(initVal_gamma);
   if (gaussFrac_ != nullptr)
     gaussFrac_->setVal(initVal_gaussFrac);
   if (gaussFrac2_ != nullptr)
     gaussFrac2_->setVal(initVal_gaussFrac2);
   if (expCoeffa0_ != nullptr)
     expCoeffa0_->setVal(initVal_expCoeffa0);
   if (expCoeffa1_ != nullptr)
     expCoeffa1_->setVal(initVal_expCoeffa1);
   if (expCoeffa2_ != nullptr)
     expCoeffa2_->setVal(initVal_expCoeffa2);
   if (fsig_ != nullptr)
     fsig_->setVal(initVal_fsig);
   if (a0_ != nullptr)
     a0_->setVal(initVal_a0);
   if (a1_ != nullptr)
     a1_->setVal(initVal_a1);
   if (a2_ != nullptr)
     a2_->setVal(initVal_a2);
   if (a3_ != nullptr)
     a3_->setVal(initVal_a3);
   if (a4_ != nullptr)
     a4_->setVal(initVal_a4);
   if (a5_ != nullptr)
     a5_->setVal(initVal_a5);
   if (a6_ != nullptr)
     a6_->setVal(initVal_a6);
   if (alpha_ != nullptr)
     alpha_->setVal(initVal_alpha);
   if (n_ != nullptr)
     n_->setVal(initVal_n);
   if (fGCB_ != nullptr)
     fGCB_->setVal(initVal_fGCB);
 }

◆ sigma()

RooRealVar* FitWithRooFit::sigma ( )

inline

Definition at line 226 of file FitWithRooFit.h.

References sigma_.

Referenced by ZMassBinFit_OldTool().

226 { return sigma_; }

FitWithRooFit::sigma_

RooRealVar * sigma_

Definition: FitWithRooFit.h:260

◆ sigma2()

RooRealVar* FitWithRooFit::sigma2 ( )

inline

Definition at line 227 of file FitWithRooFit.h.

References sigma2_.

227 { return sigma2_; }

FitWithRooFit::sigma2_

RooRealVar * sigma2_

Definition: FitWithRooFit.h:261

◆ sigma3()

RooRealVar* FitWithRooFit::sigma3 ( )

inline

Definition at line 228 of file FitWithRooFit.h.

References sigma3_.

228 { return sigma3_; }

FitWithRooFit::sigma3_

RooRealVar * sigma3_

Definition: FitWithRooFit.h:262

Member Data Documentation

◆ a0_

RooRealVar* FitWithRooFit::a0_

protected

Definition at line 270 of file FitWithRooFit.h.

Referenced by a0(), buildBackgroundModel(), initA0(), and reinitializeParameters().

◆ a1_

RooRealVar* FitWithRooFit::a1_

protected

Definition at line 271 of file FitWithRooFit.h.

Referenced by a1(), buildBackgroundModel(), initA1(), and reinitializeParameters().

◆ a2_

RooRealVar* FitWithRooFit::a2_

protected

Definition at line 272 of file FitWithRooFit.h.

Referenced by a2(), buildBackgroundModel(), initA2(), and reinitializeParameters().

◆ a3_

RooRealVar* FitWithRooFit::a3_

protected

Definition at line 273 of file FitWithRooFit.h.

Referenced by a3(), buildBackgroundModel(), initA3(), and reinitializeParameters().

◆ a4_

RooRealVar* FitWithRooFit::a4_

protected

Definition at line 274 of file FitWithRooFit.h.

Referenced by a4(), buildBackgroundModel(), initA4(), and reinitializeParameters().

◆ a5_

RooRealVar* FitWithRooFit::a5_

protected

Definition at line 275 of file FitWithRooFit.h.

Referenced by a5(), buildBackgroundModel(), initA5(), and reinitializeParameters().

◆ a6_

RooRealVar* FitWithRooFit::a6_

protected

Definition at line 276 of file FitWithRooFit.h.

Referenced by a6(), buildBackgroundModel(), initA6(), and reinitializeParameters().

◆ alpha_

RooRealVar* FitWithRooFit::alpha_

protected

Definition at line 277 of file FitWithRooFit.h.

Referenced by alpha(), buildSignalModel(), initAlpha(), and reinitializeParameters().

◆ expCoeffa0_

RooRealVar* FitWithRooFit::expCoeffa0_

protected

Definition at line 266 of file FitWithRooFit.h.

Referenced by buildBackgroundModel(), expCoeffa0(), initExpCoeffA0(), and reinitializeParameters().

◆ expCoeffa1_

RooRealVar* FitWithRooFit::expCoeffa1_

protected

Definition at line 267 of file FitWithRooFit.h.

Referenced by buildBackgroundModel(), expCoeffa1(), initExpCoeffA1(), and reinitializeParameters().

◆ expCoeffa2_

RooRealVar* FitWithRooFit::expCoeffa2_

protected

Definition at line 268 of file FitWithRooFit.h.

Referenced by buildBackgroundModel(), expCoeffa2(), initExpCoeffA2(), and reinitializeParameters().

◆ fGCB_

RooRealVar* FitWithRooFit::fGCB_

protected

Definition at line 279 of file FitWithRooFit.h.

Referenced by buildSignalModel(), fGCB(), initFGCB(), and reinitializeParameters().

◆ fsig_

RooRealVar* FitWithRooFit::fsig_

protected

Definition at line 269 of file FitWithRooFit.h.

Referenced by buildBackgroundModel(), buildModel(), fsig(), initFsig(), and reinitializeParameters().

◆ gamma_

RooRealVar* FitWithRooFit::gamma_

protected

Definition at line 263 of file FitWithRooFit.h.

Referenced by buildSignalModel(), gamma(), initGamma(), and reinitializeParameters().

◆ gaussFrac2_

RooRealVar* FitWithRooFit::gaussFrac2_

protected

Definition at line 265 of file FitWithRooFit.h.

Referenced by buildSignalModel(), gaussFrac2(), initGaussFrac2(), and reinitializeParameters().

◆ gaussFrac_

RooRealVar* FitWithRooFit::gaussFrac_

protected

Definition at line 264 of file FitWithRooFit.h.

Referenced by buildSignalModel(), gaussFrac(), initGaussFrac(), and reinitializeParameters().

◆ initVal_a0

double FitWithRooFit::initVal_a0

protected

Definition at line 295 of file FitWithRooFit.h.

Referenced by initA0(), and reinitializeParameters().

◆ initVal_a1

double FitWithRooFit::initVal_a1

protected

Definition at line 296 of file FitWithRooFit.h.

Referenced by initA1(), and reinitializeParameters().

◆ initVal_a2

double FitWithRooFit::initVal_a2

protected

Definition at line 297 of file FitWithRooFit.h.

Referenced by initA2(), and reinitializeParameters().

◆ initVal_a3

double FitWithRooFit::initVal_a3

protected

Definition at line 298 of file FitWithRooFit.h.

Referenced by initA3(), and reinitializeParameters().

◆ initVal_a4

double FitWithRooFit::initVal_a4

protected

Definition at line 299 of file FitWithRooFit.h.

Referenced by initA4(), and reinitializeParameters().

◆ initVal_a5

double FitWithRooFit::initVal_a5

protected

Definition at line 300 of file FitWithRooFit.h.

Referenced by initA5(), and reinitializeParameters().

◆ initVal_a6

double FitWithRooFit::initVal_a6

protected

Definition at line 301 of file FitWithRooFit.h.

Referenced by initA6(), and reinitializeParameters().

◆ initVal_alpha

double FitWithRooFit::initVal_alpha

protected

Definition at line 302 of file FitWithRooFit.h.

Referenced by initAlpha(), and reinitializeParameters().

◆ initVal_expCoeffa0

double FitWithRooFit::initVal_expCoeffa0

protected

Definition at line 291 of file FitWithRooFit.h.

Referenced by initExpCoeffA0(), and reinitializeParameters().

◆ initVal_expCoeffa1

double FitWithRooFit::initVal_expCoeffa1

protected

Definition at line 292 of file FitWithRooFit.h.

Referenced by initExpCoeffA1(), and reinitializeParameters().

◆ initVal_expCoeffa2

double FitWithRooFit::initVal_expCoeffa2

protected

Definition at line 293 of file FitWithRooFit.h.

Referenced by initExpCoeffA2(), and reinitializeParameters().

◆ initVal_fGCB

double FitWithRooFit::initVal_fGCB

protected

Definition at line 304 of file FitWithRooFit.h.

Referenced by initFGCB(), and reinitializeParameters().

◆ initVal_fsig

double FitWithRooFit::initVal_fsig

protected

Definition at line 294 of file FitWithRooFit.h.

Referenced by initFsig(), and reinitializeParameters().

◆ initVal_gamma

double FitWithRooFit::initVal_gamma

protected

Definition at line 288 of file FitWithRooFit.h.

Referenced by initGamma(), and reinitializeParameters().

◆ initVal_gaussFrac

double FitWithRooFit::initVal_gaussFrac

protected

Definition at line 289 of file FitWithRooFit.h.

Referenced by initGaussFrac(), and reinitializeParameters().

◆ initVal_gaussFrac2

double FitWithRooFit::initVal_gaussFrac2

protected

Definition at line 290 of file FitWithRooFit.h.

Referenced by initGaussFrac2(), and reinitializeParameters().

◆ initVal_mean

double FitWithRooFit::initVal_mean

protected

Definition at line 282 of file FitWithRooFit.h.

Referenced by initMean(), and reinitializeParameters().

◆ initVal_mean2

double FitWithRooFit::initVal_mean2

protected

Definition at line 283 of file FitWithRooFit.h.

Referenced by initMean2(), and reinitializeParameters().

◆ initVal_mean3

double FitWithRooFit::initVal_mean3

protected

Definition at line 284 of file FitWithRooFit.h.

Referenced by initMean3(), and reinitializeParameters().

◆ initVal_n

double FitWithRooFit::initVal_n

protected

Definition at line 303 of file FitWithRooFit.h.

Referenced by initN(), and reinitializeParameters().

◆ initVal_sigma

double FitWithRooFit::initVal_sigma

protected

Definition at line 285 of file FitWithRooFit.h.

Referenced by initSigma(), and reinitializeParameters().

◆ initVal_sigma2

double FitWithRooFit::initVal_sigma2

protected

Definition at line 286 of file FitWithRooFit.h.

Referenced by initSigma2(), and reinitializeParameters().

◆ initVal_sigma3

double FitWithRooFit::initVal_sigma3

protected

Definition at line 287 of file FitWithRooFit.h.

Referenced by initSigma3(), and reinitializeParameters().

◆ mean2_

RooRealVar* FitWithRooFit::mean2_

protected

Definition at line 258 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initMean2(), mean2(), and reinitializeParameters().

◆ mean3_

RooRealVar* FitWithRooFit::mean3_

protected

Definition at line 259 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initMean3(), mean3(), and reinitializeParameters().

◆ mean_

RooRealVar* FitWithRooFit::mean_

protected

Definition at line 257 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initMean(), mean(), and reinitializeParameters().

◆ n_

RooRealVar* FitWithRooFit::n_

protected

Definition at line 278 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initN(), n(), and reinitializeParameters().

◆ sigma2_

RooRealVar* FitWithRooFit::sigma2_

protected

Definition at line 261 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initSigma2(), reinitializeParameters(), and sigma2().

◆ sigma3_

RooRealVar* FitWithRooFit::sigma3_

protected

Definition at line 262 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initSigma3(), reinitializeParameters(), and sigma3().

◆ sigma_

RooRealVar* FitWithRooFit::sigma_

protected

Definition at line 260 of file FitWithRooFit.h.

Referenced by buildSignalModel(), initSigma(), reinitializeParameters(), and sigma().

◆ useChi2_

bool FitWithRooFit::useChi2_

Definition at line 253 of file FitWithRooFit.h.

Referenced by fit(), and ZMassBinFit_OldTool().

Public Member Functions

Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FitWithRooFit()

Member Function Documentation

◆ a0()

◆ a1()

◆ a2()

◆ a3()

◆ a4()

◆ a5()

◆ a6()

◆ alpha()

◆ buildBackgroundModel()

◆ buildModel()

◆ buildSignalModel()

◆ expCoeffa0()

◆ expCoeffa1()

◆ expCoeffa2()

◆ fGCB()

◆ fit()

◆ fsig()

◆ gamma()

◆ gaussFrac()

◆ gaussFrac2()

◆ importTH1()

◆ initA0()

◆ initA1()

◆ initA2()

◆ initA3()

◆ initA4()

◆ initA5()

◆ initA6()

◆ initAlpha()

◆ initExpCoeffA0()

◆ initExpCoeffA1()

◆ initExpCoeffA2()

◆ initFGCB()

◆ initFsig()

◆ initGamma()

◆ initGaussFrac()

◆ initGaussFrac2()

◆ initMean()

◆ initMean2()

◆ initMean3()

◆ initN()

◆ initSigma()

◆ initSigma2()

◆ initSigma3()

◆ mean()

◆ mean2()

◆ mean3()

◆ n()

◆ reinitializeParameters()

◆ sigma()

◆ sigma2()

◆ sigma3()

Member Data Documentation

◆ a0_

◆ a1_

◆ a2_

◆ a3_

◆ a4_

◆ a5_

◆ a6_

◆ alpha_

◆ expCoeffa0_

◆ expCoeffa1_

◆ expCoeffa2_

◆ fGCB_

◆ fsig_

◆ gamma_

◆ gaussFrac2_

◆ gaussFrac_

◆ initVal_a0

◆ initVal_a1

◆ initVal_a2

◆ initVal_a3