CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes
FlavourHistograms< T > Class Template Reference

#include <FlavourHistorgrams.h>

Inheritance diagram for FlavourHistograms< T >:
TrackIPHistograms< T >

Public Member Functions

int * arrayDimension () const
 
std::string baseNameDescription () const
 
std::string baseNameTitle () const
 
void divide (const FlavourHistograms< T > &bHD) const
 
void epsPlot (const std::string &name)
 
void fill (const int &flavour, const T &variable) 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 bool &statistics_, const bool &plotLog_, const bool &plotNormalized_, const std::string &plotFirst_, const bool &update, const std::string &folder, const bool &mc)
 
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_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 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

void fillVariable (const int &flavour, const T &var) const
 

Protected Attributes

int * theArrayDimension
 
std::string theBaseNameDescription
 
std::string theBaseNameTitle
 
TCanvas * theCanvas
 
MonitorElementtheHisto_all
 
MonitorElementtheHisto_b
 
MonitorElementtheHisto_c
 
MonitorElementtheHisto_d
 
MonitorElementtheHisto_dus
 
MonitorElementtheHisto_dusg
 
MonitorElementtheHisto_g
 
MonitorElementtheHisto_ni
 
MonitorElementtheHisto_s
 
MonitorElementtheHisto_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

bool mcPlots_
 

Detailed Description

template<class T>
class FlavourHistograms< T >

Definition at line 28 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 bool &  statistics_,
const bool &  plotLog_,
const bool &  plotNormalized_,
const std::string &  plotFirst_,
const bool &  update,
const std::string &  folder,
const bool &  mc 
)

Definition at line 145 of file FlavourHistorgrams.h.

References HistoProviderDQM::access(), 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_s, FlavourHistograms< T >::theHisto_u, FlavourHistograms< T >::theLowerBound, FlavourHistograms< T >::theNBins, FlavourHistograms< T >::thePlotFirst, FlavourHistograms< T >::theStatistics, and FlavourHistograms< T >::theUpperBound.

148  :
149  // BaseFlavourHistograms () ,
150  // theVariable ( variable_ ) ,
151  theMaxDimension(-1), theIndexToPlot(-1), theBaseNameTitle ( baseNameTitle_ ) , theBaseNameDescription ( baseNameDescription_ ) ,
152  theNBins ( nBins_ ) , theLowerBound ( lowerBound_ ) , theUpperBound ( upperBound_ ) ,
153  theStatistics ( statistics_ ) , thePlotLog ( plotLog_ ) , thePlotNormalized ( plotNormalized_ ) ,
154  thePlotFirst ( plotFirst_ ), theMin(-1.), theMax(-1.), mcPlots_(mc)
155 {
156  // defaults for array dimensions
157  theArrayDimension = 0 ;
158 
159  // check plot order string
160  if ( thePlotFirst == "l" || thePlotFirst == "c" || thePlotFirst == "b" ) {
161  // OK
162  }
163  else {
164  // not correct: print warning and set default (l)
165  std::cout << "FlavourHistograms::FlavourHistograms : thePlotFirst was not correct : " << thePlotFirst << std::endl ;
166  std::cout << "FlavourHistograms::FlavourHistograms : Set it to default value (l)! " << std::endl ;
167  thePlotFirst = "l" ;
168  }
169 
170  if (!update) {
171  // book histos
172  HistoProviderDQM prov("Btag",folder);
173  theHisto_all = (prov.book1D( theBaseNameTitle + "ALL" , theBaseNameDescription + " all jets" , theNBins , theLowerBound , theUpperBound )) ;
174  if (mcPlots_) {
175  theHisto_d = (prov.book1D ( theBaseNameTitle + "D" , theBaseNameDescription + " d-jets" , theNBins , theLowerBound , theUpperBound )) ;
176  theHisto_u = (prov.book1D ( theBaseNameTitle + "U" , theBaseNameDescription + " u-jets" , theNBins , theLowerBound , theUpperBound )) ;
177  theHisto_s = (prov.book1D ( theBaseNameTitle + "S" , theBaseNameDescription + " s-jets" , theNBins , theLowerBound , theUpperBound )) ;
178  theHisto_c = (prov.book1D ( theBaseNameTitle + "C" , theBaseNameDescription + " c-jets" , theNBins , theLowerBound , theUpperBound )) ;
179  theHisto_b = (prov.book1D ( theBaseNameTitle + "B" , theBaseNameDescription + " b-jets" , theNBins , theLowerBound , theUpperBound )) ;
180  theHisto_g = (prov.book1D ( theBaseNameTitle + "G" , theBaseNameDescription + " g-jets" , theNBins , theLowerBound , theUpperBound )) ;
181  theHisto_ni = (prov.book1D ( theBaseNameTitle + "NI" , theBaseNameDescription + " ni-jets" , theNBins , theLowerBound , theUpperBound )) ;
182  theHisto_dus = (prov.book1D ( theBaseNameTitle + "DUS" , theBaseNameDescription + " dus-jets" , theNBins , theLowerBound , theUpperBound )) ;
183  theHisto_dusg = (prov.book1D ( theBaseNameTitle + "DUSG" , theBaseNameDescription + " dusg-jets" , theNBins , theLowerBound , theUpperBound )) ;
184  }else{
185  theHisto_d = 0;
186  theHisto_u = 0;
187  theHisto_s = 0;
188  theHisto_c = 0;
189  theHisto_b = 0;
190  theHisto_g = 0;
191  theHisto_ni = 0;
192  theHisto_dus = 0;
193  theHisto_dusg = 0;
194  }
195  // statistics if requested
196  if ( theStatistics ) {
197  theHisto_all ->getTH1F()->Sumw2() ;
198  if (mcPlots_ ) {
199 
200  theHisto_d ->getTH1F()->Sumw2() ;
201  theHisto_u ->getTH1F()->Sumw2() ;
202  theHisto_s ->getTH1F()->Sumw2() ;
203  theHisto_c ->getTH1F()->Sumw2() ;
204  theHisto_b ->getTH1F()->Sumw2() ;
205  theHisto_g ->getTH1F()->Sumw2() ;
206  theHisto_ni ->getTH1F()->Sumw2() ;
207  theHisto_dus ->getTH1F()->Sumw2() ;
208  theHisto_dusg->getTH1F()->Sumw2() ;
209  }
210  }
211  } else {
212  HistoProviderDQM prov("Btag",folder);
213  theHisto_all = prov.access(theBaseNameTitle + "ALL" ) ;
214  if (mcPlots_) {
215 
216  theHisto_d = prov.access(theBaseNameTitle + "D" ) ;
217  theHisto_u = prov.access(theBaseNameTitle + "U" ) ;
218  theHisto_s = prov.access(theBaseNameTitle + "S" ) ;
219  theHisto_c = prov.access(theBaseNameTitle + "C" ) ;
220  theHisto_b =prov.access(theBaseNameTitle + "B" ) ;
221  theHisto_g =prov.access(theBaseNameTitle + "G" ) ;
222  theHisto_ni =prov.access(theBaseNameTitle + "NI" ) ;
223  theHisto_dus =prov.access(theBaseNameTitle + "DUS" ) ;
224  theHisto_dusg =prov.access(theBaseNameTitle + "DUSG") ;
225  }
226  }
227 
228  // defaults for other data members
229  theCanvas = 0 ;
230 }
MonitorElement * theHisto_g
MonitorElement * theHisto_c
MonitorElement * theHisto_u
MonitorElement * theHisto_all
MonitorElement * theHisto_d
MonitorElement * theHisto_dusg
std::string theBaseNameDescription
TH1F * getTH1F(void) const
MonitorElement * theHisto_b
std::string theBaseNameTitle
#define update(a, b)
tuple cout
Definition: gather_cfg.py:121
MonitorElement * theHisto_s
MonitorElement * theHisto_dus
MonitorElement * theHisto_ni
template<class T >
FlavourHistograms< T >::~FlavourHistograms ( )
virtual

Definition at line 234 of file FlavourHistorgrams.h.

234  {
235  // delete the canvas*/
236  delete theCanvas ;
237 }
template<class T>
FlavourHistograms< T >::FlavourHistograms ( )
inlineprivate

Definition at line 138 of file FlavourHistorgrams.h.

138 {}

Member Function Documentation

template<class T>
int* FlavourHistograms< T >::arrayDimension ( ) const
inline

Definition at line 75 of file FlavourHistorgrams.h.

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

75 { return theArrayDimension; }
template<class T>
std::string FlavourHistograms< T >::baseNameDescription ( ) const
inline

Definition at line 67 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

67 { return theBaseNameDescription ; }
std::string theBaseNameDescription
template<class T>
std::string FlavourHistograms< T >::baseNameTitle ( ) const
inline

Definition at line 66 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

66 { return theBaseNameTitle ; }
std::string theBaseNameTitle
template<class T>
void FlavourHistograms< T >::divide ( const FlavourHistograms< T > &  bHD) const

Definition at line 452 of file FlavourHistorgrams.h.

References MultipleCompare::Divide(), getTH1F(), 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_s(), and FlavourHistograms< T >::histo_u().

Referenced by python.rootplot.utilities.Hist::__div__(), and EffPurFromHistos::EffPurFromHistos().

452  {
453  // divide histos using binomial errors
454  //
455  // ATTENTION: It's the responsability of the user to make sure that the HistoDescriptions
456  // involved in this operation have been constructed with the statistics option switched on!!
457  //
458  theHisto_all ->getTH1F()-> Divide ( theHisto_all->getTH1F() , bHD.histo_all () , 1.0 , 1.0 , "b" ) ;
459  if (mcPlots_) {
460  theHisto_d ->getTH1F()-> Divide ( theHisto_d ->getTH1F() , bHD.histo_d () , 1.0 , 1.0 , "b" ) ;
461  theHisto_u ->getTH1F()-> Divide ( theHisto_u ->getTH1F() , bHD.histo_u () , 1.0 , 1.0 , "b" ) ;
462  theHisto_s ->getTH1F()-> Divide ( theHisto_s ->getTH1F() , bHD.histo_s () , 1.0 , 1.0 , "b" ) ;
463  theHisto_c ->getTH1F()-> Divide ( theHisto_c ->getTH1F() , bHD.histo_c () , 1.0 , 1.0 , "b" ) ;
464  theHisto_b ->getTH1F()-> Divide ( theHisto_b ->getTH1F() , bHD.histo_b () , 1.0 , 1.0 , "b" ) ;
465  theHisto_g ->getTH1F()-> Divide ( theHisto_g ->getTH1F() , bHD.histo_g () , 1.0 , 1.0 , "b" ) ;
466  theHisto_ni ->getTH1F()-> Divide ( theHisto_ni->getTH1F() , bHD.histo_ni () , 1.0 , 1.0 , "b" ) ;
467  theHisto_dus ->getTH1F()-> Divide ( theHisto_dus->getTH1F() , bHD.histo_dus () , 1.0 , 1.0 , "b" ) ;
468  theHisto_dusg ->getTH1F()-> Divide ( theHisto_dusg->getTH1F() , bHD.histo_dusg() , 1.0 , 1.0 , "b" ) ;
469  }
470 }
MonitorElement * theHisto_g
MonitorElement * theHisto_c
TH1F * histo_all() const
MonitorElement * theHisto_u
TH1F * histo_c() const
MonitorElement * theHisto_all
TH1F * histo_dus() const
TH1F * histo_b() const
TH1F * histo_g() const
TH1F * histo_d() const
TH1F * getTH1F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
MonitorElement * theHisto_d
MonitorElement * theHisto_dusg
TH1F * getTH1F(void) const
MonitorElement * theHisto_b
TH1F * histo_s() const
TH1F * histo_ni() const
TH1F * histo_u() const
MonitorElement * theHisto_s
TH1F * histo_dusg() const
MonitorElement * theHisto_dus
MonitorElement * theHisto_ni
template<class T >
void FlavourHistograms< T >::epsPlot ( const std::string &  name)
template<class T>
void FlavourHistograms< T >::fill ( const int &  flavour,
const T variable 
) const

Definition at line 251 of file FlavourHistorgrams.h.

Referenced by SoftLeptonTagPlotter::analyzeTag(), TrackCountingTagPlotter::analyzeTag(), TrackIPTagPlotter::analyzeTag(), JetTagPlotter::analyzeTag(), and TrackIPHistograms< T >::fill().

252 {
253  // For single variables and arrays (for arrays only a single index can be filled)
254  fillVariable ( flavour , variable ) ;
255 }
void fillVariable(const int &flavour, const T &var) const
int flavour(const Candidate &part)
Definition: pdgIdUtils.h:31
template<class T>
void FlavourHistograms< T >::fill ( const int &  flavour,
const T variable 
) const

Definition at line 258 of file FlavourHistorgrams.h.

References i.

259 {
260  if ( theArrayDimension == 0 ) {
261  // single variable
262  fillVariable ( flavour , *variable ) ;
263  } else {
264  // array
266  //
267  for ( int i = 0 ; i != iMax ; ++i ) {
268  // check if only one index to be plotted (<0: switched off -> plot all)
269  if ( ( theIndexToPlot < 0 ) || ( i == theIndexToPlot ) ) {
270  fillVariable ( flavour , *(variable+i) ) ;
271  }
272  }
273 
274  // if single index to be filled but not enough entries: fill 0.0 (convention!)
275  if ( theIndexToPlot >= iMax ) {
276  // cout << "==>> The index to be filled is too big -> fill 0.0 : " << theBaseNameTitle << " : " << theIndexToPlot << " >= " << iMax << endl ;
277  const T& theZero = static_cast<T> ( 0.0 ) ;
278  fillVariable ( flavour , theZero ) ;
279  }
280  }
281 }
int i
Definition: DBlmapReader.cc:9
void fillVariable(const int &flavour, const T &var) const
long double T
int flavour(const Candidate &part)
Definition: pdgIdUtils.h:31
template<class T>
void FlavourHistograms< T >::fillVariable ( const int &  flavour,
const T var 
) const
protected

Definition at line 474 of file FlavourHistorgrams.h.

474  {
475  // all
476  theHisto_all ->Fill ( var ) ;
477  // flavour specific
478  if (!mcPlots_) return;
479 
480  switch(flavour) {
481  case 1:
482  theHisto_d->Fill( var );
483  theHisto_dus->Fill( var );
484  theHisto_dusg->Fill( var );
485  return;
486  case 2:
487  theHisto_u->Fill( var );
488  theHisto_dus->Fill( var );
489  theHisto_dusg->Fill( var );
490  return;
491  case 3:
492  theHisto_s->Fill( var );
493  theHisto_dus->Fill( var );
494  theHisto_dusg->Fill( var );
495  return;
496  case 4:
497  theHisto_c->Fill( var );
498  return;
499  case 5:
500  theHisto_b->Fill( var );
501  return;
502  case 21:
503  theHisto_g->Fill( var );
504  theHisto_dusg->Fill( var );
505  return;
506  default:
507  theHisto_ni->Fill( var );
508  return;
509  }
510 }
MonitorElement * theHisto_g
MonitorElement * theHisto_c
MonitorElement * theHisto_u
MonitorElement * theHisto_all
void Fill(long long x)
MonitorElement * theHisto_d
MonitorElement * theHisto_dusg
MonitorElement * theHisto_b
MonitorElement * theHisto_s
int flavour(const Candidate &part)
Definition: pdgIdUtils.h:31
MonitorElement * theHisto_dus
MonitorElement * theHisto_ni
template<class T >
std::vector< TH1F * > FlavourHistograms< T >::getHistoVector ( ) const

Definition at line 513 of file FlavourHistorgrams.h.

References getTH1F().

Referenced by EffPurFromHistos::EffPurFromHistos().

514 {
515  std::vector<TH1F*> histoVector;
516  histoVector.push_back ( theHisto_all->getTH1F() );
517  if (mcPlots_) {
518  histoVector.push_back ( theHisto_d->getTH1F() );
519  histoVector.push_back ( theHisto_u->getTH1F() );
520  histoVector.push_back ( theHisto_s->getTH1F() );
521  histoVector.push_back ( theHisto_c->getTH1F() );
522  histoVector.push_back ( theHisto_b->getTH1F() );
523  histoVector.push_back ( theHisto_g ->getTH1F() );
524  histoVector.push_back ( theHisto_ni->getTH1F() );
525  histoVector.push_back ( theHisto_dus->getTH1F() );
526  histoVector.push_back ( theHisto_dusg->getTH1F());
527  }
528  return histoVector;
529 }
MonitorElement * theHisto_g
MonitorElement * theHisto_c
MonitorElement * theHisto_u
MonitorElement * theHisto_all
TH1F * getTH1F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
MonitorElement * theHisto_d
MonitorElement * theHisto_dusg
TH1F * getTH1F(void) const
MonitorElement * theHisto_b
MonitorElement * theHisto_s
MonitorElement * theHisto_dus
MonitorElement * theHisto_ni
template<class T>
TH1F* FlavourHistograms< T >::histo_all ( ) const
inline

Definition at line 80 of file FlavourHistorgrams.h.

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

80 { return theHisto_all->getTH1F() ; }
MonitorElement * theHisto_all
TH1F * getTH1F(void) const
template<class T>
TH1F* FlavourHistograms< T >::histo_b ( ) const
inline

Definition at line 85 of file FlavourHistorgrams.h.

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

85 { return theHisto_b->getTH1F() ; }
TH1F * getTH1F(void) const
MonitorElement * theHisto_b
template<class T>
TH1F* FlavourHistograms< T >::histo_c ( ) const
inline

Definition at line 84 of file FlavourHistorgrams.h.

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

84 { return theHisto_c->getTH1F() ; }
MonitorElement * theHisto_c
TH1F * getTH1F(void) const
template<class T>
TH1F* FlavourHistograms< T >::histo_d ( ) const
inline

Definition at line 81 of file FlavourHistorgrams.h.

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

81 { return theHisto_d ->getTH1F() ; }
MonitorElement * theHisto_d
TH1F * getTH1F(void) const
template<class T>
TH1F* FlavourHistograms< T >::histo_dus ( ) const
inline

Definition at line 88 of file FlavourHistorgrams.h.

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

88 { return theHisto_dus->getTH1F() ; }
TH1F * getTH1F(void) const
MonitorElement * theHisto_dus
template<class T>
TH1F* FlavourHistograms< T >::histo_dusg ( ) const
inline

Definition at line 89 of file FlavourHistorgrams.h.

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

89 { return theHisto_dusg->getTH1F() ; }
MonitorElement * theHisto_dusg
TH1F * getTH1F(void) const
template<class T>
TH1F* FlavourHistograms< T >::histo_g ( ) const
inline

Definition at line 86 of file FlavourHistorgrams.h.

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

86 { return theHisto_g->getTH1F() ; }
MonitorElement * theHisto_g
TH1F * getTH1F(void) const
template<class T>
TH1F* FlavourHistograms< T >::histo_ni ( ) const
inline

Definition at line 87 of file FlavourHistorgrams.h.

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

87 { return theHisto_ni->getTH1F() ; }
TH1F * getTH1F(void) const
MonitorElement * theHisto_ni
template<class T>
TH1F* FlavourHistograms< T >::histo_s ( ) const
inline

Definition at line 83 of file FlavourHistorgrams.h.

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

83 { return theHisto_s->getTH1F() ; }
TH1F * getTH1F(void) const
MonitorElement * theHisto_s
template<class T>
TH1F* FlavourHistograms< T >::histo_u ( ) const
inline

Definition at line 82 of file FlavourHistorgrams.h.

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

82 { return theHisto_u->getTH1F() ; }
MonitorElement * theHisto_u
TH1F * getTH1F(void) const
template<class T>
int FlavourHistograms< T >::indexToPlot ( ) const
inline

Definition at line 77 of file FlavourHistorgrams.h.

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

77 {return theIndexToPlot; }
template<class T>
double FlavourHistograms< T >::lowerBound ( ) const
inline

Definition at line 69 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

69 { return theLowerBound ; }
template<class T>
int FlavourHistograms< T >::maxDimension ( ) const
inline

Definition at line 76 of file FlavourHistorgrams.h.

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

76 {return theMaxDimension; }
template<class T>
int FlavourHistograms< T >::nBins ( ) const
inline

Definition at line 68 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

68 { return theNBins ; }
template<class T >
void FlavourHistograms< T >::plot ( TPad *  theCanvas = 0)

Definition at line 303 of file FlavourHistorgrams.h.

References getTH1F(), MonitorElement::getTH1F(), interpolateCardsSimple::histo, i, max(), and RecoBTag::setTDRStyle().

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

303  {
304 
305 //fixme:
306  bool btppNI = false;
307  bool btppColour = true;
308 
309  if (theCanvas)
310  theCanvas->cd();
311 
312  RecoBTag::setTDRStyle()->cd();
313  gPad->UseCurrentStyle();
314 // if ( !btppTitle ) gStyle->SetOptTitle ( 0 ) ;
315 //
316 // // here: plot histograms in a canvas
317 // theCanvas = new TCanvas ( "C" + theBaseNameTitle , "C" + theBaseNameDescription ,
318 // btppXCanvas , btppYCanvas ) ;
319 // theCanvas->SetFillColor ( 0 ) ;
320 // theCanvas->cd ( 1 ) ;
321  gPad->SetLogy ( 0 ) ;
322  if ( thePlotLog ) gPad->SetLogy ( 1 ) ;
323  gPad->SetGridx ( 0 ) ;
324  gPad->SetGridy ( 0 ) ;
325  gPad->SetTitle ( 0 ) ;
326 
327  MonitorElement * histo[4];
328  int col[4], lineStyle[4], markerStyle[4];
329  int lineWidth = 1 ;
330 
331  const double markerSize = gPad->GetWh() * gPad->GetHNDC() / 500.;
332 
333  // default (l)
334  histo[0] = theHisto_dusg ;
335  //CW histo_1 = theHisto_dus ;
336  histo[1] = theHisto_b ;
337  histo[2] = theHisto_c ;
338  histo[3]= 0 ;
339 
340  double max = theMax;
341  if (theMax<=0.) {
342  max = theHisto_dusg->getTH1F()->GetMaximum();
343  if (theHisto_b->getTH1F()->GetMaximum() > max) max = theHisto_b->getTH1F()->GetMaximum();
344  if (theHisto_c->getTH1F()->GetMaximum() > max) max = theHisto_c->getTH1F()->GetMaximum();
345  }
346 
347  if (btppNI) {
348  histo[3] = theHisto_ni ;
349  if (theHisto_ni->getTH1F()->GetMaximum() > max) max = theHisto_ni->getTH1F()->GetMaximum();
350  }
351 
352  if ( btppColour ) { // print colours
353  col[0] = 4 ;
354  col[1] = 2 ;
355  col[2] = 6 ;
356  col[3] = 3 ;
357  lineStyle[0] = 1 ;
358  lineStyle[1] = 1 ;
359  lineStyle[2] = 1 ;
360  lineStyle[3] = 1 ;
361  markerStyle[0] = 20 ;
362  markerStyle[1] = 21 ;
363  markerStyle[2] = 22 ;
364  markerStyle[3] = 23 ;
365  lineWidth = 1 ;
366  }
367  else { // different marker/line styles
368  col[1] = 1 ;
369  col[2] = 1 ;
370  col[3] = 1 ;
371  col[0] = 1 ;
372  lineStyle[0] = 2 ;
373  lineStyle[1] = 1 ;
374  lineStyle[2] = 3 ;
375  lineStyle[3] = 4 ;
376  markerStyle[0] = 20 ;
377  markerStyle[1] = 21 ;
378  markerStyle[2] = 22 ;
379  markerStyle[3] = 23 ;
380  }
381 
382  // if changing order (NI stays always last)
383 
384  // c to plot first
385  if ( thePlotFirst == "c" ) {
386  histo[0] = theHisto_c ;
387  if ( btppColour ) col[0] = 6 ;
388  if ( !btppColour ) lineStyle[0] = 3 ;
389  histo[2] = theHisto_dusg ;
390  if ( btppColour ) col[2] = 4 ;
391  if ( !btppColour ) lineStyle[2] = 2 ;
392  }
393 
394  // b to plot first
395  if ( thePlotFirst == "b" ) {
396  histo[0] = theHisto_b ;
397  if ( btppColour ) col[0] = 2 ;
398  if ( !btppColour ) lineStyle[0] = 1 ;
399  histo[1] = theHisto_dusg ;
400  if ( btppColour ) col[1] = 4 ;
401  if ( !btppColour ) lineStyle[1] = 2 ;
402  }
403 
404 
405  histo[0] ->getTH1F()->GetXaxis()->SetTitle ( theBaseNameDescription.c_str() ) ;
406  histo[0] ->getTH1F()->GetYaxis()->SetTitle ( "Arbitrary Units" ) ;
407  histo[0] ->getTH1F()->GetYaxis()->SetTitleOffset(1.25) ;
408 
409  for (int i=0; i != 4; ++i) {
410  if (histo[i]== 0 ) continue;
411  histo[i] ->getTH1F()->SetStats ( false ) ;
412  histo[i] ->getTH1F()->SetLineStyle ( lineStyle[i] ) ;
413  histo[i] ->getTH1F()->SetLineWidth ( lineWidth ) ;
414  histo[i] ->getTH1F()->SetLineColor ( col[i] ) ;
415  histo[i] ->getTH1F()->SetMarkerStyle ( markerStyle[i] ) ;
416  histo[i] ->getTH1F()->SetMarkerColor ( col[i] ) ;
417  histo[i] ->getTH1F()->SetMarkerSize ( markerSize ) ;
418  }
419 
420  if ( thePlotNormalized ) {
421 // cout <<histo[0]->GetEntries() <<" "<< histo[1]->GetEntries() <<" "<< histo[2]->GetEntries()<<" "<<histo[3]->GetEntries()<<endl;
422  if (histo[0]->getTH1F()->GetEntries() != 0) {histo[0]->getTH1F() ->DrawNormalized() ;} else { histo[0]->getTH1F() ->SetMaximum(1.0);
423 histo[0] ->getTH1F()->Draw() ;}
424  if (histo[1]->getTH1F()->GetEntries() != 0) histo[1] ->getTH1F()->DrawNormalized("Same") ;
425  if (histo[2]->getTH1F()->GetEntries() != 0) histo[2]->getTH1F() ->DrawNormalized("Same") ;
426  if ((histo[3] != 0) && (histo[3]->getTH1F()->GetEntries() != 0)) histo[3] ->getTH1F()->DrawNormalized("Same") ;
427  }
428  else {
429  histo[0]->getTH1F()->SetMaximum(max*1.05);
430  if (theMin!=-1.) histo[0]->getTH1F()->SetMinimum(theMin);
431  histo[0]->getTH1F()->Draw() ;
432  histo[1]->getTH1F()->Draw("Same") ;
433  histo[2]->getTH1F()->Draw("Same") ;
434  if ( histo[3] != 0 ) histo[3]->getTH1F()->Draw("Same") ;
435  }
436 
437 }
int i
Definition: DBlmapReader.cc:9
MonitorElement * theHisto_c
TStyle * setTDRStyle()
Definition: Tools.cc:287
const T & max(const T &a, const T &b)
TH1F * getTH1F(std::string name, std::string process, std::string rootfolder, DQMStore *dbe_, bool verb, bool clone)
MonitorElement * theHisto_dusg
std::string theBaseNameDescription
TH1F * getTH1F(void) const
MonitorElement * theHisto_b
MonitorElement * theHisto_ni
template<class T>
std::string FlavourHistograms< T >::plotFirst ( ) const
inline

Definition at line 74 of file FlavourHistorgrams.h.

74 { return thePlotFirst ; }
template<class T>
bool FlavourHistograms< T >::plotLog ( ) const
inline

Definition at line 72 of file FlavourHistorgrams.h.

72 { return thePlotLog ; }
template<class T>
bool FlavourHistograms< T >::plotNormalized ( ) const
inline

Definition at line 73 of file FlavourHistorgrams.h.

73 { return thePlotNormalized ; }
template<class T>
void FlavourHistograms< T >::SetMaximum ( const double &  max)
inline

Definition at line 61 of file FlavourHistorgrams.h.

61 { theMax = max;}
const T & max(const T &a, const T &b)
template<class T>
void FlavourHistograms< T >::SetMinimum ( const double &  min)
inline

Definition at line 62 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

62 { theMin = min;}
#define min(a, b)
Definition: mlp_lapack.h:161
template<class T >
void FlavourHistograms< T >::settitle ( const char *  title)

Definition at line 285 of file FlavourHistorgrams.h.

Referenced by JetTagPlotter::JetTagPlotter(), and TrackIPHistograms< T >::settitle().

285  {
287  if (mcPlots_) {
288 
298  }
299 }
MonitorElement * theHisto_g
MonitorElement * theHisto_c
MonitorElement * theHisto_u
MonitorElement * theHisto_all
MonitorElement * theHisto_d
MonitorElement * theHisto_dusg
MonitorElement * theHisto_b
MonitorElement * theHisto_s
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
MonitorElement * theHisto_dus
MonitorElement * theHisto_ni
template<class T>
bool FlavourHistograms< T >::statistics ( ) const
inline

Definition at line 71 of file FlavourHistorgrams.h.

71 { return theStatistics ; }
template<class T>
double FlavourHistograms< T >::upperBound ( ) const
inline

Definition at line 70 of file FlavourHistorgrams.h.

Referenced by EffPurFromHistos::EffPurFromHistos().

70 { return theUpperBound ; }

Member Data Documentation

template<class T>
bool FlavourHistograms< T >::mcPlots_
private

Definition at line 140 of file FlavourHistorgrams.h.

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

template<class T>
int* FlavourHistograms< T >::theArrayDimension
protected
template<class T>
std::string FlavourHistograms< T >::theBaseNameDescription
protected
template<class T>
std::string FlavourHistograms< T >::theBaseNameTitle
protected
template<class T>
TCanvas* FlavourHistograms< T >::theCanvas
protected

Definition at line 136 of file FlavourHistorgrams.h.

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

template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_all
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_b
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_c
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_d
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_dus
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_dusg
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_g
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_ni
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_s
protected
template<class T>
MonitorElement* FlavourHistograms< T >::theHisto_u
protected
template<class T>
int FlavourHistograms< T >::theIndexToPlot
protected

Definition at line 107 of file FlavourHistorgrams.h.

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

template<class T>
double FlavourHistograms< T >::theLowerBound
protected
template<class T>
double FlavourHistograms< T >::theMax
protected

Definition at line 118 of file FlavourHistorgrams.h.

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

template<class T>
int FlavourHistograms< T >::theMaxDimension
protected

Definition at line 106 of file FlavourHistorgrams.h.

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

template<class T>
double FlavourHistograms< T >::theMin
protected

Definition at line 118 of file FlavourHistorgrams.h.

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

template<class T>
int FlavourHistograms< T >::theNBins
protected
template<class T>
std::string FlavourHistograms< T >::thePlotFirst
protected
template<class T>
bool FlavourHistograms< T >::thePlotLog
protected

Definition at line 115 of file FlavourHistorgrams.h.

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

template<class T>
bool FlavourHistograms< T >::thePlotNormalized
protected

Definition at line 116 of file FlavourHistorgrams.h.

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

template<class T>
bool FlavourHistograms< T >::theStatistics
protected
template<class T>
double FlavourHistograms< T >::theUpperBound
protected