CDCMCTrackLookUp Class Reference

Specialisation of the lookup for the truth values of reconstructed tracks. More...

#include <CDCMCTrackLookUp.h>

Inheritance diagram for CDCMCTrackLookUp:

Public Member Functions
	CDCMCTrackLookUp ()=default
	Default constructor, needs to be public for initialization in CDCMCManager.
	CDCMCTrackLookUp (CDCMCTrackLookUp &)=delete
	Singleton: Delete copy constructor and assignment operator.
CDCMCTrackLookUp &	operator= (const CDCMCTrackLookUp &)=delete
	Operator =.
void	clear ()
	Clears all Monte Carlo information left from the last event.
MCTrackIdPurityPair	getHighestPurity (const CDCTrack &hits) const
	Get the track id with the highest corresponding purity.
ITrackType	getMCTrackId (const CDCTrack *ptrHits) const
	Getter for the Monte Carlo track id matched to this collection of hits.
int	getCorrectRLVote (const CDCTrack *ptrHits) const
	Getter for the difference of correct versus incorrect right left passage information.
double	getRLPurity (const CDCTrack *ptrHits) const
	Getter for the right left passge purity which respects the forward backward reconstruction.
const MCParticle *	getMCParticle (const CDCTrack *ptrHits) const
	Getter for the mc particle matched to this collection of hits.
const CDCHit *	getFirstHit (const CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *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 CDCTrack *ptrHits) const
	Returns the orientation of the collection of hits relative to its matched track.
EForwardBackward	areAlignedInMCTrack (const CDCTrack *ptrFromHits, const CDCTrack *ptrToHits) const
	Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track.
EForwardBackward	areAlignedInMCTrackWithRLCheck (const CDCTrack *ptrFromHits, const CDCTrack *ptrToHits) const
	Returns if the second collection of hits follows the first collection of hits in their common Monte Carlo track.
CDCTrajectory3D	getTrajectory3D (const CDCTrack *ptrHits) const
	Returns the trajectory of the collection of hits.

Static Public Member Functions
static const CDCMCTrackLookUp &	getInstance ()
	Getter for the singletone instance.

Private Member Functions
std::map< ITrackType, size_t >	getHitCountByMCTrackId (const CDCTrack &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
	Threshold for the purity that must be exceeded to be considered as a match.
float	m_minimalRLPurity
	Threshold for the correct fraction of right left passage information to be considered a match.

Detailed Description

Specialisation of the lookup for the truth values of reconstructed tracks.

Definition at line 23 of file CDCMCTrackLookUp.h.

Member Function Documentation

areAlignedInMCTrack()

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

inherited

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 168 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 CDCTrack *  ptrFromHits, const CDCTrack *  ptrToHits  ) const

inherited

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

inherited

Clears all Monte Carlo information left from the last event.

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

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

getCorrectRLVote()

int getCorrectRLVote ( const CDCTrack *  ptrHits ) const

inherited

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

Definition at line 95 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 CDCTrack *  ptrHits ) const

inherited

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 105 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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack &  hits ) const

inherited

Get the track id with the highest corresponding purity.

Definition at line 85 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);
      // cppcheck-suppress ignoredReturnValue
      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 CDCTrack &  hits ) const

privateinherited

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

Definition at line 81 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;
    }

getInstance()

const CDCMCTrackLookUp & getInstance ( )

static

Getter for the singletone instance.

Definition at line 21 of file CDCMCTrackLookUp.cc.

{
  return CDCMCManager::getMCTrackLookUp();
}

getLastHit()

const CDCHit * getLastHit ( const CDCTrack *  ptrHits ) const

inherited

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 108 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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inlineinherited

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 CDCTrack *  ptrHits ) const

inherited

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

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

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

getMCTrackId()

ITrackType getMCTrackId ( const CDCTrack *  ptrHits ) const

inherited

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 CDCTrack *  ptrHits ) const

inherited

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

Definition at line 98 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 CDCTrack *  ptrHits ) const

inherited

Returns the trajectory of the collection of hits.

Definition at line 180 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 CDCTrack *  ptrHits ) const

inherited

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

privateinherited

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

privateinherited

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/CDCMCTrackLookUp.h
• tracking/trackFindingCDC/mclookup/src/CDCMCTrackLookUp.cc