#include <HLTExoticaPlotter.h>

Public Member Functions
void	analyze (const bool &isPassTrigger, const std::string &source, const std::vector< reco::LeafCandidate > &matches, std::map< int, double > theSumEt, std::vector< float > &dxys)

void	beginJob ()

void	beginRun (const edm::Run &, const edm::EventSetup &)

const std::string	gethltpath () const

	HLTExoticaPlotter (const edm::ParameterSet &pset, const std::string &hltPath, const std::vector< unsigned int > &objectsType)

void	plotterBookHistos (DQMStore::IBooker &iBooker, const edm::Run &iRun, const edm::EventSetup &iSetup)

	~HLTExoticaPlotter ()

Private Member Functions
void	bookHist (DQMStore::IBooker &iBooker, const std::string &source, const std::string &objType, const std::string &variable)

void	fillHist (const bool &passTrigger, const std::string &source, const std::string &objType, const std::string &var, const float &value)

Private Attributes
std::map< std::string, MonitorElement * >	_elements

std::string	_hltPath

std::string	_hltProcessName

unsigned int	_nObjects

std::set< unsigned int >	_objectsType

std::vector< double >	_parametersDxy

std::vector< double >	_parametersEta

std::vector< double >	_parametersPhi

std::vector< double >	_parametersTurnOn

std::vector< double >	_parametersTurnOnSumEt

Detailed Description

Generate histograms for trigger efficiencies Exotica related Documentation available on the CMS TWiki: https://twiki.cern.ch/twiki/bin/view/CMS/EXOTICATriggerValidation

Author: Thiago R. Fernandez Perez Tomei Based and adapted from: J. Duarte Campderros code from HLTriggerOffline/Higgs J. Klukas, M. Vander Donckt and J. Alcaraz code from the HLTriggerOffline/Muon package.

Definition at line 44 of file HLTExoticaPlotter.h.

Constructor & Destructor Documentation

HLTExoticaPlotter::HLTExoticaPlotter	(	const edm::ParameterSet &	pset,
		const std::string &	hltPath,
		const std::vector< unsigned int > &	objectsType
	)

Definition at line 20 of file HLTExoticaPlotter.cc.

References LogDebug.

                                                                                   :
     _hltPath(hltPath),
     _hltProcessName(pset.getParameter<std::string>("hltProcessName")),
     _objectsType(std::set<unsigned int>(objectsType.begin(), objectsType.end())),
     _nObjects(objectsType.size()),
     _parametersEta(pset.getParameter<std::vector<double> >("parametersEta")),
     _parametersPhi(pset.getParameter<std::vector<double> >("parametersPhi")),
     _parametersTurnOn(pset.getParameter<std::vector<double> >("parametersTurnOn")),
     _parametersTurnOnSumEt(pset.getParameter<std::vector<double> >("parametersTurnOnSumEt")),
     _parametersDxy(pset.getParameter<std::vector<double> >("parametersDxy"))
 {
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::constructor()";
 }

HLTExoticaPlotter::~HLTExoticaPlotter ( )

Definition at line 36 of file HLTExoticaPlotter.cc.

37 {

38 }

Member Function Documentation

void HLTExoticaPlotter::analyze	(	const bool &	isPassTrigger,
		const std::string &	source,
		const std::vector< reco::LeafCandidate > &	matches,
		std::map< int, double >	theSumEt,
		std::vector< float > &	dxys
	)

Definition at line 107 of file HLTExoticaPlotter.cc.

References _objectsType, EVTColContainer::ELEC, eta, fillHist(), EVTColContainer::getTypeString(), j, LogDebug, EVTColContainer::MUON, phi, EnergyCorrector::pt, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::analyze()";
     if (!isPassTrigger) {
         return;
     }
 
     std::map<unsigned int, int> countobjects;
     // Initializing the count of the used object
     for (std::set<unsigned int>::iterator co = _objectsType.begin();
          co != _objectsType.end(); ++co) {
         countobjects[*co] = 0;
     }
 
     int counttotal = 0;
 
     // 3 : pt1, pt2, pt3
     const int totalobjectssize3 = 3 * countobjects.size();
     // Fill the histos if pass the trigger (just the two with higher pt)
     for (size_t j = 0; j < matches.size(); ++j) {
         // Is this object owned by this trigger? If not we are not interested...
         if (_objectsType.find(matches[j].pdgId()) == _objectsType.end()) {
             continue;
         }
 
         const unsigned int objType = matches[j].pdgId();
         const std::string objTypeStr = EVTColContainer::getTypeString(objType);
 
         float pt  =   matches[j].pt();
         float eta =   matches[j].eta();
         float phi =   matches[j].phi();
 
     if ( !( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) {
           this->fillHist(isPassTrigger, source, objTypeStr, "Eta", eta);
           this->fillHist(isPassTrigger, source, objTypeStr, "Phi", phi);
     }
     else if( source!="gen" ) {
       if(theSumEt[objType]>=0 && countobjects[objType] == 0) {
         this->fillHist(isPassTrigger, source, objTypeStr, "SumEt", theSumEt[objType]);
       }
     }
 
         if ( dxys.size() && (objType == EVTColContainer::ELEC || objType == EVTColContainer::MUON ) ) 
            this->fillHist(isPassTrigger, source, objTypeStr, "Dxy", dxys[j] );
 
         if (countobjects[objType] == 0) {
       if ( !( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) || source!="gen" ) {
             this->fillHist(isPassTrigger, source, objTypeStr, "MaxPt1", pt);
       }
       // Filled the high pt ...
       ++(countobjects[objType]);
       ++counttotal;
         } 
     else if (countobjects[objType] == 1) {
       if( !( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) {
             this->fillHist(isPassTrigger, source, objTypeStr, "MaxPt2", pt);
       }
       // Filled the second high pt ...
       ++(countobjects[objType]);
       ++counttotal;
         } 
     else if (countobjects[objType] == 2) {
       if( !( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) {
             this->fillHist(isPassTrigger, source, objTypeStr, "MaxPt3", pt);
       }
       // Filled the third highest pt ...
       ++(countobjects[objType]);
       ++counttotal;
         } 
     else {
       if (counttotal == totalobjectssize3) {
         break;
       }
         }
 
     } // end loop over matches
 }

void HLTExoticaPlotter::beginJob ( void )

Definition at line 41 of file HLTExoticaPlotter.cc.

42 {

43 }

void HLTExoticaPlotter::beginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

void HLTExoticaPlotter::bookHist	(	DQMStore::IBooker &	iBooker,
		const std::string &	source,
		const std::string &	objType,
		const std::string &	variable
	)

private

Definition at line 190 of file HLTExoticaPlotter.cc.

References _elements, _hltPath, _parametersDxy, _parametersEta, _parametersPhi, _parametersTurnOn, _parametersTurnOnSumEt, DQMStore::IBooker::book1D(), prof2calltree::edges, h, i, LogDebug, bookConverter::max, min(), mergeVDriftHistosByStation::name, source, AlCaHLTBitMon_QueryRunRegistry::string, and indexGen::title.

Referenced by plotterBookHistos().

 {
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::bookHist()";
     std::string sourceUpper = source;
     sourceUpper[0] = std::toupper(sourceUpper[0]);
     std::string name = source + objType + variable + "_" + _hltPath;
     TH1F * h = 0;
 
     if (variable.find("SumEt") != std::string::npos) {
         std::string title = "Sum ET of " + sourceUpper + " " + objType;
         const size_t nBins = _parametersTurnOnSumEt.size() - 1;
         float * edges = new float[nBins + 1];
         for (size_t i = 0; i < nBins + 1; i++) {
             edges[i] = _parametersTurnOnSumEt[i];
         }
         h = new TH1F(name.c_str(), title.c_str(), nBins, edges);
         delete[] edges;
     }
     else if (variable.find("Dxy") != std::string::npos) {
       std::string title = "Dxy " + sourceUpper + " " + objType;
       int    nBins = _parametersDxy[0];
       double min   = _parametersDxy[1];
       double max   = _parametersDxy[2];
       h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
     }
     else if (variable.find("MaxPt") != std::string::npos) {
         std::string desc = (variable == "MaxPt1") ? "Leading" : 
                            (variable == "MaxPt2") ? "Next-to-Leading" : "Next-to-next-to-Leading";
         std::string title = "pT of " + desc + " " + sourceUpper + " " + objType + " "
                             "where event pass the " + _hltPath;
         const size_t nBins = _parametersTurnOn.size() - 1;
         float * edges = new float[nBins + 1];
         for (size_t i = 0; i < nBins + 1; i++) {
             edges[i] = _parametersTurnOn[i];
         }
         h = new TH1F(name.c_str(), title.c_str(), nBins, edges);
         delete [] edges;
     } 
 
     else {
         std::string symbol = (variable == "Eta") ? "#eta" : "#phi";
         std::string title  = symbol + " of " + sourceUpper + " " + objType + " " +
                              "where event pass the " + _hltPath;
         std::vector<double> params = (variable == "Eta") ? _parametersEta : _parametersPhi;
 
         int    nBins = (int)params[0];
         double min   = params[1];
         double max   = params[2];
         h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
     }
 
     h->Sumw2();
     _elements[name] = iBooker.book1D(name, h);
     //    LogDebug("ExoticaValidation") << "                        booked histo with name " << name << "\n"
     //                << "                        at location " << (unsigned long int)_elements[name];
     delete h;
 }

void HLTExoticaPlotter::fillHist	(	const bool &	passTrigger,
		const std::string &	source,
		const std::string &	objType,
		const std::string &	var,
		const float &	value
	)

private

Definition at line 251 of file HLTExoticaPlotter.cc.

References _elements, _hltPath, LogDebug, mergeVDriftHistosByStation::name, source, AlCaHLTBitMon_QueryRunRegistry::string, and relativeConstraints::value.

Referenced by analyze().

 {
     std::string sourceUpper = source;
     sourceUpper[0] = toupper(sourceUpper[0]);
     std::string name = source + objType + variable + "_" + _hltPath;
 
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::fillHist()" << name << " " << value;
     _elements[name]->Fill(value);
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::fillHist()" << name << " worked";
 }

const std::string HLTExoticaPlotter::gethltpath ( ) const

inline

Definition at line 57 of file HLTExoticaPlotter.h.

References _hltPath.

     {
         return _hltPath;
     }

void HLTExoticaPlotter::plotterBookHistos	(	DQMStore::IBooker &	iBooker,
		const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

Definition at line 47 of file HLTExoticaPlotter.cc.

References _objectsType, bookHist(), EVTColContainer::ELEC, EVTColContainer::getTypeString(), i, LogDebug, EVTColContainer::MUON, source, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
     LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::plotterBookHistos()";
     for (std::set<unsigned int>::iterator it = _objectsType.begin();
          it != _objectsType.end(); ++it) {
         std::vector<std::string> sources(2);
         sources[0] = "gen";
         sources[1] = "rec";
 
         const std::string objTypeStr = EVTColContainer::getTypeString(*it);
 
         for (size_t i = 0; i < sources.size(); i++) {
             std::string source = sources[i];
 
             if ( source == "gen" ) {
               if ( TString(objTypeStr).Contains("MET") ||
                    TString(objTypeStr).Contains("MHT") ||
                    TString(objTypeStr).Contains("Jet")    ) {
                 continue;
               } else {
                 bookHist(iBooker, source, objTypeStr, "MaxPt1");
                 bookHist(iBooker, source, objTypeStr, "MaxPt2");
                 bookHist(iBooker, source, objTypeStr, "MaxPt3");
                 bookHist(iBooker, source, objTypeStr, "Eta");
                 bookHist(iBooker, source, objTypeStr, "Phi");
  
                 // If the target is electron or muon,
                 // we will add Dxy plots.
                 if ( *it == EVTColContainer::ELEC ||
                      *it == EVTColContainer::MUON    ) { 
                   bookHist(iBooker, source, objTypeStr, "Dxy");
                 }
               }
             } else { // reco
               if ( TString(objTypeStr).Contains("MET") ||
                    TString(objTypeStr).Contains("MHT")    ) {
                 bookHist(iBooker, source, objTypeStr, "MaxPt1");
                 bookHist(iBooker, source, objTypeStr, "SumEt");
               } else {
                 bookHist(iBooker, source, objTypeStr, "MaxPt1");
                 bookHist(iBooker, source, objTypeStr, "MaxPt2");
                 bookHist(iBooker, source, objTypeStr, "MaxPt3");
                 bookHist(iBooker, source, objTypeStr, "Eta");
                 bookHist(iBooker, source, objTypeStr, "Phi");
 
                 // If the target is electron or muon,
                 // we will add Dxy plots.
                 if ( *it == EVTColContainer::ELEC ||
                      *it == EVTColContainer::MUON    ) { 
                   bookHist(iBooker, source, objTypeStr, "Dxy");
                 }
               }
             }       
 
         }
     }
 }

Member Data Documentation

std::map<std::string, MonitorElement *> HLTExoticaPlotter::_elements

private

Definition at line 81 of file HLTExoticaPlotter.h.

Referenced by bookHist(), and fillHist().

std::string HLTExoticaPlotter::_hltPath

private

Definition at line 68 of file HLTExoticaPlotter.h.

Referenced by bookHist(), fillHist(), and gethltpath().

std::string HLTExoticaPlotter::_hltProcessName

private

Definition at line 69 of file HLTExoticaPlotter.h.

unsigned int HLTExoticaPlotter::_nObjects

private

Definition at line 73 of file HLTExoticaPlotter.h.

std::set<unsigned int> HLTExoticaPlotter::_objectsType

private

Definition at line 71 of file HLTExoticaPlotter.h.

Referenced by analyze(), and plotterBookHistos().

std::vector<double> HLTExoticaPlotter::_parametersDxy

private

Definition at line 79 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersEta

private

Definition at line 75 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersPhi

private

Definition at line 76 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersTurnOn

private

Definition at line 77 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersTurnOnSumEt

private

Definition at line 78 of file HLTExoticaPlotter.h.

Referenced by bookHist().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation