CMS 3D CMS Logo

Functions
muonid Namespace Reference

Functions

reco::TrackRef findOppositeTrack (const edm::Handle< reco::TrackCollection > &collection, const reco::Track &muon, double angleMatch=0.01, double momentumMatch=0.05)
 
std::pair< double, double > matchTracks (const reco::Track &ref, const reco::Track &probe)
 

Function Documentation

◆ findOppositeTrack()

reco::TrackRef muonid::findOppositeTrack ( const edm::Handle< reco::TrackCollection > &  collection,
const reco::Track muon,
double  angleMatch = 0.01,
double  momentumMatch = 0.05 
)

Definition at line 28 of file MuonCosmicsId.cc.

31  {
32  for (unsigned int i = 0; i < tracks->size(); ++i) {
33  // When both tracks are reconstructed as outside going, sign is -1
34  // otherwise it's +1. There is also a crazy case of both are outside
35  // going, then sign is -1 as well.
36  const std::pair<double, double>& match = matchTracks(muonTrack, tracks->at(i));
37  if (match.first < angleMatch && match.second < momentumMatch)
38  return reco::TrackRef(tracks, i);
39  }
40  return reco::TrackRef();
41 }

References mps_fire::i, match(), matchTracks(), and PDWG_EXOHSCP_cff::tracks.

Referenced by MuonCosmicCompatibilityFiller::backToBack2LegCosmic(), CosmicsMuonIdProducer::produce(), and TrajectorySeedFromMuonProducer::produce().

◆ matchTracks()

std::pair< double, double > muonid::matchTracks ( const reco::Track ref,
const reco::Track probe 
)

Definition at line 11 of file MuonCosmicsId.cc.

11  {
12  std::pair<double, double> result(0, 0);
13  // When both tracks are reconstructed as outside going, sign is -1
14  // otherwise it's +1. There is also a crazy case of both are outside
15  // going, then sign is -1 as well.
16  int match_sign = directionAlongMomentum(ref) == directionAlongMomentum(probe) ? -1 : +1;
17  double sprod = ref.px() * probe.px() + ref.py() * probe.py() + ref.pz() * probe.pz();
18  double argCos = match_sign * (sprod / ref.p() / probe.p());
19  if (argCos < -1.0)
20  argCos = -1.0;
21  if (argCos > 1.0)
22  argCos = 1.0;
23  result.first = acos(argCos);
24  result.second = fabs(probe.pt() - ref.pt()) / sqrt(ref.pt() * probe.pt()); //SK: take a geom-mean pt
25  return result;
26 }

References directionAlongMomentum(), reco::TrackBase::p(), reco::TrackBase::pt(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), mps_fire::result, and mathSSE::sqrt().

Referenced by findOppositeTrack().

PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
mps_fire.i
i
Definition: mps_fire.py:355
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:605
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:611
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:608
reco::TrackBase::py
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:614
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
reco::TrackRef
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:20
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
muonid::matchTracks
std::pair< double, double > matchTracks(const reco::Track &ref, const reco::Track &probe)
Definition: MuonCosmicsId.cc:11
directionAlongMomentum
bool directionAlongMomentum(const reco::Track &track)
Definition: MuonCosmicsId.cc:4
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:617
mps_fire.result
result
Definition: mps_fire.py:303