Trk::FitParameters Class Reference

#include <FitParameters.h>

Collaboration diagram for Trk::FitParameters:

Public Member Functions
	FitParameters (const Perigee &perigee)
	FitParameters (double d0, double z0, double cosPhi, double sinPhi, double cotTheta, double ptInv0, const PerigeeSurface &surface)
	FitParameters (const FitParameters &parameters)=default
	FitParameters (FitParameters &&)=default
FitParameters &	operator= (const FitParameters &)=default
FitParameters &	operator= (FitParameters &&)=default
	~FitParameters (void)=default
void	addAlignment (bool constrained, double localAngle, double localOffset)
void	addScatterer (double phi, double theta)
double	alignmentAngle (int alignment) const
double	alignmentAngleConstraint (int alignment) const
double	alignmentOffset (int alignment) const
double	alignmentOffsetConstraint (int alignment) const
const Surface *	associatedSurface (void) const
double	cosPhi (void) const
double	cosTheta (void) const
double	cotTheta (void) const
void	covariance (Amg::MatrixX *finalCovariance, const Amg::MatrixX *fullCovariance)
double	d0 (void) const
void	d0 (double value)
double	difference (int param) const
const Amg::VectorX &	differences (void) const
Amg::Vector3D	direction (void) const
bool	extremeMomentum (void) const
void	extremeMomentum (bool value)
const Amg::MatrixX *	finalCovariance (void) const
int	firstAlignmentParameter (void) const
void	firstAlignmentParameter (int value)
int	firstScatteringParameter (void) const
void	firstScatteringParameter (int value)
bool	fitEnergyDeposit (void) const
void	fitEnergyDeposit (double minEnergyDeposit)
bool	fitMomentum (void) const
void	fitMomentum (bool value)
const Amg::MatrixX *	fullCovariance (void) const
TrackSurfaceIntersection	intersection (void) const
int	numberAlignments (void) const
void	numberAlignments (int numberAlignments)
int	numberOscillations (void) const
int	numberParameters (void) const
void	numberParameters (int numberParameters)
int	numberScatterers (void) const
void	numberScatterers (int numberScatterers)
void	performCutStep (double cutStep)
bool	phiInstability (void) const
double	ptInv0 (void) const
const Amg::MatrixX	parameterDifference (const Amg::VectorX &parameters) const
Perigee *	perigee (void) const
const Amg::Vector3D &	position (void) const
void	print (MsgStream &log) const
void	printCovariance (MsgStream &log) const
void	printVerbose (MsgStream &log) const
double	qOverP (void) const
void	qOverP (double value)
double	qOverP1 (void) const
void	qOverP1 (double value)
void	reset (const FitParameters &parameters)
void	resetOscillations (void)
double	scattererPhi (int scatterer) const
double	scattererTheta (int scatterer) const
ScatteringAngles	scatteringAngles (const FitMeasurement &fitMeasurement, int scatterer=-1) const
void	setPhiInstability (void)
double	sinPhi (void) const
double	sinTheta (void) const
Perigee *	startingPerigee (void) const
TrackParameters *	trackParameters (MsgStream &log, const FitMeasurement &measurement, bool withCovariance=false)
void	update (const Amg::VectorX &differences)
void	update (Amg::Vector3D position, Amg::Vector3D direction, double qOverP, const Amg::MatrixX &leadingCovariance)
const Amg::Vector3D &	vertex (void) const
double	z0 (void) const

Private Attributes
std::vector< double >	m_alignmentAngle
std::vector< double >	m_alignmentAngleConstraint
std::vector< double >	m_alignmentOffset
std::vector< double >	m_alignmentOffsetConstraint
double	m_cosPhi
double	m_cosPhi1
double	m_cosTheta
double	m_cosTheta1
double	m_cotTheta
double	m_d0
Amg::VectorX	m_differences {}
bool	m_eigen {}
bool	m_extremeMomentum
Amg::MatrixX *	m_finalCovariance
int	m_firstAlignmentParameter
int	m_firstScatteringParameter
bool	m_fitEnergyDeposit
bool	m_fitMomentum
const Amg::MatrixX *	m_fullCovariance
double	m_minEnergyDeposit
int	m_numberAlignments
int	m_numberOscillations
int	m_numberParameters
int	m_numberScatterers
double	m_oldDifference
const Perigee *	m_perigee
bool	m_phiInstability
Amg::Vector3D	m_position
double	m_qOverP
double	m_qOverP1
std::vector< double >	m_scattererPhi
std::vector< double >	m_scattererTheta
double	m_sinPhi
double	m_sinPhi1
double	m_sinTheta
double	m_sinTheta1
const Surface *	m_surface
Amg::Vector3D	m_vertex
double	m_z0

Detailed Description

Definition at line 29 of file FitParameters.h.

Constructor & Destructor Documentation

FitParameters() [1/4]

Trk::FitParameters::FitParameters

(

const Perigee &

perigee

)

Definition at line 26 of file FitParameters.cxx.

    : m_cosPhi1(0.),
      m_cosTheta1(0.),
      m_d0(perigee.parameters()[Trk::d0]),
      m_differences{},
      m_extremeMomentum(false),
      m_finalCovariance(nullptr),
      m_firstAlignmentParameter(0),
      m_firstScatteringParameter(0),
      m_fitEnergyDeposit(false),
      m_fitMomentum(true),
      m_fullCovariance(nullptr),
      m_minEnergyDeposit(0.),
      m_numberAlignments(0),
      m_numberOscillations(0),
      m_numberParameters(0),
      m_numberScatterers(0),
      m_oldDifference(0.),
      m_perigee(&perigee),
      m_phiInstability(false),
      m_qOverP1(0.),
      m_sinPhi1(0.),
      m_sinTheta1(0.),
      m_surface(&perigee.associatedSurface()),
      m_vertex(perigee.associatedSurface().center()),
      m_z0(perigee.position().z()) {
  Amg::Vector3D momentum = perigee.momentum();
  const double ptInv0 = 1. / momentum.perp();
  m_cosPhi = ptInv0 * momentum.x();
  m_sinPhi = ptInv0 * momentum.y();
  m_cotTheta = ptInv0 * momentum.z();
  m_sinTheta = 1. / std::sqrt(1. + m_cotTheta * m_cotTheta);
  m_cosTheta = m_cotTheta * m_sinTheta;
  m_qOverP = perigee.charge() * ptInv0 * m_sinTheta;
  m_position = Amg::Vector3D(m_vertex.x() - m_d0 * m_sinPhi,
                             m_vertex.y() + m_d0 * m_cosPhi, m_z0);
}

FitParameters() [2/4]

Trk::FitParameters::FitParameters	(	double	d0,
		double	z0,
		double	cosPhi,
		double	sinPhi,
		double	cotTheta,
		double	ptInv0,
		const PerigeeSurface &	surface )

Definition at line 64 of file FitParameters.cxx.

    : m_cosPhi(cosPhi),
      m_cosPhi1(0.),
      m_cosTheta1(0.),
      m_cotTheta(cotTheta),
      m_d0(d0),
      m_extremeMomentum(false),
      m_finalCovariance(nullptr),
      m_firstAlignmentParameter(0),
      m_firstScatteringParameter(0),
      m_fitEnergyDeposit(false),
      m_fitMomentum(true),
      m_fullCovariance(nullptr),
      m_minEnergyDeposit(0.),
      m_numberAlignments(0),
      m_numberOscillations(0),
      m_numberParameters(0),
      m_numberScatterers(0),
      m_oldDifference(0.),
      m_perigee(nullptr),
      m_phiInstability(false),
      m_qOverP1(0.),
      m_sinPhi(sinPhi),
      m_sinPhi1(0.),
      m_sinTheta1(0.),
      m_surface(&surface),
      m_vertex(surface.center()),
      m_z0(z0) {
  m_sinTheta = 1. / std::sqrt(1. + m_cotTheta * m_cotTheta);
  m_cosTheta = m_cotTheta * m_sinTheta;
  m_qOverP = ptInv0 * m_sinTheta;
  m_z0 = z0;
  m_position = Amg::Vector3D(m_vertex.x() - m_d0 * m_sinPhi,
                             m_vertex.y() + m_d0 * m_cosPhi, m_z0);
}

FitParameters() [3/4]

Trk::FitParameters::FitParameters ( const FitParameters & parameters )

default

FitParameters() [4/4]

Trk::FitParameters::FitParameters ( FitParameters && )

default

~FitParameters()

Trk::FitParameters::~FitParameters ( void )

default

Member Function Documentation

addAlignment()

void Trk::FitParameters::addAlignment	(	bool	constrained,
		double	localAngle,
		double	localOffset )

Definition at line 103 of file FitParameters.cxx.

                                                {
  m_alignmentAngle[m_numberAlignments] = angle;
  m_alignmentOffset[m_numberAlignments] = offset;
  if (constrained) {
    m_alignmentAngleConstraint[m_numberAlignments] = angle;
    m_alignmentOffsetConstraint[m_numberAlignments] = offset;
  }
  ++m_numberAlignments;
}

addScatterer()

void Trk::FitParameters::addScatterer	(	double	phi,
		double	theta )

Definition at line 114 of file FitParameters.cxx.

                                                         {
  m_scattererPhi[m_numberScatterers] = phi;
  m_scattererTheta[m_numberScatterers] = theta;
  ++m_numberScatterers;
}

alignmentAngle()

double Trk::FitParameters::alignmentAngle ( int alignment ) const

inline

Definition at line 156 of file FitParameters.h.

                                                               {
  return m_alignmentAngle[alignment];
}

alignmentAngleConstraint()

double Trk::FitParameters::alignmentAngleConstraint ( int alignment ) const

inline

Definition at line 160 of file FitParameters.h.

                                                                         {
  return m_alignmentAngleConstraint[alignment];
}

alignmentOffset()

double Trk::FitParameters::alignmentOffset ( int alignment ) const

inline

Definition at line 164 of file FitParameters.h.

                                                                {
  return m_alignmentOffset[alignment];
}

alignmentOffsetConstraint()

double Trk::FitParameters::alignmentOffsetConstraint ( int alignment ) const

inline

Definition at line 168 of file FitParameters.h.

                                                                          {
  return m_alignmentOffsetConstraint[alignment];
}

associatedSurface()

const Surface * Trk::FitParameters::associatedSurface

(

void

)

const

Definition at line 120 of file FitParameters.cxx.

                                                          {
  if (!m_perigee)
    return nullptr;
  return &m_perigee->associatedSurface();
}

cosPhi()

double Trk::FitParameters::cosPhi ( void ) const

inline

Definition at line 172 of file FitParameters.h.

                                              {
  return m_cosPhi;
}

cosTheta()

double Trk::FitParameters::cosTheta ( void ) const

inline

Definition at line 176 of file FitParameters.h.

                                                {
  return m_cosTheta;
}

cotTheta()

double Trk::FitParameters::cotTheta ( void ) const

inline

Definition at line 180 of file FitParameters.h.

                                                {
  return m_cotTheta;
}

covariance()

void Trk::FitParameters::covariance	(	Amg::MatrixX *	finalCovariance,
		const Amg::MatrixX *	fullCovariance )

Definition at line 126 of file FitParameters.cxx.

                                                                 {
  m_finalCovariance = finalCovariance;
  m_fullCovariance = fullCovariance;
}

d0() [1/2]

void Trk::FitParameters::d0

(

double

value

)

Definition at line 132 of file FitParameters.cxx.

                                   {
  m_d0 = value;
  m_position = Amg::Vector3D(m_vertex.x() - m_d0 * m_sinPhi,
                             m_vertex.y() + m_d0 * m_cosPhi, m_z0);
}

d0() [2/2]

double Trk::FitParameters::d0 ( void ) const

inline

Definition at line 184 of file FitParameters.h.

                                          {
  return m_d0;
}

difference()

double Trk::FitParameters::difference ( int param ) const

inline

Use as_const to avoid compiler warning

Definition at line 188 of file FitParameters.h.

                                                       {
  if (m_differences.size() == 0) {
    return 0.;
  }
  return std::as_const(m_differences)(param);
}

differences()

const Amg::VectorX & Trk::FitParameters::differences ( void ) const

inline

Definition at line 196 of file FitParameters.h.

                                                              {
  return m_differences;
}

direction()

Amg::Vector3D Trk::FitParameters::direction ( void ) const

inline

Definition at line 200 of file FitParameters.h.

                                                      {
  return Amg::Vector3D(m_cosPhi * m_sinTheta, m_sinPhi * m_sinTheta,
                       m_cosTheta);
}

extremeMomentum() [1/2]

void Trk::FitParameters::extremeMomentum

(

bool

value

)

Definition at line 138 of file FitParameters.cxx.

                                              {
  m_extremeMomentum = value;
  m_fitEnergyDeposit = !value;
}

extremeMomentum() [2/2]

bool Trk::FitParameters::extremeMomentum ( void ) const

inline

Definition at line 205 of file FitParameters.h.

                                                     {
  return m_extremeMomentum;
}

finalCovariance()

const Amg::MatrixX * Trk::FitParameters::finalCovariance ( void ) const

inline

Definition at line 209 of file FitParameters.h.

                                                                  {
  return m_finalCovariance;
}

firstAlignmentParameter() [1/2]

void Trk::FitParameters::firstAlignmentParameter

(

int

value

)

firstAlignmentParameter() [2/2]

int Trk::FitParameters::firstAlignmentParameter ( void ) const

inline

Definition at line 213 of file FitParameters.h.

                                                            {
  return m_firstAlignmentParameter;
}

firstScatteringParameter() [1/2]

void Trk::FitParameters::firstScatteringParameter

(

int

value

)

firstScatteringParameter() [2/2]

int Trk::FitParameters::firstScatteringParameter ( void ) const

inline

Definition at line 217 of file FitParameters.h.

                                                             {
  return m_firstScatteringParameter;
}

fitEnergyDeposit() [1/2]

void Trk::FitParameters::fitEnergyDeposit

(

double

minEnergyDeposit

)

Definition at line 143 of file FitParameters.cxx.

                                                            {
  m_fitEnergyDeposit = true;
  m_minEnergyDeposit = minEnergyDeposit;
}

fitEnergyDeposit() [2/2]

bool Trk::FitParameters::fitEnergyDeposit ( void ) const

inline

Definition at line 221 of file FitParameters.h.

                                                      {
  return m_fitEnergyDeposit;
}

fitMomentum() [1/2]

void Trk::FitParameters::fitMomentum

(

bool

value

)

Definition at line 148 of file FitParameters.cxx.

                                          {
  m_fitMomentum = value;
}

fitMomentum() [2/2]

bool Trk::FitParameters::fitMomentum ( void ) const

inline

Definition at line 225 of file FitParameters.h.

                                                 {
  return m_fitMomentum;
}

fullCovariance()

const Amg::MatrixX * Trk::FitParameters::fullCovariance ( void ) const

inline

Definition at line 229 of file FitParameters.h.

                                                                 {
  return m_fullCovariance;
}

intersection()

TrackSurfaceIntersection Trk::FitParameters::intersection

(

void

)

const

Definition at line 152 of file FitParameters.cxx.

                                                               {
  return TrackSurfaceIntersection(
      m_position,
      Amg::Vector3D(m_cosPhi * m_sinTheta, m_sinPhi * m_sinTheta, m_cosTheta),
      0.);
}

numberAlignments() [1/2]

void Trk::FitParameters::numberAlignments

(

int

numberAlignments

)

Definition at line 159 of file FitParameters.cxx.

                                                         {
  m_numberAlignments = 0;
  if (!numberAlignments)
    return;
  m_alignmentAngle = std::vector<double>(numberAlignments, 0.);
  m_alignmentAngleConstraint = std::vector<double>(numberAlignments, 0.);
  m_alignmentOffset = std::vector<double>(numberAlignments, 0.);
  m_alignmentOffsetConstraint = std::vector<double>(numberAlignments, 0.);
}

numberAlignments() [2/2]

int Trk::FitParameters::numberAlignments ( void ) const

inline

Definition at line 233 of file FitParameters.h.

                                                     {
  return m_numberAlignments;
}

numberOscillations()

int Trk::FitParameters::numberOscillations ( void ) const

inline

Definition at line 237 of file FitParameters.h.

                                                       {
  return m_numberOscillations;
}

numberParameters() [1/2]

void Trk::FitParameters::numberParameters

(

int

numberParameters

)

Definition at line 169 of file FitParameters.cxx.

                                                         {
  m_numberParameters = numberParameters;
  m_firstAlignmentParameter = numberParameters - 2 * m_numberAlignments - 1;
  m_firstScatteringParameter =
      m_firstAlignmentParameter - 2 * m_numberScatterers;
}

numberParameters() [2/2]

int Trk::FitParameters::numberParameters ( void ) const

inline

Definition at line 241 of file FitParameters.h.

                                                     {
  return m_numberParameters;
}

numberScatterers() [1/2]

void Trk::FitParameters::numberScatterers

(

int

numberScatterers

)

Definition at line 176 of file FitParameters.cxx.

                                                         {
  m_numberScatterers = 0;
  if (!numberScatterers)
    return;
  m_scattererPhi = std::vector<double>(numberScatterers, 0.);
  m_scattererTheta = std::vector<double>(numberScatterers, 0.);
}

numberScatterers() [2/2]

int Trk::FitParameters::numberScatterers ( void ) const

inline

Definition at line 245 of file FitParameters.h.

                                                     {
  return m_numberScatterers;
}

operator=() [1/2]

FitParameters & Trk::FitParameters::operator= ( const FitParameters & )

default

operator=() [2/2]

FitParameters & Trk::FitParameters::operator= ( FitParameters && )

default

parameterDifference()

const Amg::MatrixX Trk::FitParameters::parameterDifference

(

const Amg::VectorX &

parameters

)

const

Definition at line 184 of file FitParameters.cxx.

                                        {
  Amg::MatrixX difference(1, 5);
  difference(0, 0) = parameters(0) - m_d0;
  difference(0, 1) = parameters(1) - m_z0;
  difference(0, 2) = parameters(3) * m_cosPhi - parameters(2) * m_sinPhi;
 
  // FIXME: diff is now delta(theta) : check sign
  // difference[0][3] = parameters[4] - m_cotTheta;
  difference(0, 3) = std::sin(parameters(4)) * m_cosTheta -
                     std::cos(parameters(4)) * m_sinTheta;
  difference(0, 4) = (parameters(5) - m_qOverP) * Gaudi::Units::TeV;
  return difference;
}

performCutStep()

void Trk::FitParameters::performCutStep

(

double

cutStep

)

Definition at line 199 of file FitParameters.cxx.

                                                 {
  // revert parameters to previous parameter change with cutStep*value
  // i.e. 0 < cutstep < 1 such that cutStep = 0 gives complete reversion
 
  Amg::VectorX cutDifferences(m_differences);
  cutDifferences *= (cutStep - 1.);
  Amg::VectorX oldDifferences(m_differences);
  oldDifferences *= cutStep;
 
  // leave phi alone when unstable
  if (m_phiInstability) {
    cutDifferences(2) = 0.;
    oldDifferences(2) = (m_differences)(2);
  }
 
  // apply cut
  update(cutDifferences);
  m_differences = Amg::VectorX(oldDifferences);
 
  m_numberOscillations = 0;
  m_oldDifference = 0.;
  // std::cout << " after cutstep " << std::endl;
}

perigee()

Perigee * Trk::FitParameters::perigee

(

void

)

const

Definition at line 223 of file FitParameters.cxx.

                                          {
  // copy 'final' covariance
  AmgSymMatrix(5) covMatrix = AmgSymMatrix(5)(*m_finalCovariance);
  const double pT = std::abs(m_sinTheta / m_qOverP);
  double charge = 1.;
  if (m_qOverP < 0.)
    charge = -1.;
  const Amg::Vector3D momentum(pT * m_cosPhi, pT * m_sinPhi, pT * m_cotTheta);
 
  if (m_surface) {
    return new Perigee(m_position, momentum, charge,
                       dynamic_cast<const Trk::PerigeeSurface&>(*m_surface),
                       std::move(covMatrix));
  } else {
    return new Perigee(m_position, momentum, charge, m_vertex,
                       std::move(covMatrix));
  }
}

phiInstability()

bool Trk::FitParameters::phiInstability

(

void

)

const

Definition at line 242 of file FitParameters.cxx.

                                             {
  return m_phiInstability;
}

position()

const Amg::Vector3D & Trk::FitParameters::position ( void ) const

inline

Definition at line 253 of file FitParameters.h.

                                                            {
  return m_position;
}

print()

void Trk::FitParameters::print

(

MsgStream &

log

)

const

Definition at line 246 of file FitParameters.cxx.

                                              {
  log << std::setiosflags(std::ios::fixed | std::ios::right) << std::setw(16)
      << std::setprecision(1) << m_position.perp() << " perigee radius"
      << std::setw(10) << std::setprecision(3) << m_d0 << " d0" << std::setw(11)
      << std::setprecision(3) << m_position.z() << " z0" << std::setw(11)
      << std::setprecision(6) << std::atan2(m_sinPhi, m_cosPhi) << " phi"
      << std::setw(11) << std::setprecision(6) << std::acos(m_cosTheta)
      << " theta" << std::setw(13) << std::setprecision(3)
      << m_sinTheta / (m_qOverP * Gaudi::Units::GeV) << " pT (GeV)";
}

printCovariance()

void Trk::FitParameters::printCovariance

(

MsgStream &

log

)

const

Definition at line 257 of file FitParameters.cxx.

                                                        {
  double error00 = 0.;
  const Amg::MatrixX& cov = *m_finalCovariance;
  if (cov(0, 0) > 0.)
    error00 = std::sqrt(cov(0, 0));
  double error11 = 0.;
  if (cov(1, 1) > 0.)
    error11 = std::sqrt(cov(1, 1));
  double error22 = 0.;
  if (cov(2, 2) > 0.)
    error22 = std::sqrt(cov(2, 2));
  double error33 = 0.;
  if (cov(3, 3) > 0.)
    error33 = std::sqrt(cov(3, 3));
  double error44 = 0.;
  if (cov(4, 4) > 0.)
    error44 = std::sqrt(cov(4, 4));
  double correl02 = 0.;
  const double denom02 = cov(0, 0) * cov(2, 2);
  if (denom02 > 0.)
    correl02 = cov(0, 2) / std::sqrt(denom02);
  double correl13 = 0.;
  const double denom13 = cov(1, 1) * cov(3, 3);
  if (denom13 > 0.)
    correl13 = cov(1, 3) / std::sqrt(denom13);
  log << std::setiosflags(std::ios::fixed | std::ios::right) << std::setw(10)
      << std::setprecision(3) << error00 << std::setw(14)
      << std::setprecision(3) << error11 << std::setw(14)
      << std::setprecision(6) << error22 << std::setw(15)
      << std::setprecision(6) << error33 << std::setw(19)
      << std::setprecision(3)
      << (error44 * m_sinTheta) / (m_qOverP * m_qOverP * Gaudi::Units::GeV)
      << " (covariance)" << std::setw(9) << std::setprecision(3) << correl02
      << " Cd0phi" << std::setw(8) << std::setprecision(3) << correl13
      << " Cz0theta" << std::endl;
}

printVerbose()

void Trk::FitParameters::printVerbose

(

MsgStream &

log

)

const

Definition at line 294 of file FitParameters.cxx.

                                                     {
  log << std::endl;
 
  if (m_differences.size() != 0) {
    const Amg::VectorX& differences = m_differences;
    log << "      dParams ====" << std::setiosflags(std::ios::fixed)
        << std::setw(10) << std::setprecision(4) << differences(0) << " (0) "
        << std::setw(10) << std::setprecision(4) << differences(1) << " (1) "
        << std::setw(10) << std::setprecision(5) << differences(2) << " (2) "
        << std::setw(10) << std::setprecision(5) << differences(3) << " (3) "
        << std::setw(13) << std::setprecision(9)
        << differences(4) * Gaudi::Units::GeV / Gaudi::Units::TeV << " (4) ";
    if (m_fitEnergyDeposit)
      log << std::setiosflags(std::ios::fixed) << std::setw(13)
          << std::setprecision(9)
          << differences(5) * Gaudi::Units::GeV / Gaudi::Units::TeV << " (5) ";
    log << std::endl;
 
    if (m_numberAlignments) {
      log << "       dAlign ==== ";
      unsigned param = m_firstAlignmentParameter;
      for (int scat = 0; scat < m_numberAlignments; ++scat) {
        ++param;
        if (scat % 5 == 0 && scat > 0)
          log << std::endl << "                   ";
        log << std::setiosflags(std::ios::fixed) << std::setw(10)
            << std::setprecision(6) << differences(param);
        ++param;
        log << std::setiosflags(std::ios::fixed) << std::setw(10)
            << std::setprecision(6) << differences(param) << " ("
            << std::setw(2) << scat << "A)     ";
      }
      log << std::endl;
    }
 
    if (m_numberScatterers) {
      log << "        dScat ==== ";
      unsigned param = m_firstScatteringParameter;
      for (int scat = 0; scat < m_numberScatterers; ++scat) {
        ++param;
        if (scat % 5 == 0 && scat > 0)
          log << std::endl << "                   ";
        log << std::setiosflags(std::ios::fixed) << std::setw(10)
            << std::setprecision(6) << differences(param);
        ++param;
        log << std::setiosflags(std::ios::fixed) << std::setw(10)
            << std::setprecision(6) << differences(param) << " ("
            << std::setw(2) << scat << "S)     ";
      }
      log << std::endl;
    }
  }
 
  log << std::setiosflags(std::ios::fixed | std::ios::right)
      << "   parameters:   " << std::setw(12) << std::setprecision(4) << m_d0
      << " transverse impact  " << std::setw(10) << std::setprecision(4)
      << m_position.z() << " z0  " << std::setw(10) << std::setprecision(6)
      << std::atan2(m_sinPhi, m_cosPhi) << std::setw(10) << std::setprecision(6)
      << m_cotTheta << " phi,cotTheta  " << std::setw(13)
      << std::setprecision(9) << m_qOverP / m_sinTheta << " inverse pT  "
      << std::setw(12) << std::setprecision(6)
      << m_sinTheta / (m_qOverP * Gaudi::Units::GeV) << " signed pT  ";
  if (m_fitEnergyDeposit) {
    // TODO: should give fitted energy loss
    log << std::endl
        << "    E before/after energy deposit" << std::setw(12)
        << std::setprecision(3) << 1. / std::abs(m_qOverP * Gaudi::Units::GeV)
        << std::setw(12) << std::setprecision(3)
        << 1. / std::abs(m_qOverP1 * Gaudi::Units::GeV);
  }
  if (m_numberAlignments) {
    log << std::endl << "   alignment number, angle, offset:  ";
    for (int align = 0; align < m_numberAlignments; ++align) {
      log << std::setiosflags(std::ios::fixed) << std::setw(6) << align
          << std::setw(10) << std::setprecision(6) << m_alignmentAngle[align]
          << std::setw(10) << std::setprecision(6) << m_alignmentOffset[align];
      if ((align + 1) % 5 == 0)
        log << std::endl << "                                                ";
    }
  }
  if (m_numberScatterers) {
    log << std::endl << "   scatterer number, delta(phi), delta(theta):  ";
    for (int scat = 0; scat < m_numberScatterers; ++scat) {
      log << std::setiosflags(std::ios::fixed) << std::setw(6) << scat
          << std::setw(10) << std::setprecision(6) << m_scattererPhi[scat]
          << std::setw(10) << std::setprecision(6) << m_scattererTheta[scat];
      if ((scat + 1) % 5 == 0)
        log << std::endl << "                                                ";
    }
  }
  log << std::endl;
}

ptInv0()

double Trk::FitParameters::ptInv0 ( void ) const

inline

Definition at line 249 of file FitParameters.h.

                                              {
  return m_qOverP / m_sinTheta;
}

qOverP() [1/2]

void Trk::FitParameters::qOverP

(

double

value

)

Definition at line 387 of file FitParameters.cxx.

                                       {
  m_qOverP = value;
}

qOverP() [2/2]

double Trk::FitParameters::qOverP ( void ) const

inline

Definition at line 257 of file FitParameters.h.

                                              {
  return m_qOverP;
}

qOverP1() [1/2]

void Trk::FitParameters::qOverP1

(

double

value

)

Definition at line 391 of file FitParameters.cxx.

                                        {
  m_qOverP1 = value;
}

qOverP1() [2/2]

double Trk::FitParameters::qOverP1 ( void ) const

inline

Definition at line 261 of file FitParameters.h.

                                               {
  return m_qOverP1;
}

reset()

void Trk::FitParameters::reset

(

const FitParameters &

parameters

)

Definition at line 395 of file FitParameters.cxx.

                                                         {
  // method is needed to complement copy in places where design uses
  // a reference to a FitParameter pointer
  // essentially a copy, with history of previous iteration removed
  m_cosPhi = parameters.m_cosPhi;
  m_cosPhi1 = parameters.m_cosPhi1;
  m_cosTheta = parameters.m_cosTheta;
  m_cosTheta1 = parameters.m_cosTheta1;
  m_cotTheta = parameters.m_cotTheta;
  m_d0 = parameters.m_d0;
  m_extremeMomentum = parameters.m_extremeMomentum;
  m_finalCovariance = parameters.m_finalCovariance;
  m_firstScatteringParameter = parameters.m_firstScatteringParameter;
  m_firstScatteringParameter = parameters.m_firstScatteringParameter;
  m_fitEnergyDeposit = parameters.m_fitEnergyDeposit;
  m_fitMomentum = parameters.m_fitMomentum;
  m_fullCovariance = parameters.m_fullCovariance;
  m_minEnergyDeposit = parameters.m_minEnergyDeposit;
  m_numberAlignments = parameters.m_numberAlignments;
  m_numberOscillations = parameters.m_numberOscillations;
  m_numberParameters = parameters.m_numberParameters;
  m_numberScatterers = parameters.m_numberScatterers;
  m_oldDifference = parameters.m_oldDifference;
  m_perigee = parameters.m_perigee;
  m_phiInstability = parameters.m_phiInstability;
  m_position = parameters.m_position;
  m_qOverP = parameters.m_qOverP;
  m_qOverP1 = parameters.m_qOverP1;
  m_sinPhi = parameters.m_sinPhi;
  m_sinPhi1 = parameters.m_sinPhi1;
  m_sinTheta = parameters.m_sinTheta;
  m_sinTheta1 = parameters.m_sinTheta1;
  m_vertex = parameters.vertex();
  m_z0 = parameters.m_z0;
  for (int s = 0; s != m_numberAlignments; ++s) {
    m_alignmentAngle[s] = parameters.m_alignmentAngle[s];
    m_alignmentAngleConstraint[s] = parameters.m_alignmentAngleConstraint[s];
    m_alignmentOffset[s] = parameters.m_alignmentOffset[s];
    m_alignmentOffsetConstraint[s] = parameters.m_alignmentOffsetConstraint[s];
  }
  for (int s = 0; s != m_numberScatterers; ++s) {
    m_scattererPhi[s] = parameters.m_scattererPhi[s];
    m_scattererTheta[s] = parameters.m_scattererTheta[s];
  }
 
  // restore difference history
  if (parameters.m_differences.size() != 0) {
    m_differences = Amg::VectorX(parameters.m_differences);
  } else {
    m_differences = Amg::VectorX(m_numberParameters);
    m_differences.setZero();
  }
 
  m_numberOscillations = 0;
  m_oldDifference = 0.;
}

resetOscillations()

void Trk::FitParameters::resetOscillations ( void )

inline

Definition at line 281 of file FitParameters.h.

                                                 {
  m_numberOscillations = 0;
}

scattererPhi()

double Trk::FitParameters::scattererPhi ( int scatterer ) const

inline

Definition at line 265 of file FitParameters.h.

                                                             {
  return m_scattererPhi[scatterer];
}

scattererTheta()

double Trk::FitParameters::scattererTheta ( int scatterer ) const

inline

Definition at line 269 of file FitParameters.h.

                                                               {
  return m_scattererTheta[scatterer];
}

scatteringAngles()

ScatteringAngles Trk::FitParameters::scatteringAngles	(	const FitMeasurement &	fitMeasurement,
		int	scatterer = -1 ) const

Definition at line 452 of file FitParameters.cxx.

                                                               {
  // scattering sigma used in chi2 computation
  const double scattererSigmaTheta = 1. / fitMeasurement.weight();
  const double scattererSigmaPhi =
      scattererSigmaTheta /
      fitMeasurement.intersection(FittedTrajectory).direction().perp();
  if (scatterer < 0) {
    return {0., 0., scattererSigmaPhi, scattererSigmaTheta};
  } else {
    return {m_scattererPhi[scatterer], m_scattererTheta[scatterer],
            scattererSigmaPhi, scattererSigmaTheta};
  }
}

setPhiInstability()

void Trk::FitParameters::setPhiInstability

(

void

)

Definition at line 467 of file FitParameters.cxx.

                                          {
  m_phiInstability = true;
}

sinPhi()

double Trk::FitParameters::sinPhi ( void ) const

inline

Definition at line 273 of file FitParameters.h.

                                              {
  return m_sinPhi;
}

sinTheta()

double Trk::FitParameters::sinTheta ( void ) const

inline

Definition at line 277 of file FitParameters.h.

                                                {
  return m_sinTheta;
}

startingPerigee()

Perigee * Trk::FitParameters::startingPerigee

(

void

)

const

Definition at line 471 of file FitParameters.cxx.

                                                  {
  // create momentum
  const double pT = std::abs(m_sinTheta / m_qOverP);
  double charge = 1.;
  if (m_qOverP < 0.)
    charge = -1.;
  const Amg::Vector3D momentum(pT * m_cosPhi, pT * m_sinPhi, pT * m_cotTheta);
 
  return new Perigee(m_position, momentum, charge, m_vertex);
}

trackParameters()

TrackParameters * Trk::FitParameters::trackParameters	(	MsgStream &	log,
		const FitMeasurement &	measurement,
		bool	withCovariance = false )

Definition at line 482 of file FitParameters.cxx.

                                                                            {
  // make checks necessary for the TrackParameters to be computed
  //   1) a Surface is required
  if (!measurement.surface()) {
    log << MSG::WARNING
        << "FitParameters::trackParameters - measurement lacks Surface"
        << endmsg;
    return nullptr;
  }
 
  //   2) a SurfaceIntersection is required
  if (!measurement.hasIntersection(FittedTrajectory)) {
    log << MSG::WARNING
        << "FitParameters::trackParameters - invalid measurement" << endmsg;
    return nullptr;
  }
 
  //   3) the intersection position has to lie sufficiently close to the Surface
  const TrackSurfaceIntersection& intersection =
      measurement.intersection(FittedTrajectory);
  Amg::Vector2D localPos;
  if (!measurement.surface()->globalToLocal(
          intersection.position(), intersection.direction(), localPos)) {
    log << MSG::WARNING
        << "FitParameters::trackParameters - globalToLocal failure" << endmsg;
    return nullptr;
  }
 
  // cache parameters at EnergyDeposit
  if (measurement.isEnergyDeposit()) {
    m_sinTheta1 = intersection.direction().perp();
    m_cosPhi1 = intersection.direction().x() / m_sinTheta1;
    m_sinPhi1 = intersection.direction().y() / m_sinTheta1;
    m_cosTheta1 = intersection.direction().z();
    m_qOverP1 = measurement.qOverP();
  }
 
  // propagate full covariance to form localCovariance
  std::optional<AmgSymMatrix(5)> covMatrix = std::nullopt;
  if (withCovariance && (measurement.isDrift() || measurement.isCluster() ||
                         measurement.isPerigee())) {
    Amg::Vector3D direction = intersection.direction();
    const double sigma = 1. / measurement.weight();
    const double sigma2 = 1. / measurement.weight2();
    int const lastParameter = measurement.lastParameter();
    Amg::MatrixX jacobian = Amg::MatrixX::Zero(5, lastParameter);
    for (int i = 0; i != lastParameter; ++i) {
      jacobian(0, i) = sigma * measurement.derivative(i);
      jacobian(1, i) = sigma2 * measurement.derivative2(i);
    }
 
    // phi,theta derivs into jac
    jacobian(2, 2) = 1.;
    jacobian(3, 3) = 1.;
    int i = m_firstScatteringParameter;
    while (++i < lastParameter) {
      jacobian(2, i) = 1.;
      ++i;
      jacobian(3, i) = 1.;
    }
 
    // only if fit to curvature
    if (m_fitMomentum && m_qOverP) {
      const double sinTheta = direction.perp();
      if (m_fitEnergyDeposit && measurement.afterCalo()) {
        const double deltaPhi =
            (direction.y() * m_cosPhi1 - direction.x() * m_sinPhi1) / sinTheta;
        const double deltaTheta =
            (sinTheta * m_cosTheta1 - direction.z() * m_sinTheta1);
        jacobian(0, 5) *= Gaudi::Units::TeV;
        jacobian(1, 5) *= Gaudi::Units::TeV;
        jacobian(2, 5) = deltaPhi / measurement.qOverP();
        jacobian(3, 5) = deltaTheta / measurement.qOverP();
        jacobian(4, 5) = measurement.qOverP() / m_qOverP1;
      } else {
        const double deltaPhi =
            (direction.y() * m_cosPhi - direction.x() * m_sinPhi) / sinTheta;
        const double deltaTheta =
            (sinTheta * m_cosTheta - direction.z() * m_sinTheta);
        jacobian(0, 4) *= Gaudi::Units::TeV;
        jacobian(1, 4) *= Gaudi::Units::TeV;
        jacobian(2, 4) = deltaPhi / measurement.qOverP();
        jacobian(3, 4) = deltaTheta / measurement.qOverP();
        jacobian(4, 4) = measurement.qOverP() / m_qOverP;
      }
    } else {
      jacobian(4, 4) = 1.;
    }
 
    // similarity transform
    covMatrix = AmgSymMatrix(5)(
        jacobian * m_fullCovariance->block(0, 0, lastParameter, lastParameter) *
        jacobian.transpose());
  }
 
  double phi = intersection.direction().phi();
  double theta = intersection.direction().theta();
  if (measurement.isFlipped()) {
    if (phi > 0.) {
      phi -= M_PI;
    } else {
      phi += M_PI;
    }
    theta = M_PI - theta;
  }
 
  if (!measurement.surface()) {
    log << MSG::WARNING
        << "FitParameters::trackParameters - unrecognized surface" << endmsg;
    return nullptr;
  }
  // finally can create the appropriate TrackParameters based on the surface
  return measurement.surface()
      ->createUniqueTrackParameters(localPos[locR], localPos[locZ], phi, theta,
                                    measurement.qOverP(), std::move(covMatrix))
      .release();
}

update() [1/2]

void Trk::FitParameters::update	(	Amg::Vector3D	position,
		Amg::Vector3D	direction,
		double	qOverP,
		const Amg::MatrixX &	leadingCovariance )

Definition at line 670 of file FitParameters.cxx.

                                                                {
  // update parameters after leading material corrections
  m_position = position;
  m_sinTheta = direction.perp();
  const double sinThetaInv = 1. / m_sinTheta;
  m_cotTheta = sinThetaInv * m_cosTheta;
  m_cosPhi = sinThetaInv * direction.x();
  m_sinPhi = sinThetaInv * direction.y();
  m_cotTheta = sinThetaInv * direction.z();
  m_cosTheta = direction.z();
  m_qOverP = qOverP;
  m_z0 = position.z();
  m_d0 = (m_vertex.x() - position.x()) * m_sinPhi -
         (m_vertex.y() - position.y()) * m_cosPhi;
 
  // update covariance
  // (*m_finalCovariance) += leadingCovariance; // ??
  for (int i = 0; i != 5; ++i) {
    for (int j = 0; j != 5; ++j) {
      (*m_finalCovariance)(i, j) += leadingCovariance(i, j);
    }
  }
}

update() [2/2]

void Trk::FitParameters::update

(

const Amg::VectorX &

differences

)

Definition at line 601 of file FitParameters.cxx.

                                                        {
  // keep update values in case of cutStep procedure
  if (m_numberOscillations && m_oldDifference * differences(4) < 0.) {
    ++m_numberOscillations;
  } else {
    m_numberOscillations = 1;
  }
  m_differences = Amg::VectorX(differences);
  m_oldDifference = differences(4);
 
  // misalignment parameters
  std::vector<double>::iterator a = m_alignmentAngle.begin();
  std::vector<double>::iterator o = m_alignmentOffset.begin();
  int align = m_firstAlignmentParameter;
  for (int i = 0; i != m_numberAlignments; ++i) {
    (*a++) += differences(++align);
    (*o++) += differences(++align);
  }
 
  // scattering angles
  std::vector<double>::iterator p = m_scattererPhi.begin();
  std::vector<double>::iterator t = m_scattererTheta.begin();
  int scat = m_firstScatteringParameter;
  for (int i = 0; i != m_numberScatterers; ++i) {
    (*p++) += differences(++scat);
    (*t++) += differences(++scat);
  }
 
  // qOverP, cotTheta
  if (m_fitMomentum)
    m_qOverP += differences(4) / Gaudi::Units::TeV;
  m_cotTheta -= differences(3) / (m_sinTheta * m_sinTheta);
 
  // impose charge conservation and decreasing energy
  if (m_fitEnergyDeposit) {
    m_qOverP1 += differences(5) / Gaudi::Units::TeV;
    const double deposit = 1. / std::abs(m_qOverP) - 1. / std::abs(m_qOverP1);
    if (std::abs(deposit) < std::abs(m_minEnergyDeposit) ||
        deposit * m_minEnergyDeposit < 0. || m_qOverP * m_qOverP1 < 0.) {
      m_qOverP = 1. / (1. / std::abs(m_qOverP1) + m_minEnergyDeposit);
      if (m_qOverP1 < 0.)
        m_qOverP = -m_qOverP;
    }
  }
 
  // protect phi against some rounding instabilities
  double sinDPhi = differences(2);
  double cosDPhi = 0.;
  if (std::abs(sinDPhi) < 1.0) {
    cosDPhi = std::sqrt(1. - sinDPhi * sinDPhi);
  } else {
    if (sinDPhi > 0.) {
      sinDPhi = 1.0;
    } else {
      sinDPhi = -1.0;
    }
  }
 
  const double cosPhi = m_cosPhi * cosDPhi - m_sinPhi * sinDPhi;
  m_sinPhi = m_sinPhi * cosDPhi + m_cosPhi * sinDPhi;
  m_cosPhi = cosPhi;
  m_z0 += differences(1);
  m_d0 += differences(0);
  m_sinTheta = 1. / std::sqrt(1. + m_cotTheta * m_cotTheta);
  m_cosTheta = m_cotTheta * m_sinTheta;
  m_position = Amg::Vector3D(m_vertex.x() - m_d0 * m_sinPhi,
                             m_vertex.y() + m_d0 * m_cosPhi, m_z0);
}

vertex()

const Amg::Vector3D & Trk::FitParameters::vertex ( void ) const

inline

Definition at line 285 of file FitParameters.h.

                                                          {
  return m_vertex;
}

z0()

double Trk::FitParameters::z0 ( void ) const

inline

Definition at line 289 of file FitParameters.h.

                                          {
  return m_position.z();
}

Member Data Documentation

m_alignmentAngle

std::vector<double> Trk::FitParameters::m_alignmentAngle

private

Definition at line 113 of file FitParameters.h.

m_alignmentAngleConstraint

std::vector<double> Trk::FitParameters::m_alignmentAngleConstraint

private

Definition at line 114 of file FitParameters.h.

m_alignmentOffset

std::vector<double> Trk::FitParameters::m_alignmentOffset

private

Definition at line 115 of file FitParameters.h.

m_alignmentOffsetConstraint

std::vector<double> Trk::FitParameters::m_alignmentOffsetConstraint

private

Definition at line 116 of file FitParameters.h.

m_cosPhi

double Trk::FitParameters::m_cosPhi

private

Definition at line 117 of file FitParameters.h.

m_cosPhi1

double Trk::FitParameters::m_cosPhi1

private

Definition at line 118 of file FitParameters.h.

m_cosTheta

double Trk::FitParameters::m_cosTheta

private

Definition at line 119 of file FitParameters.h.

m_cosTheta1

double Trk::FitParameters::m_cosTheta1

private

Definition at line 120 of file FitParameters.h.

m_cotTheta

double Trk::FitParameters::m_cotTheta

private

Definition at line 121 of file FitParameters.h.

m_d0

double Trk::FitParameters::m_d0

private

Definition at line 122 of file FitParameters.h.

m_differences

Amg::VectorX Trk::FitParameters::m_differences {}

private

Definition at line 123 of file FitParameters.h.

{};

m_eigen

bool Trk::FitParameters::m_eigen {}

private

Definition at line 124 of file FitParameters.h.

{};

m_extremeMomentum

bool Trk::FitParameters::m_extremeMomentum

private

Definition at line 125 of file FitParameters.h.

m_finalCovariance

Amg::MatrixX* Trk::FitParameters::m_finalCovariance

private

Definition at line 126 of file FitParameters.h.

m_firstAlignmentParameter

int Trk::FitParameters::m_firstAlignmentParameter

private

Definition at line 127 of file FitParameters.h.

m_firstScatteringParameter

int Trk::FitParameters::m_firstScatteringParameter

private

Definition at line 128 of file FitParameters.h.

m_fitEnergyDeposit

bool Trk::FitParameters::m_fitEnergyDeposit

private

Definition at line 129 of file FitParameters.h.

m_fitMomentum

bool Trk::FitParameters::m_fitMomentum

private

Definition at line 130 of file FitParameters.h.

m_fullCovariance

const Amg::MatrixX* Trk::FitParameters::m_fullCovariance

private

Definition at line 131 of file FitParameters.h.

m_minEnergyDeposit

double Trk::FitParameters::m_minEnergyDeposit

private

Definition at line 132 of file FitParameters.h.

m_numberAlignments

int Trk::FitParameters::m_numberAlignments

private

Definition at line 133 of file FitParameters.h.

m_numberOscillations

int Trk::FitParameters::m_numberOscillations

private

Definition at line 134 of file FitParameters.h.

m_numberParameters

int Trk::FitParameters::m_numberParameters

private

Definition at line 135 of file FitParameters.h.

m_numberScatterers

int Trk::FitParameters::m_numberScatterers

private

Definition at line 136 of file FitParameters.h.

m_oldDifference

double Trk::FitParameters::m_oldDifference

private

Definition at line 137 of file FitParameters.h.

m_perigee

const Perigee* Trk::FitParameters::m_perigee

private

Definition at line 138 of file FitParameters.h.

m_phiInstability

bool Trk::FitParameters::m_phiInstability

private

Definition at line 139 of file FitParameters.h.

m_position

Amg::Vector3D Trk::FitParameters::m_position

private

Definition at line 140 of file FitParameters.h.

m_qOverP

double Trk::FitParameters::m_qOverP

private

Definition at line 141 of file FitParameters.h.

m_qOverP1

double Trk::FitParameters::m_qOverP1

private

Definition at line 142 of file FitParameters.h.

m_scattererPhi

std::vector<double> Trk::FitParameters::m_scattererPhi

private

Definition at line 143 of file FitParameters.h.

m_scattererTheta

std::vector<double> Trk::FitParameters::m_scattererTheta

private

Definition at line 144 of file FitParameters.h.

m_sinPhi

double Trk::FitParameters::m_sinPhi

private

Definition at line 145 of file FitParameters.h.

m_sinPhi1

double Trk::FitParameters::m_sinPhi1

private

Definition at line 146 of file FitParameters.h.

m_sinTheta

double Trk::FitParameters::m_sinTheta

private

Definition at line 147 of file FitParameters.h.

m_sinTheta1

double Trk::FitParameters::m_sinTheta1

private

Definition at line 148 of file FitParameters.h.

m_surface

const Surface* Trk::FitParameters::m_surface

private

Definition at line 149 of file FitParameters.h.

m_vertex

Amg::Vector3D Trk::FitParameters::m_vertex

private

Definition at line 150 of file FitParameters.h.

m_z0

double Trk::FitParameters::m_z0

private

Definition at line 151 of file FitParameters.h.

---

The documentation for this class was generated from the following files:

• FitParameters.h
• FitParameters.cxx