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 29 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(), EgHLTOffHistBins_cfi::et, SpecificPFMETData::MuonFraction, SpecificPFMETData::NeutralEMFraction, SpecificPFMETData::NeutralHadFraction, reco::Candidate::pdgId(), funct::sin(), timingPdfMaker::specific, 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;
}