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

Definition at line 21 of file Strip1DMeasurementTransformator.h.

Definition at line 20 of file Strip1DMeasurementTransformator.h.

Definition at line 22 of file Strip1DMeasurementTransformator.h.

Definition at line 19 of file Strip1DMeasurementTransformator.h.

Constructor & Destructor Documentation

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

Member Function Documentation

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; }
double Strip1DMeasurementTransformator::hitError ( ) const

Definition at line 29 of file Strip1DMeasurementTransformator.cc.

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

Referenced by KFStrip1DUpdator::update(), and ~Strip1DMeasurementTransformator().

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

17  {
18  return topology()->measurementPosition(hit().localPosition()).x();
19 }
virtual MeasurementPoint measurementPosition(const LocalPoint &) const =0
const StripTopology * topology() const
return real topology, taking sensor deformations into account
T x() const
Definition: PV2DBase.h:43
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; }
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 =0
virtual const GeomDetType & type() const
Definition: GeomDet.cc:69
virtual const Topology & topology() const
Definition: GeomDet.cc:67
virtual const GeomDetUnit * detUnit() const
double Strip1DMeasurementTransformator::projectedTrajectoryError ( ) const

Definition at line 37 of file Strip1DMeasurementTransformator.cc.

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

Referenced by KFStrip1DUpdator::update(), and ~Strip1DMeasurementTransformator().

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

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

Definition at line 41 of file Strip1DMeasurementTransformator.cc.

References MaterialEffects_cfi::A, funct::cos(), SimDataFormats::CaloAnalysis::cp, class-composition::H, hit(), idealTopology(), StripTopology::localPitch(), TrackingRecHit::localPosition(), phi, indexGen::s2, funct::sin(), mathSSE::sqrt(), state(), digitizers_cfi::strip, StripTopology::stripAngle(), OrderedSet::t, createJobs::tmp, topology(), PV3DBase< T, PVType, FrameType >::x(), and y.

Referenced by KFStrip1DUpdator::update(), and ~Strip1DMeasurementTransformator().

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
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
const StripTopology * idealTopology() const
return ideal topology, as if the sensor were flat
virtual float stripAngle(float strip) const =0
T sqrt(T t)
Definition: SSEVec.h:19
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
virtual LocalPoint localPosition() const =0
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
virtual float localPitch(const LocalPoint &) const =0
tmp
align.sh
Definition: createJobs.py:716
T x() const
Definition: PV3DBase.h:59
const TSOS& Strip1DMeasurementTransformator::state ( ) const
inline
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; }
const AlgebraicSymMatrix55 & Strip1DMeasurementTransformator::trajectoryError ( ) const

Definition at line 33 of file Strip1DMeasurementTransformator.cc.

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

Referenced by KFStrip1DUpdator::update(), and ~Strip1DMeasurementTransformator().

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

Member Data Documentation

const StripTopology * Strip1DMeasurementTransformator::theIdealTopology
private

Definition at line 47 of file Strip1DMeasurementTransformator.h.

Referenced by idealTopology(), and init().

const TrackingRecHit& Strip1DMeasurementTransformator::theRecHit
private

Definition at line 45 of file Strip1DMeasurementTransformator.h.

Referenced by hit().

TSOS Strip1DMeasurementTransformator::theState
private

Definition at line 46 of file Strip1DMeasurementTransformator.h.

Referenced by state().

const StripTopology* Strip1DMeasurementTransformator::theTopology
private

Definition at line 47 of file Strip1DMeasurementTransformator.h.

Referenced by init(), and topology().