RecoTrackUtil Struct Reference

Structure to summarize utility function to output a list of hits into a RecoTrack. More...

#include <RecoTrackUtil.h>

Static Public Member Functions
static RecoTrack *	storeInto (const TrackingUtilities::CDCTrack &track, StoreArray< RecoTrack > &recoTracks, const double momentumSeedMagnitude)
	For magnetic monopoles; estimates charge sign from all stereo hits, momentum direction from hits in closest superlayer, fixing magnitude with given value. Ignores building CDCTrajectory3D since it is based on a helix.
static RecoTrack *	storeInto (const TrackingUtilities::CDCTrack &track, StoreArray< RecoTrack > &recoTracks)
	Copies the track information to the RecoTrack.
static RecoTrack *	storeInto (const TrackingUtilities::CDCTrajectory3D &traj3D, StoreArray< RecoTrack > &recoTracks)
	Copies the trajectory information to the RecoTrack.
static RecoTrack *	storeInto (const TrackingUtilities::CDCTrajectory3D &traj3D, const double bZ, StoreArray< RecoTrack > &recoTracks)
	Copies the trajectory information to the RecoTrack.
template<class ARLHitHolderRange>
static void	fill (const ARLHitHolderRange &rlWireHitHolders, RecoTrack &recoTrack)
	Translates a range of hits and inserts them in the reco track.

Detailed Description

Structure to summarize utility function to output a list of hits into a RecoTrack.

Definition at line 26 of file RecoTrackUtil.h.

Member Function Documentation

fill()

template<class ARLHitHolderRange>

void fill ( const ARLHitHolderRange & rlWireHitHolders, RecoTrack & recoTrack )

static

Translates a range of hits and inserts them in the reco track.

fill() specialisation for a class that holds right-left wire-hit information

List of explicit specialisations - at the moment only CDCTrack.

Currently only instantiated for RecoTrack.

Parameters

	rlWireHitHolders	A range of hits, which elements support ->getRLWireHit() to access a wire hit including a right left passage hypotheses.
[out]	recoTrack	RecoTrack to be filled.

Definition at line 135 of file RecoTrackUtil.cc.

{
  int sortingParameter = -1;
  for (const auto& rlWireHitHolder : rlWireHitHolders) {
    ++sortingParameter;
 
    const CDCRLWireHit rlWireHit = rlWireHitHolder.getRLWireHit();
    const CDCWireHit& wireHit = rlWireHit.getWireHit();
    const CDCHit* cdcHit = wireHit.getHit();
 
    // Right left ambiguity resolution
    ERightLeft rlInfo = rlWireHit.getRLInfo();
    using RightLeftInformation = RecoHitInformation::RightLeftInformation;
    if (rlInfo == ERightLeft::c_Left) {
      recoTrack.addCDCHit(cdcHit,
                          sortingParameter,
                          RightLeftInformation::c_left,
                          RecoHitInformation::c_CDCTrackFinder);
    } else if (rlInfo == ERightLeft::c_Right) {
      recoTrack.addCDCHit(cdcHit,
                          sortingParameter,
                          RightLeftInformation::c_right,
                          RecoHitInformation::c_CDCTrackFinder);
    } else if (rlInfo == ERightLeft::c_Invalid) {
      recoTrack.addCDCHit(cdcHit,
                          sortingParameter,
                          RightLeftInformation::c_invalidRightLeftInformation,
                          RecoHitInformation::c_CDCTrackFinder);
    } else {
      recoTrack.addCDCHit(cdcHit,
                          sortingParameter,
                          RightLeftInformation::c_undefinedRightLeftInformation,
                          RecoHitInformation::c_CDCTrackFinder);
    }
  }
}

storeInto() [1/4]

RecoTrack * storeInto ( const TrackingUtilities::CDCTrack & track, StoreArray< RecoTrack > & recoTracks )

static

Copies the track information to the RecoTrack.

Definition at line 87 of file RecoTrackUtil.cc.

{
  const CDCTrajectory3D& traj3D = track.getStartTrajectory3D();
  RecoTrack* newRecoTrack = storeInto(traj3D, recoTracks);
  if (not newRecoTrack) return nullptr;
 
  RecoTrackUtil::fill(track, *newRecoTrack);
  return newRecoTrack;
}

storeInto() [2/4]

RecoTrack * storeInto ( const TrackingUtilities::CDCTrack & track, StoreArray< RecoTrack > & recoTracks, const double momentumSeedMagnitude )

static

For magnetic monopoles; estimates charge sign from all stereo hits, momentum direction from hits in closest superlayer, fixing magnitude with given value. Ignores building CDCTrajectory3D since it is based on a helix.

Definition at line 42 of file RecoTrackUtil.cc.

{
  CDCTrack firstHits;
  ISuperLayer closestLayer = track.getStartISuperLayer();
  if (ISuperLayerUtil::isAxial(closestLayer)) closestLayer = ISuperLayerUtil::getNextOutwards(closestLayer);
 
  for (const CDCRecoHit3D& hit : track) {
    if (hit.isAxial()) continue;
    firstHits.push_back(hit);
    if ((hit.getISuperLayer() != closestLayer) and (firstHits.size() > 3)) break;
  }
 
  const CDCSZFitter& szFitter = CDCSZFitter::getFitter();
  const CDCTrajectorySZ& szTrajectory = szFitter.fitWithStereoHits(firstHits);
 
  ROOT::Math::XYZVector position(track.getStartTrajectory3D().getSupport().x(),
                                 track.getStartTrajectory3D().getSupport().y(),
                                 szTrajectory.getZ0());
 
  ROOT::Math::XYZVector momentum(track.getStartTrajectory3D().getFlightDirection3DAtSupport());
  const double z0 = szTrajectory.getZ0();
  const double lambda = std::atan(szTrajectory.getTanLambda());
  momentum *= (std::cos(lambda));
  momentum.SetZ(std::sin(lambda));
  momentum *= momentumSeedMagnitude;
 
  ESign charge;
  double votes(0);
  int nHits(0);
  closestLayer = track.getEndISuperLayer();
  if (ISuperLayerUtil::isAxial(closestLayer)) closestLayer = ISuperLayerUtil::getNextInwards(closestLayer);
  for (auto it = track.rbegin(); it != track.rend(); it++) {
    CDCRecoHit3D hit = *it;
    if (hit.isAxial()) continue;
    votes += (hit.getRecoZ() - z0 - std::tan(lambda) * hit.getArcLength2D());
    nHits++;
    if ((hit.getISuperLayer() != closestLayer) and (nHits > 3)) break;
  }
  votes > 0 ? charge = ESign::c_Plus : charge = ESign::c_Minus;
 
  RecoTrack* newRecoTrack = recoTracks.appendNew(position, momentum, charge);
  RecoTrackUtil::fill(track, *newRecoTrack);
  return newRecoTrack;
}

storeInto() [3/4]

RecoTrack * storeInto ( const TrackingUtilities::CDCTrajectory3D & traj3D, const double bZ, StoreArray< RecoTrack > & recoTracks )

static

Copies the trajectory information to the RecoTrack.

Definition at line 104 of file RecoTrackUtil.cc.

{
  // Set the start parameters
  ROOT::Math::XYZVector position = traj3D.getSupport();
  ROOT::Math::XYZVector momentum =
    bZ == 0 ? traj3D.getFlightDirection3DAtSupport() : traj3D.getMom3DAtSupport(bZ);
  ESign charge = traj3D.getChargeSign();
 
  // Do not propagate invalid fits, signal that the fit is invalid to the caller.
  if (not ESignUtil::isValid(charge) or VectorUtil::hasNAN(momentum) or VectorUtil::hasNAN(position)) {
    return nullptr;
  }
 
  RecoTrack* newRecoTrack = recoTracks.appendNew(position, momentum, charge);
  if (std::isfinite(traj3D.getFlightTime())) {
    newRecoTrack->setTimeSeed(traj3D.getFlightTime());
  }
 
  const CovarianceMatrix<6>& cov6 = traj3D.getCartesianCovariance(bZ);
  TMatrixDSym covSeed = TMatrixConversion::toTMatrix(cov6);
  newRecoTrack->setSeedCovariance(covSeed);
  return newRecoTrack;
}

storeInto() [4/4]

RecoTrack * storeInto ( const TrackingUtilities::CDCTrajectory3D & traj3D, StoreArray< RecoTrack > & recoTracks )

static

Copies the trajectory information to the RecoTrack.

Definition at line 98 of file RecoTrackUtil.cc.

{
  ROOT::Math::XYZVector position = traj3D.getSupport();
  return storeInto(traj3D, CDCBFieldUtil::getBFieldZ(position), recoTracks);
}

---

The documentation for this struct was generated from the following files:

• tracking/trackFindingCDC/eventdata/utils/include/RecoTrackUtil.h
• tracking/trackFindingCDC/eventdata/utils/src/RecoTrackUtil.cc