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 "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/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 "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/Higgs/interface/HLTHiggsPlotter.h"


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

class EVTColContainer;

class HLTHiggsSubAnalysis 
{   
        public:
        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);

        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 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);
        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::string _hltProcessName;
        
        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;

        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::PFTauCollection> _recLabelsPFTau;

        
            std::vector<double> _parametersEta;
            std::vector<double> _parametersPhi;
            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::Muon>        * _recMuonSelector;
            StringCutObjectSelector<reco::GsfElectron> * _recElecSelector;
            StringCutObjectSelector<reco::CaloMET>     * _recCaloMETSelector;
            StringCutObjectSelector<reco::PFTau>       * _recPFTauSelector;
            StringCutObjectSelector<reco::Photon>      * _recPhotonSelector;
            StringCutObjectSelector<reco::Track>       * _recTrackSelector;
        

        // The plotters: managers of each hlt path where the plots are done
        std::vector<HLTHiggsPlotter> _analyzers;
        
        HLTConfigProvider _hltConfig;
        
            DQMStore* _dbe;
            std::map<std::string, MonitorElement *> _elements;      
};


#endif