HLTExoticaSubAnalysis Class Reference

#include <HLTExoticaSubAnalysis.h>

Public Types
enum	Level { Level::GEN = 98, Level::RECO = 99 }

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iEventSetup, EVTColContainer *cols)
void	beginJob ()
void	beginRun (const edm::Run &iRun, const edm::EventSetup &iEventSetup)
void	endRun ()
	HLTExoticaSubAnalysis (const edm::ParameterSet &pset, const std::string &analysisname, edm::ConsumesCollector &&consCollector)
	Constructor. More...
void	subAnalysisBookHistos (DQMStore::IBooker &iBooker, const edm::Run &iRun, const edm::EventSetup &iSetup)
	~HLTExoticaSubAnalysis ()
	End Constructor. More...

Private Member Functions
void	bookHist (DQMStore::IBooker &iBooker, const std::string &source, const std::string &objType, const std::string &variable)
	The internal functions to book and fill histograms. More...
void	fillHist (const std::string &source, const std::string &objType, const std::string &variable, const float &value)
void	getHandlesToObjects (const edm::Event &iEvent, EVTColContainer *col)
	Gets the collections themselves. More...
void	getNamesOfObjects (const edm::ParameterSet &anpset)
	Creates the maps that map which collection should come from which label. More...
const std::vector< unsigned int >	getObjectsType (const std::string &hltpath) const
	closes analyze method More...
void	initSelector (const unsigned int &objtype)
	Initializes the selectors of the objects based on which object it is. More...
void	insertCandidates (const unsigned int &objtype, const EVTColContainer *col, std::vector< reco::LeafCandidate > *matches, std::map< int, double > &theSumEt, std::map< int, std::vector< const reco::Track * > > &trkObjs)
void	registerConsumes (edm::ConsumesCollector &consCollector)
	Registers consumption of objects. More...

Private Attributes
std::string	_analysisname
	The name of this sub-analysis. More...
edm::InputTag	_beamSpotLabel
edm::EDGetTokenT< reco::BeamSpot >	_bsToken
bool	_drop_pt2
bool	_drop_pt3
std::map< std::string, MonitorElement * >	_elements
	Structure of the MonitorElements. More...
std::map< unsigned int, std::string >	_genCut
	gen/rec objects cuts More...
std::map< unsigned int, std::string >	_genCut_leading
	gen/rec pt-leading objects cuts More...
StringCutObjectSelector < reco::GenMET > *	_genMETSelector
edm::InputTag	_genParticleLabel
edm::EDGetTokenT < reco::GenParticleCollection >	_genParticleToken
	And also the tokens to get the object collections. More...
std::map< unsigned int, StringCutObjectSelector < reco::GenParticle > * >	_genSelectorMap
HLTConfigProvider	_hltConfig
	Interface to the HLT information. More...
std::set< std::string >	_hltPaths
	The hlt paths found in the hltConfig. More...
std::vector< std::string >	_hltPathsToCheck
	The hlt paths to check for. More...
std::string	_hltProcessName
	The labels of the object collections to be used in this analysis. More...
StringCutObjectSelector < l1extra::L1EtMissParticle > *	_l1METSelector
unsigned int	_minCandidates
	The minimum number of reco/gen candidates needed by the analysis. More...
std::vector< double >	_parametersDxy
std::vector< double >	_parametersEta
	Some kinematical parameters. More...
std::vector< double >	_parametersPhi
std::vector< double >	_parametersTurnOn
std::vector< double >	_parametersTurnOnSumEt
std::vector< HLTExoticaPlotter >	_plotters
	The plotters: managers of each hlt path where the plots are done. More...
edm::ParameterSet	_pset
	Internal, working copy of the PSet passed from above. More...
StringCutObjectSelector < reco::CaloJet > *	_recCaloJetSelector
StringCutObjectSelector < reco::CaloMET > *	_recCaloMETSelector
StringCutObjectSelector < reco::CaloMET > *	_recCaloMHTSelector
std::map< unsigned int, std::string >	_recCut
std::map< unsigned int, std::string >	_recCut_leading
StringCutObjectSelector < reco::GsfElectron > *	_recElecSelector
std::map< unsigned int, edm::InputTag >	_recLabels
StringCutObjectSelector < reco::MET > *	_recMETSelector
StringCutObjectSelector < reco::Muon > *	_recMuonSelector
StringCutObjectSelector < reco::Track > *	_recMuonTrkSelector
StringCutObjectSelector < reco::PFJet > *	_recPFJetSelector
StringCutObjectSelector < reco::PFMET > *	_recPFMETSelector
StringCutObjectSelector < reco::PFMET > *	_recPFMHTSelector
StringCutObjectSelector < reco::PFTau > *	_recPFTauSelector
StringCutObjectSelector < reco::Photon > *	_recPhotonSelector
StringCutObjectSelector < reco::Track > *	_recTrackSelector
std::map< std::string, std::string >	_shortpath2long
	Relation between the short and long versions of the path. More...
std::map< unsigned int, edm::EDGetToken >	_tokens
std::map< std::string, int >	_triggerCounter
	counting HLT passed events More...
edm::InputTag	_trigResultsLabel
edm::EDGetTokenT < edm::TriggerResults >	_trigResultsToken

Detailed Description

This class is the main workhorse of the package. It makes the histograms for one given analysis, taking care of all HLT paths related to that analysis.

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 66 of file HLTExoticaSubAnalysis.h.

Member Enumeration Documentation

enum HLTExoticaSubAnalysis::Level

strong

Enumerator
GEN
RECO

Definition at line 68 of file HLTExoticaSubAnalysis.h.

                    {
        GEN = 98,
        RECO = 99
    };

Constructor & Destructor Documentation

HLTExoticaSubAnalysis::HLTExoticaSubAnalysis	(	const edm::ParameterSet &	pset,
		const std::string &	analysisname,
		edm::ConsumesCollector &&	consCollector
	)

Constructor.

Get the vector of paths to check, for this particular analysis.

Get the minimum candidates, for this particular analysis.

Definition at line 28 of file HLTExoticaSubAnalysis.cc.

References _drop_pt2, _drop_pt3, _genCut, _genCut_leading, _hltPathsToCheck, _minCandidates, _parametersDxy, _parametersEta, _parametersPhi, _parametersTurnOn, _parametersTurnOnSumEt, _pset, _recCut, _recCut_leading, _recLabels, edm::ParameterSet::exists(), edm::ParameterSet::existsAs(), getNamesOfObjects(), edm::ParameterSet::getParameter(), EVTColContainer::getTypeString(), edm::ParameterSet::getUntrackedParameter(), edm::ParameterSet::insert(), LogDebug, registerConsumes(), edm::ParameterSet::retrieve(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                :
    _pset(pset),
    _analysisname(analysisname),
    _minCandidates(0),
    _hltProcessName(pset.getParameter<std::string>("hltProcessName")),
    _genParticleLabel(pset.getParameter<std::string>("genParticleLabel")),
    _trigResultsLabel("TriggerResults", "", _hltProcessName),
    _beamSpotLabel(pset.getParameter<std::string>("beamSpotLabel")),
    _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),
    _recMuonSelector(0),
    _recMuonTrkSelector(0),
    _recTrackSelector(0),
    _recElecSelector(0),
    _recMETSelector(0),
    _recPFMETSelector(0),
    _recPFMHTSelector(0),
    _genMETSelector(0),
    _recCaloMETSelector(0),
    _recCaloMHTSelector(0),
    _l1METSelector(0),
    _recPFTauSelector(0),
    _recPhotonSelector(0),
    _recPFJetSelector(0),
    _recCaloJetSelector(0)
{

    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::constructor()";

    // Specific parameters for this analysis
    edm::ParameterSet anpset = pset.getParameter<edm::ParameterSet>(analysisname);

    // If this analysis has a particular set of binnings, use it.
    // (Taken from the analysis-specific parameter set, of course)
    // The "true" in the beginning of _pset.insert() means
    // "overwrite the parameter if need be".
    if (anpset.exists("parametersTurnOn")) {
        _parametersTurnOn = anpset.getParameter<std::vector<double> >("parametersTurnOn");
        _pset.insert(true, "parametersTurnOn", anpset.retrieve("parametersTurnOn"));
    }
    if (anpset.exists("parametersEta")) {
        _parametersEta = anpset.getParameter<std::vector<double> >("parametersEta");
        _pset.insert(true, "parametersEta", anpset.retrieve("parametersEta"));
    }
    if (anpset.exists("parametersPhi")) {
        _parametersPhi = anpset.getParameter<std::vector<double> >("parametersPhi");
        _pset.insert(true, "parametersPhi", anpset.retrieve("parametersPhi"));
    }
    if (anpset.exists("parametersDxy")) {
        _parametersDxy = anpset.getParameter<std::vector<double> >("parametersDxy");
        _pset.insert(true, "parametersDxy", anpset.retrieve("parametersDxy"));
    }
    if (anpset.exists("parametersTurnOnSumEt")) {
        _parametersTurnOnSumEt = anpset.getParameter<std::vector<double> >("parametersTurnOnSumEt");
        _pset.insert(true, "parametersTurnOnSumEt", anpset.retrieve("parametersTurnOnSumEt"));
    }
    if (anpset.exists("dropPt2")) {
        _drop_pt2 = anpset.getParameter<bool>("dropPt2");
        _pset.insert(true, "dropPt2", anpset.retrieve("dropPt2"));
    }
    if (anpset.exists("dropPt3")) {
        _drop_pt3 = anpset.getParameter<bool>("dropPt3");
        _pset.insert(true, "dropPt3", anpset.retrieve("dropPt3"));
    }

    // Get names of objects that we may want to get from the event.
    // Notice that genParticles are dealt with separately.
    this->getNamesOfObjects(anpset);
    
    // Since now we have the names, we should register the consumption
    // of objects.
    this->registerConsumes(consCollector);

    // Generic objects: Initialization of basic phase space cuts.
    for (std::map<unsigned int, edm::InputTag>::const_iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
    const std::string objStr = EVTColContainer::getTypeString(it->first);
        _genCut[it->first] = pset.getParameter<std::string>(objStr + "_genCut");
        _recCut[it->first] = pset.getParameter<std::string>(objStr + "_recCut");
        auto const genCutParam = objStr + "_genCut_leading";
        if (pset.exists(genCutParam)) {
          _genCut_leading[it->first] = pset.getParameter<std::string>(genCutParam);
        } else {
          _genCut_leading[it->first] = "pt>0"; // no cut
        }
        auto const recCutParam = objStr + "_recCut_leading";
        if (pset.exists(recCutParam)) {
          _recCut_leading[it->first] = pset.getParameter<std::string>(recCutParam);
        } else {
          _recCut_leading[it->first] = "pt>0"; // no cut
        }
    }

    //--- Updating parameters if has to be modified for this particular specific analysis
    for (std::map<unsigned int, edm::InputTag>::const_iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
    const std::string objStr = EVTColContainer::getTypeString(it->first);

        auto const genCutParam = objStr + "_genCut";
        if(anpset.existsAs<std::string>(genCutParam,false) ) {
            _genCut[it->first] = anpset.getUntrackedParameter<std::string>(genCutParam);
        }

        auto const recCutParam = objStr + "_recCut";
        if(anpset.existsAs<std::string>(recCutParam,false) ) {
           _recCut[it->first] = anpset.getUntrackedParameter<std::string>(recCutParam);
        } 

    }

    _hltPathsToCheck = anpset.getParameter<std::vector<std::string> >("hltPathsToCheck");
    _minCandidates = anpset.getParameter<unsigned int>("minCandidates");

} 

HLTExoticaSubAnalysis::~HLTExoticaSubAnalysis

(

)

End Constructor.

Definition at line 152 of file HLTExoticaSubAnalysis.cc.

References _genMETSelector, _genSelectorMap, _l1METSelector, _recCaloJetSelector, _recCaloMETSelector, _recCaloMHTSelector, _recElecSelector, _recMETSelector, _recMuonSelector, _recMuonTrkSelector, _recPFJetSelector, _recPFMETSelector, _recPFMHTSelector, _recPFTauSelector, _recPhotonSelector, _recTrackSelector, and python.multivaluedict::map().

{
    for (std::map<unsigned int, StringCutObjectSelector<reco::GenParticle>* >::iterator it = _genSelectorMap.begin();
         it != _genSelectorMap.end(); ++it) {
        delete it->second;
        it->second = 0;
    }
    delete _recMuonSelector;
    _recMuonSelector = 0;
    delete _recMuonTrkSelector;
    _recMuonTrkSelector = 0;
    delete _recTrackSelector;
    _recTrackSelector = 0;
    delete _recElecSelector;
    _recElecSelector = 0;
    delete _recPhotonSelector;
    _recPhotonSelector = 0;
    delete _recMETSelector;
    _recMETSelector = 0;
    delete _recPFMETSelector;
    _recPFMETSelector = 0;
    delete _recPFMHTSelector;
    _recPFMHTSelector = 0;
    delete _genMETSelector;
    _genMETSelector = 0;
    delete _recCaloMETSelector;
    _recCaloMETSelector = 0;
    delete _recCaloMHTSelector;
    _recCaloMHTSelector = 0;
    delete _l1METSelector;
    _l1METSelector = 0;
    delete _recPFTauSelector;
    _recPFTauSelector = 0;
    delete _recPFJetSelector;
    _recPFJetSelector = 0;
    delete _recCaloJetSelector;
    _recCaloJetSelector = 0;
}

Member Function Documentation

void HLTExoticaSubAnalysis::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iEventSetup,
		EVTColContainer *	cols
	)

Method to fill all relevant histograms. Notice that we update the EVTColContaner to point to the collections we want.

We are going to make a fake reco::LeafCandidate, with our particleType as the pdgId. This is an alternative to the older implementation with MatchStruct.

Filling the histograms if pass the minimum amount of candidates needed by the analysis:

GEN CASE ///

Close GEN case

RECO CASE ///

Debugging.

Close RECO case

Definition at line 383 of file HLTExoticaSubAnalysis.cc.

References _drop_pt2, _drop_pt3, _genCut, _genCut_leading, _genSelectorMap, _minCandidates, _plotters, _recCut_leading, _recLabels, _shortpath2long, _triggerCounter, edm::HLTGlobalStatus::accept(), EVTColContainer::bs, funct::cos(), EVTColContainer::ELEC, eta, fillHist(), EVTColContainer::genParticles, getHandlesToObjects(), EVTColContainer::getTypeString(), i, initSelector(), insertCandidates(), j, findQualityFiles::jj, LogDebug, EVTColContainer::MUON, reco::Candidate::p4(), phi, reco::BeamSpot::position(), EnergyCorrector::pt, edm::second(), benchmark_cfg::select, funct::sin(), findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), EVTColContainer::triggerResults, trigNames, findQualityFiles::v, reco::Candidate::vertex(), reco::BeamSpot::x0(), and reco::BeamSpot::y0().

{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::analyze()";

    // Loop over _recLabels to make sure everything is alright.
    /*
    std::cout << "Now printing the _recLabels" << std::endl;
    for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
    std::cout <<  "Number: " << it->first << "\t" << "Label: " << it->second.label() << std::endl;  
    }
    */
    
    // Initialize the collection (the ones which have not been initialiazed yet)
    //std::cout << "Setting handles to objects..." << std::endl;
    this->getHandlesToObjects(iEvent, cols);

    // Utility map, mapping kinds of objects (GEN, RECO) to strings ("gen","rec")
    //std::map<Level, std::string> u2str;
    //u2str[Level::GEN] = "gen";
    //u2str[Level::RECO] = "rec";

    // Extract the match structure containing the gen/reco candidates (electron, muons,...). This part is common to all the SubAnalyses
    std::vector<reco::LeafCandidate> matchesGen; matchesGen.clear();
    std::vector<reco::LeafCandidate> matchesReco; matchesReco.clear();
    std::map<int , double> theSumEt; // map< pdgId ; SumEt > in order to keep track of the MET type
    std::map<int, std::vector<const reco::Track*> > trkObjs;

    // --- deal with GEN objects first.
    // Make each good GEN object into the base cand for a MatchStruct
    // Our definition of "good" is "passes the selector" defined in the config.py
    // Save all the MatchStructs in the "matchesGen" vector.

    for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
        // Here we are filling the vector of StringCutObjectSelector<reco::GenParticle>
        // with objects constructed from the strings saved in _genCut.
        // Initialize selectors when first event

    //std::cout << "Loop over the kinds of objects: objects of kind " << it->first << std::endl;
    

    if (!_genSelectorMap[it->first]) {
            _genSelectorMap[it->first] = new StringCutObjectSelector<reco::GenParticle>(_genCut[it->first]);
        }

        const std::string objTypeStr = EVTColContainer::getTypeString(it->first);
        // genAnyMET doesn't make sense. No need their matchesGens
        if ( TString(objTypeStr).Contains("MET") || 
             TString(objTypeStr).Contains("MHT") ||
             TString(objTypeStr).Contains("Jet")   ) continue;

        // Now loop over the genParticles, and apply the operator() over each of them.
        // Fancy syntax: for objects X and Y, X.operator()(Y) is the same as X(Y).
        for (size_t i = 0; i < cols->genParticles->size(); ++i) {
        //std::cout << "Now matchesGen.size() is " << matchesGen.size() << std::endl;
            if (_genSelectorMap[it->first]->operator()(cols->genParticles->at(i))) {
                const reco::Candidate* cand = &(cols->genParticles->at(i));
        //std::cout << "Found good cand: cand->pt() = " << cand->pt() << std::endl;
        //matchesGen.push_back(MatchStruct(cand, it->first));
        reco::LeafCandidate v(0,cand->p4(),cand->vertex(),it->first,0,true);

                matchesGen.push_back(v);
            }
        }
    }

    // Sort the matches by pT for later filling of turn-on curve
    //std::cout << "Before sorting: matchesGen.size() = " << matchesGen.size() << std::endl;

    // GreaterByPt<reco::LeafCandidate> comparator;
    // std::sort(matchesGen.begin(), 
    //        matchesGen.end(), 
    //        comparator);

    // --- same for RECO objects
    // Extraction of the objects candidates
    if(verbose>0) LogDebug("ExoticaValidation") << "-- enter loop over recLabels" ;
    for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
      //std::cout << "Filling RECO \"matchesReco\" vector for particle kind it->first = "
      //      << it->first << ", which means " << it->second.label() << std::endl;
      // Reco selectors (the function takes into account if it was instantiated
      // before or not) ### Thiago ---> Then why don't we put it in the beginRun???
      this->initSelector(it->first);
      // -- Storing the matchesReco
      this->insertCandidates(it->first, cols, &matchesReco, theSumEt, trkObjs);
      if(verbose>0) LogDebug("ExoticaValidation") << "--- " << EVTColContainer::getTypeString(it->first) 
                  << " sumEt=" << theSumEt[it->first] ;
    }

    // std::sort(matchesReco.begin(), 
    //        matchesReco.end(), 
    //        comparator);

    // -- Trigger Results
    const edm::TriggerNames trigNames = iEvent.triggerNames(*(cols->triggerResults));

    // counting HLT passed events for debugging
    for (std::vector<HLTExoticaPlotter>::iterator an = _plotters.begin(); an != _plotters.end(); ++an) {
      const std::string hltPath = _shortpath2long[an->gethltpath()];
      const bool ispassTrigger =  cols->triggerResults->accept(trigNames.triggerIndex(hltPath));
      if (ispassTrigger) _triggerCounter.find(an->gethltpath())->second++;
    }


    //for (std::map<unsigned int, std::vector<MatchStruct> >::iterator it = sourceMatchMap.begin(); it != sourceMatchMap.end(); ++it) {
    // it->first: gen/reco   it->second: HLT matches (std::vector<MatchStruct>)

    //if (it->second.size() < _minCandidates) {  // FIXME: A bug is potentially here: what about the mixed channels?
    //continue;
    //}

    if( matchesGen.size() >= _minCandidates) {  // FIXME: A bug is potentially here: what about the mixed channels?
      // Okay, there are enough candidates. Move on!

      // Filling the gen/reco objects (eff-denominators):
      // Just the first two different ones, if there are more
      // The countobjects maps uints (object types, really) --> integers.
      // Example:
      // | uint | int |
      // |  0   |  1  | --> 1 muon used
      // |  1   |  2  | --> 2 electrons used  

      // Initializing the count of the used objects.
      std::map<unsigned int, int> countobjects;
      for (std::map<unsigned int, edm::InputTag>::iterator co = _recLabels.begin();
       co != _recLabels.end(); ++co) {
    //countobjects->insert(std::pair<unsigned int, int>(co->first, 0));
    countobjects.insert(std::pair<unsigned int, int>(co->first, 0));
      }
    
      int counttotal = 0;

      // 3 : pt1, pt2, pt3
      int totalobjectssize = 1;
      if (!_drop_pt2) totalobjectssize++;
      if (!_drop_pt3) totalobjectssize++;
      totalobjectssize *= countobjects.size();


      bool isPassedLeadingCut = true;
      // We will proceed only when cuts for the pt-leading are satisified.
      for (size_t j = 0; j != matchesGen.size(); ++j) {
    const unsigned int objType = matchesGen[j].pdgId();
        // Cut for the pt-leading object 
        StringCutObjectSelector<reco::LeafCandidate> select( _genCut_leading[objType] );
        if ( !select( matchesGen[j] ) ) { // No interest case
          isPassedLeadingCut = false;     // Will skip the following matchesGen loop
          matchesGen.clear();
          break;
        }
      }
 
      std::vector<float> dxys; dxys.clear();

      for (size_t j = 0; ( j != matchesGen.size() ) && isPassedLeadingCut; ++j) {
    const unsigned int objType = matchesGen[j].pdgId();
    //std::cout << "(4) Gonna call with " << objType << std::endl;
    const std::string objTypeStr = EVTColContainer::getTypeString(objType);
      
    float pt  = matchesGen[j].pt();

    if (countobjects[objType] == 0) {
      this->fillHist("gen", objTypeStr, "MaxPt1", pt);
      ++(countobjects[objType]);
      ++counttotal;
    } 
    else if (countobjects[objType] == 1 && !_drop_pt2) {
      this->fillHist("gen", objTypeStr, "MaxPt2", pt);
      ++(countobjects[objType]);
      ++counttotal;
    } 
    else if (countobjects[objType] == 2 && !_drop_pt3) {
      this->fillHist("gen", objTypeStr, "MaxPt3", pt);
      ++(countobjects[objType]);
      ++counttotal;
    } 
    else {
      // Already the minimum three objects has been filled, get out...
      if (counttotal == totalobjectssize) {
        size_t max_size = matchesGen.size();
        for ( size_t jj = j; jj < max_size; jj++ ) {
          matchesGen.erase(matchesGen.end());
        }
        break;
      }
    }
        
    float eta = matchesGen[j].eta();
    float phi = matchesGen[j].phi();

    this->fillHist("gen", objTypeStr, "Eta", eta);
    this->fillHist("gen", objTypeStr, "Phi", phi);

        // If the target is electron or muon,

        if ( objType == EVTColContainer::MUON ||
             objType == EVTColContainer::ELEC   ) {
          const math::XYZPoint & vtx = matchesGen[j].vertex();
          float momphi = matchesGen[j].momentum().phi();
          float dxyGen = (-(vtx.x()-cols->bs->x0())*sin(momphi)+(vtx.y()-cols->bs->y0())*cos(momphi));
          dxys.push_back(dxyGen);
          this->fillHist("gen", objTypeStr, "Dxy", dxyGen);
        }

      } // Closes loop in gen
    
      // Calling to the plotters analysis (where the evaluation of the different trigger paths are done)
      //const std::string source = "gen";
      for (std::vector<HLTExoticaPlotter>::iterator an = _plotters.begin(); an != _plotters.end(); ++an) {
    const std::string hltPath = _shortpath2long[an->gethltpath()];
    const bool ispassTrigger =  cols->triggerResults->accept(trigNames.triggerIndex(hltPath));
    LogDebug("ExoticaValidation") << "                        preparing to call the plotters analysis";
        an->analyze(ispassTrigger, "gen", matchesGen, theSumEt, dxys);
    LogDebug("ExoticaValidation") << "                        called the plotter";
      }
    } 


    {
    if(matchesReco.size() < _minCandidates) return; // FIXME: A bug is potentially here: what about the mixed channels?

    // Okay, there are enough candidates. Move on!
    
    // Filling the gen/reco objects (eff-denominators):
    // Just the first two different ones, if there are more
    // The countobjects maps uints (object types, really) --> integers.
    // Example:
    // | uint | int |
    // |  0   |  1  | --> 1 muon used
    // |  1   |  2  | --> 2 electrons used  
    // Initializing the count of the used objects.
    //std::map<unsigned int, int> * countobjects = new std::map<unsigned int, int>;
    std::map<unsigned int, int> countobjects;
    for (std::map<unsigned int, edm::InputTag>::iterator co = _recLabels.begin();
         co != _recLabels.end(); ++co) {
      countobjects.insert(std::pair<unsigned int, int>(co->first, 0));
    }
    
    int counttotal = 0;

        // 3 : pt1, pt2, pt3
        int totalobjectssize = 1;
        if (!_drop_pt2) totalobjectssize++;
        if (!_drop_pt3) totalobjectssize++;
        totalobjectssize *= countobjects.size();
    
    //std::cout << "Our RECO vector has matchesReco.size() = " << matchesReco.size() << std::endl;

        std::vector<float> dxys; dxys.clear();

        bool isPassedLeadingCut = true;
        // We will proceed only when cuts for the pt-leading are satisified.
        for (size_t j = 0; j != matchesReco.size(); ++j) {
          const unsigned int objType = matchesReco[j].pdgId();
          // Cut for the pt-leading object 
          StringCutObjectSelector<reco::LeafCandidate> select( _recCut_leading[objType] );
          if ( !select( matchesReco[j] ) ) { // No interest case 
            isPassedLeadingCut = false;      // Will skip the following matchesReco loop
            matchesReco.clear();
            break;
          }
        }

    for (size_t j = 0; ( j != matchesReco.size() ) && isPassedLeadingCut; ++j) {
        const unsigned int objType = matchesReco[j].pdgId();
        //std::cout << "(4) Gonna call with " << objType << std::endl;
        const std::string objTypeStr = EVTColContainer::getTypeString(objType);
        
        float pt  = matchesReco[j].pt();

        if (countobjects[objType] == 0) {
        this->fillHist("rec", objTypeStr, "MaxPt1", pt);
        ++(countobjects[objType]);
        ++counttotal;
        } 
        else if (countobjects[objType] == 1 && !_drop_pt2) {
          if( ! ( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) {
        this->fillHist("rec", objTypeStr, "MaxPt2", pt);
          } 
          ++(countobjects[objType]);
          ++counttotal;
        } 
        else if (countobjects[objType] == 2 && !_drop_pt3) {
          if( ! ( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) {
        this->fillHist("rec", objTypeStr, "MaxPt3", pt);
          } 
          ++(countobjects[objType]);
          ++counttotal;
        } 
        else {
          // Already the minimum three objects has been filled, get out...
          if (counttotal == totalobjectssize) {
        size_t max_size = matchesReco.size();
        for ( size_t jj = j; jj < max_size; jj++ ) {
          matchesReco.erase(matchesReco.end());
        }
        break;
          }
        }

        float eta = matchesReco[j].eta();
        float phi = matchesReco[j].phi();
    
            if ( !( TString(objTypeStr).Contains("MET") || TString(objTypeStr).Contains("MHT") ) ) { 
          this->fillHist("rec", objTypeStr, "Eta", eta);
          this->fillHist("rec", objTypeStr, "Phi", phi);
            }
        else {
          this->fillHist("rec", objTypeStr, "SumEt", theSumEt[objType]);
        }
          
            if (trkObjs[objType].size()>=j+1) {
              float dxyRec = trkObjs[objType].at(j)->dxy(cols->bs->position());
              this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
              dxys.push_back(dxyRec);
            }

    } // Closes loop in reco

    //LogDebug("ExoticaValidation") << "                        deleting countobjects";
    //delete countobjects;
    
    // Calling to the plotters analysis (where the evaluation of the different trigger paths are done)
    //const std::string source = "reco";
    for (std::vector<HLTExoticaPlotter>::iterator an = _plotters.begin(); an != _plotters.end(); ++an) {
        const std::string hltPath = _shortpath2long[an->gethltpath()];
        const bool ispassTrigger =  cols->triggerResults->accept(trigNames.triggerIndex(hltPath));
        LogDebug("ExoticaValidation") << "                        preparing to call the plotters analysis";
            an->analyze(ispassTrigger, "rec", matchesReco, theSumEt, dxys);
        LogDebug("ExoticaValidation") << "                        called the plotter";
    }
    } 

} 

void HLTExoticaSubAnalysis::beginJob

(

void

)

Definition at line 192 of file HLTExoticaSubAnalysis.cc.

{
}

void HLTExoticaSubAnalysis::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iEventSetup
	)

Construct the plotters right here. For that we need to create the _hltPaths vector.

Definition at line 276 of file HLTExoticaSubAnalysis.cc.

References _analysisname, _hltConfig, _hltPaths, _hltPathsToCheck, _hltProcessName, _plotters, _pset, _recLabels, _shortpath2long, _triggerCounter, newFWLiteAna::found, i, HLTConfigProvider::init(), j, LogDebug, LogTrace, fed_dqm_sourceclient-live_cfg::path, listBenchmarks::pattern, AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerNames().

{

    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::beginRun()";
    

    // Initialize the HLT config.
    bool changedConfig(true);
    if (!_hltConfig.init(iRun, iSetup, _hltProcessName, changedConfig)) {
        edm::LogError("ExoticaValidation") << "HLTExoticaSubAnalysis::constructor(): "
                                           << "Initialization of HLTConfigProvider failed!";
    }

    // Parse the input paths to get them if they are in the table and associate
    // them to the last filter of the path (in order to extract the objects).
    _hltPaths.clear();
    for (size_t i = 0; i < _hltPathsToCheck.size(); ++i) {
        bool found = false;
        TPRegexp pattern(_hltPathsToCheck[i]);
    
    // Loop over triggerNames from _hltConfig
    for (size_t j = 0 ; j < _hltConfig.triggerNames().size(); ++j) {
            std::string thetriggername = _hltConfig.triggerNames()[j];
            if (TString(thetriggername).Contains(pattern)) {
                _hltPaths.insert(thetriggername);
                found = true;
            }
        if(verbose>2 && i==0) 
          LogDebug("ExoticaValidation") << "--- TRIGGER PATH : " << thetriggername;
        }
    
    // Oh dear, the path we wanted seems to not be available
        if (! found && verbose>2) {
            edm::LogWarning("ExoticaValidation") << "HLTExoticaSubAnalysis::constructor(): In "
                                                 << _analysisname << " subfolder NOT found the path: '"
                                                 << _hltPathsToCheck[i] << "*'" ;
        }
    } // Close loop over paths to check.
    
    // At this point, _hltpaths contains the names of the paths to check
    // that were found. Let's log it at trace level.
    LogTrace("ExoticaValidation") << "SubAnalysis: " << _analysisname 
                  << "\nHLT Trigger Paths found >>>";
    for (std::set<std::string>::const_iterator iter = _hltPaths.begin();
     iter != _hltPaths.end(); ++iter) {
        LogTrace("ExoticaValidation") << (*iter) << "\n";
    }

    // Initialize the plotters (analysers for each trigger path)
    _plotters.clear();
    for (std::set<std::string>::iterator iPath = _hltPaths.begin();
         iPath != _hltPaths.end(); ++iPath) {
        // Avoiding the dependence of the version number for the trigger paths
        std::string path = * iPath;
        std::string shortpath = path;
        if (path.rfind("_v") < path.length()) {
            shortpath = path.substr(0, path.rfind("_v"));
        }
        _shortpath2long[shortpath] = path;

        // Objects needed by the HLT path
    // Thiago: instead of trying to decode the objects from the path,
    // put the burden on the user to tell us which objects are needed.
        //const std::vector<unsigned int> objsNeedHLT = this->getObjectsType(shortpath);
    std::vector<unsigned int> objsNeedHLT;
        for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin() ;
             it != _recLabels.end(); ++it) {
            objsNeedHLT.push_back(it->first);
        }

        /*std::vector<unsigned int> userInstantiate;
        // Sanity check: the object needed by a trigger path should be
        // introduced by the user via config python (_recLabels datamember)
        for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin() ;
             it != _recLabels.end(); ++it) {
            userInstantiate.push_back(it->first);
        }
        for (std::vector<unsigned int>::const_iterator it = objsNeedHLT.begin(); it != objsNeedHLT.end();
             ++it) {
            if (std::find(userInstantiate.begin(), userInstantiate.end(), *it) ==
                userInstantiate.end()) {
                edm::LogError("ExoticaValidation") << "In HLTExoticaSubAnalysis::beginRun, "
                                                   << "Incoherence found in the python configuration file!!\nThe SubAnalysis '"
                                                   << _analysisname << "' has been asked to evaluate the trigger path '"
                                                   << shortpath << "' (found it in 'hltPathsToCheck') BUT this path"
                                                   << " needs a '" << EVTColContainer::getTypeString(*it)
                                                   << "' which has not been instantiated ('recVariableLabels'"
                                                   << ")" ;
                exit(-1); // This should probably throw an exception...
            }
        }
    */
        LogTrace("ExoticaValidation") << " --- " << shortpath;

        // The hlt path, the objects (electrons, muons, photons, ...)
        // needed to evaluate the path are the argumens of the plotter
        HLTExoticaPlotter analyzer(_pset, shortpath, objsNeedHLT);
        _plotters.push_back(analyzer);
        // counting HLT passed events for debug
        _triggerCounter.insert(std::map<std::string,int>::value_type(shortpath,0));
    }// Okay, at this point we have prepared all the plotters.

}

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

private

The internal functions to book and fill histograms.

Definition at line 1069 of file HLTExoticaSubAnalysis.cc.

References _elements, _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 subAnalysisBookHistos().

{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::bookHist()";
    std::string sourceUpper = source;
    sourceUpper[0] = std::toupper(sourceUpper[0]);
    std::string name = source + objType + variable ;
    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" : "Next-to-Leading";
      std::string title = "pT of " + desc + " " + 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 {
      std::string symbol = (variable == "Eta") ? "#eta" : "#phi";
      std::string title  = symbol + " of " + sourceUpper + " " + objType;
      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();
    // This is the trick, that takes a normal TH1F and puts it in in the DQM
    // machinery. Seems to be easy!
    // Updated to use the new iBooker machinery.
    _elements[name] = iBooker.book1D(name, h);
    delete h;
}

void HLTExoticaSubAnalysis::endRun

(

void

)

Definition at line 1337 of file HLTExoticaSubAnalysis.cc.

References _analysisname, _triggerCounter, cmsBatch::log, and LogDebug.

                                   {
    // Dump trigger results
    std::stringstream log;
    log << std::endl;
    log << "===========================================================================" << std::endl;
    log << "          Trigger Results ( " << _analysisname << " )                      " << std::endl;
    log << "===========================================================================" << std::endl;
    log << std::setw(18) << "# of passed events : HLT path names"                        << std::endl;
    log << "-------------------:-------------------------------------------------------" << std::endl;
    for (std::map<std::string,int>::iterator it = _triggerCounter.begin();
         it != _triggerCounter.end(); ++it) {
      log << std::setw(18) << it->second <<  " : " << it->first << std::endl;
    }
    log << "===========================================================================" << std::endl;
    LogDebug("ExoticaValidation") << log.str().data();
}

void HLTExoticaSubAnalysis::fillHist ( const std::string &  source, const std::string &  objType, const std::string &  variable, const float &  value  )

private

Definition at line 1130 of file HLTExoticaSubAnalysis.cc.

References _elements, 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 ;

    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::fillHist() " << name << " " << value;
    _elements[name]->Fill(value);
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::fillHist() " << name << " worked";
}

void HLTExoticaSubAnalysis::getHandlesToObjects ( const edm::Event &  iEvent, EVTColContainer *  col  )

private

Gets the collections themselves.

Definition at line 953 of file HLTExoticaSubAnalysis.cc.

References _bsToken, _genParticleToken, _tokens, _trigResultsToken, EVTColContainer::bs, EVTColContainer::CALOJET, EVTColContainer::CALOMET, EVTColContainer::CALOMHT, EVTColContainer::ELEC, EVTColContainer::GENMET, EVTColContainer::genParticles, edm::Event::getByToken(), EVTColContainer::isCommonInit(), edm::HandleBase::isValid(), EVTColContainer::L1MET, LogDebug, EVTColContainer::MET, EVTColContainer::MUON, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::PFMET, EVTColContainer::PFMHT, EVTColContainer::PFTAU, EVTColContainer::PHOTON, edm::Handle< T >::product(), EVTColContainer::set(), EVTColContainer::setCaloMHT(), EVTColContainer::setPFMHT(), EVTColContainer::TRACK, and EVTColContainer::triggerResults.

Referenced by analyze().

{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::getHandlesToObjects()";

    if (! col->isCommonInit()) {
        // Extract the trigger results (path info, pass,...)
        edm::Handle<edm::TriggerResults> trigResults;
        iEvent.getByToken(_trigResultsToken, trigResults);
        if (trigResults.isValid()) {
            col->triggerResults = trigResults.product();
        LogDebug("ExoticaValidation") << "Added handle to triggerResults";
        }

        // Extract the genParticles
        edm::Handle<reco::GenParticleCollection> genPart;
        iEvent.getByToken(_genParticleToken, genPart);
        if (genPart.isValid()) {
            col->genParticles = genPart.product();
        LogDebug("ExoticaValidation") << "Added handle to genParticles";
        }

        // BeamSpot for dxy 
        edm::Handle<reco::BeamSpot> bsHandle;
        iEvent.getByToken(_bsToken, bsHandle);
        if (bsHandle.isValid()) {
            col->bs = bsHandle.product();
        }
    }

    // Loop over the tokens and extract all other objects
    LogDebug("ExoticaValidation") << "We have got " << _tokens.size() << "tokens";
    for (std::map<unsigned int, edm::EDGetToken>::iterator it = _tokens.begin();
         it != _tokens.end(); ++it) {
        if (it->first == EVTColContainer::MUON) {
            edm::Handle<reco::MuonCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::MUTRK) {
            edm::Handle<reco::TrackCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else if (it->first == EVTColContainer::TRACK) {
            edm::Handle<reco::TrackCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else if (it->first == EVTColContainer::ELEC) {
            edm::Handle<reco::GsfElectronCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::PHOTON) {
            edm::Handle<reco::PhotonCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::MET) {
            edm::Handle<reco::METCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::PFMET) {
      edm::Handle<reco::PFMETCollection> theHandle;
      iEvent.getByToken(it->second, theHandle);
      if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::PFMHT) {
      edm::Handle<reco::PFMETCollection> theHandle;
      iEvent.getByToken(it->second, theHandle);
      if (theHandle.isValid()) col->setPFMHT(theHandle.product());
        } 
    else if (it->first == EVTColContainer::GENMET) {
            edm::Handle<reco::GenMETCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
        else if (it->first == EVTColContainer::CALOMET) {
            edm::Handle<reco::CaloMETCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else if (it->first == EVTColContainer::CALOMHT) {
      edm::Handle<reco::CaloMETCollection> theHandle;
      iEvent.getByToken(it->second, theHandle);
      if (theHandle.isValid()) col->setCaloMHT(theHandle.product());
        } 
        else if (it->first == EVTColContainer::L1MET) {
            edm::Handle<l1extra::L1EtMissParticleCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else if (it->first == EVTColContainer::PFTAU) {
            edm::Handle<reco::PFTauCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        } 
    else if (it->first == EVTColContainer::PFJET) {
            edm::Handle<reco::PFJetCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else if (it->first == EVTColContainer::CALOJET) {
            edm::Handle<reco::CaloJetCollection> theHandle;
            iEvent.getByToken(it->second, theHandle);
            if (theHandle.isValid()) col->set(theHandle.product());
        }
    else {
            edm::LogError("ExoticaValidation") << "HLTExoticaSubAnalysis::getHandlesToObjects "
                                               << " NOT IMPLEMENTED (yet) ERROR: '" << it->first << "'";
        }
    }
}

void HLTExoticaSubAnalysis::getNamesOfObjects ( const edm::ParameterSet &  anpset )

private

Creates the maps that map which collection should come from which label.

Definition at line 772 of file HLTExoticaSubAnalysis.cc.

References _analysisname, _genSelectorMap, _recLabels, EVTColContainer::CALOJET, EVTColContainer::CALOMET, EVTColContainer::CALOMHT, EVTColContainer::ELEC, edm::ParameterSet::exists(), EVTColContainer::GENMET, edm::ParameterSet::getParameter(), EVTColContainer::L1MET, LogDebug, EVTColContainer::MET, EVTColContainer::MUON, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::PFMET, EVTColContainer::PFMHT, EVTColContainer::PFTAU, EVTColContainer::PHOTON, and EVTColContainer::TRACK.

Referenced by HLTExoticaSubAnalysis().

{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::getNamesOfObjects()";

    if (anpset.exists("recMuonLabel")) {
        _recLabels[EVTColContainer::MUON] = anpset.getParameter<edm::InputTag>("recMuonLabel");
        _genSelectorMap[EVTColContainer::MUON] = 0 ;
    }
    if (anpset.exists("recMuonTrkLabel")) {
        _recLabels[EVTColContainer::MUTRK] = anpset.getParameter<edm::InputTag>("recMuonTrkLabel");
        _genSelectorMap[EVTColContainer::MUTRK] = 0 ;
    }
    if (anpset.exists("recTrackLabel")) {
        _recLabels[EVTColContainer::TRACK] = anpset.getParameter<edm::InputTag>("recTrackLabel");
        _genSelectorMap[EVTColContainer::TRACK] = 0 ;
    }
    if (anpset.exists("recElecLabel")) {
        _recLabels[EVTColContainer::ELEC] = anpset.getParameter<edm::InputTag>("recElecLabel");
        _genSelectorMap[EVTColContainer::ELEC] = 0 ;
    }
    if (anpset.exists("recPhotonLabel")) {
        _recLabels[EVTColContainer::PHOTON] = anpset.getParameter<edm::InputTag>("recPhotonLabel");
        _genSelectorMap[EVTColContainer::PHOTON] = 0 ;
    }
    if (anpset.exists("recMETLabel")) {
        _recLabels[EVTColContainer::MET] = anpset.getParameter<edm::InputTag>("recMETLabel");
        _genSelectorMap[EVTColContainer::MET] = 0 ;
    }
    if (anpset.exists("recPFMETLabel")) {
        _recLabels[EVTColContainer::PFMET] = anpset.getParameter<edm::InputTag>("recPFMETLabel");
        _genSelectorMap[EVTColContainer::PFMET] = 0 ;
    }
    if (anpset.exists("recPFMHTLabel")) {
        _recLabels[EVTColContainer::PFMHT] = anpset.getParameter<edm::InputTag>("recPFMHTLabel");
        _genSelectorMap[EVTColContainer::PFMHT] = 0 ;
    }
    if (anpset.exists("genMETLabel")) {
        _recLabels[EVTColContainer::GENMET] = anpset.getParameter<edm::InputTag>("genMETLabel");
        _genSelectorMap[EVTColContainer::GENMET] = 0 ;
    }
    if (anpset.exists("recCaloMETLabel")) {
        _recLabels[EVTColContainer::CALOMET] = anpset.getParameter<edm::InputTag>("recCaloMETLabel");
        _genSelectorMap[EVTColContainer::CALOMET] = 0 ;
    }
    if (anpset.exists("recCaloMHTLabel")) {
        _recLabels[EVTColContainer::CALOMHT] = anpset.getParameter<edm::InputTag>("recCaloMHTLabel");
        _genSelectorMap[EVTColContainer::CALOMHT] = 0 ;
    }
    if (anpset.exists("hltMETLabel")) {
        _recLabels[EVTColContainer::CALOMET] = anpset.getParameter<edm::InputTag>("hltMETLabel");
        _genSelectorMap[EVTColContainer::CALOMET] = 0 ;
    }
    if (anpset.exists("l1METLabel")) {
        _recLabels[EVTColContainer::L1MET] = anpset.getParameter<edm::InputTag>("l1METLabel");
        _genSelectorMap[EVTColContainer::L1MET] = 0 ;
    }
    if (anpset.exists("recPFTauLabel")) {
        _recLabels[EVTColContainer::PFTAU] = anpset.getParameter<edm::InputTag>("recPFTauLabel");
        _genSelectorMap[EVTColContainer::PFTAU] = 0 ;
    }
    if (anpset.exists("recPFJetLabel")) {
        _recLabels[EVTColContainer::PFJET] = anpset.getParameter<edm::InputTag>("recPFJetLabel");
        _genSelectorMap[EVTColContainer::PFJET] = 0 ;
    }
    if (anpset.exists("recCaloJetLabel")) {
        _recLabels[EVTColContainer::CALOJET] = anpset.getParameter<edm::InputTag>("recCaloJetLabel");
        _genSelectorMap[EVTColContainer::CALOJET] = 0 ;
    }

    if (_recLabels.size() < 1) {
        edm::LogError("ExoticaValidation") << "HLTExoticaSubAnalysis::getNamesOfObjects, "
                                           << "Not included any object (recMuonLabel, recElecLabel, ...)  "
                                           << "in the analysis " << _analysisname;
        return;
    }
}

const std::vector< unsigned int > HLTExoticaSubAnalysis::getObjectsType ( const std::string &  hltpath ) const

private

closes analyze method

Return the objects (muons,electrons,photons,...) needed by a HLT path. Will in general return: 0 for muon, 1 for electron, 2 for photon, 3 for PFMET, 4 for PFTau, 5 for Jet. Notice that this function is really based on a parsing of the name of the path; any incongruences there may lead to problems.

Definition at line 732 of file HLTExoticaSubAnalysis.cc.

References EVTColContainer::CALOJET, EVTColContainer::CALOMET, EVTColContainer::CALOMHT, EVTColContainer::ELEC, EVTColContainer::GENMET, EVTColContainer::getTypeString(), i, EVTColContainer::L1MET, LogDebug, EVTColContainer::MET, EVTColContainer::MUON, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::PFMET, EVTColContainer::PFMHT, EVTColContainer::PFTAU, EVTColContainer::PHOTON, AlCaHLTBitMon_QueryRunRegistry::string, and EVTColContainer::TRACK.

{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::getObjectsType()";

    static const unsigned int objSize = 15;
    static const unsigned int objtriggernames[] = {
        EVTColContainer::MUON,
        EVTColContainer::MUTRK,
        EVTColContainer::TRACK,
        EVTColContainer::ELEC,
        EVTColContainer::PHOTON,
        EVTColContainer::MET,
        EVTColContainer::PFMET,
        EVTColContainer::PFMHT,
        EVTColContainer::GENMET,
        EVTColContainer::CALOMET,
        EVTColContainer::CALOMHT,
        EVTColContainer::L1MET,
        EVTColContainer::PFTAU,
        EVTColContainer::PFJET,
        EVTColContainer::CALOJET
    };

    std::set<unsigned int> objsType;
    // The object to deal has to be entered via the config .py
    for (unsigned int i = 0; i < objSize; ++i) {
    //std::cout << "(5) Gonna call with " << objtriggernames[i] << std::endl;
        std::string objTypeStr = EVTColContainer::getTypeString(objtriggernames[i]);
        // Check if it is needed this object for this trigger
        if (! TString(hltPath).Contains(objTypeStr)) {
            continue;
        }

        objsType.insert(objtriggernames[i]);
    }

    return std::vector<unsigned int>(objsType.begin(), objsType.end());
}

void HLTExoticaSubAnalysis::initSelector ( const unsigned int &  objtype )

private

Initializes the selectors of the objects based on which object it is.

Definition at line 1145 of file HLTExoticaSubAnalysis.cc.

References _genMETSelector, _l1METSelector, _recCaloJetSelector, _recCaloMETSelector, _recCaloMHTSelector, _recCut, _recElecSelector, _recMETSelector, _recMuonSelector, _recMuonTrkSelector, _recPFJetSelector, _recPFMETSelector, _recPFMHTSelector, _recPFTauSelector, _recPhotonSelector, _recTrackSelector, EVTColContainer::CALOJET, EVTColContainer::CALOMET, EVTColContainer::CALOMHT, EVTColContainer::ELEC, EVTColContainer::GENMET, EVTColContainer::L1MET, LogDebug, EVTColContainer::MET, EVTColContainer::MUON, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::PFMET, EVTColContainer::PFMHT, EVTColContainer::PFTAU, EVTColContainer::PHOTON, and EVTColContainer::TRACK.

Referenced by analyze().

{

    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::initSelector()";

    if (objtype == EVTColContainer::MUON && _recMuonSelector == 0) {
        _recMuonSelector = new StringCutObjectSelector<reco::Muon>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::MUTRK && _recMuonTrkSelector == 0) {
        _recMuonTrkSelector = new StringCutObjectSelector<reco::Track>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::TRACK && _recTrackSelector == 0) {
        _recTrackSelector = new StringCutObjectSelector<reco::Track>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::ELEC && _recElecSelector == 0) {
        _recElecSelector = new StringCutObjectSelector<reco::GsfElectron>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::PHOTON && _recPhotonSelector == 0) {
        _recPhotonSelector = new StringCutObjectSelector<reco::Photon>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::MET && _recMETSelector == 0) {
        _recMETSelector = new StringCutObjectSelector<reco::MET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::PFMET && _recPFMETSelector == 0) {
        _recPFMETSelector = new StringCutObjectSelector<reco::PFMET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::PFMHT && _recPFMHTSelector == 0) {
        _recPFMHTSelector = new StringCutObjectSelector<reco::PFMET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::GENMET && _genMETSelector == 0) {
        _genMETSelector = new StringCutObjectSelector<reco::GenMET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::CALOMET && _recCaloMETSelector == 0) {
        _recCaloMETSelector = new StringCutObjectSelector<reco::CaloMET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::CALOMHT && _recCaloMHTSelector == 0) {
        _recCaloMHTSelector = new StringCutObjectSelector<reco::CaloMET>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::L1MET && _l1METSelector == 0) {
        _l1METSelector = new StringCutObjectSelector<l1extra::L1EtMissParticle>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::PFTAU && _recPFTauSelector == 0) {
        _recPFTauSelector = new StringCutObjectSelector<reco::PFTau>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::PFJET && _recPFJetSelector == 0) {
        _recPFJetSelector = new StringCutObjectSelector<reco::PFJet>(_recCut[objtype]);
    } else if (objtype == EVTColContainer::CALOJET && _recCaloJetSelector == 0) {
        _recCaloJetSelector = new StringCutObjectSelector<reco::CaloJet>(_recCut[objtype]);
    }
    /* else
    {
    FIXME: ERROR NOT IMPLEMENTED
    }*/
}

void HLTExoticaSubAnalysis::insertCandidates ( const unsigned int &  objtype, const EVTColContainer *  col, std::vector< reco::LeafCandidate > *  matches, std::map< int, double > &  theSumEt, std::map< int, std::vector< const reco::Track * > > &  trkObjs  )

private

This function applies the selectors initialized previously to the objects, and matches the passing objects to HLT objects.

This is a special case. Passing a PFMET* to the constructor of MatchStruct will trigger the usage of the special constructor which also sets the sumEt member.

Definition at line 1188 of file HLTExoticaSubAnalysis.cc.

References _genMETSelector, _l1METSelector, _recCaloJetSelector, _recCaloMETSelector, _recCaloMHTSelector, _recElecSelector, _recMuonSelector, _recMuonTrkSelector, _recPFJetSelector, _recPFMETSelector, _recPFMHTSelector, _recPFTauSelector, _recPhotonSelector, _recTrackSelector, EVTColContainer::CALOJET, EVTColContainer::caloJets, EVTColContainer::CALOMET, EVTColContainer::caloMETs, EVTColContainer::CALOMHT, EVTColContainer::caloMHTs, EVTColContainer::ELEC, EVTColContainer::electrons, EVTColContainer::GENMET, EVTColContainer::genMETs, i, EVTColContainer::L1MET, EVTColContainer::l1METs, LogDebug, visualization-live-secondInstance_cfg::m, EVTColContainer::MUON, EVTColContainer::muons, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::pfJets, EVTColContainer::PFMET, EVTColContainer::pfMETs, EVTColContainer::PFMHT, EVTColContainer::pfMHTs, EVTColContainer::PFTAU, EVTColContainer::pfTaus, EVTColContainer::PHOTON, EVTColContainer::photons, EVTColContainer::TRACK, and EVTColContainer::tracks.

Referenced by analyze().

{
    
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::insertCandidates()"; 

    theSumEt[objType] = -1;
    
    if (objType == EVTColContainer::MUON) {
        for (size_t i = 0; i < cols->muons->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting muon " << i ;
        if (_recMuonSelector->operator()(cols->muons->at(i))) {
        reco::LeafCandidate m(0, cols->muons->at(i).p4(), cols->muons->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        
                // for making dxy plots
                trkObjs[objType].push_back(cols->muons->at(i).bestTrack());
        }
        }
    } else if (objType == EVTColContainer::MUTRK) {
        for (size_t i = 0; i < cols->tracks->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting muonTrack " << i ;
        if (_recMuonTrkSelector->operator()(cols->tracks->at(i))) {
        ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>> mom4;
                ROOT::Math::XYZVector mom3 = cols->tracks->at(i).momentum();
                mom4.SetXYZT(mom3.x(),mom3.y(),mom3.z(),mom3.r());
        reco::LeafCandidate m(0, mom4, cols->tracks->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
            }
        }
    } else if (objType == EVTColContainer::TRACK) {
        for (size_t i = 0; i < cols->tracks->size(); i++) {
            LogDebug("ExoticaValidation") << "Inserting Track " << i ;
            if (_recTrackSelector->operator()(cols->tracks->at(i))) {
                ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>> mom4;
                ROOT::Math::XYZVector mom3 = cols->tracks->at(i).momentum();
                mom4.SetXYZT(mom3.x(),mom3.y(),mom3.z(),mom3.r());
                reco::LeafCandidate m(0, mom4, cols->tracks->at(i).vertex(), objType, 0, true);
                matches->push_back(m);
            }
        }
    } else if (objType == EVTColContainer::ELEC) {
        for (size_t i = 0; i < cols->electrons->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting electron " << i ;
            if (_recElecSelector->operator()(cols->electrons->at(i))) {
        reco::LeafCandidate m(0, cols->electrons->at(i).p4(), cols->electrons->at(i).vertex(), objType, 0, true);
        matches->push_back(m);

                // for making dxy plots
                trkObjs[objType].push_back(cols->electrons->at(i).bestTrack());
        }
        }
    } else if (objType == EVTColContainer::PHOTON) {
        for (size_t i = 0; i < cols->photons->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting photon " << i ;
            if (_recPhotonSelector->operator()(cols->photons->at(i))) {
        reco::LeafCandidate m(0, cols->photons->at(i).p4(), cols->photons->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
            }
        }
    } else if (objType == EVTColContainer::PFMET) {
        for (size_t i = 0; i < cols->pfMETs->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting PFMET " << i ;
        if (_recPFMETSelector->operator()(cols->pfMETs->at(i))) {
        reco::LeafCandidate m(0, cols->pfMETs->at(i).p4(), cols->pfMETs->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->pfMETs->at(i).sumEt();
        }
        }
    } else if (objType == EVTColContainer::PFMHT) {
        for (size_t i = 0; i < cols->pfMHTs->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting PFMHT " << i ;
        if (_recPFMHTSelector->operator()(cols->pfMHTs->at(i))) {
        reco::LeafCandidate m(0, cols->pfMHTs->at(i).p4(), cols->pfMHTs->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->pfMHTs->at(i).sumEt();
        }
        }
    } else if (objType == EVTColContainer::GENMET) {
        for (size_t i = 0; i < cols->genMETs->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting GENMET " << i ;
        if (_genMETSelector->operator()(cols->genMETs->at(i))) {
        reco::LeafCandidate m(0, cols->genMETs->at(i).p4(), cols->genMETs->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->genMETs->at(i).sumEt();
        }
        }
    } else if (objType == EVTColContainer::CALOMET) {
        for (size_t i = 0; i < cols->caloMETs->size(); i++) {
            LogDebug("ExoticaValidation") << "Inserting CALOMET " << i ;
            if (_recCaloMETSelector->operator()(cols->caloMETs->at(i))) {
                reco::LeafCandidate m(0, cols->caloMETs->at(i).p4(), cols->caloMETs->at(i).vertex(), objType, 0, true);
                matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->caloMETs->at(i).sumEt();
            }
        }
    } else if (objType == EVTColContainer::CALOMHT) {
        for (size_t i = 0; i < cols->caloMHTs->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting CaloMHT " << i ;
        if (_recCaloMHTSelector->operator()(cols->caloMHTs->at(i))) {
        reco::LeafCandidate m(0, cols->caloMHTs->at(i).p4(), cols->caloMHTs->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->caloMHTs->at(i).sumEt();
        }
        }
    } else if (objType == EVTColContainer::L1MET) {
        for (size_t i = 0; i < cols->l1METs->size(); i++) {
            LogDebug("ExoticaValidation") << "Inserting L1MET " << i ;
            if (_l1METSelector->operator()(cols->l1METs->at(i))) {
                reco::LeafCandidate m(0, cols->l1METs->at(i).p4(), cols->l1METs->at(i).vertex(), objType, 0, true);
                matches->push_back(m);
        if(i==0) theSumEt[objType] = cols->l1METs->at(i).etTotal();
            }
        }
    } else if (objType == EVTColContainer::PFTAU) {
        for (size_t i = 0; i < cols->pfTaus->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting PFtau " << i ;
        if (_recPFTauSelector->operator()(cols->pfTaus->at(i))) {
        reco::LeafCandidate m(0, cols->pfTaus->at(i).p4(), cols->pfTaus->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
        }
        }
    } else if (objType == EVTColContainer::PFJET) {
        for (size_t i = 0; i < cols->pfJets->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting jet " << i ;
            if (_recPFJetSelector->operator()(cols->pfJets->at(i))) {
        reco::LeafCandidate m(0, cols->pfJets->at(i).p4(), cols->pfJets->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
            }
        }
    } else if (objType == EVTColContainer::CALOJET) {
        for (size_t i = 0; i < cols->caloJets->size(); i++) {
        LogDebug("ExoticaValidation") << "Inserting jet " << i ;
            if (_recCaloJetSelector->operator()(cols->caloJets->at(i))) {
        reco::LeafCandidate m(0, cols->caloJets->at(i).p4(), cols->caloJets->at(i).vertex(), objType, 0, true);
        matches->push_back(m);
            }
        }
    }

    /* else
    {
    FIXME: ERROR NOT IMPLEMENTED
    }*/
}

void HLTExoticaSubAnalysis::registerConsumes ( edm::ConsumesCollector &  consCollector )

private

Registers consumption of objects.

Definition at line 851 of file HLTExoticaSubAnalysis.cc.

References _beamSpotLabel, _bsToken, _genParticleLabel, _genParticleToken, _recLabels, _tokens, _trigResultsLabel, _trigResultsToken, EVTColContainer::CALOJET, EVTColContainer::CALOMET, EVTColContainer::CALOMHT, edm::ConsumesCollector::consumes(), EVTColContainer::ELEC, EVTColContainer::GENMET, EVTColContainer::L1MET, LogDebug, EVTColContainer::MET, EVTColContainer::MUON, EVTColContainer::MUTRK, EVTColContainer::PFJET, EVTColContainer::PFMET, EVTColContainer::PFMHT, EVTColContainer::PFTAU, EVTColContainer::PHOTON, and EVTColContainer::TRACK.

Referenced by HLTExoticaSubAnalysis().

{
    // Register that we are getting genParticles
    _genParticleToken = iC.consumes<reco::GenParticleCollection>(_genParticleLabel);    
   
    // Register that we are getting the trigger results
    _trigResultsToken = iC.consumes<edm::TriggerResults>(_trigResultsLabel);

    // Register beamspot 
    _bsToken = iC.consumes<reco::BeamSpot>(_beamSpotLabel);

    // Loop over _recLabels, see what we need, and register.
    // Then save the registered token in _tokens.
    // Remember: _recLabels is a map<uint, edm::InputTag>
    // Remember: _tokens    is a map<uint, edm::EDGetToken>
    LogDebug("ExoticaValidation") << "We have got " << _recLabels.size() << "recLabels";
    for (std::map<unsigned int, edm::InputTag>::iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
        if (it->first == EVTColContainer::MUON) {
        edm::EDGetTokenT<reco::MuonCollection> particularToken = iC.consumes<reco::MuonCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::MUTRK) {
        edm::EDGetTokenT<reco::TrackCollection> particularToken = iC.consumes<reco::TrackCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
        }
    else if (it->first == EVTColContainer::TRACK) {
        edm::EDGetTokenT<reco::TrackCollection> particularToken = iC.consumes<reco::TrackCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
        }
    else if (it->first == EVTColContainer::ELEC) {
            edm::EDGetTokenT<reco::GsfElectronCollection> particularToken = iC.consumes<reco::GsfElectronCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::PHOTON) {
            edm::EDGetTokenT<reco::PhotonCollection> particularToken = iC.consumes<reco::PhotonCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::MET) {
            edm::EDGetTokenT<reco::METCollection> particularToken = iC.consumes<reco::METCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::PFMET) {
            edm::EDGetTokenT<reco::PFMETCollection> particularToken = iC.consumes<reco::PFMETCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::PFMHT) {
            edm::EDGetTokenT<reco::PFMETCollection> particularToken = iC.consumes<reco::PFMETCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::GENMET) {
            edm::EDGetTokenT<reco::GenMETCollection> particularToken = iC.consumes<reco::GenMETCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::CALOMET) {
            edm::EDGetTokenT<reco::CaloMETCollection> particularToken = iC.consumes<reco::CaloMETCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::CALOMHT) {
            edm::EDGetTokenT<reco::CaloMETCollection> particularToken = iC.consumes<reco::CaloMETCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::L1MET) {
            edm::EDGetTokenT<l1extra::L1EtMissParticleCollection> particularToken = iC.consumes<l1extra::L1EtMissParticleCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::PFTAU) {
            edm::EDGetTokenT<reco::PFTauCollection> particularToken = iC.consumes<reco::PFTauCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::PFJET) {
            edm::EDGetTokenT<reco::PFJetCollection> particularToken = iC.consumes<reco::PFJetCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else if (it->first == EVTColContainer::CALOJET) {
            edm::EDGetTokenT<reco::CaloJetCollection> particularToken = iC.consumes<reco::CaloJetCollection>(it->second);
        edm::EDGetToken token(particularToken);
        _tokens[it->first] = token;
    } 
    else {
            edm::LogError("ExoticaValidation") << "HLTExoticaSubAnalysis::registerConsumes"
                                               << " NOT IMPLEMENTED (yet) ERROR: '" << it->second.label() << "'";
        }
    }
    
}

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

Method to book all relevant histograms in the DQMStore. Uses the IBooker interface for thread safety. Intended to be called from master object.

Definition at line 203 of file HLTExoticaSubAnalysis.cc.

References _analysisname, _drop_pt2, _drop_pt3, _plotters, _recLabels, bookHist(), EVTColContainer::ELEC, EVTColContainer::getTypeString(), i, LogDebug, EVTColContainer::MUON, DQMStore::IBooker::setCurrentFolder(), source, and AlCaHLTBitMon_QueryRunRegistry::string.

{
    
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::subAnalysisBookHistos()";

    // Create the folder structure inside HLT/Exotica
    std::string baseDir = "HLT/Exotica/" + _analysisname + "/";
    iBooker.setCurrentFolder(baseDir);

    // Book the gen/reco analysis-dependent histograms (denominators)
    for (std::map<unsigned int, edm::InputTag>::const_iterator it = _recLabels.begin();
         it != _recLabels.end(); ++it) {
    const std::string objStr = EVTColContainer::getTypeString(it->first);
        std::vector<std::string> sources(2);
        sources[0] = "gen";
        sources[1] = "rec";
    
        for (size_t i = 0; i < sources.size(); i++) {
      std::string source = sources[i];

          if ( source == "gen" ) {
            if ( TString(objStr).Contains("MET") ||
                 TString(objStr).Contains("MHT") ||
                 TString(objStr).Contains("Jet")    ) {
              continue;
            } else {
              bookHist(iBooker, source, objStr, "MaxPt1");
              if (!_drop_pt2) bookHist(iBooker, source, objStr, "MaxPt2");
              if (!_drop_pt3) bookHist(iBooker, source, objStr, "MaxPt3");
              bookHist(iBooker, source, objStr, "Eta");
              bookHist(iBooker, source, objStr, "Phi");

              // If the target is electron or muon,
              // we will add Dxy plots.
              if ( it->first == EVTColContainer::ELEC ||
                   it->first == EVTColContainer::MUON    ) { 
                 bookHist(iBooker, source, objStr, "Dxy");
              }
            }
          } else { // reco
            if ( TString(objStr).Contains("MET") ||
                 TString(objStr).Contains("MHT")    ) {
              bookHist(iBooker, source, objStr, "MaxPt1");
              bookHist(iBooker, source, objStr, "SumEt");
            } else {
              bookHist(iBooker, source, objStr, "MaxPt1");
              if (!_drop_pt2) bookHist(iBooker, source, objStr, "MaxPt2");
              if (!_drop_pt3) bookHist(iBooker, source, objStr, "MaxPt3");
              bookHist(iBooker, source, objStr, "Eta");
              bookHist(iBooker, source, objStr, "Phi");

              // If the target is electron or muon,
              // we will add Dxy plots.
              if ( it->first == EVTColContainer::ELEC ||
                   it->first == EVTColContainer::MUON    ) { 
                 bookHist(iBooker, source, objStr, "Dxy");
              }
            }
          }

        }
    } // closes loop in _recLabels

    // Call the plotterBookHistos() (which books all the path dependent histograms)
    LogDebug("ExoticaValidation") << "                        number of plotters = " << _plotters.size();
    for (std::vector<HLTExoticaPlotter>::iterator it = _plotters.begin();
         it != _plotters.end(); ++it) {
    it->plotterBookHistos(iBooker, iRun, iSetup);
    }
}

Member Data Documentation

std::string HLTExoticaSubAnalysis::_analysisname

private

The name of this sub-analysis.

Definition at line 123 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun(), endRun(), getNamesOfObjects(), and subAnalysisBookHistos().

edm::InputTag HLTExoticaSubAnalysis::_beamSpotLabel

private

Definition at line 139 of file HLTExoticaSubAnalysis.h.

Referenced by registerConsumes().

edm::EDGetTokenT<reco::BeamSpot> HLTExoticaSubAnalysis::_bsToken

private

Definition at line 144 of file HLTExoticaSubAnalysis.h.

Referenced by getHandlesToObjects(), and registerConsumes().

bool HLTExoticaSubAnalysis::_drop_pt2

private

Definition at line 155 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), HLTExoticaSubAnalysis(), and subAnalysisBookHistos().

bool HLTExoticaSubAnalysis::_drop_pt3

private

Definition at line 156 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), HLTExoticaSubAnalysis(), and subAnalysisBookHistos().

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

private

Structure of the MonitorElements.

Definition at line 194 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and fillHist().

std::map<unsigned int, std::string> HLTExoticaSubAnalysis::_genCut

private

gen/rec objects cuts

Definition at line 159 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), and HLTExoticaSubAnalysis().

std::map<unsigned int, std::string> HLTExoticaSubAnalysis::_genCut_leading

private

gen/rec pt-leading objects cuts

Definition at line 162 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), and HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::GenMET>* HLTExoticaSubAnalysis::_genMETSelector

private

Definition at line 175 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

edm::InputTag HLTExoticaSubAnalysis::_genParticleLabel

private

Definition at line 137 of file HLTExoticaSubAnalysis.h.

Referenced by registerConsumes().

edm::EDGetTokenT<reco::GenParticleCollection> HLTExoticaSubAnalysis::_genParticleToken

private

And also the tokens to get the object collections.

Definition at line 142 of file HLTExoticaSubAnalysis.h.

Referenced by getHandlesToObjects(), and registerConsumes().

std::map<unsigned int, StringCutObjectSelector<reco::GenParticle> *> HLTExoticaSubAnalysis::_genSelectorMap

private

The concrete String selectors (use the string cuts introduced via the config python)

Definition at line 167 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), getNamesOfObjects(), and ~HLTExoticaSubAnalysis().

HLTConfigProvider HLTExoticaSubAnalysis::_hltConfig

private

Interface to the HLT information.

Definition at line 191 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun().

std::set<std::string> HLTExoticaSubAnalysis::_hltPaths

private

The hlt paths found in the hltConfig.

Definition at line 131 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun().

std::vector<std::string> HLTExoticaSubAnalysis::_hltPathsToCheck

private

The hlt paths to check for.

Definition at line 129 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun(), and HLTExoticaSubAnalysis().

std::string HLTExoticaSubAnalysis::_hltProcessName

private

The labels of the object collections to be used in this analysis.

Definition at line 136 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun().

StringCutObjectSelector<l1extra::L1EtMissParticle>* HLTExoticaSubAnalysis::_l1METSelector

private

Definition at line 178 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

unsigned int HLTExoticaSubAnalysis::_minCandidates

private

The minimum number of reco/gen candidates needed by the analysis.

Definition at line 126 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), and HLTExoticaSubAnalysis().

std::vector<double> HLTExoticaSubAnalysis::_parametersDxy

private

Definition at line 152 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and HLTExoticaSubAnalysis().

std::vector<double> HLTExoticaSubAnalysis::_parametersEta

private

Some kinematical parameters.

Definition at line 148 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and HLTExoticaSubAnalysis().

std::vector<double> HLTExoticaSubAnalysis::_parametersPhi

private

Definition at line 149 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and HLTExoticaSubAnalysis().

std::vector<double> HLTExoticaSubAnalysis::_parametersTurnOn

private

Definition at line 150 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and HLTExoticaSubAnalysis().

std::vector<double> HLTExoticaSubAnalysis::_parametersTurnOnSumEt

private

Definition at line 151 of file HLTExoticaSubAnalysis.h.

Referenced by bookHist(), and HLTExoticaSubAnalysis().

std::vector<HLTExoticaPlotter> HLTExoticaSubAnalysis::_plotters

private

The plotters: managers of each hlt path where the plots are done.

Definition at line 185 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), beginRun(), and subAnalysisBookHistos().

edm::ParameterSet HLTExoticaSubAnalysis::_pset

private

Internal, working copy of the PSet passed from above.

Definition at line 120 of file HLTExoticaSubAnalysis.h.

Referenced by beginRun(), and HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::CaloJet>* HLTExoticaSubAnalysis::_recCaloJetSelector

private

Definition at line 182 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::CaloMET>* HLTExoticaSubAnalysis::_recCaloMETSelector

private

Definition at line 176 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::CaloMET>* HLTExoticaSubAnalysis::_recCaloMHTSelector

private

Definition at line 177 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

std::map<unsigned int, std::string> HLTExoticaSubAnalysis::_recCut

private

Definition at line 160 of file HLTExoticaSubAnalysis.h.

Referenced by HLTExoticaSubAnalysis(), and initSelector().

std::map<unsigned int, std::string> HLTExoticaSubAnalysis::_recCut_leading

private

Definition at line 163 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), and HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::GsfElectron>* HLTExoticaSubAnalysis::_recElecSelector

private

Definition at line 171 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

std::map<unsigned int, edm::InputTag> HLTExoticaSubAnalysis::_recLabels

private

Definition at line 140 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), beginRun(), getNamesOfObjects(), HLTExoticaSubAnalysis(), registerConsumes(), and subAnalysisBookHistos().

StringCutObjectSelector<reco::MET>* HLTExoticaSubAnalysis::_recMETSelector

private

Definition at line 172 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::Muon>* HLTExoticaSubAnalysis::_recMuonSelector

private

Definition at line 168 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::Track>* HLTExoticaSubAnalysis::_recMuonTrkSelector

private

Definition at line 169 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::PFJet>* HLTExoticaSubAnalysis::_recPFJetSelector

private

Definition at line 181 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::PFMET>* HLTExoticaSubAnalysis::_recPFMETSelector

private

Definition at line 173 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::PFMET>* HLTExoticaSubAnalysis::_recPFMHTSelector

private

Definition at line 174 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::PFTau>* HLTExoticaSubAnalysis::_recPFTauSelector

private

Definition at line 179 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::Photon>* HLTExoticaSubAnalysis::_recPhotonSelector

private

Definition at line 180 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

StringCutObjectSelector<reco::Track>* HLTExoticaSubAnalysis::_recTrackSelector

private

Definition at line 170 of file HLTExoticaSubAnalysis.h.

Referenced by initSelector(), insertCandidates(), and ~HLTExoticaSubAnalysis().

std::map<std::string, std::string> HLTExoticaSubAnalysis::_shortpath2long

private

Relation between the short and long versions of the path.

Definition at line 133 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), and beginRun().

std::map<unsigned int, edm::EDGetToken> HLTExoticaSubAnalysis::_tokens

private

Definition at line 145 of file HLTExoticaSubAnalysis.h.

Referenced by getHandlesToObjects(), and registerConsumes().

std::map<std::string,int> HLTExoticaSubAnalysis::_triggerCounter

private

counting HLT passed events

Definition at line 188 of file HLTExoticaSubAnalysis.h.

Referenced by analyze(), beginRun(), and endRun().

edm::InputTag HLTExoticaSubAnalysis::_trigResultsLabel

private

Definition at line 138 of file HLTExoticaSubAnalysis.h.

Referenced by registerConsumes().

edm::EDGetTokenT<edm::TriggerResults> HLTExoticaSubAnalysis::_trigResultsToken

private

Definition at line 143 of file HLTExoticaSubAnalysis.h.

Referenced by getHandlesToObjects(), and registerConsumes().