HLTExoticaSubAnalysis.cc

Go to the documentation of this file.

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

#include "CommonTools/Utils/interface/PtComparator.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "HLTriggerOffline/Exotica/interface/HLTExoticaSubAnalysis.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 <algorithm>
#include <set>

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;
      }
    }

    int jel = 0;
    int jmu = 0;
    int jmutrk = 0;

    // jel, jmu and jmutrk are being used as a dedicated counters to avoid getting
    // non-existent elements inside trkObjs[11], trkObjs[13] and trkObjs[130], respectively
    // more information in the issue https://github.com/cms-sw/cmssw/issues/32550

    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 (objType == 11) {
        float dxyRec = trkObjs[objType].at(jel)->dxy(cols->bs->position());
        this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
        dxys.push_back(dxyRec);
        ++jel;
      }

      if (objType == 13) {
        float dxyRec = trkObjs[objType].at(jmu)->dxy(cols->bs->position());
        this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
        dxys.push_back(dxyRec);
        ++jmu;
      }

      if (objType == 130) {
        float dxyRec = trkObjs[objType].at(jmutrk)->dxy(cols->bs->position());
        this->fillHist("rec", objTypeStr, "Dxy", dxyRec);
        dxys.push_back(dxyRec);
        ++jmutrk;
      }

    }  // 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();
}