EmDQMReco.h

Go to the documentation of this file.

#ifndef HLTriggerOffline_Egamma_EmDQMReco_H
#define HLTriggerOffline_Egamma_EmDQMReco_H

// Base Class Headers
#include "CommonTools/Utils/interface/PtComparator.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Common/interface/AssociationMap.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/RefToBase.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
#include "HepMC/GenParticle.h"
#include "TDirectory.h"

#include <memory>
#include <vector>

class EmDQMReco;

template <class T>
class HistoFillerReco {
public:
  HistoFillerReco(EmDQMReco *d) : dqm(d){};
  ~HistoFillerReco(){};

  void fillHistos(edm::Handle<trigger::TriggerEventWithRefs> &triggerObj,
                  const edm::Event &iEvent,
                  unsigned int n,
                  std::vector<reco::Particle> &sortedReco,
                  bool plotReco,
                  bool plotMonpath);
  std::vector<edm::EDGetTokenT<edm::AssociationMap<edm::OneToValue<T, float>>>> isoNameTokens_;

private:
  EmDQMReco *dqm;
};

class EmDQMReco : public DQMEDAnalyzer {
  //----------------------------------------

  class FourVectorMonitorElements {
  public:
    FourVectorMonitorElements(EmDQMReco *_parent,
                              DQMStore::IBooker &iBooker,
                              const std::string &histogramNameTemplate,
                              const std::string &histogramTitleTemplate);

    void fill(const math::XYZTLorentzVector &momentum);

  private:
    EmDQMReco *parent;

    MonitorElement *etMonitorElement;
    MonitorElement *etaMonitorElement;
    MonitorElement *phiMonitorElement;
  };
  //----------------------------------------

public:
  friend class HistoFillerReco<reco::ElectronCollection>;
  friend class HistoFillerReco<reco::RecoEcalCandidateCollection>;
  friend class HistoFillerReco<l1extra::L1EmParticleCollection>;

  explicit EmDQMReco(const edm::ParameterSet &pset);

  ~EmDQMReco() override;

  // Operations
  void analyze(const edm::Event &event, const edm::EventSetup &) override;
  void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;

private:
  // Input from cfg file

  std::vector<edm::InputTag> theHLTCollectionLabels;

  unsigned int numOfHLTCollectionLabels;  // Will be size of above vector

  bool useHumanReadableHistTitles;
  std::vector<std::string> theHLTCollectionHumanNames;  // Human-readable names for the collections
  // edm::InputTag theL1Seed;
  std::vector<int> theHLTOutputTypes;
  std::vector<bool> plotiso;
  std::vector<std::vector<edm::InputTag>> isoNames;  // there has to be a better solution
  std::vector<std::pair<double, double>> plotBounds;
  std::string theHltName;
  HLTConfigProvider hltConfig_;
  bool isHltConfigInitialized_;

  //          Read from configuration file                  //
  // parameters for generator study
  unsigned int reqNum;
  int pdgGen;
  double recoEtaAcc;
  double recoEtAcc;
  // plotting paramters
  double plotEtaMax;
  double plotPtMin;
  double plotPtMax;
  double plotPhiMax;

  unsigned int plotBins;

  // preselction cuts
  // edm::InputTag recocutCollection_;
  unsigned int recocut_;

  std::string triggerNameRecoMonPath;

  std::string processNameRecoMonPath;

  edm::EDGetTokenT<reco::GsfElectronCollection> recoElectronsInput;
  edm::EDGetTokenT<std::vector<reco::SuperCluster>> recoObjectsEBT;
  edm::EDGetTokenT<std::vector<reco::SuperCluster>> recoObjectsEET;
  edm::EDGetTokenT<edm::TriggerResults> hltResultsT;
  edm::EDGetTokenT<trigger::TriggerEventWithRefs> triggerObjT;
  //          Create Histograms                             //

  // std::vector<MonitorElement*> etahist;
  //  std::vector<MonitorElement*> ethist;
  //  std::vector<MonitorElement*> phiHist;

  std::vector<std::unique_ptr<FourVectorMonitorElements>> standardHist;

  //  std::vector<MonitorElement*> etahistmatchreco;
  //  std::vector<MonitorElement*> ethistmatchreco;
  //  std::vector<MonitorElement*> phiHistMatchReco;
  std::vector<std::unique_ptr<FourVectorMonitorElements>> histMatchReco;

  //  std::vector<MonitorElement*> etahistmatchrecomonpath;
  //  std::vector<MonitorElement*> ethistmatchrecomonpath;
  //  std::vector<MonitorElement*> phiHistMatchRecoMonPath;
  std::vector<std::unique_ptr<FourVectorMonitorElements>> histMatchRecoMonPath;

  //  std::vector<MonitorElement*> histEtOfHltObjMatchToReco;
  //  std::vector<MonitorElement*> histEtaOfHltObjMatchToReco;
  //  std::vector<MonitorElement*> histPhiOfHltObjMatchToReco;
  std::vector<std::unique_ptr<FourVectorMonitorElements>> histHltObjMatchToReco;

  std::vector<MonitorElement *> etahistiso;
  std::vector<MonitorElement *> ethistiso;
  std::vector<MonitorElement *> phiHistIso;

  std::vector<MonitorElement *> etahistisomatchreco;
  std::vector<MonitorElement *> ethistisomatchreco;
  std::vector<MonitorElement *> phiHistIsoMatchReco;

  std::vector<MonitorElement *> histEtIsoOfHltObjMatchToReco;
  std::vector<MonitorElement *> histEtaIsoOfHltObjMatchToReco;
  std::vector<MonitorElement *> histPhiIsoOfHltObjMatchToReco;
  MonitorElement *totalreco;
  MonitorElement *totalmatchreco;

  //  MonitorElement* etreco;
  //  MonitorElement* etareco;
  //  MonitorElement* phiReco;
  std::unique_ptr<FourVectorMonitorElements> histReco;

  //  MonitorElement* etrecomonpath;
  //  MonitorElement* etarecomonpath;
  //  MonitorElement* phiRecoMonPath;
  std::unique_ptr<FourVectorMonitorElements> histRecoMonpath;

  //  MonitorElement* etahistmonpath;
  //  MonitorElement* ethistmonpath;
  //  MonitorElement* phiHistMonPath;
  std::unique_ptr<FourVectorMonitorElements> histMonpath;
  int eventnum;
  // int prescale;

  // interface to DQM framework
  std::string dirname_;

  HistoFillerReco<reco::ElectronCollection> *histoFillerEle;
  HistoFillerReco<reco::RecoEcalCandidateCollection> *histoFillerClu;
  HistoFillerReco<l1extra::L1EmParticleCollection> *histoFillerL1NonIso;
  HistoFillerReco<reco::RecoEcalCandidateCollection> *histoFillerPho;
  HistoFillerReco<l1extra::L1EmParticleCollection> *histoFillerL1Iso;

  // template <class T> void
  // fillHistos(edm::Handle<trigger::TriggerEventWithRefs>&,const edm::Event&
  // ,unsigned int, std::vector<reco::Particle>&, bool, bool);
  GreaterByPt<reco::Particle> pTComparator_;
  GreaterByPt<reco::GsfElectron> pTRecoComparator_;

  //----------------------------------------
};
#endif