CMS 3D CMS Logo

List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
Strip1DMeasurementTransformator Class Reference

#include <Strip1DMeasurementTransformator.h>

Public Member Functions

const TrackingRecHithit () const
 
double hitError () const
 
double hitParameters () const
 
const StripTopologyidealTopology () const
 return ideal topology, as if the sensor were flat More...
 
double projectedTrajectoryError () const
 
double projectedTrajectoryParameters () const
 
AM15 projectionMatrix () const
 
const TSOSstate () const
 
 Strip1DMeasurementTransformator (const TSOS &aTsos, const TrackingRecHit &aHit)
 
const StripTopologytopology () const
 return real topology, taking sensor deformations into account More...
 
const ASM55trajectoryError () const
 
AV5 trajectoryParameters () const
 
 ~Strip1DMeasurementTransformator ()
 

Private Types

typedef AlgebraicMatrix15 AM15
 
typedef AlgebraicSymMatrix55 ASM55
 
typedef AlgebraicVector5 AV5
 
typedef TrajectoryStateOnSurface TSOS
 

Private Member Functions

void init ()
 

Private Attributes

const StripTopologytheIdealTopology
 
const TrackingRecHittheRecHit
 
TSOS theState
 
const StripTopologytheTopology
 

Detailed Description

Helper class for accessing the RecHit and the TrajectoryState parameters and errors in the measurement frame. The RecHits are treated as 1-dimensional, the second coordinate is ignored. Ported from ORCA.

Author
todorov, cerati

Definition at line 17 of file Strip1DMeasurementTransformator.h.

Member Typedef Documentation

◆ AM15

Definition at line 21 of file Strip1DMeasurementTransformator.h.

◆ ASM55

Definition at line 20 of file Strip1DMeasurementTransformator.h.

◆ AV5

Definition at line 22 of file Strip1DMeasurementTransformator.h.

◆ TSOS

Definition at line 19 of file Strip1DMeasurementTransformator.h.

Constructor & Destructor Documentation

◆ Strip1DMeasurementTransformator()

Strip1DMeasurementTransformator::Strip1DMeasurementTransformator ( const TSOS aTsos,
const TrackingRecHit aHit 
)

◆ ~Strip1DMeasurementTransformator()

Strip1DMeasurementTransformator::~Strip1DMeasurementTransformator ( )
inline

Definition at line 27 of file Strip1DMeasurementTransformator.h.

27 {}

Member Function Documentation

◆ hit()

const TrackingRecHit& Strip1DMeasurementTransformator::hit ( ) const
inline

Definition at line 37 of file Strip1DMeasurementTransformator.h.

References theRecHit.

Referenced by hitError(), hitParameters(), init(), and projectionMatrix().

37 { return theRecHit; }

◆ hitError()

double Strip1DMeasurementTransformator::hitError ( ) const

Definition at line 29 of file Strip1DMeasurementTransformator.cc.

References hit(), Topology::measurementError(), topology(), and MeasurementError::uu().

Referenced by KFStrip1DUpdator::update().

29  {
30  return topology()->measurementError(hit().localPosition(), hit().localPositionError()).uu();
31 }
float uu() const
virtual MeasurementError measurementError(const LocalPoint &, const LocalError &) const =0
const StripTopology * topology() const
return real topology, taking sensor deformations into account

◆ hitParameters()

double Strip1DMeasurementTransformator::hitParameters ( ) const

Definition at line 17 of file Strip1DMeasurementTransformator.cc.

References hit(), Topology::measurementPosition(), topology(), and PV2DBase< T, PVType, FrameType >::x().

Referenced by KFStrip1DUpdator::update().

17  {
18  return topology()->measurementPosition(hit().localPosition()).x();
19 }
T x() const
Definition: PV2DBase.h:43
virtual MeasurementPoint measurementPosition(const LocalPoint &) const =0
const StripTopology * topology() const
return real topology, taking sensor deformations into account

◆ idealTopology()

const StripTopology* Strip1DMeasurementTransformator::idealTopology ( ) const
inline

return ideal topology, as if the sensor were flat

Definition at line 40 of file Strip1DMeasurementTransformator.h.

References theIdealTopology.

Referenced by projectionMatrix().

40 { return theIdealTopology; }

◆ init()

void Strip1DMeasurementTransformator::init ( void  )
private

Definition at line 12 of file Strip1DMeasurementTransformator.cc.

References TrackingRecHit::detUnit(), hit(), theIdealTopology, theTopology, GeomDetType::topology(), GeomDet::topology(), and GeomDet::type().

Referenced by Strip1DMeasurementTransformator().

12  {
13  theTopology = dynamic_cast<const StripTopology*>(&(hit().detUnit()->topology()));
14  theIdealTopology = dynamic_cast<const StripTopology*>(&(hit().detUnit()->type().topology()));
15 }
virtual const Topology & topology() const
Definition: GeomDet.cc:67
virtual const GeomDetType & type() const
Definition: GeomDet.cc:69
virtual const GeomDetUnit * detUnit() const
virtual const Topology & topology() const =0

◆ projectedTrajectoryError()

double Strip1DMeasurementTransformator::projectedTrajectoryError ( ) const

Definition at line 37 of file Strip1DMeasurementTransformator.cc.

References Topology::measurementError(), state(), topology(), and MeasurementError::uu().

Referenced by KFStrip1DUpdator::update().

37  {
38  return topology()->measurementError(state().localPosition(), state().localError().positionError()).uu();
39 }
float uu() const
virtual MeasurementError measurementError(const LocalPoint &, const LocalError &) const =0
const StripTopology * topology() const
return real topology, taking sensor deformations into account

◆ projectedTrajectoryParameters()

double Strip1DMeasurementTransformator::projectedTrajectoryParameters ( ) const

Definition at line 25 of file Strip1DMeasurementTransformator.cc.

References Topology::measurementPosition(), state(), topology(), and PV2DBase< T, PVType, FrameType >::x().

Referenced by KFStrip1DUpdator::update().

25  {
26  return topology()->measurementPosition(state().localPosition()).x();
27 }
T x() const
Definition: PV2DBase.h:43
virtual MeasurementPoint measurementPosition(const LocalPoint &) const =0
const StripTopology * topology() const
return real topology, taking sensor deformations into account

◆ projectionMatrix()

AlgebraicMatrix15 Strip1DMeasurementTransformator::projectionMatrix ( ) const

Definition at line 41 of file Strip1DMeasurementTransformator.cc.

References A, funct::cos(), data-class-funcs::H, hit(), idealTopology(), StripTopology::localPitch(), TrackingRecHit::localPosition(), phi, funct::sin(), mathSSE::sqrt(), state(), nano_mu_digi_cff::strip, StripTopology::stripAngle(), submitPVValidationJobs::t, createJobs::tmp, topology(), PV3DBase< T, PVType, FrameType >::x(), and y.

Referenced by KFStrip1DUpdator::update().

41  {
42  // H(measurement <- local)
43  // m_meas = H*x_local + c
45  if (const RadialStripTopology* tmp = dynamic_cast<const RadialStripTopology*>(idealTopology())) {
46  double yHitToInter = tmp->yDistanceToIntersection(hit().localPosition().y());
47  double t = tmp->yAxisOrientation() * hit().localPosition().x() / yHitToInter;
48  double c2 = 1. / (1. + t * t); // cos(angle)**2
49  //double cs = t*c2; // sin(angle)*cos(angle); tan carries sign of sin!
50  double s2 = 1. - c2; // sin(angle)**2
51  double A = tmp->angularWidth();
52  // D is distance from intersection of edges to hit on strip
53  double D2 = hit().localPosition().x() * hit().localPosition().x() + yHitToInter * yHitToInter;
54  double D = std::sqrt(D2);
55 
56  double cp = std::sqrt(c2);
57  double sp;
58  if (t > 0) {
59  sp = std::sqrt(s2);
60  } else {
61  sp = -std::sqrt(s2);
62  }
63  H(0, 3) = cp / (D * A);
64  H(0, 4) = -sp / (D * A);
65  } else {
66  double phi = topology()->stripAngle(topology()->strip(state().localPosition()));
67  double pitch = topology()->localPitch(state().localPosition());
68  H(0, 3) = cos(phi) / pitch;
69  H(0, 4) = sin(phi) / pitch;
70  }
71  return H;
72 }
Divides< B, C > D2
Definition: Factorize.h:137
virtual float stripAngle(float strip) const =0
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
T x() const
Definition: PV3DBase.h:59
virtual float localPitch(const LocalPoint &) const =0
const StripTopology * idealTopology() const
return ideal topology, as if the sensor were flat
T sqrt(T t)
Definition: SSEVec.h:19
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
ROOT::Math::SMatrix< double, 1, 5, ROOT::Math::MatRepStd< double, 1, 5 > > AlgebraicMatrix15
const StripTopology * topology() const
return real topology, taking sensor deformations into account
DecomposeProduct< arg, typename Div::arg > D
Definition: Factorize.h:141
Definition: APVGainStruct.h:7
tmp
align.sh
Definition: createJobs.py:716
virtual LocalPoint localPosition() const =0

◆ state()

const TSOS& Strip1DMeasurementTransformator::state ( ) const
inline

◆ topology()

const StripTopology* Strip1DMeasurementTransformator::topology ( ) const
inline

return real topology, taking sensor deformations into account

Definition at line 42 of file Strip1DMeasurementTransformator.h.

References theTopology.

Referenced by hitError(), hitParameters(), projectedTrajectoryError(), projectedTrajectoryParameters(), and projectionMatrix().

42 { return theTopology; }

◆ trajectoryError()

const AlgebraicSymMatrix55 & Strip1DMeasurementTransformator::trajectoryError ( ) const

Definition at line 33 of file Strip1DMeasurementTransformator.cc.

References TrajectoryStateOnSurface::localError(), LocalTrajectoryError::matrix(), and state().

Referenced by KFStrip1DUpdator::update().

33  {
34  return state().localError().matrix();
35 }
const LocalTrajectoryError & localError() const
const AlgebraicSymMatrix55 & matrix() const

◆ trajectoryParameters()

AlgebraicVector5 Strip1DMeasurementTransformator::trajectoryParameters ( ) const

Definition at line 21 of file Strip1DMeasurementTransformator.cc.

References TrajectoryStateOnSurface::localParameters(), state(), and LocalTrajectoryParameters::vector().

Referenced by KFStrip1DUpdator::update().

21  {
22  return state().localParameters().vector();
23 }
const LocalTrajectoryParameters & localParameters() const
AlgebraicVector5 vector() const

Member Data Documentation

◆ theIdealTopology

const StripTopology * Strip1DMeasurementTransformator::theIdealTopology
private

Definition at line 47 of file Strip1DMeasurementTransformator.h.

Referenced by idealTopology(), and init().

◆ theRecHit

const TrackingRecHit& Strip1DMeasurementTransformator::theRecHit
private

Definition at line 45 of file Strip1DMeasurementTransformator.h.

Referenced by hit().

◆ theState

TSOS Strip1DMeasurementTransformator::theState
private

Definition at line 46 of file Strip1DMeasurementTransformator.h.

Referenced by state().

◆ theTopology

const StripTopology* Strip1DMeasurementTransformator::theTopology
private

Definition at line 47 of file Strip1DMeasurementTransformator.h.

Referenced by init(), and topology().