AlignableEKLMRecoHit.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#include <alignment/reconstruction/AlignableEKLMRecoHit.h>
 
#include <alignment/GlobalDerivatives.h>
#include <alignment/GlobalLabel.h>
#include <klm/dataobjects/eklm/EKLMElementNumbers.h>
#include <klm/dataobjects/KLMDigit.h>
#include <klm/dataobjects/KLMElementNumbers.h>
#include <klm/dataobjects/KLMHit2d.h>
#include <klm/dbobjects/eklm/EKLMAlignment.h>
#include <klm/eklm/geometry/GeometryData.h>
#include <klm/eklm/geometry/TransformDataGlobalAligned.h>
 
using namespace Belle2;
using namespace alignment;
 
AlignableEKLMRecoHit::AlignableEKLMRecoHit()
{
}
 
AlignableEKLMRecoHit::AlignableEKLMRecoHit(
  const EKLMAlignmentHit* hit, const genfit::TrackCandHit* trackCandHit) :
  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()
{
}
 
std::pair<std::vector<int>, TMatrixD> AlignableEKLMRecoHit::globalDerivatives(const genfit::StateOnPlane* sop)
{
  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);
 
}
 
genfit::AbsMeasurement* AlignableEKLMRecoHit::clone() const
{
  return new AlignableEKLMRecoHit(*this);
}