Belle II Software development
DriftLengthBasedEventTimeExtractor Class Referencefinalabstract

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

#include <DriftLengthBasedEventTimeExtractor.h>

Inheritance diagram for DriftLengthBasedEventTimeExtractor:
BaseEventTimeExtractor< RecoTrack * > Findlet< AIOTypes... > CompositeProcessingSignalListener ProcessingSignalListener

Public Types

using IOTypes
 Types that should be served to apply on invocation.
 
using IOTypes
 Types that should be served to apply on invocation.
 
using IOVectors
 Vector types that should be served to apply on invocation.
 
using IOVectors
 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.
 
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 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 initialize () override
 Initialize the event t0 store obj ptr.
 
virtual void beginEvent () override
 Create 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 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
 Short hand for ToRangeImpl.
 
using ToVector
 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 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< EventT0m_eventT0
 Pointer to the storage of the eventwise T0 estimation in the data store.
 
StoreObjPtr< EventT0m_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.
 
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.
 
bool m_param_overwriteExistingEstimation
 Parameter if it is fine to overwrite the current EventT0.
 
std::optional< EventT0::EventT0Componentm_eventT0Before
 Storage for the event t0 before this module ran.
 
std::optional< EventT0::EventT0Componentm_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
 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.
 

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 [1/2]

using IOTypes
inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOTypes [2/2]

using IOTypes
inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOVectors [1/2]

using IOVectors
inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ IOVectors [2/2]

using IOVectors
inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ Super

Type of the base class.

Definition at line 33 of file DriftLengthBasedEventTimeExtractor.h.

◆ ToVector [1/2]

using ToVector
protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

◆ ToVector [2/2]

using ToVector
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 53 of file CompositeProcessingSignalListener.cc.

56{
58}
Class to extract the event t0.

◆ apply()

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

Timing extraction for this findlet.

Definition at line 74 of file DriftLengthBasedEventTimeExtractor.cc.

75{
76 TrackFitter trackFitter;
77 m_wasSuccessful = false;
78
79
80 // TODO: in principle this should also be possible for CDCTracks
81 // Collect all event t0 hypothesis with their weights
82 std::vector<std::pair<double, double>> eventT0WithWeights;
83
84 for (RecoTrack* recoTrack : recoTracks) {
85 if (not trackFitter.fit(*recoTrack)) {
86 continue;
87 }
88
89 const std::vector<CDCHit*>& cdcHits = recoTrack->getSortedCDCHitList();
90 for (CDCHit* cdcHit : cdcHits) {
91 RecoHitInformation* recoHitInformation = recoTrack->getRecoHitInformation(cdcHit);
92 try {
93 const genfit::MeasuredStateOnPlane& mSoP = recoTrack->getMeasuredStateOnPlaneFromRecoHit(recoHitInformation);
94 addEventT0AndWeight(mSoP, *cdcHit, *recoTrack, recoHitInformation, eventT0WithWeights);
95 } catch (...) {
96 continue;
97 }
98 }
99 }
100
101 // calculate the weighted median
102 std::sort(eventT0WithWeights.begin(), eventT0WithWeights.end(), [](const auto & lhs, const auto & rhs) {
103 return lhs.first < rhs.first;
104 });
105
106 double extractedEventT0 = NAN;
107
108 double weightSum = 0;
109 for (const auto& pair : eventT0WithWeights) {
110 weightSum += pair.second;
111 }
112
113 double cumSum = 0;
114 for (const auto& pair : eventT0WithWeights) {
115 cumSum += pair.second;
116 if (cumSum > 0.5 * weightSum) {
117 extractedEventT0 = pair.first;
118 break;
119 }
120 }
121
122 if (not std::isnan(extractedEventT0)) {
123 EventT0::EventT0Component eventT0Component(extractedEventT0, NAN, Const::CDC, "drift length");
124 m_eventT0->setEventT0(eventT0Component);
125 m_wasSuccessful = true;
126 B2DEBUG(25, "Drift length gave a result of " << extractedEventT0);
127 } else {
128 B2DEBUG(25, "Extracted event t0 is nan.");
129 }
130}
bool fit(RecoTrack &recoTrack, genfit::AbsTrackRep *trackRepresentation, bool resortHits=false) const
Fit a reco track with a given non-default track representation.

◆ beginEvent() [1/2]

virtual void beginEvent ( )
overridevirtualinherited

Create the event t0 store obj ptr.

Reimplemented from ProcessingSignalListener.

◆ beginEvent() [2/2]

void beginEvent ( )
overridevirtualinherited

Create the event t0 store obj ptr.

Reimplemented from ProcessingSignalListener.

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

46 {
47 Super::beginEvent();
48
49 m_wasSuccessful = false;
50
51 if (not m_eventT0.isValid()) {
52 m_eventT0.create();
53 }
54
55 m_eventT0Before = m_eventT0->getEventT0Component();
56 }

◆ beginRun()

void beginRun ( )
overrideinherited

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

Definition at line 33 of file CompositeProcessingSignalListener.cc.

◆ endRun()

void endRun ( )
overrideinherited

Receive and dispatch signal for the end of the run.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

◆ exposeParameters() [1/2]

◆ exposeParameters() [2/2]

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

◆ getDescription() [1/2]

virtual std::string getDescription ( )
inlinevirtualinherited

Brief description of the purpose of the concrete findlet.

Definition at line 60 of file Findlet.h.

61 {
62 return "(no description)";
63 }

◆ getDescription() [2/2]

virtual std::string getDescription ( )
inlinevirtualinherited

Brief description of the purpose of the concrete findlet.

Definition at line 60 of file Findlet.h.

61 {
62 return "(no description)";
63 }

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )
protectedinherited

Get the number of currently registered listeners.

Definition at line 56 of file CompositeProcessingSignalListener.cc.

61{
63}

◆ initialize() [1/2]

virtual void initialize ( )
overridevirtualinherited

Initialize the event t0 store obj ptr.

Reimplemented from ProcessingSignalListener.

◆ initialize() [2/2]

void initialize ( )
overridevirtualinherited

Initialize the event t0 store obj ptr.

Reimplemented from ProcessingSignalListener.

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

39 {
40 Super::initialize();
41 m_eventT0.registerInDataStore();
42 }

◆ resetEventT0()

void resetEventT0 ( ) const
protectedinherited

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

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

60 {
61 if (m_eventT0Before) {
62 m_eventT0->setEventT0(*m_eventT0Before);
63 } else {
64 m_eventT0->clearEventT0();
65 }
66 }

◆ terminate()

void terminate ( )
overrideinherited

Receive and dispatch Signal for termination of the event processing.

Definition at line 42 of file CompositeProcessingSignalListener.cc.

48{
50 psl->terminate();
51 }
53}
void terminate() override
Receive and dispatch Signal for termination of the event processing.

◆ wasSuccessful()

bool wasSuccessful ( ) const
inherited

Returns true if the last run t0 extraction was successful.

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

22 {
23 return m_wasSuccessful;
24 }

Member Data Documentation

◆ m_eventT0 [1/2]

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_eventT0 [2/2]

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 [1/2]

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_eventT0Before [2/2]

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
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 [1/2]

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_param_overwriteExistingEstimation [2/2]

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

◆ m_wasSuccessful [1/2]

bool m_wasSuccessful
protectedinherited

Variable to show that the execution was successful.

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

◆ m_wasSuccessful [2/2]

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: