FitWithRooFit Class Reference

#include <FitWithRooFit.h>

Public Member Functions
RooRealVar *	a0 ()
RooRealVar *	a1 ()
RooRealVar *	a2 ()
RooRealVar *	a3 ()
RooRealVar *	a4 ()
RooRealVar *	a5 ()
RooRealVar *	a6 ()
RooRealVar *	alpha ()
std::unique_ptr< RooAbsPdf >	buildBackgroundModel (RooRealVar *x, const TString &backgroundType)
	Build the model for the specified background type. More...
std::unique_ptr< RooAbsPdf >	buildModel (RooRealVar *x, const TString &signalType, const TString &backgroundType)
	Build the model to fit. More...
std::unique_ptr< 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, double xMin=0., double xMax=0., bool sumW2Error=false)
RooRealVar *	fsig ()
RooRealVar *	gamma ()
RooRealVar *	gaussFrac ()
RooRealVar *	gaussFrac2 ()
std::unique_ptr< RooDataHist >	importTH1 (TH1 *histo, double xMin, double xMax)
void	initA0 (double value, double min, double max, const TString &name="a0", const TString &title="a0")
void	initA1 (double value, double min, double max, const TString &name="a1", const TString &title="a1")
void	initA2 (double value, double min, double max, const TString &name="a2", const TString &title="a2")
void	initA3 (double value, double min, double max, const TString &name="a3", const TString &title="a3")
void	initA4 (double value, double min, double max, const TString &name="a4", const TString &title="a4")
void	initA5 (double value, double min, double max, const TString &name="a5", const TString &title="a5")
void	initA6 (double value, double min, double max, const TString &name="a6", const TString &title="a6")
void	initAlpha (double value, double min, double max, const TString &name="alpha", const TString &title="alpha")
void	initExpCoeffA0 (double value, double min, double max, const TString &name="expCoeffa0", const TString &title="expCoeffa0")
void	initExpCoeffA1 (double value, double min, double max, const TString &name="expCoeffa1", const TString &title="expCoeffa1")
void	initExpCoeffA2 (double value, double min, double max, const TString &name="expCoeffa2", const TString &title="expCoeffa2")
void	initFGCB (double value, double min, double max, const TString &name="fGCB", const TString &title="fGCB")
void	initFsig (double value, double min, double max, const TString &name="fsig", const TString &title="signal fraction")
void	initGamma (double value, double min, double max, const TString &name="gamma", const TString &title="gamma")
void	initGaussFrac (double value, double min, double max, const TString &name="GaussFrac", const TString &title="GaussFrac")
void	initGaussFrac2 (double value, double min, double max, const TString &name="GaussFrac2", const TString &title="GaussFrac2")
void	initMean (double value, double min, double max, const TString &name="mean", const TString &title="mean")
void	initMean2 (double value, double min, double max, const TString &name="mean2", const TString &title="mean2")
void	initMean3 (double value, double min, double max, const TString &name="mean3", const TString &title="mean3")
void	initN (double value, double min, double max, const TString &name="n", const TString &title="n")
void	initSigma (double value, double min, double max, const TString &name="sigma", const TString &title="sigma")
void	initSigma2 (double value, double min, double max, const TString &name="sigma2", const TString &title="sigma2")
void	initSigma3 (double value, double min, 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_ = false

Protected Attributes
std::unique_ptr< RooRealVar >	a0_
std::unique_ptr< RooRealVar >	a1_
std::unique_ptr< RooRealVar >	a2_
std::unique_ptr< RooRealVar >	a3_
std::unique_ptr< RooRealVar >	a4_
std::unique_ptr< RooRealVar >	a5_
std::unique_ptr< RooRealVar >	a6_
std::unique_ptr< RooRealVar >	alpha_
std::unique_ptr< RooRealVar >	expCoeffa0_
std::unique_ptr< RooRealVar >	expCoeffa1_
std::unique_ptr< RooRealVar >	expCoeffa2_
std::unique_ptr< RooRealVar >	fGCB_
std::unique_ptr< RooRealVar >	fsig_
std::unique_ptr< RooRealVar >	gamma_
std::unique_ptr< RooRealVar >	gaussFrac2_
std::unique_ptr< RooRealVar >	gaussFrac_
double	initVal_a0 = 0.0
double	initVal_a1 = 0.0
double	initVal_a2 = 0.0
double	initVal_a3 = 0.0
double	initVal_a4 = 0.0
double	initVal_a5 = 0.0
double	initVal_a6 = 0.0
double	initVal_alpha = 0.0
double	initVal_expCoeffa0 = 0.0
double	initVal_expCoeffa1 = 0.0
double	initVal_expCoeffa2 = 0.0
double	initVal_fGCB = 0.0
double	initVal_fsig = 0.0
double	initVal_gamma = 0.0
double	initVal_gaussFrac = 0.0
double	initVal_gaussFrac2 = 0.0
double	initVal_mean = 0.0
double	initVal_mean2 = 0.0
double	initVal_mean3 = 0.0
double	initVal_n = 0.0
double	initVal_sigma = 0.0
double	initVal_sigma2 = 0.0
double	initVal_sigma3 = 0.0
std::unique_ptr< RooRealVar >	mean2_
std::unique_ptr< RooRealVar >	mean3_
std::unique_ptr< RooRealVar >	mean_
std::unique_ptr< RooRealVar >	n_
std::unique_ptr< RooRealVar >	sigma2_
std::unique_ptr< RooRealVar >	sigma3_
std::unique_ptr< RooRealVar >	sigma_

Detailed Description

This class allows to use RooFit to perform a fit on a TH1 histogram.
The currently implemented functions are:

• signal:
  – gaussian
  – double gaussian
  – voigtian
  
• background:
  – exponential
  It is possible to select any combination of signal and background.
  The fit() method receives the TH1, two strings specifying the signal and background type and the min and max x of the histogram over which to perform the fit.
  The variables of the fit must be initialized separately. For example when doing a gaussian+exponential fit the initMean, initSigma, initFSig and initExpCoeff methods must be used to create and initialize the corresponding variables.
  The methods names after the variables return the fit result.

Definition at line 46 of file FitWithRooFit.h.

Member Function Documentation

a0()

RooRealVar* FitWithRooFit::a0 ( )

inline

Definition at line 100 of file FitWithRooFit.h.

References a0_.

{ return a0_.get(); }

a1()

RooRealVar* FitWithRooFit::a1 ( )

inline

Definition at line 101 of file FitWithRooFit.h.

References a1_.

{ return a1_.get(); }

a2()

RooRealVar* FitWithRooFit::a2 ( )

inline

Definition at line 102 of file FitWithRooFit.h.

References a2_.

{ return a2_.get(); }

a3()

RooRealVar* FitWithRooFit::a3 ( )

inline

Definition at line 103 of file FitWithRooFit.h.

References a3_.

{ return a3_.get(); }

a4()

RooRealVar* FitWithRooFit::a4 ( )

inline

Definition at line 104 of file FitWithRooFit.h.

References a4_.

{ return a4_.get(); }

a5()

RooRealVar* FitWithRooFit::a5 ( )

inline

Definition at line 105 of file FitWithRooFit.h.

References a5_.

{ return a5_.get(); }

a6()

RooRealVar* FitWithRooFit::a6 ( )

inline

Definition at line 106 of file FitWithRooFit.h.

References a6_.

{ return a6_.get(); }

alpha()

RooRealVar* FitWithRooFit::alpha ( )

inline

Definition at line 107 of file FitWithRooFit.h.

References alpha_.

{ return alpha_.get(); }

buildBackgroundModel()

std::unique_ptr< RooAbsPdf > FitWithRooFit::buildBackgroundModel	(	RooRealVar *	x,
		const TString &	backgroundType
	)

Build the model for the specified background type.

Definition at line 415 of file FitWithRooFit.cc.

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

Referenced by buildModel().

                                                                                                         {
  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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<RooChebychev>(
        "6thOrderPol", "6thOrderPol", *x, RooArgList{*a0_, *a1_, *a2_, *a3_, *a4_, *a5_, *a6_});
  }
  return nullptr;
}

buildModel()

std::unique_ptr< RooAbsPdf > FitWithRooFit::buildModel	(	RooRealVar *	x,
		const TString &	signalType,
		const TString &	backgroundType
	)

Build the model to fit.

Definition at line 481 of file FitWithRooFit.cc.

References buildBackgroundModel(), buildSignalModel(), makeMuonMisalignmentScenario::components, beamvalidation::exit(), fsig_, isotrackApplyRegressor::model, eostools::move(), and x.

Referenced by fit().

                                                                                    {
  std::unique_ptr<RooAbsPdf> model;

  std::unique_ptr<RooAbsPdf> signal = buildSignalModel(x, signalType);
  std::unique_ptr<RooAbsPdf> background = buildBackgroundModel(x, backgroundType);

  if ((signal != nullptr) && (background != nullptr)) {
    // Combine signal and background pdfs
    edm::LogPrint("FitWithRooFit") << "Building model with signal and backgound";
    RooArgList components{*signal.release(), *background.release()};
    model = std::make_unique<RooAddPdf>("model", "model", components, *fsig_);
    model->addOwnedComponents(components);
  } else if (signal != nullptr) {
    edm::LogPrint("FitWithRooFit") << "Building model with signal";
    model = std::move(signal);
  } else if (background != nullptr) {
    edm::LogPrint("FitWithRooFit") << "Building model with backgound";
    model = std::move(background);
  } else {
    edm::LogWarning("FitWithRooFit") << "Nothing to fit";
    exit(0);
  }
  return model;
}

buildSignalModel()

std::unique_ptr< RooAbsPdf > FitWithRooFit::buildSignalModel	(	RooRealVar *	x,
		const TString &	signalType
	)

Build the model for the specified signal type.

Definition at line 251 of file FitWithRooFit.cc.

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

Referenced by buildModel().

                                                                                                 {
  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);
    }
    return std::make_unique<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_);
    RooArgList components{*gaussModel, *gaussModel2};
    auto out = std::make_unique<RooAddModel>("doubleGaussian", "double gaussian", components, *gaussFrac_);
    out->addOwnedComponents(components);
    return out;
  } 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_);
    RooArgList components{*gaussModel, *gaussModel2, *gaussModel3};
    auto out = std::make_unique<RooAddModel>(
        "tripleGaussian", "triple gaussian", components, RooArgList{*gaussFrac_, *gaussFrac2_});
    out->addOwnedComponents(components);
    return out;
  } 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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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);
    }
    return std::make_unique<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_);
    auto out = std::make_unique<RooFFTConvPdf>("breitWignerTimesCB", "breitWignerTimesCB", *x, *bw, *cb);
    out->addOwnedComponents({*bw, *cb});
    return out;
  } 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))",
                                          {*x, *mean_, *gamma_});
    RooAbsPdf* cb = new RooCBShape("crystalBall", "crystalBall", *x, *mean2_, *sigma_, *alpha_, *n_);
    auto out = std::make_unique<RooFFTConvPdf>("relBreitWignerTimesCB", "relBreitWignerTimesCB", *x, *bw, *cb);
    out->addOwnedComponents({*bw, *cb});
    return out;
  } 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_);
    RooArgList components{*tempCB, *tempGaussian};

    auto out = std::make_unique<RooAddPdf>("gaussianPlusCrystalBall", "gaussianPlusCrystalBall", components, *fGCB_);
    out->addOwnedComponents(components);
    return out;
  } 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_);
    RooArgList components{*tempVoigt, *tempCB};

    auto out = std::make_unique<RooAddPdf>("voigtianPlusCrystalBall", "voigtianPlusCrystalBall", components, *fGCB_);
    out->addOwnedComponents(components);
    return out;
  } 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_);
    RooArgList components{*tempCB, *tempBW};

    auto out =
        std::make_unique<RooAddPdf>("breitWignerPlusCrystalBall", "breitWignerPlusCrystalBall", components, *fGCB_);
    out->addOwnedComponents(components);
    return out;
  }

  else if (signalType != "") {
    edm::LogError("FitWithRooFit") << "Unknown signal function: " << signalType << ". Signal will not be in the model";
  }
  return nullptr;
}

expCoeffa0()

RooRealVar* FitWithRooFit::expCoeffa0 ( )

inline

Definition at line 96 of file FitWithRooFit.h.

References expCoeffa0_.

{ return expCoeffa0_.get(); }

expCoeffa1()

RooRealVar* FitWithRooFit::expCoeffa1 ( )

inline

Definition at line 97 of file FitWithRooFit.h.

References expCoeffa1_.

{ return expCoeffa1_.get(); }

expCoeffa2()

RooRealVar* FitWithRooFit::expCoeffa2 ( )

inline

Definition at line 98 of file FitWithRooFit.h.

References expCoeffa2_.

{ return expCoeffa2_.get(); }

fGCB()

RooRealVar* FitWithRooFit::fGCB ( )

inline

Definition at line 109 of file FitWithRooFit.h.

References fGCB_.

{ return fGCB_.get(); }

fit()

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

Definition at line 18 of file FitWithRooFit.cc.

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

Referenced by trackingPlots.Iteration::modules().

                                                                                                                   {
  reinitializeParameters();

  std::unique_ptr<RooDataHist> dh = importTH1(histo, xMin, xMax);
  RooRealVar x(*static_cast<RooRealVar*>(dh->get()->find("InvMass")));

  // 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
  std::unique_ptr<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_) {
    // use  RooFit::PrintLevel(-1) to reduce output to screen
    model->fitTo(*dh, RooFit::SumW2Error(sumW2Error), RooFit::PrintLevel(-1));
  }
  // 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));

  // this causes segfaults when the fit causes a Warning in <ROOT::Math::Fitter::CalculateHessErrors>: Error when calculating Hessian
  //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));

  // this causes segfaults when the fit causes a Warning in <ROOT::Math::Fitter::CalculateHessErrors>: Error when calculating Hessian
  //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 99 of file FitWithRooFit.h.

References fsig_.

{ return fsig_.get(); }

gamma()

RooRealVar* FitWithRooFit::gamma ( )

inline

Definition at line 93 of file FitWithRooFit.h.

References gamma_.

{ return gamma_.get(); }

gaussFrac()

RooRealVar* FitWithRooFit::gaussFrac ( )

inline

Definition at line 94 of file FitWithRooFit.h.

References gaussFrac_.

{ return gaussFrac_.get(); }

gaussFrac2()

RooRealVar* FitWithRooFit::gaussFrac2 ( )

inline

Definition at line 95 of file FitWithRooFit.h.

References gaussFrac2_.

{ return gaussFrac2_.get(); }

importTH1()

std::unique_ptr< RooDataHist > FitWithRooFit::importTH1	(	TH1 *	histo,
		double	xMin,
		double	xMax
	)

Definition at line 6 of file FitWithRooFit.cc.

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

Referenced by fit().

                                                                                        {
  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_unique<RooDataHist>("dh", "dh", x, RooFit::Import(*histo));
}

initA0()

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

Definition at line 169 of file FitWithRooFit.cc.

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

                                                                                                          {
  a0_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a0 = value;
}

initA1()

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

Definition at line 174 of file FitWithRooFit.cc.

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

                                                                                                          {
  a1_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a1 = value;
}

initA2()

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

Definition at line 179 of file FitWithRooFit.cc.

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

                                                                                                          {
  a2_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a2 = value;
}

initA3()

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

Definition at line 184 of file FitWithRooFit.cc.

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

                                                                                                          {
  a3_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a3 = value;
}

initA4()

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

Definition at line 189 of file FitWithRooFit.cc.

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

                                                                                                          {
  a4_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a4 = value;
}

initA5()

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

Definition at line 194 of file FitWithRooFit.cc.

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

                                                                                                          {
  a5_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a5 = value;
}

initA6()

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

Definition at line 199 of file FitWithRooFit.cc.

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

                                                                                                          {
  a6_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_a6 = value;
}

initAlpha()

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

Definition at line 204 of file FitWithRooFit.cc.

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

                                                                                                             {
  alpha_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_alpha = value;
}

initExpCoeffA0()

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

Definition at line 149 of file FitWithRooFit.cc.

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

                                                                                                                  {
  expCoeffa0_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_expCoeffa0 = value;
}

initExpCoeffA1()

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

Definition at line 154 of file FitWithRooFit.cc.

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

                                                                                                                  {
  expCoeffa1_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_expCoeffa1 = value;
}

initExpCoeffA2()

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

Definition at line 159 of file FitWithRooFit.cc.

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

                                                                                                                  {
  expCoeffa2_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_expCoeffa2 = value;
}

initFGCB()

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

Definition at line 214 of file FitWithRooFit.cc.

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

                                                                                                            {
  fGCB_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_fGCB = value;
}

initFsig()

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

Definition at line 164 of file FitWithRooFit.cc.

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

                                                                                                            {
  fsig_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_fsig = value;
}

initGamma()

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

Definition at line 134 of file FitWithRooFit.cc.

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

                                                                                                             {
  gamma_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_gamma = value;
}

initGaussFrac()

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

Definition at line 139 of file FitWithRooFit.cc.

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

                                                                                                                 {
  gaussFrac_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_gaussFrac = value;
}

initGaussFrac2()

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

Definition at line 144 of file FitWithRooFit.cc.

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

                                                                                                                  {
  gaussFrac2_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_gaussFrac2 = value;
}

initMean()

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

Definition at line 104 of file FitWithRooFit.cc.

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

                                                                                                            {
  mean_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_mean = value;
}

initMean2()

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

Definition at line 109 of file FitWithRooFit.cc.

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

                                                                                                             {
  mean2_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_mean2 = value;
}

initMean3()

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

Definition at line 114 of file FitWithRooFit.cc.

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

                                                                                                             {
  mean3_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_mean3 = value;
}

initN()

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

Definition at line 209 of file FitWithRooFit.cc.

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

                                                                                                         {
  n_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_n = value;
}

initSigma()

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

Definition at line 119 of file FitWithRooFit.cc.

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

                                                                                                             {
  sigma_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_sigma = value;
}

initSigma2()

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

Definition at line 124 of file FitWithRooFit.cc.

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

                                                                                                              {
  sigma2_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_sigma2 = value;
}

initSigma3()

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

Definition at line 129 of file FitWithRooFit.cc.

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

                                                                                                              {
  sigma3_ = std::make_unique<RooRealVar>(name, title, value, min, max);
  initVal_sigma3 = value;
}

mean()

RooRealVar* FitWithRooFit::mean ( )

inline

Definition at line 87 of file FitWithRooFit.h.

References mean_.

{ return mean_.get(); }

mean2()

RooRealVar* FitWithRooFit::mean2 ( )

inline

Definition at line 88 of file FitWithRooFit.h.

References mean2_.

{ return mean2_.get(); }

mean3()

RooRealVar* FitWithRooFit::mean3 ( )

inline

Definition at line 89 of file FitWithRooFit.h.

References mean3_.

{ return mean3_.get(); }

n()

RooRealVar* FitWithRooFit::n ( )

inline

Definition at line 108 of file FitWithRooFit.h.

References n_.

{ return n_.get(); }

reinitializeParameters()

void FitWithRooFit::reinitializeParameters

(

)

Definition at line 219 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_, sigma_, heppy_batch::val, and ALCARECOEcalPhiSym_cff::var.

Referenced by fit().

                                           {
  auto initParam = [](std::unique_ptr<RooRealVar>& var, double val) {
    if (var)
      var->setVal(val);
  };

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

sigma()

RooRealVar* FitWithRooFit::sigma ( )

inline

Definition at line 90 of file FitWithRooFit.h.

References sigma_.

{ return sigma_.get(); }

sigma2()

RooRealVar* FitWithRooFit::sigma2 ( )

inline

Definition at line 91 of file FitWithRooFit.h.

References sigma2_.

{ return sigma2_.get(); }

sigma3()

RooRealVar* FitWithRooFit::sigma3 ( )

inline

Definition at line 92 of file FitWithRooFit.h.

References sigma3_.

{ return sigma3_.get(); }

Member Data Documentation

a0_

std::unique_ptr<RooRealVar> FitWithRooFit::a0_

protected

Definition at line 134 of file FitWithRooFit.h.

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

a1_

std::unique_ptr<RooRealVar> FitWithRooFit::a1_

protected

Definition at line 135 of file FitWithRooFit.h.

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

a2_

std::unique_ptr<RooRealVar> FitWithRooFit::a2_

protected

Definition at line 136 of file FitWithRooFit.h.

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

a3_

std::unique_ptr<RooRealVar> FitWithRooFit::a3_

protected

Definition at line 137 of file FitWithRooFit.h.

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

a4_

std::unique_ptr<RooRealVar> FitWithRooFit::a4_

protected

Definition at line 138 of file FitWithRooFit.h.

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

a5_

std::unique_ptr<RooRealVar> FitWithRooFit::a5_

protected

Definition at line 139 of file FitWithRooFit.h.

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

a6_

std::unique_ptr<RooRealVar> FitWithRooFit::a6_

protected

Definition at line 140 of file FitWithRooFit.h.

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

alpha_

std::unique_ptr<RooRealVar> FitWithRooFit::alpha_

protected

Definition at line 141 of file FitWithRooFit.h.

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

expCoeffa0_

std::unique_ptr<RooRealVar> FitWithRooFit::expCoeffa0_

protected

Definition at line 130 of file FitWithRooFit.h.

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

expCoeffa1_

std::unique_ptr<RooRealVar> FitWithRooFit::expCoeffa1_

protected

Definition at line 131 of file FitWithRooFit.h.

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

expCoeffa2_

std::unique_ptr<RooRealVar> FitWithRooFit::expCoeffa2_

protected

Definition at line 132 of file FitWithRooFit.h.

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

fGCB_

std::unique_ptr<RooRealVar> FitWithRooFit::fGCB_

protected

Definition at line 143 of file FitWithRooFit.h.

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

fsig_

std::unique_ptr<RooRealVar> FitWithRooFit::fsig_

protected

Definition at line 133 of file FitWithRooFit.h.

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

gamma_

std::unique_ptr<RooRealVar> FitWithRooFit::gamma_

protected

Definition at line 127 of file FitWithRooFit.h.

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

gaussFrac2_

std::unique_ptr<RooRealVar> FitWithRooFit::gaussFrac2_

protected

Definition at line 129 of file FitWithRooFit.h.

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

gaussFrac_

std::unique_ptr<RooRealVar> FitWithRooFit::gaussFrac_

protected

Definition at line 128 of file FitWithRooFit.h.

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

initVal_a0

double FitWithRooFit::initVal_a0 = 0.0

protected

Definition at line 159 of file FitWithRooFit.h.

Referenced by initA0(), and reinitializeParameters().

initVal_a1

double FitWithRooFit::initVal_a1 = 0.0

protected

Definition at line 160 of file FitWithRooFit.h.

Referenced by initA1(), and reinitializeParameters().

initVal_a2

double FitWithRooFit::initVal_a2 = 0.0

protected

Definition at line 161 of file FitWithRooFit.h.

Referenced by initA2(), and reinitializeParameters().

initVal_a3

double FitWithRooFit::initVal_a3 = 0.0

protected

Definition at line 162 of file FitWithRooFit.h.

Referenced by initA3(), and reinitializeParameters().

initVal_a4

double FitWithRooFit::initVal_a4 = 0.0

protected

Definition at line 163 of file FitWithRooFit.h.

Referenced by initA4(), and reinitializeParameters().

initVal_a5

double FitWithRooFit::initVal_a5 = 0.0

protected

Definition at line 164 of file FitWithRooFit.h.

Referenced by initA5(), and reinitializeParameters().

initVal_a6

double FitWithRooFit::initVal_a6 = 0.0

protected

Definition at line 165 of file FitWithRooFit.h.

Referenced by initA6(), and reinitializeParameters().

initVal_alpha

double FitWithRooFit::initVal_alpha = 0.0

protected

Definition at line 166 of file FitWithRooFit.h.

Referenced by initAlpha(), and reinitializeParameters().

initVal_expCoeffa0

double FitWithRooFit::initVal_expCoeffa0 = 0.0

protected

Definition at line 155 of file FitWithRooFit.h.

Referenced by initExpCoeffA0(), and reinitializeParameters().

initVal_expCoeffa1

double FitWithRooFit::initVal_expCoeffa1 = 0.0

protected

Definition at line 156 of file FitWithRooFit.h.

Referenced by initExpCoeffA1(), and reinitializeParameters().

initVal_expCoeffa2

double FitWithRooFit::initVal_expCoeffa2 = 0.0

protected

Definition at line 157 of file FitWithRooFit.h.

Referenced by initExpCoeffA2(), and reinitializeParameters().

initVal_fGCB

double FitWithRooFit::initVal_fGCB = 0.0

protected

Definition at line 168 of file FitWithRooFit.h.

Referenced by initFGCB(), and reinitializeParameters().

initVal_fsig

double FitWithRooFit::initVal_fsig = 0.0

protected

Definition at line 158 of file FitWithRooFit.h.

Referenced by initFsig(), and reinitializeParameters().

initVal_gamma

double FitWithRooFit::initVal_gamma = 0.0

protected

Definition at line 152 of file FitWithRooFit.h.

Referenced by initGamma(), and reinitializeParameters().

initVal_gaussFrac

double FitWithRooFit::initVal_gaussFrac = 0.0

protected

Definition at line 153 of file FitWithRooFit.h.

Referenced by initGaussFrac(), and reinitializeParameters().

initVal_gaussFrac2

double FitWithRooFit::initVal_gaussFrac2 = 0.0

protected

Definition at line 154 of file FitWithRooFit.h.

Referenced by initGaussFrac2(), and reinitializeParameters().

initVal_mean

double FitWithRooFit::initVal_mean = 0.0

protected

Definition at line 146 of file FitWithRooFit.h.

Referenced by initMean(), and reinitializeParameters().

initVal_mean2

double FitWithRooFit::initVal_mean2 = 0.0

protected

Definition at line 147 of file FitWithRooFit.h.

Referenced by initMean2(), and reinitializeParameters().

initVal_mean3

double FitWithRooFit::initVal_mean3 = 0.0

protected

Definition at line 148 of file FitWithRooFit.h.

Referenced by initMean3(), and reinitializeParameters().

initVal_n

double FitWithRooFit::initVal_n = 0.0

protected

Definition at line 167 of file FitWithRooFit.h.

Referenced by initN(), and reinitializeParameters().

initVal_sigma

double FitWithRooFit::initVal_sigma = 0.0

protected

Definition at line 149 of file FitWithRooFit.h.

Referenced by initSigma(), and reinitializeParameters().

initVal_sigma2

double FitWithRooFit::initVal_sigma2 = 0.0

protected

Definition at line 150 of file FitWithRooFit.h.

Referenced by initSigma2(), and reinitializeParameters().

initVal_sigma3

double FitWithRooFit::initVal_sigma3 = 0.0

protected

Definition at line 151 of file FitWithRooFit.h.

Referenced by initSigma3(), and reinitializeParameters().

mean2_

std::unique_ptr<RooRealVar> FitWithRooFit::mean2_

protected

Definition at line 122 of file FitWithRooFit.h.

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

mean3_

std::unique_ptr<RooRealVar> FitWithRooFit::mean3_

protected

Definition at line 123 of file FitWithRooFit.h.

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

mean_

std::unique_ptr<RooRealVar> FitWithRooFit::mean_

protected

Definition at line 121 of file FitWithRooFit.h.

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

n_

std::unique_ptr<RooRealVar> FitWithRooFit::n_

protected

Definition at line 142 of file FitWithRooFit.h.

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

sigma2_

std::unique_ptr<RooRealVar> FitWithRooFit::sigma2_

protected

Definition at line 125 of file FitWithRooFit.h.

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

sigma3_

std::unique_ptr<RooRealVar> FitWithRooFit::sigma3_

protected

Definition at line 126 of file FitWithRooFit.h.

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

sigma_

std::unique_ptr<RooRealVar> FitWithRooFit::sigma_

protected

Definition at line 124 of file FitWithRooFit.h.

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

useChi2_

bool FitWithRooFit::useChi2_ = false

Definition at line 117 of file FitWithRooFit.h.

Referenced by fit().