AlignableBKLMRecoHit Class Reference

Alignable BKLM hit. More...

#include <AlignableBKLMRecoHit.h>

Inheritance diagram for AlignableBKLMRecoHit:

Public Member Functions
	AlignableBKLMRecoHit ()
	Constructor.
	AlignableBKLMRecoHit (const KLMHit2d *hit, const genfit::TrackCandHit *trackCandHit)
	Constructor.
virtual	~AlignableBKLMRecoHit ()
	Destructor.
genfit::AbsMeasurement *	clone () const override
	Creating a deep copy of this hit.
virtual std::vector< genfit::MeasurementOnPlane * >	constructMeasurementsOnPlane (const genfit::StateOnPlane &state) const override
	Measurement construction.
virtual const genfit::AbsHMatrix *	constructHMatrix (const genfit::AbsTrackRep *) const override
	Get the genfit projection H matrix (to U,V).
virtual std::pair< std::vector< int >, TMatrixD >	globalDerivatives (const genfit::StateOnPlane *sop) override
	Labels and derivatives of residuals (local measurement coordinates) w.r.t.

Private Types
enum	{ HIT_DIMENSIONS = 2 }

Private Member Functions
	ClassDefOverride (AlignableBKLMRecoHit, 1)
	Class version.

Private Attributes
uint16_t	m_KLMModule
	KLM module number.
int	m_Layer
	Hit layer.
const bklm::Module *	m_Module
	Module used to get geometry information.
CLHEP::Hep3Vector	global
	not streamed
const double	halfheight_sci = 0.5
	not streamed

Detailed Description

Alignable BKLM hit.

Definition at line 30 of file AlignableBKLMRecoHit.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
HIT_DIMENSIONS	sensitive Dimensions of the Hit

Definition at line 111 of file AlignableBKLMRecoHit.h.

{ HIT_DIMENSIONS = 2  };

Constructor & Destructor Documentation

AlignableBKLMRecoHit() [1/2]

AlignableBKLMRecoHit ( )

inline

Constructor.

Definition at line 38 of file AlignableBKLMRecoHit.h.

    {
    }

AlignableBKLMRecoHit() [2/2]

AlignableBKLMRecoHit ( const KLMHit2d *  hit, const genfit::TrackCandHit *  trackCandHit  )

explicit

Constructor.

the position (in global coordinates) of this module's sensitive-volume origin

the directioin (in global coordinates) of this module's U axis

the directioin (in global coordinates) of this module's V axis

Definition at line 25 of file AlignableBKLMRecoHit.cc.

                                                                                        :
  genfit::PlanarMeasurement(HIT_DIMENSIONS)
{
  int section = hit->getSection();
  int sector = hit->getSector();
  m_Layer = hit->getLayer();
  const KLMElementNumbers* elementNumbers = &(KLMElementNumbers::Instance());
  m_KLMModule = elementNumbers->moduleNumberBKLM(section, sector, m_Layer);
 
  bklm::GeometryPar* m_GeoPar = Belle2::bklm::GeometryPar::instance();
  m_Module = m_GeoPar->findModule(section, sector, m_Layer);
 
  //+++ global coordinates of the hit
  global[0] = hit->getPositionX();
  global[1] = hit->getPositionY();
  global[2] = hit->getPositionZ();
 
  //+++ local coordinates of the hit
  CLHEP::Hep3Vector local = m_Module->globalToLocal(global);
 
  //+++ local coordinates in KLM m_Layer of the hit
  double localU = local[1]; //phi
  double localV = local[2]; //z
  double errorU = m_Module->getPhiStripWidth() / sqrt(12);
  double errorV = m_Module->getZStripWidth() / sqrt(12);
 
  if (hit->inRPC()) {
    //+++ scale localU and localV based on measured-in-Belle resolution
    int nStrips = hit->getPhiStripMax() - hit->getPhiStripMin() + 1;
    double dn = nStrips - 1.5;
    double factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.60;//measured-in-Belle resolution
    errorU = errorU * sqrt(factor);
 
    nStrips = hit->getZStripMax() - hit->getZStripMin() + 1;
    dn = nStrips - 1.5;
    factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.55;//measured-in-Belle resolution
    errorV = errorV * sqrt(factor);
  }
 
  //Set positions
  rawHitCoords_(0) = localU;
  rawHitCoords_(1) = localV;
  //Set the error covariance matrix
  rawHitCov_(0, 0) = errorU * errorU; // error in U, squared
  rawHitCov_(0, 1) = 0.;
  rawHitCov_(1, 0) = 0.;
  rawHitCov_(1, 1) = errorV * errorV; // error in V, squared
 
  //+++ Construct vectors o, u, v of m_Layer in global coords.
  CLHEP::Hep3Vector gOrigin = m_Module->getGlobalOrigin();
  CLHEP::Hep3Vector lOrigin = m_Module->globalToLocal(gOrigin) + m_Module->getLocalReconstructionShift();
  CLHEP::Hep3Vector gOrigin_midway = m_Module->localToGlobal(lOrigin);
  CLHEP::Hep3Vector uAxis(0, 1, 0);
  CLHEP::Hep3Vector vAxis(0, 0, 1);
  CLHEP::Hep3Vector gUaxis = m_Module->localToGlobal(uAxis) - gOrigin;
  CLHEP::Hep3Vector gVaxis = m_Module->localToGlobal(vAxis) - gOrigin;
 
  B2Vector3D origin_mid(gOrigin_midway[0], gOrigin_midway[1], gOrigin_midway[2]);
 
  B2Vector3D uGlobal(gUaxis[0], gUaxis[1], gUaxis[2]);
  B2Vector3D vGlobal(gVaxis[0], gVaxis[1], gVaxis[2]);
 
  genfit::SharedPlanePtr detPlane(new genfit::DetPlane(origin_mid, uGlobal, vGlobal, 0));
  setPlane(detPlane, m_KLMModule);
}

~AlignableBKLMRecoHit()

virtual ~AlignableBKLMRecoHit ( )

inlinevirtual

Destructor.

Definition at line 51 of file AlignableBKLMRecoHit.h.

    {
    }

Member Function Documentation

clone()

genfit::AbsMeasurement * clone ( ) const

override

Creating a deep copy of this hit.

Definition at line 94 of file AlignableBKLMRecoHit.cc.

{
  return new AlignableBKLMRecoHit(*this);
}

constructHMatrix()

virtual const genfit::AbsHMatrix * constructHMatrix ( const genfit::AbsTrackRep *  ) const

inlineoverridevirtual

Get the genfit projection H matrix (to U,V).

Definition at line 70 of file AlignableBKLMRecoHit.h.

    {
      return new genfit::HMatrixUV();
    };

constructMeasurementsOnPlane()

std::vector< genfit::MeasurementOnPlane * > constructMeasurementsOnPlane ( const genfit::StateOnPlane &  state ) const

overridevirtual

Measurement construction.

Definition at line 99 of file AlignableBKLMRecoHit.cc.

{
  TVectorD predFglo = state.get6DState();
  TVectorD correctedLocal(2);
  CLHEP::Hep3Vector localmom(predFglo[3], predFglo[4], predFglo[5]);
  localmom = m_Module->RotateToLocal(localmom);
 
  //do correction for scintillators
  if (m_Layer == 1 || m_Layer == 2) {
    CLHEP::Hep3Vector global_shift_z(0, 0, predFglo[5]*halfheight_sci / sqrt(predFglo[3]*predFglo[3] + predFglo[4]*predFglo[4]));
    CLHEP::Hep3Vector local_corrected_z = m_Module->globalToLocal(global - global_shift_z);
    double local_shift_u = localmom[1] / localmom[0] * halfheight_sci;
    correctedLocal[0] = local_corrected_z[1] + local_shift_u;
    //correctedLocal[0] = local_corrected_z[1];
    correctedLocal[1] = local_corrected_z[2];
  } else {
    correctedLocal[0] = rawHitCoords_[0];
    correctedLocal[1] = rawHitCoords_[1];
  }
 
  // return std::vector<genfit::MeasurementOnPlane*>(1, new genfit::MeasurementOnPlane(rawHitCoords_, rawHitCov_, state.getPlane(),
  return std::vector<genfit::MeasurementOnPlane*>(1, new genfit::MeasurementOnPlane(correctedLocal, rawHitCov_, state.getPlane(),
                                                  state.getRep(), this->constructHMatrix(state.getRep())));
}

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 124 of file AlignableBKLMRecoHit.cc.

{
  std::vector<int> labGlobal;
 
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaU));
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaV));
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaW));
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaAlpha));
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaBeta));
  labGlobal.push_back(GlobalLabel::construct<BKLMAlignment>(m_KLMModule, KLMAlignmentData::c_DeltaGamma));
 
  // Matrix of global derivatives
  TMatrixD derGlobal(2, 6);
  derGlobal.Zero();
 
  // track u-slope in local sensor system
  double uSlope = sop->getState()[1];
  // track v-slope in local sensor system
  double vSlope = sop->getState()[2];
  // Predicted track u-position in local sensor system
  double uPos = sop->getState()[3];
  // Predicted track v-position in local sensor system
  double vPos = sop->getState()[4];
 
  // Global derivatives for alignment in module local coordinates
  derGlobal(0, 0) = 1.0;
  derGlobal(0, 1) = 0.0;
  derGlobal(0, 2) = - uSlope;
  derGlobal(0, 3) = vPos * uSlope;
  derGlobal(0, 4) = -uPos * uSlope;
  derGlobal(0, 5) = vPos;
 
  derGlobal(1, 0) = 0.0;
  derGlobal(1, 1) = 1.0;
  derGlobal(1, 2) = - vSlope;
  derGlobal(1, 3) = vPos * vSlope;
  derGlobal(1, 4) = -uPos * vSlope;
  derGlobal(1, 5) = -uPos;
 
  return alignment::GlobalDerivatives(labGlobal, derGlobal);
 
}

Member Data Documentation

global

CLHEP::Hep3Vector global

private

not streamed

global coordiante of the hit.

Definition at line 123 of file AlignableBKLMRecoHit.h.

halfheight_sci

const double halfheight_sci = 0.5

private

not streamed

half height of scintillator module

Definition at line 126 of file AlignableBKLMRecoHit.h.

m_KLMModule

uint16_t m_KLMModule

private

KLM module number.

Definition at line 114 of file AlignableBKLMRecoHit.h.

m_Layer

int m_Layer

private

Hit layer.

Definition at line 117 of file AlignableBKLMRecoHit.h.

m_Module

const bklm::Module* m_Module

private

Module used to get geometry information.

Definition at line 120 of file AlignableBKLMRecoHit.h.

---

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

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