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

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 29 of file MuonCosmicsId.cc.

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

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

33 {
34  for (unsigned int i=0; i<tracks->size(); ++i){
35  // When both tracks are reconstructed as outside going, sign is -1
36  // otherwise it's +1. There is also a crazy case of both are outside
37  // going, then sign is -1 as well.
38  const std::pair<double,double>& match = matchTracks(muonTrack,tracks->at(i));
39  if ( match.first < angleMatch && match.second < momentumMatch )
40  return reco::TrackRef(tracks,i);
41  }
42  return reco::TrackRef();
43 
44 }
std::pair< double, double > matchTracks(const reco::Track &ref, const reco::Track &probe)
edm::Ref< TrackCollection > TrackRef
persistent reference to a Track
Definition: TrackFwd.h:21
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
std::pair< double, double > muonid::matchTracks ( const reco::Track ref,
const reco::Track probe 
)

Definition at line 12 of file MuonCosmicsId.cc.

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(), and JetPlusTrackCorrector::getSumPtForBeta().

13 {
14  std::pair<double,double> result(0,0);
15  // When both tracks are reconstructed as outside going, sign is -1
16  // otherwise it's +1. There is also a crazy case of both are outside
17  // going, then sign is -1 as well.
18  int match_sign = directionAlongMomentum(ref)==directionAlongMomentum(probe) ? -1 : +1;
19  double sprod = ref.px()*probe.px() + ref.py()*probe.py() + ref.pz()*probe.pz();
20  double argCos = match_sign*(sprod/ref.p()/probe.p());
21  if (argCos < -1.0) argCos = -1.0;
22  if (argCos > 1.0) 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 }
double p() const
momentum vector magnitude
Definition: TrackBase.h:615
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:627
bool directionAlongMomentum(const reco::Track &track)
Definition: MuonCosmicsId.cc:4
T sqrt(T t)
Definition: SSEVec.h:18
double pt() const
track transverse momentum
Definition: TrackBase.h:621
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:639
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:633