CDCMCHitCollectionLookUp< ACDCHitCollection > Class Template Reference

Interface class to the Monte Carlo information for collections of hits. More...

#include <CDCMCHitCollectionLookUp.h>

Public Member Functions
void	clear ()
	Clears all Monte Carlo information left from the last event.
	CDCMCHitCollectionLookUp ()
	Dummy constructor to silence the intel-compiler warnings (NOTE: using "=default" will make the intel compilation fail)
MCTrackIdPurityPair	getHighestPurity (const ACDCHitCollection &hits) const
	Get the track id with the highest corresponding purity.
ITrackType	getMCTrackId (const ACDCHitCollection *ptrHits) const
	Getter for the Monte Carlo track id matched to this collection of hits. More...
int	getCorrectRLVote (const ACDCHitCollection *ptrHits) const
	Getter for the difference of correct versus incorrect right left passage informations.
double	getRLPurity (const ACDCHitCollection *ptrHits) const
	Getter for the right left passge purity which respects the forward backward reconstruction.
const MCParticle *	getMCParticle (const ACDCHitCollection *ptrHits) const
	Getter for the mc particle matched to this collection of hits.
const CDCHit *	getFirstHit (const ACDCHitCollection *ptrHits) const
	Getter for the first hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
const CDCHit *	getLastHit (const ACDCHitCollection *ptrHits) const
	Getter for the last hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getFirstInTrackId (const ACDCHitCollection *ptrHits) const
	Getter for in track id for the first hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getLastInTrackId (const ACDCHitCollection *ptrHits) const
	Getter for in track id for the last hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getFirstInTrackSegmentId (const ACDCHitCollection *ptrHits) const
	Getter for in track collection of hits id for the first hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getLastInTrackSegmentId (const ACDCHitCollection *ptrHits) const
	Getter for in track collection of hits id for the last hit in the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getFirstNPassedSuperLayers (const ACDCHitCollection *ptrHits) const
	Getter for number of passed superlayer till the first hit the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getLastNPassedSuperLayers (const ACDCHitCollection *ptrHits) const
	Getter for number of passed superlayer till the last hit the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getFirstNLoops (const ACDCHitCollection *ptrHits) const
	Getter for number of loops till the first hit the collection of hits which has the Monte Carlo track id matched to this collection of hits.
Index	getLastNLoops (const ACDCHitCollection *ptrHits) const
	Getter for number of loops till the last hit the collection of hits which has the Monte Carlo track id matched to this collection of hits.
EForwardBackward	isForwardOrBackwardToMCTrack (const ACDCHitCollection *ptrHits) const
	Returns the orientation of the collection of hits relative to its matched track. More...
EForwardBackward	areAlignedInMCTrack (const ACDCHitCollection *ptrFromHits, const ACDCHitCollection *ptrToHits) const
	Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track. More...
EForwardBackward	areAlignedInMCTrackWithRLCheck (const ACDCHitCollection *ptrFromHits, const ACDCHitCollection *ptrToHits) const
	Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track. More...
CDCTrajectory3D	getTrajectory3D (const ACDCHitCollection *ptrHits) const
	Returns the trajectory of the collection of hits.

Private Member Functions
std::map< ITrackType, size_t >	getHitCountByMCTrackId (const ACDCHitCollection &hits) const
	Fill a map with the number of hits for each track id contained in the given hit range.

Private Attributes
float	m_minimalMatchPurity = 0.5
	Threshold for the purity that must be exceeded to be considered as a match.
float	m_minimalRLPurity = 0.5
	Threshold for the correct fraction of right left passage informations to be considered a match.

Detailed Description

template<typename ACDCHitCollection>
class Belle2::TrackFindingCDC::CDCMCHitCollectionLookUp< ACDCHitCollection >

Interface class to the Monte Carlo information for collections of hits.

Definition at line 70 of file CDCMCHitCollectionLookUp.h.

Member Function Documentation

areAlignedInMCTrack()

EForwardBackward areAlignedInMCTrack	(	const ACDCHitCollection *	ptrFromHits,
		const ACDCHitCollection *	ptrToHits
	)		const

Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track.

• Returns EForwardBackward::c_Invalid if the collection of hits is not matched in the same Monte Carlo track.
• Returns EForwardBackward::c_Forward if both collections of hits are coaligned with the Monte Carlo track and the second collection of hits comes after the first.
• Returns EForwardBackward::c_Backward if both collections of hits are anticoaligned with the Monte Carlo track and the second collection of hits comes after the first in the reversed sense.

Check if the segments are in the same track

Definition at line 202 of file CDCMCHitCollectionLookUp.icc.h.

    {
      ITrackType fromMCTrackId = getMCTrackId(ptrFromHits);
      if (fromMCTrackId == INVALID_ITRACK) return EForwardBackward::c_Invalid;
 
      ITrackType toMCTrackId = getMCTrackId(ptrToHits);
      if (toMCTrackId == INVALID_ITRACK) return EForwardBackward::c_Invalid;
 
      if (fromMCTrackId != toMCTrackId) return EForwardBackward::c_Invalid;
 
      // Check if the segments are sensable on their own
      EForwardBackward fromFBInfo = isForwardOrBackwardToMCTrack(ptrFromHits);
      if (fromFBInfo == EForwardBackward::c_Invalid) return EForwardBackward::c_Invalid;
 
      EForwardBackward toFBInfo = isForwardOrBackwardToMCTrack(ptrToHits);
      if (toFBInfo == EForwardBackward::c_Invalid) return EForwardBackward::c_Invalid;
 
      if (fromFBInfo != toFBInfo) return EForwardBackward::c_Invalid;
 
      {
        // Now check if hits are aligned within their common track
        // Index firstNPassedSuperLayersOfFromHits = getFirstNPassedSuperLayers(ptrFromHits);
        Index lastNPassedSuperLayersOfFromHits = getLastNPassedSuperLayers(ptrFromHits);
        // if (firstNPassedSuperLayersOfFromHits == c_InvalidIndex) return
        // EForwardBackward::c_Invalid;
        if (lastNPassedSuperLayersOfFromHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        Index firstNPassedSuperLayersOfToHits = getFirstNPassedSuperLayers(ptrToHits);
        // Index lastNPassedSuperLayersOfToHits = getLastNPassedSuperLayers(ptrToHits);
        if (firstNPassedSuperLayersOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
        // if (lastNPassedSuperLayersOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        if (lastNPassedSuperLayersOfFromHits < firstNPassedSuperLayersOfToHits) {
          if (fromFBInfo == EForwardBackward::c_Forward and
              toFBInfo == EForwardBackward::c_Forward) {
            return EForwardBackward::c_Forward;
          } else {
            return EForwardBackward::c_Invalid;
          }
        } else if (firstNPassedSuperLayersOfToHits < lastNPassedSuperLayersOfFromHits) {
          if (fromFBInfo == EForwardBackward::c_Backward and
              toFBInfo == EForwardBackward::c_Backward) {
            return EForwardBackward::c_Backward;
          } else {
            return EForwardBackward::c_Invalid;
          }
        }
      }
 
      {
        // Now we are in the same true segment with both segments
        // Index firstInTrackSegmentIdOfFromHits = getFirstInTrackSegmentId(ptrFromHits);
        Index lastInTrackSegmentIdOfFromHits = getLastInTrackSegmentId(ptrFromHits);
        // if (firstInTrackSegmentIdOfFromHits == c_InvalidIndex) return
        // EForwardBackward::c_Invalid;
        if (lastInTrackSegmentIdOfFromHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        Index firstInTrackSegmentIdOfToHits = getFirstInTrackSegmentId(ptrToHits);
        // Index lastInTrackSegmentIdOfToHits = getLastInTrackSegmentId(ptrToHits);
        if (firstInTrackSegmentIdOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
        // if (lastInTrackSegmentIdOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        if (lastInTrackSegmentIdOfFromHits < firstInTrackSegmentIdOfToHits) {
          if (fromFBInfo == EForwardBackward::c_Forward and
              toFBInfo == EForwardBackward::c_Forward) {
            return EForwardBackward::c_Forward;
          } else {
            return EForwardBackward::c_Invalid;
          }
        } else if (firstInTrackSegmentIdOfToHits < lastInTrackSegmentIdOfFromHits) {
          // Test if to segment lies before in the mc track
          // Hence the whole pair of segments is reverse to the track direction of flight
          if (fromFBInfo == EForwardBackward::c_Backward and
              toFBInfo == EForwardBackward::c_Backward) {
            return EForwardBackward::c_Backward;
          } else {
            return EForwardBackward::c_Invalid;
          }
        }
      }
 
      {
        // Now we are in the same true segment with both of the hits
        // Index firstInTrackIdOfFromHits = getFirstInTrackId(ptrFromHits);
        Index lastInTrackIdOfFromHits = getLastInTrackId(ptrFromHits);
        // if (firstInTrackIdOfFromHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
        if (lastInTrackIdOfFromHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        Index firstInTrackIdOfToHits = getFirstInTrackId(ptrToHits);
        // Index lastInTrackIdOfToHits = getLastInTrackId(ptrToHits);
        if (firstInTrackIdOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
        // if (lastInTrackIdOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        // Relax conditions somewhat such that segments may overlap at the borders.
 
        if (lastInTrackIdOfFromHits - 1 < firstInTrackIdOfToHits + 1) {
          if (fromFBInfo == EForwardBackward::c_Forward and
              toFBInfo == EForwardBackward::c_Forward) {
            return EForwardBackward::c_Forward;
          }
        }
 
        if (firstInTrackIdOfToHits - 1 < lastInTrackIdOfFromHits + 1) {
          if (fromFBInfo == EForwardBackward::c_Backward and
              toFBInfo == EForwardBackward::c_Backward) {
            return EForwardBackward::c_Backward;
          }
        }
      }
      // FIXME: Handle intertwined hits that are not cleanly consecutive along the track?
      return EForwardBackward::c_Invalid;
    }

areAlignedInMCTrackWithRLCheck()

EForwardBackward areAlignedInMCTrackWithRLCheck	(	const ACDCHitCollection *	ptrFromHits,
		const ACDCHitCollection *	ptrToHits
	)		const

Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track.

Also checks that the majority of the right left passage informations agrees with the forward backward information

Definition at line 319 of file CDCMCHitCollectionLookUp.icc.h.

getMCTrackId()

ITrackType getMCTrackId

(

const ACDCHitCollection *

ptrHits

)

const

Getter for the Monte Carlo track id matched to this collection of hits.

On first encounter of a collection of hits this evaluates the purities for the contained track ids. The match is valid if the highest purity exceeds the s_minimalMatchPurity threshold. In case the highest purity is to low to result will be INVALID_ITRACK. To result is stored in member map object for fast look up for repeated calls.

Definition at line 93 of file CDCMCHitCollectionLookUp.icc.h.

isForwardOrBackwardToMCTrack()

EForwardBackward isForwardOrBackwardToMCTrack

(

const ACDCHitCollection *

ptrHits

)

const

Returns the orientation of the collection of hits relative to its matched track.

• Returns EForwardBackward::c_Invalid if the collection of hits is not matched to any track.
• Returns EForwardBackward::c_Forward if the collection of hits is coaligned with the matched track.
• Returns EForwardBackward::c_Backward if the collection of hits is coaligned with the matched track interpreted in reverse.
• Returns EForwardBackward::c_Unknown if the matter cannot be decided.

Definition at line 184 of file CDCMCHitCollectionLookUp.icc.h.

---

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

• tracking/trackFindingCDC/mclookup/include/CDCMCHitCollectionLookUp.h
• tracking/trackFindingCDC/mclookup/include/CDCMCHitCollectionLookUp.icc.h