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.

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

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

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

matchTracks()

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

Definition at line 11 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().

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