PFMuonAlgo Class Reference

#include <PFMuonAlgo.h>

Classes
class	IndexPtComparator
class	TrackMETComparator
class	TrackPtErrorSorter

Public Member Functions
void	addMissingMuons (edm::Handle< reco::MuonCollection >, reco::PFCandidateCollection *cands)
void	changeTrack (reco::PFCandidate &, const MuonTrackTypePair &)
	PFMuonAlgo (edm::ParameterSet const &, bool postMuonCleaning)
	constructor More...
void	postClean (reco::PFCandidateCollection *)
bool	reconstructMuon (reco::PFCandidate &, const reco::MuonRef &, bool allowLoose=false)
void	setInputsForCleaning (reco::VertexCollection const &)
void	setVetoes (const reco::PFCandidateCollection &vetoes)
std::unique_ptr< reco::PFCandidateCollection >	transferAddedMuonCandidates ()
std::unique_ptr< reco::PFCandidateCollection >	transferCleanedCosmicCandidates ()
std::unique_ptr< reco::PFCandidateCollection >	transferCleanedFakeCandidates ()
std::unique_ptr< reco::PFCandidateCollection >	transferCleanedTrackerAndGlobalCandidates ()
std::unique_ptr< reco::PFCandidateCollection >	transferPunchThroughCleanedHadronCandidates ()
std::unique_ptr< reco::PFCandidateCollection >	transferPunchThroughCleanedMuonCandidates ()

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &iDesc)
static bool	hasValidTrack (const reco::MuonRef &muonRef, bool loose, double maxDPtOPt)
static bool	isGlobalLooseMuon (const reco::PFBlockElement &elt)
static bool	isGlobalLooseMuon (const reco::MuonRef &muonRef)
static bool	isGlobalTightMuon (const reco::PFBlockElement &elt)
static bool	isGlobalTightMuon (const reco::MuonRef &muonRef)
static bool	isIsolatedMuon (const reco::PFBlockElement &elt)
static bool	isIsolatedMuon (const reco::MuonRef &muonRef)
static bool	isLooseMuon (const reco::PFBlockElement &elt)
static bool	isLooseMuon (const reco::MuonRef &muonRef)
static bool	isMuon (const reco::PFBlockElement &elt)
static bool	isMuon (const reco::MuonRef &muonRef)
static bool	isTightMuonPOG (const reco::MuonRef &muonRef)
static bool	isTrackerLooseMuon (const reco::PFBlockElement &elt)
static bool	isTrackerLooseMuon (const reco::MuonRef &muonRef)
static bool	isTrackerTightMuon (const reco::PFBlockElement &elt)
static bool	isTrackerTightMuon (const reco::MuonRef &muonRef)
static int	muAssocToTrack (const reco::TrackRef &trackref, const reco::MuonCollection &muons)
static std::vector< reco::Muon::MuonTrackTypePair >	muonTracks (const reco::MuonRef &muon, double maxDPtOPt=1e+9, bool includeSA=false)
static void	printMuonProperties (const reco::MuonRef &muonRef)

Private Types
typedef reco::Muon::MuonTrackType	MuonTrackType
typedef reco::Muon::MuonTrackTypePair	MuonTrackTypePair

Private Member Functions
bool	cleanMismeasured (reco::PFCandidate &, unsigned int)
bool	cleanPunchThroughAndFakes (reco::PFCandidate &, reco::PFCandidateCollection *, unsigned int)
void	estimateEventQuantities (const reco::PFCandidateCollection *)
std::pair< double, double >	getMinMaxMET2 (const reco::PFCandidate &)
MuonTrackTypePair	getTrackWithSmallestError (const std::vector< MuonTrackTypePair > &)
void	removeDeadCandidates (reco::PFCandidateCollection *, const std::vector< unsigned int > &)
std::vector< MuonTrackTypePair >	tracksPointingAtMET (const std::vector< MuonTrackTypePair > &)
std::vector< MuonTrackTypePair >	tracksWithBetterMET (const std::vector< MuonTrackTypePair > &, const reco::PFCandidate &)

Private Attributes
const double	cosmicRejDistance_
const double	errorCompScale_
const double	eventFactorCosmics_
const double	eventFractionCleaning_
const double	eventFractionRejection_
std::vector< unsigned int >	maskedIndices_
const double	maxDPtOPt_
const double	metFactorCleaning_
const double	metFactorFake_
const double	metFactorHighEta_
const double	metFactorRejection_
const double	metSigForCleaning_
const double	metSigForRejection_
double	METX_
double	METY_
const double	minPostCleaningPt_
const double	minPunchThroughEnergy_
const double	minPunchThroughMomentum_
std::unique_ptr< reco::PFCandidateCollection >	pfAddedMuonCandidates_
	the collection of added muon candidates More...
std::unique_ptr< reco::PFCandidateCollection >	pfCleanedTrackerAndGlobalMuonCandidates_
	the collection of tracker/global cleaned muon candidates More...
std::unique_ptr< reco::PFCandidateCollection >	pfCosmicsMuonCleanedCandidates_
	the collection of cosmics cleaned muon candidates More...
std::unique_ptr< reco::PFCandidateCollection >	pfFakeMuonCleanedCandidates_
	the collection of fake cleaned muon candidates More...
std::unique_ptr< reco::PFCandidateCollection >	pfPunchThroughHadronCleanedCandidates_
	the collection of punch-through cleaned neutral hadron candidates More...
std::unique_ptr< reco::PFCandidateCollection >	pfPunchThroughMuonCleanedCandidates_
	the collection of punch-through cleaned muon candidates More...
const bool	postCleaning_
const double	ptFactorHighEta_
const double	punchThroughFactor_
const double	punchThroughMETFactor_
double	sumet_
double	sumetPU_
const reco::TrackBase::TrackQuality	trackQuality_
const reco::VertexCollection *	vertices_
const reco::PFCandidateCollection *	vetoes_ = nullptr

Detailed Description

Definition at line 13 of file PFMuonAlgo.h.

Member Typedef Documentation

MuonTrackType

typedef reco::Muon::MuonTrackType PFMuonAlgo::MuonTrackType

private

Definition at line 15 of file PFMuonAlgo.h.

MuonTrackTypePair

typedef reco::Muon::MuonTrackTypePair PFMuonAlgo::MuonTrackTypePair

private

Definition at line 14 of file PFMuonAlgo.h.

Constructor & Destructor Documentation

PFMuonAlgo()

PFMuonAlgo::PFMuonAlgo	(	edm::ParameterSet const &	iConfig,
		bool	postMuonCleaning
	)

constructor

Definition at line 17 of file PFMuonAlgo.cc.

    : pfCosmicsMuonCleanedCandidates_(std::make_unique<reco::PFCandidateCollection>()),
      pfCleanedTrackerAndGlobalMuonCandidates_(std::make_unique<reco::PFCandidateCollection>()),
      pfFakeMuonCleanedCandidates_(std::make_unique<reco::PFCandidateCollection>()),
      pfPunchThroughMuonCleanedCandidates_(std::make_unique<reco::PFCandidateCollection>()),
      pfPunchThroughHadronCleanedCandidates_(std::make_unique<reco::PFCandidateCollection>()),
      pfAddedMuonCandidates_(std::make_unique<reco::PFCandidateCollection>()),

      maxDPtOPt_(iConfig.getParameter<double>("maxDPtOPt")),
      trackQuality_(reco::TrackBase::qualityByName(iConfig.getParameter<std::string>("trackQuality"))),
      errorCompScale_(iConfig.getParameter<double>("ptErrorScale")),

      postCleaning_(postMuonCleaning),  // disable by default (for HLT)
      eventFractionCleaning_(iConfig.getParameter<double>("eventFractionForCleaning")),
      minPostCleaningPt_(iConfig.getParameter<double>("minPtForPostCleaning")),
      eventFactorCosmics_(iConfig.getParameter<double>("eventFactorForCosmics")),
      metSigForCleaning_(iConfig.getParameter<double>("metSignificanceForCleaning")),
      metSigForRejection_(iConfig.getParameter<double>("metSignificanceForRejection")),
      metFactorCleaning_(iConfig.getParameter<double>("metFactorForCleaning")),
      eventFractionRejection_(iConfig.getParameter<double>("eventFractionForRejection")),
      metFactorRejection_(iConfig.getParameter<double>("metFactorForRejection")),
      metFactorHighEta_(iConfig.getParameter<double>("metFactorForHighEta")),
      ptFactorHighEta_(iConfig.getParameter<double>("ptFactorForHighEta")),
      metFactorFake_(iConfig.getParameter<double>("metFactorForFakes")),
      minPunchThroughMomentum_(iConfig.getParameter<double>("minMomentumForPunchThrough")),
      minPunchThroughEnergy_(iConfig.getParameter<double>("minEnergyForPunchThrough")),
      punchThroughFactor_(iConfig.getParameter<double>("punchThroughFactor")),
      punchThroughMETFactor_(iConfig.getParameter<double>("punchThroughMETFactor")),
      cosmicRejDistance_(iConfig.getParameter<double>("cosmicRejectionDistance")) {}

Member Function Documentation

addMissingMuons()

void PFMuonAlgo::addMissingMuons	(	edm::Handle< reco::MuonCollection >	muons,
		reco::PFCandidateCollection *	cands
	)

Definition at line 724 of file PFMuonAlgo.cc.

References HLT_2022v12_cff::cands, changeTrack(), ALCARECOTkAlJpsiMuMu_cff::charge, cosmicRejDistance_, HLT_2022v12_cff::dEta, HLT_2022v12_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, mps_fire::i, edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), pfDeepBoostedJetPreprocessParams_cfi::lower_bound, maxDPtOPt_, METX_, METY_, reco::PFCandidate::mu, reco::PFCandidate::muonRef(), PDWG_BPHSkim_cff::muons, muonTracks(), pfAddedMuonCandidates_, pfElectronTranslator_cfi::PFCandidate, postCleaning_, jetUpdater_cfi::sort, mathSSE::sqrt(), reco::PFCandidate::trackRef(), tracks, tracksPointingAtMET(), vertices_, PbPb_ZMuSkimMuonDPG_cff::veto, and vetoes_.

                                                                                                      {
  if (!postCleaning_)
    return;

  if (pfAddedMuonCandidates_.get())
    pfAddedMuonCandidates_->clear();
  else
    pfAddedMuonCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  // do not recover muons if they are in the veto list
  typedef std::pair<edm::ProductID, unsigned> TrackProdIDKey;
  std::vector<TrackProdIDKey> vetoed;
  if (vetoes_) {
    for (const auto& veto : *vetoes_) {
      if (veto.trackRef().isNull())
        continue;
      vetoed.emplace_back(veto.trackRef().id(), veto.trackRef().key());
    }
    std::sort(vetoed.begin(), vetoed.end());
  }

  for (unsigned imu = 0; imu < muons->size(); ++imu) {
    reco::MuonRef muonRef(muons, imu);
    bool used = false;
    bool hadron = false;
    for (unsigned i = 0; i < cands->size(); i++) {
      const PFCandidate& pfc = cands->at(i);
      if (!pfc.trackRef().isNonnull())
        continue;
      if (pfc.trackRef().isNonnull() && pfc.trackRef() == muonRef->track())
        hadron = true;
      if (!pfc.muonRef().isNonnull())
        continue;

      if (pfc.muonRef()->innerTrack() == muonRef->innerTrack())
        used = true;
      else {
        // Check if the stand-alone muon is not a spurious copy of an existing muon
        // (Protection needed for HLT)
        if (pfc.muonRef()->isStandAloneMuon() && muonRef->isStandAloneMuon()) {
          double dEta = pfc.muonRef()->standAloneMuon()->eta() - muonRef->standAloneMuon()->eta();
          double dPhi = pfc.muonRef()->standAloneMuon()->phi() - muonRef->standAloneMuon()->phi();
          double dR = sqrt(dEta * dEta + dPhi * dPhi);
          if (dR < 0.005) {
            used = true;
          }
        }
      }

      if (used)
        break;
    }

    if (!vetoed.empty()) {
      TrackProdIDKey trk = std::make_pair(muonRef->track().id(), muonRef->track().key());
      auto lower = std::lower_bound(vetoed.begin(), vetoed.end(), trk);
      bool inVetoList = (lower != vetoed.end() && *lower == trk);
      if (inVetoList)
        used = true;
    }

    if (used || hadron || (!muonRef.isNonnull()))
      continue;

    TrackMETComparator comparator(METX_, METY_);
    //Low pt dont need to be cleaned

    std::vector<reco::Muon::MuonTrackTypePair> tracks = muonTracks(muonRef, maxDPtOPt_, true);
    //If there is at least 1 track choice  try to change the track
    if (!tracks.empty()) {
      //Find tracks that change dramatically MET or Pt
      std::vector<reco::Muon::MuonTrackTypePair> tracksThatChangeMET = tracksPointingAtMET(tracks);
      //From those tracks get the one with smallest MET
      if (!tracksThatChangeMET.empty()) {
        reco::Muon::MuonTrackTypePair bestTrackType =
            *std::min_element(tracksThatChangeMET.begin(), tracksThatChangeMET.end(), comparator);

        //Make sure it is not cosmic
        if ((vertices_->empty()) || bestTrackType.first->dz(vertices_->at(0).position()) < cosmicRejDistance_) {
          //make a pfcandidate
          int charge = bestTrackType.first->charge() > 0 ? 1 : -1;
          math::XYZTLorentzVector momentum(bestTrackType.first->px(),
                                           bestTrackType.first->py(),
                                           bestTrackType.first->pz(),
                                           sqrt(bestTrackType.first->p() * bestTrackType.first->p() + 0.1057 * 0.1057));

          cands->push_back(PFCandidate(charge, momentum, reco::PFCandidate::mu));

          changeTrack(cands->back(), bestTrackType);

          if (muonRef->track().isNonnull())
            cands->back().setTrackRef(muonRef->track());

          cands->back().setMuonRef(muonRef);

          pfAddedMuonCandidates_->push_back(cands->back());
        }
      }
    }
  }
}

changeTrack()

void PFMuonAlgo::changeTrack	(	reco::PFCandidate &	candidate,
		const MuonTrackTypePair &	track
	)

Definition at line 596 of file PFMuonAlgo.cc.

References HCALHighEnergyHPDFilter_cfi::energy, reco::PFCandidate::mu, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, reco::LeafCandidate::setCharge(), reco::PFCandidate::setMuonTrackType(), reco::LeafCandidate::setP4(), reco::PFCandidate::setParticleType(), reco::LeafCandidate::setVertex(), mathSSE::sqrt(), HLT_2022v12_cff::track, and PbPb_ZMuSkimMuonDPG_cff::trackType.

Referenced by addMissingMuons(), cleanMismeasured(), and reconstructMuon().

                                                                                       {
  using namespace reco;
  reco::TrackRef bestTrack = track.first;
  MuonTrackType trackType = track.second;
  //OK Now redefine the canddiate with that track
  double px = bestTrack->px();
  double py = bestTrack->py();
  double pz = bestTrack->pz();
  double energy = sqrt(bestTrack->p() * bestTrack->p() + 0.1057 * 0.1057);

  candidate.setCharge(bestTrack->charge() > 0 ? 1 : -1);
  candidate.setP4(math::XYZTLorentzVector(px, py, pz, energy));
  candidate.setParticleType(reco::PFCandidate::mu);
  //    candidate.setTrackRef( bestTrack );
  candidate.setMuonTrackType(trackType);
  candidate.setVertex(bestTrack->vertex());
}

cleanMismeasured()

bool PFMuonAlgo::cleanMismeasured ( reco::PFCandidate &  pfc, unsigned int  i  )

private

Definition at line 845 of file PFMuonAlgo.cc.

References changeTrack(), eventFractionRejection_, mps_fire::i, maskedIndices_, maxDPtOPt_, metFactorRejection_, metSigForRejection_, METX_, METY_, minPostCleaningPt_, reco::PFCandidate::muonRef(), muonTracks(), pfCleanedTrackerAndGlobalMuonCandidates_, pfFakeMuonCleanedCandidates_, reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), sumet_, sumetPU_, tracks, and tracksWithBetterMET().

Referenced by postClean().

                                                                      {
  using namespace std;
  using namespace reco;
  bool cleaned = false;

  //First define the MET without this guy
  double METNOX = METX_ - pfc.px();
  double METNOY = METY_ - pfc.py();
  double SUMETNO = sumet_ - pfc.pt();

  TrackMETComparator comparator(METNOX, METNOY);
  //Low pt dont need to be cleaned
  if (pfc.pt() < minPostCleaningPt_)
    return false;
  std::vector<reco::Muon::MuonTrackTypePair> tracks = muonTracks(pfc.muonRef(), maxDPtOPt_, false);

  //If there is more than 1 track choice  try to change the track
  if (tracks.size() > 1) {
    //Find tracks that change dramatically MET or Pt
    std::vector<reco::Muon::MuonTrackTypePair> tracksThatChangeMET = tracksWithBetterMET(tracks, pfc);
    //From those tracks get the one with smallest MET
    if (!tracksThatChangeMET.empty()) {
      reco::Muon::MuonTrackTypePair bestTrackType =
          *std::min_element(tracksThatChangeMET.begin(), tracksThatChangeMET.end(), comparator);
      changeTrack(pfc, bestTrackType);

      pfCleanedTrackerAndGlobalMuonCandidates_->push_back(pfc);
      //update eventquantities
      METX_ = METNOX + pfc.px();
      METY_ = METNOY + pfc.py();
      sumet_ = SUMETNO + pfc.pt();
    }
  }

  //Now attempt to kill it
  if (!(pfc.muonRef()->isGlobalMuon() && pfc.muonRef()->isTrackerMuon())) {
    //define MET significance and SUM ET
    double MET2 = METX_ * METX_ + METY_ * METY_;
    double newMET2 = METNOX * METNOX + METNOY * METNOY;
    double METSig = sqrt(MET2) / sqrt(sumet_ - sumetPU_);
    if (METSig > metSigForRejection_)
      if ((newMET2 < MET2 / metFactorRejection_) &&
          ((SUMETNO - sumetPU_) / (sumet_ - sumetPU_) < eventFractionRejection_)) {
        pfFakeMuonCleanedCandidates_->push_back(pfc);
        maskedIndices_.push_back(i);
        METX_ = METNOX;
        METY_ = METNOY;
        sumet_ = SUMETNO;
        cleaned = true;
      }
  }
  return cleaned;
}

cleanPunchThroughAndFakes()

bool PFMuonAlgo::cleanPunchThroughAndFakes ( reco::PFCandidate &  pfc, reco::PFCandidateCollection *  cands, unsigned int  imu  )

private

Definition at line 943 of file PFMuonAlgo.cc.

References HLT_2022v12_cff::cands, reco::PFCandidate::elementsInBlocks(), reco::LeafCandidate::eta(), DQMScaleToClient_cfi::factor, getMinMaxMET2(), reco::PFCandidate::h, reco::PFCandidate::h0, mps_fire::i, isIsolatedMuon(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), maskedIndices_, SiStripPI::max, metFactorFake_, metFactorHighEta_, METX_, METY_, minPostCleaningPt_, minPunchThroughEnergy_, minPunchThroughMomentum_, reco::PFCandidate::muonRef(), reco::LeafCandidate::p(), reco::PFCandidate::particleId(), pfFakeMuonCleanedCandidates_, pfPunchThroughHadronCleanedCandidates_, pfPunchThroughMuonCleanedCandidates_, reco::LeafCandidate::pt(), ptFactorHighEta_, punchThroughFactor_, punchThroughMETFactor_, reco::PFCandidate::rawEcalEnergy(), reco::PFCandidate::rawHcalEnergy(), edm::second(), sumet_, and sumetPU_.

Referenced by postClean().

                                                             {
  using namespace reco;

  bool cleaned = false;

  if (pfc.pt() < minPostCleaningPt_)
    return false;

  double METXNO = METX_ - pfc.pt();
  double METYNO = METY_ - pfc.pt();
  double MET2NO = METXNO * METXNO + METYNO * METYNO;
  double MET2 = METX_ * METX_ + METY_ * METY_;
  bool fake1 = false;

  std::pair<double, double> met2 = getMinMaxMET2(pfc);

  //Check for Fakes at high pseudorapidity
  if (pfc.muonRef()->standAloneMuon().isNonnull())
    fake1 = fabs(pfc.eta()) > 2.15 && met2.first < met2.second / 2 && MET2NO < MET2 / metFactorHighEta_ &&
            pfc.muonRef()->standAloneMuon()->pt() < pfc.pt() / ptFactorHighEta_;

  double factor = std::max(2., 2000. / (sumet_ - pfc.pt() - sumetPU_));
  bool fake2 =
      (pfc.pt() / (sumet_ - sumetPU_) < 0.25 && MET2NO < MET2 / metFactorFake_ && met2.first < met2.second / factor);

  bool punchthrough = pfc.p() > minPunchThroughMomentum_ && pfc.rawHcalEnergy() > minPunchThroughEnergy_ &&
                      pfc.rawEcalEnergy() + pfc.rawHcalEnergy() > pfc.p() / punchThroughFactor_ &&
                      !isIsolatedMuon(pfc.muonRef()) && MET2NO < MET2 / punchThroughMETFactor_;

  if (fake1 || fake2 || punchthrough) {
    // Find the block of the muon
    const PFCandidate::ElementsInBlocks& eleInBlocks = pfc.elementsInBlocks();
    if (!eleInBlocks.empty()) {
      PFBlockRef blockRefMuon = eleInBlocks[0].first;
      unsigned indexMuon = eleInBlocks[0].second;
      if (eleInBlocks.size() > 1)
        indexMuon = eleInBlocks[1].second;

      // Check if the muon gave rise to a neutral hadron
      double iHad = 1E9;
      bool hadron = false;
      for (unsigned i = imu + 1; i < cands->size(); ++i) {
        const PFCandidate& pfcn = cands->at(i);
        const PFCandidate::ElementsInBlocks& ele = pfcn.elementsInBlocks();
        if (ele.empty()) {
          continue;
        }
        PFBlockRef blockRefHadron = ele[0].first;
        unsigned indexHadron = ele[0].second;
        // We are out of the block -> exit the loop
        if (blockRefHadron.key() != blockRefMuon.key())
          break;
        // Check that this particle is a neutral hadron
        if (indexHadron == indexMuon && pfcn.particleId() == reco::PFCandidate::h0) {
          iHad = i;
          hadron = true;
        }
        if (hadron)
          break;
      }

      if (hadron) {
        double rescaleFactor = cands->at(iHad).p() / cands->at(imu).p();
        METX_ -= cands->at(imu).px() + cands->at(iHad).px();
        METY_ -= cands->at(imu).py() + cands->at(iHad).py();
        sumet_ -= cands->at(imu).pt();
        cands->at(imu).rescaleMomentum(rescaleFactor);
        maskedIndices_.push_back(iHad);
        pfPunchThroughHadronCleanedCandidates_->push_back(cands->at(iHad));
        cands->at(imu).setParticleType(reco::PFCandidate::h);
        pfPunchThroughMuonCleanedCandidates_->push_back(cands->at(imu));
        METX_ += cands->at(imu).px();
        METY_ += cands->at(imu).py();
        sumet_ += cands->at(imu).pt();

      } else if (fake1 || fake2) {
        METX_ -= cands->at(imu).px();
        METY_ -= cands->at(imu).py();
        sumet_ -= cands->at(imu).pt();
        maskedIndices_.push_back(imu);
        pfFakeMuonCleanedCandidates_->push_back(cands->at(imu));
        cleaned = true;
      }
    }
  }
  return cleaned;
}

estimateEventQuantities()

void PFMuonAlgo::estimateEventQuantities ( const reco::PFCandidateCollection *  pfc )

private

Definition at line 620 of file PFMuonAlgo.cc.

References mps_fire::i, METX_, METY_, sumet_, and sumetPU_.

Referenced by postClean().

                                                                             {
  //SUM ET
  sumetPU_ = 0.0;
  METX_ = 0.;
  METY_ = 0.;
  sumet_ = 0.0;
  for (reco::PFCandidateCollection::const_iterator i = pfc->begin(); i != pfc->end(); ++i) {
    sumet_ += i->pt();
    METX_ += i->px();
    METY_ += i->py();
  }
}

fillPSetDescription()

void PFMuonAlgo::fillPSetDescription ( edm::ParameterSetDescription &  iDesc )

static

Definition at line 1043 of file PFMuonAlgo.cc.

References edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PFTICLProducer::fillDescriptions(), and PFProducer::fillDescriptions().

                                                                      {
  // Muon ID and post cleaning parameters
  iDesc.add<double>("maxDPtOPt", 1.0);
  iDesc.add<std::string>("trackQuality", "highPurity");
  iDesc.add<double>("ptErrorScale", 8.0);

  iDesc.add<double>("eventFractionForCleaning", 0.5);
  iDesc.add<double>("minPtForPostCleaning", 20.0);
  iDesc.add<double>("eventFactorForCosmics", 10.0);
  iDesc.add<double>("metSignificanceForCleaning", 3.0);
  iDesc.add<double>("metSignificanceForRejection", 4.0);
  iDesc.add<double>("metFactorForCleaning", 4.0);
  iDesc.add<double>("eventFractionForRejection", 0.8);
  iDesc.add<double>("metFactorForRejection", 4.0);
  iDesc.add<double>("metFactorForHighEta", 25.0);
  iDesc.add<double>("ptFactorForHighEta", 2.0);
  iDesc.add<double>("metFactorForFakes", 4.0);
  iDesc.add<double>("minMomentumForPunchThrough", 100.0);
  iDesc.add<double>("minEnergyForPunchThrough", 100.0);
  iDesc.add<double>("punchThroughFactor", 3.0);
  iDesc.add<double>("punchThroughMETFactor", 4.0);
  iDesc.add<double>("cosmicRejectionDistance", 1.0);
}

getMinMaxMET2()

std::pair< double, double > PFMuonAlgo::getMinMaxMET2 ( const reco::PFCandidate &  pfc )

private

Definition at line 826 of file PFMuonAlgo.cc.

References mps_fire::i, maxDPtOPt_, METX_, METY_, reco::PFCandidate::muonRef(), muonTracks(), funct::pow(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and tracks.

Referenced by cleanPunchThroughAndFakes().

                                                                            {
  std::vector<reco::Muon::MuonTrackTypePair> tracks = muonTracks((pfc.muonRef()), maxDPtOPt_, true);

  double METXNO = METX_ - pfc.px();
  double METYNO = METY_ - pfc.py();
  std::vector<double> met2;
  for (unsigned int i = 0; i < tracks.size(); ++i) {
    met2.push_back(pow(METXNO + tracks.at(i).first->px(), 2) + pow(METYNO + tracks.at(i).first->py(), 2));
  }

  //PROTECT for cases of only one track. If there is only one track it will crash .
  //Has never happened but could likely happen!

  if (tracks.size() > 1)
    return std::make_pair(*std::min_element(met2.begin(), met2.end()), *std::max_element(met2.begin(), met2.end()));
  else
    return std::make_pair(0, 0);
}

getTrackWithSmallestError()

reco::Muon::MuonTrackTypePair PFMuonAlgo::getTrackWithSmallestError ( const std::vector< MuonTrackTypePair > &  tracks )

private

Definition at line 614 of file PFMuonAlgo.cc.

References MCScenario_CRAFT1_22X::sorter(), and tracks.

Referenced by reconstructMuon().

                                                        {
  TrackPtErrorSorter sorter;
  return *std::min_element(tracks.begin(), tracks.end(), sorter);
}

hasValidTrack()

bool PFMuonAlgo::hasValidTrack ( const reco::MuonRef &  muonRef, bool  loose, double  maxDPtOPt  )

static

Definition at line 348 of file PFMuonAlgo.cc.

References loose, particleFlowBlock_cfi::maxDPtOPt, and muonTracks().

Referenced by GeneralTracksImporter::importToBlock().

                                                                                       {
  if (loose)
    return !muonTracks(muonRef).empty();
  else
    return !muonTracks(muonRef, maxDPtOPt).empty();
}

isGlobalLooseMuon() [1/2]

bool PFMuonAlgo::isGlobalLooseMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 76 of file PFMuonAlgo.cc.

References cms::cuda::assert(), and reco::PFBlockElementTrack::muonRef().

Referenced by isLooseMuon().

                                                                {
  const reco::PFBlockElementTrack* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);
  reco::MuonRef muonRef = eltTrack->muonRef();

  return isGlobalLooseMuon(muonRef);
}

isGlobalLooseMuon() [2/2]

bool PFMuonAlgo::isGlobalLooseMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 197 of file PFMuonAlgo.cc.

References muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), SiStripPI::min, TrackCollections2monitor_cff::normalizedChi2, quality, mps_fire::result, and muon::TMLastStationAngTight.

                                                             {
  if (!muonRef.isNonnull())
    return false;
  if (!muonRef->isGlobalMuon())
    return false;
  if (!muonRef->isStandAloneMuon())
    return false;

  reco::TrackRef standAloneMu = muonRef->standAloneMuon();
  reco::TrackRef combinedMu = muonRef->combinedMuon();
  reco::TrackRef trackerMu = muonRef->track();

  unsigned nMuonHits =
      standAloneMu->hitPattern().numberOfValidMuonDTHits() + 2 * standAloneMu->hitPattern().numberOfValidMuonCSCHits();

  bool quality = false;

  if (muonRef->isTrackerMuon()) {
    bool result = combinedMu->normalizedChi2() < 100.;

    bool laststation = muon::isGoodMuon(*muonRef, muon::TMLastStationAngTight);

    int nMatches = muonRef->numberOfMatches();

    quality = laststation && nMuonHits > 12 && nMatches > 1;

    return result && quality;

  } else {
    // Check the quality of the stand-alone muon :
    // good chi**2 and large number of hits and good pt error
    if (nMuonHits <= 15 || standAloneMu->normalizedChi2() > 10. ||
        standAloneMu->ptError() / standAloneMu->pt() > 0.20) {
      quality = false;
    } else {
      // If the stand-alone muon is good, check the global muon
      if (combinedMu->normalizedChi2() > standAloneMu->normalizedChi2()) {
        // If the combined muon is worse than the stand-alone, it
        // means that either the corresponding tracker track was not
        // reconstructed, or that the sta muon comes from a late
        // pion decay (hence with a momentum smaller than the track)
        // Take the stand-alone muon only if its momentum is larger
        // than that of the track

        // Note that here we even take the standAlone if it has a smaller pT, in contrast to GlobalTight
        if (standAloneMu->pt() > trackerMu->pt() || combinedMu->normalizedChi2() < 5.)
          quality = true;
      } else {
        // If the combined muon is better (and good enough), take the
        // global muon
        if (combinedMu->ptError() / combinedMu->pt() < std::min(0.20, standAloneMu->ptError() / standAloneMu->pt()))
          quality = true;
      }
    }
  }

  return quality;
}

isGlobalTightMuon() [1/2]

bool PFMuonAlgo::isGlobalTightMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 67 of file PFMuonAlgo.cc.

References cms::cuda::assert(), and reco::PFBlockElementTrack::muonRef().

Referenced by isMuon().

                                                                {
  const reco::PFBlockElementTrack* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);
  reco::MuonRef muonRef = eltTrack->muonRef();

  return isGlobalTightMuon(muonRef);
}

isGlobalTightMuon() [2/2]

bool PFMuonAlgo::isGlobalTightMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 111 of file PFMuonAlgo.cc.

References muon::GlobalMuonPromptTight, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), SiStripPI::min, quality, mps_fire::result, and muon::TM2DCompatibilityTight.

                                                             {
  if (!muonRef.isNonnull())
    return false;

  if (!muonRef->isGlobalMuon())
    return false;
  if (!muonRef->isStandAloneMuon())
    return false;

  if (muonRef->isTrackerMuon()) {
    bool result = muon::isGoodMuon(*muonRef, muon::GlobalMuonPromptTight);

    bool isTM2DCompatibilityTight = muon::isGoodMuon(*muonRef, muon::TM2DCompatibilityTight);
    int nMatches = muonRef->numberOfMatches();
    bool quality = nMatches > 2 || isTM2DCompatibilityTight;

    return result && quality;

  } else {
    reco::TrackRef standAloneMu = muonRef->standAloneMuon();

    // No tracker muon -> Request a perfect stand-alone muon, or an even better global muon
    bool result = false;

    // Check the quality of the stand-alone muon :
    // good chi**2 and large number of hits and good pt error
    if ((standAloneMu->hitPattern().numberOfValidMuonDTHits() < 22 &&
         standAloneMu->hitPattern().numberOfValidMuonCSCHits() < 15) ||
        standAloneMu->normalizedChi2() > 10. || standAloneMu->ptError() / standAloneMu->pt() > 0.20) {
      result = false;
    } else {
      reco::TrackRef combinedMu = muonRef->combinedMuon();
      reco::TrackRef trackerMu = muonRef->track();

      // If the stand-alone muon is good, check the global muon
      if (combinedMu->normalizedChi2() > standAloneMu->normalizedChi2()) {
        // If the combined muon is worse than the stand-alone, it
        // means that either the corresponding tracker track was not
        // reconstructed, or that the sta muon comes from a late
        // pion decay (hence with a momentum smaller than the track)
        // Take the stand-alone muon only if its momentum is larger
        // than that of the track
        result = standAloneMu->pt() > trackerMu->pt();
      } else {
        // If the combined muon is better (and good enough), take the
        // global muon
        result =
            combinedMu->ptError() / combinedMu->pt() < std::min(0.20, standAloneMu->ptError() / standAloneMu->pt());
      }
    }

    return result;
  }

  return false;
}

isIsolatedMuon() [1/2]

bool PFMuonAlgo::isIsolatedMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 94 of file PFMuonAlgo.cc.

References cms::cuda::assert(), and reco::PFBlockElementTrack::muonRef().

Referenced by cleanPunchThroughAndFakes(), PFAlgo::createCandidatesHCAL(), GreedyMuonPFCandidateFilter::filter(), and isMuon().

                                                             {
  const reco::PFBlockElementTrack* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);
  reco::MuonRef muonRef = eltTrack->muonRef();

  return isIsolatedMuon(muonRef);
}

isIsolatedMuon() [2/2]

bool PFMuonAlgo::isIsolatedMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 279 of file PFMuonAlgo.cc.

References edm::Ref< C, T, F >::isNonnull(), and objects.autophobj::relIso.

                                                          {
  if (!muonRef.isNonnull())
    return false;
  if (!muonRef->isIsolationValid())
    return false;

  // Isolated Muons which are missed by standard cuts are nearly always global+tracker
  if (!muonRef->isGlobalMuon())
    return false;

  // If it's not a tracker muon, only take it if there are valid muon hits

  reco::TrackRef standAloneMu = muonRef->standAloneMuon();

  if (!muonRef->isTrackerMuon()) {
    if (standAloneMu->hitPattern().numberOfValidMuonDTHits() == 0 &&
        standAloneMu->hitPattern().numberOfValidMuonCSCHits() == 0)
      return false;
  }

  // for isolation, take the smallest pt available to reject fakes

  reco::TrackRef combinedMu = muonRef->combinedMuon();
  double smallestMuPt = combinedMu->pt();

  if (standAloneMu->pt() < smallestMuPt)
    smallestMuPt = standAloneMu->pt();

  if (muonRef->isTrackerMuon()) {
    reco::TrackRef trackerMu = muonRef->track();
    if (trackerMu->pt() < smallestMuPt)
      smallestMuPt = trackerMu->pt();
  }

  double sumPtR03 = muonRef->isolationR03().sumPt;

  double relIso = sumPtR03 / smallestMuPt;

  if (relIso < 0.1)
    return true;
  else
    return false;
}

isLooseMuon() [1/2]

bool PFMuonAlgo::isLooseMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 57 of file PFMuonAlgo.cc.

References cms::cuda::assert(), and reco::PFBlockElementTrack::muonRef().

Referenced by PFAlgo::createCandidatesHCAL(), GeneralTracksImporter::importToBlock(), and reconstructMuon().

                                                          {
  const reco::PFBlockElementTrack* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);

  reco::MuonRef muonRef = eltTrack->muonRef();

  return isLooseMuon(muonRef);
}

isLooseMuon() [2/2]

bool PFMuonAlgo::isLooseMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 107 of file PFMuonAlgo.cc.

References isGlobalLooseMuon(), and isTrackerLooseMuon().

                                                       {
  return (isGlobalLooseMuon(muonRef) || isTrackerLooseMuon(muonRef));
}

isMuon() [1/2]

bool PFMuonAlgo::isMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 48 of file PFMuonAlgo.cc.

References cms::cuda::assert().

Referenced by PFAlgo::createCandidatesHCAL(), GeneralTracksImporter::importToBlock(), PFEGammaAlgo::isMuon(), reconstructMuon(), and PFAlgo::recoTracksNotHCAL().

                                                     {
  const auto* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);
  reco::MuonRef muonRef = eltTrack->muonRef();

  return isMuon(muonRef);
}

isMuon() [2/2]

bool PFMuonAlgo::isMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 103 of file PFMuonAlgo.cc.

References isGlobalTightMuon(), isIsolatedMuon(), and isTrackerTightMuon().

                                                  {
  return isGlobalTightMuon(muonRef) || isTrackerTightMuon(muonRef) || isIsolatedMuon(muonRef);
}

isTightMuonPOG()

bool PFMuonAlgo::isTightMuonPOG ( const reco::MuonRef &  muonRef )

static

Definition at line 323 of file PFMuonAlgo.cc.

References muon::GlobalMuonPromptTight, and muon::isGoodMuon().

                                                          {
  if (!muon::isGoodMuon(*muonRef, muon::GlobalMuonPromptTight))
    return false;

  if (!muonRef->isTrackerMuon())
    return false;

  if (muonRef->numberOfMatches() < 2)
    return false;

  //const reco::TrackRef& combinedMuon = muonRef->combinedMuon();
  const reco::TrackRef& combinedMuon = muonRef->globalTrack();

  if (combinedMuon->hitPattern().numberOfValidTrackerHits() < 11)
    return false;

  if (combinedMuon->hitPattern().numberOfValidPixelHits() == 0)
    return false;

  if (combinedMuon->hitPattern().numberOfValidMuonHits() == 0)
    return false;

  return true;
}

isTrackerLooseMuon() [1/2]

static bool PFMuonAlgo::isTrackerLooseMuon ( const reco::PFBlockElement &  elt )

static

Referenced by isLooseMuon().

isTrackerLooseMuon() [2/2]

bool PFMuonAlgo::isTrackerLooseMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 256 of file PFMuonAlgo.cc.

References muon::AllArbitrated, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), quality, and muon::TMLastStationAngTight.

                                                              {
  if (!muonRef.isNonnull())
    return false;
  if (!muonRef->isTrackerMuon())
    return false;

  reco::TrackRef trackerMu = muonRef->track();

  if (trackerMu->ptError() / trackerMu->pt() > 0.20)
    return false;

  // this doesn't seem to be necessary on the small samples looked at, but keep it around as insurance
  if (trackerMu->pt() > 20.)
    return false;

  bool isAllArbitrated = muon::isGoodMuon(*muonRef, muon::AllArbitrated);
  bool isTMLastStationAngTight = muon::isGoodMuon(*muonRef, muon::TMLastStationAngTight);

  bool quality = isAllArbitrated && isTMLastStationAngTight;

  return quality;
}

isTrackerTightMuon() [1/2]

bool PFMuonAlgo::isTrackerTightMuon ( const reco::PFBlockElement &  elt )

static

Definition at line 85 of file PFMuonAlgo.cc.

References cms::cuda::assert(), and reco::PFBlockElementTrack::muonRef().

Referenced by isMuon().

                                                                 {
  const reco::PFBlockElementTrack* eltTrack = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);

  assert(eltTrack);
  reco::MuonRef muonRef = eltTrack->muonRef();

  return isTrackerTightMuon(muonRef);
}

isTrackerTightMuon() [2/2]

bool PFMuonAlgo::isTrackerTightMuon ( const reco::MuonRef &  muonRef )

static

Definition at line 168 of file PFMuonAlgo.cc.

References muon::AllArbitrated, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), muon::TM2DCompatibilityTight, and HLT_2022v12_cff::track.

                                                              {
  if (!muonRef.isNonnull())
    return false;

  if (!muonRef->isTrackerMuon())
    return false;

  reco::TrackRef trackerMu = muonRef->track();
  const reco::Track& track = *trackerMu;

  unsigned nTrackerHits = track.hitPattern().numberOfValidTrackerHits();

  if (nTrackerHits <= 12)
    return false;

  bool isAllArbitrated = muon::isGoodMuon(*muonRef, muon::AllArbitrated);

  bool isTM2DCompatibilityTight = muon::isGoodMuon(*muonRef, muon::TM2DCompatibilityTight);

  if (!isAllArbitrated || !isTM2DCompatibilityTight)
    return false;

  if ((trackerMu->ptError() / trackerMu->pt() > 0.10)) {
    //std::cout<<" PT ERROR > 10 % "<< trackerMu->pt() <<std::endl;
    return false;
  }
  return true;
}

muAssocToTrack()

int PFMuonAlgo::muAssocToTrack ( const reco::TrackRef &  trackref, const reco::MuonCollection &  muons  )

static

Definition at line 479 of file PFMuonAlgo.cc.

References HLT_2022v12_cff::distance, visualization-live-secondInstance_cfg::m, HLT_2022v12_cff::muon, and PDWG_BPHSkim_cff::muons.

Referenced by GeneralTracksImporter::importToBlock(), and PFTICLProducer::produce().

                                                                                            {
  auto muon = std::find_if(muons.cbegin(), muons.cend(), [&](const reco::Muon& m) {
    return (m.track().isNonnull() && m.track() == trackref);
  });
  return (muon != muons.cend() ? std::distance(muons.cbegin(), muon) : -1);
}

muonTracks()

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::muonTracks ( const reco::MuonRef &  muon, double  maxDPtOPt = 1e+9, bool  includeSA = false  )

static

Definition at line 486 of file PFMuonAlgo.cc.

References reco::Muon::CombinedTrack, reco::Muon::DYT, reco::Muon::InnerTrack, particleFlowBlock_cfi::maxDPtOPt, MillePedeFileConverter_cfg::out, reco::Muon::OuterTrack, reco::Muon::Picky, and reco::Muon::TPFMS.

Referenced by addMissingMuons(), cleanMismeasured(), getMinMaxMET2(), hasValidTrack(), and reconstructMuon().

                                                                                  {
  std::vector<reco::Muon::MuonTrackTypePair> out;

  if (muon->globalTrack().isNonnull() && muon->globalTrack()->pt() > 0)
    if (muon->globalTrack()->ptError() / muon->globalTrack()->pt() < maxDPtOPt)
      out.emplace_back(muon->globalTrack(), reco::Muon::CombinedTrack);

  if (muon->innerTrack().isNonnull() && muon->innerTrack()->pt() > 0)
    if (muon->innerTrack()->ptError() / muon->innerTrack()->pt() < maxDPtOPt)  //Here Loose!@
      out.emplace_back(muon->innerTrack(), reco::Muon::InnerTrack);

  bool pickyExists = false;
  double pickyDpt = 99999.;
  if (muon->pickyTrack().isNonnull() && muon->pickyTrack()->pt() > 0) {
    pickyDpt = muon->pickyTrack()->ptError() / muon->pickyTrack()->pt();
    if (pickyDpt < maxDPtOPt)
      out.emplace_back(muon->pickyTrack(), reco::Muon::Picky);
    pickyExists = true;
  }

  bool dytExists = false;
  double dytDpt = 99999.;
  if (muon->dytTrack().isNonnull() && muon->dytTrack()->pt() > 0) {
    dytDpt = muon->dytTrack()->ptError() / muon->dytTrack()->pt();
    if (dytDpt < maxDPtOPt)
      out.emplace_back(muon->dytTrack(), reco::Muon::DYT);
    dytExists = true;
  }

  //Magic: TPFMS is not a really good track especially under misalignment
  //IT is kind of crap because if mu system is displaced it can make a change
  //So allow TPFMS if there is no picky or the error of tpfms is better than picky
  //AND if there is no DYT or the error of tpfms is better than DYT
  if (muon->tpfmsTrack().isNonnull() && muon->tpfmsTrack()->pt() > 0) {
    double tpfmsDpt = muon->tpfmsTrack()->ptError() / muon->tpfmsTrack()->pt();
    if (((pickyExists && tpfmsDpt < pickyDpt) || (!pickyExists)) &&
        ((dytExists && tpfmsDpt < dytDpt) || (!dytExists)) && tpfmsDpt < maxDPtOPt)
      out.emplace_back(muon->tpfmsTrack(), reco::Muon::TPFMS);
  }

  if (includeSA && muon->outerTrack().isNonnull())
    if (muon->outerTrack()->ptError() / muon->outerTrack()->pt() < maxDPtOPt)
      out.emplace_back(muon->outerTrack(), reco::Muon::OuterTrack);

  return out;
}

postClean()

void PFMuonAlgo::postClean

(

reco::PFCandidateCollection *

cands

)

Definition at line 633 of file PFMuonAlgo.cc.

References HLT_2022v12_cff::cands, cleanMismeasured(), cleanPunchThroughAndFakes(), cosmicRejDistance_, ecal_dqm_sourceclient-live_cfg::cosmics, estimateEventQuantities(), eventFactorCosmics_, mps_fire::i, maskedIndices_, METX_, METY_, reco::PFCandidate::mu, reco::PFCandidate::muonRef(), PDWG_BPHSkim_cff::muons, pfCleanedTrackerAndGlobalMuonCandidates_, pfCosmicsMuonCleanedCandidates_, pfFakeMuonCleanedCandidates_, pfPunchThroughHadronCleanedCandidates_, pfPunchThroughMuonCleanedCandidates_, postCleaning_, reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), removeDeadCandidates(), jetUpdater_cfi::sort, sumet_, sumetPU_, and vertices_.

                                                           {
  using namespace std;
  using namespace reco;
  if (!postCleaning_)
    return;

  //Initialize vectors

  if (pfCosmicsMuonCleanedCandidates_.get())
    pfCosmicsMuonCleanedCandidates_->clear();
  else
    pfCosmicsMuonCleanedCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  if (pfCleanedTrackerAndGlobalMuonCandidates_.get())
    pfCleanedTrackerAndGlobalMuonCandidates_->clear();
  else
    pfCleanedTrackerAndGlobalMuonCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  if (pfFakeMuonCleanedCandidates_.get())
    pfFakeMuonCleanedCandidates_->clear();
  else
    pfFakeMuonCleanedCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  if (pfPunchThroughMuonCleanedCandidates_.get())
    pfPunchThroughMuonCleanedCandidates_->clear();
  else
    pfPunchThroughMuonCleanedCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  if (pfPunchThroughHadronCleanedCandidates_.get())
    pfPunchThroughHadronCleanedCandidates_->clear();
  else
    pfPunchThroughHadronCleanedCandidates_ = std::make_unique<reco::PFCandidateCollection>();

  pfPunchThroughHadronCleanedCandidates_->clear();

  maskedIndices_.clear();

  //Estimate MET and SumET

  estimateEventQuantities(cands);

  std::vector<int> muons;
  std::vector<int> cosmics;
  //get the muons
  for (unsigned int i = 0; i < cands->size(); ++i) {
    const reco::PFCandidate& cand = (*cands)[i];
    if (cand.particleId() == reco::PFCandidate::mu)
      muons.push_back(i);
  }
  //Then sort the muon indicess by decsending pt

  IndexPtComparator comparator(cands);
  std::sort(muons.begin(), muons.end(), comparator);

  //first kill cosmics
  double METXCosmics = 0;
  double METYCosmics = 0;
  double SUMETCosmics = 0.0;

  for (unsigned int i = 0; i < muons.size(); ++i) {
    const PFCandidate& pfc = cands->at(muons[i]);
    double origin = 0.0;
    if (!vertices_->empty() && vertices_->at(0).isValid() && !vertices_->at(0).isFake())
      origin = pfc.muonRef()->muonBestTrack()->dxy(vertices_->at(0).position());

    if (origin > cosmicRejDistance_) {
      cosmics.push_back(muons[i]);
      METXCosmics += pfc.px();
      METYCosmics += pfc.py();
      SUMETCosmics += pfc.pt();
    }
  }
  double MET2Cosmics = METXCosmics * METXCosmics + METYCosmics * METYCosmics;

  if (SUMETCosmics > (sumet_ - sumetPU_) / eventFactorCosmics_ && MET2Cosmics < METX_ * METX_ + METY_ * METY_)
    for (unsigned int i = 0; i < cosmics.size(); ++i) {
      maskedIndices_.push_back(cosmics[i]);
      pfCosmicsMuonCleanedCandidates_->push_back(cands->at(cosmics[i]));
    }

  //Loop on the muons candidates and clean
  for (unsigned int i = 0; i < muons.size(); ++i) {
    if (cleanMismeasured(cands->at(muons[i]), muons[i]))
      continue;
    cleanPunchThroughAndFakes(cands->at(muons[i]), cands, muons[i]);
  }

  //OK Now do the hard job ->remove the candidates that were cleaned
  removeDeadCandidates(cands, maskedIndices_);
}

printMuonProperties()

void PFMuonAlgo::printMuonProperties ( const reco::MuonRef &  muonRef )

static

Definition at line 355 of file PFMuonAlgo.cc.

References gather_cfg::cout, dumpMFGeometry_cfg::delta, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), SiStripPI::max, SiStripPI::min, particleFlowDisplacedVertex_cfi::ratio, mathSSE::sqrt(), egammaDQMOffline_cff::standAlone, muon::TM2DCompatibilityLoose, muon::TM2DCompatibilityTight, muon::TMLastStationAngLoose, muon::TMLastStationAngTight, muon::TMLastStationLoose, muon::TMLastStationOptimizedBarrelLowPtLoose, muon::TMLastStationOptimizedBarrelLowPtTight, muon::TMLastStationOptimizedLowPtLoose, muon::TMLastStationOptimizedLowPtTight, muon::TMLastStationTight, muon::TMOneStationLoose, muon::TMOneStationTight, HLT_2022v12_cff::track, and PbPb_ZMuSkimMuonDPG_cff::tracker.

                                                               {
  if (!muonRef.isNonnull())
    return;

  bool isGL = muonRef->isGlobalMuon();
  bool isTR = muonRef->isTrackerMuon();
  bool isST = muonRef->isStandAloneMuon();

  std::cout << " GL: " << isGL << " TR: " << isTR << " ST: " << isST << std::endl;
  std::cout << " nMatches " << muonRef->numberOfMatches() << std::endl;

  if (muonRef->isGlobalMuon()) {
    reco::TrackRef combinedMu = muonRef->combinedMuon();
    std::cout << " GL,  pt: " << combinedMu->pt() << " +/- " << combinedMu->ptError() / combinedMu->pt()
              << " chi**2 GBL : " << combinedMu->normalizedChi2() << std::endl;
    std::cout << " Total Muon Hits : " << combinedMu->hitPattern().numberOfValidMuonHits() << "/"
              << combinedMu->hitPattern().numberOfLostMuonHits()
              << " DT Hits : " << combinedMu->hitPattern().numberOfValidMuonDTHits() << "/"
              << combinedMu->hitPattern().numberOfLostMuonDTHits()
              << " CSC Hits : " << combinedMu->hitPattern().numberOfValidMuonCSCHits() << "/"
              << combinedMu->hitPattern().numberOfLostMuonCSCHits()
              << " RPC Hits : " << combinedMu->hitPattern().numberOfValidMuonRPCHits() << "/"
              << combinedMu->hitPattern().numberOfLostMuonRPCHits() << std::endl;

    std::cout << "  # of Valid Tracker Hits " << combinedMu->hitPattern().numberOfValidTrackerHits() << std::endl;
    std::cout << "  # of Valid Pixel Hits " << combinedMu->hitPattern().numberOfValidPixelHits() << std::endl;
  }
  if (muonRef->isStandAloneMuon()) {
    reco::TrackRef standAloneMu = muonRef->standAloneMuon();
    std::cout << " ST,  pt: " << standAloneMu->pt() << " +/- " << standAloneMu->ptError() / standAloneMu->pt()
              << " eta : " << standAloneMu->eta()
              << " DT Hits : " << standAloneMu->hitPattern().numberOfValidMuonDTHits() << "/"
              << standAloneMu->hitPattern().numberOfLostMuonDTHits()
              << " CSC Hits : " << standAloneMu->hitPattern().numberOfValidMuonCSCHits() << "/"
              << standAloneMu->hitPattern().numberOfLostMuonCSCHits()
              << " RPC Hits : " << standAloneMu->hitPattern().numberOfValidMuonRPCHits() << "/"
              << standAloneMu->hitPattern().numberOfLostMuonRPCHits()
              << " chi**2 STA : " << standAloneMu->normalizedChi2() << std::endl;
  }

  if (muonRef->isTrackerMuon()) {
    reco::TrackRef trackerMu = muonRef->track();
    const reco::Track& track = *trackerMu;
    std::cout << " TR,  pt: " << trackerMu->pt() << " +/- " << trackerMu->ptError() / trackerMu->pt()
              << " chi**2 TR : " << trackerMu->normalizedChi2() << std::endl;
    std::cout << " nTrackerHits " << track.hitPattern().numberOfValidTrackerHits() << std::endl;
    std::cout << "TMLastStationAngLoose               " << muon::isGoodMuon(*muonRef, muon::TMLastStationAngLoose)
              << std::endl
              << "TMLastStationAngTight               " << muon::isGoodMuon(*muonRef, muon::TMLastStationAngTight)
              << std::endl
              << "TMLastStationLoose               " << muon::isGoodMuon(*muonRef, muon::TMLastStationLoose)
              << std::endl
              << "TMLastStationTight               " << muon::isGoodMuon(*muonRef, muon::TMLastStationTight)
              << std::endl
              << "TMOneStationLoose                " << muon::isGoodMuon(*muonRef, muon::TMOneStationLoose) << std::endl
              << "TMOneStationTight                " << muon::isGoodMuon(*muonRef, muon::TMOneStationTight) << std::endl
              << "TMLastStationOptimizedLowPtLoose "
              << muon::isGoodMuon(*muonRef, muon::TMLastStationOptimizedLowPtLoose) << std::endl
              << "TMLastStationOptimizedLowPtTight "
              << muon::isGoodMuon(*muonRef, muon::TMLastStationOptimizedLowPtTight) << std::endl
              << "TMLastStationOptimizedBarrelLowPtLoose "
              << muon::isGoodMuon(*muonRef, muon::TMLastStationOptimizedBarrelLowPtLoose) << std::endl
              << "TMLastStationOptimizedBarrelLowPtTight "
              << muon::isGoodMuon(*muonRef, muon::TMLastStationOptimizedBarrelLowPtTight) << std::endl
              << std::endl;
  }

  std::cout << "TM2DCompatibilityLoose           " << muon::isGoodMuon(*muonRef, muon::TM2DCompatibilityLoose)
            << std::endl
            << "TM2DCompatibilityTight           " << muon::isGoodMuon(*muonRef, muon::TM2DCompatibilityTight)
            << std::endl;

  if (muonRef->isGlobalMuon() && muonRef->isTrackerMuon() && muonRef->isStandAloneMuon()) {
    reco::TrackRef combinedMu = muonRef->combinedMuon();
    reco::TrackRef trackerMu = muonRef->track();
    reco::TrackRef standAloneMu = muonRef->standAloneMuon();

    double sigmaCombined = combinedMu->ptError() / (combinedMu->pt() * combinedMu->pt());
    double sigmaTracker = trackerMu->ptError() / (trackerMu->pt() * trackerMu->pt());
    double sigmaStandAlone = standAloneMu->ptError() / (standAloneMu->pt() * standAloneMu->pt());

    bool combined = combinedMu->ptError() / combinedMu->pt() < 0.20;
    bool tracker = trackerMu->ptError() / trackerMu->pt() < 0.20;
    bool standAlone = standAloneMu->ptError() / standAloneMu->pt() < 0.20;

    double delta1 = combined && tracker ? fabs(1. / combinedMu->pt() - 1. / trackerMu->pt()) /
                                              sqrt(sigmaCombined * sigmaCombined + sigmaTracker * sigmaTracker)
                                        : 100.;
    double delta2 = combined && standAlone ? fabs(1. / combinedMu->pt() - 1. / standAloneMu->pt()) /
                                                 sqrt(sigmaCombined * sigmaCombined + sigmaStandAlone * sigmaStandAlone)
                                           : 100.;
    double delta3 = standAlone && tracker ? fabs(1. / standAloneMu->pt() - 1. / trackerMu->pt()) /
                                                sqrt(sigmaStandAlone * sigmaStandAlone + sigmaTracker * sigmaTracker)
                                          : 100.;

    double delta =
        standAloneMu->hitPattern().numberOfValidMuonDTHits() + standAloneMu->hitPattern().numberOfValidMuonCSCHits() > 0
            ? std::min(delta3, std::min(delta1, delta2))
            : std::max(delta3, std::max(delta1, delta2));

    std::cout << "delta = " << delta << " delta1 " << delta1 << " delta2 " << delta2 << " delta3 " << delta3
              << std::endl;

    double ratio = combinedMu->ptError() / combinedMu->pt() / (trackerMu->ptError() / trackerMu->pt());
    //if ( ratio > 2. && delta < 3. ) std::cout << "ALARM ! " << ratio << ", " << delta << std::endl;
    std::cout << " ratio " << ratio << " combined mu pt " << combinedMu->pt() << std::endl;
    //bool quality3 =  ( combinedMu->pt() < 50. || ratio < 2. ) && delta <  3.;
  }

  double sumPtR03 = muonRef->isolationR03().sumPt;
  double emEtR03 = muonRef->isolationR03().emEt;
  double hadEtR03 = muonRef->isolationR03().hadEt;
  double relIsoR03 = (sumPtR03 + emEtR03 + hadEtR03) / muonRef->pt();
  double sumPtR05 = muonRef->isolationR05().sumPt;
  double emEtR05 = muonRef->isolationR05().emEt;
  double hadEtR05 = muonRef->isolationR05().hadEt;
  double relIsoR05 = (sumPtR05 + emEtR05 + hadEtR05) / muonRef->pt();
  std::cout << " 0.3 Radion Rel Iso: " << relIsoR03 << " sumPt " << sumPtR03 << " emEt " << emEtR03 << " hadEt "
            << hadEtR03 << std::endl;
  std::cout << " 0.5 Radion Rel Iso: " << relIsoR05 << " sumPt " << sumPtR05 << " emEt " << emEtR05 << " hadEt "
            << hadEtR05 << std::endl;
  return;
}

reconstructMuon()

bool PFMuonAlgo::reconstructMuon	(	reco::PFCandidate &	candidate,
		const reco::MuonRef &	muon,
		bool	allowLoose = false
	)

Definition at line 537 of file PFMuonAlgo.cc.

References changeTrack(), errorCompScale_, getTrackWithSmallestError(), reco::Muon::InnerTrack, isLooseMuon(), jetsAK4_CHS_cff::isMu, isMuon(), maxDPtOPt_, muonTracks(), reco::PFCandidate::setMuonRef(), muon::tevOptimized(), trackQuality_, and PbPb_ZMuSkimMuonDPG_cff::trackType.

                                                                                                     {
  using namespace std;
  using namespace reco;

  if (!muon.isNonnull())
    return false;

  bool isMu = false;

  if (allowLoose)
    isMu = isMuon(muon) || isLooseMuon(muon);
  else
    isMu = isMuon(muon);

  if (!isMu)
    return false;

  //get the valid tracks(without standalone except we allow loose muons)
  //MIKE: Here we need to be careful. If we have a muon inside a dense
  //jet environment often the track is badly measured. In this case
  //we should not apply Dpt/Pt<1

  std::vector<reco::Muon::MuonTrackTypePair> validTracks = muonTracks(muon, maxDPtOPt_);
  if (!allowLoose)
    validTracks = muonTracks(muon, maxDPtOPt_);
  else
    validTracks = muonTracks(muon);

  if (validTracks.empty())
    return false;

  //check what is the track used.Rerun TuneP
  reco::Muon::MuonTrackTypePair bestTrackPair = muon::tevOptimized(*muon);

  TrackRef bestTrack = bestTrackPair.first;
  MuonTrackType trackType = bestTrackPair.second;

  MuonTrackTypePair trackPairWithSmallestError = getTrackWithSmallestError(validTracks);
  TrackRef trackWithSmallestError = trackPairWithSmallestError.first;

  if (trackType == reco::Muon::InnerTrack &&
      (!bestTrack->quality(trackQuality_) ||
       bestTrack->ptError() / bestTrack->pt() >
           errorCompScale_ * trackWithSmallestError->ptError() / trackWithSmallestError->pt())) {
    bestTrack = trackWithSmallestError;
    trackType = trackPairWithSmallestError.second;
  } else if (trackType != reco::Muon::InnerTrack &&
             bestTrack->ptError() / bestTrack->pt() >
                 errorCompScale_ * trackWithSmallestError->ptError() / trackWithSmallestError->pt()) {
    bestTrack = trackWithSmallestError;
    trackType = trackPairWithSmallestError.second;
  }

  changeTrack(candidate, std::make_pair(bestTrack, trackType));
  candidate.setMuonRef(muon);

  return true;
}

removeDeadCandidates()

void PFMuonAlgo::removeDeadCandidates ( reco::PFCandidateCollection *  obj, const std::vector< unsigned int > &  indices  )

private

Definition at line 1033 of file PFMuonAlgo.cc.

References mps_fire::i, dqmdumpme::indices, N, and getGTfromDQMFile::obj.

Referenced by postClean().

                                                                                                            {
  size_t N = indices.size();
  size_t collSize = obj->size();

  for (size_t i = 0; i < N; ++i)
    obj->at(indices.at(i)) = obj->at(collSize - i - 1);

  obj->resize(collSize - indices.size());
}

setInputsForCleaning()

void PFMuonAlgo::setInputsForCleaning

(

reco::VertexCollection const &

vertices

)

Definition at line 941 of file PFMuonAlgo.cc.

References AlignmentTracksFromVertexSelector_cfi::vertices, and vertices_.

{ vertices_ = &vertices; }

setVetoes()

void PFMuonAlgo::setVetoes ( const reco::PFCandidateCollection &  vetoes )

inline

Definition at line 54 of file PFMuonAlgo.h.

References vetoes_.

Referenced by PFProducer::produce().

{ vetoes_ = &vetoes; }

tracksPointingAtMET()

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::tracksPointingAtMET ( const std::vector< MuonTrackTypePair > &  tracks )

private

Definition at line 924 of file PFMuonAlgo.cc.

References mps_fire::i, metFactorCleaning_, METX_, METY_, funct::pow(), and tracks.

Referenced by addMissingMuons().

                                                        {
  std::vector<reco::Muon::MuonTrackTypePair> outputTracks;

  double newMET2 = 0.0;

  for (unsigned int i = 0; i < tracks.size(); ++i) {
    //calculate new SUM ET and MET2
    newMET2 = pow(METX_ + tracks.at(i).first->px(), 2) + pow(METY_ + tracks.at(i).first->py(), 2);

    if (newMET2 < (METX_ * METX_ + METY_ * METY_) / metFactorCleaning_)
      outputTracks.push_back(tracks.at(i));
  }

  return outputTracks;
}

tracksWithBetterMET()

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::tracksWithBetterMET ( const std::vector< MuonTrackTypePair > &  tracks, const reco::PFCandidate &  pfc  )

private

Definition at line 899 of file PFMuonAlgo.cc.

References eventFractionCleaning_, mps_fire::i, metFactorCleaning_, metSigForCleaning_, METX_, METY_, funct::pow(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), sumet_, sumetPU_, and tracks.

Referenced by cleanMismeasured().

                                                                                    {
  std::vector<reco::Muon::MuonTrackTypePair> outputTracks;

  double METNOX = METX_ - pfc.px();
  double METNOY = METY_ - pfc.py();
  double SUMETNO = sumet_ - pfc.pt();
  double MET2 = METX_ * METX_ + METY_ * METY_;
  double newMET2 = 0.0;
  double newSUMET = 0.0;
  double METSIG = sqrt(MET2) / sqrt(sumet_ - sumetPU_);

  if (METSIG > metSigForCleaning_)
    for (unsigned int i = 0; i < tracks.size(); ++i) {
      //calculate new SUM ET and MET2
      newSUMET = SUMETNO + tracks.at(i).first->pt() - sumetPU_;
      newMET2 = pow(METNOX + tracks.at(i).first->px(), 2) + pow(METNOY + tracks.at(i).first->py(), 2);

      if (newSUMET / (sumet_ - sumetPU_) > eventFractionCleaning_ && newMET2 < MET2 / metFactorCleaning_)
        outputTracks.push_back(tracks.at(i));
    }

  return outputTracks;
}

transferAddedMuonCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferAddedMuonCandidates ( )

inline

Definition at line 76 of file PFMuonAlgo.h.

References eostools::move(), and pfAddedMuonCandidates_.

                                                                         {
    return std::move(pfAddedMuonCandidates_);
  }

transferCleanedCosmicCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferCleanedCosmicCandidates ( )

inline

Definition at line 56 of file PFMuonAlgo.h.

References eostools::move(), and pfCosmicsMuonCleanedCandidates_.

                                                                             {
    return std::move(pfCosmicsMuonCleanedCandidates_);
  }

transferCleanedFakeCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferCleanedFakeCandidates ( )

inline

Definition at line 64 of file PFMuonAlgo.h.

References eostools::move(), and pfFakeMuonCleanedCandidates_.

                                                                           {
    return std::move(pfFakeMuonCleanedCandidates_);
  }

transferCleanedTrackerAndGlobalCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferCleanedTrackerAndGlobalCandidates ( )

inline

Definition at line 60 of file PFMuonAlgo.h.

References eostools::move(), and pfCleanedTrackerAndGlobalMuonCandidates_.

                                                                                       {
    return std::move(pfCleanedTrackerAndGlobalMuonCandidates_);
  }

transferPunchThroughCleanedHadronCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferPunchThroughCleanedHadronCandidates ( )

inline

Definition at line 72 of file PFMuonAlgo.h.

References eostools::move(), and pfPunchThroughHadronCleanedCandidates_.

                                                                                         {
    return std::move(pfPunchThroughHadronCleanedCandidates_);
  }

transferPunchThroughCleanedMuonCandidates()

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::transferPunchThroughCleanedMuonCandidates ( )

inline

Definition at line 68 of file PFMuonAlgo.h.

References eostools::move(), and pfPunchThroughMuonCleanedCandidates_.

                                                                                       {
    return std::move(pfPunchThroughMuonCleanedCandidates_);
  }

Member Data Documentation

cosmicRejDistance_

const double PFMuonAlgo::cosmicRejDistance_

private

Definition at line 148 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), and postClean().

errorCompScale_

const double PFMuonAlgo::errorCompScale_

private

Definition at line 130 of file PFMuonAlgo.h.

Referenced by reconstructMuon().

eventFactorCosmics_

const double PFMuonAlgo::eventFactorCosmics_

private

Definition at line 135 of file PFMuonAlgo.h.

Referenced by postClean().

eventFractionCleaning_

const double PFMuonAlgo::eventFractionCleaning_

private

Definition at line 133 of file PFMuonAlgo.h.

Referenced by tracksWithBetterMET().

eventFractionRejection_

const double PFMuonAlgo::eventFractionRejection_

private

Definition at line 139 of file PFMuonAlgo.h.

Referenced by cleanMismeasured().

maskedIndices_

std::vector<unsigned int> PFMuonAlgo::maskedIndices_

private

Definition at line 118 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), cleanPunchThroughAndFakes(), and postClean().

maxDPtOPt_

const double PFMuonAlgo::maxDPtOPt_

private

Definition at line 127 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), cleanMismeasured(), getMinMaxMET2(), and reconstructMuon().

metFactorCleaning_

const double PFMuonAlgo::metFactorCleaning_

private

Definition at line 138 of file PFMuonAlgo.h.

Referenced by tracksPointingAtMET(), and tracksWithBetterMET().

metFactorFake_

const double PFMuonAlgo::metFactorFake_

private

Definition at line 143 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

metFactorHighEta_

const double PFMuonAlgo::metFactorHighEta_

private

Definition at line 141 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

metFactorRejection_

const double PFMuonAlgo::metFactorRejection_

private

Definition at line 140 of file PFMuonAlgo.h.

Referenced by cleanMismeasured().

metSigForCleaning_

const double PFMuonAlgo::metSigForCleaning_

private

Definition at line 136 of file PFMuonAlgo.h.

Referenced by tracksWithBetterMET().

metSigForRejection_

const double PFMuonAlgo::metSigForRejection_

private

Definition at line 137 of file PFMuonAlgo.h.

Referenced by cleanMismeasured().

METX_

double PFMuonAlgo::METX_

private

Definition at line 152 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), cleanMismeasured(), cleanPunchThroughAndFakes(), estimateEventQuantities(), getMinMaxMET2(), postClean(), tracksPointingAtMET(), and tracksWithBetterMET().

METY_

double PFMuonAlgo::METY_

private

Definition at line 153 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), cleanMismeasured(), cleanPunchThroughAndFakes(), estimateEventQuantities(), getMinMaxMET2(), postClean(), tracksPointingAtMET(), and tracksWithBetterMET().

minPostCleaningPt_

const double PFMuonAlgo::minPostCleaningPt_

private

Definition at line 134 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), and cleanPunchThroughAndFakes().

minPunchThroughEnergy_

const double PFMuonAlgo::minPunchThroughEnergy_

private

Definition at line 145 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

minPunchThroughMomentum_

const double PFMuonAlgo::minPunchThroughMomentum_

private

Definition at line 144 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

pfAddedMuonCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfAddedMuonCandidates_

private

the collection of added muon candidates

Definition at line 116 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), and transferAddedMuonCandidates().

pfCleanedTrackerAndGlobalMuonCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfCleanedTrackerAndGlobalMuonCandidates_

private

the collection of tracker/global cleaned muon candidates

Definition at line 108 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), postClean(), and transferCleanedTrackerAndGlobalCandidates().

pfCosmicsMuonCleanedCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfCosmicsMuonCleanedCandidates_

private

the collection of cosmics cleaned muon candidates

Definition at line 106 of file PFMuonAlgo.h.

Referenced by postClean(), and transferCleanedCosmicCandidates().

pfFakeMuonCleanedCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfFakeMuonCleanedCandidates_

private

the collection of fake cleaned muon candidates

Definition at line 110 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), cleanPunchThroughAndFakes(), postClean(), and transferCleanedFakeCandidates().

pfPunchThroughHadronCleanedCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfPunchThroughHadronCleanedCandidates_

private

the collection of punch-through cleaned neutral hadron candidates

Definition at line 114 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes(), postClean(), and transferPunchThroughCleanedHadronCandidates().

pfPunchThroughMuonCleanedCandidates_

std::unique_ptr<reco::PFCandidateCollection> PFMuonAlgo::pfPunchThroughMuonCleanedCandidates_

private

the collection of punch-through cleaned muon candidates

Definition at line 112 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes(), postClean(), and transferPunchThroughCleanedMuonCandidates().

postCleaning_

const bool PFMuonAlgo::postCleaning_

private

Definition at line 132 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), and postClean().

ptFactorHighEta_

const double PFMuonAlgo::ptFactorHighEta_

private

Definition at line 142 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

punchThroughFactor_

const double PFMuonAlgo::punchThroughFactor_

private

Definition at line 146 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

punchThroughMETFactor_

const double PFMuonAlgo::punchThroughMETFactor_

private

Definition at line 147 of file PFMuonAlgo.h.

Referenced by cleanPunchThroughAndFakes().

sumet_

double PFMuonAlgo::sumet_

private

Definition at line 151 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), cleanPunchThroughAndFakes(), estimateEventQuantities(), postClean(), and tracksWithBetterMET().

sumetPU_

double PFMuonAlgo::sumetPU_

private

Definition at line 150 of file PFMuonAlgo.h.

Referenced by cleanMismeasured(), cleanPunchThroughAndFakes(), estimateEventQuantities(), postClean(), and tracksWithBetterMET().

trackQuality_

const reco::TrackBase::TrackQuality PFMuonAlgo::trackQuality_

private

Definition at line 129 of file PFMuonAlgo.h.

Referenced by reconstructMuon().

vertices_

const reco::VertexCollection* PFMuonAlgo::vertices_

private

Definition at line 123 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), postClean(), and setInputsForCleaning().

vetoes_

const reco::PFCandidateCollection* PFMuonAlgo::vetoes_ = nullptr

private

Definition at line 120 of file PFMuonAlgo.h.

Referenced by addMissingMuons(), and setVetoes().