#include <HLTExoticaPlotter.h>

Public Types
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

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
bool	_drop_pt2

bool	_drop_pt3

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 42 of file HLTExoticaPlotter.h.

Member Typedef Documentation

◆ DQMStore

typedef dqm::legacy::DQMStore HLTExoticaPlotter::DQMStore

Definition at line 44 of file HLTExoticaPlotter.h.

◆ MonitorElement

typedef dqm::legacy::MonitorElement HLTExoticaPlotter::MonitorElement

Definition at line 45 of file HLTExoticaPlotter.h.

Constructor & Destructor Documentation

◆ HLTExoticaPlotter()

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

Definition at line 17 of file HLTExoticaPlotter.cc.

References _drop_pt2, _drop_pt3, LogDebug, and muonDTDigis_cfi::pset.

     : _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")),
       _drop_pt2(false),
       _drop_pt3(false) {
   if (pset.exists("dropPt2")) {
     _drop_pt2 = pset.getParameter<bool>("dropPt2");
   }
   if (pset.exists("dropPt3")) {
     _drop_pt3 = pset.getParameter<bool>("dropPt3");
   }
   LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::constructor()";
 }

◆ ~HLTExoticaPlotter()

HLTExoticaPlotter::~HLTExoticaPlotter ( )

Definition at line 40 of file HLTExoticaPlotter.cc.

40 {}

Member Function Documentation

◆ analyze()

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 101 of file HLTExoticaPlotter.cc.

References _drop_pt2, _drop_pt3, _objectsType, cms::cuda::co, EVTColContainer::ELEC, PVValHelper::eta, fillHist(), EVTColContainer::getTypeString(), dqmiolumiharvest::j, LogDebug, oniaPATMuonsWithTrigger_cff::matches, EVTColContainer::MUON, EVTColContainer::MUTRK, electrons_cff::objType, phi, DiDispStaMuonMonitor_cfi::pt, source, 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
   int totalobjectssize = 1;
   if (!_drop_pt2)
     totalobjectssize++;
   if (!_drop_pt3)
     totalobjectssize++;
   totalobjectssize *= countobjects.size();
   // Fill the histos if pass the trigger (just the two with higher pt)
   unsigned int jaux = 0;
   // jaux is being used as a dedicated counter to avoid getting
   // a non-existent element inside dxys
   // more information in the issue https://github.com/cms-sw/cmssw/issues/32550
   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()) {
       ++jaux;
       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.empty() &&
         (objType == EVTColContainer::ELEC || objType == EVTColContainer::MUON || objType == EVTColContainer::MUTRK)) {
       this->fillHist(isPassTrigger, source, objTypeStr, "Dxy", dxys[jaux]);
       ++jaux;
     }
 
     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 && !_drop_pt2) {
       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 && !_drop_pt3) {
       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 == totalobjectssize) {
         break;
       }
     }
 
   }  // end loop over matches
 }

◆ beginJob()

void HLTExoticaPlotter::beginJob ( void )

Definition at line 42 of file HLTExoticaPlotter.cc.

42 {}

◆ beginRun()

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

◆ bookHist()

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, dqm::implementation::IBooker::book1D(), submitPVResolutionJobs::desc, SelectiveReadoutTask_cfi::edges, h, mps_fire::i, createfilelist::int, LogDebug, SiStripPI::max, SiStripPI::min, Skims_PA_cff::name, seedmultiplicitymonitor_newtracking_cfi::nBins, electrons_cff::objType, submitPVValidationJobs::params, source, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and HPSPFTaus_cff::variable.

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 = nullptr;
 
   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";
     if (variable == "MaxPt1") {
       desc = "Leading";
     } else if (variable == "MaxPt2") {
       desc = "Next-to-Leading";
     } else {
       desc = "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();
 
   if (source == "gen") {
     if (objType != "refittedStandAloneMuons") {
       _elements[name] = iBooker.book1D(name, h);
     }
   } else {
     _elements[name] = iBooker.book1D(name, h);
   }
 
   //    LogDebug("ExoticaValidation") << "                        booked histo
   //    with name " << name << "\n"
   //                  << "                        at location " <<
   //(unsigned long int)_elements[name];
   delete h;
 }

◆ fillHist()

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 266 of file HLTExoticaPlotter.cc.

References _elements, _hltPath, LogDebug, Skims_PA_cff::name, electrons_cff::objType, source, AlCaHLTBitMon_QueryRunRegistry::string, relativeConstraints::value, and HPSPFTaus_cff::variable.

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;
 
   if (source == "gen") {
     if (objType != "refittedStandAloneMuons") {
       _elements[name]->Fill(value);
     }
   } else {
     _elements[name]->Fill(value);
   }
 
   LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::fillHist()" << name << " worked";
 }

◆ gethltpath()

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

inline

Definition at line 60 of file HLTExoticaPlotter.h.

References _hltPath.

60 { return _hltPath; }

HLTExoticaPlotter::_hltPath

std::string _hltPath

Definition: HLTExoticaPlotter.h:73

◆ plotterBookHistos()

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

Definition at line 44 of file HLTExoticaPlotter.cc.

References _drop_pt2, _drop_pt3, _objectsType, bookHist(), EVTColContainer::ELEC, EVTColContainer::getTypeString(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogDebug, EVTColContainer::MUON, EVTColContainer::MUTRK, source, CalibrationSummaryClient_cfi::sources, 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");
           if (!_drop_pt2)
             bookHist(iBooker, source, objTypeStr, "MaxPt2");
           if (!_drop_pt3)
             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 || *it == EVTColContainer::MUTRK) {
             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");
           if (!_drop_pt2)
             bookHist(iBooker, source, objTypeStr, "MaxPt2");
           if (!_drop_pt3)
             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 || *it == EVTColContainer::MUTRK) {
             bookHist(iBooker, source, objTypeStr, "Dxy");
           }
         }
       }
     }
   }
 }

Member Data Documentation

◆ _drop_pt2

bool HLTExoticaPlotter::_drop_pt2

private

Definition at line 87 of file HLTExoticaPlotter.h.

Referenced by analyze(), HLTExoticaPlotter(), and plotterBookHistos().

◆ _drop_pt3

bool HLTExoticaPlotter::_drop_pt3

private

Definition at line 88 of file HLTExoticaPlotter.h.

Referenced by analyze(), HLTExoticaPlotter(), and plotterBookHistos().

◆ _elements

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

private

Definition at line 90 of file HLTExoticaPlotter.h.

Referenced by bookHist(), and fillHist().

◆ _hltPath

std::string HLTExoticaPlotter::_hltPath

private

Definition at line 73 of file HLTExoticaPlotter.h.

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

◆ _hltProcessName

std::string HLTExoticaPlotter::_hltProcessName

private

Definition at line 74 of file HLTExoticaPlotter.h.

◆ _nObjects

unsigned int HLTExoticaPlotter::_nObjects

private

Definition at line 78 of file HLTExoticaPlotter.h.

◆ _objectsType

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

private

Definition at line 76 of file HLTExoticaPlotter.h.

Referenced by analyze(), and plotterBookHistos().

◆ _parametersDxy

std::vector<double> HLTExoticaPlotter::_parametersDxy

private

Definition at line 84 of file HLTExoticaPlotter.h.

Referenced by bookHist().

◆ _parametersEta

std::vector<double> HLTExoticaPlotter::_parametersEta

private

Definition at line 80 of file HLTExoticaPlotter.h.

Referenced by bookHist().

◆ _parametersPhi

std::vector<double> HLTExoticaPlotter::_parametersPhi

private

Definition at line 81 of file HLTExoticaPlotter.h.

Referenced by bookHist().

◆ _parametersTurnOn

std::vector<double> HLTExoticaPlotter::_parametersTurnOn

private

Definition at line 82 of file HLTExoticaPlotter.h.

Referenced by bookHist().

◆ _parametersTurnOnSumEt

std::vector<double> HLTExoticaPlotter::_parametersTurnOnSumEt

private

Definition at line 83 of file HLTExoticaPlotter.h.

Referenced by bookHist().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ HLTExoticaPlotter()

◆ ~HLTExoticaPlotter()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ bookHist()

◆ fillHist()

◆ gethltpath()

◆ plotterBookHistos()

Member Data Documentation

◆ _drop_pt2

◆ _drop_pt3

◆ _elements

◆ _hltPath

◆ _hltProcessName

◆ _nObjects

◆ _objectsType

◆ _parametersDxy

◆ _parametersEta

◆ _parametersPhi

◆ _parametersTurnOn

◆ _parametersTurnOnSumEt