mkfit::TrackState Struct Reference

#include <Track.h>

Public Member Functions
void	convertFromCartesianToCCS ()
void	convertFromCCSToCartesian ()
void	convertFromCCSToGlbCurvilinear ()
void	convertFromGlbCurvilinearToCCS ()
float	einvpT () const
float	emomEta () const
float	emomPhi () const
float	eposEta () const
float	eposPhi () const
float	eposR () const
float	epT () const
float	epxpx () const
float	epypy () const
float	epzpz () const
float	etheta () const
float	exx () const
float	exy () const
float	exz () const
float	eyy () const
float	eyz () const
float	ezz () const
float	invpT () const
SMatrix66	jacobianCartesianToCCS (float px, float py, float pz) const
SMatrix66	jacobianCCSToCartesian (float invpt, float phi, float theta) const
SMatrix66	jacobianCCSToCurvilinear (float invpt, float cosP, float sinP, float cosT, float sinT, short charge) const
SMatrix66	jacobianCurvilinearToCCS (float px, float py, float pz, short charge) const
float	momEta () const
float	momPhi () const
float	p () const
float	posEta () const
SVector3	position () const
float	posPhi () const
float	posR () const
float	posRsq () const
float	pT () const
float	px () const
float	py () const
float	pz () const
float	theta () const
	TrackState ()
	TrackState (int charge, const SVector3 &pos, const SVector3 &mom, const SMatrixSym66 &err)
float	x () const
float	y () const
float	z () const

Public Attributes
short	charge
SMatrixSym66	errors
SVector6	parameters
bool	valid

Detailed Description

Definition at line 50 of file Track.h.

Constructor & Destructor Documentation

TrackState() [1/2]

mkfit::TrackState::TrackState ( )

inline

Definition at line 53 of file Track.h.

: valid(true) {}

TrackState() [2/2]

mkfit::TrackState::TrackState ( int  charge, const SVector3 &  pos, const SVector3 &  mom, const SMatrixSym66 &  err  )

inline

Definition at line 54 of file Track.h.

        : parameters(SVector6(pos.At(0), pos.At(1), pos.At(2), mom.At(0), mom.At(1), mom.At(2))),
          errors(err),
          charge(charge),
          valid(true) {}

Member Function Documentation

convertFromCartesianToCCS()

void mkfit::TrackState::convertFromCartesianToCCS

(

)

Definition at line 13 of file Track.cc.

References mkfit::getPhi(), mkfit::getTheta(), jacobianCartesianToCCS(), parameters, DiDispStaMuonMonitor_cfi::pt, px(), py(), pz(), mathSSE::sqrt(), and theta().

                                             {
    //assume we are currently in cartesian coordinates and want to move to ccs
    const float px = parameters.At(3);
    const float py = parameters.At(4);
    const float pz = parameters.At(5);
    const float pt = std::sqrt(px * px + py * py);
    const float phi = getPhi(px, py);
    const float theta = getTheta(pt, pz);
    parameters.At(3) = 1.f / pt;
    parameters.At(4) = phi;
    parameters.At(5) = theta;
    SMatrix66 jac = jacobianCartesianToCCS(px, py, pz);
    errors = ROOT::Math::Similarity(jac, errors);
  }

convertFromCCSToCartesian()

void mkfit::TrackState::convertFromCCSToCartesian

(

)

Definition at line 28 of file Track.cc.

References funct::cos(), jacobianCCSToCartesian(), parameters, DiDispStaMuonMonitor_cfi::pt, funct::sin(), and theta().

                                             {
    //assume we are currently in ccs coordinates and want to move to cartesian
    const float invpt = parameters.At(3);
    const float phi = parameters.At(4);
    const float theta = parameters.At(5);
    const float pt = 1.f / invpt;
    float cosP = std::cos(phi);
    float sinP = std::sin(phi);
    float cosT = std::cos(theta);
    float sinT = std::sin(theta);
    parameters.At(3) = cosP * pt;
    parameters.At(4) = sinP * pt;
    parameters.At(5) = cosT * pt / sinT;
    SMatrix66 jac = jacobianCCSToCartesian(invpt, phi, theta);
    errors = ROOT::Math::Similarity(jac, errors);
  }

convertFromCCSToGlbCurvilinear()

void mkfit::TrackState::convertFromCCSToGlbCurvilinear

(

)

Definition at line 92 of file Track.cc.

References charge, funct::cos(), jacobianCCSToCurvilinear(), parameters, DiDispStaMuonMonitor_cfi::pt, funct::sin(), and theta().

                                                  {
    //assume we are currently in ccs coordinates and want to move to global state with cartesian error
    const float invpt = parameters.At(3);
    const float phi = parameters.At(4);
    const float theta = parameters.At(5);
    const float pt = 1.f / invpt;
    float cosP = std::cos(phi);
    float sinP = std::sin(phi);
    float cosT = std::cos(theta);
    float sinT = std::sin(theta);
    parameters.At(3) = cosP * pt;
    parameters.At(4) = sinP * pt;
    parameters.At(5) = cosT * pt / sinT;
    SMatrix66 jac = jacobianCCSToCurvilinear(invpt, cosP, sinP, cosT, sinT, charge);
    errors = ROOT::Math::Similarity(jac, errors);
  }

convertFromGlbCurvilinearToCCS()

void mkfit::TrackState::convertFromGlbCurvilinearToCCS

(

)

Definition at line 77 of file Track.cc.

References charge, mkfit::getPhi(), mkfit::getTheta(), jacobianCurvilinearToCCS(), parameters, DiDispStaMuonMonitor_cfi::pt, px(), py(), pz(), mathSSE::sqrt(), and theta().

                                                  {
    //assume we are currently in global state with curvilinear error and want to move to ccs
    const float px = parameters.At(3);
    const float py = parameters.At(4);
    const float pz = parameters.At(5);
    const float pt = std::sqrt(px * px + py * py);
    const float phi = getPhi(px, py);
    const float theta = getTheta(pt, pz);
    parameters.At(3) = 1.f / pt;
    parameters.At(4) = phi;
    parameters.At(5) = theta;
    SMatrix66 jac = jacobianCurvilinearToCCS(px, py, pz, charge);
    errors = ROOT::Math::Similarity(jac, errors);
  }

einvpT()

float mkfit::TrackState::einvpT ( ) const

inline

Definition at line 107 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(3, 3)); }

emomEta()

float mkfit::TrackState::emomEta ( ) const

inline

Definition at line 111 of file Track.h.

References parameters, funct::sin(), and mathSSE::sqrt().

Referenced by mkfit::TrackBase::emomEta().

{ return std::sqrt(errors.At(5, 5)) / std::sin(parameters.At(5)); }

emomPhi()

float mkfit::TrackState::emomPhi ( ) const

inline

Definition at line 108 of file Track.h.

References mathSSE::sqrt().

Referenced by mkfit::TrackBase::emomPhi().

{ return std::sqrt(errors.At(4, 4)); }

eposEta()

float mkfit::TrackState::eposEta ( ) const

inline

Definition at line 84 of file Track.h.

References mkfit::getEtaErr2(), mathSSE::sqrt(), x(), y(), and z().

                          {
      return std::sqrt(getEtaErr2(x(),
                                  y(),
                                  z(),
                                  errors.At(0, 0),
                                  errors.At(1, 1),
                                  errors.At(2, 2),
                                  errors.At(0, 1),
                                  errors.At(0, 2),
                                  errors.At(1, 2)));
    }

eposPhi()

float mkfit::TrackState::eposPhi ( ) const

inline

Definition at line 83 of file Track.h.

References mkfit::getPhiErr2(), mathSSE::sqrt(), x(), and y().

{ return std::sqrt(getPhiErr2(x(), y(), errors.At(0, 0), errors.At(1, 1), errors.At(0, 1))); }

eposR()

float mkfit::TrackState::eposR ( ) const

inline

Definition at line 82 of file Track.h.

References mkfit::getRadErr2(), mathSSE::sqrt(), x(), and y().

{ return std::sqrt(getRadErr2(x(), y(), errors.At(0, 0), errors.At(1, 1), errors.At(0, 1))); }

epT()

float mkfit::TrackState::epT ( ) const

inline

Definition at line 110 of file Track.h.

References parameters, and mathSSE::sqrt().

Referenced by mkfit::TrackBase::epT().

{ return std::sqrt(errors.At(3, 3)) / (parameters.At(3) * parameters.At(3)); }

epxpx()

float mkfit::TrackState::epxpx ( ) const

inline

Definition at line 112 of file Track.h.

References mkfit::getPxPxErr2(), invpT(), momPhi(), and mathSSE::sqrt().

{ return std::sqrt(getPxPxErr2(invpT(), momPhi(), errors.At(3, 3), errors.At(4, 4))); }

epypy()

float mkfit::TrackState::epypy ( ) const

inline

Definition at line 113 of file Track.h.

References mkfit::getPyPyErr2(), invpT(), momPhi(), and mathSSE::sqrt().

{ return std::sqrt(getPyPyErr2(invpT(), momPhi(), errors.At(3, 3), errors.At(4, 4))); }

epzpz()

float mkfit::TrackState::epzpz ( ) const

inline

Definition at line 114 of file Track.h.

References mkfit::getPyPyErr2(), invpT(), mathSSE::sqrt(), and theta().

{ return std::sqrt(getPyPyErr2(invpT(), theta(), errors.At(3, 3), errors.At(5, 5))); }

etheta()

float mkfit::TrackState::etheta ( ) const

inline

Definition at line 109 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(5, 5)); }

exx()

float mkfit::TrackState::exx ( ) const

inline

Definition at line 75 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(0, 0)); }

exy()

float mkfit::TrackState::exy ( ) const

inline

Definition at line 78 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(0, 1)); }

exz()

float mkfit::TrackState::exz ( ) const

inline

Definition at line 79 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(0, 2)); }

eyy()

float mkfit::TrackState::eyy ( ) const

inline

Definition at line 76 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(1, 1)); }

eyz()

float mkfit::TrackState::eyz ( ) const

inline

Definition at line 80 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(1, 2)); }

ezz()

float mkfit::TrackState::ezz ( ) const

inline

Definition at line 77 of file Track.h.

References mathSSE::sqrt().

{ return std::sqrt(errors.At(2, 2)); }

invpT()

float mkfit::TrackState::invpT ( ) const

inline

Definition at line 97 of file Track.h.

References parameters.

Referenced by epxpx(), epypy(), epzpz(), and mkfit::TrackBase::invpT().

{ return parameters.At(3); }

jacobianCartesianToCCS()

SMatrix66 mkfit::TrackState::jacobianCartesianToCCS	(	float	px,
		float	py,
		float	pz
	)		const

Definition at line 62 of file Track.cc.

References SiStripOfflineCRack_cfg::p2, DiDispStaMuonMonitor_cfi::pt, px(), py(), pz(), and mathSSE::sqrt().

Referenced by convertFromCartesianToCCS().

                                                                                 {
    //arguments are passed so that the function can be used both starting from ccs and from cartesian
    SMatrix66 jac = ROOT::Math::SMatrixIdentity();
    const float pt = std::sqrt(px * px + py * py);
    const float p2 = px * px + py * py + pz * pz;
    jac(3, 3) = -px / (pt * pt * pt);
    jac(3, 4) = -py / (pt * pt * pt);
    jac(4, 3) = -py / (pt * pt);
    jac(4, 4) = px / (pt * pt);
    jac(5, 3) = px * pz / (pt * p2);
    jac(5, 4) = py * pz / (pt * p2);
    jac(5, 5) = -pt / p2;
    return jac;
  }

jacobianCCSToCartesian()

SMatrix66 mkfit::TrackState::jacobianCCSToCartesian	(	float	invpt,
		float	phi,
		float	theta
	)		const

Definition at line 45 of file Track.cc.

References funct::cos(), DiDispStaMuonMonitor_cfi::pt, funct::sin(), and theta().

Referenced by convertFromCCSToCartesian().

                                                                                        {
    //arguments are passed so that the function can be used both starting from ccs and from cartesian
    SMatrix66 jac = ROOT::Math::SMatrixIdentity();
    float cosP = std::cos(phi);
    float sinP = std::sin(phi);
    float cosT = std::cos(theta);
    float sinT = std::sin(theta);
    const float pt = 1.f / invpt;
    jac(3, 3) = -cosP * pt * pt;
    jac(3, 4) = -sinP * pt;
    jac(4, 3) = -sinP * pt * pt;
    jac(4, 4) = cosP * pt;
    jac(5, 3) = -cosT * pt * pt / sinT;
    jac(5, 5) = -pt / (sinT * sinT);
    return jac;
  }

jacobianCCSToCurvilinear()

SMatrix66 mkfit::TrackState::jacobianCCSToCurvilinear	(	float	invpt,
		float	cosP,
		float	sinP,
		float	cosT,
		float	sinT,
		short	charge
	)		const

Definition at line 109 of file Track.cc.

References charge.

Referenced by convertFromCCSToGlbCurvilinear().

                                                                                       {
    SMatrix66 jac;
    jac(3, 0) = -sinP;
    jac(4, 0) = -cosP * cosT;
    jac(3, 1) = cosP;
    jac(4, 1) = -sinP * cosT;
    jac(4, 2) = sinT;
    jac(0, 3) = charge * sinT;
    jac(0, 5) = charge * cosT * invpt;
    jac(1, 5) = -1.f;
    jac(2, 4) = 1.f;

    return jac;
  }

jacobianCurvilinearToCCS()

SMatrix66 mkfit::TrackState::jacobianCurvilinearToCCS	(	float	px,
		float	py,
		float	pz,
		short	charge
	)		const

Definition at line 125 of file Track.cc.

References charge, l1tSlwPFPuppiJets_cfi::cosPhi, DiDispStaMuonMonitor_cfi::pt, HLT_2022v12_cff::pt2, px(), py(), pz(), l1tSlwPFPuppiJets_cfi::sinPhi, and mathSSE::sqrt().

Referenced by convertFromGlbCurvilinearToCCS().

                                                                                                 {
    const float pt2 = px * px + py * py;
    const float pt = sqrt(pt2);
    const float invpt2 = 1.f / pt2;
    const float invpt = 1.f / pt;
    const float invp = 1.f / sqrt(pt2 + pz * pz);
    const float sinPhi = py * invpt;
    const float cosPhi = px * invpt;
    const float sinLam = pz * invp;
    const float cosLam = pt * invp;

    SMatrix66 jac;
    jac(0, 3) = -sinPhi;
    jac(0, 4) = -sinLam * cosPhi;
    jac(1, 3) = cosPhi;
    jac(1, 4) = -sinLam * sinPhi;
    jac(2, 4) = cosLam;
    jac(3, 0) = charge / cosLam;  //assumes |charge|==1 ; else 1.f/charge here
    jac(3, 1) = pz * invpt2;
    jac(4, 2) = 1.f;
    jac(5, 1) = -1.f;

    return jac;
  }

momEta()

float mkfit::TrackState::momEta ( ) const

inline

Definition at line 104 of file Track.h.

References mkfit::getEta(), and theta().

Referenced by mkfit::TTreeValidation::fillEfficiencyTree(), mkfit::TTreeValidation::fillFakeRateTree(), and mkfit::TrackBase::momEta().

{ return getEta(theta()); }

momPhi()

float mkfit::TrackState::momPhi ( ) const

inline

Definition at line 98 of file Track.h.

References parameters.

Referenced by epxpx(), epypy(), mkfit::TTreeValidation::fillEfficiencyTree(), mkfit::TTreeValidation::fillFakeRateTree(), and mkfit::TrackBase::momPhi().

{ return parameters.At(4); }

p()

float mkfit::TrackState::p ( ) const

inline

Definition at line 105 of file Track.h.

References parameters, pT(), and funct::sin().

Referenced by mkfit::TrackBase::p(), and Electron.Electron::ptErr().

{ return pT() / std::sin(parameters.At(5)); }

posEta()

float mkfit::TrackState::posEta ( ) const

inline

Definition at line 72 of file Track.h.

References mkfit::getEta(), posR(), and z().

{ return getEta(posR(), z()); }

position()

SVector3 mkfit::TrackState::position ( ) const

inline

Definition at line 59 of file Track.h.

References parameters.

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

posPhi()

float mkfit::TrackState::posPhi ( ) const

inline

Definition at line 71 of file Track.h.

References mkfit::getPhi(), x(), and y().

{ return getPhi(x(), y()); }

posR()

float mkfit::TrackState::posR ( ) const

inline

Definition at line 69 of file Track.h.

References mkfit::getHypot(), x(), and y().

Referenced by posEta().

{ return getHypot(x(), y()); }

posRsq()

float mkfit::TrackState::posRsq ( ) const

inline

Definition at line 70 of file Track.h.

References x(), and y().

Referenced by mkfit::TrackBase::posRsq().

{ return x() * x() + y() * y(); }

pT()

float mkfit::TrackState::pT ( ) const

inline

Definition at line 100 of file Track.h.

References funct::abs(), f, and parameters.

Referenced by mkfit::TTreeValidation::fillEfficiencyTree(), mkfit::TTreeValidation::fillFakeRateTree(), p(), mkfit::TrackBase::pT(), px(), py(), and pz().

{ return std::abs(1.f / parameters.At(3)); }

px()

float mkfit::TrackState::px ( ) const

inline

Definition at line 101 of file Track.h.

References funct::cos(), parameters, and pT().

Referenced by convertFromCartesianToCCS(), convertFromGlbCurvilinearToCCS(), jacobianCartesianToCCS(), jacobianCurvilinearToCCS(), and mkfit::TrackBase::px().

{ return pT() * std::cos(parameters.At(4)); }

py()

float mkfit::TrackState::py ( ) const

inline

Definition at line 102 of file Track.h.

References parameters, pT(), and funct::sin().

Referenced by convertFromCartesianToCCS(), convertFromGlbCurvilinearToCCS(), jacobianCartesianToCCS(), jacobianCurvilinearToCCS(), and mkfit::TrackBase::py().

{ return pT() * std::sin(parameters.At(4)); }

pz()

float mkfit::TrackState::pz ( ) const

inline

Definition at line 103 of file Track.h.

References parameters, pT(), and funct::tan().

Referenced by convertFromCartesianToCCS(), convertFromGlbCurvilinearToCCS(), jacobianCartesianToCCS(), jacobianCurvilinearToCCS(), and mkfit::TrackBase::pz().

{ return pT() / std::tan(parameters.At(5)); }

theta()

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

inline

Definition at line 99 of file Track.h.

References parameters.

Referenced by convertFromCartesianToCCS(), convertFromCCSToCartesian(), convertFromCCSToGlbCurvilinear(), convertFromGlbCurvilinearToCCS(), epzpz(), jacobianCCSToCartesian(), momEta(), mkfit::TrackBase::theta(), and Tau.Tau::zImpact().

{ return parameters.At(5); }

x()

float mkfit::TrackState::x ( ) const

inline

Definition at line 66 of file Track.h.

References parameters.

Referenced by svgfig.Curve.Sample::__repr__(), svgfig.Ellipse::__repr__(), eposEta(), eposPhi(), eposR(), geometryXMLparser.Alignable::pos(), posPhi(), posR(), posRsq(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

{ return parameters.At(0); }

y()

float mkfit::TrackState::y ( ) const

inline

Definition at line 67 of file Track.h.

References parameters.

Referenced by svgfig.Ellipse::__repr__(), eposEta(), eposPhi(), eposR(), geometryXMLparser.Alignable::pos(), posPhi(), posR(), posRsq(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

{ return parameters.At(1); }

z()

float mkfit::TrackState::z ( ) const

inline

Definition at line 68 of file Track.h.

References parameters.

Referenced by eposEta(), geometryXMLparser.Alignable::pos(), posEta(), and ntupleDataFormat._HitObject::r3D().

{ return parameters.At(2); }

Member Data Documentation

charge

short mkfit::TrackState::charge

Definition at line 62 of file Track.h.

Referenced by mkfit::TrackBase::charge(), convertFromCCSToGlbCurvilinear(), convertFromGlbCurvilinearToCCS(), jacobianCCSToCurvilinear(), jacobianCurvilinearToCCS(), and mkfit::TrackBase::setCharge().

errors

SMatrixSym66 mkfit::TrackState::errors

Definition at line 61 of file Track.h.

Referenced by mkfit::TrackBase::errArray(), mkfit::TrackBase::errors(), mkfit::TrackBase::errors_nc(), mkfit::TrackBase::hasNanNSillyValues(), and mkfit::TrackBase::hasSillyValues().

parameters

SVector6 mkfit::TrackState::parameters

Definition at line 60 of file Track.h.

Referenced by mkfit::conformalFitMPlex(), convertFromCartesianToCCS(), convertFromCCSToCartesian(), convertFromCCSToGlbCurvilinear(), convertFromGlbCurvilinearToCCS(), emomEta(), epT(), mkfit::TTreeValidation::fillEfficiencyTree(), mkfit::TTreeValidation::fillFakeRateTree(), invpT(), mkfit::TrackBase::momentum(), momPhi(), p(), mkfit::TrackBase::parameters(), mkfit::TrackBase::parameters_nc(), mkfit::TrackBase::posArray(), mkfit::TrackBase::posEta(), position(), mkfit::TrackBase::position(), mkfit::TrackBase::posPhi(), mkfit::TrackBase::posR(), pT(), px(), py(), pz(), theta(), x(), mkfit::TrackBase::x(), y(), mkfit::TrackBase::y(), z(), and mkfit::TrackBase::z().

valid

bool mkfit::TrackState::valid

Definition at line 63 of file Track.h.