HLTExoticaSubAnalysis.cc

Go to the documentation of this file.

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Utilities/interface/TypeID.h"
#include "DataFormats/Common/interface/Handle.h"

#include "HLTriggerOffline/Exotica/interface/HLTExoticaSubAnalysis.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "CommonTools/Utils/interface/PtComparator.h"
#include "HLTriggerOffline/Exotica/src/EVTColContainer.cc"

#include "TPRegexp.h"
#include "TString.h"

#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Common/interface/TriggerNames.h"

#include <set>
#include <algorithm>

static constexpr int verbose=0;

HLTExoticaSubAnalysis::HLTExoticaSubAnalysis(const edm::ParameterSet & pset,
                                             const std::string & analysisname,
                         edm::ConsumesCollector && consCollector) :
    _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(nullptr),
    _recMuonTrkSelector(nullptr),
    _recTrackSelector(nullptr),
    _recElecSelector(nullptr),
    _recMETSelector(nullptr),
    _recPFMETSelector(nullptr),
    _recPFMHTSelector(nullptr),
    _genMETSelector(nullptr),
    _recCaloMETSelector(nullptr),
    _recCaloMHTSelector(nullptr),
    _l1METSelector(nullptr),
    _recPFTauSelector(nullptr),
    _recPhotonSelector(nullptr),
    _recPFJetSelector(nullptr),
    _recCaloJetSelector(nullptr)
{

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


void HLTExoticaSubAnalysis::beginJob()
{
}


// 2014-02-03 -- Thiago
// Due to the fact that the DQM has to be thread safe now, we have to do things differently:
// 1) Implement the bookHistograms() method in the container class
// 2) Make the iBooker from above be known to this class
// 3) Separate all booking histograms routines in this and any auxiliary classe to be called
// from bookHistograms() in the container class
void HLTExoticaSubAnalysis::subAnalysisBookHistos(DQMStore::IBooker &iBooker,
                                                  const edm::Run & iRun,
                                                  const edm::EventSetup & iSetup)
{
    
    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 ||  
                   it->first == EVTColContainer::MUTRK   ) { 
                 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 ||  
                   it->first == EVTColContainer::MUTRK    ) { 
                 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);
    }
}

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

    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::analyze(const edm::Event & iEvent, const edm::EventSetup & iSetup, EVTColContainer * cols)
{
    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::MUTRK ||
             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";
    }
    } 

} 


// Return the objects (muons,electrons,photons,...) needed by a hlt path.
const std::vector<unsigned int> HLTExoticaSubAnalysis::getObjectsType(const std::string & hltPath) const
{
    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());
}

// Booking the maps: recLabels and genParticle selectors
void HLTExoticaSubAnalysis::getNamesOfObjects(const edm::ParameterSet & anpset)
{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::getNamesOfObjects()";

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

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

// Register consumption of objects.
// I have chosen to centralize all consumes() calls here.
void HLTExoticaSubAnalysis::registerConsumes(edm::ConsumesCollector & iC)
{
    // 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() << "'";
        }
    }
    
}

// Setting the collections of objects in EVTColContainer
void HLTExoticaSubAnalysis::getHandlesToObjects(const edm::Event & iEvent, EVTColContainer * col)
{
    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 << "'";
        }
    }
}

// Booking the histograms, and putting them in DQM
void HLTExoticaSubAnalysis::bookHist(DQMStore::IBooker & iBooker,
                     const std::string & source,
                                     const std::string & objType, 
                     const std::string & variable)
{
    LogDebug("ExoticaValidation") << "In HLTExoticaSubAnalysis::bookHist()";
    std::string sourceUpper = source;
    sourceUpper[0] = std::toupper(sourceUpper[0]);
    std::string name = source + objType + variable ;
    TH1F * h = nullptr;
    
    if (variable.find("SumEt") != std::string::npos) {
      std::string title = "Sum ET of " + sourceUpper + " " + objType;
      const size_t nBins = _parametersTurnOnSumEt.size() - 1;
      float * edges = new float[nBins + 1];
      for (size_t i = 0; i < nBins + 1; i++) {
    edges[i] = _parametersTurnOnSumEt[i];
      }
      h = new TH1F(name.c_str(), title.c_str(), nBins, edges);
      delete[] edges;
    }

    else if (variable.find("Dxy") != std::string::npos) {
      std::string title = "Dxy " + sourceUpper + " " + objType;
      int    nBins = _parametersDxy[0];
      double min   = _parametersDxy[1];
      double max   = _parametersDxy[2];
      h = new TH1F(name.c_str(), title.c_str(), nBins, min, max);
    }

    else if (variable.find("MaxPt") != std::string::npos) {
      std::string desc = (variable == "MaxPt1") ? "Leading" : "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;
}

// Fill the histograms
void HLTExoticaSubAnalysis::fillHist(const std::string & source, 
                     const std::string & objType, 
                     const std::string & variable, 
                     const float & value)
{
    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";
}

// Initialize the selectors
void HLTExoticaSubAnalysis::initSelector(const unsigned int & objtype)
{

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

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

// Insert the HLT candidates
void HLTExoticaSubAnalysis::insertCandidates(const unsigned int & objType, const EVTColContainer * cols, std::vector<reco::LeafCandidate> * matches, 
                                             std::map<int,double> & theSumEt, std::map<int,std::vector<const reco::Track*> > & trkObjs )
{
    
    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);

                // for making dxy plots
                trkObjs[objType].push_back(&cols->tracks->at(i));
            }
        }
    } 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::endRun() {
    // 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();
}