FlavourHistograms< T > Class Template Reference

#include <FlavourHistorgrams.h>

Inheritance diagram for FlavourHistograms< T >:

Public Types
typedef dqm::legacy::DQMStore	DQMStore
typedef dqm::legacy::MonitorElement	MonitorElement

Public Member Functions
int *	arrayDimension () const
std::string	baseNameDescription () const
std::string	baseNameTitle () const
void	divide (const FlavourHistograms< T > &bHD)
void	epsPlot (const std::string &name)
void	fill (const int &flavour, const T &variable) const
void	fill (const int &flavour, const T &variable, const T &w) const
void	fill (const int &flavour, const T *variable) const
	FlavourHistograms (const std::string &baseNameTitle_, const std::string &baseNameDescription_, const int &nBins_, const double &lowerBound_, const double &upperBound_, const std::string &plotFirst_, const std::string &folder, const unsigned int &mc, DQMStore::IGetter &iget)
	FlavourHistograms (const std::string &baseNameTitle_, const std::string &baseNameDescription_, const int &nBins_, const double &lowerBound_, const double &upperBound_, const bool &statistics_, const bool &plotLog_, const bool &plotNormalized_, const std::string &plotFirst_, const std::string &folder, const unsigned int &mc, DQMStore::IBooker &ibook)
std::vector< TH1F * >	getHistoVector () const
TH1F *	histo_all () const
TH1F *	histo_b () const
TH1F *	histo_c () const
TH1F *	histo_d () const
TH1F *	histo_dus () const
TH1F *	histo_dusg () const
TH1F *	histo_g () const
TH1F *	histo_ni () const
TH1F *	histo_pu () const
TH1F *	histo_s () const
TH1F *	histo_u () const
int	indexToPlot () const
double	lowerBound () const
int	maxDimension () const
int	nBins () const
void	plot (TPad *theCanvas=nullptr)
std::string	plotFirst () const
bool	plotLog () const
bool	plotNormalized () const
void	setEfficiencyFlag ()
void	SetMaximum (const double &max)
void	SetMinimum (const double &min)
void	settitle (const char *title)
bool	statistics () const
double	upperBound () const
virtual	~FlavourHistograms ()

Protected Member Functions
double	ClopperPearsonUnc (double num, double den)
void	ComputeEfficiency (TH1F *num, TH1F *den, int bin)
void	fillVariable (const int &flavour, const T &var, const T &w) const

Protected Attributes
int *	theArrayDimension
std::string	theBaseNameDescription
std::string	theBaseNameTitle
TCanvas *	theCanvas
MonitorElement *	theHisto_all
MonitorElement *	theHisto_b
MonitorElement *	theHisto_c
MonitorElement *	theHisto_d
MonitorElement *	theHisto_dus
MonitorElement *	theHisto_dusg
MonitorElement *	theHisto_g
MonitorElement *	theHisto_ni
MonitorElement *	theHisto_pu
MonitorElement *	theHisto_s
MonitorElement *	theHisto_u
int	theIndexToPlot
double	theLowerBound
double	theMax
int	theMaxDimension
double	theMin
int	theNBins
std::string	thePlotFirst
bool	thePlotLog
bool	thePlotNormalized
bool	theStatistics
double	theUpperBound

Private Member Functions
	FlavourHistograms ()

Private Attributes
unsigned int	mcPlots_

Detailed Description

template<class T>
class FlavourHistograms< T >

Definition at line 26 of file FlavourHistorgrams.h.

Member Typedef Documentation

template<class T>

typedef dqm::legacy::DQMStore FlavourHistograms< T >::DQMStore

Definition at line 28 of file FlavourHistorgrams.h.

template<class T>

typedef dqm::legacy::MonitorElement FlavourHistograms< T >::MonitorElement

Definition at line 29 of file FlavourHistorgrams.h.

Constructor & Destructor Documentation

template<class T >

FlavourHistograms< T >::FlavourHistograms	(	const std::string &	baseNameTitle_,
		const std::string &	baseNameDescription_,
		const int &	nBins_,
		const double &	lowerBound_,
		const double &	upperBound_,
		const std::string &	plotFirst_,
		const std::string &	folder,
		const unsigned int &	mc,
		DQMStore::IGetter &	iget
	)

Definition at line 157 of file FlavourHistorgrams.h.

References gather_cfg::cout, dqm::implementation::IGetter::get(), FlavourHistograms< T >::mcPlots_, FlavourHistograms< T >::theArrayDimension, FlavourHistograms< T >::theBaseNameTitle, FlavourHistograms< T >::theCanvas, FlavourHistograms< T >::theHisto_all, FlavourHistograms< T >::theHisto_b, FlavourHistograms< T >::theHisto_c, FlavourHistograms< T >::theHisto_d, FlavourHistograms< T >::theHisto_dus, FlavourHistograms< T >::theHisto_dusg, FlavourHistograms< T >::theHisto_g, FlavourHistograms< T >::theHisto_ni, FlavourHistograms< T >::theHisto_pu, FlavourHistograms< T >::theHisto_s, FlavourHistograms< T >::theHisto_u, and FlavourHistograms< T >::thePlotFirst.

    : theMaxDimension(-1),
      theIndexToPlot(-1),
      theBaseNameTitle(baseNameTitle_),
      theBaseNameDescription(baseNameDescription_),
      theNBins(nBins_),
      theLowerBound(lowerBound_),
      theUpperBound(upperBound_),
      theStatistics(false),
      thePlotLog(false),
      thePlotNormalized(false),
      thePlotFirst(plotFirst_),
      theMin(-1.),
      theMax(-1.),
      mcPlots_(mc) {
  // defaults for array dimensions
  theArrayDimension = nullptr;

  // check plot order string
  if (thePlotFirst == "l" || thePlotFirst == "c" || thePlotFirst == "b") {
    // OK
  } else {
    // not correct: print warning and set default (l)
    std::cout << "FlavourHistograms::FlavourHistograms : thePlotFirst was not correct : " << thePlotFirst << std::endl;
    std::cout << "FlavourHistograms::FlavourHistograms : Set it to default value (l)! " << std::endl;
    thePlotFirst = "l";
  }

  if (mcPlots_ % 2 == 0)
    theHisto_all = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "ALL");
  else
    theHisto_all = nullptr;
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      theHisto_d = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "D");
      theHisto_u = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "U");
      theHisto_s = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "S");
      theHisto_g = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "G");
      theHisto_dus = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "DUS");
    } else {
      theHisto_d = nullptr;
      theHisto_u = nullptr;
      theHisto_s = nullptr;
      theHisto_g = nullptr;
      theHisto_dus = nullptr;
    }
    theHisto_c = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "C");
    theHisto_b = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "B");
    theHisto_ni = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "NI");
    theHisto_dusg = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "DUSG");
    theHisto_pu = iget.get("Btag/" + folder + "/" + theBaseNameTitle + "PU");
  } else {
    theHisto_d = nullptr;
    theHisto_u = nullptr;
    theHisto_s = nullptr;
    theHisto_c = nullptr;
    theHisto_b = nullptr;
    theHisto_g = nullptr;
    theHisto_ni = nullptr;
    theHisto_dus = nullptr;
    theHisto_dusg = nullptr;
    theHisto_pu = nullptr;
  }
  // defaults for other data members
  theCanvas = nullptr;
}

template<class T >

FlavourHistograms< T >::FlavourHistograms	(	const std::string &	baseNameTitle_,
		const std::string &	baseNameDescription_,
		const int &	nBins_,
		const double &	lowerBound_,
		const double &	upperBound_,
		const bool &	statistics_,
		const bool &	plotLog_,
		const bool &	plotNormalized_,
		const std::string &	plotFirst_,
		const std::string &	folder,
		const unsigned int &	mc,
		DQMStore::IBooker &	ibook
	)

Definition at line 233 of file FlavourHistorgrams.h.

References HistoProviderDQM::book1D(), gather_cfg::cout, dqm::impl::MonitorElement::enableSumw2(), FlavourHistograms< T >::mcPlots_, FlavourHistograms< T >::theArrayDimension, FlavourHistograms< T >::theBaseNameDescription, FlavourHistograms< T >::theBaseNameTitle, FlavourHistograms< T >::theCanvas, FlavourHistograms< T >::theHisto_all, FlavourHistograms< T >::theHisto_b, FlavourHistograms< T >::theHisto_c, FlavourHistograms< T >::theHisto_d, FlavourHistograms< T >::theHisto_dus, FlavourHistograms< T >::theHisto_dusg, FlavourHistograms< T >::theHisto_g, FlavourHistograms< T >::theHisto_ni, FlavourHistograms< T >::theHisto_pu, FlavourHistograms< T >::theHisto_s, FlavourHistograms< T >::theHisto_u, FlavourHistograms< T >::theLowerBound, FlavourHistograms< T >::theNBins, FlavourHistograms< T >::thePlotFirst, FlavourHistograms< T >::theStatistics, and FlavourHistograms< T >::theUpperBound.

    : theMaxDimension(-1),
      theIndexToPlot(-1),
      theBaseNameTitle(baseNameTitle_),
      theBaseNameDescription(baseNameDescription_),
      theNBins(nBins_),
      theLowerBound(lowerBound_),
      theUpperBound(upperBound_),
      theStatistics(statistics_),
      thePlotLog(plotLog_),
      thePlotNormalized(plotNormalized_),
      thePlotFirst(plotFirst_),
      theMin(-1.),
      theMax(-1.),
      mcPlots_(mc) {
  // defaults for array dimensions
  theArrayDimension = nullptr;

  // check plot order string
  if (thePlotFirst == "l" || thePlotFirst == "c" || thePlotFirst == "b") {
    // OK
  } else {
    // not correct: print warning and set default (l)
    std::cout << "FlavourHistograms::FlavourHistograms : thePlotFirst was not correct : " << thePlotFirst << std::endl;
    std::cout << "FlavourHistograms::FlavourHistograms : Set it to default value (l)! " << std::endl;
    thePlotFirst = "l";
  }

  // book histos
  HistoProviderDQM prov("Btag", folder, ibook);
  if (mcPlots_ % 2 == 0)
    theHisto_all = (prov.book1D(
        theBaseNameTitle + "ALL", theBaseNameDescription + " all jets", theNBins, theLowerBound, theUpperBound));
  else
    theHisto_all = nullptr;
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      theHisto_d = (prov.book1D(
          theBaseNameTitle + "D", theBaseNameDescription + " d-jets", theNBins, theLowerBound, theUpperBound));
      theHisto_u = (prov.book1D(
          theBaseNameTitle + "U", theBaseNameDescription + " u-jets", theNBins, theLowerBound, theUpperBound));
      theHisto_s = (prov.book1D(
          theBaseNameTitle + "S", theBaseNameDescription + " s-jets", theNBins, theLowerBound, theUpperBound));
      theHisto_g = (prov.book1D(
          theBaseNameTitle + "G", theBaseNameDescription + " g-jets", theNBins, theLowerBound, theUpperBound));
      theHisto_dus = (prov.book1D(
          theBaseNameTitle + "DUS", theBaseNameDescription + " dus-jets", theNBins, theLowerBound, theUpperBound));
    } else {
      theHisto_d = nullptr;
      theHisto_u = nullptr;
      theHisto_s = nullptr;
      theHisto_g = nullptr;
      theHisto_dus = nullptr;
    }
    theHisto_c = (prov.book1D(
        theBaseNameTitle + "C", theBaseNameDescription + " c-jets", theNBins, theLowerBound, theUpperBound));
    theHisto_b = (prov.book1D(
        theBaseNameTitle + "B", theBaseNameDescription + " b-jets", theNBins, theLowerBound, theUpperBound));
    theHisto_ni = (prov.book1D(
        theBaseNameTitle + "NI", theBaseNameDescription + " ni-jets", theNBins, theLowerBound, theUpperBound));
    theHisto_dusg = (prov.book1D(
        theBaseNameTitle + "DUSG", theBaseNameDescription + " dusg-jets", theNBins, theLowerBound, theUpperBound));
    theHisto_pu = (prov.book1D(
        theBaseNameTitle + "PU", theBaseNameDescription + " pu-jets", theNBins, theLowerBound, theUpperBound));
  } else {
    theHisto_d = nullptr;
    theHisto_u = nullptr;
    theHisto_s = nullptr;
    theHisto_c = nullptr;
    theHisto_b = nullptr;
    theHisto_g = nullptr;
    theHisto_ni = nullptr;
    theHisto_dus = nullptr;
    theHisto_dusg = nullptr;
    theHisto_pu = nullptr;
  }

  // statistics if requested
  if (theStatistics) {
    if (theHisto_all)
      theHisto_all->enableSumw2();
    if (mcPlots_) {
      if (mcPlots_ > 2) {
        theHisto_d->enableSumw2();
        theHisto_u->enableSumw2();
        theHisto_s->enableSumw2();
        theHisto_g->enableSumw2();
        theHisto_dus->enableSumw2();
      }
      theHisto_c->enableSumw2();
      theHisto_b->enableSumw2();
      theHisto_ni->enableSumw2();
      theHisto_dusg->enableSumw2();
      theHisto_pu->enableSumw2();
    }
  }
  // defaults for other data members
  theCanvas = nullptr;
}

template<class T >

FlavourHistograms< T >::~FlavourHistograms ( )

virtual

Definition at line 345 of file FlavourHistorgrams.h.

                                         {
  // delete the canvas*/
  delete theCanvas;
}

template<class T>

FlavourHistograms< T >::FlavourHistograms ( )

inlineprivate

Definition at line 151 of file FlavourHistorgrams.h.

{}

Member Function Documentation

template<class T>

int* FlavourHistograms< T >::arrayDimension ( ) const

inline

Definition at line 86 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theArrayDimension.

Referenced by TrackIPHistograms< T >::fill().

{ return theArrayDimension; }

template<class T>

std::string FlavourHistograms< T >::baseNameDescription ( ) const

inline

Definition at line 78 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theBaseNameDescription.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theBaseNameDescription; }

template<class T>

std::string FlavourHistograms< T >::baseNameTitle ( ) const

inline

Definition at line 77 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theBaseNameTitle.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theBaseNameTitle; }

template<class T >

double FlavourHistograms< T >::ClopperPearsonUnc ( double  num, double  den  )

protected

Definition at line 570 of file FlavourHistorgrams.h.

References SiStripPI::max.

                                                                     {
  double effVal = num / den;
  double errLo = TEfficiency::ClopperPearson(static_cast<int>(den), static_cast<int>(num), 0.683, false);
  double errUp = TEfficiency::ClopperPearson(static_cast<int>(den), static_cast<int>(num), 0.683, true);
  return std::max(effVal - errLo, errUp - effVal);
}

template<class T >

void FlavourHistograms< T >::ComputeEfficiency ( TH1F *  num, TH1F *  den, int  bin  )

protected

Definition at line 578 of file FlavourHistorgrams.h.

                                                                          {
  double effVal = 1.;
  double errVal = 0.;
  double numVal = num->GetBinContent(bin);
  double denVal = den->GetBinContent(bin);
  if (denVal > 0 && numVal <= denVal) {
    effVal = numVal / denVal;
    errVal = ClopperPearsonUnc(numVal, denVal);
  }
  num->SetBinContent(bin, effVal);
  num->SetBinError(bin, errVal);
}

template<class T >

void FlavourHistograms< T >::divide

(

const FlavourHistograms< T > &

bHD

)

Definition at line 592 of file FlavourHistorgrams.h.

References newFWLiteAna::bin, FlavourHistograms< T >::histo_all(), FlavourHistograms< T >::histo_b(), FlavourHistograms< T >::histo_c(), FlavourHistograms< T >::histo_d(), FlavourHistograms< T >::histo_dus(), FlavourHistograms< T >::histo_dusg(), FlavourHistograms< T >::histo_g(), FlavourHistograms< T >::histo_ni(), FlavourHistograms< T >::histo_pu(), FlavourHistograms< T >::histo_s(), and FlavourHistograms< T >::histo_u().

Referenced by python.rootplot.utilities.Hist::__div__().

                                                                 {
  for (int bin = 0; bin < theNBins + 2; bin++) {
    if (theHisto_all)
      ComputeEfficiency(theHisto_all->getTH1F(), bHD.histo_all(), bin);
    if (mcPlots_) {
      if (mcPlots_ > 2) {
        ComputeEfficiency(theHisto_d->getTH1F(), bHD.histo_d(), bin);
        ComputeEfficiency(theHisto_u->getTH1F(), bHD.histo_u(), bin);
        ComputeEfficiency(theHisto_s->getTH1F(), bHD.histo_s(), bin);
        ComputeEfficiency(theHisto_g->getTH1F(), bHD.histo_g(), bin);
        ComputeEfficiency(theHisto_dus->getTH1F(), bHD.histo_dus(), bin);
      }
      ComputeEfficiency(theHisto_c->getTH1F(), bHD.histo_c(), bin);
      ComputeEfficiency(theHisto_b->getTH1F(), bHD.histo_b(), bin);
      ComputeEfficiency(theHisto_ni->getTH1F(), bHD.histo_ni(), bin);
      ComputeEfficiency(theHisto_dusg->getTH1F(), bHD.histo_dusg(), bin);
      ComputeEfficiency(theHisto_pu->getTH1F(), bHD.histo_pu(), bin);
    }
  }
}

template<class T >

void FlavourHistograms< T >::epsPlot

(

const std::string &

name

)

Definition at line 562 of file FlavourHistorgrams.h.

References bigModule::plot().

                                                        {
  TCanvas tc(theBaseNameTitle.c_str(), theBaseNameDescription.c_str());

  plot(&tc);
  tc.Print((name + theBaseNameTitle + ".eps").c_str());
}

template<class T >

void FlavourHistograms< T >::fill	(	const int &	flavour,
		const T &	variable
	)		const

Definition at line 352 of file FlavourHistorgrams.h.

Referenced by TrackIPHistograms< T >::fill().

                                                                           {
  // For single variables and arrays (for arrays only a single index can be filled)
  fillVariable(flavour, variable, 1.);
}

template<class T >

void FlavourHistograms< T >::fill	(	const int &	flavour,
		const T &	variable,
		const T &	w
	)		const

Definition at line 358 of file FlavourHistorgrams.h.

                                                                                       {
  // For single variables and arrays (for arrays only a single index can be filled)
  fillVariable(flavour, variable, w);
}

template<class T >

void FlavourHistograms< T >::fill	(	const int &	flavour,
		const T *	variable
	)		const

Definition at line 364 of file FlavourHistorgrams.h.

References mps_fire::i.

                                                                           {
  if (theArrayDimension == nullptr) {
    // single variable
    fillVariable(flavour, *variable, 1.);
  } else {
    // array
    int iMax = (*theArrayDimension > theMaxDimension) ? theMaxDimension : *theArrayDimension;
    //
    for (int i = 0; i != iMax; ++i) {
      // check if only one index to be plotted (<0: switched off -> plot all)
      if ((theIndexToPlot < 0) || (i == theIndexToPlot)) {
        fillVariable(flavour, *(variable + i), 1.);
      }
    }

    // if single index to be filled but not enough entries: fill 0.0 (convention!)
    if (theIndexToPlot >= iMax) {
      // cout << "==>> The index to be filled is too big -> fill 0.0 : " << theBaseNameTitle << " : " << theIndexToPlot << " >= " << iMax << endl ;
      const T& theZero = static_cast<T>(0.0);
      fillVariable(flavour, theZero, 1.);
    }
  }
}

template<class T >

void FlavourHistograms< T >::fillVariable ( const int &  flavour, const T &  var, const T &  w  ) const

protected

Definition at line 634 of file FlavourHistorgrams.h.

                                                                                          {
  // all, except for the Jet Multiplicity which is not filled for each jets but for each events
  if ((theBaseNameDescription != "Jet Multiplicity" || flavour == -1) && theHisto_all)
    theHisto_all->Fill(var, w);

  // flavour specific
  if (!mcPlots_ || (theBaseNameDescription == "Jet Multiplicity" && flavour == -1))
    return;

  switch (flavour) {
    case 1:
      if (mcPlots_ > 2) {
        theHisto_d->Fill(var, w);
        if (theBaseNameDescription != "Jet Multiplicity")
          theHisto_dus->Fill(var, w);
      }
      if (theBaseNameDescription != "Jet Multiplicity")
        theHisto_dusg->Fill(var, w);
      return;
    case 2:
      if (mcPlots_ > 2) {
        theHisto_u->Fill(var, w);
        if (theBaseNameDescription != "Jet Multiplicity")
          theHisto_dus->Fill(var, w);
      }
      if (theBaseNameDescription != "Jet Multiplicity")
        theHisto_dusg->Fill(var, w);
      return;
    case 3:
      if (mcPlots_ > 2) {
        theHisto_s->Fill(var, w);
        if (theBaseNameDescription != "Jet Multiplicity")
          theHisto_dus->Fill(var, w);
      }
      if (theBaseNameDescription != "Jet Multiplicity")
        theHisto_dusg->Fill(var, w);
      return;
    case 4:
      theHisto_c->Fill(var, w);
      return;
    case 5:
      theHisto_b->Fill(var, w);
      return;
    case 21:
      if (mcPlots_ > 2)
        theHisto_g->Fill(var, w);
      if (theBaseNameDescription != "Jet Multiplicity")
        theHisto_dusg->Fill(var, w);
      return;
    case 123:
      if (mcPlots_ > 2 && theBaseNameDescription == "Jet Multiplicity")
        theHisto_dus->Fill(var, w);
      return;
    case 12321:
      if (theBaseNameDescription == "Jet Multiplicity")
        theHisto_dusg->Fill(var, w);
      return;
    case 20:
      theHisto_pu->Fill(var, w);
      return;
    default:
      theHisto_ni->Fill(var, w);
      return;
  }
}

template<class T >

std::vector< TH1F * > FlavourHistograms< T >::getHistoVector

(

)

const

Definition at line 701 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

                                                            {
  std::vector<TH1F*> histoVector;
  if (theHisto_all)
    histoVector.push_back(theHisto_all->getTH1F());
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      histoVector.push_back(theHisto_d->getTH1F());
      histoVector.push_back(theHisto_u->getTH1F());
      histoVector.push_back(theHisto_s->getTH1F());
      histoVector.push_back(theHisto_g->getTH1F());
      histoVector.push_back(theHisto_dus->getTH1F());
    }
    histoVector.push_back(theHisto_c->getTH1F());
    histoVector.push_back(theHisto_b->getTH1F());
    histoVector.push_back(theHisto_ni->getTH1F());
    histoVector.push_back(theHisto_dusg->getTH1F());
    histoVector.push_back(theHisto_pu->getTH1F());
  }
  return histoVector;
}

template<class T>

TH1F* FlavourHistograms< T >::histo_all ( ) const

inline

Definition at line 91 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_all.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_all->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_b ( ) const

inline

Definition at line 96 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_b.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_b->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_c ( ) const

inline

Definition at line 95 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_c.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_c->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_d ( ) const

inline

Definition at line 92 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_d.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_d->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_dus ( ) const

inline

Definition at line 99 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_dus.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_dus->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_dusg ( ) const

inline

Definition at line 100 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_dusg.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_dusg->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_g ( ) const

inline

Definition at line 97 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_g.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_g->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_ni ( ) const

inline

Definition at line 98 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_ni.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_ni->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_pu ( ) const

inline

Definition at line 101 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_pu.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_pu->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_s ( ) const

inline

Definition at line 94 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_s.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_s->getTH1F(); }

template<class T>

TH1F* FlavourHistograms< T >::histo_u ( ) const

inline

Definition at line 93 of file FlavourHistorgrams.h.

References dqm::legacy::MonitorElement::getTH1F(), and FlavourHistograms< T >::theHisto_u.

Referenced by FlavourHistograms< T >::divide().

{ return theHisto_u->getTH1F(); }

template<class T>

int FlavourHistograms< T >::indexToPlot ( ) const

inline

Definition at line 88 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theIndexToPlot.

Referenced by TrackIPHistograms< T >::fill().

{ return theIndexToPlot; }

template<class T>

double FlavourHistograms< T >::lowerBound ( ) const

inline

Definition at line 80 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theLowerBound.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theLowerBound; }

template<class T>

int FlavourHistograms< T >::maxDimension ( ) const

inline

Definition at line 87 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theMaxDimension.

Referenced by TrackIPHistograms< T >::fill().

{ return theMaxDimension; }

template<class T>

int FlavourHistograms< T >::nBins ( ) const

inline

Definition at line 79 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theNBins.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theNBins; }

template<class T >

void FlavourHistograms< T >::plot

(

TPad *

theCanvas = nullptr

)

Definition at line 409 of file FlavourHistorgrams.h.

References cuy::col, dqm::legacy::MonitorElement::getTH1F(), timingPdfMaker::histo, mps_fire::i, SiStripPI::max, dqm::impl::MonitorElement::setAxisTitle(), and RecoBTag::setTDRStyle().

                                                 {
  //fixme:
  bool btppNI = false;
  bool btppColour = true;

  if (theCanvas)
    theCanvas->cd();

  RecoBTag::setTDRStyle()->cd();
  gPad->UseCurrentStyle();
  //   if ( !btppTitle ) gStyle->SetOptTitle ( 0 ) ;
  //
  //   // here: plot histograms in a canvas
  //   theCanvas = new TCanvas ( "C" + theBaseNameTitle , "C" + theBaseNameDescription ,
  //                 btppXCanvas , btppYCanvas ) ;
  //   theCanvas->SetFillColor ( 0 ) ;
  //   theCanvas->cd  ( 1 ) ;
  gPad->SetLogy(0);
  if (thePlotLog)
    gPad->SetLogy(1);
  gPad->SetGridx(0);
  gPad->SetGridy(0);
  gPad->SetTitle(nullptr);

  MonitorElement* histo[4];
  int col[4], lineStyle[4], markerStyle[4];
  int lineWidth = 1;

  const double markerSize = gPad->GetWh() * gPad->GetHNDC() / 500.;

  // default (l)
  histo[0] = theHisto_dusg;
  //CW histo_1 = theHisto_dus ;
  histo[1] = theHisto_b;
  histo[2] = theHisto_c;
  histo[3] = nullptr;

  double max = theMax;
  if (theMax <= 0.) {
    max = theHisto_dusg->getTH1F()->GetMaximum();
    if (theHisto_b->getTH1F()->GetMaximum() > max)
      max = theHisto_b->getTH1F()->GetMaximum();
    if (theHisto_c->getTH1F()->GetMaximum() > max)
      max = theHisto_c->getTH1F()->GetMaximum();
  }

  if (btppNI) {
    histo[3] = theHisto_ni;
    if (theHisto_ni->getTH1F()->GetMaximum() > max)
      max = theHisto_ni->getTH1F()->GetMaximum();
  }

  if (btppColour) {  // print colours
    col[0] = 4;
    col[1] = 2;
    col[2] = 6;
    col[3] = 3;
    lineStyle[0] = 1;
    lineStyle[1] = 1;
    lineStyle[2] = 1;
    lineStyle[3] = 1;
    markerStyle[0] = 20;
    markerStyle[1] = 21;
    markerStyle[2] = 22;
    markerStyle[3] = 23;
    lineWidth = 1;
  } else {  // different marker/line styles
    col[1] = 1;
    col[2] = 1;
    col[3] = 1;
    col[0] = 1;
    lineStyle[0] = 2;
    lineStyle[1] = 1;
    lineStyle[2] = 3;
    lineStyle[3] = 4;
    markerStyle[0] = 20;
    markerStyle[1] = 21;
    markerStyle[2] = 22;
    markerStyle[3] = 23;
  }

  // if changing order (NI stays always last)

  // c to plot first
  if (thePlotFirst == "c") {
    histo[0] = theHisto_c;
    if (btppColour)
      col[0] = 6;
    if (!btppColour)
      lineStyle[0] = 3;
    histo[2] = theHisto_dusg;
    if (btppColour)
      col[2] = 4;
    if (!btppColour)
      lineStyle[2] = 2;
  }

  // b to plot first
  if (thePlotFirst == "b") {
    histo[0] = theHisto_b;
    if (btppColour)
      col[0] = 2;
    if (!btppColour)
      lineStyle[0] = 1;
    histo[1] = theHisto_dusg;
    if (btppColour)
      col[1] = 4;
    if (!btppColour)
      lineStyle[1] = 2;
  }

  histo[0]->setAxisTitle(theBaseNameDescription);
  histo[0]->getTH1F()->GetYaxis()->SetTitle("Arbitrary Units");
  histo[0]->getTH1F()->GetYaxis()->SetTitleOffset(1.25);

  for (int i = 0; i != 4; ++i) {
    if (histo[i] == nullptr)
      continue;
    histo[i]->getTH1F()->SetStats(false);
    histo[i]->getTH1F()->SetLineStyle(lineStyle[i]);
    histo[i]->getTH1F()->SetLineWidth(lineWidth);
    histo[i]->getTH1F()->SetLineColor(col[i]);
    histo[i]->getTH1F()->SetMarkerStyle(markerStyle[i]);
    histo[i]->getTH1F()->SetMarkerColor(col[i]);
    histo[i]->getTH1F()->SetMarkerSize(markerSize);
  }

  if (thePlotNormalized) {
    if (histo[0]->getTH1F()->GetEntries() != 0) {
      histo[0]->getTH1F()->DrawNormalized();
    } else {
      histo[0]->getTH1F()->SetMaximum(1.0);
      histo[0]->getTH1F()->Draw();
    }
    if (histo[1]->getTH1F()->GetEntries() != 0)
      histo[1]->getTH1F()->DrawNormalized("Same");
    if (histo[2]->getTH1F()->GetEntries() != 0)
      histo[2]->getTH1F()->DrawNormalized("Same");
    if ((histo[3] != nullptr) && (histo[3]->getTH1F()->GetEntries() != 0))
      histo[3]->getTH1F()->DrawNormalized("Same");
  } else {
    histo[0]->getTH1F()->SetMaximum(max * 1.05);
    if (theMin != -1.)
      histo[0]->getTH1F()->SetMinimum(theMin);
    histo[0]->getTH1F()->Draw();
    histo[1]->getTH1F()->Draw("Same");
    histo[2]->getTH1F()->Draw("Same");
    if (histo[3] != nullptr)
      histo[3]->getTH1F()->Draw("Same");
  }
}

template<class T>

std::string FlavourHistograms< T >::plotFirst ( ) const

inline

Definition at line 85 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::thePlotFirst.

{ return thePlotFirst; }

template<class T>

bool FlavourHistograms< T >::plotLog ( ) const

inline

Definition at line 83 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::thePlotLog.

{ return thePlotLog; }

template<class T>

bool FlavourHistograms< T >::plotNormalized ( ) const

inline

Definition at line 84 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::thePlotNormalized.

{ return thePlotNormalized; }

template<class T >

void FlavourHistograms< T >::setEfficiencyFlag

(

)

Definition at line 614 of file FlavourHistorgrams.h.

                                             {
  if (theHisto_all)
    theHisto_all->setEfficiencyFlag();
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      theHisto_d->setEfficiencyFlag();
      theHisto_u->setEfficiencyFlag();
      theHisto_s->setEfficiencyFlag();
      theHisto_g->setEfficiencyFlag();
      theHisto_dus->setEfficiencyFlag();
    }
    theHisto_c->setEfficiencyFlag();
    theHisto_b->setEfficiencyFlag();
    theHisto_ni->setEfficiencyFlag();
    theHisto_dusg->setEfficiencyFlag();
    theHisto_pu->setEfficiencyFlag();
  }
}

template<class T>

void FlavourHistograms< T >::SetMaximum ( const double &  max )

inline

Definition at line 73 of file FlavourHistorgrams.h.

References SiStripPI::max, and FlavourHistograms< T >::theMax.

{ theMax = max; }

template<class T>

void FlavourHistograms< T >::SetMinimum ( const double &  min )

inline

Definition at line 74 of file FlavourHistorgrams.h.

References SiStripPI::min, and FlavourHistograms< T >::theMin.

{ theMin = min; }

template<class T >

void FlavourHistograms< T >::settitle

(

const char *

title

)

Definition at line 389 of file FlavourHistorgrams.h.

Referenced by TrackIPHistograms< T >::settitle().

                                                     {
  if (theHisto_all)
    theHisto_all->setAxisTitle(title);
  if (mcPlots_) {
    if (mcPlots_ > 2) {
      theHisto_d->setAxisTitle(title);
      theHisto_u->setAxisTitle(title);
      theHisto_s->setAxisTitle(title);
      theHisto_g->setAxisTitle(title);
      theHisto_dus->setAxisTitle(title);
    }
    theHisto_c->setAxisTitle(title);
    theHisto_b->setAxisTitle(title);
    theHisto_ni->setAxisTitle(title);
    theHisto_dusg->setAxisTitle(title);
    theHisto_pu->setAxisTitle(title);
  }
}

template<class T>

bool FlavourHistograms< T >::statistics ( ) const

inline

Definition at line 82 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theStatistics.

{ return theStatistics; }

template<class T>

double FlavourHistograms< T >::upperBound ( ) const

inline

Definition at line 81 of file FlavourHistorgrams.h.

References FlavourHistograms< T >::theUpperBound.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theUpperBound; }

Member Data Documentation

template<class T>

unsigned int FlavourHistograms< T >::mcPlots_

private

Definition at line 153 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms().

template<class T>

int* FlavourHistograms< T >::theArrayDimension

protected

Definition at line 117 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::arrayDimension(), and FlavourHistograms< T >::FlavourHistograms().

template<class T>

std::string FlavourHistograms< T >::theBaseNameDescription

protected

Definition at line 122 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::baseNameDescription(), and FlavourHistograms< T >::FlavourHistograms().

template<class T>

std::string FlavourHistograms< T >::theBaseNameTitle

protected

Definition at line 121 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::baseNameTitle(), and FlavourHistograms< T >::FlavourHistograms().

template<class T>

TCanvas* FlavourHistograms< T >::theCanvas

protected

Definition at line 148 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_all

protected

Definition at line 133 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_all().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_b

protected

Definition at line 138 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_b().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_c

protected

Definition at line 137 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_c().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_d

protected

Definition at line 134 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_d().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_dus

protected

Definition at line 141 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_dus().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_dusg

protected

Definition at line 142 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_dusg().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_g

protected

Definition at line 139 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_g().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_ni

protected

Definition at line 140 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_ni().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_pu

protected

Definition at line 143 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_pu().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_s

protected

Definition at line 136 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_s().

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_u

protected

Definition at line 135 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::histo_u().

template<class T>

int FlavourHistograms< T >::theIndexToPlot

protected

Definition at line 119 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::indexToPlot().

template<class T>

double FlavourHistograms< T >::theLowerBound

protected

Definition at line 124 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::lowerBound().

template<class T>

double FlavourHistograms< T >::theMax

protected

Definition at line 130 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::SetMaximum().

template<class T>

int FlavourHistograms< T >::theMaxDimension

protected

Definition at line 118 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::maxDimension().

template<class T>

double FlavourHistograms< T >::theMin

protected

Definition at line 130 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::SetMinimum().

template<class T>

int FlavourHistograms< T >::theNBins

protected

Definition at line 123 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::nBins().

template<class T>

std::string FlavourHistograms< T >::thePlotFirst

protected

Definition at line 129 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::plotFirst().

template<class T>

bool FlavourHistograms< T >::thePlotLog

protected

Definition at line 127 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::plotLog().

template<class T>

bool FlavourHistograms< T >::thePlotNormalized

protected

Definition at line 128 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::plotNormalized().

template<class T>

bool FlavourHistograms< T >::theStatistics

protected

Definition at line 126 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::statistics().

template<class T>

double FlavourHistograms< T >::theUpperBound

protected

Definition at line 125 of file FlavourHistorgrams.h.

Referenced by FlavourHistograms< T >::FlavourHistograms(), and FlavourHistograms< T >::upperBound().