SimplePXDStateFilter Class Reference

A very simple filter for all space points. More...

#include <SimplePXDStateFilter.h>

Inheritance diagram for SimplePXDStateFilter:

Collaboration diagram for SimplePXDStateFilter:

Public Types
using	Object = AObject
	Type of the object to be analysed.
using	Interface = Filter< AObject >
	Mark this class as the basic interface.

Public Member Functions
void	beginRun () final
	Set the cached B field.
	SimplePXDStateFilter ()=default
	default constructor to silence the cppcheck warnings
TrackFindingCDC::Weight	operator() (const BasePXDStateFilter::Object &pair) final
	Function to evaluate the object. More...
virtual void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix)
	Expose the set of parameters of the filter to the module parameter list. More...
virtual bool	needsTruthInformation ()
	Indicates if the filter requires Monte Carlo information.
Weight	operator() (const Object *obj)
	Function to evaluate the object. More...
void	initialize () override
	Receive and dispatch signal before the start of the event processing.
void	beginEvent () override
	Receive and dispatch signal for the start of a new event.
void	endRun () override
	Receive and dispatch signal for the end of the run.
void	terminate () override
	Receive and dispatch Signal for termination of the event processing.

Protected Member Functions
void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.
int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Private Types
using	MaximalValueArray = double[2][3]
	Shortcut for a 2x3 array.
using	Super = ProcessingSignalListener
	Type of the base class.

Private Attributes
double	m_cachedBField
	Cache for the B field at the IP.
PXDKalmanStepper	m_kalmanStepper
	Kalman stepper (CKF) for PXD.
std::vector< ProcessingSignalListener * >	m_subordinaryProcessingSignalListeners
	References to subordinary signal processing listener contained in this findlet.
bool	m_initialized = false
	Flag to keep track whether initialization happend before.
bool	m_terminated = false
	Flag to keep track whether termination happend before.
std::string	m_initializedAs
	Name of the type during initialisation.

Static Private Attributes
static constexpr const MaximalValueArray	m_param_maximumHelixDistanceXY
	Maximum distance calculated with helix extrapolation in filter 1. Numbers calculated on MC. More...
static constexpr const MaximalValueArray	m_param_maximumResidual
	Maximum distance calculated with normal extrapolation in filter 2. Numbers calculated on MC. More...
static constexpr const MaximalValueArray	m_param_maximumChi2
	Maximum chi^2 in filter 3. Numbers calculated on MC. More...

Detailed Description

A very simple filter for all space points.

Definition at line 29 of file SimplePXDStateFilter.h.

Member Function Documentation

exposeParameters()

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

virtualinherited

Expose the set of parameters of the filter to the module parameter list.

Note that not all filters have yet exposed their parameters in this way.

This method is deprecated as the exposeParams below uses a less compile heavy equivalent.

Reimplemented in ChoosableFromVarSet< Filter< ATruthVarSet::Object > >, ChoosableFromVarSet< Filter< PXDResultVarSet ::Object > >, ChoosableFromVarSet< Filter< SVDResultVarSet ::Object > >, Recording< Filter< AVarSet::Object > >, MVA< Filter< MVAFeasibleSegmentPairVarSet ::Object > >, MVA< BaseAxialSegmentPairFilter >, MVA< Filter< MVARealisticSegmentPairVarSet ::Object > >, MVA< BaseSegmentPairRelationFilter >, MVA< BaseTrackRelationFilter >, MVA< BaseSegmentRelationFilter >, MVA< Filter< BasicClusterVarSet ::Object > >, MVA< BaseFacetRelationFilter >, MVA< BaseFacetFilter >, MVA< Filter< AVarSet::Object > >, AdvanceFilter< AState, AnAdvancer >, DuplicateCDCPathPairFilter, BridgingWireHitRelationFilter, Chi2FacetRelationFilter, SimpleFacetRelationFilter, MCSymmetric< BaseAxialSegmentPairFilter >, MCSymmetric< BaseSegmentPairFilter >, MCSymmetric< BaseSegmentTripleRelationFilter >, MCSymmetric< BaseSegmentPairRelationFilter >, MCSymmetric< BaseTrackRelationFilter >, MCSymmetric< BaseSegmentTripleFilter >, MCSymmetric< BaseSegmentRelationFilter >, MCSymmetric< BaseFacetRelationFilter >, MCSymmetric< BaseFacetFilter >, LayerPXDRelationFilter< AFilter, APrefilter >, ExtrapolateAndUpdateCDCStateFilter, PhiTrackRelationFilter, PhiRecoTrackRelationFilter, RoughCDCfromEclStateFilter, RoughCDCStateFilter, ResidualSVDStateFilter, DistanceCDCStateFilter, LoosePXDPairFilter, SensorPXDPairFilter, LooseSVDPairFilter, SensorSVDPairFilter, Chi2FacetFilter, FeasibleRLFacetFilter, RealisticFacetFilter, SimpleFacetFilter, MCFacetRelationFilter, LayerSVDRelationFilter< AFilter, APrefilter >, MCSegmentTripleFilter, SectorMapBasedSVDPairFilter, WholeWireHitRelationFilter, MCAxialSegmentPairFilter, MCSegmentPairRelationFilter, MCSegmentRelationFilter, MCSegmentTripleRelationFilter, and MCSegmentPairFilter.

Definition at line 42 of file Filter.icc.h.

operator()() [1/2]

Weight operator() ( const BasePXDStateFilter::Object &  obj )

finalvirtual

Function to evaluate the object.

Base implementation accepts all objects.

Parameters

obj	The object to be accepted or rejected.

Returns

A finite float value if the object is accepted. NAN if the object is rejected.

Reimplemented from Filter< AObject >.

Definition at line 48 of file SimplePXDStateFilter.cc.

{
  const std::vector<TrackFindingCDC::WithWeight<const CKFToPXDState*>>& previousStates = pair.first;
  CKFToPXDState* currentState = pair.second;
 
  const auto* spacePoint = currentState->getHit();
 
  genfit::MeasuredStateOnPlane firstMeasurement;
  if (currentState->mSoPSet()) {
    firstMeasurement = currentState->getMeasuredStateOnPlane();
  } else {
    firstMeasurement = previousStates.back()->getMeasuredStateOnPlane();
  }
 
  Vector3D position = Vector3D(firstMeasurement.getPos());
  Vector3D momentum = Vector3D(firstMeasurement.getMom());
 
  const Vector3D hitPosition = static_cast<Vector3D>(spacePoint->getPosition());
 
  const bool sameHemisphere = fabs(position.phi() - hitPosition.phi()) < TMath::PiOver2();
  if (not sameHemisphere) {
    return NAN;
  }
 
  double valueToCheck;
  const MaximalValueArray* maximumValues;
 
  if (not currentState->isFitted() and not currentState->mSoPSet()) {
    // Filter 1
    const RecoTrack* cdcTrack = previousStates.front()->getSeed();
 
    B2ASSERT("A path without a seed?", cdcTrack);
 
    const CDCTrajectory3D trajectory(position, 0, momentum, cdcTrack->getChargeSeed(), m_cachedBField);
 
    const double arcLength = trajectory.calcArcLength2D(hitPosition);
    const Vector2D& trackPositionAtHit2D = trajectory.getTrajectory2D().getPos2DAtArcLength2D(arcLength);
    const double trackPositionAtHitZ = trajectory.getTrajectorySZ().mapSToZ(arcLength);
    const Vector3D trackPositionAtHit(trackPositionAtHit2D, trackPositionAtHitZ);
    const Vector3D differenceHelix = trackPositionAtHit - hitPosition;
 
    valueToCheck = differenceHelix.xy().norm();
    maximumValues = &m_param_maximumHelixDistanceXY;
  } else if (not currentState->isFitted()) {
    // Filter 2
    const double residual = m_kalmanStepper.calculateResidual(firstMeasurement, *currentState);
 
    valueToCheck = residual;
    maximumValues = &m_param_maximumResidual;
  } else {
    // Filter 3
    valueToCheck = currentState->getChi2();
    maximumValues = &m_param_maximumChi2;
  }
 
  const unsigned int layer = currentState->getGeometricalLayer();
  if (valueToCheck > (*maximumValues)[layer - 1][getPTRange(momentum)]) {
    return NAN;
  }
 
  return valueToCheck;
}

operator()() [2/2]

Weight operator() ( const Object *  obj )

inherited

Function to evaluate the object.

Base implementation accepts all objects, except nullptr.

Parameters

obj	The object to be accepted or rejected.

Returns

A finit float value if the object is accepted. NAN if the object is rejected. Nullptr is always rejected.

Definition at line 60 of file Filter.icc.h.

Member Data Documentation

m_param_maximumChi2

constexpr const SimplePXDStateFilter::MaximalValueArray m_param_maximumChi2

staticconstexprprivate

Initial value:

= {
      {10, 50, 125},
      {20, 60, 85}
    }

Maximum chi^2 in filter 3. Numbers calculated on MC.

Definition at line 62 of file SimplePXDStateFilter.h.

m_param_maximumHelixDistanceXY

constexpr const SimplePXDStateFilter::MaximalValueArray m_param_maximumHelixDistanceXY

staticconstexprprivate

Initial value:

= {
      {0.295060, 0.295060, 0.295060},
      {0.295060, 0.295060, 0.295060}
    }

Maximum distance calculated with helix extrapolation in filter 1. Numbers calculated on MC.

Definition at line 50 of file SimplePXDStateFilter.h.

m_param_maximumResidual

constexpr const SimplePXDStateFilter::MaximalValueArray m_param_maximumResidual

staticconstexprprivate

Initial value:

= {
      {0.679690, 0.677690, 0.427060},
      {1.006930, 0.926120, 0.687030}
    }

Maximum distance calculated with normal extrapolation in filter 2. Numbers calculated on MC.

Definition at line 56 of file SimplePXDStateFilter.h.

---

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

• tracking/ckf/pxd/filters/states/include/SimplePXDStateFilter.h
• tracking/ckf/pxd/filters/states/src/SimplePXDStateFilter.cc