SegmentCreatorMCTruth Class Reference

Findlet that generates segments from wire hits using the mc truth information. More...

#include <SegmentCreatorMCTruth.h>

Inheritance diagram for SegmentCreatorMCTruth:

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	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) final
	Expose the parameters to a module.
std::string	getDescription () final
	Short description of the findlet.
void	initialize () final
	Initialize the Module before event processing.
void	beginEvent () final
	Start processing the current event.
void	apply (const std::vector< CDCWireHit > &inputWireHits, std::vector< CDCSegment2D > &outputSegments) final
	Main function of the segment finding by the cellular automaton.
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	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.
int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Private Types
using	Super = Findlet< const CDCWireHit, CDCSegment2D >
	Type of the base class.

Private Attributes
bool	m_param_reconstructedDriftLength = false
	Parameter : Setup the drift length as it can be estimated from two dimensional information.
bool	m_param_reconstructedPositions = false
	Parameter : Switch to reconstruct the positions in the segments imitating the facet ca picking up all correct hits.
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

Findlet that generates segments from wire hits using the mc truth information.

Definition at line 26 of file SegmentCreatorMCTruth.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 = Findlet<const CDCWireHit, CDCSegment2D>

private

Type of the base class.

Definition at line 30 of file SegmentCreatorMCTruth.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 ( const std::vector< CDCWireHit > &  inputWireHits, std::vector< CDCSegment2D > &  outputSegments  )

final

Main function of the segment finding by the cellular automaton.

Definition at line 67 of file SegmentCreatorMCTruth.cc.

{
  const CDCMCTrackStore& mcTrackStore = CDCMCTrackStore::getInstance();
  const CDCSimHitLookUp& simHitLookUp = CDCSimHitLookUp::getInstance();
 
  using CDCHitVector = CDCMCTrackStore::CDCHitVector;
 
  const std::map<ITrackType, std::vector<CDCHitVector>>& mcSegmentsByMCParticleIdx =
                                                       mcTrackStore.getMCSegmentsByMCParticleIdx();
 
  std::size_t nSegments = 0;
  for (const std::pair<ITrackType, std::vector<CDCHitVector>> mcSegmentsAndMCParticleIdx : mcSegmentsByMCParticleIdx) {
    const std::vector<CDCHitVector>& mcSegments = mcSegmentsAndMCParticleIdx.second;
    nSegments += mcSegments.size();
  }
 
  outputSegments.reserve(outputSegments.size() + nSegments);
  for (const std::pair<ITrackType, std::vector<CDCHitVector>> mcSegmentsAndMCParticleIdx : mcSegmentsByMCParticleIdx) {
 
    const std::vector<CDCHitVector>& mcSegments = mcSegmentsAndMCParticleIdx.second;
    for (const CDCHitVector& mcSegment : mcSegments) {
      outputSegments.push_back(CDCSegment2D());
      CDCSegment2D& segment2D = outputSegments.back();
      for (const CDCHit* ptrHit : mcSegment) {
        const CDCWireHit* wireHit = simHitLookUp.getWireHit(ptrHit, inputWireHits);
        if (not wireHit) continue;
 
        CDCRecoHit2D recoHit2D = simHitLookUp.getClosestPrimaryRecoHit2D(ptrHit, inputWireHits);
        segment2D.push_back(recoHit2D);
      }
      if (segment2D.size() < 3) outputSegments.pop_back();
    }
  }
 
  CDC::RealisticTDCCountTranslator tdcCountTranslator;
  for (CDCSegment2D& segment : outputSegments) {
    for (CDCRecoHit2D& recoHit2D : segment) {
      Vector2D flightDirection = recoHit2D.getFlightDirection2D();
      Vector2D recoPos2D = recoHit2D.getRecoPos2D();
      double alpha = recoPos2D.angleWith(flightDirection);
 
      const CDCWire& wire = recoHit2D.getWire();
      const CDCHit* hit = recoHit2D.getWireHit().getHit();
      const bool rl = recoHit2D.getRLInfo() == ERightLeft::c_Right;
 
      double driftLength = recoHit2D.getRefDriftLength();
      if (m_param_reconstructedDriftLength) {
        // Setup the drift length such that only information
        // that would be present in two dimensional reconstruction is used
        const double beta = 1;
        double flightTimeEstimate = 0;
        FlightTimeEstimator::instance().getFlightTime2D(recoPos2D, alpha, beta);
 
        driftLength =
          tdcCountTranslator.getDriftLength(hit->getTDCCount(),
                                            wire.getWireID(),
                                            flightTimeEstimate,
                                            rl,
                                            wire.getRefZ(),
                                            alpha);
      } else {
        // In case the true drift length should be kept at least smear it with its variance.
        double driftLengthVariance = tdcCountTranslator.getDriftLengthResolution(driftLength,
                                     wire.getWireID(),
                                     rl,
                                     wire.getRefZ(),
                                     alpha);
 
        driftLength += gRandom->Gaus(0, std::sqrt(driftLengthVariance));
      }
      bool snapRecoPos = true;
      recoHit2D.setRefDriftLength(driftLength, snapRecoPos);
    }
  }
 
  if (m_param_reconstructedPositions) {
    for (CDCSegment2D& segment : outputSegments) {
      if (segment.size() > 1) {
        CDCRLWireHitSegment rlWireHitSegment = segment.getRLWireHitSegment();
        segment = CDCSegment2D::reconstructUsingFacets(rlWireHitSegment);
      }
    }
  }
 
  for (CDCSegment2D& segment : outputSegments) {
    segment.receiveISuperCluster();
  }
  std::sort(outputSegments.begin(), outputSegments.end());
}

beginEvent()

void beginEvent ( )

finalvirtual

Start processing the current event.

Reimplemented from ProcessingSignalListener.

Definition at line 61 of file SegmentCreatorMCTruth.cc.

{
  CDCMCManager::getInstance().fill();
  Super::beginEvent();
}

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  )

finalvirtual

Expose the parameters to a module.

Reimplemented from Findlet< const CDCWireHit, CDCSegment2D >.

Definition at line 33 of file SegmentCreatorMCTruth.cc.

{
  moduleParamList->addParameter(prefixed(prefix, "reconstructedDriftLength"),
                                m_param_reconstructedDriftLength,
                                "Switch to assign the reconstructed drift length to each hit, "
                                "as it can be estimated from two dimensional information only.",
                                m_param_reconstructedDriftLength);
 
  moduleParamList->addParameter(prefixed(prefix, "reconstructedPositions"),
                                m_param_reconstructedPositions,
                                "Switch to reconstruct the positions in the segments "
                                "imitating the facet ca picking up all correct hits.",
                                m_param_reconstructedPositions);
}

getDescription()

std::string getDescription ( )

finalvirtual

Short description of the findlet.

Reimplemented from Findlet< const CDCWireHit, CDCSegment2D >.

Definition at line 48 of file SegmentCreatorMCTruth.cc.

{
  return "Constructs segments from wire hits using the mc truth information.";
}

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 ( )

finalvirtual

Initialize the Module before event processing.

Reimplemented from ProcessingSignalListener.

Definition at line 54 of file SegmentCreatorMCTruth.cc.

{
  CDCMCManager::getInstance().requireTruthInformation();
  Super::initialize();
}

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();
}

Member Data Documentation

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_reconstructedDriftLength

bool m_param_reconstructedDriftLength = false

private

Parameter : Setup the drift length as it can be estimated from two dimensional information.

Definition at line 51 of file SegmentCreatorMCTruth.h.

m_param_reconstructedPositions

bool m_param_reconstructedPositions = false

private

Parameter : Switch to reconstruct the positions in the segments imitating the facet ca picking up all correct hits.

Definition at line 54 of file SegmentCreatorMCTruth.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.

---

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

• tracking/trackFindingCDC/findlets/minimal/include/SegmentCreatorMCTruth.h
• tracking/trackFindingCDC/findlets/minimal/src/SegmentCreatorMCTruth.cc