PATMETProducer.h

Go to the documentation of this file.

//
//

#ifndef PhysicsTools_PatAlgos_PATMETProducer_h
#define PhysicsTools_PatAlgos_PATMETProducer_h

#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "CommonTools/Utils/interface/EtComparator.h"

#include "DataFormats/PatCandidates/interface/MET.h"
#include "PhysicsTools/PatAlgos/interface/EfficiencyLoader.h"
#include "PhysicsTools/PatAlgos/interface/KinResolutionsLoader.h"

#include "DataFormats/PatCandidates/interface/UserData.h"
#include "PhysicsTools/PatAlgos/interface/PATUserDataHelper.h"

#include "DataFormats/Candidate/interface/Candidate.h"
#include "RecoMET/METAlgorithms/interface/METSignificance.h"

namespace pat {

  class PATMETProducer : public edm::stream::EDProducer<> {
  public:
    explicit PATMETProducer(const edm::ParameterSet& iConfig);
    ~PATMETProducer() override;

    void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;

    static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

  private:
    // configurables
    edm::InputTag metSrc_;
    edm::EDGetTokenT<edm::View<reco::MET> > metToken_;
    bool addGenMET_;
    edm::EDGetTokenT<edm::View<reco::GenMET> > genMETToken_;
    bool addResolutions_;
    pat::helper::KinResolutionsLoader resolutionLoader_;
    bool addMuonCorr_;
    edm::InputTag muonSrc_;
    // tools
    GreaterByEt<MET> eTComparator_;

    bool addEfficiencies_;
    pat::helper::EfficiencyLoader efficiencyLoader_;

    bool useUserData_;
    pat::PATUserDataHelper<pat::MET> userDataHelper_;

    //MET Significance
    bool calculateMETSignificance_;
    metsig::METSignificance* metSigAlgo_;
    edm::EDGetTokenT<edm::View<reco::Jet> > jetToken_;
    edm::EDGetTokenT<edm::View<reco::Candidate> > pfCandToken_;
    std::vector<edm::EDGetTokenT<edm::View<reco::Candidate> > > lepTokens_;
    edm::EDGetTokenT<double> rhoToken_;
    std::string jetResPtType_;
    std::string jetResPhiType_;
    std::string jetSFType_;

    const reco::METCovMatrix getMETCovMatrix(const edm::Event& event,
                                             const edm::EventSetup& iSetup,
                                             double& sumPtUnclustered) const;
  };

}  // namespace pat

#endif