Type1PFMET Class Reference

Computes the Type-1 corrections for pfMET. A specific version of the Type1MET class from the JetMETCorrections/Type1MET package. More...

#include <Type1PFMET.h>

Inheritance diagram for Type1PFMET:

Public Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)
	Type1PFMET (const edm::ParameterSet &)
	Type1PFMET ()
virtual	~Type1PFMET ()
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	run (const reco::METCollection &uncorMET, const JetCorrector &corrector, const reco::PFJetCollection &uncorJet, double jetPTthreshold, double jetEMfracLimit, double jetMufracLimit, reco::METCollection *corMET)

Private Attributes
std::string	correctorLabel
edm::InputTag	inputUncorJetsTag
std::string	inputUncorMetLabel
double	jetEMfracLimit
double	jetMufracLimit
double	jetPTthreshold

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

Computes the Type-1 corrections for pfMET. A specific version of the Type1MET class from the JetMETCorrections/Type1MET package.

Author

Michal Bluj

Date

February 2009

Definition at line 31 of file Type1PFMET.h.

Constructor & Destructor Documentation

Type1PFMET::Type1PFMET ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 14 of file Type1PFMET.cc.

References edm::ParameterSet::getParameter().

{
  inputUncorMetLabel  = iConfig.getParameter<std::string>("inputUncorMetLabel");
  inputUncorJetsTag = iConfig.getParameter<edm::InputTag>("inputUncorJetsTag");
  correctorLabel   = iConfig.getParameter<std::string>("corrector");
  jetPTthreshold      = iConfig.getParameter<double>("jetPTthreshold");
  jetEMfracLimit      = iConfig.getParameter<double>("jetEMfracLimit");
  jetMufracLimit      = iConfig.getParameter<double>("jetMufracLimit");
  produces<METCollection>();
}

Type1PFMET::Type1PFMET ( )

explicit

Definition at line 24 of file Type1PFMET.cc.

{}

Type1PFMET::~Type1PFMET ( )

virtual

Definition at line 27 of file Type1PFMET.cc.

{}

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 30 of file Type1PFMET.cc.

References edm::Event::getByLabel(), JetCorrector::getJetCorrector(), convertSQLitetoXML_cfg::output, edm::Event::put(), and DTTTrigCorrFirst::run.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

{
  using namespace edm;
  Handle<PFJetCollection> inputUncorJets;
  iEvent.getByLabel( inputUncorJetsTag, inputUncorJets );
  const JetCorrector* corrector = JetCorrector::getJetCorrector (correctorLabel, iSetup);
  Handle<METCollection> inputUncorMet;                     //Define Inputs
  iEvent.getByLabel( inputUncorMetLabel,  inputUncorMet );     //Get Inputs
  std::auto_ptr<METCollection> output( new METCollection() );  //Create empty output
  run( *(inputUncorMet.product()), *corrector, *(inputUncorJets.product()), 
       jetPTthreshold, jetEMfracLimit, jetMufracLimit,
       &*output );                                         //Invoke the algorithm
  iEvent.put( output );                                        //Put output into Event
}

void Type1PFMET::run ( const reco::METCollection &  uncorMET, const JetCorrector &  corrector, const reco::PFJetCollection &  uncorJet, double  jetPTthreshold, double  jetEMfracLimit, double  jetMufracLimit, reco::METCollection *  corMET  )

private

Definition at line 45 of file Type1PFMET.cc.

References dtNoiseDBValidation_cfg::cerr, corr, JetCorrector::correction(), delta, metsig::jet, reco::MET::mEtCorr(), CorrMETData::mex, CorrMETData::mey, reco::LeafCandidate::px(), reco::LeafCandidate::py(), query::result, mathSSE::sqrt(), CorrMETData::sumet, reco::MET::sumEt(), and reco::LeafCandidate::vertex().

Referenced by Types.LuminosityBlockID::cppID().

{
  if (!corMET) {
    std::cerr << "Type1METAlgo_run-> undefined output MET collection. Stop. " << std::endl;
    return;
  }

  double DeltaPx = 0.0;
  double DeltaPy = 0.0;
  double DeltaSumET = 0.0;
  // ---------------- Calculate jet corrections, but only for those uncorrected jets
  // ---------------- which are above the given threshold.  This requires that the
  // ---------------- uncorrected jets be matched with the corrected jets.
  for( PFJetCollection::const_iterator jet = uncorJet.begin(); jet != uncorJet.end(); ++jet) {
    if( jet->pt() > jetPTthreshold ) {
      double emEFrac = 
    jet->chargedEmEnergyFraction() + jet->neutralEmEnergyFraction();
      double muEFrac = jet->chargedMuEnergyFraction();
      if( emEFrac < jetEMfracLimit
      && muEFrac < jetMufracLimit ) {
    double corr = corrector.correction (*jet) - 1.; // correction itself
    DeltaPx +=  jet->px() * corr;
    DeltaPy +=  jet->py() * corr;
    DeltaSumET += jet->et() * corr;
      }
    }
  }
  //----------------- Calculate and set deltas for new MET correction
  CorrMETData delta;
  delta.mex   =  - DeltaPx;    //correction to MET (from Jets) is negative,    
  delta.mey   =  - DeltaPy;    //since MET points in direction opposite of jets
  delta.sumet =  DeltaSumET; 
  //----------------- Fill holder with corrected MET (= uncorrected + delta) values
  const MET* u = &(uncorMET.front());
  double corrMetPx = u->px()+delta.mex;
  double corrMetPy = u->py()+delta.mey;
  MET::LorentzVector correctedMET4vector( corrMetPx, corrMetPy, 0., 
                      sqrt (corrMetPx*corrMetPx + corrMetPy*corrMetPy)
                      );
  //----------------- get previous corrections and push into new corrections 
  std::vector<CorrMETData> corrections = u->mEtCorr();
  corrections.push_back( delta );
  //----------------- Push onto MET Collection
  MET result = MET(u->sumEt()+delta.sumet, 
           corrections, 
           correctedMET4vector, 
           u->vertex() );
  corMET->push_back(result);
  
  return;
}

Member Data Documentation

std::string Type1PFMET::correctorLabel

private

Definition at line 41 of file Type1PFMET.h.

edm::InputTag Type1PFMET::inputUncorJetsTag

private

Definition at line 40 of file Type1PFMET.h.

std::string Type1PFMET::inputUncorMetLabel

private

Definition at line 39 of file Type1PFMET.h.

double Type1PFMET::jetEMfracLimit

private

Definition at line 43 of file Type1PFMET.h.

double Type1PFMET::jetMufracLimit

private

Definition at line 44 of file Type1PFMET.h.

double Type1PFMET::jetPTthreshold

private

Definition at line 42 of file Type1PFMET.h.