GlobalMuonTrackMatcher Class Reference

#include <GlobalMuonTrackMatcher.h>

Public Types
typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand

Public Member Functions
std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface >	convertToTSOSMuHit (const TrackCand &, const TrackCand &) const
std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface >	convertToTSOSTk (const TrackCand &, const TrackCand &) const
std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface >	convertToTSOSTkHit (const TrackCand &, const TrackCand &) const
	GlobalMuonTrackMatcher (const edm::ParameterSet &, const MuonServiceProxy *)
	constructor More...
double	match (const TrackCand &sta, const TrackCand &track, int matchOption=0, int surfaceOption=1) const
std::vector< TrackCand >	match (const TrackCand &sta, const std::vector< TrackCand > &tracks) const
	choose all tracks with a matching-chi2 less than Chi2Cut More...
double	match_Chi2 (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
double	match_D (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
double	match_d (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
double	match_dist (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
double	match_R_IP (const TrackCand &, const TrackCand &) const
double	match_Rmom (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
double	match_Rpos (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
std::vector< TrackCand >::const_iterator	matchOne (const TrackCand &sta, const std::vector< TrackCand > &tracks) const
	choose the one tracker track which best matches a muon track More...
bool	matchTight (const TrackCand &sta, const TrackCand &track) const
	check if two tracks are compatible (less than Chi2Cut, DeltaDCut, DeltaRCut) More...
bool	samePlane (const TrajectoryStateOnSurface &, const TrajectoryStateOnSurface &) const
virtual	~GlobalMuonTrackMatcher ()
	destructor More...

Private Attributes
double	theChi2_1
double	theChi2_2
double	theChi2_3
double	theDeltaD_1
double	theDeltaD_2
double	theDeltaD_3
double	theDeltaR_1
double	theDeltaR_2
double	theDeltaR_3
double	theEta_threshold
double	theLocChi2
double	theMinP
double	theMinPt
std::string	theOutPropagatorName
double	thePt_threshold1
double	thePt_threshold2
double	theQual_1
double	theQual_2
double	theQual_3
const MuonServiceProxy *	theService

Detailed Description

Match a standalone muon track with the most compatible tracker tracks in a TrackCollection. GlobalMuonTrackMatcher is used during global muon reconstruction to check the compatability of a tracker track with a standalone muon track. The compatability is determined on a chi2 comparison, of the local parameters of the two corresponding TrajectoryStates on the surface of the innermost muon measurement for momentum below a threshold and above this, with the position and direction parameters on the mentioned surface. If the comparison of local parameters fails to yield any matches, then it makes a comparison of the TSOS local direction.

Author

Edwin Antillon Purdue University

Chang Liu Purdue University

Adam Everett Purdue University

Norbert Neumeister Purdue University

Definition at line 39 of file GlobalMuonTrackMatcher.h.

Member Typedef Documentation

typedef std::pair<const Trajectory*, reco::TrackRef> GlobalMuonTrackMatcher::TrackCand

Definition at line 41 of file GlobalMuonTrackMatcher.h.

Constructor & Destructor Documentation

GlobalMuonTrackMatcher::GlobalMuonTrackMatcher	(	const edm::ParameterSet &	par,
		const MuonServiceProxy *	service
	)

constructor

Definition at line 49 of file GlobalMuonTrackMatcher.cc.

References edm::ParameterSet::getParameter(), theChi2_1, theChi2_2, theChi2_3, theDeltaD_1, theDeltaD_2, theDeltaD_3, theDeltaR_1, theDeltaR_2, theDeltaR_3, theEta_threshold, theLocChi2, theMinP, theMinPt, theOutPropagatorName, thePt_threshold1, thePt_threshold2, theQual_1, theQual_2, and theQual_3.

    : theService(service) {
  theMinP = par.getParameter<double>("MinP");
  theMinPt = par.getParameter<double>("MinPt");
  thePt_threshold1 = par.getParameter<double>("Pt_threshold1");
  thePt_threshold2 = par.getParameter<double>("Pt_threshold2");
  theEta_threshold = par.getParameter<double>("Eta_threshold");
  theChi2_1 = par.getParameter<double>("Chi2Cut_1");
  theChi2_2 = par.getParameter<double>("Chi2Cut_2");
  theChi2_3 = par.getParameter<double>("Chi2Cut_3");
  theLocChi2 = par.getParameter<double>("LocChi2Cut");
  theDeltaD_1 = par.getParameter<double>("DeltaDCut_1");
  theDeltaD_2 = par.getParameter<double>("DeltaDCut_2");
  theDeltaD_3 = par.getParameter<double>("DeltaDCut_3");
  theDeltaR_1 = par.getParameter<double>("DeltaRCut_1");
  theDeltaR_2 = par.getParameter<double>("DeltaRCut_2");
  theDeltaR_3 = par.getParameter<double>("DeltaRCut_3");
  theQual_1 = par.getParameter<double>("Quality_1");
  theQual_2 = par.getParameter<double>("Quality_2");
  theQual_3 = par.getParameter<double>("Quality_3");
  theOutPropagatorName = par.getParameter<string>("Propagator");
}

GlobalMuonTrackMatcher::~GlobalMuonTrackMatcher ( )

virtual

destructor

Definition at line 75 of file GlobalMuonTrackMatcher.cc.

{}

Member Function Documentation

std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface > GlobalMuonTrackMatcher::convertToTSOSMuHit	(	const TrackCand &	staCand,
		const TrackCand &	tkCand
	)		const

Definition at line 417 of file GlobalMuonTrackMatcher.cc.

References taus_updatedMVAIds_cff::category, relativeConstraints::empty, TrajectoryStateOnSurface::freeState(), TrajectoryStateOnSurface::globalPosition(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), LogTrace, mag(), trajectoryStateTransform::outerStateOnSurface(), theMinP, theMinPt, theOutPropagatorName, and theService.

Referenced by match(), matchOne(), and matchTight().

                                                             {
  const string category = "GlobalMuonTrackMatcher";
  TrajectoryStateOnSurface empty;
  TransientTrack muTT(*staCand.second, &*theService->magneticField(), theService->trackingGeometry());
  TrajectoryStateOnSurface innerMuTSOS = muTT.innermostMeasurementState();
  TrajectoryStateOnSurface outerTkTsos, innerTkTsos;
  if (tkCand.second.isNonnull()) {
    // make sure the tracker track has enough momentum to reach the muon chambers
    if (!(tkCand.second->p() < theMinP || tkCand.second->pt() < theMinPt)) {
      TrajectoryStateOnSurface innerTkTsos;

      outerTkTsos = trajectoryStateTransform::outerStateOnSurface(
          *tkCand.second, *theService->trackingGeometry(), &*theService->magneticField());
      innerTkTsos = trajectoryStateTransform::innerStateOnSurface(
          *tkCand.second, *theService->trackingGeometry(), &*theService->magneticField());
      // for cosmics, outer-most referst to last traversed layer
      if ((innerMuTSOS.globalPosition() - outerTkTsos.globalPosition()).mag() >
          (innerMuTSOS.globalPosition() - innerTkTsos.globalPosition()).mag())
        outerTkTsos = innerTkTsos;
    }
  }

  if (!innerMuTSOS.isValid() || !outerTkTsos.isValid()) {
    LogTrace(category) << "A TSOS validity problem! MuTSOS " << innerMuTSOS.isValid() << " TkTSOS "
                       << outerTkTsos.isValid();
    return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);
  }

  const Surface& refSurface = innerMuTSOS.surface();
  TrajectoryStateOnSurface tkAtMu =
      theService->propagator(theOutPropagatorName)->propagate(*outerTkTsos.freeState(), refSurface);

  if (!tkAtMu.isValid()) {
    LogTrace(category) << "Could not propagate Muon and Tracker track to the same muon hit surface!";
    return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);
  }

  return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(innerMuTSOS, tkAtMu);
}

std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface > GlobalMuonTrackMatcher::convertToTSOSTk	(	const TrackCand &	staCand,
		const TrackCand &	tkCand
	)		const

Definition at line 359 of file GlobalMuonTrackMatcher.cc.

References taus_updatedMVAIds_cff::category, relativeConstraints::empty, TrajectoryStateOnSurface::isValid(), LogTrace, trajectoryStateTransform::outerStateOnSurface(), samePlane(), TrajectoryStateOnSurface::surface(), theMinP, theMinPt, theOutPropagatorName, and theService.

                                                             {
  const string category = "GlobalMuonTrackMatcher";

  TrajectoryStateOnSurface empty;

  TransientTrack muTT(*staCand.second, &*theService->magneticField(), theService->trackingGeometry());
  TrajectoryStateOnSurface impactMuTSOS = muTT.impactPointState();

  TrajectoryStateOnSurface outerTkTsos;
  if (tkCand.second.isNonnull()) {
    // make sure the tracker track has enough momentum to reach the muon chambers
    if (!(tkCand.second->p() < theMinP || tkCand.second->pt() < theMinPt)) {
      outerTkTsos = trajectoryStateTransform::outerStateOnSurface(
          *tkCand.second, *theService->trackingGeometry(), &*theService->magneticField());
    }
  }

  if (!impactMuTSOS.isValid() || !outerTkTsos.isValid())
    return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);

  // define StateOnTrackerBound object
  StateOnTrackerBound fromInside(&*theService->propagator(theOutPropagatorName));

  // extrapolate to outer tracker surface
  TrajectoryStateOnSurface tkTsosFromMu = fromInside(impactMuTSOS);
  TrajectoryStateOnSurface tkTsosFromTk = fromInside(outerTkTsos);

  if (!samePlane(tkTsosFromMu, tkTsosFromTk)) {
    // propagate tracker track to same surface as muon
    bool same1, same2;
    TrajectoryStateOnSurface newTkTsosFromTk, newTkTsosFromMu;
    if (tkTsosFromMu.isValid())
      newTkTsosFromTk = theService->propagator(theOutPropagatorName)->propagate(outerTkTsos, tkTsosFromMu.surface());
    same1 = samePlane(newTkTsosFromTk, tkTsosFromMu);
    LogTrace(category) << "Propagating to same tracker surface (Mu):" << same1;
    if (!same1) {
      if (tkTsosFromTk.isValid())
        newTkTsosFromMu = theService->propagator(theOutPropagatorName)->propagate(impactMuTSOS, tkTsosFromTk.surface());
      same2 = samePlane(newTkTsosFromMu, tkTsosFromTk);
      LogTrace(category) << "Propagating to same tracker surface (Tk):" << same2;
    }
    if (same1)
      tkTsosFromTk = newTkTsosFromTk;
    else if (same2)
      tkTsosFromMu = newTkTsosFromMu;
    else {
      LogTrace(category) << "Could not propagate Muon and Tracker track to the same tracker bound!";
      return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);
    }
  }

  return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(tkTsosFromMu, tkTsosFromTk);
}

std::pair< TrajectoryStateOnSurface, TrajectoryStateOnSurface > GlobalMuonTrackMatcher::convertToTSOSTkHit	(	const TrackCand &	staCand,
		const TrackCand &	tkCand
	)		const

Definition at line 461 of file GlobalMuonTrackMatcher.cc.

References taus_updatedMVAIds_cff::category, relativeConstraints::empty, TrajectoryStateOnSurface::globalPosition(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), LogTrace, mag(), trajectoryStateTransform::outerStateOnSurface(), TrajectoryStateOnSurface::surface(), theMinP, theMinPt, theOutPropagatorName, and theService.

Referenced by match().

                                                             {
  const string category = "GlobalMuonTrackMatcher";

  TrajectoryStateOnSurface empty;

  TransientTrack muTT(*staCand.second, &*theService->magneticField(), theService->trackingGeometry());
  TrajectoryStateOnSurface impactMuTSOS = muTT.impactPointState();
  TrajectoryStateOnSurface innerMuTSOS = muTT.innermostMeasurementState();

  TrajectoryStateOnSurface outerTkTsos, innerTkTsos;
  if (tkCand.second.isNonnull()) {
    // make sure the tracker track has enough momentum to reach the muon chambers
    if (!(tkCand.second->p() < theMinP || tkCand.second->pt() < theMinPt)) {
      outerTkTsos = trajectoryStateTransform::outerStateOnSurface(
          *tkCand.second, *theService->trackingGeometry(), &*theService->magneticField());
      innerTkTsos = trajectoryStateTransform::innerStateOnSurface(
          *tkCand.second, *theService->trackingGeometry(), &*theService->magneticField());

      // for cosmics, outer-most referst to last traversed layer
      if ((innerMuTSOS.globalPosition() - outerTkTsos.globalPosition()).mag() >
          (innerMuTSOS.globalPosition() - innerTkTsos.globalPosition()).mag())
        outerTkTsos = innerTkTsos;
    }
  }

  if (!impactMuTSOS.isValid() || !outerTkTsos.isValid()) {
    LogTrace(category) << "A TSOS validity problem! MuTSOS " << impactMuTSOS.isValid() << " TkTSOS "
                       << outerTkTsos.isValid();
    return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);
  }

  const Surface& refSurface = outerTkTsos.surface();
  TrajectoryStateOnSurface muAtTk =
      theService->propagator(theOutPropagatorName)->propagate(*impactMuTSOS.freeState(), refSurface);

  if (!muAtTk.isValid()) {
    LogTrace(category) << "Could not propagate Muon and Tracker track to the same tracker hit surface!";
    return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(empty, empty);
  }

  return pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface>(muAtTk, outerTkTsos);
}

double GlobalMuonTrackMatcher::match	(	const TrackCand &	sta,
		const TrackCand &	track,
		int	matchOption = 0,
		int	surfaceOption = 1
	)		const

check if two tracks are compatible matchOption: 0 = chi2, 1 = distance, 2 = deltaR surfaceOption: 0 = outermost tracker surface, 1 = innermost muon system surface

Definition at line 100 of file GlobalMuonTrackMatcher.cc.

References convertToTSOSMuHit(), convertToTSOSTkHit(), match_Chi2(), match_d(), match_dist(), and match_Rpos().

Referenced by GlobalCosmicMuonTrajectoryBuilder::match(), GlobalTrackQualityProducer::produce(), GlobalMuonTrajectoryBuilder::trajectories(), and L3MuonTrajectoryBuilder::trajectories().

                                                              {
  std::pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface> tsosPair;
  if (surfaceOption == 0)
    tsosPair = convertToTSOSMuHit(sta, track);
  if (surfaceOption == 1)
    tsosPair = convertToTSOSTkHit(sta, track);
  if (surfaceOption != 0 && surfaceOption != 1)
    return -1.0;

  if (matchOption == 0) {
    // chi^2
    return match_Chi2(tsosPair.first, tsosPair.second);
  } else if (matchOption == 1) {
    // distance
    return match_d(tsosPair.first, tsosPair.second);
  } else if (matchOption == 2) {
    // deltaR
    return match_Rpos(tsosPair.first, tsosPair.second);
  } else if (matchOption == 3) {
    return match_dist(tsosPair.first, tsosPair.second);
  } else {
    return -1.0;
  }
}

vector< GlobalMuonTrackMatcher::TrackCand > GlobalMuonTrackMatcher::match	(	const TrackCand &	sta,
		const std::vector< TrackCand > &	tracks
	)		const

choose all tracks with a matching-chi2 less than Chi2Cut

Definition at line 160 of file GlobalMuonTrackMatcher.cc.

References PV3DBase< T, PVType, FrameType >::barePhi(), HLT_2018_cff::cands, taus_updatedMVAIds_cff::category, hltPixelTracks_cff::chi2, convertToTSOSMuHit(), SiPixelRawToDigiRegional_cfi::deltaPhi, PbPb_ZMuSkimMuonDPG_cff::deltaR, HLT_2018_cff::distance, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), cuy::ii, createfilelist::int, TrajectoryStateOnSurface::isValid(), findQualityFiles::jj, LogTrace, match_Chi2(), match_d(), match_dist(), match_Rpos(), mps_fire::result, theChi2_1, theChi2_2, theDeltaD_1, theDeltaD_2, theDeltaR_1, theDeltaR_2, theEta_threshold, theLocChi2, thePt_threshold1, thePt_threshold2, theQual_1, theQual_2, and theQual_3.

                                                                                                               {
  const string category = "GlobalMuonTrackMatcher";

  vector<TrackCand> result;

  if (tracks.empty())
    return result;

  typedef std::pair<TrackCand, TrajectoryStateOnSurface> TrackCandWithTSOS;
  vector<TrackCandWithTSOS> cands;
  int iiTk = 1;
  TrajectoryStateOnSurface muonTSOS;

  LogTrace(category) << "   ***" << endl << "STA Muon pT " << sta.second->pt();
  LogTrace(category) << "   Tk in Region " << tracks.size() << endl;

  for (vector<TrackCand>::const_iterator is = tracks.begin(); is != tracks.end(); ++is, iiTk++) {
    // propagate to a common surface
    std::pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface> tsosPair = convertToTSOSMuHit(sta, *is);
    LogTrace(category) << "    Tk " << iiTk << " of " << tracks.size() << "  ConvertToMuHitSurface muon isValid "
                       << tsosPair.first.isValid() << " tk isValid " << tsosPair.second.isValid() << endl;
    if (tsosPair.first.isValid())
      muonTSOS = tsosPair.first;
    cands.push_back(TrackCandWithTSOS(*is, tsosPair.second));
  }

  // initialize variables
  double min_chisq = 999999;
  double min_d = 999999;
  double min_de = 999999;
  double min_r_pos = 999999;
  std::vector<bool> passes(cands.size(), false);
  int jj = 0;

  int iTkCand = 1;
  for (vector<TrackCandWithTSOS>::const_iterator ii = cands.begin(); ii != cands.end(); ++ii, jj++, iTkCand++) {
    // tracks that are able not able propagate to a common surface
    if (!muonTSOS.isValid() || !(*ii).second.isValid())
      continue;

    // calculate matching variables
    double distance = match_d(muonTSOS, (*ii).second);
    double chi2 = match_Chi2(muonTSOS, (*ii).second);
    double loc_chi2 = match_dist(muonTSOS, (*ii).second);
    double deltaR = match_Rpos(muonTSOS, (*ii).second);

    LogTrace(category) << "   iTk " << iTkCand << " of " << cands.size() << " eta "
                       << (*ii).second.globalPosition().eta() << " phi " << (*ii).second.globalPosition().phi() << endl;
    LogTrace(category) << "    distance " << distance << " distance cut "
                       << " " << endl;
    LogTrace(category) << "    chi2 " << chi2 << " chi2 cut "
                       << " " << endl;
    LogTrace(category) << "    loc_chi2 " << loc_chi2 << " locChi2 cut "
                       << " " << endl;
    LogTrace(category) << "    deltaR " << deltaR << " deltaR cut "
                       << " " << endl;

    if ((*ii).second.globalMomentum().perp() < thePt_threshold1) {
      LogTrace(category) << "    Enters  a1" << endl;

      if ((chi2 > 0 && fabs((*ii).second.globalMomentum().eta()) < theEta_threshold && chi2 < theChi2_1) ||
          (distance > 0 && distance / (*ii).first.second->pt() < theDeltaD_1 && loc_chi2 > 0 &&
           loc_chi2 < theLocChi2)) {
        LogTrace(category) << "    Passes a1" << endl;
        result.push_back((*ii).first);
        passes[jj] = true;
      }
    }
    if ((passes[jj] == false) && (*ii).second.globalMomentum().perp() < thePt_threshold2) {
      LogTrace(category) << "    Enters a2" << endl;
      if ((chi2 > 0 && chi2 < theChi2_2) || (distance > 0 && distance < theDeltaD_2)) {
        LogTrace(category) << "    Passes a2" << endl;
        result.push_back((*ii).first);
        passes[jj] = true;
      }
    } else {
      LogTrace(category) << "    Enters a3" << endl;
      if (distance > 0 && distance < theDeltaD_3 && deltaR > 0 && deltaR < theDeltaR_1) {
        LogTrace(category) << "    Passes a3" << endl;
        result.push_back((*ii).first);
        passes[jj] = true;
      }
    }

    if (passes[jj]) {
      if (distance < min_d)
        min_d = distance;
      if (loc_chi2 < min_de)
        min_de = loc_chi2;
      if (deltaR < min_r_pos)
        min_r_pos = deltaR;
      if (chi2 < min_chisq)
        min_chisq = chi2;
    }
  }

  // re-initialize mask counter
  jj = 0;

  if (result.empty()) {
    LogTrace(category) << "   Stage 1 returned 0 results";
    for (vector<TrackCandWithTSOS>::const_iterator is = cands.begin(); is != cands.end(); ++is, jj++) {
      double deltaR = match_Rpos(muonTSOS, (*is).second);

      if (muonTSOS.isValid() && (*is).second.isValid()) {
        // check matching between tracker and muon tracks using dEta cut looser then dPhi cut
        LogTrace(category) << "    Stage 2 deltaR " << deltaR << " deltaEta "
                           << fabs((*is).second.globalPosition().eta() - muonTSOS.globalPosition().eta() <
                                   1.5 * theDeltaR_2)
                           << " deltaPhi "
                           << (fabs(deltaPhi((*is).second.globalPosition().barePhi(),
                                             muonTSOS.globalPosition().barePhi())) < theDeltaR_2)
                           << endl;

        if (fabs((*is).second.globalPosition().eta() - muonTSOS.globalPosition().eta()) < 1.5 * theDeltaR_2 &&
            fabs(deltaPhi((*is).second.globalPosition().barePhi(), muonTSOS.globalPosition().barePhi())) <
                theDeltaR_2) {
          result.push_back((*is).first);
          passes[jj] = true;
        }
      }

      if (passes[jj]) {
        double distance = match_d(muonTSOS, (*is).second);
        double chi2 = match_Chi2(muonTSOS, (*is).second);
        double loc_chi2 = match_dist(muonTSOS, (*is).second);
        if (distance < min_d)
          min_d = distance;
        if (loc_chi2 < min_de)
          min_de = loc_chi2;
        if (deltaR < min_r_pos)
          min_r_pos = deltaR;
        if (chi2 < min_chisq)
          min_chisq = chi2;
      }
    }
  }

  for (vector<TrackCand>::const_iterator iTk = result.begin(); iTk != result.end(); ++iTk) {
    LogTrace(category) << "   -----" << endl
                       << "selected pt " << iTk->second->pt() << " eta " << iTk->second->eta() << " phi "
                       << iTk->second->phi() << endl;
  }

  if (result.size() < 2)
    return result;
  else
    result.clear();

  LogTrace(category) << "   Cleaning matched candiates" << endl;

  // re-initialize mask counter
  jj = 0;

  for (vector<TrackCandWithTSOS>::const_iterator is = cands.begin(); is != cands.end(); ++is, jj++) {
    if (!passes[jj])
      continue;

    double distance = match_d(muonTSOS, (*is).second);
    double chi2 = match_Chi2(muonTSOS, (*is).second);
    //unused    double loc_chi2 = match_dist(muonTSOS,(*is).second);
    double deltaR = match_Rpos(muonTSOS, (*is).second);

    // compute quality as the relative ratio to the minimum found for each variable

    int qual = (int)(chi2 / min_chisq + distance / min_d + deltaR / min_r_pos);
    int n_min =
        ((chi2 / min_chisq == 1) ? 1 : 0) + ((distance / min_d == 1) ? 1 : 0) + ((deltaR / min_r_pos == 1) ? 1 : 0);

    if (n_min == 3) {
      result.push_back((*is).first);
    }

    if (n_min == 2 && qual < theQual_1) {
      result.push_back((*is).first);
    }

    if (n_min == 1 && qual < theQual_2) {
      result.push_back((*is).first);
    }

    if (n_min == 0 && qual < theQual_3) {
      result.push_back((*is).first);
    }
  }

  for (vector<TrackCand>::const_iterator iTk = result.begin(); iTk != result.end(); ++iTk) {
    LogTrace(category) << "   -----" << endl
                       << "selected pt " << iTk->second->pt() << " eta " << iTk->second->eta() << " phi "
                       << iTk->second->phi() << endl;
  }

  return result;
}

double GlobalMuonTrackMatcher::match_Chi2	(	const TrajectoryStateOnSurface &	tsos1,
		const TrajectoryStateOnSurface &	tsos2
	)		const

Definition at line 533 of file GlobalMuonTrackMatcher.cc.

References taus_updatedMVAIds_cff::category, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), visualization-live-secondInstance_cfg::m, LocalTrajectoryError::matrix(), findQualityFiles::v, and LocalTrajectoryParameters::vector().

Referenced by match(), matchOne(), and matchTight().

                                                                                       {
  const string category = "GlobalMuonTrackMatcher";
  //LogTrace(category) << "match_Chi2 sanity check: " << tsos1.isValid() << " " << tsos2.isValid();
  if (!tsos1.isValid() || !tsos2.isValid())
    return -1.;

  AlgebraicVector5 v(tsos1.localParameters().vector() - tsos2.localParameters().vector());
  AlgebraicSymMatrix55 m(tsos1.localError().matrix() + tsos2.localError().matrix());

  //LogTrace(category) << "match_Chi2 vector v " << v;

  bool ierr = !m.Invert();

  if (ierr) {
    edm::LogInfo(category) << "Error inverting covariance matrix";
    return -1;
  }

  double est = ROOT::Math::Similarity(v, m);

  //LogTrace(category) << "Chi2 " << est;

  return est;
}

double GlobalMuonTrackMatcher::match_D	(	const TrajectoryStateOnSurface &	sta,
		const TrajectoryStateOnSurface &	tk
	)		const

Definition at line 596 of file GlobalMuonTrackMatcher.cc.

References TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), and mag().

Referenced by samePlane().

                                                                                                                    {
  if (!sta.isValid() || !tk.isValid())
    return -1;
  return (sta.globalPosition() - tk.globalPosition()).mag();
}

double GlobalMuonTrackMatcher::match_d	(	const TrajectoryStateOnSurface &	sta,
		const TrajectoryStateOnSurface &	tk
	)		const

Definition at line 605 of file GlobalMuonTrackMatcher.cc.

References TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), and mag().

Referenced by match(), matchTight(), and samePlane().

                                                                                                                    {
  if (!sta.isValid() || !tk.isValid())
    return -1;
  return (sta.localPosition() - tk.localPosition()).mag();
}

double GlobalMuonTrackMatcher::match_dist	(	const TrajectoryStateOnSurface &	sta,
		const TrajectoryStateOnSurface &	tk
	)		const

Definition at line 615 of file GlobalMuonTrackMatcher.cc.

References taus_updatedMVAIds_cff::category, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), LogTrace, visualization-live-secondInstance_cfg::m, LocalTrajectoryError::positionError(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by match().

                                                                                    {
  const string category = "GlobalMuonTrackMatcher";

  if (!sta.isValid() || !tk.isValid())
    return -1;

  AlgebraicMatrix22 m;
  m(0, 0) = tk.localError().positionError().xx() + sta.localError().positionError().xx();
  m(1, 0) = m(0, 1) = tk.localError().positionError().xy() + sta.localError().positionError().xy();
  m(1, 1) = tk.localError().positionError().yy() + sta.localError().positionError().yy();

  AlgebraicVector2 v;
  v[0] = tk.localPosition().x() - sta.localPosition().x();
  v[1] = tk.localPosition().y() - sta.localPosition().y();

  if (!m.Invert()) {
    LogTrace(category) << "Error inverting local matrix ";
    return -1;
  }

  return ROOT::Math::Similarity(v, m);
}

double GlobalMuonTrackMatcher::match_R_IP	(	const TrackCand &	staCand,
		const TrackCand &	tkCand
	)		const

Definition at line 562 of file GlobalMuonTrackMatcher.cc.

References HGC3DClusterGenMatchSelector_cfi::dR.

                                                                                                 {
  double dR = 99.0;
  if (tkCand.second.isNonnull()) {
    dR = (deltaR<double>(staCand.second->eta(), staCand.second->phi(), tkCand.second->eta(), tkCand.second->phi()));
  }

  return dR;
}

double GlobalMuonTrackMatcher::match_Rmom	(	const TrajectoryStateOnSurface &	sta,
		const TrajectoryStateOnSurface &	tk
	)		const

Definition at line 574 of file GlobalMuonTrackMatcher.cc.

References PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::isValid(), and PV3DBase< T, PVType, FrameType >::phi().

                                                                                    {
  if (!sta.isValid() || !tk.isValid())
    return -1;
  return (deltaR<double>(
      sta.globalMomentum().eta(), sta.globalMomentum().phi(), tk.globalMomentum().eta(), tk.globalMomentum().phi()));
}

double GlobalMuonTrackMatcher::match_Rpos	(	const TrajectoryStateOnSurface &	sta,
		const TrajectoryStateOnSurface &	tk
	)		const

Definition at line 585 of file GlobalMuonTrackMatcher.cc.

References PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), and PV3DBase< T, PVType, FrameType >::phi().

Referenced by match().

                                                                                    {
  if (!sta.isValid() || !tk.isValid())
    return -1;
  return (deltaR<double>(
      sta.globalPosition().eta(), sta.globalPosition().phi(), tk.globalPosition().eta(), tk.globalPosition().phi()));
}

std::vector< GlobalMuonTrackMatcher::TrackCand >::const_iterator GlobalMuonTrackMatcher::matchOne	(	const TrackCand &	sta,
		const std::vector< TrackCand > &	tracks
	)		const

choose the one tracker track which best matches a muon track

Definition at line 132 of file GlobalMuonTrackMatcher.cc.

References hltPixelTracks_cff::chi2, convertToTSOSMuHit(), match_Chi2(), and mps_fire::result.

                                                                    {
  if (tracks.empty())
    return tracks.end();

  double minChi2 = 1000.0;
  vector<TrackCand>::const_iterator result = tracks.end();
  for (vector<TrackCand>::const_iterator is = tracks.begin(); is != tracks.end(); ++is) {
    // propagate to common surface
    std::pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface> tsosPair = convertToTSOSMuHit(sta, *is);

    // calculate chi^2 of local track parameters
    double chi2 = match_Chi2(tsosPair.first, tsosPair.second);
    if (chi2 > 0. && chi2 <= minChi2) {
      minChi2 = chi2;
      result = is;
    }
  }

  return result;
}

bool GlobalMuonTrackMatcher::matchTight	(	const TrackCand &	sta,
		const TrackCand &	track
	)		const

check if two tracks are compatible (less than Chi2Cut, DeltaDCut, DeltaRCut)

Definition at line 80 of file GlobalMuonTrackMatcher.cc.

References hltPixelTracks_cff::chi2, convertToTSOSMuHit(), HLT_2018_cff::distance, match_Chi2(), match_d(), theChi2_2, and theDeltaD_2.

Referenced by GlobalTrackQualityProducer::produce().

                                                                                          {
  std::pair<TrajectoryStateOnSurface, TrajectoryStateOnSurface> tsosPair = convertToTSOSMuHit(sta, track);

  double chi2 = match_Chi2(tsosPair.first, tsosPair.second);
  if (chi2 > 0. && chi2 < theChi2_2)
    return true;

  double distance = match_d(tsosPair.first, tsosPair.second);
  if (distance > 0. && distance < theDeltaD_2)
    return true;

  //double deltaR = match_Rpos(tsosPair.first,tsosPair.second);
  //if ( deltaR > 0. && deltaR < theDeltaR_3 ) return true;

  return false;
}

bool GlobalMuonTrackMatcher::samePlane	(	const TrajectoryStateOnSurface &	tsos1,
		const TrajectoryStateOnSurface &	tsos2
	)		const

Definition at line 508 of file GlobalMuonTrackMatcher.cc.

References TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), match_D(), match_d(), p1, p2, TrajectoryStateOnSurface::surface(), and Surface::tangentPlane().

Referenced by convertToTSOSTk().

                                                                                    {
  if (!tsos1.isValid() || !tsos2.isValid())
    return false;

  if (fabs(match_D(tsos1, tsos2) - match_d(tsos1, tsos2)) > 0.1)
    return false;

  const float maxtilt = 0.999;
  const float maxdist = 0.01;  // in cm

  auto p1(tsos1.surface().tangentPlane(tsos1.localPosition()));
  auto p2(tsos2.surface().tangentPlane(tsos2.localPosition()));

  bool returnValue = ((fabs(p1->normalVector().dot(p2->normalVector())) > maxtilt) ||
                      (fabs((p1->toLocal(p2->position())).z()) < maxdist))
                         ? true
                         : false;

  return returnValue;
}

Member Data Documentation

double GlobalMuonTrackMatcher::theChi2_1

private

Definition at line 85 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theChi2_2

private

Definition at line 86 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), match(), and matchTight().

double GlobalMuonTrackMatcher::theChi2_3

private

Definition at line 87 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher().

double GlobalMuonTrackMatcher::theDeltaD_1

private

Definition at line 89 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theDeltaD_2

private

Definition at line 90 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), match(), and matchTight().

double GlobalMuonTrackMatcher::theDeltaD_3

private

Definition at line 91 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher().

double GlobalMuonTrackMatcher::theDeltaR_1

private

Definition at line 92 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theDeltaR_2

private

Definition at line 93 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theDeltaR_3

private

Definition at line 94 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher().

double GlobalMuonTrackMatcher::theEta_threshold

private

Definition at line 84 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theLocChi2

private

Definition at line 88 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theMinP

private

Definition at line 80 of file GlobalMuonTrackMatcher.h.

Referenced by convertToTSOSMuHit(), convertToTSOSTk(), convertToTSOSTkHit(), and GlobalMuonTrackMatcher().

double GlobalMuonTrackMatcher::theMinPt

private

Definition at line 81 of file GlobalMuonTrackMatcher.h.

Referenced by convertToTSOSMuHit(), convertToTSOSTk(), convertToTSOSTkHit(), and GlobalMuonTrackMatcher().

std::string GlobalMuonTrackMatcher::theOutPropagatorName

private

Definition at line 100 of file GlobalMuonTrackMatcher.h.

Referenced by convertToTSOSMuHit(), convertToTSOSTk(), convertToTSOSTkHit(), and GlobalMuonTrackMatcher().

double GlobalMuonTrackMatcher::thePt_threshold1

private

Definition at line 82 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::thePt_threshold2

private

Definition at line 83 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theQual_1

private

Definition at line 95 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theQual_2

private

Definition at line 96 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

double GlobalMuonTrackMatcher::theQual_3

private

Definition at line 97 of file GlobalMuonTrackMatcher.h.

Referenced by GlobalMuonTrackMatcher(), and match().

const MuonServiceProxy* GlobalMuonTrackMatcher::theService

private

Definition at line 99 of file GlobalMuonTrackMatcher.h.

Referenced by convertToTSOSMuHit(), convertToTSOSTk(), and convertToTSOSTkHit().