AlignableEKLMRecoHit Class Reference

Alignable EKLM hit. More...

#include <AlignableEKLMRecoHit.h>

Inheritance diagram for AlignableEKLMRecoHit:

Public Member Functions
	AlignableEKLMRecoHit ()
	Constructor.
	AlignableEKLMRecoHit (const EKLMAlignmentHit *hit, const genfit::TrackCandHit *trackCandHit)
	Constructor.
	~AlignableEKLMRecoHit ()
	Destructor.
virtual std::pair< std::vector< int >, TMatrixD >	globalDerivatives (const genfit::StateOnPlane *sop) override
	Labels and derivatives of residuals (local measurement coordinates) w.r.t.
genfit::AbsMeasurement *	clone () const override
	Clone.

Private Member Functions
	ClassDefOverride (AlignableEKLMRecoHit, 2)
	Needed to make objects storable.

Private Attributes
int	m_Section
	Section number.
int	m_Sector
	Sector number.
int	m_Layer
	Layer number.
uint16_t	m_KLMModule
	KLM module number.
uint16_t	m_Segment
	Segment number.
B2Vector3D	m_StripV
	V direction.

Detailed Description

Alignable EKLM hit.

Definition at line 29 of file AlignableEKLMRecoHit.h.

Constructor & Destructor Documentation

AlignableEKLMRecoHit() [1/2]

AlignableEKLMRecoHit

(

)

Constructor.

Definition at line 25 of file AlignableEKLMRecoHit.cc.

{
}

AlignableEKLMRecoHit() [2/2]

AlignableEKLMRecoHit	(	const EKLMAlignmentHit *	hit,
		const genfit::TrackCandHit *	trackCandHit
	)

Constructor.

Definition at line 29 of file AlignableEKLMRecoHit.cc.

                                                                       :
  genfit::PlanarMeasurement(1)
{
  (void)trackCandHit;
  int digit, plane, segment, strip;
  const HepGeom::Transform3D* t;
  CLHEP::HepRotation r;
  CLHEP::Hep3Vector origin;
  CLHEP::Hep3Vector u(1, 0, 0);
  CLHEP::Hep3Vector v(0, 1, 0);
  B2Vector3D origin2, u2, v2;
  const EKLM::GeometryData* geoDat = &(EKLM::GeometryData::Instance());
  const KLMElementNumbers* elementNumbers = &(KLMElementNumbers::Instance());
  const EKLMElementNumbers* eklmElementNumbers =
    &(EKLMElementNumbers::Instance());
  const EKLM::TransformDataGlobalAligned* transformData =
    &(EKLM::TransformDataGlobalAligned::Instance());
  RelationVector<KLMHit2d> hit2ds = hit->getRelationsTo<KLMHit2d>();
  if (hit2ds.size() != 1)
    B2FATAL("Incorrect number of related KLMHit2ds.");
  RelationVector<KLMDigit> eklmDigits = hit2ds[0]->getRelationsTo<KLMDigit>();
  if (eklmDigits.size() != 2)
    B2FATAL("Incorrect number of related KLMDigits.");
  digit = hit->getDigitIdentifier();
  m_Section = eklmDigits[digit]->getSection();
  m_Layer = eklmDigits[digit]->getLayer();
  m_Sector = eklmDigits[digit]->getSector();
  plane = eklmDigits[digit]->getPlane();
  segment = (eklmDigits[digit]->getStrip() - 1) /
            eklmElementNumbers->getNStripsSegment() + 1;
  strip = (segment - 1) * eklmElementNumbers->getNStripsSegment() + 1;
  m_KLMModule = elementNumbers->moduleNumberEKLM(m_Section, m_Sector, m_Layer);
  m_Segment = eklmElementNumbers->segmentNumber(
                m_Section, m_Layer, m_Sector, plane, segment);
  t = transformData->getStripTransform(
        m_Section, m_Layer, m_Sector, plane, strip);
  origin = t->getTranslation();
  origin2.SetX(origin.x() / CLHEP::cm * Unit::cm);
  origin2.SetY(origin.y() / CLHEP::cm * Unit::cm);
  origin2.SetZ(origin.z() / CLHEP::cm * Unit::cm);
  r = t->getRotation();
  u = r * u;
  v = r * v;
  u2.SetX(u.x());
  u2.SetY(u.y());
  u2.SetZ(u.z());
  v2.SetX(v.x());
  v2.SetY(v.y());
  v2.SetZ(v.z());
  t = transformData->getSectorTransform(m_Section, m_Layer, m_Sector);
  r = t->getRotation().inverse();
  v = r * v;
  m_StripV.SetX(v.unit().x());
  m_StripV.SetY(v.unit().y());
  genfit::SharedPlanePtr detPlane(new genfit::DetPlane(origin2, u2, v2, 0));
  setPlane(detPlane, m_Segment);
  rawHitCoords_[0] = geoDat->getStripGeometry()->getWidth() *
                     ((eklmDigits[digit]->getStrip() - 1) %
                      eklmElementNumbers->getNStripsSegment()) / CLHEP::cm * Unit::cm;
  rawHitCov_[0][0] = pow(geoDat->getStripGeometry()->getWidth() /
                         CLHEP::cm * Unit::cm, 2) / 12;
  setStripV();
}

~AlignableEKLMRecoHit()

~AlignableEKLMRecoHit

(

)

Destructor.

Definition at line 94 of file AlignableEKLMRecoHit.cc.

{
}

Member Function Documentation

clone()

genfit::AbsMeasurement * clone ( ) const

override

Clone.

Definition at line 142 of file AlignableEKLMRecoHit.cc.

{
  return new AlignableEKLMRecoHit(*this);
}

globalDerivatives()

std::pair< std::vector< int >, TMatrixD > globalDerivatives ( const genfit::StateOnPlane *  sop )

overridevirtual

Labels and derivatives of residuals (local measurement coordinates) w.r.t.

alignment/calibration parameters Matrix "G" of derivatives valid for given prediction of track state:

G(i, j) = d_residual_i/d_parameter_j

For 2D measurement (u,v):

G = ( du/da du/db du/dc ... ) ( dv/da dv/db dv/dc ... )

for calibration parameters a, b, c.

For 1D measurement:

G = ( 0 0 0 ... ) ( dv/da dv/db dv/dc ... ) for V-strip,

G = ( du/da du/db du/dc ... ) ( 0 0 0 ... ) for U-strip,

Measurements with more dimesions (slopes, curvature) should provide full 4-5Dx(n params) matrix (state as (q/p, u', v', u, v) or (u', v', u, v))

Parameters

sop	Predicted state of the track as linearization point around which derivatives of alignment/calibration parameters shall be computed

Returns

pair<vector<int>, TMatrixD> With matrix with number of rows = dimension of residual, number of columns = number of parameters. number of columns must match vector<int>.size().

Definition at line 98 of file AlignableEKLMRecoHit.cc.

{
  std::vector<int> labGlobal;
  labGlobal.push_back(GlobalLabel::construct<EKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaU)); // dx
  labGlobal.push_back(GlobalLabel::construct<EKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaV)); // dy
  labGlobal.push_back(GlobalLabel::construct<EKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaGamma)); // drot
 
  /* Local parameters. */
  const double dalpha = 0;
  const double dxs = 0;
  const double dys = 0;
  const double sinda = sin(dalpha);
  const double cosda = cos(dalpha);
  /* Local position in sector coordinates. */
  HepGeom::Point3D<double> globalPos;
  HepGeom::Transform3D t;
  const EKLM::TransformDataGlobalAligned* transformData =
    &(EKLM::TransformDataGlobalAligned::Instance());
  t = (*transformData->getSectorTransform(m_Section, m_Layer, m_Sector)).inverse();
  globalPos.setX(sop->getPos().X() / Unit::cm * CLHEP::cm);
  globalPos.setY(sop->getPos().Y() / Unit::cm * CLHEP::cm);
  globalPos.setZ(sop->getPos().Z() / Unit::cm * CLHEP::cm);
  globalPos = t * globalPos;
  double x = globalPos.x() / CLHEP::cm * Unit::cm;
  double y = globalPos.y() / CLHEP::cm * Unit::cm;
  /*
   * Matrix of global derivatives of local coordinate (second dimension is
   * added because of resizing in GblFitterInfo::constructGblPoint()) by
   * variations of the detector element position.
   * The sign is inverted to match the convention for other detectors.
   */
  TMatrixD derGlobal(2, 3);
  derGlobal(0, 0) = 0;
  derGlobal(0, 1) = 0;
  derGlobal(0, 2) = 0;
  derGlobal(1, 0) = cosda * m_StripV.X() - sinda * m_StripV.Y();
  derGlobal(1, 1) = sinda * m_StripV.X() + cosda * m_StripV.Y();
  derGlobal(1, 2) =
    -((x - dxs) * (-sinda * m_StripV.X() - cosda * m_StripV.Y()) +
      (y - dys) * (cosda * m_StripV.X() - sinda * m_StripV.Y()));
  return alignment::GlobalDerivatives(labGlobal, derGlobal);
 
}

Member Data Documentation

m_KLMModule

uint16_t m_KLMModule

private

KLM module number.

Definition at line 100 of file AlignableEKLMRecoHit.h.

m_Layer

int m_Layer

private

Layer number.

Definition at line 97 of file AlignableEKLMRecoHit.h.

m_Section

int m_Section

private

Section number.

Definition at line 91 of file AlignableEKLMRecoHit.h.

m_Sector

int m_Sector

private

Sector number.

Definition at line 94 of file AlignableEKLMRecoHit.h.

m_Segment

uint16_t m_Segment

private

Segment number.

Definition at line 103 of file AlignableEKLMRecoHit.h.

m_StripV

B2Vector3D m_StripV

private

V direction.

Definition at line 106 of file AlignableEKLMRecoHit.h.

---

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

• alignment/reconstruction/include/AlignableEKLMRecoHit.h
• alignment/reconstruction/src/AlignableEKLMRecoHit.cc