PFSpecificAlgo Class Reference

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

Public Member Functions
	PFSpecificAlgo ()
SpecificPFMETData	run (const edm::View< reco::Candidate > &pfCands, edm::ValueMap< float > const *weights=nullptr)

Detailed Description

Description: Adds Particle Flow specific information to MET

Implementation: [Notes on implementation]

Definition at line 28 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,
		edm::ValueMap< float > const *	weights = nullptr
	)

Definition at line 15 of file PFSpecificAlgo.cc.

References funct::abs(), SpecificPFMETData::ChargedEMFraction, SpecificPFMETData::ChargedHadFraction, alignCSCRings::e, edm::View< T >::empty(), reco::Candidate::energy(), reco::Candidate::get(), SpecificPFMETData::MuonFraction, SpecificPFMETData::NeutralEMFraction, SpecificPFMETData::NeutralHadFraction, reco::Candidate::pdgId(), edm::View< T >::ptrs(), funct::sin(), timingPdfMaker::specific, reco::Candidate::theta(), theta(), SpecificPFMETData::Type6Fraction, SpecificPFMETData::Type7Fraction, and histoStyle::weight.

Referenced by 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 (auto const& pfCandPtr : pfCands.ptrs()) {
    const reco::Candidate* pfCand = pfCandPtr.get();
    float weight = (weights != nullptr) ? (*weights)[pfCandPtr] : 1.0;
    if (!pfCand)
      continue;
    const double theta = pfCand->theta();
    const double e = pfCand->energy() * weight;
    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;
}