Event time extraction based on the principle of the CDC drift time calculation. More...

#include <DriftLengthBasedEventTimeExtractor.h>

Inheritance diagram for DriftLengthBasedEventTimeExtractor:

Public Types
using	IOTypes = std::tuple< AIOTypes... >
	Types that should be served to apply on invocation.

using	IOVectors = std::tuple< std::vector< AIOTypes >... >
	Vector types that should be served to apply on invocation.

Public Member Functions
void	apply (std::vector< RecoTrack * > &) override final
	Timing extraction for this findlet.

bool	wasSuccessful () const
	Returns true if the last run t0 extraction was successful.

virtual void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) override
	Expose our parameters to the super module.

virtual void	initialize () override
	Initialize the event t0 store obj ptr.

virtual void	beginEvent () override
	Create the event t0 store obj ptr.

virtual std::string	getDescription ()
	Brief description of the purpose of the concrete findlet.

virtual void	apply (ToVector< AIOTypes > &... ioVectors)=0
	Main function executing the algorithm.

void	beginRun () override
	Receive and dispatch signal for the beginning of a new run.

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 Types
using	ToVector = typename ToVectorImpl< T >::Type
	Short hand for ToRangeImpl.

Protected Member Functions
void	resetEventT0 () const
	Reset the t0 value to cached value if it exists or clear it otherwise.

void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.

int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Protected Attributes
StoreObjPtr< EventT0 >	m_eventT0
	Pointer to the storage of the eventwise T0 estimation in the data store.

bool	m_wasSuccessful
	Variable to show that the execution was successful.

Private Types
using	Super = BaseEventTimeExtractor< RecoTrack * >
	Type of the base class.

Private Attributes
bool	m_param_overwriteExistingEstimation
	Parameter if it is fine to overwrite the current EventT0.

std::optional< EventT0::EventT0Component >	m_eventT0Before
	Storage for the event t0 before this module ran.

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 happened before.

bool	m_terminated = false
	Flag to keep track whether termination happened before.

std::string	m_initializedAs
	Name of the type during initialisation.

Detailed Description

Event time extraction based on the principle of the CDC drift time calculation.

The measured time is given by: T_meas = T_drift + T_prop + T_0 + T_flight + T_walk As all terms except T_0 are known from the track fit, it is possible to use this information to extract the T0 from this.

Will set a single candidate and the event t0, if successful. Assumes the tracks to already be fitted with the current event0. Will assume an eventt0 of 0, if none is set.

Definition at line 30 of file DriftLengthBasedEventTimeExtractor.h.

Member Typedef Documentation

◆ IOTypes

using IOTypes = std::tuple<AIOTypes...>

inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOVectors

using IOVectors = std::tuple< std::vector<AIOTypes>... >

inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ Super

using Super = BaseEventTimeExtractor<RecoTrack*>

private

Type of the base class.

Definition at line 33 of file DriftLengthBasedEventTimeExtractor.h.

◆ ToVector

using ToVector = typename ToVectorImpl<T>::Type

protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Member Function Documentation

◆ addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener * psl )

protectedinherited

Register a processing signal listener to be notified.

Definition at line 55 of file CompositeProcessingSignalListener.cc.

{
  m_subordinaryProcessingSignalListeners.push_back(psl);
}

◆ apply()

void apply ( std::vector< RecoTrack * > & recoTracks )

finaloverride

Timing extraction for this findlet.

Definition at line 74 of file DriftLengthBasedEventTimeExtractor.cc.

{
  TrackFitter trackFitter;
  m_wasSuccessful = false;
 
 
  // TODO: in principle this should also be possible for CDCTracks
  // Collect all event t0 hypothesis with their weights
  std::vector<std::pair<double, double>> eventT0WithWeights;
 
  for (RecoTrack* recoTrack : recoTracks) {
    if (not trackFitter.fit(*recoTrack)) {
      continue;
    }
 
    const std::vector<CDCHit*>& cdcHits = recoTrack->getSortedCDCHitList();
    for (CDCHit* cdcHit : cdcHits) {
      RecoHitInformation* recoHitInformation = recoTrack->getRecoHitInformation(cdcHit);
      try {
        const genfit::MeasuredStateOnPlane& mSoP = recoTrack->getMeasuredStateOnPlaneFromRecoHit(recoHitInformation);
        addEventT0AndWeight(mSoP, *cdcHit, *recoTrack, recoHitInformation, eventT0WithWeights);
      } catch (...) {
        continue;
      }
    }
  }
 
  // calculate the weighted median
  std::sort(eventT0WithWeights.begin(), eventT0WithWeights.end(), [](const auto & lhs, const auto & rhs) {
    return lhs.first < rhs.first;
  });
 
  double extractedEventT0 = NAN;
 
  double weightSum = 0;
  for (const auto& pair : eventT0WithWeights) {
    weightSum += pair.second;
  }
 
  double cumSum = 0;
  for (const auto& pair : eventT0WithWeights) {
    cumSum += pair.second;
    if (cumSum > 0.5 * weightSum) {
      extractedEventT0 = pair.first;
      break;
    }
  }
 
  if (not std::isnan(extractedEventT0)) {
    EventT0::EventT0Component eventT0Component(extractedEventT0, NAN, Const::CDC, "drift length");
    m_eventT0->setEventT0(eventT0Component);
    m_wasSuccessful = true;
    B2DEBUG(25, "Drift length gave a result of " << extractedEventT0);
  } else {
    B2DEBUG(25, "Extracted event t0 is nan.");
  }
}

◆ beginEvent()

void beginEvent

overridevirtualinherited

Create the event t0 store obj ptr.

Reimplemented from CompositeProcessingSignalListener.

Definition at line 39 of file BaseEventTimeExtractor.icc.h.

  {
    Super::beginEvent();
 
    m_wasSuccessful = false;
 
    if (not m_eventT0.isValid()) {
      m_eventT0.create();
    }
 
    m_eventT0Before = m_eventT0->getEventT0Component();
  }

◆ beginRun()

void beginRun ( )

overridevirtualinherited

Receive and dispatch signal for the beginning of a new run.

Reimplemented from ProcessingSignalListener.

Reimplemented in LayerRelationFilter< AFilter >, FourHitFilter, QualityIndicatorFilter, ThreeHitFilter, TwoHitVirtualIPFilter, TwoHitVirtualIPQIFilter, RecoTrackStorer, ROIFinder, SpacePointLoaderAndPreparer, and TrackCandidateResultRefiner.

Definition at line 23 of file CompositeProcessingSignalListener.cc.

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

◆ endRun()

void endRun ( )

overridevirtualinherited

Receive and dispatch signal for the end of the run.

Reimplemented from ProcessingSignalListener.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

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

◆ exposeParameters()

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

overridevirtualinherited

Expose our parameters to the super module.

Reimplemented from Findlet< AIOTypes... >.

Reimplemented in IterativeEventTimeExtractor< AFindlet >, IterativeEventTimeExtractor< Chi2BasedEventTimeExtractor >, IterativeEventTimeExtractor< DriftLengthBasedEventTimeExtractor >, Chi2BasedEventTimeExtractor, FullGridChi2TrackTimeExtractor, FullGridDriftLengthTrackTimeExtractor, GridEventTimeExtractor< AFindlet >, GridEventTimeExtractor< Belle2::Chi2BasedEventTimeExtractor >, GridEventTimeExtractor< Belle2::DriftLengthBasedEventTimeExtractor >, IterativeChi2BasedEventTimeExtractor, and IterativeDriftLengthBasedEventTimeExtractor.

Definition at line 33 of file BaseEventTimeExtractor.icc.h.

  {
    moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "overwriteExistingEstimation"),
                                  m_param_overwriteExistingEstimation,
                                  "Is it fine to overwrite the current EventT0?",
                                  m_param_overwriteExistingEstimation);
 
    Super::exposeParameters(moduleParamList, prefix);
  }

◆ getDescription()

virtual std::string getDescription ( )

inlinevirtualinherited

Brief description of the purpose of the concrete findlet.

Reimplemented in HitBasedT0Extractor.

Definition at line 60 of file Findlet.h.

      {
        return "(no description)";
      }

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )

protectedinherited

Get the number of currently registered listeners.

Definition at line 60 of file CompositeProcessingSignalListener.cc.

{
  return m_subordinaryProcessingSignalListeners.size();
}

◆ initialize()

void initialize

overridevirtualinherited

Initialize the event t0 store obj ptr.

Reimplemented from CompositeProcessingSignalListener.

Definition at line 36 of file BaseEventTimeExtractor.icc.h.

  {
    Super::initialize();
    m_eventT0.registerInDataStore();
  }

◆ resetEventT0()

void resetEventT0

protectedinherited

Reset the t0 value to cached value if it exists or clear it otherwise.

Definition at line 43 of file BaseEventTimeExtractor.icc.h.

  {
    if (m_eventT0Before) {
      m_eventT0->setEventT0(*m_eventT0Before);
    } else {
      m_eventT0->clearEventT0();
    }
  }

◆ terminate()

void terminate ( )

overridevirtualinherited

Receive and dispatch Signal for termination of the event processing.

Reimplemented from ProcessingSignalListener.

Reimplemented in StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::HyperHough >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::QuadraticLegendre >, and StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::Z0TanLambdaLegendre >.

Definition at line 47 of file CompositeProcessingSignalListener.cc.

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

◆ wasSuccessful()

bool wasSuccessful

inherited

Returns true if the last run t0 extraction was successful.

Definition at line 30 of file BaseEventTimeExtractor.icc.h.

  {
    return m_wasSuccessful;
  }

Member Data Documentation

◆ m_eventT0

StoreObjPtr<EventT0> m_eventT0

protectedinherited

Pointer to the storage of the eventwise T0 estimation in the data store.

Definition at line 46 of file BaseEventTimeExtractor.dcl.h.

◆ m_eventT0Before

std::optional<EventT0::EventT0Component> m_eventT0Before

privateinherited

Storage for the event t0 before this module ran.

Definition at line 56 of file BaseEventTimeExtractor.dcl.h.

◆ m_initialized

bool m_initialized = false

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_param_overwriteExistingEstimation

bool m_param_overwriteExistingEstimation

privateinherited

Parameter if it is fine to overwrite the current EventT0.

Definition at line 53 of file BaseEventTimeExtractor.dcl.h.

◆ m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners

privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 52 of file CompositeProcessingSignalListener.h.

◆ m_terminated

bool m_terminated = false

privateinherited

Flag to keep track whether termination happened before.

Definition at line 55 of file ProcessingSignalListener.h.

◆ m_wasSuccessful

bool m_wasSuccessful

protectedinherited

Variable to show that the execution was successful.

Definition at line 49 of file BaseEventTimeExtractor.dcl.h.

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

tracking/eventTimeExtraction/findlets/include/DriftLengthBasedEventTimeExtractor.h
tracking/eventTimeExtraction/findlets/src/DriftLengthBasedEventTimeExtractor.cc

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Types

Private Attributes

Detailed Description

Member Typedef Documentation

◆ IOTypes

◆ IOVectors

◆ Super

◆ ToVector

Member Function Documentation

◆ addProcessingSignalListener()

◆ apply()

◆ beginEvent()

◆ beginRun()

◆ endRun()

◆ exposeParameters()

◆ getDescription()

◆ getNProcessingSignalListener()

◆ initialize()

◆ resetEventT0()

◆ terminate()

◆ wasSuccessful()

Member Data Documentation

◆ m_eventT0

◆ m_eventT0Before

◆ m_initialized

◆ m_initializedAs

◆ m_param_overwriteExistingEstimation

◆ m_subordinaryProcessingSignalListeners

◆ m_terminated

◆ m_wasSuccessful