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

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

inline

Definition at line 346 of file GenericTnPFitter.h.

References dqmTnP::AbstractFitter::simPdf.

        : 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");
    };

dqmTnP::VoigtianPlusExponentialFitter::~VoigtianPlusExponentialFitter ( )

overridedefault

Member Function Documentation

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);
      chi2 = (RooChi2Var("chi2Fail", "chi2Fail", pdfFail, dataFail, DataError(RooAbsData::Poisson)).getVal() +
              RooChi2Var("chi2Pass", "chi2Pass", pdfPass, dataPass, DataError(RooAbsData::Poisson)).getVal()) /
             (2 * all->GetNbinsX() - 8);
      if (chi2 > 3) {
        efficiency.setVal(0.5);
        efficiency.setError(0.5);
      }
    }

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, and dqmTnP::AbstractFitter::mean.

Referenced by DQMGenericTnPClient::calculateEfficiency().

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

Member Data Documentation

RooExponential dqmTnP::VoigtianPlusExponentialFitter::exponentialFail

protected

Definition at line 339 of file GenericTnPFitter.h.

RooExponential dqmTnP::VoigtianPlusExponentialFitter::exponentialPass

protected

Definition at line 341 of file GenericTnPFitter.h.

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfFail

protected

Definition at line 342 of file GenericTnPFitter.h.

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfPass

protected

Definition at line 343 of file GenericTnPFitter.h.

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopeFail

protected

Definition at line 338 of file GenericTnPFitter.h.

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopePass

protected

Definition at line 340 of file GenericTnPFitter.h.

RooVoigtian dqmTnP::VoigtianPlusExponentialFitter::voigtian

protected

Definition at line 337 of file GenericTnPFitter.h.

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(), and python.rootplot.root2matplotlib.Hist::barh().