SimpleSVDStateFilter Class Reference

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

#include <SimpleSVDStateFilter.h>

Inheritance diagram for SimpleSVDStateFilter:

Public Types
using	Object
	Type of the object to be analysed.
using	Object
	Type of the object to be analysed.
using	Interface
	Mark this class as the basic interface.
using	Interface
	Mark this class as the basic interface.

Public Member Functions
void	beginRun () final
	Set the cached B field.
	SimpleSVDStateFilter ()=default
	default constructor to silence cppcheck warnings
TrackFindingCDC::Weight	operator() (const BaseSVDStateFilter::Object &pair) final
	Function to evaluate the object.
virtual bool	needsTruthInformation ()
	Indicates if the filter requires Monte Carlo information.
virtual bool	needsTruthInformation ()
	Indicates if the filter requires Monte Carlo information.
Weight	operator() (const Object *obj)
	Function to evaluate the object.
virtual std::vector< float >	operator() (const std::vector< Object * > &objs)
	Function to evaluate a vector of objects Base implementation applies the function to each object.
virtual Weight	operator() (const Object &obj)
	Function to evaluate the object.
Weight	operator() (const Object *obj)
	Function to evaluate the object.
virtual std::vector< float >	operator() (const std::vector< Object * > &objs)
	Function to evaluate a vector of objects Base implementation applies the function to each object.
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.
virtual void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix)
	Expose the set of parameters of the filter to the module parameter list.

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[4][3]
	Shortcut for a 4x3 array.
using	Super
	Type of the base class.

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

Static Private Attributes
static constexpr const MaximalValueArray	m_param_maximumHelixDistance
	Maximum distance calculated with helix extrapolation in filter 1. Numbers calculated on MC.
static constexpr const MaximalValueArray	m_param_maximumResidual
	Maximum distance calculated with normal extrapolation in filter 3. Numbers calculated on MC.
static constexpr const MaximalValueArray	m_param_maximumResidual2
	Maximum distance calculated with normal extrapolation in filter 2. Numbers calculated on MC.

Detailed Description

A very simple filter for all space points.

Definition at line 19 of file SimpleSVDStateFilter.h.

Member Typedef Documentation

Interface [1/2]

using Interface

inherited

Mark this class as the basic interface.

Definition at line 38 of file Filter.dcl.h.

Interface [2/2]

using Interface

inherited

Mark this class as the basic interface.

Definition at line 38 of file Filter.dcl.h.

MaximalValueArray

using MaximalValueArray = double[4][3]

private

Shortcut for a 4x3 array.

Definition at line 37 of file SimpleSVDStateFilter.h.

Object [1/2]

using Object

inherited

Type of the object to be analysed.

Definition at line 35 of file Filter.dcl.h.

Object [2/2]

using Object

inherited

Type of the object to be analysed.

Definition at line 35 of file Filter.dcl.h.

Super

using Super

privateinherited

Type of the base class.

Definition at line 26 of file CompositeProcessingSignalListener.h.

Member Function Documentation

addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener * psl )

protectedinherited

Register a processing signal listener to be notified.

Definition at line 53 of file CompositeProcessingSignalListener.cc.

{
  m_subordinaryProcessingSignalListeners.push_back(psl);
}

beginEvent()

void beginEvent ( )

overrideinherited

Receive and dispatch signal for the start of a new event.

Definition at line 36 of file CompositeProcessingSignalListener.cc.

{
  Super::beginEvent();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginEvent();
  }
}

beginRun()

void beginRun ( )

finalvirtual

Set the cached B field.

Reimplemented from CompositeProcessingSignalListener.

Definition at line 41 of file SimpleSVDStateFilter.cc.

{
  m_cachedBField = TrackFindingCDC::CDCBFieldUtil::getBFieldZ();
}

endRun()

void endRun ( )

overrideinherited

Receive and dispatch signal for the end of the run.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->endRun();
  }
  Super::endRun();
}

exposeParameters()

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

inlinevirtualinherited

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

Parameters

moduleParamList	Module parameter list to add the new parameters to
prefix	Prefix (or sometimes postfix) to be added to all nely added parameters.

Definition at line 47 of file CompositeProcessingSignalListener.h.

      {};

getNProcessingSignalListener()

int getNProcessingSignalListener ( )

protectedinherited

Get the number of currently registered listeners.

Definition at line 56 of file CompositeProcessingSignalListener.cc.

{
  return m_subordinaryProcessingSignalListeners.size();
}

initialize()

void initialize ( )

overrideinherited

Receive and dispatch signal before the start of the event processing.

Definition at line 30 of file CompositeProcessingSignalListener.cc.

{
  Super::initialize();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->initialize();
  }
}

needsTruthInformation()

bool needsTruthInformation ( )

virtualinherited

Indicates if the filter requires Monte Carlo information.

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

    {
      return false;
    }

operator()() [1/6]

Weight operator() ( const BaseSVDStateFilter::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< std::pair< const std::vector< TrackFindingCDC::WithWeight< const CKFToSVDState * > >, CKFToSVDState * > >.

Definition at line 46 of file SimpleSVDStateFilter.cc.

{
  const std::vector<TrackFindingCDC::WithWeight<const CKFToSVDState*>>& previousStates = pair.first;
  CKFToSVDState* 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.norm();
    maximumValues = &m_param_maximumHelixDistance;
  } else {
    // Filter 2 + 3
    const double residual = m_kalmanStepper.calculateResidual(firstMeasurement, *currentState);
    valueToCheck = residual;
    if (currentState->isFitted()) {
      maximumValues = &m_param_maximumResidual2;
    } else {
      maximumValues = &m_param_maximumResidual;
    }
  }
 
  const unsigned int layer = currentState->getGeometricalLayer();
  if (valueToCheck > (*maximumValues)[layer - 3][getPTRange(momentum)]) {
    return NAN;
  }
 
  return valueToCheck;
}

operator()() [2/6]

Weight operator() ( const Object & obj )

virtualinherited

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.

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

    {
      return 1;
    }

operator()() [3/6]

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.

operator()() [4/6]

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 69 of file Filter.icc.h.

    {
      return obj ? operator()(*obj) : NAN;
    }

operator()() [5/6]

virtual std::vector< float > operator() ( const std::vector< Object * > & objs )

virtualinherited

Function to evaluate a vector of objects Base implementation applies the function to each object.

Can be optimized for MVA filters

Parameters

objs	A vector of pointers to objects

Returns

A vector of float or NAN values. See above

operator()() [6/6]

std::vector< float > operator() ( const std::vector< Object * > & objs )

virtualinherited

Function to evaluate a vector of objects Base implementation applies the function to each object.

Can be optimized for MVA filters

Parameters

objs	A vector of pointers to objects

Returns

A vector of float or NAN values. See above

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

    {
      std::vector<float> out;
      for (const auto& obj : objs) {
        out.push_back(operator()(obj));
      }
      return out;
    }

terminate()

void terminate ( )

overrideinherited

Receive and dispatch Signal for termination of the event processing.

Definition at line 42 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->terminate();
  }
  Super::terminate();
}

Member Data Documentation

m_cachedBField

double m_cachedBField

private

Cache for the B field at the IP.

Definition at line 31 of file SimpleSVDStateFilter.h.

m_initialized

bool m_initialized

privateinherited

Flag to keep track whether initialization happened before.

Definition at line 52 of file ProcessingSignalListener.h.

m_initializedAs

std::string m_initializedAs

privateinherited

Name of the type during initialisation.

Definition at line 58 of file ProcessingSignalListener.h.

m_kalmanStepper

SVDKalmanStepper m_kalmanStepper

private

Kalman stepper (CKF) for SVD.

Definition at line 34 of file SimpleSVDStateFilter.h.

m_param_maximumHelixDistance

const double m_param_maximumHelixDistance

staticconstexprprivate

Initial value:

= {
      
      {1.38532, 1.44346, 1.53386}, 
      {1.96701, 2.17819, 2.19784}, 
      {2.33064, 2.67768, 2.67174}, 
      {6, 6, 6}, 
    }

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

Definition at line 40 of file SimpleSVDStateFilter.h.

                                                                            {
      // high pt, middle pt, low pt
      {1.38532, 1.44346, 1.53386}, // layer 3
      {1.96701, 2.17819, 2.19784}, // layer 4
      {2.33064, 2.67768, 2.67174}, // layer 5
      {6, 6, 6}, // layer 6
    };

m_param_maximumResidual

const double m_param_maximumResidual

staticconstexprprivate

Initial value:

= {
      {1.52328, 1.57714, 1.63069}, 
      {1.61301, 1.78301, 1.75973}, 
      {1.85188, 2.19571, 2.38167}, 
      {5.5, 5, 5}, 
    }

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

Definition at line 49 of file SimpleSVDStateFilter.h.

                                                                       {
      {1.52328, 1.57714, 1.63069}, // layer 3
      {1.61301, 1.78301, 1.75973}, // layer 4
      {1.85188, 2.19571, 2.38167}, // layer 5
      {5.5, 5, 5}, // layer 6
    };

m_param_maximumResidual2

const double m_param_maximumResidual2

staticconstexprprivate

Initial value:

= {
      {0.69713, 0.39872, 0.41854}, 
      {0.72737, 0.40816, 0.53138}, 
      {0.35624, 0.12846, 0.22516}, 
      {0.00975, 0.06432, 0.5969}, 
    }

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

Definition at line 57 of file SimpleSVDStateFilter.h.

                                                                        {
      {0.69713, 0.39872, 0.41854}, // layer 3
      {0.72737, 0.40816, 0.53138}, // layer 4
      {0.35624, 0.12846, 0.22516}, // layer 5
      {0.00975, 0.06432, 0.5969}, // layer 6
    };

m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners

privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 60 of file CompositeProcessingSignalListener.h.

m_terminated

bool m_terminated

privateinherited

Flag to keep track whether termination happened before.

Definition at line 55 of file ProcessingSignalListener.h.

---

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

• tracking/ckf/svd/filters/states/include/SimpleSVDStateFilter.h
• tracking/ckf/svd/filters/states/src/SimpleSVDStateFilter.cc