#include <Type0PFMETcorrInputProducer.h>

Inheritance diagram for Type0PFMETcorrInputProducer:

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

	~Type0PFMETcorrInputProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

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

Private Attributes
TFormula *	correction_

double	minDz_

std::string	moduleLabel_

edm::InputTag	srcHardScatterVertex_

edm::InputTag	srcPFCandidateToVertexAssociations_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

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

Version

Revision:: 1.1

Id:: Type0PFMETcorrInputProducer.h,v 1.1 2012/04/19 17:59:40 veelken Exp

Definition at line 27 of file Type0PFMETcorrInputProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file Type0PFMETcorrInputProducer.cc.

References correction_, edm::ParameterSet::getParameter(), minDz_, moduleLabel_, srcHardScatterVertex_, and srcPFCandidateToVertexAssociations_.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label")),
     correction_(0)
 {
   srcPFCandidateToVertexAssociations_ = cfg.getParameter<edm::InputTag>("srcPFCandidateToVertexAssociations");
   srcHardScatterVertex_ = 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 ( )

Definition at line 39 of file Type0PFMETcorrInputProducer.cc.

References correction_.

 {
   delete correction_;
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDProducer.

Definition at line 44 of file Type0PFMETcorrInputProducer.cc.

References correction_, reco::PFCandidate::e, edm::Event::getByLabel(), reco::PFCandidate::h, minDz_, reco::PFCandidate::mu, reco::LeafCandidate::p4(), reco::PFCandidate::particleId(), pfMETCorrectionType0_cfi::pfCandidateToVertexAssociation, phi, edm::Event::put(), srcHardScatterVertex_, and srcPFCandidateToVertexAssociations_.

 {
   edm::Handle<reco::VertexCollection> hardScatterVertex;
   evt.getByLabel(srcHardScatterVertex_, hardScatterVertex);
 
   edm::Handle<PFCandToVertexAssMap> pfCandidateToVertexAssociations;
   evt.getByLabel(srcPFCandidateToVertexAssociations_, pfCandidateToVertexAssociations);
 
   std::auto_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(pfMEtCorrection);
 }