8#include <tracking/ckf/cdc/filters/states/ExtrapolateAndUpdateCDCStateFilter.h>
10#include <tracking/trackFindingCDC/numerics/ERightLeft.h>
11#include <cdc/dataobjects/CDCRecoHit.h>
13#include <tracking/ckf/cdc/entities/CDCCKFState.h>
14#include <tracking/ckf/cdc/entities/CDCCKFPath.h>
16#include <tracking/trackFindingCDC/utilities/StringManipulation.h>
17#include <framework/core/ModuleParamList.h>
22 TrackFindingCDC::ERightLeft setRLInfo(
const genfit::MeasuredStateOnPlane& mSoP,
CDCCKFState& state)
25 const auto& wire = state.getWireHit()->getWire();
28 const auto& hitPosition = wire.getWirePos3DAtZ(trackPosition.z());
31 TrackFindingCDC::ERightLeft rlInfo = trackPosToWire.
xy().
isRightOrLeftOf(mom.xy());
33 state.setRLinfo(rlInfo);
38ExtrapolateAndUpdateCDCStateFilter::ExtrapolateAndUpdateCDCStateFilter()
47 moduleParamList->
getParameter<std::string>(TrackFindingCDC::prefixed(prefix,
"direction")).setDefaultValue(
"forward");
56 genfit::MeasuredStateOnPlane mSoP = lastState.
getTrackState();
59 CDCRecoHit recoHit(wireHit->getHit(),
nullptr);
68 B2ASSERT(
"Should be exactly two measurements", measurements.size() == 2);
70 const auto& rightLeft = setRLInfo(mSoP, state);
72 if (rightLeft == TrackFindingCDC::ERightLeft::c_Right) {
73 state.setChi2(
m_updater.kalmanStep(mSoP, *(measurements[1])));
75 state.setChi2(
m_updater.kalmanStep(mSoP, *(measurements[0])));
78 delete measurements[0];
79 delete measurements[1];
81 state.setTrackState(mSoP);
82 setRLInfo(mSoP, state);
84 return 1. / state.getChi2();
85 }
catch (
const genfit::Exception& e) {
Define states for CKF algorithm, which can be seed track or CDC wire hit.
const genfit::MeasuredStateOnPlane & getTrackState() const
Get genfit track state (but first check if already present)
This class is used to transfer CDC information to the track fit.
std::vector< genfit::MeasurementOnPlane * > constructMeasurementsOnPlane(const genfit::StateOnPlane &state) const override
build MeasurementsOnPlane
genfit::SharedPlanePtr constructPlane(const genfit::StateOnPlane &state) const override
Methods that actually interface to Genfit.
TrackFindingCDC::Weight operator()(const BaseCDCStateFilter::Object &pair) final
Extrapolate along the path (pair.first) to the CDC wireHit-state (pair.second). Return 1/chi2 if Ok,...
KalmanStepper< 1 > m_updater
Kalman filter updater.
void exposeParameters(ModuleParamList *moduleParamList, const std::string &prefix) override
Expose the parameters of the sub findlets.
Advancer m_extrapolator
Kalman filter extrapolator.
The Module parameter list class.
Class representing a hit wire in the central drift chamber.
void addProcessingSignalListener(ProcessingSignalListener *psl)
std::pair< const CDCCKFPath *, CDCCKFState * > Object
ERightLeft isRightOrLeftOf(const Vector2D &rhs) const
Indicates if the given vector is more left or more right if you looked in the direction of this vecto...
A three dimensional vector.
const Vector2D & xy() const
Getter for the xy projected vector ( reference ! )
ModuleParam< T > & getParameter(const std::string &name) const
Returns a reference to a parameter.
std::vector< CDCCKFState > CDCCKFPath
Shortcut for the collection of CDC CKF-algorithm states.
Abstract base class for different kinds of events.