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.

{ }

Member Function Documentation

SpecificPFMETData PFSpecificAlgo::run

(

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

pfCands

)

Definition at line 15 of file PFSpecificAlgo.cc.

References edm::View< T >::begin(), SpecificPFMETData::ChargedEMFraction, SpecificPFMETData::ChargedHadFraction, alignCSCRings::e, edm::View< T >::end(), reco::LeafCandidate::energy(), SpecificPFMETData::MuonFraction, SpecificPFMETData::NeutralEMFraction, SpecificPFMETData::NeutralHadFraction, reco::PFCandidate::particleId(), funct::sin(), edm::View< T >::size(), timingPdfMaker::specific, reco::LeafCandidate::theta(), theta(), SpecificPFMETData::Type6Fraction, and SpecificPFMETData::Type7Fraction.

Referenced by cms::PFMETProducer::produce(), and cms::METProducer::produce_PFMET().

{
  if(!pfCands.size()) 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::PFCandidate* pfCand = dynamic_cast<const reco::PFCandidate*> (&(*iPfCand));
      if (!pfCand) continue;
      const double theta = pfCand->theta();
      const double e     = pfCand->energy();
      const double et    = e*sin(theta);

      if (pfCand->particleId() == 1) ChargedHadEt += et;
      if (pfCand->particleId() == 2) ChargedEMEt += et;
      if (pfCand->particleId() == 3) MuonEt += et;
      if (pfCand->particleId() == 4) NeutralEMEt += et;
      if (pfCand->particleId() == 5) NeutralHadEt += et;
      if (pfCand->particleId() == 6) type6Et += et;
      if (pfCand->particleId() == 7) type7Et += et;
    }

  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;
}