HLTHiggsSubAnalysis.h

Go to the documentation of this file.

#ifndef HLTriggerOffline_Higgs_HLTHiggsSubAnalysis_H
#define HLTriggerOffline_Higgs_HLTHiggsSubAnalysis_H

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"

#include "SimDataFormats/PileupSummaryInfo/interface/PileupSummaryInfo.h"
#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETFwd.h"
#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/TauReco/interface/PFTauFwd.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/BTauReco/interface/JetTag.h"
#include "CommonTools/Utils/interface/StringCutObjectSelector.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"

#include "HLTHiggsPlotter.h"

#include <vector>
#include <set>
#include <map>
#include <cstring>

struct EVTColContainer;

class HLTHiggsSubAnalysis {
public:
  typedef dqm::legacy::DQMStore DQMStore;
  typedef dqm::legacy::MonitorElement MonitorElement;

  enum { GEN, RECO };

  HLTHiggsSubAnalysis(const edm::ParameterSet &pset, const std::string &analysisname, edm::ConsumesCollector &&iC);
  ~HLTHiggsSubAnalysis();
  void beginJob();
  void beginRun(const edm::Run &iRun, const edm::EventSetup &iEventSetup);
  void analyze(const edm::Event &iEvent, const edm::EventSetup &iEventSetup, EVTColContainer *cols);
  void bookHistograms(DQMStore::IBooker &);

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

private:
  void bookobjects(const edm::ParameterSet &anpset, edm::ConsumesCollector &iC);
  void initobjects(const edm::Event &iEvent, EVTColContainer *col);
  void InitSelector(const unsigned int &objtype);
  void initAndInsertJets(const edm::Event &iEvent, EVTColContainer *cols, std::vector<MatchStruct> *matches);
  void passJetCuts(std::vector<MatchStruct> *matches,
                   std::map<std::string, bool> &jetCutResult,
                   float &dEtaqq,
                   float &mqq,
                   float &dPhibb,
                   float &CSV1,
                   float &CSV2,
                   float &CSV3);
  void passOtherCuts(const std::vector<MatchStruct> &matches, std::map<std::string, bool> &jetCutResult);
  void insertcandidates(const unsigned int &objtype, const EVTColContainer *col, std::vector<MatchStruct> *matches);

  void bookHist(const std::string &source, const std::string &objType, const std::string &variable, DQMStore::IBooker &);
  void fillHist(const std::string &source, const std::string &objType, const std::string &variable, const float &value);

  edm::ParameterSet _pset;

  std::string _analysisname;

  unsigned int _minCandidates;

  double _HtJetPtMin;
  double _HtJetEtaMax;

  bool _bookHtPlots;

  std::string _hltProcessName;
  std::string _histDirectory;

  std::vector<std::string> _hltPathsToCheck;
  std::set<std::string> _hltPaths;

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

  // The name of the object collections to be used in this
  // analysis.
  edm::EDGetTokenT<reco::GenParticleCollection> _genParticleLabel;
  edm::EDGetTokenT<reco::GenJetCollection> _genJetLabel;
  edm::EDGetTokenT<reco::PFJetCollection> _recoHtJetLabel;

  std::map<unsigned int, std::string> _recLabels;
  edm::EDGetTokenT<reco::MuonCollection> _recLabelsMuon;
  edm::EDGetTokenT<reco::GsfElectronCollection> _recLabelsElec;
  edm::EDGetTokenT<reco::PhotonCollection> _recLabelsPhoton;
  edm::EDGetTokenT<reco::CaloMETCollection> _recLabelsCaloMET;
  edm::EDGetTokenT<reco::PFMETCollection> _recLabelsPFMET;
  edm::EDGetTokenT<reco::PFTauCollection> _recLabelsPFTau;
  edm::EDGetTokenT<reco::PFJetCollection> _recLabelsPFJet;
  edm::EDGetTokenT<reco::JetTagCollection> _recTagPFJet;
  edm::EDGetTokenT<std::vector<PileupSummaryInfo> > _puSummaryInfo;

  std::vector<double> _parametersEta;
  std::vector<double> _parametersPhi;
  std::vector<double> _parametersPu;
  std::vector<double> _parametersHt;
  std::vector<double> _parametersTurnOn;
  edm::EDGetTokenT<edm::TriggerResults> _trigResultsTag;

  std::map<unsigned int, double> _cutMinPt;
  std::map<unsigned int, double> _cutMaxEta;
  std::map<unsigned int, unsigned int> _cutMotherId;     //TO BE DEPRECATED (HLTMATCH)
  std::map<unsigned int, std::vector<double> > _cutsDr;  // TO BE DEPRECATED (HLTMATCH)
  std::map<unsigned int, std::string> _genCut;
  std::map<unsigned int, std::string> _recCut;

  std::map<unsigned int, StringCutObjectSelector<reco::GenParticle> *> _genSelectorMap;
  StringCutObjectSelector<reco::GenJet> *_genJetSelector;
  StringCutObjectSelector<reco::Muon> *_recMuonSelector;
  StringCutObjectSelector<reco::GsfElectron> *_recElecSelector;
  StringCutObjectSelector<reco::CaloMET> *_recCaloMETSelector;
  StringCutObjectSelector<reco::PFMET> *_recPFMETSelector;
  StringCutObjectSelector<reco::PFTau> *_recPFTauSelector;
  StringCutObjectSelector<reco::Photon> *_recPhotonSelector;
  StringCutObjectSelector<reco::PFJet> *_recPFJetSelector;
  StringCutObjectSelector<reco::Track> *_recTrackSelector;

  //N-1 cut values
  std::vector<double> _NminOneCuts;
  bool _useNminOneCuts;
  unsigned int NptPlots;

  // The plotters: managers of each hlt path where the plots are done
  std::vector<HLTHiggsPlotter> _analyzers;

  HLTConfigProvider _hltConfig;
  std::map<std::string, MonitorElement *> _elements;
};
#endif