PFSpecificAlgo Class Reference

#include <RecoMET/METAlgorithms/interface/PFSpecificAlgo.h>

Public Member Functions
	PFSpecificAlgo ()
SpecificPFMETData	run (const edm::View< reco::Candidate > &pfCands)

Detailed Description

Description: Adds Particle Flow specific information to MET

Implementation: [Notes on implementation]

Definition at line 27 of file PFSpecificAlgo.h.

Constructor & Destructor Documentation

PFSpecificAlgo::PFSpecificAlgo ( )

inline

Definition at line 30 of file PFSpecificAlgo.h.

References run().

{ }

Member Function Documentation

SpecificPFMETData PFSpecificAlgo::run

(

const edm::View< reco::Candidate > &

pfCands

)

Definition at line 15 of file PFSpecificAlgo.cc.

References funct::abs(), edm::View< T >::begin(), SpecificPFMETData::ChargedEMFraction, SpecificPFMETData::ChargedHadFraction, MillePedeFileConverter_cfg::e, edm::View< T >::empty(), edm::View< T >::end(), reco::Candidate::energy(), stringResolutionProvider_cfi::et, SpecificPFMETData::MuonFraction, SpecificPFMETData::NeutralEMFraction, SpecificPFMETData::NeutralHadFraction, reco::Candidate::pdgId(), funct::sin(), reco::Candidate::theta(), theta(), SpecificPFMETData::Type6Fraction, and SpecificPFMETData::Type7Fraction.

Referenced by PFSpecificAlgo(), reco::PFMETProducerMVA::produce(), and cms::PFMETProducer::produce().

{
  if(pfCands.empty()) return SpecificPFMETData();

  double NeutralEMEt = 0.0;
  double NeutralHadEt = 0.0;
  double ChargedEMEt = 0.0;
  double ChargedHadEt = 0.0;
  double MuonEt = 0.0;
  double type6Et = 0.0;
  double type7Et = 0.0;

  for( edm::View<reco::Candidate>::const_iterator iPfCand = pfCands.begin(); iPfCand != pfCands.end(); ++iPfCand)
    {   
      const reco::Candidate* pfCand = dynamic_cast<const reco::Candidate*> (&(*iPfCand));
      if (!pfCand) continue;
      const double theta = pfCand->theta();
      const double e     = pfCand->energy();
      const double et    = e*sin(theta);
      switch (abs(pfCand->pdgId())) {
        case 211: ChargedHadEt += et; break;
        case  11: ChargedEMEt += et; break;
        case  13: MuonEt += et; break;
        case  22: NeutralEMEt += et; break;
        case 130: NeutralHadEt += et; break;
        case   1: type6Et += et; break;
        case   2: type7Et += et; break;
      }
    }

  const double Et_total = NeutralEMEt + NeutralHadEt + ChargedEMEt + ChargedHadEt + MuonEt + type6Et + type7Et;
  SpecificPFMETData specific;
  if (Et_total!=0.0)
  {
    specific.NeutralEMFraction = NeutralEMEt/Et_total;
    specific.NeutralHadFraction = NeutralHadEt/Et_total;
    specific.ChargedEMFraction = ChargedEMEt/Et_total;
    specific.ChargedHadFraction = ChargedHadEt/Et_total;
    specific.MuonFraction = MuonEt/Et_total;
    specific.Type6Fraction = type6Et/Et_total;
    specific.Type7Fraction = type7Et/Et_total;
  }
  return specific;
}