HLTExoticaPlotter Class Reference

#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)
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 >	_parametersEta
std::vector< double >	_parametersPhi
std::vector< double >	_parametersTurnOn

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 19 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"))
{
    LogDebug("ExoticaValidation") << "In HLTExoticaPlotter::constructor()";
}

HLTExoticaPlotter::~HLTExoticaPlotter

(

)

Definition at line 33 of file HLTExoticaPlotter.cc.

{
}

Member Function Documentation

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

Definition at line 88 of file HLTExoticaPlotter.cc.

References _objectsType, eta(), fillHist(), EVTColContainer::getTypeString(), j, LogDebug, 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;
    const int totalobjectssize2 = 2 * 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 (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 (counttotal == totalobjectssize2) {
        break;
      }
        }

    } // end loop over matches
}

void HLTExoticaPlotter::beginJob

(

void

)

Definition at line 38 of file HLTExoticaPlotter.cc.

{
}

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

References _elements, _hltPath, _parametersEta, _parametersPhi, _parametersTurnOn, 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 = _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 if (variable.find("MaxPt") != std::string::npos) {
        std::string desc = (variable == "MaxPt1") ? "Leading" : "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 211 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 56 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 44 of file HLTExoticaPlotter.cc.

References _objectsType, bookHist(), EVTColContainer::getTypeString(), i, LogDebug, 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, "Eta");
                bookHist(iBooker, source, objTypeStr, "Phi");
              }
            } 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, "Eta");
                bookHist(iBooker, source, objTypeStr, "Phi");
              }
            }       

        }
    }
}

Member Data Documentation

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

private

Definition at line 78 of file HLTExoticaPlotter.h.

Referenced by bookHist(), and fillHist().

std::string HLTExoticaPlotter::_hltPath

private

Definition at line 67 of file HLTExoticaPlotter.h.

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

std::string HLTExoticaPlotter::_hltProcessName

private

Definition at line 68 of file HLTExoticaPlotter.h.

unsigned int HLTExoticaPlotter::_nObjects

private

Definition at line 72 of file HLTExoticaPlotter.h.

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

private

Definition at line 70 of file HLTExoticaPlotter.h.

Referenced by analyze(), and plotterBookHistos().

std::vector<double> HLTExoticaPlotter::_parametersEta

private

Definition at line 74 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersPhi

private

Definition at line 75 of file HLTExoticaPlotter.h.

Referenced by bookHist().

std::vector<double> HLTExoticaPlotter::_parametersTurnOn

private

Definition at line 76 of file HLTExoticaPlotter.h.

Referenced by bookHist().