#include <Type0PFMETcorrInputProducer.h>

Inheritance diagram for Type0PFMETcorrInputProducer:

Public Member Functions
	Type0PFMETcorrInputProducer (const edm::ParameterSet &)

	~Type0PFMETcorrInputProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
TFormula *	correction_

edm::EDGetTokenT< reco::VertexCollection >	hardScatterVertexToken_

double	minDz_

std::string	moduleLabel_

edm::EDGetTokenT< PFCandToVertexAssMap >	pfCandidateToVertexAssociationsToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Compute Type 0 (PF)MET corrections ( https://indico.cern.ch/getFile.py/access?contribId=4&resId=1&materialId=slides&confId=161159 )

Author: Christian Veelken, LLR

Definition at line 29 of file Type0PFMETcorrInputProducer.h.

Constructor & Destructor Documentation

Type0PFMETcorrInputProducer::Type0PFMETcorrInputProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 13 of file Type0PFMETcorrInputProducer.cc.

References correction_, edm::ParameterSet::getParameter(), hardScatterVertexToken_, minDz_, moduleLabel_, pfCandidateToVertexAssociationsToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label")),
     correction_(nullptr)
 {
   pfCandidateToVertexAssociationsToken_ = consumes<PFCandToVertexAssMap>(cfg.getParameter<edm::InputTag>("srcPFCandidateToVertexAssociations"));
   hardScatterVertexToken_ = consumes<reco::VertexCollection>(cfg.getParameter<edm::InputTag>("srcHardScatterVertex"));
 
   edm::ParameterSet cfgCorrection_function = cfg.getParameter<edm::ParameterSet>("correction");
   std::string corrFunctionName = std::string(moduleLabel_).append("correction");
   std::string corrFunctionFormula = cfgCorrection_function.getParameter<std::string>("formula");
   correction_ = new TFormula(corrFunctionName.data(), corrFunctionFormula.data());
   int numParameter = correction_->GetNpar();
   for ( int iParameter = 0; iParameter < numParameter; ++iParameter ) {
     std::string parName = Form("par%i", iParameter);
     double parValue = cfgCorrection_function.getParameter<double>(parName);
     correction_->SetParameter(iParameter, parValue);
   }
 
   minDz_ = cfg.getParameter<double>("minDz");
 
   produces<CorrMETData>();
 }

Type0PFMETcorrInputProducer::~Type0PFMETcorrInputProducer ( )

override

Definition at line 36 of file Type0PFMETcorrInputProducer.cc.

References correction_.

 {
   delete correction_;
 }

Member Function Documentation

void Type0PFMETcorrInputProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

overrideprivate

Definition at line 41 of file Type0PFMETcorrInputProducer.cc.

References Abs(), edm::AssociationMap< Tag >::begin(), correction_, DEFINE_FWK_MODULE, reco::PFCandidate::e, edm::AssociationMap< Tag >::end(), edm::Event::getByToken(), reco::PFCandidate::h, hardScatterVertexToken_, edm::Ref< C, T, F >::key(), minDz_, eostools::move(), reco::PFCandidate::mu, reco::LeafCandidate::p4(), reco::PFCandidate::particleId(), pfCandidateToVertexAssociationsToken_, phi, EnergyCorrector::pt, and edm::Event::put().

 {
   edm::Handle<reco::VertexCollection> hardScatterVertex;
   evt.getByToken(hardScatterVertexToken_, hardScatterVertex);
 
   edm::Handle<PFCandToVertexAssMap> pfCandidateToVertexAssociations;
   evt.getByToken(pfCandidateToVertexAssociationsToken_, pfCandidateToVertexAssociations);
 
   std::unique_ptr<CorrMETData> pfMEtCorrection(new CorrMETData());
   
   for ( PFCandToVertexAssMap::const_iterator pfCandidateToVertexAssociation = pfCandidateToVertexAssociations->begin();
     pfCandidateToVertexAssociation != pfCandidateToVertexAssociations->end(); ++pfCandidateToVertexAssociation ) {
     reco::VertexRef vertex = pfCandidateToVertexAssociation->key;
     const PFCandQualityPairVector& pfCandidates_vertex = pfCandidateToVertexAssociation->val;
     
     bool isHardScatterVertex = false;
     for ( reco::VertexCollection::const_iterator hardScatterVertex_i = hardScatterVertex->begin();
       hardScatterVertex_i != hardScatterVertex->end(); ++hardScatterVertex_i ) {
       if ( TMath::Abs(vertex->position().z() - hardScatterVertex_i->position().z()) < minDz_ ) {
     isHardScatterVertex = true;
     break;
       }
     }
   
     if ( !isHardScatterVertex ) {
       reco::Candidate::LorentzVector sumChargedPFCandP4_vertex;
       for ( PFCandQualityPairVector::const_iterator pfCandidate_vertex = pfCandidates_vertex.begin();
         pfCandidate_vertex != pfCandidates_vertex.end(); ++pfCandidate_vertex ) {
     const reco::PFCandidate& pfCandidate = (*pfCandidate_vertex->first);
     if ( pfCandidate.particleId() == reco::PFCandidate::h  ||
          pfCandidate.particleId() == reco::PFCandidate::e  ||
          pfCandidate.particleId() == reco::PFCandidate::mu ) {
       sumChargedPFCandP4_vertex += pfCandidate.p4();
     }
       }
       
       double pt = sumChargedPFCandP4_vertex.pt();
       double phi = sumChargedPFCandP4_vertex.phi();
       double ptCorr = correction_->Eval(pt);
       double pxCorr = TMath::Cos(phi)*ptCorr;
       double pyCorr = TMath::Sin(phi)*ptCorr;
 
       pfMEtCorrection->mex += pxCorr;
       pfMEtCorrection->mey += pyCorr;
     }
   }
 
   evt.put(std::move(pfMEtCorrection));
 }