#include <CaloJetMETcorrInputProducerT.h>

Inheritance diagram for CaloJetMETcorrInputProducerT< T, Textractor >:

Public Member Functions
	CaloJetMETcorrInputProducerT (const edm::ParameterSet &cfg)

	~CaloJetMETcorrInputProducerT ()

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 &evt, const edm::EventSetup &es)

Private Attributes
double	jetCorrEtaMax_

Textractor	jetCorrExtractor_

std::string	jetCorrLabel_

std::string	moduleLabel_

std::string	offsetCorrLabel_

bool	skipEM_

double	skipEMfractionThreshold_

edm::InputTag	src_

edm::InputTag	srcMET_

double	type1JetPtThreshold_

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

template<typename T, typename Textractor>
class CaloJetMETcorrInputProducerT< T, Textractor >

Definition at line 60 of file CaloJetMETcorrInputProducerT.h.

Constructor & Destructor Documentation

template<typename T , typename Textractor >

CaloJetMETcorrInputProducerT< T, Textractor >::CaloJetMETcorrInputProducerT ( const edm::ParameterSet & cfg )

inlineexplicit

Definition at line 64 of file CaloJetMETcorrInputProducerT.h.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label"))
   {
     src_ = cfg.getParameter<edm::InputTag>("src");
 
     offsetCorrLabel_ = ( cfg.exists("offsetCorrLabel") ) ?
       cfg.getParameter<std::string>("offsetCorrLabel") : "";
     jetCorrLabel_ = cfg.getParameter<std::string>("jetCorrLabel");
 
     jetCorrEtaMax_ = ( cfg.exists("jetCorrEtaMax") ) ?
       cfg.getParameter<double>("jetCorrEtaMax") : 9.9;
 
     type1JetPtThreshold_ = cfg.getParameter<double>("type1JetPtThreshold");
 
     skipEM_ = cfg.getParameter<bool>("skipEM");
     if ( skipEM_ ) {
       skipEMfractionThreshold_ = cfg.getParameter<double>("skipEMfractionThreshold");
     }
 
     if ( cfg.exists("srcMET") ) {
       srcMET_ = cfg.getParameter<edm::InputTag>("srcMET");
     }
 
     produces<CorrMETData>("type1");
     produces<CorrMETData>("type2");
     produces<CorrMETData>("offset");
   }

template<typename T , typename Textractor >

CaloJetMETcorrInputProducerT< T, Textractor >::~CaloJetMETcorrInputProducerT ( )

inline

Definition at line 91 of file CaloJetMETcorrInputProducerT.h.

91 {}

Member Function Documentation

template<typename T , typename Textractor >

void CaloJetMETcorrInputProducerT< T, Textractor >::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

inlineprivatevirtual

Implements edm::EDProducer.

Definition at line 95 of file CaloJetMETcorrInputProducerT.h.

   {
     std::auto_ptr<CorrMETData> type1Correction(new CorrMETData());
     std::auto_ptr<CorrMETData> unclEnergySum(new CorrMETData());
     std::auto_ptr<CorrMETData> offsetEnergySum(new CorrMETData());
 
     typedef std::vector<T> JetCollection;
     edm::Handle<JetCollection> jets;
     evt.getByLabel(src_, jets);
 
     typedef edm::View<reco::MET> METView;
     edm::Handle<METView> met;
     if ( srcMET_.label() != "" ) {
       evt.getByLabel(srcMET_, met);
       if ( met->size() != 1 )
     throw cms::Exception("CaloJetMETcorrInputProducer::produce") 
       << "Failed to find unique MET in the event, src = " << srcMET_.label() << " !!\n";
 
 //--- compute "unclustered energy" by sutracting from the reconstructed MET 
 //   (i.e. from the negative vectorial sum of all particles reconstructed in the event)
 //    the momenta of (high Pt) jets which enter Type 1 MET corrections 
 //
 //    NOTE: MET = -(jets + muons + "unclustered energy"),
 //          so "unclustered energy" = -(MET + jets + muons),
 //          i.e. (high Pt) jets enter the sum of "unclustered energy" with negative sign.
 //
       unclEnergySum->mex   = -met->front().px();
       unclEnergySum->mey   = -met->front().py();
       unclEnergySum->sumet =  met->front().sumEt();
     }
 
     int numJets = jets->size();
     for ( int jetIndex = 0; jetIndex < numJets; ++jetIndex ) {
       const T& rawJet = jets->at(jetIndex);
 
       static CaloJetMETcorrInputProducer_namespace::InputTypeCheckerT<T> checkInputType;
       checkInputType(rawJet);
 
       static CaloJetMETcorrInputProducer_namespace::RawJetExtractorT<T> rawJetExtractor;
       reco::Candidate::LorentzVector rawJetP4 = rawJetExtractor(rawJet);
       
       reco::Candidate::LorentzVector corrJetP4 = jetCorrExtractor_(rawJet, jetCorrLabel_, &evt, &es, jetCorrEtaMax_);
 
       if ( corrJetP4.pt() > type1JetPtThreshold_ ) {
     
     unclEnergySum->mex   -= rawJetP4.px();
     unclEnergySum->mey   -= rawJetP4.py();
     unclEnergySum->sumet -= rawJetP4.Et();
     
     if ( skipEM_ && rawJet.emEnergyFraction() > skipEMfractionThreshold_ ) continue;
     
     reco::Candidate::LorentzVector rawJetP4offsetCorr = rawJetP4;
     if ( offsetCorrLabel_ != "" ) {
       rawJetP4offsetCorr = jetCorrExtractor_(rawJet, offsetCorrLabel_, &evt, &es, jetCorrEtaMax_);
       
       offsetEnergySum->mex   += (rawJetP4.px() - rawJetP4offsetCorr.px());
       offsetEnergySum->mey   += (rawJetP4.py() - rawJetP4offsetCorr.py());
       offsetEnergySum->sumet += (rawJetP4.Et() - rawJetP4offsetCorr.Et());
     }
 
 //--- MET balances momentum of reconstructed particles,
 //    hence correction to jets and corresponding Type 1 MET correction are of opposite sign
     type1Correction->mex   -= (corrJetP4.px() - rawJetP4offsetCorr.px());
     type1Correction->mey   -= (corrJetP4.py() - rawJetP4offsetCorr.py());
     type1Correction->sumet += (corrJetP4.Et() - rawJetP4offsetCorr.Et());
       } 
     }
 
 //--- add 
 //     o Type 1 MET correction                (difference corrected-uncorrected jet energy for jets of (corrected) Pt > 20 GeV)
 //     o momentum sum of "unclustered energy" (jets of (corrected) Pt < 20 GeV)
 //     o momentum sum of "offset energy"      (sum of energy attributed to pile-up/underlying event)
 //    to the event
     evt.put(type1Correction, "type1");
     evt.put(unclEnergySum,   "type2");
     evt.put(offsetEnergySum, "offset");
   }