HLTExoticaSubAnalysis.h

Go to the documentation of this file.

#ifndef HLTriggerOffline_Exotica_HLTExoticaSubAnalysis_H
#define HLTriggerOffline_Exotica_HLTExoticaSubAnalysis_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/Framework/interface/ConsumesCollector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/Candidate/interface/LeafCandidate.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"
#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/METReco/interface/GenMETCollection.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/L1Trigger/interface/L1EtMissParticle.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/TauReco/interface/PFTauFwd.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "CommonTools/Utils/interface/StringCutObjectSelector.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"

#include "HLTriggerOffline/Exotica/interface/HLTExoticaPlotter.h"

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

class EVTColContainer;

class HLTExoticaSubAnalysis {
public:
    enum class Level{
        GEN = 98,
        RECO = 99
    };

    HLTExoticaSubAnalysis(const edm::ParameterSet & pset,
                          const std::string & analysisname,
              edm::ConsumesCollector && consCollector);
    ~HLTExoticaSubAnalysis();
    void beginJob();
    void beginRun(const edm::Run & iRun, const edm::EventSetup & iEventSetup);
    void endRun();

    void subAnalysisBookHistos(DQMStore::IBooker &iBooker, const edm::Run & iRun, const edm::EventSetup & iSetup);

    void analyze(const edm::Event & iEvent, const edm::EventSetup & iEventSetup, EVTColContainer * cols);

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

    void getNamesOfObjects(const edm::ParameterSet & anpset);
    void registerConsumes(edm::ConsumesCollector & consCollector);
    void getHandlesToObjects(const edm::Event & iEvent, EVTColContainer * col);
    void initSelector(const unsigned int & objtype);
    void insertCandidates(const unsigned int & objtype, const EVTColContainer * col,
                          std::vector<reco::LeafCandidate> * matches, 
                          std::map<int,double> & theSumEt,
                          std::map<int, std::vector<const reco::Track*> > & trkObjs);

    void bookHist(DQMStore::IBooker &iBooker, const std::string & source, const std::string & objType,
                  const std::string & variable);
    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;

    std::vector<std::string> _hltPathsToCheck;
    std::set<std::string> _hltPaths;
    std::map<std::string, std::string> _shortpath2long;

    std::string _hltProcessName;
    edm::InputTag _genParticleLabel;
    edm::InputTag _trigResultsLabel;
    edm::InputTag _beamSpotLabel;
    std::map<unsigned int, edm::InputTag> _recLabels;
    edm::EDGetTokenT<reco::GenParticleCollection> _genParticleToken;
    edm::EDGetTokenT<edm::TriggerResults> _trigResultsToken;
    edm::EDGetTokenT<reco::BeamSpot> _bsToken;
    std::map<unsigned int, edm::EDGetToken> _tokens;

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

    // flag to switch off
    bool _drop_pt2;
    bool _drop_pt3;

    std::map<unsigned int, std::string> _genCut;
    std::map<unsigned int, std::string> _recCut;
    std::map<unsigned int, std::string> _genCut_leading;
    std::map<unsigned int, std::string> _recCut_leading;

    std::map<unsigned int, StringCutObjectSelector<reco::GenParticle> *> _genSelectorMap;
    StringCutObjectSelector<reco::Muon>        * _recMuonSelector;
    StringCutObjectSelector<reco::Track>       * _recMuonTrkSelector;
    StringCutObjectSelector<reco::Track>       * _recTrackSelector;
    StringCutObjectSelector<reco::GsfElectron> * _recElecSelector;
    StringCutObjectSelector<reco::MET>         * _recMETSelector;
    StringCutObjectSelector<reco::PFMET>       * _recPFMETSelector;
    StringCutObjectSelector<reco::PFMET>       * _recPFMHTSelector;
    StringCutObjectSelector<reco::GenMET>      * _genMETSelector;
    StringCutObjectSelector<reco::CaloMET>     * _recCaloMETSelector;
    StringCutObjectSelector<reco::CaloMET>     * _recCaloMHTSelector;
    StringCutObjectSelector<l1extra::L1EtMissParticle> * _l1METSelector;
    StringCutObjectSelector<reco::PFTau>       * _recPFTauSelector;
    StringCutObjectSelector<reco::Photon>      * _recPhotonSelector;
    StringCutObjectSelector<reco::PFJet>       * _recPFJetSelector;
    StringCutObjectSelector<reco::CaloJet>     * _recCaloJetSelector;

    std::vector<HLTExoticaPlotter> _plotters;

    std::map<std::string,int> _triggerCounter;

    HLTConfigProvider _hltConfig;

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


#endif