DriftLengthEstimator Struct Reference

Helper construct implementing the (re)estimation of the drift length for various hit objects. More...

#include <DriftLengthEstimator.h>

Public Member Functions
void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix)
	Add the parameters of the estimator to the module.
double	updateDriftLength (CDCRecoHit2D &recoHit2D)
	Update the drift length of the reconstructed hit in place.
double	updateDriftLength (CDCRecoHit3D &recoHit3D, double tanLambda)
	Update the drift length of the reconstructed hit in place.
void	updateDriftLength (CDCFacet &facet)
	Reestimate the drift length of all three contained drift circles.
void	updateDriftLength (CDCSegment2D &segment)
	Update the drift length of the contained reconstructed hit in place.
void	updateDriftLength (CDCSegment3D &segment, double tanLambda)
	Update the drift length of the contained reconstructed hit in place.
void	updateDriftLength (CDCTrack &track, double tanLambda)
	Update the drift length of the contained reconstructed hit in place.

Public Attributes
bool	m_param_useAlphaInDriftLength = true
	Parameter : Switch to serve the alpha angle to the drift length translator.
double	m_param_tofMassScale = NAN
	Parameter : Mass to estimate the velocity in the flight time to the hit.

Detailed Description

Helper construct implementing the (re)estimation of the drift length for various hit objects.

Definition at line 36 of file DriftLengthEstimator.h.

Member Function Documentation

exposeParameters()

void exposeParameters	(	ModuleParamList *	moduleParamList,
		const std::string &	prefix
	)

Add the parameters of the estimator to the module.

Definition at line 43 of file DriftLengthEstimator.cc.

{
  moduleParamList->addParameter(prefixed(prefix, "useAlphaInDriftLength"),
                                m_param_useAlphaInDriftLength,
                                "Switch to serve the alpha angle to the drift length translator",
                                m_param_useAlphaInDriftLength);
 
  moduleParamList->addParameter(prefixed(prefix, "tofMassScale"),
                                m_param_tofMassScale,
                                "Mass to estimate the velocity in the flight time to the hit",
                                m_param_tofMassScale);
}

updateDriftLength() [1/6]

void updateDriftLength

(

CDCFacet &

facet

)

Reestimate the drift length of all three contained drift circles.

Using the additional flight direction information the accuracy of the drift length can be increased alot helping the filters following this step

Definition at line 88 of file DriftLengthEstimator.cc.

{
  CDC::RealisticTDCCountTranslator tdcCountTranslator;
  const FlightTimeEstimator& flightTimeEstimator = FlightTimeEstimator::instance();
 
  const UncertainParameterLine2D& line = facet.getFitLine();
  Vector2D flightDirection = line->tangential();
  Vector2D centralPos2D = line->closest(facet.getMiddleWire().getRefPos2D());
  double alpha = centralPos2D.angleWith(flightDirection);
  if (not m_param_useAlphaInDriftLength) {
    alpha = 0;
  }
 
  auto doUpdate = [&](CDCRLWireHit & rlWireHit, Vector2D recoPos2D) {
    const CDCWire& wire = rlWireHit.getWire();
    const CDCHit* hit = rlWireHit.getWireHit().getHit();
    const bool rl = rlWireHit.getRLInfo() == ERightLeft::c_Right;
    const double beta = 1;
    double flightTimeEstimate = flightTimeEstimator.getFlightTime2D(recoPos2D, alpha, beta);
    double driftLength = tdcCountTranslator.getDriftLength(hit->getTDCCount(),
                                                           wire.getWireID(),
                                                           flightTimeEstimate,
                                                           rl,
                                                           wire.getRefZ(),
                                                           alpha);
    rlWireHit.setRefDriftLength(driftLength);
  };
 
  doUpdate(facet.getStartRLWireHit(), facet.getStartRecoPos2D());
  doUpdate(facet.getMiddleRLWireHit(), facet.getMiddleRecoPos2D());
  doUpdate(facet.getEndRLWireHit(), facet.getEndRecoPos2D());
 
  // More accurate implementation
  // double startDriftLength = updateDriftLength(facet.getStartRecoHit2D());
  // facet.getStartRLWireHit().setRefDriftLength(startDriftLength);
 
  // double middleDriftLength = updateDriftLength(facet.getMiddleRecoHit2D());
  // facet.getMiddleRLWireHit().setRefDriftLength(middleDriftLength);
 
  // double endDriftLength = updateDriftLength(facet.getEndRecoHit2D());
  // facet.getEndRLWireHit().setRefDriftLength(endDriftLength);
}

updateDriftLength() [2/6]

double updateDriftLength

(

CDCRecoHit2D &

recoHit2D

)

Update the drift length of the reconstructed hit in place.

Definition at line 56 of file DriftLengthEstimator.cc.

{
  CDC::RealisticTDCCountTranslator tdcCountTranslator;
  const FlightTimeEstimator& flightTimeEstimator = FlightTimeEstimator::instance();
 
  Vector2D flightDirection = recoHit2D.getFlightDirection2D();
  Vector2D recoPos2D = recoHit2D.getRecoPos2D();
  double alpha = recoPos2D.angleWith(flightDirection);
  const double beta = 1;
  double flightTimeEstimate = flightTimeEstimator.getFlightTime2D(recoPos2D, alpha, beta);
 
  const CDCWire& wire = recoHit2D.getWire();
  const CDCHit* hit = recoHit2D.getWireHit().getHit();
  const bool rl = recoHit2D.getRLInfo() == ERightLeft::c_Right;
 
  if (not m_param_useAlphaInDriftLength) {
    alpha = 0;
  }
 
  double driftLength = tdcCountTranslator.getDriftLength(hit->getTDCCount(),
                                                         wire.getWireID(),
                                                         flightTimeEstimate,
                                                         rl,
                                                         wire.getRefZ(),
                                                         alpha);
  if (driftLength > -2 and driftLength < 16) {
    bool snapRecoPos = true;
    recoHit2D.setRefDriftLength(driftLength, snapRecoPos);
  }
  return driftLength;
}

updateDriftLength() [3/6]

double updateDriftLength	(	CDCRecoHit3D &	recoHit3D,
		double	tanLambda
	)

Update the drift length of the reconstructed hit in place.

Definition at line 138 of file DriftLengthEstimator.cc.

{
  CDC::RealisticTDCCountTranslator tdcCountTranslator;
  const FlightTimeEstimator& flightTimeEstimator = FlightTimeEstimator::instance();
 
  Vector2D flightDirection = recoHit3D.getFlightDirection2D();
  const Vector3D& recoPos3D = recoHit3D.getRecoPos3D();
  const Vector2D& recoPos2D = recoPos3D.xy();
  double alpha = recoPos2D.angleWith(flightDirection);
  const double beta = 1;
  double flightTimeEstimate = flightTimeEstimator.getFlightTime2D(recoPos2D, alpha, beta);
 
  if (std::isnan(tanLambda)) {
    tanLambda = recoPos3D.z() / recoPos3D.cylindricalR();
  }
  const double theta = M_PI / 2 - std::atan(tanLambda);
  flightTimeEstimate *= hypot2(1, tanLambda);
 
  const CDCWire& wire = recoHit3D.getWire();
  const CDCHit* hit = recoHit3D.getWireHit().getHit();
  const bool rl = recoHit3D.getRLInfo() == ERightLeft::c_Right;
  double driftLength =
    tdcCountTranslator.getDriftLength(hit->getTDCCount(),
                                      wire.getWireID(),
                                      flightTimeEstimate,
                                      rl,
                                      recoPos3D.z(),
                                      alpha,
                                      theta,
                                      hit->getADCCount());
  if (driftLength > -2 and driftLength < 16) {
    bool snapRecoPos = true;
    recoHit3D.setRecoDriftLength(driftLength, snapRecoPos);
  }
  return driftLength;
}

updateDriftLength() [4/6]

void updateDriftLength

(

CDCSegment2D &

segment

)

Update the drift length of the contained reconstructed hit in place.

Definition at line 131 of file DriftLengthEstimator.cc.

{
  for (CDCRecoHit2D& recoHit2D : segment) {
    updateDriftLength(recoHit2D);
  }
}

updateDriftLength() [5/6]

void updateDriftLength	(	CDCSegment3D &	segment,
		double	tanLambda
	)

Update the drift length of the contained reconstructed hit in place.

Definition at line 177 of file DriftLengthEstimator.cc.

{
  for (CDCRecoHit3D& recoHit3D : segment3D) {
    updateDriftLength(recoHit3D, tanLambda);
  }
}

updateDriftLength() [6/6]

void updateDriftLength	(	CDCTrack &	track,
		double	tanLambda
	)

Update the drift length of the contained reconstructed hit in place.

Definition at line 185 of file DriftLengthEstimator.cc.

{
  for (CDCRecoHit3D& recoHit3D : track) {
    updateDriftLength(recoHit3D, tanLambda);
  }
}

Member Data Documentation

m_param_tofMassScale

double m_param_tofMassScale = NAN

Parameter : Mass to estimate the velocity in the flight time to the hit.

Definition at line 67 of file DriftLengthEstimator.h.

m_param_useAlphaInDriftLength

bool m_param_useAlphaInDriftLength = true

Parameter : Switch to serve the alpha angle to the drift length translator.

Definition at line 64 of file DriftLengthEstimator.h.

---

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

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