Trk::FitMeasurement Class Reference

#include <FitMeasurement.h>

Collaboration diagram for Trk::FitMeasurement:

Public Member Functions
	FitMeasurement (int hitIndex, HitOnTrack *hitOnTrack, const MeasurementBase *measurementBase)
	FitMeasurement (const MaterialEffectsBase *materialEffects, double particleMass, const Amg::Vector3D &position, double qOverP=0., bool calo=false)
	FitMeasurement (double radiationThickness, double deltaE, double particleMass, const Amg::Vector3D &position, const Amg::Vector3D &direction, double qOverP, const Surface *surface=0)
	FitMeasurement (const AlignmentEffectsOnTrack *alignmentEffects, const Amg::Vector3D &direction, const Amg::Vector3D &position)
	FitMeasurement (const TrackSurfaceIntersection &intersection, double shift)
	FitMeasurement (const TrackStateOnSurface &TSOS)
	FitMeasurement (int hitIndex, HitOnTrack *hitOnTrack, const Amg::Vector3D &position, double sigma, double sigma2, double sinStereo, int status, const Surface *surface, MeasurementType type)
	FitMeasurement (int hitIndex, HitOnTrack *hitOnTrack, const Amg::Vector3D &position, double sigma, double signedDriftDistance, double sinStereo, const Surface *surface)
	FitMeasurement (const TrackParameters &perigee)
	FitMeasurement (double d0, const Amg::Vector3D &position, double sigma)
	~FitMeasurement (void)
bool	afterCalo (void) const
double	alignmentAngle (void) const
void	alignmentAngle (double value)
const AlignmentEffectsOnTrack *	alignmentEffects (void) const
double	alignmentOffset (void) const
void	alignmentOffset (double value)
unsigned	alignmentParameter (void) const
void	alignmentParameter (unsigned value)
unsigned	alignmentParameter2 (void) const
void	alignmentParameter2 (unsigned value)
double	derivative (int param) const
void	derivative (int param, double value)
void	derivative (double *pointer)
double	derivative2 (int param) const
void	derivative2 (int param, double value)
void	derivative2 (double *pointer)
int	derivativeRow (void) const
void	derivativeRow (int row)
double	d0 (void) const
double	energyLoss (void) const
double	energyLossSigma (void) const
unsigned	firstParameter (void) const
void	firstParameter (unsigned value)
void	flipDriftDirection (void)
int	hitIndex (void) const
const HitOnTrack *	hitOnTrack (void) const
HitOnTrack *	hitOnTrack (void)
bool	hasIntersection (ExtrapolationType type) const
const TrackSurfaceIntersection &	intersection (ExtrapolationType type) const
void	intersection (ExtrapolationType type, const std::optional< TrackSurfaceIntersection > &value)
bool	isAlignment (void) const
bool	isCluster (void) const
bool	isDrift (void) const
bool	isEnergyDeposit (void) const
bool	isFlipped (void) const
bool	isMaterialDelimiter (void) const
bool	isOutlier (void) const
bool	isPassive (void) const
bool	isPerigee (void) const
bool	isPositionMeasurement (void) const
bool	isPseudo (void) const
bool	isScatterer (void) const
bool	isTrapezoidCluster (void) const
bool	isVertex (void) const
bool	is2Dimensional (void) const
unsigned	lastParameter (void) const
void	lastParameter (unsigned value, bool afterCalo=false)
const MaterialEffectsBase *	materialEffects (void) const
const MeasurementBase *	measurementBase (void) const
double	minEnergyDeposit (void) const
const Amg::Vector3D &	minimizationDirection (void) const
void	minimizationDirection (const Amg::Vector3D &value)
const Amg::Vector3D &	normal (void) const
int	numberDoF (void) const
void	numberDoF (int value)
bool	numericalDerivative (void) const
const Amg::VectorX &	perigee (void) const
const Amg::MatrixX &	perigeeWeight (void) const
const Amg::Vector3D &	position (void) const
void	print (MsgStream &log) const
double	qOverP (void) const
void	qOverP (double value)
double	radiationThickness (void) const
double	residual (void) const
void	residual (double value)
void	residual (std::vector< double >::iterator pointer)
double	residual2 (void) const
void	residual2 (double value)
double	scattererPhi (void) const
void	scattererPhi (double value)
double	scattererTheta (void) const
void	scattererTheta (double value)
void	scatteringAngle (double angle, double totalRadiationThickness)
const Amg::Vector3D &	sensorDirection (void) const
void	setEnergyGain (void)
void	setMaterialEffectsOwner (void)
void	setNumericalDerivative (void)
void	setOutlier (void)
void	setSigma (void)
void	setSigmaMinus (void)
void	setSigmaPlus (void)
void	setSigmaSymmetric (void)
double	sigma (void) const
double	sigma2 (void) const
double	signedDriftDistance (void) const
int	status (void) const
const Surface *	surface (void) const
MeasurementType	type (void) const
void	unsetOutlier (void)
double	weight (void) const
double	weight2 (void) const

Static Public Member Functions
static void	printHeading (MsgStream &log)

Private Member Functions
	FitMeasurement (const FitMeasurement &)
FitMeasurement &	operator= (const FitMeasurement &)

Private Attributes
bool	m_afterCalo
const AlignmentEffectsOnTrack *	m_alignmentEffects
unsigned	m_alignmentParameter
unsigned	m_alignmentParameter2
double	m_betaSquared
double *	m_derivative
double *	m_derivative2
int	m_derivativeRow
double	m_d0
double	m_energyLoss
unsigned	m_firstParameter
bool	m_flippedDriftDistance
int	m_hitIndex
HitOnTrack *	m_hitOnTrack
std::array< std::optional< TrackSurfaceIntersection >, ExtrapolationTypes >	m_intersection
unsigned	m_lastParameter
const MaterialEffectsBase *	m_materialEffects
bool	m_materialEffectsOwner
const MeasurementBase *	m_measurementBase
double	m_minEnergyDeposit
Amg::Vector3D	m_minimizationDirection
Amg::Vector3D	m_normal
int	m_numberDoF
bool	m_numericalDerivative
bool	m_outlier
double	m_particleMassSquared
Amg::VectorX	m_perigee
Amg::MatrixX	m_perigeeWeight
Amg::Vector3D	m_position
double	m_qOverP
double	m_radiationThickness
std::vector< double >::iterator	m_residual
double	m_scaleFactor
double	m_scatterPhi
double	m_scatterTheta
double	m_scatteringAngle
double	m_scatteringAngleOffSet
double	m_secondResidual
Amg::Vector3D	m_sensorDirection
double	m_sigma
double	m_sigmaMinus
double	m_sigmaPlus
double	m_signedDriftDistance
int	m_status
const Surface *	m_surface
MeasurementType	m_type
double	m_weight
double	m_weight2

Detailed Description

Definition at line 40 of file FitMeasurement.h.

Constructor & Destructor Documentation

FitMeasurement() [1/11]

Trk::FitMeasurement::FitMeasurement	(	int	hitIndex,
		HitOnTrack *	hitOnTrack,
		const MeasurementBase *	measurementBase )

Definition at line 41 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(hitIndex),
      m_hitOnTrack(hitOnTrack),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(measurementBase),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(measurementBase->localCovariance().cols()),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(measurementBase->associatedSurface().center()),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(&measurementBase->associatedSurface()),
      m_type(stripCluster),
      m_weight(1.),
      m_weight2(1.) {
  double sigma = 0.;
  if (m_numberDoF > 0)
    sigma = Amg::error(measurementBase->localCovariance(), locX);
  double sigma2 = 0.;
  if (m_numberDoF > 1)
    sigma2 = Amg::error(measurementBase->localCovariance(), locY);
 
  // remaining data according to surface (normal, sensor, minimization
  // directions etc)
  if (dynamic_cast<const PlaneSurface*>(m_surface)) {
    m_normal = Amg::Vector3D(m_surface->transform().rotation().col(2));
    const Amg::Vector3D posptr =
        m_surface->localToGlobal(measurementBase->localParameters());
    m_position = posptr;
 
    // special case to get sensor for endcap trapezoids (discs represented as
    // planes) thus sensor direction rotates according to localX - rather than
    // being parallel and there is only 1 'real' degree of freedom for
    // covariance
    if (m_numberDoF == 2 &&
        dynamic_cast<const TrapezoidBounds*>(&m_surface->bounds())) {
      m_numberDoF = 1;
      m_type = trapezoidCluster;
      const double aa = measurementBase->localCovariance()(locX, locX);
      const double ab = measurementBase->localCovariance()(locX, locY);
      const double bb = measurementBase->localCovariance()(locY, locY);
      const double sum = aa + bb;
      const double diff = std::sqrt(sum * sum - 4. * (aa * bb - ab * ab));
      // used by obsolete scaling
      // double lengthSq    = 0.5*(sum + diff);
      const double widthSq = 0.5 * (sum - diff);
      sigma = std::sqrt(widthSq);
      const double term = 0.5 * (aa - bb) / diff;
      const double cosStereo = std::sqrt(0.5 - term);
      double sinStereo = 0.0;
      if (term > -0.5) {
        sinStereo = std::sqrt(0.5 + term);
        if (ab * m_normal.z() < 0.)
          sinStereo = -sinStereo;
      }
 
      const Amg::Vector3D& axis = m_surface->transform().rotation().col(1);
      m_sensorDirection =
          Amg::Vector3D(axis(0) * cosStereo + axis(1) * sinStereo,
                        axis(1) * cosStereo - axis(0) * sinStereo, axis(2));
    } else {
      m_sensorDirection =
          Amg::Vector3D(m_surface->transform().rotation().col(1));
    }
 
    m_minimizationDirection = Amg::Vector3D(m_sensorDirection.cross(m_normal));
    if (m_numberDoF == 2)
      m_type = pixelCluster;
  } else if (dynamic_cast<const StraightLineSurface*>(m_surface)) {
    // StraightLines can be driftCircles or pseudoMeasurements along the wire
    if (!measurementBase->localParameters().contains(locY)) {
      // driftCircles will eventually have sagged surfaces.
      // then may need something like:
      //    get z-along wire from (globPos-center).dot(sensor)
      //    then loc3D = (0,0,zalongwire)
      //    position = surface->transform() * loc3D
      // but probably just same as using centre with wire direction ...
      m_sensorDirection =
          Amg::Vector3D(m_surface->transform().rotation().col(2));
      m_signedDriftDistance = measurementBase->localParameters()[driftRadius];
      m_type = driftCircle;
    } else {
      // pseudomeasurement - minimize along wire direction
      m_minimizationDirection =
          Amg::Vector3D(m_surface->transform().rotation().col(2));
      const double mag = m_surface->center().mag();
      m_normal = mag > 1e-6 ? Amg::Vector3D(m_surface->center() / mag)
                            : Amg::Vector3D(m_surface->normal());
      m_position = measurementBase->globalPosition();
      m_sensorDirection =
          Amg::Vector3D(m_normal.cross(m_minimizationDirection));
      m_signedDriftDistance = 0.;
      m_type = pseudoMeasurement;
    }
  }
 
  const PerigeeSurface* perigee =
      dynamic_cast<const PerigeeSurface*>(m_surface);
  if (perigee) {
    m_position = m_surface->center();
    m_sensorDirection = Amg::Vector3D(0., 0., 1.);
    m_type = transverseVertex;
    if (m_numberDoF == 2)
      m_type = vertex;
  }
 
  // there are no measurements from Atlas detectors of cylinder or disc type
  //     const CylinderSurface* cylinder        = dynamic_cast<const
  //     CylinderSurface*>(m_surface); if (cylinder)
  //     {
  //     }
 
  //     const DiscSurface* disc            = dynamic_cast<const
  //     DiscSurface*>(m_surface); if (disc)
  //     {
  //     }
 
  // add protection against junk input
  if (sigma > 0.)
    m_weight = 1. / sigma;
  if (sigma2 > 0.)
    m_weight2 = 1. / sigma2;
}

FitMeasurement() [2/11]

Trk::FitMeasurement::FitMeasurement	(	const MaterialEffectsBase *	materialEffects,
		double	particleMass,
		const Amg::Vector3D &	position,
		double	qOverP = 0.,
		bool	calo = false )

Definition at line 200 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(materialEffects),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(0),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(particleMass * particleMass),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(qOverP),
      m_radiationThickness(materialEffects->thicknessInX0()),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(&materialEffects->associatedSurface()),
      m_type(endcapScatterer),
      m_weight(0.),
      m_weight2(0.) {
  if (dynamic_cast<const CylinderSurface*>(m_surface) ||
      std::abs(m_surface->normal()(2)) < 0.5)
    m_type = barrelScatterer;
 
  if (calo)
    m_type = calorimeterScatterer;
 
  // set any energy loss
  const EnergyLoss* energyLoss = nullptr;
  const ScatteringAngles* scattering = nullptr;
  const MaterialEffectsOnTrack* meot =
      dynamic_cast<const MaterialEffectsOnTrack*>(materialEffects);
  if (meot) {
    energyLoss = meot->energyLoss();
    scattering = meot->scatteringAngles();
  }
 
  if (energyLoss) {
    // note: EDM defines energy loss to be negative
    //       but here take opposite as this convention is not accepted by
    //       CaloEnergy clients (exception to take verbatim for CaloEnergy)
    // m_energyLoss = -energyLoss->deltaE();
    m_energyLoss = std::abs(energyLoss->deltaE());
 
    // calo energy deposit treated as a single pure energy loss measurement,
    // with fit taking error into account
    if (calo && !scattering && energyLoss->sigmaDeltaE() > 0.) {
      m_energyLoss = energyLoss->deltaE();
      m_numberDoF = 1;
      m_sigma = energyLoss->sigmaDeltaE();
      m_sigmaMinus = energyLoss->sigmaMinusDeltaE();
      m_sigmaPlus = energyLoss->sigmaPlusDeltaE();
      m_type = energyDeposit;
      m_weight = 1. / m_sigma;
      m_minEnergyDeposit = 0.5 * m_energyLoss;
      if (m_energyLoss > 0. && m_energyLoss > 5.0 * m_sigmaMinus) {
        m_minEnergyDeposit = m_energyLoss - 2.5 * m_sigmaMinus;
      } else if (m_energyLoss < 0. && m_energyLoss < -5.0 * m_sigmaMinus) {
        m_minEnergyDeposit = m_energyLoss + 2.5 * m_sigmaMinus;
      }
    }
  }
 
  // initialize scattering angles
  if (scattering) {
    m_scatterPhi = scattering->deltaPhi();
    m_scatterTheta = scattering->deltaTheta();
    if (calo)
      m_scatteringAngleOffSet = scattering->sigmaDeltaTheta();
    //        if(calo) std::cout << " Calorimeter scaterrer with sigmaDeltaPhi "
    //        << scattering->sigmaDeltaPhi() << " sigmaDeltaTheta " <<
    //        scattering->sigmaDeltaTheta() << std::endl;
  }
}

FitMeasurement() [3/11]

Trk::FitMeasurement::FitMeasurement	(	double	radiationThickness,
		double	deltaE,
		double	particleMass,
		const Amg::Vector3D &	position,
		const Amg::Vector3D &	direction,
		double	qOverP,
		const Surface *	surface = 0 )

Definition at line 307 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(-deltaE),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(true),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(0),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(particleMass * particleMass),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(qOverP),
      m_radiationThickness(radiationThickness),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(surface),
      m_type(endcapScatterer),
      m_weight(0.),
      m_weight2(0.) {
  // plane surface with normal along input direction
  if (m_surface && (dynamic_cast<const CylinderSurface*>(m_surface) ||
                    std::abs(m_surface->normal()(2)) < 0.5))
    m_type = barrelScatterer;
  else if (std::abs(direction(2)) < 0.5)
    m_type = barrelScatterer;
  if (!m_surface) {
    const CurvilinearUVT uvt(direction);
    m_surface = new PlaneSurface(position, uvt);
  }
 
  // create MaterialEffects
  std::bitset<MaterialEffectsBase::NumberOfMaterialEffectsTypes> typeMaterial;
  if (deltaE != 0.)
    typeMaterial.set(MaterialEffectsBase::EnergyLossEffects);
  auto energyLoss = std::make_unique<EnergyLoss>(deltaE, 0., 0., 0.);
  m_materialEffects = new MaterialEffectsOnTrack(
      radiationThickness, std::move(energyLoss), *m_surface, typeMaterial);
  if (!surface)
    delete m_surface;
  m_surface = &m_materialEffects->associatedSurface();
 
  m_intersection[FittedTrajectory] =
      std::make_optional<TrackSurfaceIntersection>(position, direction, 0.);
}

FitMeasurement() [4/11]

Trk::FitMeasurement::FitMeasurement	(	const AlignmentEffectsOnTrack *	alignmentEffects,
		const Amg::Vector3D &	direction,
		const Amg::Vector3D &	position )

Definition at line 386 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(alignmentEffects),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(2),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(alignmentEffects->deltaAngle()),
      m_scatterTheta(alignmentEffects->deltaTranslation()),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(alignmentEffects->sigmaDeltaAngle()),
      m_sigmaPlus(alignmentEffects->sigmaDeltaTranslation()),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(&alignmentEffects->associatedSurface()),
      m_type(alignment),
      m_weight(1.),
      m_weight2(1.) {
  // set weights
  if (m_sigmaMinus)
    m_weight = 1. / m_sigmaMinus;
  if (m_sigmaPlus)
    m_weight2 = 1. / m_sigmaPlus;
 
  m_intersection[FittedTrajectory] =
      std::make_optional<TrackSurfaceIntersection>(position, direction, 0.);
}

FitMeasurement() [5/11]

Trk::FitMeasurement::FitMeasurement	(	const TrackSurfaceIntersection &	intersection,
		double	shift )

Definition at line 447 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(0),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(intersection.position()),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(nullptr),
      m_type(materialDelimiter),
      m_weight(0.),
      m_weight2(0.) {
  // plane surface with normal along input direction and shift wrt position
  const Amg::Vector3D offset = intersection.direction() * shift;
  m_position += offset;
  const CurvilinearUVT uvt(intersection.direction());
  m_surface = new PlaneSurface(m_position, uvt);
 
  m_intersection[FittedTrajectory] = std::make_optional<TrackSurfaceIntersection>(
      m_position, intersection.direction(), 0.);
}

FitMeasurement() [6/11]

Trk::FitMeasurement::FitMeasurement

(

const TrackStateOnSurface &

TSOS

)

Definition at line 507 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(0),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(TSOS.trackParameters()->position()),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(&TSOS.trackParameters()->associatedSurface()),
      m_type(hole),
      m_weight(0.),
      m_weight2(0.) {
  // set type if necessary (hole = default)
  if (TSOS.type(TrackStateOnSurface::Outlier)) {
    m_outlier = true;
  }
  //     else if (TSOS.type(TrackStateOnSurface::Perigee))
  //     {
  //    m_type = MSperigee;
  //     }
}

FitMeasurement() [7/11]

Trk::FitMeasurement::FitMeasurement	(	int	hitIndex,
		HitOnTrack *	hitOnTrack,
		const Amg::Vector3D &	position,
		double	sigma,
		double	sigma2,
		double	sinStereo,
		int	status,
		const Surface *	surface,
		MeasurementType	type )

Definition at line 566 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(hitIndex),
      m_hitOnTrack(hitOnTrack),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_numberDoF(1),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(status),
      m_surface(surface),
      m_type(type),
      m_weight(1.),
      m_weight2(1.) {
  // pixel has 2-D measurement
  if (type == pixelCluster) {
    m_numberDoF = 2;
  }
 
  // special treatment for projective trapezoidal chambers in the endcap
  // take sensorDirection from stereo angle as I don't understand Surface axes
  // in this case
  m_normal = Amg::Vector3D(m_surface->transform().rotation().col(2));
  if (m_numberDoF == 1 && std::abs(m_normal.z()) > 0.99 &&
      std::abs(sinStereo) < 0.5)  // end-cap projective geometry
  {
    const double cosStereo = std::sqrt(1. - sinStereo * sinStereo);
    m_sensorDirection =
        Amg::Vector3D(position(0) * cosStereo + position(1) * sinStereo,
                      -position(0) * sinStereo + position(1) * cosStereo, 0.);
    (m_sensorDirection) /= m_sensorDirection.perp();
  } else  // otherwise chambers have parallel strips with sensor direction =
          // appropriate module axis
  {
    m_sensorDirection = Amg::Vector3D(m_surface->transform().rotation().col(1));
  }
  m_minimizationDirection = Amg::Vector3D(m_sensorDirection.cross(m_normal));
 
  // add protection against junk input
  if (sigma > 0.) {
    m_weight = 1. / sigma;
  } else {
    m_numberDoF = 0;
    m_outlier = true;
  }
  if (m_numberDoF == 2) {
    if (sigma2 > 0.) {
      m_weight2 = 1. / sigma2;
    } else {
      m_numberDoF = 0;
      m_outlier = true;
    }
  }
}

FitMeasurement() [8/11]

Trk::FitMeasurement::FitMeasurement	(	int	hitIndex,
		HitOnTrack *	hitOnTrack,
		const Amg::Vector3D &	position,
		double	sigma,
		double	signedDriftDistance,
		double	sinStereo,
		const Surface *	surface )

Definition at line 656 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(hitIndex),
      m_hitOnTrack(hitOnTrack),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(1),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(signedDriftDistance),
      m_status(0),
      m_surface(surface),
      m_type(driftCircle),
      m_weight(1.),
      m_weight2(1.) {
  m_sensorDirection = Amg::Vector3D(m_surface->transform().rotation().col(2));
 
  // add protection against junk input
  if (sigma > 0.) {
    m_weight = 1. / sigma;
  } else {
    m_numberDoF = 0;
    m_outlier = true;
  }
}

FitMeasurement() [9/11]

Trk::FitMeasurement::FitMeasurement

(

const TrackParameters &

perigee

)

Definition at line 718 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(0.),
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(0),
      m_numericalDerivative(false),
      m_outlier(true),  // use base class for additional trackParameters at
                        // detector boundary
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(perigee.associatedSurface().center()),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection{},
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(&perigee.associatedSurface()),
      m_type(perigeeParameters),
      m_weight(1.),
      m_weight2(1.) {
  // perigee axis needed for propagation of fitted parameters
  m_sensorDirection = Amg::Vector3D(m_surface->transform().rotation().col(2));
 
  // is this perigee to be used as a measurement?
  if (perigee.covariance() && !m_outlier) {
    // use as measurement
    m_numberDoF = 5;
    m_outlier = false;
 
    // perigeeParameters as HepVector
    Amg::Vector3D momentum = perigee.momentum();
    double ptInv0 = 1. / momentum.perp();
    const double cosPhi = ptInv0 * momentum(0);
    const double sinPhi = ptInv0 * momentum(1);
    const double cotTheta = ptInv0 * momentum(2);
    ptInv0 *= perigee.charge();
 
    Amg::VectorX parameters(6);
    parameters(0) = perigee.parameters()[Trk::d0];
    parameters(1) = perigee.parameters()[Trk::z0];
    parameters(2) = cosPhi;
    parameters(3) = sinPhi;
    parameters(4) = cotTheta;
    parameters(5) = ptInv0;
    m_perigee = Amg::VectorX(parameters);
 
    // weight = inverse covariance
    AmgSymMatrix(5) covariance(*perigee.covariance());
 
    // convert to internal units (TeV) to avoid rounding
    for (int row = 0; row < 5; ++row) {
      covariance(4, row) *= Gaudi::Units::TeV;
      covariance(row, 4) = covariance(4, row);
    }
    covariance(4, 4) *= Gaudi::Units::TeV;
    covariance.inverse();
 
    JacobianCotThetaPtToThetaP jacobian(cotTheta, ptInv0);
    m_perigeeWeight =
        Amg::MatrixX(jacobian * covariance * jacobian.transpose());
  }
}

FitMeasurement() [10/11]

Trk::FitMeasurement::FitMeasurement	(	double	d0,
		const Amg::Vector3D &	position,
		double	sigma )

Definition at line 811 of file FitMeasurement.cxx.

    : m_afterCalo(false),
      m_alignmentEffects(nullptr),
      m_alignmentParameter(0),
      m_alignmentParameter2(0),
      m_betaSquared(1.),
      m_derivative(nullptr),
      m_derivative2(nullptr),
      m_derivativeRow(-1),
      m_d0(d0),  // FIXME:: kept for cache tag as d0 is never used anywhere
      m_energyLoss(0.),
      m_firstParameter(0),
      m_flippedDriftDistance(false),
      m_hitIndex(0),
      m_hitOnTrack(nullptr),
      m_lastParameter(0),
      m_materialEffects(nullptr),
      m_materialEffectsOwner(false),
      m_measurementBase(nullptr),
      m_minEnergyDeposit(0.),
      m_minimizationDirection{},
      m_normal{},
      m_numberDoF(1),
      m_numericalDerivative(false),
      m_outlier(false),
      m_particleMassSquared(0.),
      m_perigee{},
      m_perigeeWeight{},
      m_position(position),
      m_qOverP(0.),
      m_radiationThickness(0.),
      m_residual(0),
      m_scaleFactor(0.),
      m_scatterPhi(0.),
      m_scatterTheta(0.),
      m_scatteringAngle(0.),
      m_scatteringAngleOffSet(0.),
      m_secondResidual(0.),
      m_sensorDirection(Amg::Vector3D(0., 0., 1.)),
      m_sigma(0.),
      m_sigmaMinus(0.),
      m_sigmaPlus(0.),
      m_signedDriftDistance(0.),
      m_status(0),
      m_surface(new const Trk::PerigeeSurface(position)),
      m_type(transverseVertex),
      m_weight(1. / sigma),
      m_weight2(1.) {}

~FitMeasurement()

Trk::FitMeasurement::~FitMeasurement

(

void

)

Definition at line 862 of file FitMeasurement.cxx.

                                    {
  if ((m_type == materialDelimiter || m_type == transverseVertex) &&
      !m_measurementBase)
    delete m_surface;
  if (m_materialEffects && m_materialEffectsOwner)
    delete m_materialEffects;
}

FitMeasurement() [11/11]

Trk::FitMeasurement::FitMeasurement ( const FitMeasurement & )

private

Member Function Documentation

afterCalo()

bool Trk::FitMeasurement::afterCalo ( void ) const

inline

Definition at line 243 of file FitMeasurement.h.

                                                {
  return m_afterCalo;
}

alignmentAngle() [1/2]

void Trk::FitMeasurement::alignmentAngle ( double value )

inline

Definition at line 251 of file FitMeasurement.h.

                                                       {
  m_scatterPhi = value;
}

alignmentAngle() [2/2]

double Trk::FitMeasurement::alignmentAngle ( void ) const

inline

Definition at line 247 of file FitMeasurement.h.

                                                       {
  return m_scatterPhi;
}

alignmentEffects()

const AlignmentEffectsOnTrack * Trk::FitMeasurement::alignmentEffects ( void ) const

inline

Definition at line 255 of file FitMeasurement.h.

                {
  return m_alignmentEffects;
}

alignmentOffset() [1/2]

void Trk::FitMeasurement::alignmentOffset ( double value )

inline

Definition at line 264 of file FitMeasurement.h.

                                                        {
  m_scatterTheta = value;
}

alignmentOffset() [2/2]

double Trk::FitMeasurement::alignmentOffset ( void ) const

inline

Definition at line 260 of file FitMeasurement.h.

                                                        {
  return m_scatterTheta;
}

alignmentParameter() [1/2]

void Trk::FitMeasurement::alignmentParameter ( unsigned value )

inline

Definition at line 272 of file FitMeasurement.h.

                                                             {
  m_alignmentParameter = value;
}

alignmentParameter() [2/2]

unsigned Trk::FitMeasurement::alignmentParameter ( void ) const

inline

Definition at line 268 of file FitMeasurement.h.

                                                             {
  return m_alignmentParameter;
}

alignmentParameter2() [1/2]

void Trk::FitMeasurement::alignmentParameter2 ( unsigned value )

inline

Definition at line 280 of file FitMeasurement.h.

                                                              {
  m_alignmentParameter2 = value;
}

alignmentParameter2() [2/2]

unsigned Trk::FitMeasurement::alignmentParameter2 ( void ) const

inline

Definition at line 276 of file FitMeasurement.h.

                                                              {
  return m_alignmentParameter2;
}

d0()

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

inline

Definition at line 316 of file FitMeasurement.h.

                                           {
  return m_d0;
}

derivative() [1/3]

void Trk::FitMeasurement::derivative ( double * pointer )

inline

Definition at line 292 of file FitMeasurement.h.

                                                      {
  m_derivative = pointer;
}

derivative() [2/3]

double Trk::FitMeasurement::derivative ( int param ) const

inline

Definition at line 284 of file FitMeasurement.h.

                                                        {
  return *(m_derivative + param);
}

derivative() [3/3]

void Trk::FitMeasurement::derivative ( int param, double value )

inline

Definition at line 288 of file FitMeasurement.h.

                                                              {
  *(m_derivative + param) = value;
}

derivative2() [1/3]

void Trk::FitMeasurement::derivative2 ( double * pointer )

inline

Definition at line 304 of file FitMeasurement.h.

                                                       {
  m_derivative2 = pointer;
}

derivative2() [2/3]

double Trk::FitMeasurement::derivative2 ( int param ) const

inline

Definition at line 296 of file FitMeasurement.h.

                                                         {
  return *(m_derivative2 + param);
}

derivative2() [3/3]

void Trk::FitMeasurement::derivative2 ( int param, double value )

inline

Definition at line 300 of file FitMeasurement.h.

                                                               {
  *(m_derivative2 + param) = value;
}

derivativeRow() [1/2]

void Trk::FitMeasurement::derivativeRow ( int row )

inline

Definition at line 312 of file FitMeasurement.h.

                                                 {
  m_derivativeRow = row;
}

derivativeRow() [2/2]

int Trk::FitMeasurement::derivativeRow ( void ) const

inline

Definition at line 308 of file FitMeasurement.h.

                                                   {
  return m_derivativeRow;
}

energyLoss()

double Trk::FitMeasurement::energyLoss ( void ) const

inline

Definition at line 320 of file FitMeasurement.h.

                                                   {
  return m_energyLoss;
}

energyLossSigma()

double Trk::FitMeasurement::energyLossSigma ( void ) const

inline

Definition at line 324 of file FitMeasurement.h.

                                                        {
  return m_sigma;
}

firstParameter() [1/2]

void Trk::FitMeasurement::firstParameter ( unsigned value )

inline

Definition at line 332 of file FitMeasurement.h.

                                                         {
  m_firstParameter = value;
}

firstParameter() [2/2]

unsigned Trk::FitMeasurement::firstParameter ( void ) const

inline

Definition at line 328 of file FitMeasurement.h.

                                                         {
  return m_firstParameter;
}

flipDriftDirection()

void Trk::FitMeasurement::flipDriftDirection ( void )

inline

Definition at line 336 of file FitMeasurement.h.

                                                   {
  if (m_type == driftCircle) {
    m_flippedDriftDistance = true;
    m_signedDriftDistance = -m_signedDriftDistance;
  }
}

hasIntersection()

bool Trk::FitMeasurement::hasIntersection ( ExtrapolationType type ) const

inline

Definition at line 355 of file FitMeasurement.h.

                                                                        {
  return (bool)m_intersection[type];
}

hitIndex()

int Trk::FitMeasurement::hitIndex ( void ) const

inline

Definition at line 343 of file FitMeasurement.h.

                                              {
  return m_hitIndex;
}

hitOnTrack() [1/2]

HitOnTrack * Trk::FitMeasurement::hitOnTrack ( void )

inline

Definition at line 351 of file FitMeasurement.h.

                                                  {
  return m_hitOnTrack;
}

hitOnTrack() [2/2]

const HitOnTrack * Trk::FitMeasurement::hitOnTrack ( void ) const

inline

Definition at line 347 of file FitMeasurement.h.

                                                              {
  return m_hitOnTrack;
}

intersection() [1/2]

const TrackSurfaceIntersection & Trk::FitMeasurement::intersection ( ExtrapolationType type ) const

inline

Definition at line 359 of file FitMeasurement.h.

                                  {
  return *m_intersection[type];
}

intersection() [2/2]

void Trk::FitMeasurement::intersection	(	ExtrapolationType	type,
		const std::optional< TrackSurfaceIntersection > &	value )

Definition at line 870 of file FitMeasurement.cxx.

                                                                                      {
  m_intersection[type] = value;
}

is2Dimensional()

bool Trk::FitMeasurement::is2Dimensional ( void ) const

inline

Definition at line 423 of file FitMeasurement.h.

                                                     {
  return (m_numberDoF == 2);
}

isAlignment()

bool Trk::FitMeasurement::isAlignment ( void ) const

inline

Definition at line 364 of file FitMeasurement.h.

                                                  {
  return (m_type == alignment);
}

isCluster()

bool Trk::FitMeasurement::isCluster ( void ) const

inline

Definition at line 368 of file FitMeasurement.h.

                                                {
  return (m_type == pixelCluster || m_type == stripCluster ||
          m_type == trapezoidCluster);
}

isDrift()

bool Trk::FitMeasurement::isDrift ( void ) const

inline

Definition at line 373 of file FitMeasurement.h.

                                              {
  return (m_type == driftCircle);
}

isEnergyDeposit()

bool Trk::FitMeasurement::isEnergyDeposit ( void ) const

inline

Definition at line 377 of file FitMeasurement.h.

                                                      {
  return (m_type == energyDeposit);
}

isFlipped()

bool Trk::FitMeasurement::isFlipped ( void ) const

inline

Definition at line 381 of file FitMeasurement.h.

                                                {
  return (m_type == driftCircle && m_flippedDriftDistance);
}

isMaterialDelimiter()

bool Trk::FitMeasurement::isMaterialDelimiter ( void ) const

inline

Definition at line 385 of file FitMeasurement.h.

                                                          {
  return (m_type == materialDelimiter);
}

isOutlier()

bool Trk::FitMeasurement::isOutlier ( void ) const

inline

Definition at line 389 of file FitMeasurement.h.

                                                {
  return m_outlier;
}

isPassive()

bool Trk::FitMeasurement::isPassive ( void ) const

inline

Definition at line 393 of file FitMeasurement.h.

                                                {
  return (m_type > bremPoint);
}

isPerigee()

bool Trk::FitMeasurement::isPerigee ( void ) const

inline

Definition at line 397 of file FitMeasurement.h.

                                                {
  return (m_type == perigeeParameters);
}

isPositionMeasurement()

bool Trk::FitMeasurement::isPositionMeasurement ( void ) const

inline

Definition at line 401 of file FitMeasurement.h.

                                                            {
  return (m_type < barrelScatterer);
}

isPseudo()

bool Trk::FitMeasurement::isPseudo ( void ) const

inline

Definition at line 405 of file FitMeasurement.h.

                                               {
  return (m_type == pseudoMeasurement);
}

isScatterer()

bool Trk::FitMeasurement::isScatterer ( void ) const

inline

Definition at line 409 of file FitMeasurement.h.

                                                  {
  return (m_type == barrelScatterer || m_type == endcapScatterer ||
          m_type == calorimeterScatterer || m_type == barrelInert ||
          m_type == endcapInert);
}

isTrapezoidCluster()

bool Trk::FitMeasurement::isTrapezoidCluster ( void ) const

inline

Definition at line 415 of file FitMeasurement.h.

                                                         {
  return (m_type == trapezoidCluster);
}

isVertex()

bool Trk::FitMeasurement::isVertex ( void ) const

inline

Definition at line 419 of file FitMeasurement.h.

                                               {
  return (m_type == vertex || m_type == transverseVertex);
}

lastParameter() [1/2]

void Trk::FitMeasurement::lastParameter ( unsigned value, bool afterCalo = false )

inline

Definition at line 431 of file FitMeasurement.h.

                                                                        {
  m_afterCalo = afterCalo;
  m_lastParameter = value;
}

lastParameter() [2/2]

unsigned Trk::FitMeasurement::lastParameter ( void ) const

inline

Definition at line 427 of file FitMeasurement.h.

                                                        {
  return m_lastParameter;
}

materialEffects()

const MaterialEffectsBase * Trk::FitMeasurement::materialEffects ( void ) const

inline

Definition at line 436 of file FitMeasurement.h.

                                                                            {
  return m_materialEffects;
}

measurementBase()

const MeasurementBase * Trk::FitMeasurement::measurementBase ( void ) const

inline

Definition at line 440 of file FitMeasurement.h.

                                                                        {
  return m_measurementBase;
}

minEnergyDeposit()

double Trk::FitMeasurement::minEnergyDeposit ( void ) const

inline

Definition at line 444 of file FitMeasurement.h.

                                                         {
  return m_minEnergyDeposit;
}

minimizationDirection() [1/2]

void Trk::FitMeasurement::minimizationDirection ( const Amg::Vector3D & value )

inline

Definition at line 452 of file FitMeasurement.h.

                                                                          {
  m_minimizationDirection = value;
}

minimizationDirection() [2/2]

const Amg::Vector3D & Trk::FitMeasurement::minimizationDirection ( void ) const

inline

Definition at line 448 of file FitMeasurement.h.

                                                                          {
  return m_minimizationDirection;
}

normal()

const Amg::Vector3D & Trk::FitMeasurement::normal ( void ) const

inline

Definition at line 456 of file FitMeasurement.h.

                                                           {
  return m_normal;
}

numberDoF() [1/2]

void Trk::FitMeasurement::numberDoF ( int value )

inline

Definition at line 466 of file FitMeasurement.h.

                                               {
  m_numberDoF = value;
}

numberDoF() [2/2]

int Trk::FitMeasurement::numberDoF ( void ) const

inline

Definition at line 460 of file FitMeasurement.h.

                                               {
  if (m_outlier)
    return 0;
  return m_numberDoF;
}

numericalDerivative()

bool Trk::FitMeasurement::numericalDerivative ( void ) const

inline

Definition at line 470 of file FitMeasurement.h.

                                                          {
  return m_numericalDerivative;
}

operator=()

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

private

perigee()

const Amg::VectorX & Trk::FitMeasurement::perigee ( void ) const

inline

Definition at line 474 of file FitMeasurement.h.

                                                           {
  return m_perigee;
}

perigeeWeight()

const Amg::MatrixX & Trk::FitMeasurement::perigeeWeight ( void ) const

inline

Definition at line 478 of file FitMeasurement.h.

                                                                 {
  return m_perigeeWeight;
}

position()

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

inline

Definition at line 482 of file FitMeasurement.h.

                                                             {
  return m_position;
}

print()

void Trk::FitMeasurement::print

(

MsgStream &

log

)

const

Definition at line 882 of file FitMeasurement.cxx.

                                               {
  Amg::Vector3D position = m_position;
  log << m_type << std::setiosflags(std::ios::fixed);
  if (numberDoF()) {
    log << std::setw(9) << std::setprecision(3) << *m_residual;
    if (m_numberDoF > 1) {
      log << std::setw(9) << std::setprecision(3) << *(m_residual + 1);
    } else if (m_alignmentParameter2) {
      log << " A" << std::setw(1) << m_alignmentParameter << " A"
          << std::setw(1) << m_alignmentParameter2 << "   ";
    } else if (m_alignmentParameter) {
      log << " A" << std::setw(1) << m_alignmentParameter << "      ";
    } else {
      log << "         ";
    }
  } else {
    if (isPositionMeasurement()) {
      log << std::setw(9) << std::setprecision(3) << *m_residual << " outlier ";
    } else {
      log << "                  ";
    }
    if (m_intersection[FittedTrajectory])
      position = m_intersection[FittedTrajectory]->position();
  }
  log << std::setw(10) << std::setprecision(1) << position.perp()
      << std::setw(9) << std::setprecision(4) << position.phi() << std::setw(10)
      << std::setprecision(1) << position(2);
 
  if (isPositionMeasurement()) {
    log << std::setw(13) << std::setprecision(3) << 1. / m_weight;
    if (m_numberDoF == 2) {
      log << std::setw(8) << std::setprecision(3) << 1. / m_weight2;
    } else if (isDrift()) {
      log << "(" << std::setw(7) << std::setprecision(3)
          << m_signedDriftDistance << ")";
    } else {
      log << "        ";
    }
  } else if (isScatterer()) {
    log << std::setw(33) << std::setprecision(3)
        << 1. / std::abs(m_qOverP * Gaudi::Units::GeV) << std::setw(12)
        << std::setprecision(4) << m_energyLoss / Gaudi::Units::GeV;
    if (m_type < barrelInert || m_scatteringAngle > 0.) {
      const double totScat = sqrt(m_scatteringAngle * std::abs(m_qOverP) *
                                m_scatteringAngle * std::abs(m_qOverP) +
                            m_scatteringAngleOffSet * m_scatteringAngleOffSet);
      log << std::setw(16) << std::setprecision(6) << totScat << std::setw(13)
          << std::setprecision(3) << m_radiationThickness;
    }
  } else if (isAlignment()) {
    log << std::setw(13) << std::setprecision(3) << 1. / m_weight;
    if (m_numberDoF == 2) {
      log << std::setw(8) << std::setprecision(3) << 1. / m_weight2;
    }
  } else if (isEnergyDeposit()) {
    if (m_numberDoF) {
      log << std::setw(13) << std::setprecision(3)
          << 1. / (m_weight * Gaudi::Units::GeV) << std::setw(32)
          << std::setprecision(4) << m_energyLoss / Gaudi::Units::GeV;
    } else {
      log << std::setw(33) << std::setprecision(3)
          << 1. / std::abs(m_qOverP * Gaudi::Units::GeV);
      log << std::setw(12) << std::setprecision(4)
          << m_energyLoss / Gaudi::Units::GeV;
    }
  }
  log << std::endl;
}

printHeading()

void Trk::FitMeasurement::printHeading ( MsgStream & log )

static

Definition at line 875 of file FitMeasurement.cxx.

                                                {
  log << "                 residual 1........2         r      phi         z"
      << "      sigma 1.......2      energy  energyLoss scatteringAngle  "
         "integral X0"
      << std::endl;
}

qOverP() [1/2]

void Trk::FitMeasurement::qOverP

(

double

value

)

Definition at line 951 of file FitMeasurement.cxx.

                                        {
  m_qOverP = value;
 
  // for scatterer measurements: correct the weight for the local momentum value
  if ((m_type == barrelScatterer || m_type == endcapScatterer ||
       m_type == calorimeterScatterer) &&
      m_scatteringAngle > 0.) {
    const double pSquare = 1. / (value * value);
    m_betaSquared = pSquare / (pSquare + m_particleMassSquared);
    m_weight = std::sqrt(m_betaSquared * pSquare) / m_scatteringAngle;
  }
 
  if (m_type == calorimeterScatterer && m_scatteringAngleOffSet > 0) {
    if (m_weight > 0) {
      m_weight = sqrt(1. / (1. / m_weight / m_weight +
                            m_scatteringAngleOffSet * m_scatteringAngleOffSet));
    } else {
      m_weight = 1. / m_scatteringAngleOffSet;
    }
  }
}

qOverP() [2/2]

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

inline

Definition at line 486 of file FitMeasurement.h.

                                               {
  return m_qOverP;
}

radiationThickness()

double Trk::FitMeasurement::radiationThickness ( void ) const

inline

Definition at line 490 of file FitMeasurement.h.

                                                           {
  return m_radiationThickness;
}

residual() [1/3]

void Trk::FitMeasurement::residual ( double value )

inline

Definition at line 498 of file FitMeasurement.h.

                                                 {
  *m_residual = value;
}

residual() [2/3]

void Trk::FitMeasurement::residual ( std::vector< double >::iterator pointer )

inline

Definition at line 502 of file FitMeasurement.h.

                                                                      {
  m_residual = pointer;
}

residual() [3/3]

double Trk::FitMeasurement::residual ( void ) const

inline

Definition at line 494 of file FitMeasurement.h.

                                                 {
  return *m_residual;
}

residual2() [1/2]

void Trk::FitMeasurement::residual2 ( double value )

inline

Definition at line 510 of file FitMeasurement.h.

                                                  {
  *(m_residual + 1) = value;
}

residual2() [2/2]

double Trk::FitMeasurement::residual2 ( void ) const

inline

Definition at line 506 of file FitMeasurement.h.

                                                  {
  return *(m_residual + 1);
}

scattererPhi() [1/2]

void Trk::FitMeasurement::scattererPhi ( double value )

inline

Definition at line 518 of file FitMeasurement.h.

                                                     {
  m_scatterPhi = value;
}

scattererPhi() [2/2]

double Trk::FitMeasurement::scattererPhi ( void ) const

inline

Definition at line 514 of file FitMeasurement.h.

                                                     {
  return m_scatterPhi;
}

scattererTheta() [1/2]

void Trk::FitMeasurement::scattererTheta ( double value )

inline

Definition at line 526 of file FitMeasurement.h.

                                                       {
  m_scatterPhi = value;
}

scattererTheta() [2/2]

double Trk::FitMeasurement::scattererTheta ( void ) const

inline

Definition at line 522 of file FitMeasurement.h.

                                                       {
  return m_scatterTheta;
}

scatteringAngle()

void Trk::FitMeasurement::scatteringAngle	(	double	angle,
		double	totalRadiationThickness )

Definition at line 973 of file FitMeasurement.cxx.

                                                                     {
  // update the m_scatteringAngleOffSet at initialisation
 
  if (m_scatteringAngle == 0. && m_scatteringAngleOffSet > 0 &&
      m_qOverP != 0.) {
    //
    const double angle_iPat = angle * std::abs(m_qOverP);
 
    //      std::cout << " scatteringAngle type " <<  m_type << " angle_iPat "
    //      << angle_iPat << " m_scatteringAngleOffSet " <<
    //      m_scatteringAngleOffSet;
 
    if (angle_iPat < m_scatteringAngleOffSet) {
      m_scatteringAngleOffSet =
          sqrt(m_scatteringAngleOffSet * m_scatteringAngleOffSet -
               angle_iPat * angle_iPat);
    } else {
      m_scatteringAngleOffSet = 0.;
    }
    //      std::cout << " corrected m_scatteringAngleOffSet " <<
    //      m_scatteringAngleOffSet << std::endl;
  }
 
  m_scatteringAngle = angle;
  m_radiationThickness = totalRadiationThickness;
  if (m_type == barrelInert) {
    m_type = barrelScatterer;
  } else if (m_type == endcapInert) {
    m_type = endcapScatterer;
  }
 
  if (m_qOverP != 0.) {
    const double pSquare = 1. / (m_qOverP * m_qOverP);
    m_betaSquared = pSquare / (pSquare + m_particleMassSquared);
    m_weight = std::sqrt(m_betaSquared * pSquare) / m_scatteringAngle;
    if (m_type == calorimeterScatterer && m_scatteringAngleOffSet > 0) {
      if (m_weight > 0) {
        m_weight =
            sqrt(1. / (1. / m_weight / m_weight +
                       m_scatteringAngleOffSet * m_scatteringAngleOffSet));
      } else {
        m_weight = 1. / m_scatteringAngleOffSet;
      }
    }
  }
}

sensorDirection()

const Amg::Vector3D & Trk::FitMeasurement::sensorDirection ( void ) const

inline

Definition at line 530 of file FitMeasurement.h.

                                                                    {
  return m_sensorDirection;
}

setEnergyGain()

void Trk::FitMeasurement::setEnergyGain ( void )

inline

Definition at line 534 of file FitMeasurement.h.

                                              {
  m_energyLoss = -fabs(m_energyLoss);
}

setMaterialEffectsOwner()

void Trk::FitMeasurement::setMaterialEffectsOwner ( void )

inline

Definition at line 538 of file FitMeasurement.h.

                                                        {
  m_materialEffectsOwner = true;
}

setNumericalDerivative()

void Trk::FitMeasurement::setNumericalDerivative ( void )

inline

Definition at line 542 of file FitMeasurement.h.

                                                       {
  m_numericalDerivative = true;
}

setOutlier()

void Trk::FitMeasurement::setOutlier ( void )

inline

Definition at line 546 of file FitMeasurement.h.

                                           {
  m_numberDoF = 0;
  m_outlier = true;
  m_weight = 1.0;
  m_weight2 = 1.0;
}

setSigma()

void Trk::FitMeasurement::setSigma ( void )

inline

Definition at line 553 of file FitMeasurement.h.

                                         {
  m_weight = 1. / m_sigma;
}

setSigmaMinus()

void Trk::FitMeasurement::setSigmaMinus ( void )

inline

Definition at line 557 of file FitMeasurement.h.

                                              {
  m_weight = 1. / m_sigmaMinus;
}

setSigmaPlus()

void Trk::FitMeasurement::setSigmaPlus ( void )

inline

Definition at line 561 of file FitMeasurement.h.

                                             {
  m_weight = 1. / m_sigmaPlus;
}

setSigmaSymmetric()

void Trk::FitMeasurement::setSigmaSymmetric

(

void

)

Definition at line 1021 of file FitMeasurement.cxx.

                                           {
  const double sigma = std::sqrt(
      0.5 * (m_sigmaPlus * m_sigmaPlus + m_sigmaMinus * m_sigmaMinus));
  m_weight = 1. / sigma;
}

sigma()

double Trk::FitMeasurement::sigma ( void ) const

inline

Definition at line 565 of file FitMeasurement.h.

                                              {
  if (!m_weight)
    return 0;
  return 1. / m_weight;
}

sigma2()

double Trk::FitMeasurement::sigma2 ( void ) const

inline

Definition at line 571 of file FitMeasurement.h.

                                               {
  if (!m_weight2)
    return 0;
  return 1. / m_weight2;
}

signedDriftDistance()

double Trk::FitMeasurement::signedDriftDistance ( void ) const

inline

Definition at line 577 of file FitMeasurement.h.

                                                            {
  return m_signedDriftDistance;
}

status()

int Trk::FitMeasurement::status ( void ) const

inline

Definition at line 581 of file FitMeasurement.h.

                                            {
  return m_status;
}

surface()

const Surface * Trk::FitMeasurement::surface ( void ) const

inline

Definition at line 585 of file FitMeasurement.h.

                                                        {
  return m_surface;
}

type()

MeasurementType Trk::FitMeasurement::type ( void ) const

inline

Definition at line 589 of file FitMeasurement.h.

                                                      {
  return m_type;
}

unsetOutlier()

void Trk::FitMeasurement::unsetOutlier ( void )

inline

Definition at line 593 of file FitMeasurement.h.

                                             {
  m_outlier = false;
}

weight()

double Trk::FitMeasurement::weight ( void ) const

inline

Definition at line 597 of file FitMeasurement.h.

                                               {
  return m_weight;
}

weight2()

double Trk::FitMeasurement::weight2 ( void ) const

inline

Definition at line 601 of file FitMeasurement.h.

                                                {
  return m_weight2;
}

Member Data Documentation

m_afterCalo

bool Trk::FitMeasurement::m_afterCalo

private

Definition at line 190 of file FitMeasurement.h.

m_alignmentEffects

const AlignmentEffectsOnTrack* Trk::FitMeasurement::m_alignmentEffects

private

Definition at line 191 of file FitMeasurement.h.

m_alignmentParameter

unsigned Trk::FitMeasurement::m_alignmentParameter

private

Definition at line 192 of file FitMeasurement.h.

m_alignmentParameter2

unsigned Trk::FitMeasurement::m_alignmentParameter2

private

Definition at line 193 of file FitMeasurement.h.

m_betaSquared

double Trk::FitMeasurement::m_betaSquared

private

Definition at line 194 of file FitMeasurement.h.

m_d0

double Trk::FitMeasurement::m_d0

private

Definition at line 198 of file FitMeasurement.h.

m_derivative

double* Trk::FitMeasurement::m_derivative

private

Definition at line 195 of file FitMeasurement.h.

m_derivative2

double* Trk::FitMeasurement::m_derivative2

private

Definition at line 196 of file FitMeasurement.h.

m_derivativeRow

int Trk::FitMeasurement::m_derivativeRow

private

Definition at line 197 of file FitMeasurement.h.

m_energyLoss

double Trk::FitMeasurement::m_energyLoss

private

Definition at line 199 of file FitMeasurement.h.

m_firstParameter

unsigned Trk::FitMeasurement::m_firstParameter

private

Definition at line 200 of file FitMeasurement.h.

m_flippedDriftDistance

bool Trk::FitMeasurement::m_flippedDriftDistance

private

Definition at line 201 of file FitMeasurement.h.

m_hitIndex

int Trk::FitMeasurement::m_hitIndex

private

Definition at line 202 of file FitMeasurement.h.

m_hitOnTrack

HitOnTrack* Trk::FitMeasurement::m_hitOnTrack

private

Definition at line 203 of file FitMeasurement.h.

m_intersection

std::array<std::optional<TrackSurfaceIntersection>, ExtrapolationTypes> Trk::FitMeasurement::m_intersection

private

Definition at line 205 of file FitMeasurement.h.

m_lastParameter

unsigned Trk::FitMeasurement::m_lastParameter

private

Definition at line 206 of file FitMeasurement.h.

m_materialEffects

const MaterialEffectsBase* Trk::FitMeasurement::m_materialEffects

private

Definition at line 207 of file FitMeasurement.h.

m_materialEffectsOwner

bool Trk::FitMeasurement::m_materialEffectsOwner

private

Definition at line 208 of file FitMeasurement.h.

m_measurementBase

const MeasurementBase* Trk::FitMeasurement::m_measurementBase

private

Definition at line 209 of file FitMeasurement.h.

m_minEnergyDeposit

double Trk::FitMeasurement::m_minEnergyDeposit

private

Definition at line 210 of file FitMeasurement.h.

m_minimizationDirection

Amg::Vector3D Trk::FitMeasurement::m_minimizationDirection

private

Definition at line 211 of file FitMeasurement.h.

m_normal

Amg::Vector3D Trk::FitMeasurement::m_normal

private

Definition at line 212 of file FitMeasurement.h.

m_numberDoF

int Trk::FitMeasurement::m_numberDoF

private

Definition at line 213 of file FitMeasurement.h.

m_numericalDerivative

bool Trk::FitMeasurement::m_numericalDerivative

private

Definition at line 214 of file FitMeasurement.h.

m_outlier

bool Trk::FitMeasurement::m_outlier

private

Definition at line 215 of file FitMeasurement.h.

m_particleMassSquared

double Trk::FitMeasurement::m_particleMassSquared

private

Definition at line 216 of file FitMeasurement.h.

m_perigee

Amg::VectorX Trk::FitMeasurement::m_perigee

private

Definition at line 217 of file FitMeasurement.h.

m_perigeeWeight

Amg::MatrixX Trk::FitMeasurement::m_perigeeWeight

private

Definition at line 218 of file FitMeasurement.h.

m_position

Amg::Vector3D Trk::FitMeasurement::m_position

private

Definition at line 219 of file FitMeasurement.h.

m_qOverP

double Trk::FitMeasurement::m_qOverP

private

Definition at line 220 of file FitMeasurement.h.

m_radiationThickness

double Trk::FitMeasurement::m_radiationThickness

private

Definition at line 221 of file FitMeasurement.h.

m_residual

std::vector<double>::iterator Trk::FitMeasurement::m_residual

private

Definition at line 222 of file FitMeasurement.h.

m_scaleFactor

double Trk::FitMeasurement::m_scaleFactor

private

Definition at line 223 of file FitMeasurement.h.

m_scatteringAngle

double Trk::FitMeasurement::m_scatteringAngle

private

Definition at line 226 of file FitMeasurement.h.

m_scatteringAngleOffSet

double Trk::FitMeasurement::m_scatteringAngleOffSet

private

Definition at line 227 of file FitMeasurement.h.

m_scatterPhi

double Trk::FitMeasurement::m_scatterPhi

private

Definition at line 224 of file FitMeasurement.h.

m_scatterTheta

double Trk::FitMeasurement::m_scatterTheta

private

Definition at line 225 of file FitMeasurement.h.

m_secondResidual

double Trk::FitMeasurement::m_secondResidual

private

Definition at line 228 of file FitMeasurement.h.

m_sensorDirection

Amg::Vector3D Trk::FitMeasurement::m_sensorDirection

private

Definition at line 229 of file FitMeasurement.h.

m_sigma

double Trk::FitMeasurement::m_sigma

private

Definition at line 230 of file FitMeasurement.h.

m_sigmaMinus

double Trk::FitMeasurement::m_sigmaMinus

private

Definition at line 231 of file FitMeasurement.h.

m_sigmaPlus

double Trk::FitMeasurement::m_sigmaPlus

private

Definition at line 232 of file FitMeasurement.h.

m_signedDriftDistance

double Trk::FitMeasurement::m_signedDriftDistance

private

Definition at line 233 of file FitMeasurement.h.

m_status

int Trk::FitMeasurement::m_status

private

Definition at line 234 of file FitMeasurement.h.

m_surface

const Surface* Trk::FitMeasurement::m_surface

private

Definition at line 235 of file FitMeasurement.h.

m_type

MeasurementType Trk::FitMeasurement::m_type

private

Definition at line 236 of file FitMeasurement.h.

m_weight

double Trk::FitMeasurement::m_weight

private

Definition at line 237 of file FitMeasurement.h.

m_weight2

double Trk::FitMeasurement::m_weight2

private

Definition at line 238 of file FitMeasurement.h.

---

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

• FitMeasurement.h
• FitMeasurement.cxx