dqmTnP::VoigtianPlusExponentialFitter Class Reference

#include <GenericTnPFitter.h>

Inheritance diagram for dqmTnP::VoigtianPlusExponentialFitter:

Public Member Functions
void	fit (TH1 *pass, TH1 *all)
void	setup (double expectedMean_, double massLow, double massHigh, double expectedSigma_, double width_)
	VoigtianPlusExponentialFitter (bool verbose=false)
	~VoigtianPlusExponentialFitter ()
Public Member Functions inherited from dqmTnP::AbstractFitter
	AbstractFitter (bool verbose_=false)
TString	calculateEfficiency (TH3 *pass, TH3 *all, int massDimension, TProfile2D *&eff, TProfile2D *&effChi2, TString plotName="")
TString	calculateEfficiency (TH2 *pass, TH2 *all, int massDimension, TProfile *&eff, TProfile *&effChi2, 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 ()

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 301 of file GenericTnPFitter.h.

Constructor & Destructor Documentation

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

inline

Definition at line 312 of file GenericTnPFitter.h.

References dqmTnP::AbstractFitter::simPdf, and create_public_lumi_plots::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");
  };

dqmTnP::VoigtianPlusExponentialFitter::~VoigtianPlusExponentialFitter ( )

inline

Definition at line 327 of file GenericTnPFitter.h.

{};

Member Function Documentation

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

inlinevirtual

Implements dqmTnP::AbstractFitter.

Definition at line 335 of file GenericTnPFitter.h.

References cond::ecalcond::all, dqmTnP::AbstractFitter::category, dqmTnP::AbstractFitter::chi2, dqmTnP::AbstractFitter::data, dqmTnP::AbstractFitter::efficiency, dqmTnP::AbstractFitter::expectedMean, dqmTnP::AbstractFitter::expectedSigma, cmsPerfPublish::fail(), dqmTnP::AbstractFitter::mass, dqmTnP::AbstractFitter::mean, dqmTnP::AbstractFitter::nBackgroundFail, dqmTnP::AbstractFitter::nBackgroundPass, dqmTnP::AbstractFitter::nSignalAll, dqmTnP::AbstractFitter::sigma, dqmTnP::AbstractFitter::simPdf, and EcalTangentSkim_cfg::Verbose.

                               {
    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 328 of file GenericTnPFitter.h.

References dqmTnP::AbstractFitter::expectedMean, dqmTnP::AbstractFitter::expectedSigma, dqmTnP::AbstractFitter::mass, dqmTnP::AbstractFitter::mean, and create_public_lumi_plots::width.

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 306 of file GenericTnPFitter.h.

RooExponential dqmTnP::VoigtianPlusExponentialFitter::exponentialPass

protected

Definition at line 308 of file GenericTnPFitter.h.

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfFail

protected

Definition at line 309 of file GenericTnPFitter.h.

RooAddPdf dqmTnP::VoigtianPlusExponentialFitter::pdfPass

protected

Definition at line 310 of file GenericTnPFitter.h.

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopeFail

protected

Definition at line 305 of file GenericTnPFitter.h.

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::slopePass

protected

Definition at line 307 of file GenericTnPFitter.h.

RooVoigtian dqmTnP::VoigtianPlusExponentialFitter::voigtian

protected

Definition at line 304 of file GenericTnPFitter.h.

RooRealVar dqmTnP::VoigtianPlusExponentialFitter::width

protected

Definition at line 303 of file GenericTnPFitter.h.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), Vispa.Main.MainWindow.MainWindow::_saveIni(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.rootplot.root2matplotlib.Hist::bar(), python.rootplot.root2matplotlib.Hist::barh(), Vispa.Gui.VispaWidget.VispaWidget::boundingRect(), Vispa.Gui.VispaWidget.VispaWidget::contentRect(), Vispa.Gui.VispaWidget.VispaWidget::defineArrowBackgroundShape(), Vispa.Gui.VispaWidget.VispaWidget::defineCircleBackgroundShape(), Vispa.Gui.VispaWidget.VispaWidget::defineRoundRectBackgroundShape(), Vispa.Gui.VispaWidget.VispaWidget::drawHeaderBackground(), Vispa.Gui.ConnectableWidget.ConnectableWidget::drawPortNames(), Vispa.Gui.VispaWidget.VispaWidget::imageRectF(), Vispa.Gui.ConnectableWidget.ConnectableWidget::positionizeMenuWidget(), Vispa.Views.PropertyView.PropertyView::resizeEvent(), Vispa.Views.PropertyView.PropertyView::sectionResized(), Vispa.Gui.VispaWidget.VispaWidget::setZoom(), and Vispa.Main.MainWindow.MainWindow::updateStartupScreenGeometry().