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.
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.
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.
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.
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.
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 70 of file CDCMCHitCollectionLookUp.h.

Constructor & Destructor Documentation

CDCMCHitCollectionLookUp()

CDCMCHitCollectionLookUp ( )

inline

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

Definition at line 77 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 201 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 meaningful 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 information agrees with the forward backward information

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

    {
      EForwardBackward result = areAlignedInMCTrack(ptrFromHits, ptrToHits);
      if (result == EForwardBackward::c_Invalid) return result;
 
      int fromCorrectRLVote = getCorrectRLVote(ptrFromHits);
      int toCorrectRLVote = getCorrectRLVote(ptrToHits);
 
      if (result == 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 EForwardBackward::c_Invalid;
    }

clear()

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

int getCorrectRLVote

(

const ACDCHitCollection *

ptrHits

)

const

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

Definition at line 105 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()

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

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 113 of file CDCMCHitCollectionLookUp.h.

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

getFirstInTrackSegmentId()

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 123 of file CDCMCHitCollectionLookUp.h.

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

getFirstNLoops()

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 141 of file CDCMCHitCollectionLookUp.h.

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

getFirstNPassedSuperLayers()

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 133 of file CDCMCHitCollectionLookUp.h.

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

getHighestPurity()

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(), LessOf<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;
        ERightLeft mcRLInfo = mcHitLookUp.getRLInfo(hit);
        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()

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

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 165 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 : reversedRange(hits)) {
        if (mcTrackId == mcHitLookUp.getMCTrackId(hit)) return hit;
      }
      return nullptr;
    }

getLastInTrackId()

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 117 of file CDCMCHitCollectionLookUp.h.

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

getLastInTrackSegmentId()

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 127 of file CDCMCHitCollectionLookUp.h.

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

getLastNLoops()

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 145 of file CDCMCHitCollectionLookUp.h.

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

getLastNPassedSuperLayers()

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 137 of file CDCMCHitCollectionLookUp.h.

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

getMCParticle()

const MCParticle * getMCParticle

(

const ACDCHitCollection *

ptrHits

)

const

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

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

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

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

double getRLPurity

(

const ACDCHitCollection *

ptrHits

)

const

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

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

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

getTrajectory3D()

CDCTrajectory3D getTrajectory3D

(

const ACDCHitCollection *

ptrHits

)

const

Returns the trajectory of the collection of hits.

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

    {
      if (not ptrHits) {
        B2WARNING("Segment is nullptr. Could not get fit.");
        return 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 CDCTrajectory3D();
        }
      }
 
      const CDCSimHit& primarySimHit = *ptrPrimarySimHit;
 
      Vector3D mom3D{primarySimHit.getMomentum()};
      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 CDCTrajectory3D();
      }
 
      const TParticlePDG& tPDGParticle = *ptrTPDGParticle;
 
      double charge = tPDGParticle.Charge() / 3.0;
 
      ESign chargeSign = sign(charge);
 
      CDCTrajectory3D trajectory3D(pos3D, time, mom3D, charge);
 
      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()

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 183 of file CDCMCHitCollectionLookUp.icc.h.

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

Member Data Documentation

m_minimalMatchPurity

float m_minimalMatchPurity = 0.5

private

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

Definition at line 184 of file CDCMCHitCollectionLookUp.h.

m_minimalRLPurity

float m_minimalRLPurity = 0.5

private

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

Definition at line 187 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