mkfit::Track Class Reference

#include <Track.h>

Inheritance diagram for mkfit::Track:

Public Member Functions
void	addHitIdx (int hitIdx, int hitLyr, float chi2)
void	addHitIdx (const HitOnTrack &hot, float chi2)
const HitOnTrack *	beginHitsOnTrack () const
HitOnTrack *	beginHitsOnTrack_nc ()
bool	canReachRadius (float R) const
void	countAndSetNFoundHits ()
const HitOnTrack *	endHitsOnTrack () const
std::vector< int >	foundLayers () const
int	getHitIdx (int posHitIdx) const
int	getHitLyr (int posHitIdx) const
HitOnTrack	getHitOnTrack (int posHitIdx) const
const HitOnTrack *	getHitsOnTrackArray () const
int	getLastFoundHitIdx () const
int	getLastFoundHitLyr () const
HitOnTrack	getLastFoundHitOnTrack () const
int	getLastFoundHitPos () const
int	getLastFoundMCHitID (const std::vector< HitVec > &globalHitVec) const
int	getLastHitIdx () const
int	getLastHitLyr () const
HitOnTrack	getLastHitOnTrack () const
int	getMCHitIDFromLayer (const std::vector< HitVec > &globalHitVec, int layer) const
HitVec	hitsVector (const std::vector< HitVec > &globalHitVec) const
float	maxReachRadius () const
void	mcHitIDsVec (const std::vector< HitVec > &globalHitVec, const MCHitInfoVec &globalMCHitInfo, std::vector< int > &mcHitIDs) const
int	nFoundHits () const
int	nInsideMinusOneHits () const
int	nOverlapHits () const
int	nTailMinusOneHits () const
int	nTotalHits () const
int	nUniqueLayers () const
float	rAtZ (float Z) const
void	reserveHits (int nHits)
void	resetHits ()
void	resizeHits (int nHits, int nFoundHits)
void	resizeHitsForInput ()
void	setHitIdx (int posHitIdx, int newIdx)
void	setHitIdxAtPos (int pos, const HitOnTrack &hot)
void	setHitIdxLyr (int posHitIdx, int newIdx, int newLyr)
void	setNOverlapHits (int n)
void	sortHitsByLayer ()
float	swimPhiToR (const float x, const float y) const
	Track ()
	Track (const TrackBase &base)
	Track (const TrackState &state, float chi2, int label, int nHits, const HitOnTrack *hits)
	Track (int charge, const SVector3 &position, const SVector3 &momentum, const SMatrixSym66 &errors, float chi2)
	Track (const Track &t)
float	zAtR (float R, float *r_reached=nullptr) const
Public Member Functions inherited from mkfit::TrackBase
int	algoint () const
TrackAlgorithm	algorithm () const
int	charge () const
float	chi2 () const
float	d0BeamSpot (const float x_bs, const float y_bs, bool linearize=false) const
float	emomEta () const
float	emomPhi () const
float	epT () const
const float *	errArray () const
const SMatrixSym66 &	errors () const
SMatrixSym66 &	errors_nc ()
bool	getDuplicateValue () const
int	getEtaRegion () const
int	getNSeedHits () const
Status	getStatus () const
bool	hasNanNSillyValues () const
bool	hasSillyValues (bool dump, bool fix, const char *pref="")
float	invpT () const
bool	isFindable () const
bool	isNotFindable () const
int	label () const
SVector3	momentum () const
float	momEta () const
float	momPhi () const
float	p () const
const SVector6 &	parameters () const
SVector6 &	parameters_nc ()
const float *	posArray () const
float	posEta () const
SVector3	position () const
float	posPhi () const
float	posR () const
float	posRsq () const
ProdType	prodType () const
float	pT () const
float	px () const
float	py () const
float	pz () const
float	score () const
void	setAlgoint (int algo)
void	setAlgorithm (TrackAlgorithm algo)
void	setCharge (int chg)
void	setChi2 (float chi2)
void	setDuplicateValue (bool d)
void	setEtaRegion (int r)
void	setLabel (int lbl)
void	setNotFindable ()
void	setNSeedHits (int n)
void	setProdType (ProdType ptyp)
void	setScore (float s)
void	setState (const TrackState &newState)
void	setStatus (Status s)
const TrackState &	state () const
TrackState &	state_nc ()
float	theta () const
	TrackBase ()
	TrackBase (const TrackState &state, float chi2, int label)
	TrackBase (int charge, const SVector3 &position, const SVector3 &momentum, const SMatrixSym66 &errors, float chi2)
float	x () const
float	y () const
float	z () const

Private Attributes
std::vector< HitOnTrack >	hitsOnTrk_

Additional Inherited Members
Public Types inherited from mkfit::TrackBase
enum	ProdType { ProdType::NotSet = 0, ProdType::Signal = 1, ProdType::InTimePU = 2, ProdType::OutOfTimePU = 3 }
enum	TrackAlgorithm {   TrackAlgorithm::undefAlgorithm = 0, TrackAlgorithm::ctf = 1, TrackAlgorithm::duplicateMerge = 2, TrackAlgorithm::cosmics = 3,   TrackAlgorithm::initialStep = 4, TrackAlgorithm::lowPtTripletStep = 5, TrackAlgorithm::pixelPairStep = 6, TrackAlgorithm::detachedTripletStep = 7,   TrackAlgorithm::mixedTripletStep = 8, TrackAlgorithm::pixelLessStep = 9, TrackAlgorithm::tobTecStep = 10, TrackAlgorithm::jetCoreRegionalStep = 11,   TrackAlgorithm::conversionStep = 12, TrackAlgorithm::muonSeededStepInOut = 13, TrackAlgorithm::muonSeededStepOutIn = 14, TrackAlgorithm::outInEcalSeededConv = 15,   TrackAlgorithm::inOutEcalSeededConv = 16, TrackAlgorithm::nuclInter = 17, TrackAlgorithm::standAloneMuon = 18, TrackAlgorithm::globalMuon = 19,   TrackAlgorithm::cosmicStandAloneMuon = 20, TrackAlgorithm::cosmicGlobalMuon = 21, TrackAlgorithm::highPtTripletStep = 22, TrackAlgorithm::lowPtQuadStep = 23,   TrackAlgorithm::detachedQuadStep = 24, TrackAlgorithm::reservedForUpgrades1 = 25, TrackAlgorithm::reservedForUpgrades2 = 26, TrackAlgorithm::bTagGhostTracks = 27,   TrackAlgorithm::beamhalo = 28, TrackAlgorithm::gsf = 29, TrackAlgorithm::hltPixel = 30, TrackAlgorithm::hltIter0 = 31,   TrackAlgorithm::hltIter1 = 32, TrackAlgorithm::hltIter2 = 33, TrackAlgorithm::hltIter3 = 34, TrackAlgorithm::hltIter4 = 35,   TrackAlgorithm::hltIterX = 36, TrackAlgorithm::hiRegitMuInitialStep = 37, TrackAlgorithm::hiRegitMuLowPtTripletStep = 38, TrackAlgorithm::hiRegitMuPixelPairStep = 39,   TrackAlgorithm::hiRegitMuDetachedTripletStep = 40, TrackAlgorithm::hiRegitMuMixedTripletStep = 41, TrackAlgorithm::hiRegitMuPixelLessStep = 42, TrackAlgorithm::hiRegitMuTobTecStep = 43,   TrackAlgorithm::hiRegitMuMuonSeededStepInOut = 44, TrackAlgorithm::hiRegitMuMuonSeededStepOutIn = 45, TrackAlgorithm::algoSize = 46 }
	track algorithm; copy from TrackBase.h to keep in standalone builds More...
Static Public Member Functions inherited from mkfit::TrackBase
static const char *	algoint_to_cstr (int algo)
Protected Attributes inherited from mkfit::TrackBase
float	chi2_ = 0.
int	label_ = -1
short int	lastHitIdx_ = -1
short int	nFoundHits_ = 0
float	score_ = 0.
TrackState	state_
Status	status_

Detailed Description

Definition at line 365 of file Track.h.

Constructor & Destructor Documentation

Track() [1/5]

mkfit::Track::Track ( )

inline

Definition at line 367 of file Track.h.

{}

Track() [2/5]

mkfit::Track::Track ( const TrackBase &  base )

inlineexplicit

Definition at line 369 of file Track.h.

References mkfit::TrackBase::lastHitIdx_, and mkfit::TrackBase::nFoundHits_.

                                          : TrackBase(base) {
      // Reset hit counters -- caller has to initialize hits.
      lastHitIdx_ = -1;
      nFoundHits_ = 0;
    }

Track() [3/5]

mkfit::Track::Track ( const TrackState &  state, float  chi2, int  label, int  nHits, const HitOnTrack *  hits  )

inline

Definition at line 375 of file Track.h.

References addHitIdx(), f, h, hfClusterShapes_cfi::hits, phase1PixelTopology::layer, nHits, and reserveHits().

        : TrackBase(state, chi2, label) {
      reserveHits(nHits);
      for (int h = 0; h < nHits; ++h) {
        addHitIdx(hits[h].index, hits[h].layer, 0.0f);
      }
    }

Track() [4/5]

mkfit::Track::Track ( int  charge, const SVector3 &  position, const SVector3 &  momentum, const SMatrixSym66 &  errors, float  chi2  )

inline

Definition at line 383 of file Track.h.

        : TrackBase(charge, position, momentum, errors, chi2) {}

Track() [5/5]

mkfit::Track::Track ( const Track &  t )

inline

Definition at line 386 of file Track.h.

: TrackBase(t), hitsOnTrk_(t.hitsOnTrk_) {}

Member Function Documentation

addHitIdx() [1/2]

void mkfit::Track::addHitIdx ( int  hitIdx, int  hitLyr, float  chi2  )

inline

Definition at line 442 of file Track.h.

References mkfit::TrackBase::chi2(), mkfit::TrackBase::chi2_, hitsOnTrk_, mkfit::TrackBase::lastHitIdx_, and mkfit::TrackBase::nFoundHits_.

Referenced by mkfit::StdSeq::clean_cms_seedtracks_iter(), and Track().

                                                       {
      hitsOnTrk_.push_back({hitIdx, hitLyr});
      ++lastHitIdx_;
      if (hitIdx >= 0 || hitIdx == -9) {
        ++nFoundHits_;
        chi2_ += chi2;
      }
    }

addHitIdx() [2/2]

void mkfit::Track::addHitIdx ( const HitOnTrack &  hot, float  chi2  )

inline

Definition at line 451 of file Track.h.

References addHitIdx(), mkfit::TrackBase::chi2(), mkfit::HitOnTrack::index, and mkfit::HitOnTrack::layer.

Referenced by addHitIdx().

{ addHitIdx(hot.index, hot.layer, chi2); }

beginHitsOnTrack()

const HitOnTrack* mkfit::Track::beginHitsOnTrack ( ) const

inline

Definition at line 499 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::MkFinder::copy_in().

{ return hitsOnTrk_.data(); }

beginHitsOnTrack_nc()

HitOnTrack* mkfit::Track::beginHitsOnTrack_nc ( )

inline

Definition at line 502 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::MkFinder::copy_out().

{ return hitsOnTrk_.data(); }

canReachRadius()

bool mkfit::Track::canReachRadius

(

float

R

)

const

Definition at line 282 of file Track.cc.

References funct::abs(), mkfit::Config::Bfield, mkfit::TrackBase::charge(), f, dqmdumpme::k, mkfit::TrackBase::pT(), dttmaxenums::R, mkfit::Const::sol, mkfit::TrackBase::x(), and mkfit::TrackBase::y().

                                          {
    const float k = ((charge() < 0) ? 100.0f : -100.0f) / (Const::sol * Config::Bfield);
    const float ooc = 2.0f * k * pT();
    return std::abs(ooc) > R - std::hypot(x(), y());
  }

countAndSetNFoundHits()

void mkfit::Track::countAndSetNFoundHits ( )

inline

Definition at line 508 of file Track.h.

References hitsOnTrk_, mps_fire::i, mkfit::TrackBase::lastHitIdx_, and mkfit::TrackBase::nFoundHits_.

                                 {
      nFoundHits_ = 0;
      for (int i = 0; i <= lastHitIdx_; i++) {
        if (hitsOnTrk_[i].index >= 0 || hitsOnTrk_[i].index == -9)
          nFoundHits_++;
      }
    }

endHitsOnTrack()

const HitOnTrack* mkfit::Track::endHitsOnTrack ( ) const

inline

Definition at line 500 of file Track.h.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

Referenced by mkfit::MkFinder::copy_in().

{ return hitsOnTrk_.data() + (lastHitIdx_ + 1); }

foundLayers()

std::vector<int> mkfit::Track::foundLayers ( ) const

inline

Definition at line 573 of file Track.h.

References hitsOnTrk_, mkfit::TrackBase::lastHitIdx_, phase1PixelTopology::layer, and hgcalTopologyTester_cfi::layers.

                                       {
      std::vector<int> layers;
      for (int ihit = 0; ihit <= lastHitIdx_; ++ihit) {
        if (hitsOnTrk_[ihit].index >= 0 || hitsOnTrk_[ihit].index == -9) {
          layers.push_back(hitsOnTrk_[ihit].layer);
        }
      }
      return layers;
    }

getHitIdx()

int mkfit::Track::getHitIdx ( int  posHitIdx ) const

inline

Definition at line 455 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::TrackExtra::findMatchingSeedHits(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHits(), mkfit::print(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), and mkfit::TrackExtra::setMCTrackIDInfo().

{ return hitsOnTrk_[posHitIdx].index; }

getHitLyr()

int mkfit::Track::getHitLyr ( int  posHitIdx ) const

inline

Definition at line 456 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::TrackExtra::findMatchingSeedHits(), mkfit::MkFitter::inputTracksAndHits(), mkfit::print(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), and mkfit::TrackExtra::setMCTrackIDInfo().

{ return hitsOnTrk_[posHitIdx].layer; }

getHitOnTrack()

HitOnTrack mkfit::Track::getHitOnTrack ( int  posHitIdx ) const

inline

Definition at line 453 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), and mkfit::MkFitter::inputTracksAndHits().

{ return hitsOnTrk_[posHitIdx]; }

getHitsOnTrackArray()

const HitOnTrack* mkfit::Track::getHitsOnTrackArray ( ) const

inline

Definition at line 498 of file Track.h.

References hitsOnTrk_.

Referenced by mkfit::MkFinder::bkFitInputTracks(), and mkfit::MkFitter::inputTracksForFit().

{ return hitsOnTrk_.data(); }

getLastFoundHitIdx()

int mkfit::Track::getLastFoundHitIdx ( ) const

inline

Definition at line 473 of file Track.h.

References getLastFoundHitPos(), hitsOnTrk_, and mkfit::TrackBase::p().

                                   {
      int p = getLastFoundHitPos();
      return p >= 0 ? hitsOnTrk_[p].index : -1;
    }

getLastFoundHitLyr()

int mkfit::Track::getLastFoundHitLyr ( ) const

inline

Definition at line 477 of file Track.h.

References getLastFoundHitPos(), hitsOnTrk_, and mkfit::TrackBase::p().

                                   {
      int p = getLastFoundHitPos();
      return p >= 0 ? hitsOnTrk_[p].layer : -1;
    }

getLastFoundHitOnTrack()

HitOnTrack mkfit::Track::getLastFoundHitOnTrack ( ) const

inline

Definition at line 469 of file Track.h.

References getLastFoundHitPos(), hitsOnTrk_, and mkfit::TrackBase::p().

Referenced by getLastFoundMCHitID().

                                              {
      int p = getLastFoundHitPos();
      return p >= 0 ? hitsOnTrk_[p] : HitOnTrack(-1, -1);
    }

getLastFoundHitPos()

int mkfit::Track::getLastFoundHitPos ( ) const

inline

Definition at line 462 of file Track.h.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

Referenced by getLastFoundHitIdx(), getLastFoundHitLyr(), and getLastFoundHitOnTrack().

                                   {
      int hi = lastHitIdx_;
      while (hi >= 0 && hitsOnTrk_[hi].index < 0)
        --hi;
      return hi;
    }

getLastFoundMCHitID()

int mkfit::Track::getLastFoundMCHitID ( const std::vector< HitVec > &  globalHitVec ) const

inline

Definition at line 482 of file Track.h.

References getLastFoundHitOnTrack(), mkfit::HitOnTrack::index, and mkfit::HitOnTrack::layer.

                                                                         {
      HitOnTrack hot = getLastFoundHitOnTrack();
      return globalHitVec[hot.layer][hot.index].mcHitID();
    }

getLastHitIdx()

int mkfit::Track::getLastHitIdx ( ) const

inline

Definition at line 459 of file Track.h.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

{ return hitsOnTrk_[lastHitIdx_].index; }

getLastHitLyr()

int mkfit::Track::getLastHitLyr ( ) const

inline

Definition at line 460 of file Track.h.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

{ return hitsOnTrk_[lastHitIdx_].layer; }

getLastHitOnTrack()

HitOnTrack mkfit::Track::getLastHitOnTrack ( ) const

inline

Definition at line 458 of file Track.h.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

{ return hitsOnTrk_[lastHitIdx_]; }

getMCHitIDFromLayer()

int mkfit::Track::getMCHitIDFromLayer ( const std::vector< HitVec > &  globalHitVec, int  layer  ) const

inline

Definition at line 487 of file Track.h.

References hitsOnTrk_, mkfit::TrackBase::lastHitIdx_, and phase1PixelTopology::layer.

                                                                                    {
      int mcHitID = -1;
      for (int ihit = 0; ihit <= lastHitIdx_; ++ihit) {
        if (hitsOnTrk_[ihit].layer == layer) {
          mcHitID = globalHitVec[hitsOnTrk_[ihit].layer][hitsOnTrk_[ihit].index].mcHitID();
          break;
        }
      }
      return mcHitID;
    }

hitsVector()

HitVec mkfit::Track::hitsVector ( const std::vector< HitVec > &  globalHitVec ) const

inline

Definition at line 397 of file Track.h.

References hitsOnTrk_, mkfit::HitOnTrack::index, mkfit::HitOnTrack::layer, and mkfit::Config::nMaxTrkHits.

                                                                   {
      HitVec hitsVec;
      for (int ihit = 0; ihit < Config::nMaxTrkHits; ++ihit) {
        const HitOnTrack& hot = hitsOnTrk_[ihit];
        if (hot.index >= 0) {
          hitsVec.push_back(globalHitVec[hot.layer][hot.index]);
        }
      }
      return hitsVec;
    }

maxReachRadius()

float mkfit::Track::maxReachRadius

(

)

const

Definition at line 288 of file Track.cc.

References funct::abs(), mkfit::Config::Bfield, mkfit::TrackBase::charge(), f, dqmdumpme::k, mkfit::TrackBase::pT(), mkfit::TrackBase::px(), mkfit::TrackBase::py(), mkfit::Const::sol, mkfit::TrackBase::x(), and mkfit::TrackBase::y().

                                    {
    const float k = ((charge() < 0) ? 100.0f : -100.0f) / (Const::sol * Config::Bfield);
    const float abs_ooc_half = std::abs(k * pT());
    // center of helix in x,y plane
    const float x_center = x() - k * py();
    const float y_center = y() + k * px();
    return std::hypot(x_center, y_center) + abs_ooc_half;
  }

mcHitIDsVec()

void mkfit::Track::mcHitIDsVec ( const std::vector< HitVec > &  globalHitVec, const MCHitInfoVec &  globalMCHitInfo, std::vector< int > &  mcHitIDs  ) const

inline

Definition at line 408 of file Track.h.

References hitsOnTrk_, mkfit::HitOnTrack::index, mkfit::TrackBase::lastHitIdx_, and mkfit::HitOnTrack::layer.

                                                     {
      for (int ihit = 0; ihit <= lastHitIdx_; ++ihit) {
        const HitOnTrack& hot = hitsOnTrk_[ihit];
        if ((hot.index >= 0) && (static_cast<size_t>(hot.index) < globalHitVec[hot.layer].size())) {
          mcHitIDs.push_back(globalHitVec[hot.layer][hot.index].mcTrackID(globalMCHitInfo));
        } else {
          mcHitIDs.push_back(hot.index);
        }
      }
    }

nFoundHits()

int mkfit::Track::nFoundHits ( ) const

inline

Definition at line 516 of file Track.h.

References mkfit::TrackBase::nFoundHits_.

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::MkFinder::copy_in(), mkfit::getScoreCand(), mkfit::print(), mkfit::StdSeq::Quality::quality_process(), resizeHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), mkfit::TrackExtra::setMCTrackIDInfo(), and mkfit::sortByHitsChi2().

{ return nFoundHits_; }

nInsideMinusOneHits()

int mkfit::Track::nInsideMinusOneHits ( ) const

inline

Definition at line 522 of file Track.h.

References hitsOnTrk_, mps_fire::i, mkfit::TrackBase::lastHitIdx_, and dqmiodumpmetadata::n.

Referenced by mkfit::MkFinder::copy_in(), and mkfit::getScoreCand().

                                    {
      int n = 0;
      bool insideValid = false;
      for (int i = lastHitIdx_; i >= 0; --i) {
        if (hitsOnTrk_[i].index >= 0)
          insideValid = true;
        if (insideValid && hitsOnTrk_[i].index == -1)
          ++n;
      }
      return n;
    }

nOverlapHits()

int mkfit::Track::nOverlapHits ( ) const

inline

Definition at line 519 of file Track.h.

References mkfit::TrackBase::Status::n_overlaps, and mkfit::TrackBase::status_.

Referenced by mkfit::getScoreCand().

{ return status_.n_overlaps; }

nTailMinusOneHits()

int mkfit::Track::nTailMinusOneHits ( ) const

inline

Definition at line 534 of file Track.h.

References hitsOnTrk_, mps_fire::i, mkfit::TrackBase::lastHitIdx_, and dqmiodumpmetadata::n.

Referenced by mkfit::MkFinder::copy_in(), and mkfit::getScoreCand().

                                  {
      int n = 0;
      for (int i = lastHitIdx_; i >= 0; --i) {
        if (hitsOnTrk_[i].index >= 0)
          return n;
        if (hitsOnTrk_[i].index == -1)
          ++n;
      }
      return n;
    }

nTotalHits()

int mkfit::Track::nTotalHits ( ) const

inline

Definition at line 517 of file Track.h.

References mkfit::TrackBase::lastHitIdx_.

Referenced by mkfit::MkFinder::bkFitInputTracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::MkFinder::copy_in(), mkfit::TrackExtra::findMatchingSeedHits(), mkfit::print(), mkfit::MkFinder::selectHitIndices(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), and mkfit::TrackExtra::setMCTrackIDInfo().

{ return lastHitIdx_ + 1; }

nUniqueLayers()

int mkfit::Track::nUniqueLayers ( ) const

inline

Definition at line 545 of file Track.h.

References hitsOnTrk_, heavyIonCSV_trainingSettings::idx, mkfit::TrackBase::lastHitIdx_, and jetUpdater_cfi::sort.

                              {
      // make local copy in vector: sort it in place
      std::vector<HitOnTrack> tmp_hitsOnTrk(hitsOnTrk_.begin(), hitsOnTrk_.end());
      std::sort(tmp_hitsOnTrk.begin(), tmp_hitsOnTrk.end(), [](const auto& h1, const auto& h2) {
        return h1.layer < h2.layer;
      });

      // local counters
      auto lyr_cnt = 0;
      auto prev_lyr = -1;

      // loop over copy of hitsOnTrk
      for (auto ihit = 0; ihit <= lastHitIdx_; ++ihit) {
        const auto& hot = tmp_hitsOnTrk[ihit];
        const auto lyr = hot.layer;
        const auto idx = hot.index;
        if (lyr >= 0 && (idx >= 0 || idx == -9) && lyr != prev_lyr) {
          ++lyr_cnt;
          prev_lyr = lyr;
        }
      }
      return lyr_cnt;
    }

rAtZ()

float mkfit::Track::rAtZ

(

float

Z

)

const

Definition at line 374 of file Track.cc.

References alpha, mkfit::Config::Bfield, mkfit::TrackBase::charge(), funct::cos(), PVValHelper::dz, f, mkfit::TrackBase::invpT(), dqmdumpme::k, mkfit::TrackBase::px(), mkfit::TrackBase::py(), funct::sin(), mkfit::Const::sol, funct::tan(), mkfit::TrackBase::theta(), mkfit::TrackBase::x(), mkfit::TrackBase::y(), BeamSpotPI::Z, and mkfit::TrackBase::z().

                                 {
    float xc = x();
    float yc = y();
    float pxc = px();
    float pyc = py();

    const float ipt = invpT();
    const float kinv = ((charge() < 0) ? 0.01f : -0.01f) * Const::sol * Config::Bfield;
    const float k = 1.0f / kinv;

    const float dz = Z - z();
    const float alpha = dz * ipt * kinv * std::tan(theta());

    const float cosa = std::cos(alpha);
    const float sina = std::sin(alpha);

    xc += k * (pxc * sina - pyc * (1.0f - cosa));
    yc += k * (pyc * sina + pxc * (1.0f - cosa));

    // const float pxo = pxc;//copy before overwriting
    // pxc = pxc * cosa  -  pyc * sina;
    // pyc = pyc * cosa  +  pxo * sina;

    return std::hypot(xc, yc);
  }

reserveHits()

void mkfit::Track::reserveHits ( int  nHits )

inline

Definition at line 423 of file Track.h.

References hitsOnTrk_, and nHits.

Referenced by Track().

{ hitsOnTrk_.reserve(nHits); }

resetHits()

void mkfit::Track::resetHits ( )

inline

Definition at line 425 of file Track.h.

References hitsOnTrk_, mkfit::TrackBase::lastHitIdx_, and mkfit::TrackBase::nFoundHits_.

                     {
      lastHitIdx_ = -1;
      nFoundHits_ = 0;
      hitsOnTrk_.clear();
    }

resizeHits()

void mkfit::Track::resizeHits ( int  nHits, int  nFoundHits  )

inline

Definition at line 432 of file Track.h.

References hitsOnTrk_, mkfit::TrackBase::lastHitIdx_, nFoundHits(), mkfit::TrackBase::nFoundHits_, and nHits.

Referenced by mkfit::MkFinder::copy_out().

                                               {
      hitsOnTrk_.resize(nHits);
      lastHitIdx_ = nHits - 1;
      nFoundHits_ = nFoundHits;
    }

resizeHitsForInput()

void mkfit::Track::resizeHitsForInput

(

)

Definition at line 261 of file Track.cc.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

                                 {
    bzero(&hitsOnTrk_, sizeof(hitsOnTrk_));
    hitsOnTrk_.resize(lastHitIdx_ + 1);
  }

setHitIdx()

void mkfit::Track::setHitIdx ( int  posHitIdx, int  newIdx  )

inline

Definition at line 504 of file Track.h.

References hitsOnTrk_.

{ hitsOnTrk_[posHitIdx].index = newIdx; }

setHitIdxAtPos()

void mkfit::Track::setHitIdxAtPos ( int  pos, const HitOnTrack &  hot  )

inline

Definition at line 438 of file Track.h.

References hitsOnTrk_.

{ hitsOnTrk_[pos] = hot; }

setHitIdxLyr()

void mkfit::Track::setHitIdxLyr ( int  posHitIdx, int  newIdx, int  newLyr  )

inline

Definition at line 506 of file Track.h.

References hitsOnTrk_.

{ hitsOnTrk_[posHitIdx] = {newIdx, newLyr}; }

setNOverlapHits()

void mkfit::Track::setNOverlapHits ( int  n )

inline

Definition at line 520 of file Track.h.

References dqmiodumpmetadata::n, mkfit::TrackBase::Status::n_overlaps, and mkfit::TrackBase::status_.

{ status_.n_overlaps = n; }

sortHitsByLayer()

void mkfit::Track::sortHitsByLayer

(

)

Definition at line 266 of file Track.cc.

References hitsOnTrk_, and mkfit::TrackBase::lastHitIdx_.

Referenced by mkfit::StdSeq::clean_cms_seedtracks_iter().

                              {
    std::stable_sort(&hitsOnTrk_[0], &hitsOnTrk_[lastHitIdx_ + 1], [](const auto& h1, const auto& h2) {
      return h1.layer < h2.layer;
    });
  }

swimPhiToR()

float mkfit::Track::swimPhiToR	(	const float	x,
		const float	y
	)		const

Definition at line 272 of file Track.cc.

References funct::abs(), mkfit::TrackBase::charge(), HLT_2022v12_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, mkfit::getHypot(), mkfit::TrackBase::momPhi(), mkfit::Const::PI, mkfit::TrackBase::pT(), mkfit::squashPhiGeneral(), findQualityFiles::v, mkfit::TrackBase::x(), and mkfit::TrackBase::y().

                                                              {
    const float dR = getHypot(x() - x0, y() - y0);
    // XXX-ASSUMPTION-ERROR can not always reach R, should see what callers expect.
    // For now return PI to signal apex on the ohter side of the helix.
    const float v = dR / 176.f / pT() * charge();
    const float dPhi = std::abs(v) <= 1.0f ? 2.f * std::asin(v) : Const::PI;
    ;
    return squashPhiGeneral(momPhi() - dPhi);
  }

zAtR()

float mkfit::Track::zAtR	(	float	R,
		float *	r_reached = nullptr
	)		const

Definition at line 297 of file Track.cc.

References mkfit::Config::Bfield, c, mkfit::TrackBase::charge(), funct::cos(), f, mps_fire::i, triggerObjects_cff::id, mkfit::TrackBase::invpT(), dqmdumpme::k, mkfit::Config::Niter, mkfit::TrackBase::px(), mkfit::TrackBase::py(), mkfit::TrackBase::pz(), dttmaxenums::R, funct::sin(), mkfit::Const::sol, mkfit::TrackBase::x(), mkfit::TrackBase::y(), and mkfit::TrackBase::z().

                                                   {
    float xc = x();
    float yc = y();
    float pxc = px();
    float pyc = py();

    const float ipt = invpT();
    const float kinv = ((charge() < 0) ? 0.01f : -0.01f) * Const::sol * Config::Bfield;
    const float k = 1.0f / kinv;

    const float c = 0.5f * kinv * ipt;
    const float ooc = 1.0f / c;  // 2 * radius of curvature
    const float lambda = pz() * ipt;

    //printf("Track::zAtR to R=%f: k=%e, ipt=%e, c=%e, ooc=%e  -- can hit = %f (if > 1 can)\n",
    //       R, k, ipt, c, ooc, ooc / (R - std::hypot(xc,yc)));

    float D = 0;

    for (int i = 0; i < Config::Niter; ++i) {
      // compute tangental and ideal distance for the current iteration.
      // 3-rd order asin for symmetric incidence (shortest arc lenght).
      float r0 = std::hypot(xc, yc);
      float td = (R - r0) * c;
      float id = ooc * td * (1.0f + 0.16666666f * td * td);
      // This would be for line approximation:
      // float id = R - r0;
      D += id;

      //printf("%-3d r0=%f R-r0=%f td=%f id=%f id_line=%f delta_id=%g\n",
      //       i, r0, R-r0, td, id, R - r0, id - (R-r0));

      float cosa = std::cos(id * ipt * kinv);
      float sina = std::sin(id * ipt * kinv);

      //update parameters
      xc += k * (pxc * sina - pyc * (1.0f - cosa));
      yc += k * (pyc * sina + pxc * (1.0f - cosa));

      const float pxo = pxc;  //copy before overwriting
      pxc = pxc * cosa - pyc * sina;
      pyc = pyc * cosa + pxo * sina;
    }

    if (r_reached)
      *r_reached = std::hypot(xc, yc);

    return z() + lambda * D;

    // ----------------------------------------------------------------
    // Exact solution from Avery's notes ... loses precision somewhere
    // {
    //   const float a = kinv;
    //   float pT      = S.pT();

    //   float ax2y2  = a*(x*x + y*y);
    //   float T      = std::sqrt(pT*pT - 2.0f*a*(x*py - y*px) + a*ax2y2);
    //   float D0     = (T - pT) / a;
    //   float D      = (-2.0f * (x*py - y*px) + a * (x*x + y*y)) / (T + pT);

    //   float B      = c * std::sqrt((R*R - D*D) / (1.0f + 2.0f*c*D));
    //   float s1     = std::asin(B) / c;
    //   float s2     = (Const::PI - std::asin(B)) / c;

    //   printf("pt %f, invpT %f\n", pT, S.invpT());
    //   printf("lambda %f, a %f, c %f, T %f, D0 %f, D %f, B %f, s1 %f, s2 %f\n",
    //          lambda, a, c, T, D0, D, B, s1, s2);
    //   printf("%f = %f / %f\n", (R*R - D*D) / (1.0f + 2.0f*c*D), (R*R - D*D), (1.0f + 2.0f*c*D));

    //   z1 = S.z() + lambda * s1;
    //   z2 = S.z() + lambda * s2;

    //   printf("z1=%f z2=%f\n", z1, z2);
    // }
    // ----------------------------------------------------------------
  }

Member Data Documentation

hitsOnTrk_

std::vector<HitOnTrack> mkfit::Track::hitsOnTrk_

private

Definition at line 584 of file Track.h.

Referenced by addHitIdx(), beginHitsOnTrack(), beginHitsOnTrack_nc(), countAndSetNFoundHits(), endHitsOnTrack(), foundLayers(), getHitIdx(), getHitLyr(), getHitOnTrack(), getHitsOnTrackArray(), getLastFoundHitIdx(), getLastFoundHitLyr(), getLastFoundHitOnTrack(), getLastFoundHitPos(), getLastHitIdx(), getLastHitLyr(), getLastHitOnTrack(), getMCHitIDFromLayer(), hitsVector(), mcHitIDsVec(), nInsideMinusOneHits(), nTailMinusOneHits(), nUniqueLayers(), reserveHits(), resetHits(), resizeHits(), resizeHitsForInput(), setHitIdx(), setHitIdxAtPos(), setHitIdxLyr(), and sortHitsByLayer().