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.

int	getCorrectRLVote (const ACDCHitCollection *ptrHits) const
	Getter for the difference of correct versus incorrect right left passage information.

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.

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

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

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

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

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

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

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

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

TrackingUtilities::EForwardBackward	isForwardOrBackwardToMCTrack (const ACDCHitCollection *ptrHits) const
	Returns the orientation of the collection of hits relative to its matched track.

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

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

TrackingUtilities::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 information 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 73 of file CDCMCHitCollectionLookUp.h.

Constructor & Destructor Documentation

◆ CDCMCHitCollectionLookUp()

template<typename ACDCHitCollection>

CDCMCHitCollectionLookUp ( )

inline

Dummy constructor to silence the intel-compiler warnings (NOTE: using "=default" will make the intel compilation fail)

Definition at line 80 of file CDCMCHitCollectionLookUp.h.

80{};

Member Function Documentation

◆ areAlignedInMCTrack()

template<class ACDCHitCollection>

TrackingUtilities::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 TrackingUtilities::EForwardBackward::c_Invalid;
 
      ITrackType toMCTrackId = getMCTrackId(ptrToHits);
      if (toMCTrackId == INVALID_ITRACK) return TrackingUtilities::EForwardBackward::c_Invalid;
 
      if (fromMCTrackId != toMCTrackId) return TrackingUtilities::EForwardBackward::c_Invalid;
 
      // Check if the segments are meaningful on their own
      TrackingUtilities::EForwardBackward fromFBInfo = isForwardOrBackwardToMCTrack(ptrFromHits);
      if (fromFBInfo == TrackingUtilities::EForwardBackward::c_Invalid) return TrackingUtilities::EForwardBackward::c_Invalid;
 
      TrackingUtilities::EForwardBackward toFBInfo = isForwardOrBackwardToMCTrack(ptrToHits);
      if (toFBInfo == TrackingUtilities::EForwardBackward::c_Invalid) return TrackingUtilities::EForwardBackward::c_Invalid;
 
      if (fromFBInfo != toFBInfo) return TrackingUtilities::EForwardBackward::c_Invalid;
 
      {
        // Now check if hits are aligned within their common track
        // Index firstNPassedSuperLayersOfFromHits = getFirstNPassedSuperLayers(ptrFromHits);
        TrackingUtilities::Index lastNPassedSuperLayersOfFromHits = getLastNPassedSuperLayers(ptrFromHits);
        // if (firstNPassedSuperLayersOfFromHits == c_InvalidIndex) return
        // EForwardBackward::c_Invalid;
        if (lastNPassedSuperLayersOfFromHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::EForwardBackward::c_Invalid;
 
        TrackingUtilities::Index firstNPassedSuperLayersOfToHits = getFirstNPassedSuperLayers(ptrToHits);
        // Index lastNPassedSuperLayersOfToHits = getLastNPassedSuperLayers(ptrToHits);
        if (firstNPassedSuperLayersOfToHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::EForwardBackward::c_Invalid;
        // if (lastNPassedSuperLayersOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        if (lastNPassedSuperLayersOfFromHits < firstNPassedSuperLayersOfToHits) {
          if (fromFBInfo == TrackingUtilities::EForwardBackward::c_Forward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Forward) {
            return TrackingUtilities::EForwardBackward::c_Forward;
          } else {
            return TrackingUtilities::EForwardBackward::c_Invalid;
          }
        } else if (firstNPassedSuperLayersOfToHits < lastNPassedSuperLayersOfFromHits) {
          if (fromFBInfo == TrackingUtilities::EForwardBackward::c_Backward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Backward) {
            return TrackingUtilities::EForwardBackward::c_Backward;
          } else {
            return TrackingUtilities::EForwardBackward::c_Invalid;
          }
        }
      }
 
      {
        // Now we are in the same true segment with both segments
        // Index firstInTrackSegmentIdOfFromHits = getFirstInTrackSegmentId(ptrFromHits);
        TrackingUtilities::Index lastInTrackSegmentIdOfFromHits = getLastInTrackSegmentId(ptrFromHits);
        // if (firstInTrackSegmentIdOfFromHits == c_InvalidIndex) return
        // EForwardBackward::c_Invalid;
        if (lastInTrackSegmentIdOfFromHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::EForwardBackward::c_Invalid;
 
        TrackingUtilities::Index firstInTrackSegmentIdOfToHits = getFirstInTrackSegmentId(ptrToHits);
        // Index lastInTrackSegmentIdOfToHits = getLastInTrackSegmentId(ptrToHits);
        if (firstInTrackSegmentIdOfToHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::EForwardBackward::c_Invalid;
        // if (lastInTrackSegmentIdOfToHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
 
        if (lastInTrackSegmentIdOfFromHits < firstInTrackSegmentIdOfToHits) {
          if (fromFBInfo == TrackingUtilities::EForwardBackward::c_Forward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Forward) {
            return TrackingUtilities::EForwardBackward::c_Forward;
          } else {
            return TrackingUtilities::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 == TrackingUtilities::EForwardBackward::c_Backward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Backward) {
            return TrackingUtilities::EForwardBackward::c_Backward;
          } else {
            return TrackingUtilities::EForwardBackward::c_Invalid;
          }
        }
      }
 
      {
        // Now we are in the same true segment with both of the hits
        // Index firstInTrackIdOfFromHits = getFirstInTrackId(ptrFromHits);
        TrackingUtilities::Index lastInTrackIdOfFromHits = getLastInTrackId(ptrFromHits);
        // if (firstInTrackIdOfFromHits == c_InvalidIndex) return EForwardBackward::c_Invalid;
        if (lastInTrackIdOfFromHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::EForwardBackward::c_Invalid;
 
        TrackingUtilities::Index firstInTrackIdOfToHits = getFirstInTrackId(ptrToHits);
        // Index lastInTrackIdOfToHits = getLastInTrackId(ptrToHits);
        if (firstInTrackIdOfToHits == TrackingUtilities::c_InvalidIndex) return TrackingUtilities::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 == TrackingUtilities::EForwardBackward::c_Forward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Forward) {
            return TrackingUtilities::EForwardBackward::c_Forward;
          }
        }
 
        if (firstInTrackIdOfToHits - 1 < lastInTrackIdOfFromHits + 1) {
          if (fromFBInfo == TrackingUtilities::EForwardBackward::c_Backward and
              toFBInfo == TrackingUtilities::EForwardBackward::c_Backward) {
            return TrackingUtilities::EForwardBackward::c_Backward;
          }
        }
      }
      // FIXME: Handle intertwined hits that are not cleanly consecutive along the track?
      return TrackingUtilities::EForwardBackward::c_Invalid;
    }

◆ areAlignedInMCTrackWithRLCheck()

template<class ACDCHitCollection>

TrackingUtilities::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 information agrees with the forward backward information

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

    {
      TrackingUtilities::EForwardBackward result = areAlignedInMCTrack(ptrFromHits, ptrToHits);
      if (result == TrackingUtilities::EForwardBackward::c_Invalid) return result;
 
      int fromCorrectRLVote = getCorrectRLVote(ptrFromHits);
      int toCorrectRLVote = getCorrectRLVote(ptrToHits);
 
      if (result == TrackingUtilities::EForwardBackward::c_Backward) {
        fromCorrectRLVote = -fromCorrectRLVote;
        toCorrectRLVote = -toCorrectRLVote;
      }
 
      int fromNCorrectRL = (fromCorrectRLVote + ptrFromHits->size()) / 2;
      int toNCorrectRL = (toCorrectRLVote + ptrToHits->size()) / 2;
 
      float fromRLPurity = 1.0 * fromNCorrectRL / ptrFromHits->size();
      float toRLPurity = 1.0 * toNCorrectRL / ptrToHits->size();
 
      // Require the minimal rl purity and also at least 2.5 correct hits
      // (cut chosen to require all correct in single hit triplet segment)
      if (fromRLPurity > m_minimalRLPurity and toRLPurity > m_minimalRLPurity and
          fromNCorrectRL > 2.5 and toNCorrectRL > 2.5) {
        return result;
      }
 
      return TrackingUtilities::EForwardBackward::c_Invalid;
    }

◆ clear()

template<class ACDCHitCollection>

void clear ( )

Clears all Monte Carlo information left from the last event.

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

    {
      B2DEBUG(25, "Clearing CDCMCHitCollectionLookUp<ACDCHitCollection>");
    }

◆ getCorrectRLVote()

template<class ACDCHitCollection>

int getCorrectRLVote ( const ACDCHitCollection * ptrHits ) const

Getter for the difference of correct versus incorrect right left passage information.

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

    {
      if (not ptrHits) return INVALID_ITRACK;
      const ACDCHitCollection& hits = *ptrHits;
      MCTrackIdPurityPair mcTrackIdAndPurity = getHighestPurity(hits);
      if (mcTrackIdAndPurity.getPurity() >= m_minimalMatchPurity) {
        return mcTrackIdAndPurity.getCorrectRLVote();
      } else {
        return 0;
      }
    }

◆ getFirstHit()

template<class ACDCHitCollection>

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.

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

    {
      if (not ptrHits) return nullptr;
      const ACDCHitCollection& hits = *ptrHits;
 
      ITrackType mcTrackId = getMCTrackId(ptrHits);
      if (mcTrackId == INVALID_ITRACK) return nullptr;
 
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
 
      for (const CDCHit* hit : hits) {
        if (mcTrackId == mcHitLookUp.getMCTrackId(hit)) return hit;
      }
      return nullptr;
    }

◆ getFirstInTrackId()

template<typename ACDCHitCollection>

TrackingUtilities::Index getFirstInTrackId ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 116 of file CDCMCHitCollectionLookUp.h.

117 { return CDCMCHitLookUp::getInstance().getInTrackId(getFirstHit(ptrHits)); }

◆ getFirstInTrackSegmentId()

template<typename ACDCHitCollection>

TrackingUtilities::Index getFirstInTrackSegmentId ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 126 of file CDCMCHitCollectionLookUp.h.

127 { return CDCMCHitLookUp::getInstance().getInTrackSegmentId(getFirstHit(ptrHits)); }

◆ getFirstNLoops()

template<typename ACDCHitCollection>

TrackingUtilities::Index getFirstNLoops ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 144 of file CDCMCHitCollectionLookUp.h.

145 { return CDCMCHitLookUp::getInstance().getNLoops(getFirstHit(ptrHits)); }

◆ getFirstNPassedSuperLayers()

template<typename ACDCHitCollection>

TrackingUtilities::Index getFirstNPassedSuperLayers ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 136 of file CDCMCHitCollectionLookUp.h.

137 { return CDCMCHitLookUp::getInstance().getNPassedSuperLayers(getFirstHit(ptrHits)); }

◆ getHighestPurity()

template<class ACDCHitCollection>

MCTrackIdPurityPair getHighestPurity ( const ACDCHitCollection & hits ) const

Get the track id with the highest corresponding purity.

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

    {
      std::map<ITrackType, size_t> hitCountByMCTrackId = getHitCountByMCTrackId(hits);
 
      size_t nHits = 0;
      std::pair<ITrackType, size_t> highestHitCountMCTrackId(0, 0);
      static_cast<void>(std::max_element(hitCountByMCTrackId.begin(), hitCountByMCTrackId.end(),
                                         TrackingUtilities::LessOf<TrackingUtilities::Second>()));
 
      for (const auto& hitCountForMCTrackId : hitCountByMCTrackId) {
 
        nHits += hitCountForMCTrackId.second;
 
        if (highestHitCountMCTrackId.second < hitCountForMCTrackId.second) {
          highestHitCountMCTrackId = hitCountForMCTrackId;
        }
      }
 
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
 
      int correctRLVote = 0;
      for (const auto& recoHit : hits) {
        const CDCHit* hit = recoHit;
        TrackingUtilities::ERightLeft mcRLInfo = mcHitLookUp.getRLInfo(hit);
        TrackingUtilities::ERightLeft rlInfo = recoHit.getRLInfo();
        if (rlInfo == mcRLInfo) {
          ++correctRLVote;
        } else {
          --correctRLVote;
        }
      }
 
      const float purity = static_cast<float>(highestHitCountMCTrackId.second) / nHits;
      return MCTrackIdPurityPair(highestHitCountMCTrackId.first, purity, correctRLVote);
    }

◆ getHitCountByMCTrackId()

template<class ACDCHitCollection>

std::map< ITrackType, size_t > getHitCountByMCTrackId ( const ACDCHitCollection & hits ) const

private

Fill a map with the number of hits for each track id contained in the given hit range.

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

    {
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
 
      std::map<ITrackType, size_t> hitCountByMCTrackId;
      for (const CDCHit* ptrHit : hits) {
        ITrackType mcTrackId = mcHitLookUp.getMCTrackId(ptrHit);
        // cppcheck-suppress stlFindInsert
        if (hitCountByMCTrackId.count(mcTrackId) == 0) hitCountByMCTrackId[mcTrackId] = 0;
        ++(hitCountByMCTrackId[mcTrackId]);
      }
      return hitCountByMCTrackId;
    }

◆ getLastHit()

template<class ACDCHitCollection>

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.

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

    {
 
      if (not ptrHits) return nullptr;
      const ACDCHitCollection& hits = *ptrHits;
 
      ITrackType mcTrackId = getMCTrackId(ptrHits);
      if (mcTrackId == INVALID_ITRACK) return nullptr;
 
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
 
      for (const CDCHit* hit : TrackingUtilities::reversedRange(hits)) {
        if (mcTrackId == mcHitLookUp.getMCTrackId(hit)) return hit;
      }
      return nullptr;
    }

◆ getLastInTrackId()

template<typename ACDCHitCollection>

TrackingUtilities::Index getLastInTrackId ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 120 of file CDCMCHitCollectionLookUp.h.

121 { return CDCMCHitLookUp::getInstance().getInTrackId(getLastHit(ptrHits)); }

◆ getLastInTrackSegmentId()

template<typename ACDCHitCollection>

TrackingUtilities::Index getLastInTrackSegmentId ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 130 of file CDCMCHitCollectionLookUp.h.

131 { return CDCMCHitLookUp::getInstance().getInTrackSegmentId(getLastHit(ptrHits)); }

◆ getLastNLoops()

template<typename ACDCHitCollection>

TrackingUtilities::Index getLastNLoops ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 148 of file CDCMCHitCollectionLookUp.h.

149 { return CDCMCHitLookUp::getInstance().getNLoops(getLastHit(ptrHits)); }

◆ getLastNPassedSuperLayers()

template<typename ACDCHitCollection>

TrackingUtilities::Index getLastNPassedSuperLayers ( const ACDCHitCollection * ptrHits ) const

inline

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.

Definition at line 140 of file CDCMCHitCollectionLookUp.h.

141 { return CDCMCHitLookUp::getInstance().getNPassedSuperLayers(getLastHit(ptrHits)); }

◆ getMCParticle()

template<class ACDCHitCollection>

const MCParticle * getMCParticle ( const ACDCHitCollection * ptrHits ) const

Getter for the mc particle matched to this collection of hits.

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

    {
      const CDCHit* ptrHit = getFirstHit(ptrHits);
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
      return mcHitLookUp.getMCParticle(ptrHit);
    }

◆ getMCTrackId()

template<class ACDCHitCollection>

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.

    {
      if (not ptrHits) return INVALID_ITRACK;
      const ACDCHitCollection& hits = *ptrHits;
      MCTrackIdPurityPair mcTrackIdAndPurity = getHighestPurity(hits);
      if (mcTrackIdAndPurity.getPurity() >= m_minimalMatchPurity) {
        return mcTrackIdAndPurity.getMCTrackId();
      } else {
        return INVALID_ITRACK;
      }
    }

◆ getRLPurity()

template<class ACDCHitCollection>

double getRLPurity ( const ACDCHitCollection * ptrHits ) const

Getter for the right left passge purity which respects the forward backward reconstruction.

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

    {
      TrackingUtilities::EForwardBackward fbInfo = isForwardOrBackwardToMCTrack(ptrHits);
      if (fbInfo == TrackingUtilities::EForwardBackward::c_Invalid) return NAN;
 
      int correctRLVote = getCorrectRLVote(ptrHits);
 
      if (fbInfo == TrackingUtilities::EForwardBackward::c_Backward) {
        correctRLVote = -correctRLVote;
      }
 
      int nCorrectRL = (correctRLVote + ptrHits->size()) / 2;
      float rlPurity = 1.0 * nCorrectRL / ptrHits->size();
      return rlPurity;
    }

◆ getTrajectory3D()

template<class ACDCHitCollection>

TrackingUtilities::CDCTrajectory3D getTrajectory3D ( const ACDCHitCollection * ptrHits ) const

Returns the trajectory of the collection of hits.

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

    {
      if (not ptrHits) {
        B2WARNING("Segment is nullptr. Could not get fit.");
        return TrackingUtilities::CDCTrajectory3D();
      }
 
      const CDCMCHitLookUp& mcHitLookUp = CDCMCHitLookUp::getInstance();
 
      const CDCHit* ptrFirstHit = getFirstHit(ptrHits);
      const CDCSimHit* ptrPrimarySimHit = mcHitLookUp.getClosestPrimarySimHit(ptrFirstHit);
 
      if (not ptrPrimarySimHit) {
        // If there is no primary SimHit simply use the secondary simhit as reference
        ptrPrimarySimHit = mcHitLookUp.getSimHit(ptrFirstHit);
        if (not ptrPrimarySimHit) {
          return TrackingUtilities::CDCTrajectory3D();
        }
      }
 
      const CDCSimHit& primarySimHit = *ptrPrimarySimHit;
 
      TrackingUtilities::Vector3D mom3D{primarySimHit.getMomentum()};
      TrackingUtilities::Vector3D pos3D{primarySimHit.getPosTrack()};
      double time{primarySimHit.getFlightTime()};
 
      int pdgCode = primarySimHit.getPDGCode();
      const TParticlePDG* ptrTPDGParticle = TDatabasePDG::Instance()->GetParticle(pdgCode);
 
      if (not ptrTPDGParticle) {
        B2WARNING("No particle for PDG code " << pdgCode << ". Could not get fit");
        return TrackingUtilities::CDCTrajectory3D();
      }
 
      const TParticlePDG& tPDGParticle = *ptrTPDGParticle;
 
      double charge = tPDGParticle.Charge() / 3.0;
 
      TrackingUtilities::ESign chargeSign = TrackingUtilities::sign(charge);
 
      TrackingUtilities::CDCTrajectory3D trajectory3D(pos3D, time, mom3D, charge);
 
      TrackingUtilities::ESign settedChargeSign = trajectory3D.getChargeSign();
 
      if (chargeSign != settedChargeSign) {
        B2WARNING("Charge sign of mc particle is not the same as the one of the fit");
      }
 
      return trajectory3D;
    }

◆ isForwardOrBackwardToMCTrack()

template<class ACDCHitCollection>

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

    {
      TrackingUtilities::Index firstInTrackId = getFirstInTrackId(ptrHits);
      TrackingUtilities::Index lastInTrackId = getLastInTrackId(ptrHits);
      if (firstInTrackId == TrackingUtilities::c_InvalidIndex or lastInTrackId == TrackingUtilities::c_InvalidIndex) {
        return TrackingUtilities::EForwardBackward::c_Invalid;
      } else if (firstInTrackId < lastInTrackId) {
        return TrackingUtilities::EForwardBackward::c_Forward;
      } else if (firstInTrackId > lastInTrackId) {
        return TrackingUtilities::EForwardBackward::c_Backward;
      } else if (firstInTrackId == lastInTrackId) {
        return TrackingUtilities::EForwardBackward::c_Unknown;
      }
      return TrackingUtilities::EForwardBackward::c_Invalid;
    }

Member Data Documentation

◆ m_minimalMatchPurity

template<typename ACDCHitCollection>

float m_minimalMatchPurity = 0.5

private

Threshold for the purity that must be exceeded to be considered as a match.

Definition at line 187 of file CDCMCHitCollectionLookUp.h.

◆ m_minimalRLPurity

template<typename ACDCHitCollection>

float m_minimalRLPurity = 0.5

private

Threshold for the correct fraction of right left passage information to be considered a match.

Definition at line 190 of file CDCMCHitCollectionLookUp.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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CDCMCHitCollectionLookUp()

Member Function Documentation

◆ areAlignedInMCTrack()

◆ areAlignedInMCTrackWithRLCheck()

◆ clear()

◆ getCorrectRLVote()

◆ getFirstHit()

◆ getFirstInTrackId()

◆ getFirstInTrackSegmentId()

◆ getFirstNLoops()

◆ getFirstNPassedSuperLayers()

◆ getHighestPurity()

◆ getHitCountByMCTrackId()

◆ getLastHit()

◆ getLastInTrackId()

◆ getLastInTrackSegmentId()

◆ getLastNLoops()

◆ getLastNPassedSuperLayers()

◆ getMCParticle()

◆ getMCTrackId()

◆ getRLPurity()

◆ getTrajectory3D()

◆ isForwardOrBackwardToMCTrack()

Member Data Documentation

◆ m_minimalMatchPurity

◆ m_minimalRLPurity