dqmTnP::VoigtianPlusExponentialFitter Class Reference

#include <GenericTnPFitter.h>

Inheritance diagram for dqmTnP::VoigtianPlusExponentialFitter:

Public Member Functions
void	fit (TH1 *pass, TH1 *all) override
void	setup (double expectedMean_, double massLow, double massHigh, double expectedSigma_, double width_)
	VoigtianPlusExponentialFitter (bool verbose=false)
	~VoigtianPlusExponentialFitter () override=default
Public Member Functions inherited from dqmTnP::AbstractFitter
	AbstractFitter (bool verbose_=false)
TString	calculateEfficiency (TH3 *pass, TH3 *all, int massDimension, TProfile2D *&eff, TProfile2D *&effChi2, const TString &plotName="")
TString	calculateEfficiency (TH2 *pass, TH2 *all, int massDimension, TProfile *&eff, TProfile *&effChi2, const TString &plotName="")
double	getChi2 ()
double	getEfficiency ()
double	getEfficiencyError ()
void	savePlot (const TString &name)
void	setup (double expectedMean_, double massLow, double massHigh, double expectedSigma_)
virtual	~AbstractFitter ()=default

Protected Attributes
RooExponential	exponentialFail
RooExponential	exponentialPass
RooAddPdf	pdfFail
RooAddPdf	pdfPass
RooRealVar	slopeFail
RooRealVar	slopePass
RooVoigtian	voigtian
RooRealVar	width
Protected Attributes inherited from dqmTnP::AbstractFitter
RooCategory	category
double	chi2
RooDataHist *	data
RooRealVar	efficiency
double	expectedMean
double	expectedSigma
RooRealVar	mass
RooRealVar	mean
RooRealVar	nBackgroundFail
RooRealVar	nBackgroundPass
RooRealVar	nSignalAll
RooFormulaVar	nSignalFail
RooFormulaVar	nSignalPass
RooRealVar	sigma
RooSimultaneous	simPdf
bool	verbose

Detailed Description

Definition at line 334 of file GenericTnPFitter.h.

Constructor & Destructor Documentation

VoigtianPlusExponentialFitter()

dqmTnP::VoigtianPlusExponentialFitter::VoigtianPlusExponentialFitter ( bool  verbose = false )

inline

Definition at line 346 of file GenericTnPFitter.h.

References dqmTnP::AbstractFitter::simPdf, and ApeEstimator_cff::width.

        : AbstractFitter(verbose),
          width("width", "width", 2.5, "GeV"),
          voigtian("voigtian", "voigtian", mass, mean, width, sigma),
          slopeFail("slopeFail", "slopeFail", 0., -1., 0.),
          exponentialFail("linearFail", "linearFail", mass, slopeFail),
          slopePass("slopePass", "slopePass", 0., -1., 0.),
          exponentialPass("linearPass", "linearPass", mass, slopePass),
          pdfFail(
              "pdfFail", "pdfFail", RooArgList(voigtian, exponentialFail), RooArgList(nSignalFail, nBackgroundFail)),
          pdfPass(
              "pdfPass", "pdfPass", RooArgList(voigtian, exponentialPass), RooArgList(nSignalPass, nBackgroundPass)) {
      width.setConstant(kTRUE);
      simPdf.addPdf(pdfFail, "fail");
      simPdf.addPdf(pdfPass, "pass");
    };

~VoigtianPlusExponentialFitter()

dqmTnP::VoigtianPlusExponentialFitter::~VoigtianPlusExponentialFitter ( )

overridedefault

Member Function Documentation

fit()

void dqmTnP::VoigtianPlusExponentialFitter::fit ( TH1 *  pass, TH1 *  all  )

inlineoverridevirtual

Implements dqmTnP::AbstractFitter.

Definition at line 371 of file GenericTnPFitter.h.

References python.cmstools::all(), dqmTnP::AbstractFitter::category, dqmTnP::AbstractFitter::chi2, dqmTnP::AbstractFitter::data, dqmTnP::AbstractFitter::efficiency, dqmTnP::AbstractFitter::expectedMean, dqmTnP::AbstractFitter::expectedSigma, dqmTnP::AbstractFitter::mass, dqmTnP::AbstractFitter::mean, dqmTnP::AbstractFitter::nBackgroundFail, dqmTnP::AbstractFitter::nBackgroundPass, dqmTnP::AbstractFitter::nSignalAll, dqmTnP::AbstractFitter::sigma, dqmTnP::AbstractFitter::simPdf, EcalTangentSkim_cfg::Verbose, and dqmTnP::AbstractFitter::verbose.

Referenced by trackingPlots.Iteration::modules().

                                           {
      using namespace RooFit;
      all->Add(pass, -1);
      TH1*& fail = all;
      if (!data)
        delete data;
      data = new RooDataHist("data", "data", mass, Index(category), Import("fail", *fail), Import("pass", *pass));
      if (pass->Integral() + fail->Integral() < 5) {
        efficiency.setVal(0.5);
        efficiency.setError(0.5);
        chi2 = 0;
        return;
      }
      mean.setVal(expectedMean);
      sigma.setVal(expectedSigma);
      efficiency.setVal(pass->Integral() / (pass->Integral() + fail->Integral()));
      nSignalAll.setVal(0.5 * (fail->Integral() + pass->Integral()));
      nBackgroundFail.setVal(0.5 * fail->Integral());
      nBackgroundPass.setVal(0.5 * pass->Integral());
      slopeFail.setVal(0.);
      slopePass.setVal(0.);
      if (verbose) {
        simPdf.fitTo(*data);
      } else {
        simPdf.fitTo(*data, Verbose(kFALSE), PrintLevel(-1), Warnings(kFALSE), PrintEvalErrors(-1));
      }
      RooDataHist dataFail("fail", "fail", mass, fail);
      RooDataHist dataPass("pass", "pass", mass, pass);
      using AbsRealPtr = std::unique_ptr<RooAbsReal>;
      const double chi2Fail = AbsRealPtr(pdfFail.createChi2(dataFail, DataError(RooAbsData::Poisson)))->getVal();
      const double chi2Pass = AbsRealPtr(pdfPass.createChi2(dataPass, DataError(RooAbsData::Poisson)))->getVal();
      chi2 = (chi2Fail + chi2Pass) / (2 * all->GetNbinsX() - 8);
      if (chi2 > 3) {
        efficiency.setVal(0.5);
        efficiency.setError(0.5);
      }
    }

setup()

void dqmTnP::VoigtianPlusExponentialFitter::setup ( double  expectedMean_, double  massLow, double  massHigh, double  expectedSigma_, double  width_  )

inline

Definition at line 364 of file GenericTnPFitter.h.

References dqmTnP::AbstractFitter::expectedMean, dqmTnP::AbstractFitter::expectedSigma, dqmTnP::AbstractFitter::mass, EgHLTOffHistBins_cfi::massHigh, dqmTnP::AbstractFitter::mean, and ApeEstimator_cff::width.

Referenced by DQMGenericTnPClient::calculateEfficiency().

                                                                                                            {
      expectedMean = expectedMean_;
      expectedSigma = expectedSigma_;
      mass.setRange(massLow, massHigh);
      mean.setRange(massLow, massHigh);
      width.setVal(width_);
    }

Member Data Documentation

exponentialFail

RooExponential dqmTnP::VoigtianPlusExponentialFitter::exponentialFail

protected

Definition at line 339 of file GenericTnPFitter.h.

exponentialPass

RooExponential dqmTnP::VoigtianPlusExponentialFitter::exponentialPass

protected

Definition at line 341 of file GenericTnPFitter.h.

pdfFail

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfFail

protected

Definition at line 342 of file GenericTnPFitter.h.

pdfPass

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfPass

protected

Definition at line 343 of file GenericTnPFitter.h.

slopeFail

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopeFail

protected

Definition at line 338 of file GenericTnPFitter.h.

slopePass

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopePass

protected

Definition at line 340 of file GenericTnPFitter.h.

voigtian

RooVoigtian dqmTnP::VoigtianPlusExponentialFitter::voigtian

protected

Definition at line 337 of file GenericTnPFitter.h.

width

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::width

protected

Definition at line 336 of file GenericTnPFitter.h.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.root2matplotlib.Hist::bar(), python.rootplot.root2matplotlib.Hist::barh(), and python.rootplot.utilities.Hist::delete_bin().