FlavourHistograms< T > Class Template Reference

#include <FlavourHistorgrams.h>

Inheritance diagram for FlavourHistograms< T >:

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=0)
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 (TH1F *num, TH1F *den, int bin)
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 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 151 of file FlavourHistorgrams.h.

References gather_cfg::cout, DQMStore::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 = 0  ;
    
  // 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 = 0;
  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 = 0;
      theHisto_u = 0;
      theHisto_s = 0;
      theHisto_g = 0;
      theHisto_dus = 0;
    }
    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 = 0;
    theHisto_u = 0;
    theHisto_s = 0;
    theHisto_c = 0;
    theHisto_b = 0;
    theHisto_g = 0;
    theHisto_ni = 0;
    theHisto_dus = 0;
    theHisto_dusg = 0;
    theHisto_pu = 0;
  }
  // defaults for other data members
  theCanvas = 0 ;
}

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 214 of file FlavourHistorgrams.h.

References HistoProviderDQM::book1D(), gather_cfg::cout, MonitorElement::getTH1F(), 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 = 0  ;
    
  // 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 = 0;
  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 = 0;
      theHisto_u = 0;
      theHisto_s = 0;
      theHisto_g = 0;
      theHisto_dus = 0;
    }
    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 = 0;
    theHisto_u = 0;
    theHisto_s = 0;
    theHisto_c = 0;
    theHisto_b = 0;
    theHisto_g = 0;
    theHisto_ni = 0;
    theHisto_dus = 0;
    theHisto_dusg = 0;
    theHisto_pu = 0;
  }

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

template<class T >

FlavourHistograms< T >::~FlavourHistograms ( )

virtual

Definition at line 299 of file FlavourHistorgrams.h.

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

template<class T>

FlavourHistograms< T >::FlavourHistograms ( )

inlineprivate

Definition at line 144 of file FlavourHistorgrams.h.

{}

Member Function Documentation

template<class T>

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

inline

Definition at line 77 of file FlavourHistorgrams.h.

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

{ return theArrayDimension; }

template<class T>

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

inline

Definition at line 69 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theBaseNameDescription ; }

template<class T>

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

inline

Definition at line 68 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theBaseNameTitle       ; }

template<class T >

double FlavourHistograms< T >::ClopperPearsonUnc ( TH1F *  num, TH1F *  den, int  bin  )

protected

Definition at line 512 of file FlavourHistorgrams.h.

                                                                           {
  double effVal = num->GetBinContent(bin)/den->GetBinContent(bin);
  double errLo = TEfficiency::ClopperPearson((int)den->GetBinContent(bin), 
                         (int)num->GetBinContent(bin),
                         0.683,false);
  double errUp = TEfficiency::ClopperPearson((int)den->GetBinContent(bin), 
                         (int)num->GetBinContent(bin),
                         0.683,true);
  return (effVal - errLo > errUp - effVal) ? effVal - errLo : errUp - effVal; 
}

template<class T >

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

protected

Definition at line 524 of file FlavourHistorgrams.h.

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

template<class T>

void FlavourHistograms< T >::divide

(

const FlavourHistograms< T > &

bHD

)

Definition at line 536 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__(), and cuy.FindIssue::__init__().

                                                                     {
  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 503 of file FlavourHistorgrams.h.

References bigModule::plot().

Referenced by SoftLeptonTagPlotter::epsPlot(), TrackCountingTagPlotter::epsPlot(), TrackProbabilityTagPlotter::epsPlot(), and JetTagPlotter::epsPlot().

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

Referenced by SoftLeptonTagPlotter::analyzeTag(), TrackCountingTagPlotter::analyzeTag(), JetTagPlotter::analyzeTag(), and 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 313 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 320 of file FlavourHistorgrams.h.

References i.

{
  if ( theArrayDimension == 0 ) {       
    // 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 577 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 633 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 82 of file FlavourHistorgrams.h.

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

{ return theHisto_all->getTH1F()  ; }    

template<class T>

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

inline

Definition at line 87 of file FlavourHistorgrams.h.

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

{ return theHisto_b->getTH1F()    ; }

template<class T>

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

inline

Definition at line 86 of file FlavourHistorgrams.h.

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

{ return theHisto_c->getTH1F()    ; }

template<class T>

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

inline

Definition at line 83 of file FlavourHistorgrams.h.

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

{ return theHisto_d ->getTH1F()   ; }    

template<class T>

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

inline

Definition at line 90 of file FlavourHistorgrams.h.

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

{ return theHisto_dus->getTH1F()  ; }

template<class T>

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

inline

Definition at line 91 of file FlavourHistorgrams.h.

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

{ return theHisto_dusg->getTH1F() ; }

template<class T>

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

inline

Definition at line 88 of file FlavourHistorgrams.h.

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

{ return theHisto_g->getTH1F()    ; }

template<class T>

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

inline

Definition at line 89 of file FlavourHistorgrams.h.

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

{ return theHisto_ni->getTH1F()   ; }

template<class T>

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

inline

Definition at line 92 of file FlavourHistorgrams.h.

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

{ return theHisto_pu->getTH1F() ; }

template<class T>

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

inline

Definition at line 85 of file FlavourHistorgrams.h.

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

{ return theHisto_s->getTH1F()    ; }

template<class T>

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

inline

Definition at line 84 of file FlavourHistorgrams.h.

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

{ return theHisto_u->getTH1F()    ; }

template<class T>

int FlavourHistograms< T >::indexToPlot ( ) const

inline

Definition at line 79 of file FlavourHistorgrams.h.

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

{return theIndexToPlot; }

template<class T>

double FlavourHistograms< T >::lowerBound ( ) const

inline

Definition at line 71 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theLowerBound          ; } 

template<class T>

int FlavourHistograms< T >::maxDimension ( ) const

inline

Definition at line 78 of file FlavourHistorgrams.h.

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

{return theMaxDimension; }

template<class T>

int FlavourHistograms< T >::nBins ( ) const

inline

Definition at line 70 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theNBins               ; }

template<class T >

void FlavourHistograms< T >::plot

(

TPad *

theCanvas = 0

)

Definition at line 367 of file FlavourHistorgrams.h.

References cuy::col, MonitorElement::getTH1F(), timingPdfMaker::histo, i, bookConverter::max, and RecoBTag::setTDRStyle().

Referenced by cuy.FindIssue::__init__(), SoftLeptonTagPlotter::psPlot(), TrackCountingTagPlotter::psPlot(), TrackProbabilityTagPlotter::psPlot(), and JetTagPlotter::psPlot().

                                                   {

//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 ( 0 ) ;

  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]= 0 ;

  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] ->getTH1F()->GetXaxis()->SetTitle ( theBaseNameDescription.c_str() ) ;
  histo[0] ->getTH1F()->GetYaxis()->SetTitle ( "Arbitrary Units" ) ;
  histo[0] ->getTH1F()->GetYaxis()->SetTitleOffset(1.25) ;

  for (int i=0; i != 4; ++i) {
    if (histo[i]== 0 ) 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] != 0) && (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] != 0 ) histo[3]->getTH1F()->Draw("Same") ;
  }

}

template<class T>

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

inline

Definition at line 76 of file FlavourHistorgrams.h.

{ return thePlotFirst           ; }

template<class T>

bool FlavourHistograms< T >::plotLog ( ) const

inline

Definition at line 74 of file FlavourHistorgrams.h.

{ return thePlotLog             ; }

template<class T>

bool FlavourHistograms< T >::plotNormalized ( ) const

inline

Definition at line 75 of file FlavourHistorgrams.h.

{ return thePlotNormalized      ; }

template<class T >

void FlavourHistograms< T >::setEfficiencyFlag

(

)

Definition at line 557 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 63 of file FlavourHistorgrams.h.

{ theMax = max;}

template<class T>

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

inline

Definition at line 64 of file FlavourHistorgrams.h.

{ theMin = min;}

template<class T >

void FlavourHistograms< T >::settitle

(

const char *

title

)

Definition at line 347 of file FlavourHistorgrams.h.

Referenced by JetTagPlotter::JetTagPlotter(), and 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 73 of file FlavourHistorgrams.h.

{ return theStatistics          ; }

template<class T>

double FlavourHistograms< T >::upperBound ( ) const

inline

Definition at line 72 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

{ return theUpperBound          ; }

Member Data Documentation

template<class T>

unsigned int FlavourHistograms< T >::mcPlots_

private

Definition at line 146 of file FlavourHistorgrams.h.

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

template<class T>

int* FlavourHistograms< T >::theArrayDimension

protected

Definition at line 110 of file FlavourHistorgrams.h.

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

template<class T>

std::string FlavourHistograms< T >::theBaseNameDescription

protected

Definition at line 115 of file FlavourHistorgrams.h.

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

template<class T>

std::string FlavourHistograms< T >::theBaseNameTitle

protected

Definition at line 114 of file FlavourHistorgrams.h.

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

template<class T>

TCanvas* FlavourHistograms< T >::theCanvas

protected

Definition at line 142 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_all

protected

Definition at line 126 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_b

protected

Definition at line 131 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_c

protected

Definition at line 130 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_d

protected

Definition at line 127 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_dus

protected

Definition at line 134 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_dusg

protected

Definition at line 135 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_g

protected

Definition at line 132 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_ni

protected

Definition at line 133 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_pu

protected

Definition at line 136 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_s

protected

Definition at line 129 of file FlavourHistorgrams.h.

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

template<class T>

MonitorElement* FlavourHistograms< T >::theHisto_u

protected

Definition at line 128 of file FlavourHistorgrams.h.

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

template<class T>

int FlavourHistograms< T >::theIndexToPlot

protected

Definition at line 112 of file FlavourHistorgrams.h.

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

template<class T>

double FlavourHistograms< T >::theLowerBound

protected

Definition at line 117 of file FlavourHistorgrams.h.

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

template<class T>

double FlavourHistograms< T >::theMax

protected

Definition at line 123 of file FlavourHistorgrams.h.

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

template<class T>

int FlavourHistograms< T >::theMaxDimension

protected

Definition at line 111 of file FlavourHistorgrams.h.

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

template<class T>

double FlavourHistograms< T >::theMin

protected

Definition at line 123 of file FlavourHistorgrams.h.

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

template<class T>

int FlavourHistograms< T >::theNBins

protected

Definition at line 116 of file FlavourHistorgrams.h.

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

template<class T>

std::string FlavourHistograms< T >::thePlotFirst

protected

Definition at line 122 of file FlavourHistorgrams.h.

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

template<class T>

bool FlavourHistograms< T >::thePlotLog

protected

Definition at line 120 of file FlavourHistorgrams.h.

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

template<class T>

bool FlavourHistograms< T >::thePlotNormalized

protected

Definition at line 121 of file FlavourHistorgrams.h.

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

template<class T>

bool FlavourHistograms< T >::theStatistics

protected

Definition at line 119 of file FlavourHistorgrams.h.

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

template<class T>

double FlavourHistograms< T >::theUpperBound

protected

Definition at line 118 of file FlavourHistorgrams.h.

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