SoftMuonMvaRun3Estimator.cc File Reference

#include "PhysicsTools/PatAlgos/interface/SoftMuonMvaRun3Estimator.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "PhysicsTools/XGBoost/interface/XGBooster.h"

Go to the source code of this file.

Typedefs
typedef std::pair< const reco::MuonChamberMatch *, const reco::MuonSegmentMatch * >	MatchPair

Functions
float	dX (const MatchPair &match)
float	dY (const MatchPair &match)
void	fillMatchInfo (pat::XGBooster &booster, const pat::Muon &muon)
void	fillMatchInfoForStation (std::string prefix, pat::XGBooster &booster, const MatchPair &match)
const MatchPair &	getBetterMatch (const MatchPair &match1, const MatchPair &match2)
float	pullDxDz (const MatchPair &match)
float	pullDyDz (const MatchPair &match)
float	pullX (const MatchPair &match)
float	pullY (const MatchPair &match)

Typedef Documentation

MatchPair

typedef std::pair<const reco::MuonChamberMatch*, const reco::MuonSegmentMatch*> MatchPair

Definition at line 5 of file SoftMuonMvaRun3Estimator.cc.

Function Documentation

dX()

float dX

(

const MatchPair &

match

)

Definition at line 23 of file SoftMuonMvaRun3Estimator.cc.

References match().

Referenced by RPCPointVsRecHit::analyze(), RPCRecHitValid::analyze(), Quantile::calculateQ(), FWFramedTextTableCellRenderer::draw(), FWTextTableCellRenderer::draw(), FWRPZViewGeometry::estimateProjectionSizeDT(), MuonShowerDigiFiller::fill(), fillMatchInfoForStation(), getVertexD0(), muon::isGoodMuon(), MuonDetCleaner< T1, T2 >::produce(), pullX(), and LinkByRecHit::testHFEMAndHFHADByRecHit().

                                 {
  if (match.first and match.second->hasPhi())
    return (match.first->x - match.second->x);
  else
    return 9999.;
}

dY()

float dY

(

const MatchPair &

match

)

Definition at line 45 of file SoftMuonMvaRun3Estimator.cc.

References match().

Referenced by Quantile::calculateQ(), FWCollectionSummaryModelCellRenderer::draw(), FWColumnLabelCellRenderer::draw(), FWRPZViewGeometry::estimateProjectionSizeDT(), fillMatchInfoForStation(), getVertexD0(), muon::isGoodMuon(), pullY(), and LinkByRecHit::testHFEMAndHFHADByRecHit().

                                 {
  if (match.first and match.second->hasZed())
    return (match.first->y - match.second->y);
  else
    return 9999.;
}

fillMatchInfo()

void fillMatchInfo	(	pat::XGBooster &	booster,
		const pat::Muon &	muon
	)

Definition at line 76 of file SoftMuonMvaRun3Estimator.cc.

References reco::MuonSegmentMatch::BelongsToTrackByDR, reco::MuonSegmentMatch::BestInStationByDR, fillMatchInfoForStation(), dqmdumpme::first, getBetterMatch(), mps_fire::i, oniaPATMuonsWithTrigger_cff::matches, and relativeConstraints::station.

Referenced by pat::computeSoftMvaRun3().

                                                               {
  // Initiate containter for results
  const int n_stations = 2;
  std::vector<MatchPair> matches;
  for (unsigned int i = 0; i < n_stations; ++i)
    matches.push_back(std::pair(nullptr, nullptr));

  // Find best matches
  for (auto& chamberMatch : muon.matches()) {
    unsigned int station = chamberMatch.station() - 1;
    if (station >= n_stations)
      continue;

    // Find best segment match.
    // We could consider all segments, but we will restrict to segments
    // that match to this candidate better than to other muon candidates
    for (auto& segmentMatch : chamberMatch.segmentMatches) {
      if (not segmentMatch.isMask(reco::MuonSegmentMatch::BestInStationByDR) ||
          not segmentMatch.isMask(reco::MuonSegmentMatch::BelongsToTrackByDR))
        continue;

      // Multiple segment matches are possible in different
      // chambers that are either overlapping or belong to
      // different detectors. We need to select one.
      auto match_pair = MatchPair(&chamberMatch, &segmentMatch);

      if (matches[station].first)
        matches[station] = getBetterMatch(matches[station], match_pair);
      else
        matches[station] = match_pair;
    }
  }

  // Fill matching information
  fillMatchInfoForStation("match1", booster, matches[0]);
  fillMatchInfoForStation("match2", booster, matches[1]);
}

fillMatchInfoForStation()

void fillMatchInfoForStation	(	std::string	prefix,
		pat::XGBooster &	booster,
		const MatchPair &	match
	)

Definition at line 67 of file SoftMuonMvaRun3Estimator.cc.

References dX(), dY(), match(), hcallasereventfilter2012_cfi::prefix, pullDxDz(), pullDyDz(), pullX(), pullY(), and pat::XGBooster::set().

Referenced by fillMatchInfo().

                                                                                              {
  booster.set(prefix + "_dX", dX(match));
  booster.set(prefix + "_pullX", pullX(match));
  booster.set(prefix + "_pullDxDz", pullDxDz(match));
  booster.set(prefix + "_dY", dY(match));
  booster.set(prefix + "_pullY", pullY(match));
  booster.set(prefix + "_pullDyDz", pullDyDz(match));
}

getBetterMatch()

const MatchPair& getBetterMatch	(	const MatchPair &	match1,
		const MatchPair &	match2
	)

Definition at line 7 of file SoftMuonMvaRun3Estimator.cc.

References funct::abs(), and MuonSubdetId::DT.

Referenced by fillMatchInfo().

                                                                                  {
  // Prefer DT over CSC simply because it's closer to IP
  // and will have less multiple scattering (at least for
  // RB1 vs ME1/3 case). RB1 & ME1/2 overlap is tiny
  if (match2.first->detector() == MuonSubdetId::DT and match1.first->detector() != MuonSubdetId::DT)
    return match2;

  // For the rest compare local x match. We expect that
  // segments belong to the muon, so the difference in
  // local x is a reflection on how well we can measure it
  if (abs(match1.first->x - match1.second->x) > abs(match2.first->x - match2.second->x))
    return match2;

  return match1;
}

pullDxDz()

float pullDxDz

(

const MatchPair &

match

)

Definition at line 37 of file SoftMuonMvaRun3Estimator.cc.

References match(), funct::pow(), and mathSSE::sqrt().

Referenced by fillMatchInfoForStation().

                                       {
  if (match.first and match.second->hasPhi())
    return (match.first->dXdZ - match.second->dXdZ) /
           sqrt(pow(match.first->dXdZErr, 2) + pow(match.second->dXdZErr, 2));
  else
    return 9999.;
}

pullDyDz()

float pullDyDz

(

const MatchPair &

match

)

Definition at line 59 of file SoftMuonMvaRun3Estimator.cc.

References match(), funct::pow(), and mathSSE::sqrt().

Referenced by fillMatchInfoForStation().

                                       {
  if (match.first and match.second->hasZed())
    return (match.first->dYdZ - match.second->dYdZ) /
           sqrt(pow(match.first->dYdZErr, 2) + pow(match.second->dYdZErr, 2));
  else
    return 9999.;
}

pullX()

float pullX

(

const MatchPair &

match

)

Definition at line 30 of file SoftMuonMvaRun3Estimator.cc.

References dX(), match(), funct::pow(), and mathSSE::sqrt().

Referenced by ME0RecHitsValidation::analyze(), DMRChecker::analyze(), CkfDebugger::computePulls(), fillMatchInfoForStation(), muon::overlap(), and SimpleVertexTree::~SimpleVertexTree().

                                    {
  if (match.first and match.second->hasPhi())
    return dX(match) / sqrt(pow(match.first->xErr, 2) + pow(match.second->xErr, 2));
  else
    return 9999.;
}

pullY()

float pullY

(

const MatchPair &

match

)

Definition at line 52 of file SoftMuonMvaRun3Estimator.cc.

References dY(), match(), funct::pow(), and mathSSE::sqrt().

Referenced by ME0RecHitsValidation::analyze(), DMRChecker::analyze(), CkfDebugger::computePulls(), fillMatchInfoForStation(), muon::overlap(), and SimpleVertexTree::~SimpleVertexTree().

                                    {
  if (match.first and match.second->hasZed())
    return dY(match) / sqrt(pow(match.first->yErr, 2) + pow(match.second->yErr, 2));
  else
    return 9999.;
}