mkfit::TrackBase Class Reference

#include <Track.h>

Inheritance diagram for mkfit::TrackBase:

Classes
struct	Status

Public Types
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...

Public Member Functions
int	algoint () const
TrackAlgorithm	algorithm () const
bool	canReachRadius (float R) 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
float	maxReachRadius () 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	rAtZ (float Z) 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	swimPhiToR (const float x, const float y) const
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
float	zAtR (float R, float *r_reached=nullptr) const

Static Public Member Functions
static const char *	algoint_to_cstr (int algo)

Protected Attributes
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 120 of file Track.h.

Member Enumeration Documentation

ProdType

enum mkfit::TrackBase::ProdType

strong

Enumerator
NotSet
Signal
InTimePU
OutOfTimePU

Definition at line 257 of file Track.h.

{ NotSet = 0, Signal = 1, InTimePU = 2, OutOfTimePU = 3 };

TrackAlgorithm

enum mkfit::TrackBase::TrackAlgorithm

strong

track algorithm; copy from TrackBase.h to keep in standalone builds

Enumerator
undefAlgorithm
ctf
duplicateMerge
cosmics
initialStep
lowPtTripletStep
pixelPairStep
detachedTripletStep
mixedTripletStep
pixelLessStep
tobTecStep
jetCoreRegionalStep
conversionStep
muonSeededStepInOut
muonSeededStepOutIn
outInEcalSeededConv
inOutEcalSeededConv
nuclInter
standAloneMuon
globalMuon
cosmicStandAloneMuon
cosmicGlobalMuon
highPtTripletStep
lowPtQuadStep
detachedQuadStep
reservedForUpgrades1
reservedForUpgrades2
bTagGhostTracks
beamhalo
gsf
hltPixel
hltIter0
hltIter1
hltIter2
hltIter3
hltIter4
hltIterX
hiRegitMuInitialStep
hiRegitMuLowPtTripletStep
hiRegitMuPixelPairStep
hiRegitMuDetachedTripletStep
hiRegitMuMixedTripletStep
hiRegitMuPixelLessStep
hiRegitMuTobTecStep
hiRegitMuMuonSeededStepInOut
hiRegitMuMuonSeededStepOutIn
algoSize

Definition at line 272 of file Track.h.

                              {
      undefAlgorithm = 0,
      ctf = 1,
      duplicateMerge = 2,
      cosmics = 3,
      initialStep = 4,
      lowPtTripletStep = 5,
      pixelPairStep = 6,
      detachedTripletStep = 7,
      mixedTripletStep = 8,
      pixelLessStep = 9,
      tobTecStep = 10,
      jetCoreRegionalStep = 11,
      conversionStep = 12,
      muonSeededStepInOut = 13,
      muonSeededStepOutIn = 14,
      outInEcalSeededConv = 15,
      inOutEcalSeededConv = 16,
      nuclInter = 17,
      standAloneMuon = 18,
      globalMuon = 19,
      cosmicStandAloneMuon = 20,
      cosmicGlobalMuon = 21,
      // Phase1
      highPtTripletStep = 22,
      lowPtQuadStep = 23,
      detachedQuadStep = 24,
      reservedForUpgrades1 = 25,
      reservedForUpgrades2 = 26,
      bTagGhostTracks = 27,
      beamhalo = 28,
      gsf = 29,
      // HLT algo name
      hltPixel = 30,
      // steps used by PF
      hltIter0 = 31,
      hltIter1 = 32,
      hltIter2 = 33,
      hltIter3 = 34,
      hltIter4 = 35,
      // steps used by all other objects @HLT
      hltIterX = 36,
      // steps used by HI muon regional iterative tracking
      hiRegitMuInitialStep = 37,
      hiRegitMuLowPtTripletStep = 38,
      hiRegitMuPixelPairStep = 39,
      hiRegitMuDetachedTripletStep = 40,
      hiRegitMuMixedTripletStep = 41,
      hiRegitMuPixelLessStep = 42,
      hiRegitMuTobTecStep = 43,
      hiRegitMuMuonSeededStepInOut = 44,
      hiRegitMuMuonSeededStepOutIn = 45,
      algoSize = 46
    };

Constructor & Destructor Documentation

TrackBase() [1/3]

mkfit::TrackBase::TrackBase ( )

inline

Definition at line 122 of file Track.h.

{}

TrackBase() [2/3]

mkfit::TrackBase::TrackBase ( const TrackState &  state, float  chi2, int  label  )

inline

Definition at line 124 of file Track.h.

: state_(state), chi2_(chi2), label_(label) {}

TrackBase() [3/3]

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

inline

Definition at line 126 of file Track.h.

        : state_(charge, position, momentum, errors), chi2_(chi2) {}

Member Function Documentation

algoint()

int mkfit::TrackBase::algoint ( ) const

inline

Definition at line 327 of file Track.h.

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

{ return status_.algorithm; }

algoint_to_cstr()

const char * mkfit::TrackBase::algoint_to_cstr ( int  algo )

static

Definition at line 331 of file Track.cc.

References algoSize, and names.

Referenced by mkfit::Shell::Compare(), mkfit::Shell::FillByLabelMaps_CkfBase(), mkfit::Shell::PrintIterations(), mkfit::ConfigJson::save_Iterations(), and mkfit::Shell::Status().

                                                 {
    static const char* const names[] = {"undefAlgorithm",
                                        "ctf",
                                        "duplicateMerge",
                                        "cosmics",
                                        "initialStep",
                                        "lowPtTripletStep",
                                        "pixelPairStep",
                                        "detachedTripletStep",
                                        "mixedTripletStep",
                                        "pixelLessStep",
                                        "tobTecStep",
                                        "jetCoreRegionalStep",
                                        "conversionStep",
                                        "muonSeededStepInOut",
                                        "muonSeededStepOutIn",
                                        "outInEcalSeededConv",
                                        "inOutEcalSeededConv",
                                        "nuclInter",
                                        "standAloneMuon",
                                        "globalMuon",
                                        "cosmicStandAloneMuon",
                                        "cosmicGlobalMuon",
                                        "highPtTripletStep",
                                        "lowPtQuadStep",
                                        "detachedQuadStep",
                                        "reservedForUpgrades1",
                                        "reservedForUpgrades2",
                                        "bTagGhostTracks",
                                        "beamhalo",
                                        "gsf",
                                        "hltPixel",
                                        "hltIter0",
                                        "hltIter1",
                                        "hltIter2",
                                        "hltIter3",
                                        "hltIter4",
                                        "hltIterX",
                                        "hiRegitMuInitialStep",
                                        "hiRegitMuLowPtTripletStep",
                                        "hiRegitMuPixelPairStep",
                                        "hiRegitMuDetachedTripletStep",
                                        "hiRegitMuMixedTripletStep",
                                        "hiRegitMuPixelLessStep",
                                        "hiRegitMuTobTecStep",
                                        "hiRegitMuMuonSeededStepInOut",
                                        "hiRegitMuMuonSeededStepOutIn",
                                        "algoSize"};

    if (algo < 0 || algo >= (int)TrackAlgorithm::algoSize)
      return names[0];
    return names[algo];
  }

algorithm()

TrackAlgorithm mkfit::TrackBase::algorithm ( ) const

inline

Definition at line 328 of file Track.h.

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

Referenced by mkfit::Event::select_tracks_iter().

{ return TrackAlgorithm(status_.algorithm); }

canReachRadius()

bool mkfit::TrackBase::canReachRadius

(

float

R

)

const

Definition at line 213 of file Track.cc.

References funct::abs(), mkfit::Config::Bfield, charge(), f, Matriplex::hypot(), dqmdumpme::k, pT(), dttmaxenums::R, mkfit::Const::sol, x(), and 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());
  }

charge()

int mkfit::TrackBase::charge ( void  ) const

inline

Definition at line 171 of file Track.h.

References mkfit::TrackState::charge, and state_.

Referenced by mkfit::MkFinder::bkFitInputTracks(), canReachRadius(), mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::MkFinder::copy_in(), d0BeamSpot(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), mkfit::MkFitter::inputTracksAndHits(), mkfit::MkFitter::inputTracksForFit(), maxReachRadius(), rAtZ(), mkfit::MkFinder::selectHitIndices(), mkfit::MkFitter::slurpInTracksAndHits(), swimPhiToR(), and zAtR().

{ return state_.charge; }

chi2()

float mkfit::TrackBase::chi2 ( void  ) const

inline

Definition at line 172 of file Track.h.

References chi2_.

Referenced by mkfit::TrackCand::addHitIdx(), mkfit::Track::addHitIdx(), mkfit::MkFinder::bkFitOutputTracks(), mkfit::TrackCand::considerHitForOverlap(), mkfit::MkFinder::copy_in(), mkfit::getScoreCand(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), mkfit::MkFitter::inputTracksAndHits(), mkfit::MkFitter::inputTracksForFit(), mkfit::StdSeq::Quality::quality_process(), setChi2(), mkfit::MkFitter::slurpInTracksAndHits(), and mkfit::sortByHitsChi2().

{ return chi2_; }

d0BeamSpot()

float mkfit::TrackBase::d0BeamSpot	(	const float	x_bs,
		const float	y_bs,
		bool	linearize = false
	)		const

Definition at line 191 of file Track.cc.

References funct::abs(), mkfit::Config::Bfield, charge(), funct::cos(), f, Matriplex::hypot(), dqmdumpme::k, momPhi(), pT(), px(), py(), funct::sin(), mkfit::Const::sol, x(), and y().

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

                                                                                      {
    if (linearize) {
      return std::abs(std::cos(momPhi()) * (y() - y_bs) - std::sin(momPhi()) * (x() - x_bs));
    } else {
      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 - x_bs, y_center - y_bs) - abs_ooc_half;
    }
  }

emomEta()

float mkfit::TrackBase::emomEta ( ) const

inline

Definition at line 167 of file Track.h.

References mkfit::TrackState::emomEta(), and state_.

{ return state_.emomEta(); }

emomPhi()

float mkfit::TrackBase::emomPhi ( ) const

inline

Definition at line 166 of file Track.h.

References mkfit::TrackState::emomPhi(), and state_.

{ return state_.emomPhi(); }

epT()

float mkfit::TrackBase::epT ( ) const

inline

Definition at line 165 of file Track.h.

References mkfit::TrackState::epT(), and state_.

{ return state_.epT(); }

errArray()

const float* mkfit::TrackBase::errArray ( ) const

inline

Definition at line 136 of file Track.h.

References mkfit::TrackState::errors, and state_.

{ return state_.errors.Array(); }

errors()

const SMatrixSym66& mkfit::TrackBase::errors ( void  ) const

inline

Definition at line 133 of file Track.h.

References mkfit::TrackState::errors, and state_.

Referenced by mkfit::MkFinder::copy_in(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), mkfit::MkFitter::inputTracksAndHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), and mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams().

{ return state_.errors; }

errors_nc()

SMatrixSym66& mkfit::TrackBase::errors_nc ( )

inline

Definition at line 140 of file Track.h.

References mkfit::TrackState::errors, and state_.

Referenced by mkfit::MkFinder::bkFitOutputTracks(), and mkfit::MkFinder::copy_out().

{ return state_.errors; }

getDuplicateValue()

bool mkfit::TrackBase::getDuplicateValue ( ) const

inline

Definition at line 256 of file Track.h.

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

Referenced by mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), and mkfit::TrackExtra::setMCTrackIDInfo().

{ return status_.duplicate; }

getEtaRegion()

int mkfit::TrackBase::getEtaRegion ( ) const

inline

Definition at line 263 of file Track.h.

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

Referenced by mkfit::Shell::CheckMkFitLayerPlanVsReferenceHits().

{ return status_.eta_region; }

getNSeedHits()

int mkfit::TrackBase::getNSeedHits ( ) const

inline

Definition at line 261 of file Track.h.

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

Referenced by mkfit::CombCandidate::compactifyHitStorageForBestCand(), and mkfit::TrackCand::exportTrack().

{ return status_.n_seed_hits; }

getStatus()

Status mkfit::TrackBase::getStatus ( ) const

inline

Definition at line 248 of file Track.h.

References status_.

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

{ return status_; }

hasNanNSillyValues()

bool mkfit::TrackBase::hasNanNSillyValues

(

)

const

Definition at line 177 of file Track.cc.

References mkfit::TrackState::errors, mps_fire::i, mkfit::isFinite(), dqmiolumiharvest::j, mkfit::LL, and state_.

                                           {
    bool is_silly = false;
    for (int i = 0; i < LL; ++i) {
      for (int j = 0; j <= i; ++j) {
        if ((i == j && state_.errors.At(i, j) < 0) || !isFinite(state_.errors.At(i, j))) {
          is_silly = true;
          return is_silly;
        }
      }
    }
    return is_silly;
  }

hasSillyValues()

bool mkfit::TrackBase::hasSillyValues	(	bool	dump,
		bool	fix,
		const char *	pref = ""
	)

Definition at line 155 of file Track.cc.

References GCP_Ntuples_cfg::dump, mkfit::TrackState::errors, mps_fire::i, mkfit::isFinite(), dqmiolumiharvest::j, label(), mkfit::LL, pT(), and state_.

                                                                      {
    bool is_silly = false;
    for (int i = 0; i < LL; ++i) {
      for (int j = 0; j <= i; ++j) {
        if ((i == j && state_.errors.At(i, j) < 0) || !isFinite(state_.errors.At(i, j))) {
          if (!is_silly) {
            is_silly = true;
            if (dump)
              printf("%s (label=%d, pT=%f):", pref, label(), pT());
          }
          if (dump)
            printf(" (%d,%d)=%e", i, j, state_.errors.At(i, j));
          if (fix)
            state_.errors.At(i, j) = 0.00001;
        }
      }
    }
    if (is_silly && dump)
      printf("\n");
    return is_silly;
  }

invpT()

float mkfit::TrackBase::invpT ( ) const

inline

Definition at line 158 of file Track.h.

References mkfit::TrackState::invpT(), and state_.

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), rAtZ(), and zAtR().

{ return state_.invpT(); }

isFindable()

bool mkfit::TrackBase::isFindable ( ) const

inline

Definition at line 251 of file Track.h.

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

Referenced by mkfit::MkFinder::bkFitFitTracks(), and mkfit::MkFinder::selectHitIndices().

{ return !status_.not_findable; }

isNotFindable()

bool mkfit::TrackBase::isNotFindable ( ) const

inline

Definition at line 252 of file Track.h.

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

{ return status_.not_findable; }

label()

int mkfit::TrackBase::label ( ) const

inline

Definition at line 174 of file Track.h.

References label_.

Referenced by mkfit::MkFinder::bkFitFitTracks(), mkfit::Shell::CheckMkFitLayerPlanVsReferenceHits(), mkfit::MkFinder::copy_in(), mkfit::TrackCand::exportTrack(), hasSillyValues(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), mkfit::MkFitter::inputTracksAndHits(), mkfit::MkFitter::inputTracksForFit(), mkfit::print(), mkfit::StdSeq::Quality::quality_process(), mkfit::MkFinder::selectHitIndices(), mkfit::TrackExtra::setMCTrackIDInfo(), and mkfit::MkFitter::slurpInTracksAndHits().

{ return label_; }

maxReachRadius()

float mkfit::TrackBase::maxReachRadius

(

)

const

Definition at line 219 of file Track.cc.

References funct::abs(), mkfit::Config::Bfield, charge(), f, Matriplex::hypot(), dqmdumpme::k, pT(), px(), py(), mkfit::Const::sol, x(), and 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;
  }

momentum()

SVector3 mkfit::TrackBase::momentum ( ) const

inline

Definition at line 144 of file Track.h.

References mkfit::TrackState::parameters, and state_.

{ return SVector3(state_.parameters[3], state_.parameters[4], state_.parameters[5]); }

momEta()

float mkfit::TrackBase::momEta ( ) const

inline

Definition at line 161 of file Track.h.

References mkfit::TrackState::momEta(), and state_.

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), mkfit::StdSeq::Quality::quality_process(), mkfit::MkFinder::selectHitIndices(), and mkfit::sortTracksByEta().

{ return state_.momEta(); }

momPhi()

float mkfit::TrackBase::momPhi ( ) const

inline

Definition at line 160 of file Track.h.

References mkfit::TrackState::momPhi(), and state_.

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), d0BeamSpot(), mkfit::StdSeq::Quality::quality_process(), mkfit::MkFinder::selectHitIndices(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), mkfit::TrackExtra::setMCTrackIDInfo(), mkfit::sortTracksByPhi(), and swimPhiToR().

{ return state_.momPhi(); }

p()

float mkfit::TrackBase::p ( ) const

inline

Definition at line 159 of file Track.h.

References mkfit::TrackState::p(), and state_.

Referenced by mkfit::Track::getLastFoundHitIdx(), mkfit::Track::getLastFoundHitLyr(), mkfit::Track::getLastFoundHitOnTrack(), and Electron.Electron::ptErr().

{ return state_.p(); }

parameters()

const SVector6& mkfit::TrackBase::parameters ( void  ) const

inline

Definition at line 132 of file Track.h.

References mkfit::TrackState::parameters, and state_.

Referenced by mkfit::MkFinder::copy_in(), mkfit::MkFinder::findCandidates(), mkfit::MkFitter::inputSeedsTracksAndHits(), mkfit::MkFitter::inputTracksAndHitIdx(), mkfit::MkFitter::inputTracksAndHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), and mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams().

{ return state_.parameters; }

parameters_nc()

SVector6& mkfit::TrackBase::parameters_nc ( )

inline

Definition at line 139 of file Track.h.

References mkfit::TrackState::parameters, and state_.

Referenced by mkfit::MkFinder::bkFitOutputTracks(), and mkfit::MkFinder::copy_out().

{ return state_.parameters; }

posArray()

const float* mkfit::TrackBase::posArray ( ) const

inline

Definition at line 135 of file Track.h.

References mkfit::TrackState::parameters, and state_.

{ return state_.parameters.Array(); }

posEta()

float mkfit::TrackBase::posEta ( ) const

inline

Definition at line 152 of file Track.h.

References mkfit::getEta(), mkfit::TrackState::parameters, and state_.

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

{ return getEta(state_.parameters[0], state_.parameters[1], state_.parameters[2]); }

position()

SVector3 mkfit::TrackBase::position ( ) const

inline

Definition at line 143 of file Track.h.

References mkfit::TrackState::parameters, and state_.

{ return SVector3(state_.parameters[0], state_.parameters[1], state_.parameters[2]); }

posPhi()

float mkfit::TrackBase::posPhi ( ) const

inline

Definition at line 151 of file Track.h.

References mkfit::getPhi(), mkfit::TrackState::parameters, and state_.

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

{ return getPhi(state_.parameters[0], state_.parameters[1]); }

posR()

float mkfit::TrackBase::posR ( ) const

inline

Definition at line 149 of file Track.h.

References mkfit::getHypot(), mkfit::TrackState::parameters, and state_.

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

{ return getHypot(state_.parameters[0], state_.parameters[1]); }

posRsq()

float mkfit::TrackBase::posRsq ( ) const

inline

Definition at line 150 of file Track.h.

References mkfit::TrackState::posRsq(), and state_.

{ return state_.posRsq(); }

prodType()

ProdType mkfit::TrackBase::prodType ( ) const

inline

Definition at line 258 of file Track.h.

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

Referenced by mkfit::MkFinder::bkFitFitTracks(), mkfit::TrackCand::exportTrack(), and mkfit::MkFinder::selectHitIndices().

{ return ProdType(status_.prod_type); }

pT()

float mkfit::TrackBase::pT ( ) const

inline

Definition at line 157 of file Track.h.

References mkfit::TrackState::pT(), and state_.

Referenced by mkfit::MkFinder::bkFitFitTracks(), canReachRadius(), mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), d0BeamSpot(), mkfit::MkBuilder::findTracksStandard(), mkfit::getScoreCand(), hasSillyValues(), maxReachRadius(), mkfit::CandCloner::processSeedRange(), mkfit::StdSeq::Quality::quality_process(), mkfit::MkFinder::selectHitIndices(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), and swimPhiToR().

{ return state_.pT(); }

px()

float mkfit::TrackBase::px ( ) const

inline

Definition at line 154 of file Track.h.

References mkfit::TrackState::px(), and state_.

Referenced by d0BeamSpot(), maxReachRadius(), rAtZ(), and zAtR().

{ return state_.px(); }

py()

float mkfit::TrackBase::py ( ) const

inline

Definition at line 155 of file Track.h.

References mkfit::TrackState::py(), and state_.

Referenced by d0BeamSpot(), maxReachRadius(), rAtZ(), and zAtR().

{ return state_.py(); }

pz()

float mkfit::TrackBase::pz ( ) const

inline

Definition at line 156 of file Track.h.

References mkfit::TrackState::pz(), and state_.

Referenced by zAtR().

{ return state_.pz(); }

rAtZ()

float mkfit::TrackBase::rAtZ

(

float

Z

)

const

Definition at line 305 of file Track.cc.

References simBeamSpotPI::alpha, mkfit::Config::Bfield, charge(), funct::cos(), PVValHelper::dz, f, Matriplex::hypot(), invpT(), dqmdumpme::k, px(), py(), funct::sin(), mkfit::Const::sol, funct::tan(), theta(), x(), y(), beamSpotPI::Z, and 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);
  }

score()

float mkfit::TrackBase::score ( ) const

inline

Definition at line 173 of file Track.h.

References score_.

Referenced by mkfit::MkBuilder::findTracksStandard(), mkfit::CombCandidate::mergeCandsAndBestShortOne(), mkfit::CandCloner::processSeedRange(), mkfit::sortByScoreCand(), and mkfit::sortByScoreTrackCand().

{ return score_; }

setAlgoint()

void mkfit::TrackBase::setAlgoint ( int  algo )

inline

Definition at line 330 of file Track.h.

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

{ status_.algorithm = algo; }

setAlgorithm()

void mkfit::TrackBase::setAlgorithm ( TrackAlgorithm  algo )

inline

Definition at line 329 of file Track.h.

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

{ status_.algorithm = static_cast<unsigned int>(algo); }

setCharge()

void mkfit::TrackBase::setCharge ( int  chg )

inline

Definition at line 176 of file Track.h.

References mkfit::TrackState::charge, chg, and state_.

Referenced by mkfit::MkFinder::copy_out(), and mkfit::MkFinder::findCandidates().

{ state_.charge = chg; }

setChi2()

void mkfit::TrackBase::setChi2 ( float  chi2 )

inline

Definition at line 177 of file Track.h.

References chi2(), and chi2_.

Referenced by mkfit::MkFinder::bkFitOutputTracks(), and mkfit::MkFinder::copy_out().

{ chi2_ = chi2; }

setDuplicateValue()

void mkfit::TrackBase::setDuplicateValue ( bool  d )

inline

Definition at line 255 of file Track.h.

References ztail::d, mkfit::TrackBase::Status::duplicate, and status_.

{ status_.duplicate = d; }

setEtaRegion()

void mkfit::TrackBase::setEtaRegion ( int  r )

inline

Definition at line 264 of file Track.h.

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

{ status_.eta_region = r; }

setLabel()

void mkfit::TrackBase::setLabel ( int  lbl )

inline

Definition at line 179 of file Track.h.

References label_.

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

{ label_ = lbl; }

setNotFindable()

void mkfit::TrackBase::setNotFindable ( )

inline

Definition at line 253 of file Track.h.

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

{ status_.not_findable = true; }

setNSeedHits()

void mkfit::TrackBase::setNSeedHits ( int  n )

inline

Definition at line 262 of file Track.h.

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

{ status_.n_seed_hits = n; }

setProdType()

void mkfit::TrackBase::setProdType ( ProdType  ptyp )

inline

Definition at line 259 of file Track.h.

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

{ status_.prod_type = static_cast<unsigned int>(ptyp); }

setScore()

void mkfit::TrackBase::setScore ( float  s )

inline

Definition at line 178 of file Track.h.

References alignCSCRings::s, and score_.

Referenced by mkfit::MkFinder::bkFitOutputTracks(), mkfit::MkFinder::findCandidates(), mkfit::CombCandidate::mergeCandsAndBestShortOne(), mkfit::CandCloner::processSeedRange(), and mkfit::CombCandidate::reset().

{ score_ = s; }

setState()

void mkfit::TrackBase::setState ( const TrackState &  newState )

inline

Definition at line 130 of file Track.h.

References state_.

{ state_ = newState; }

setStatus()

void mkfit::TrackBase::setStatus ( Status  s )

inline

Definition at line 249 of file Track.h.

References alignCSCRings::s, and status_.

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

{ status_ = s; }

state()

const TrackState& mkfit::TrackBase::state ( ) const

inline

Definition at line 129 of file Track.h.

References state_.

Referenced by mkfit::print().

{ return state_; }

state_nc()

TrackState& mkfit::TrackBase::state_nc ( )

inline

Definition at line 141 of file Track.h.

References state_.

{ return state_; }

swimPhiToR()

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

Definition at line 203 of file Track.cc.

References funct::abs(), charge(), HGC3DClusterGenMatchSelector_cfi::dR, mkfit::getHypot(), momPhi(), mkfit::Const::PI, pT(), mkfit::squashPhiGeneral(), findQualityFiles::v, x(), and 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);
  }

theta()

float mkfit::TrackBase::theta ( void  ) const

inline

Definition at line 162 of file Track.h.

References state_, and mkfit::TrackState::theta().

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), rAtZ(), and Tau.Tau::zImpact().

{ return state_.theta(); }

x()

float mkfit::TrackBase::x ( ) const

inline

Definition at line 146 of file Track.h.

References mkfit::TrackState::parameters, and state_.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), canReachRadius(), mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), d0BeamSpot(), maxReachRadius(), geometryXMLparser.Alignable::pos(), mkfit::StdSeq::Quality::quality_process(), ntupleDataFormat._HitObject::r(), ntupleDataFormat._HitObject::r3D(), rAtZ(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), mkfit::TrackExtra::setMCTrackIDInfo(), swimPhiToR(), and zAtR().

{ return state_.parameters[0]; }

y()

float mkfit::TrackBase::y ( ) const

inline

Definition at line 147 of file Track.h.

References mkfit::TrackState::parameters, and state_.

Referenced by svgfig.Ellipse::__repr__(), canReachRadius(), mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), d0BeamSpot(), maxReachRadius(), geometryXMLparser.Alignable::pos(), mkfit::StdSeq::Quality::quality_process(), ntupleDataFormat._HitObject::r(), ntupleDataFormat._HitObject::r3D(), rAtZ(), mkfit::TrackExtra::setCMSSWTrackIDInfoByHits(), mkfit::TrackExtra::setCMSSWTrackIDInfoByTrkParams(), mkfit::TrackExtra::setMCTrackIDInfo(), swimPhiToR(), and zAtR().

{ return state_.parameters[1]; }

z()

float mkfit::TrackBase::z ( ) const

inline

Definition at line 148 of file Track.h.

References mkfit::TrackState::parameters, and state_.

Referenced by mkfit::Event::clean_cms_seedtracks(), mkfit::StdSeq::clean_cms_seedtracks_iter(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r3D(), rAtZ(), mkfit::MkFinder::selectHitIndices(), and zAtR().

{ return state_.parameters[2]; }

zAtR()

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

Definition at line 228 of file Track.cc.

References mkfit::Config::Bfield, DummyCfis::c, charge(), funct::cos(), f, Matriplex::hypot(), mps_fire::i, l1ctLayer2EG_cff::id, invpT(), dqmdumpme::k, mkfit::Config::Niter, px(), py(), pz(), dttmaxenums::R, funct::sin(), mkfit::Const::sol, x(), y(), and 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

chi2_

float mkfit::TrackBase::chi2_ = 0.

protected

Definition at line 341 of file Track.h.

Referenced by mkfit::TrackCand::addHitIdx(), mkfit::Track::addHitIdx(), chi2(), mkfit::TrackCand::popOverlap(), and setChi2().

label_

int mkfit::TrackBase::label_ = -1

protected

Definition at line 346 of file Track.h.

Referenced by Modules.SwitchProducer::__addParameter(), SequenceTypes.DummyModule::__repr__(), Modules.SwitchProducer::__setattr__(), Mixins._Labelable::_findDependencies(), label(), setLabel(), and Mixins._Labelable::setLabel().

lastHitIdx_

short int mkfit::TrackBase::lastHitIdx_ = -1

protected

Definition at line 343 of file Track.h.

Referenced by mkfit::TrackCand::addHitIdx(), mkfit::Track::addHitIdx(), mkfit::Track::countAndSetNFoundHits(), mkfit::Track::endHitsOnTrack(), mkfit::TrackCand::exportTrack(), mkfit::Track::foundLayers(), mkfit::TrackCand::getLastFoundHitLyr(), mkfit::Track::getLastFoundHitPos(), mkfit::TrackCand::getLastFoundPixelHitLyr(), mkfit::TrackCand::getLastHitIdx(), mkfit::Track::getLastHitIdx(), mkfit::TrackCand::getLastHitLyr(), mkfit::Track::getLastHitLyr(), mkfit::TrackCand::getLastHitOnTrack(), mkfit::Track::getLastHitOnTrack(), mkfit::Track::getMCHitIDFromLayer(), mkfit::TrackCand::lastCcIndex(), mkfit::Track::mcHitIDsVec(), mkfit::TrackCand::nHitsByTypeEncoded(), mkfit::Track::nInsideMinusOneHits(), mkfit::TrackCand::nLayersByTypeEncoded(), mkfit::Track::nTailMinusOneHits(), mkfit::Track::nTotalHits(), mkfit::TrackCand::nUniqueLayers(), mkfit::Track::nUniqueLayers(), mkfit::TrackCand::popOverlap(), mkfit::TrackCand::refLastHoTNode(), mkfit::Track::resetHits(), mkfit::Track::resizeHits(), mkfit::Track::resizeHitsForInput(), mkfit::TrackCand::setLastCcIndex(), mkfit::Track::sortHitsByLayer(), mkfit::Track::Track(), and mkfit::TrackCand::TrackCand().

nFoundHits_

short int mkfit::TrackBase::nFoundHits_ = 0

protected

Definition at line 344 of file Track.h.

Referenced by mkfit::TrackCand::addHitIdx(), mkfit::Track::addHitIdx(), mkfit::Track::countAndSetNFoundHits(), mkfit::TrackCand::nFoundHits(), mkfit::Track::nFoundHits(), mkfit::TrackCand::nTotalHits(), mkfit::TrackCand::popOverlap(), mkfit::Track::resetHits(), mkfit::Track::resizeHits(), mkfit::TrackCand::setNFoundHits(), mkfit::Track::Track(), and mkfit::TrackCand::TrackCand().

score_

float mkfit::TrackBase::score_ = 0.

protected

Definition at line 342 of file Track.h.

Referenced by mkfit::TrackCand::resetShortTrack(), score(), and setScore().

state_

TrackState mkfit::TrackBase::state_

protected

Definition at line 340 of file Track.h.

Referenced by charge(), emomEta(), emomPhi(), epT(), errArray(), errors(), errors_nc(), hasNanNSillyValues(), hasSillyValues(), invpT(), momentum(), momEta(), momPhi(), p(), parameters(), parameters_nc(), posArray(), posEta(), position(), posPhi(), posR(), posRsq(), pT(), px(), py(), pz(), setCharge(), setState(), state(), state_nc(), theta(), x(), y(), and z().

status_

Status mkfit::TrackBase::status_

protected

Definition at line 345 of file Track.h.

Referenced by algoint(), algorithm(), getDuplicateValue(), getEtaRegion(), getNSeedHits(), getStatus(), isFindable(), isNotFindable(), mkfit::Track::nOverlapHits(), prodType(), setAlgoint(), setAlgorithm(), setDuplicateValue(), setEtaRegion(), setNotFindable(), mkfit::Track::setNOverlapHits(), setNSeedHits(), setProdType(), and setStatus().