HLTHiggsPlotter Class Reference

#include <HLTHiggsPlotter.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< MatchStruct > &matches, const unsigned int &minCandidates)
void	analyze (const bool &isPassTrigger, const std::string &source, const std::vector< MatchStruct > &matches, std::map< std::string, bool > &nMinOne, const float &dEtaqq, const float &mqq, const float &dPhibb, const float &CSV1, const float &CSV2, const float &CSV3, const bool &passAllCuts)
void	beginJob ()
void	beginRun (const edm::Run &, const edm::EventSetup &)
void	bookHistograms (DQMStore::IBooker &, const bool &)
const std::string	gethltpath () const
	HLTHiggsPlotter (const edm::ParameterSet &pset, const std::string &hltPath, const std::vector< unsigned int > &objectsType, const unsigned int &minCandidates, const std::vector< double > &NminOneCuts)
	~HLTHiggsPlotter ()

Private Member Functions
void	bookHist (const std::string &source, const std::string &objType, const std::string &variable, DQMStore::IBooker &)
void	fillHist (const bool &passTrigger, const std::string &source, const std::string &objType, const std::string &var, const float &value)

Private Attributes
std::map< unsigned int, double >	_cutMaxEta
std::map< unsigned int, double >	_cutMinPt
std::map< unsigned int, unsigned int >	_cutMotherId
std::map< unsigned int, std::vector< double > >	_cutsDr
std::map< std::string, MonitorElement * >	_elements
std::string	_hltPath
std::string	_hltProcessName
std::vector< double >	_NminOneCuts
unsigned int	_nObjects
unsigned int	_NptPlots
std::set< unsigned int >	_objectsType
std::vector< double >	_parametersEta
std::vector< double >	_parametersPhi
std::vector< double >	_parametersTurnOn

Detailed Description

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

Author

J. Duarte Campderros (based and adapted on J. Klukas, M. Vander Donckt and J. Alcaraz code from the HLTriggerOffline/Muon package)

J. Klukas, M. Vander Donckt, J. Alcaraz

Definition at line 46 of file HLTHiggsPlotter.h.

Member Typedef Documentation

DQMStore

typedef dqm::legacy::DQMStore HLTHiggsPlotter::DQMStore

Definition at line 48 of file HLTHiggsPlotter.h.

MonitorElement

typedef dqm::legacy::MonitorElement HLTHiggsPlotter::MonitorElement

Definition at line 49 of file HLTHiggsPlotter.h.

Constructor & Destructor Documentation

HLTHiggsPlotter()

HLTHiggsPlotter::HLTHiggsPlotter	(	const edm::ParameterSet &	pset,
		const std::string &	hltPath,
		const std::vector< unsigned int > &	objectsType,
		const unsigned int &	minCandidates,
		const std::vector< double > &	NminOneCuts
	)

Definition at line 23 of file HLTHiggsPlotter.cc.

    : _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")),
      _NptPlots(NptPlots),
      _NminOneCuts(NminOneCuts) {
  for (std::set<unsigned int>::iterator it = _objectsType.begin(); it != _objectsType.end(); ++it) {
    // Some parameters extracted from the .py
    std::string objStr = EVTColContainer::getTypeString(*it);
    _cutMinPt[*it] = pset.getParameter<double>(std::string(objStr + "_cutMinPt").c_str());
    _cutMaxEta[*it] = pset.getParameter<double>(std::string(objStr + "_cutMaxEta").c_str());
  }
}

References _cutMaxEta, _cutMinPt, _objectsType, EVTColContainer::getTypeString(), muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

~HLTHiggsPlotter()

HLTHiggsPlotter::~HLTHiggsPlotter

(

)

Definition at line 45 of file HLTHiggsPlotter.cc.

{}

Member Function Documentation

analyze() [1/2]

void HLTHiggsPlotter::analyze	(	const bool &	isPassTrigger,
		const std::string &	source,
		const std::vector< MatchStruct > &	matches,
		const unsigned int &	minCandidates
	)

Definition at line 98 of file HLTHiggsPlotter.cc.

                                                                 {
  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;
  const int totalobjectssize2 = _NptPlots * 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].objType) == _objectsType.end()) {
      continue;
    }
 
    const unsigned int objType = matches[j].objType;
    const std::string objTypeStr = EVTColContainer::getTypeString(matches[j].objType);
 
    float pt = matches[j].pt;
    float eta = matches[j].eta;
    float phi = matches[j].phi;
 
    TString maxPt;
    if ((unsigned)countobjects[objType] < _NptPlots) {
      maxPt = "MaxPt";
      maxPt += (countobjects[objType] + 1);
      this->fillHist(isPassTrigger, source, objTypeStr, maxPt.Data(), pt);
      // Filled the high pt ...
      ++(countobjects[objType]);
      ++counttotal;
    } else {
      if ((unsigned)countobjects[objType] < minCandidates) {  // To get correct results for HZZ
        ++(countobjects[objType]);
        ++counttotal;
      } else
        continue;  //   Otherwise too many entries in Eta and Phi distributions
    }
 
    this->fillHist(isPassTrigger, source, objTypeStr, "Eta", eta);
    this->fillHist(isPassTrigger, source, objTypeStr, "Phi", phi);
 
    if (counttotal == totalobjectssize2) {
      break;
    }
  }
}

References _NptPlots, _objectsType, cms::cuda::co, PVValHelper::eta, fillHist(), EVTColContainer::getTypeString(), dqmiolumiharvest::j, oniaPATMuonsWithTrigger_cff::matches, L1TrackerEtMissProducer_cfi::maxPt, hltExoticaCaloHT_cff::minCandidates, electrons_cff::objType, phi, DiDispStaMuonMonitor_cfi::pt, source, and AlCaHLTBitMon_QueryRunRegistry::string.

analyze() [2/2]

void HLTHiggsPlotter::analyze	(	const bool &	isPassTrigger,
		const std::string &	source,
		const std::vector< MatchStruct > &	matches,
		std::map< std::string, bool > &	nMinOne,
		const float &	dEtaqq,
		const float &	mqq,
		const float &	dPhibb,
		const float &	CSV1,
		const float &	CSV2,
		const float &	CSV3,
		const bool &	passAllCuts
	)

Definition at line 152 of file HLTHiggsPlotter.cc.

                                                       {
  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) {
    if (!(*co == EVTColContainer::PFJET && source == "gen"))  // genJets are not there
      countobjects[*co] = 0;
  }
 
  int counttotal = 0;
  const int totalobjectssize2 = _NptPlots * 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].objType) == _objectsType.end()) {
      continue;
    }
 
    const unsigned int objType = matches[j].objType;
    const std::string objTypeStr = EVTColContainer::getTypeString(matches[j].objType);
 
    float pt = matches[j].pt;
    float eta = matches[j].eta;
    float phi = matches[j].phi;
 
    // PFMET N-1 cut
    if (objType == EVTColContainer::PFMET && _NminOneCuts[8] && !nMinOne["PFMET"])
      continue;
 
    TString maxPt;
    if ((unsigned)(countobjects)[objType] < _NptPlots) {
      maxPt = "MaxPt";
      maxPt += (countobjects[objType] + 1);
      if (objType != EVTColContainer::PFJET || nMinOne[maxPt.Data()]) {
        this->fillHist(isPassTrigger, source, objTypeStr, maxPt.Data(), pt);
      }
      ++(countobjects[objType]);
      ++counttotal;
    } else
      continue;  // if not needed (minCandidates == _NptPlots if _useNminOneCuts
    if (objType != EVTColContainer::PFJET || passAllCuts) {
      this->fillHist(isPassTrigger, source, objTypeStr, "Eta", eta);
      this->fillHist(isPassTrigger, source, objTypeStr, "Phi", phi);
    }
 
    if (counttotal == totalobjectssize2) {
      break;
    }
  }
  if (source == "rec" && _objectsType.find(EVTColContainer::PFJET) != _objectsType.end()) {
    if (_NminOneCuts[0] && nMinOne["dEtaqq"]) {
      this->fillHist(isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), "dEtaqq", dEtaqq);
    }
    if (_NminOneCuts[1] && nMinOne["mqq"]) {
      this->fillHist(isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), "mqq", mqq);
    }
    if (_NminOneCuts[2] && nMinOne["dPhibb"]) {
      this->fillHist(isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), "dPhibb", dPhibb);
    }
    if (_NminOneCuts[3]) {
      std::string nameCSVplot = "CSV1";
      if (_NminOneCuts[6])
        nameCSVplot = "maxCSV";
      if (nMinOne[nameCSVplot])
        this->fillHist(
            isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), nameCSVplot, CSV1);
    }
    if (_NminOneCuts[4] && nMinOne["CSV2"]) {
      this->fillHist(isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), "CSV2", CSV2);
    }
    if (_NminOneCuts[5] && nMinOne["CSV3"]) {
      this->fillHist(isPassTrigger, source, EVTColContainer::getTypeString(EVTColContainer::PFJET), "CSV3", CSV3);
    }
  }
}

References _NminOneCuts, _NptPlots, _objectsType, cms::cuda::co, PVValHelper::eta, fillHist(), EVTColContainer::getTypeString(), dqmiolumiharvest::j, oniaPATMuonsWithTrigger_cff::matches, L1TrackerEtMissProducer_cfi::maxPt, electrons_cff::objType, EVTColContainer::PFJET, EVTColContainer::PFMET, phi, DiDispStaMuonMonitor_cfi::pt, source, and AlCaHLTBitMon_QueryRunRegistry::string.

beginJob()

void HLTHiggsPlotter::beginJob

(

void

)

Definition at line 47 of file HLTHiggsPlotter.cc.

{}

beginRun()

void HLTHiggsPlotter::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

Definition at line 49 of file HLTHiggsPlotter.cc.

{}

bookHist()

void HLTHiggsPlotter::bookHist ( const std::string &  source, const std::string &  objType, const std::string &  variable, DQMStore::IBooker &  ibooker  )

private

Definition at line 240 of file HLTHiggsPlotter.cc.

                                                         {
  std::string sourceUpper = source;
  sourceUpper[0] = std::toupper(sourceUpper[0]);
  std::string name = source + objType + variable + "_" + _hltPath;
  TH1F *h = nullptr;
 
  if (variable.find("MaxPt") != std::string::npos) {
    std::string desc;
    if (variable == "MaxPt1")
      desc = "Leading";
    else if (variable == "MaxPt2")
      desc = "Next-to-Leading";
    else
      desc = variable.substr(5, 6) + "th Leading";
    std::string title = "pT of " + desc + " " + sourceUpper + " " + objType +
                        " "
                        "where event pass the " +
                        _hltPath;
    const size_t nBinsStandard = _parametersTurnOn.size() - 1;
    size_t nBins = nBinsStandard;
    float *edges = new float[nBinsStandard + 1];
    for (size_t i = 0; i < nBinsStandard + 1; i++) {
      edges[i] = _parametersTurnOn[i];
    }
 
    std::string jetObj = EVTColContainer::getTypeString(EVTColContainer::PFJET);
    if (objType == jetObj) {
      const size_t nBinsJets = 25;
      nBins = nBinsJets;
      delete[] edges;
      edges = new float[nBinsJets + 1];
      for (size_t i = 0; i < nBinsJets + 1; i++) {
        edges[i] = i * 10;
      }
    }
    if (objType == EVTColContainer::getTypeString(EVTColContainer::PFMET)) {
      const size_t nBinsJets = 30;
      nBins = nBinsJets;
      delete[] edges;
      edges = new float[nBinsJets + 1];
      for (size_t i = 0; i < nBinsJets + 1; i++) {
        edges[i] = i * 10;
      }
    }
    h = new TH1F(name.c_str(), title.c_str(), nBins, edges);
    delete[] edges;
  } else {
    if (variable == "dEtaqq") {
      std::string title = "#Delta #eta_{qq} of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 4.8;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "mqq") {
      std::string title = "m_{qq} of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 1000;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "dPhibb") {
      std::string title = "#Delta #phi_{bb} of " + sourceUpper + " " + objType;
      int nBins = 10;
      double min = 0;
      double max = 3.1416;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "CSV1") {
      std::string title = "CSV1 of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 1;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "CSV2") {
      std::string title = "CSV2 of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 1;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "CSV3") {
      std::string title = "CSV3 of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 1;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } else if (variable == "maxCSV") {
      std::string title = "max CSV of " + sourceUpper + " " + objType;
      int nBins = 20;
      double min = 0;
      double max = 1;
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    } 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.c_str(), h);
  delete h;
}

References _elements, _hltPath, _parametersEta, _parametersPhi, _parametersTurnOn, dqm::implementation::IBooker::book1D(), submitPVResolutionJobs::desc, SelectiveReadoutTask_cfi::edges, EVTColContainer::getTypeString(), h, mps_fire::i, createfilelist::int, SiStripPI::max, min(), Skims_PA_cff::name, seedmultiplicitymonitor_newtracking_cfi::nBins, electrons_cff::objType, CalibrationSummaryClient_cfi::params, EVTColContainer::PFJET, EVTColContainer::PFMET, source, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and taus_updatedMVAIds_cff::variable.

Referenced by bookHistograms().

bookHistograms()

void HLTHiggsPlotter::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const bool &	useNminOneCuts
	)

Definition at line 51 of file HLTHiggsPlotter.cc.

                                                                                         {
  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";
    TString maxPt;
 
    const std::string objTypeStr = EVTColContainer::getTypeString(*it);
 
    for (size_t i = 0; i < sources.size(); i++) {
      std::string source = sources[i];
 
      if (useNminOneCuts && *it == EVTColContainer::PFJET) {
        if (source == "gen")
          continue;
        else {
          // N-1 jet plots (dEtaqq, mqq, dPhibb, CSV1, maxCSV_jets, maxCSV_E, PFMET, pt1, pt2, pt3, pt4)
          if (_NminOneCuts[0])
            bookHist(source, objTypeStr, "dEtaqq", ibooker);
          if (_NminOneCuts[1])
            bookHist(source, objTypeStr, "mqq", ibooker);
          if (_NminOneCuts[2])
            bookHist(source, objTypeStr, "dPhibb", ibooker);
          if (_NminOneCuts[3]) {
            if (_NminOneCuts[6])
              bookHist(source, objTypeStr, "maxCSV", ibooker);
            else
              bookHist(source, objTypeStr, "CSV1", ibooker);
          }
          if (_NminOneCuts[4])
            bookHist(source, objTypeStr, "CSV2", ibooker);
          if (_NminOneCuts[5])
            bookHist(source, objTypeStr, "CSV3", ibooker);
        }
      }
 
      bookHist(source, objTypeStr, "Eta", ibooker);
      bookHist(source, objTypeStr, "Phi", ibooker);
      for (unsigned int i = 0; i < _NptPlots; i++) {
        maxPt = "MaxPt";
        maxPt += i + 1;
        bookHist(source, objTypeStr, maxPt.Data(), ibooker);
      }
    }
  }
}

References _NminOneCuts, _NptPlots, _objectsType, bookHist(), EVTColContainer::getTypeString(), mps_fire::i, L1TrackerEtMissProducer_cfi::maxPt, EVTColContainer::PFJET, source, CalibrationSummaryClient_cfi::sources, and AlCaHLTBitMon_QueryRunRegistry::string.

fillHist()

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

private

Definition at line 348 of file HLTHiggsPlotter.cc.

                                                   {
  std::string sourceUpper = source;
  sourceUpper[0] = toupper(sourceUpper[0]);
  std::string name = source + objType + variable + "_" + _hltPath;
 
  _elements[name]->Fill(value);
}

References _elements, _hltPath, Skims_PA_cff::name, electrons_cff::objType, source, AlCaHLTBitMon_QueryRunRegistry::string, and taus_updatedMVAIds_cff::variable.

Referenced by analyze().

gethltpath()

const std::string HLTHiggsPlotter::gethltpath ( ) const

inline

Definition at line 77 of file HLTHiggsPlotter.h.

{ return _hltPath; }

References _hltPath.

Member Data Documentation

_cutMaxEta

std::map<unsigned int, double> HLTHiggsPlotter::_cutMaxEta

private

Definition at line 100 of file HLTHiggsPlotter.h.

Referenced by HLTHiggsPlotter().

_cutMinPt

std::map<unsigned int, double> HLTHiggsPlotter::_cutMinPt

private

Definition at line 99 of file HLTHiggsPlotter.h.

Referenced by HLTHiggsPlotter().

_cutMotherId

std::map<unsigned int, unsigned int> HLTHiggsPlotter::_cutMotherId

private

Definition at line 101 of file HLTHiggsPlotter.h.

_cutsDr

std::map<unsigned int, std::vector<double> > HLTHiggsPlotter::_cutsDr

private

Definition at line 102 of file HLTHiggsPlotter.h.

_elements

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

private

Definition at line 110 of file HLTHiggsPlotter.h.

Referenced by bookHist(), and fillHist().

_hltPath

std::string HLTHiggsPlotter::_hltPath

private

Definition at line 87 of file HLTHiggsPlotter.h.

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

_hltProcessName

std::string HLTHiggsPlotter::_hltProcessName

private

Definition at line 89 of file HLTHiggsPlotter.h.

_NminOneCuts

std::vector<double> HLTHiggsPlotter::_NminOneCuts

private

Definition at line 108 of file HLTHiggsPlotter.h.

Referenced by analyze(), and bookHistograms().

_nObjects

unsigned int HLTHiggsPlotter::_nObjects

private

Definition at line 93 of file HLTHiggsPlotter.h.

_NptPlots

unsigned int HLTHiggsPlotter::_NptPlots

private

Definition at line 105 of file HLTHiggsPlotter.h.

Referenced by analyze(), and bookHistograms().

_objectsType

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

private

Definition at line 91 of file HLTHiggsPlotter.h.

Referenced by analyze(), bookHistograms(), and HLTHiggsPlotter().

_parametersEta

std::vector<double> HLTHiggsPlotter::_parametersEta

private

Definition at line 95 of file HLTHiggsPlotter.h.

Referenced by bookHist().

_parametersPhi

std::vector<double> HLTHiggsPlotter::_parametersPhi

private

Definition at line 96 of file HLTHiggsPlotter.h.

Referenced by bookHist().

_parametersTurnOn

std::vector<double> HLTHiggsPlotter::_parametersTurnOn

private

Definition at line 97 of file HLTHiggsPlotter.h.

Referenced by bookHist().