HLTHiggsPlotter.h

Go to the documentation of this file.

#ifndef HLTriggerOffline_Higgs_HLTHiggsPlotter_H
#define HLTriggerOffline_Higgs_HLTHiggsPlotter_H

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"

#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
#include "DataFormats/RecoCandidate/interface/RecoChargedCandidateFwd.h"
#include "DataFormats/RecoCandidate/interface/RecoEcalCandidate.h"
#include "DataFormats/RecoCandidate/interface/RecoEcalCandidateFwd.h"
#include "CommonTools/Utils/interface/StringCutObjectSelector.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/JetReco/interface/GenJet.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/JetReco/interface/PFJet.h"

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

#include "HLTriggerOffline/Higgs/src/MatchStruct.cc"

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

const unsigned int kNull = (unsigned int)-1;

struct EVTColContainer;

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

  HLTHiggsPlotter(const edm::ParameterSet &pset,
                  const std::string &hltPath,
                  //const std::string & lastFilter,
                  const std::vector<unsigned int> &objectsType,
                  const unsigned int &minCandidates,
                  const std::vector<double> &NminOneCuts);
  ~HLTHiggsPlotter();
  void beginJob();
  void beginRun(const edm::Run &, const edm::EventSetup &);
  void bookHistograms(DQMStore::IBooker &, const bool &);
  void analyze(const bool &isPassTrigger,
               const std::string &source,
               const std::vector<MatchStruct> &matches,
               const unsigned int &minCandidates);
  void analyze(const bool &isPassTrigger,
               const std::string &source,
               const std::vector<MatchStruct> &matches,
               std::map<std::string, bool> &nMinOne,
               const float &dEtaqq,
               const float &mqq,
               const float &dPhibb,
               const float &CSV1,
               const float &CSV2,
               const float &CSV3,
               const bool &passAllCuts);

  inline const std::string gethltpath() const { return _hltPath; }

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

  std::string _hltPath;
  //std::string _lastFilter;
  std::string _hltProcessName;

  std::set<unsigned int> _objectsType;
  // Number of objects (elec,muons, ...) needed in the hlt path
  unsigned int _nObjects;

  std::vector<double> _parametersEta;
  std::vector<double> _parametersPhi;
  std::vector<double> _parametersTurnOn;

  std::map<unsigned int, double> _cutMinPt;
  std::map<unsigned int, double> _cutMaxEta;
  std::map<unsigned int, unsigned int> _cutMotherId;
  std::map<unsigned int, std::vector<double> > _cutsDr;

  // The amount of Pt plots needed for the hlt path
  unsigned int _NptPlots;

  //N-1 cut values
  std::vector<double> _NminOneCuts;

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