AlignableSensorTRD.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#include "Align/AlignableSensorTRD.h"
#include "Align/AlignableDetectorTRD.h"
#include "Align/utils.h"
#include "Framework/Logger.h"
#include "Align/AlignmentPoint.h"
#include "DetectorsBase/GeometryManager.h"
 
namespace o2
{
namespace align
{
using namespace o2::align::utils;
using namespace TMath;
 
//_________________________________________________________
AlignableSensorTRD::AlignableSensorTRD(const char* name, int vid, int iid, int isec, Controller* ctr) : AlignableSensor(name, vid, iid, ctr), mSector(isec)
{
  // def c-tor
}
 
//____________________________________________
void AlignableSensorTRD::prepareMatrixClAlg()
{
  // prepare alignment matrix in the pseudo-LOCAL frame of TRD (account that the chamber has extra X,Y rotations
  TGeoHMatrix ma = getMatrixL2GIdeal().Inverse();
  ma *= getMatrixL2G();
  setMatrixClAlg(ma);
  //
}
 
//____________________________________________
void AlignableSensorTRD::prepareMatrixClAlgReco()
{
  // prepare alignment matrix in the pseudo-LOCAL frame of TRD (account that the chamber has extra X,Y rotations
  TGeoHMatrix ma = getMatrixL2GIdeal().Inverse();
  ma *= getMatrixL2G();
  setMatrixClAlgReco(ma);
  //
}
 
//____________________________________________
void AlignableSensorTRD::prepareMatrixL2GIdeal()
{
  TGeoHMatrix Rxy;
  Rxy.RotateX(-90);
  Rxy.RotateY(-90);
  TGeoHMatrix mtmp;
  if (!base::GeometryManager::getOriginalMatrix(getSymName(), mtmp)) {
    LOG(fatal) << "Failed to find ideal L2G matrix for " << getSymName();
  }
  mtmp *= Rxy;
  setMatrixL2GIdeal(mtmp);
}
 
//____________________________________________
void AlignableSensorTRD::prepareMatrixL2G(bool reco)
{
  TGeoHMatrix Rxy;
  Rxy.RotateX(-90);
  Rxy.RotateY(-90);
  const char* path = getSymName();
  const TGeoHMatrix* l2g = nullptr;
  if (!gGeoManager->GetAlignableEntry(path) || !(l2g = base::GeometryManager::getMatrix(path))) {
    LOGP(fatal, "Failed to find L2G matrix for {}alignable {} -> {}", gGeoManager->GetAlignableEntry(path) ? "" : "non-", path, (void*)l2g);
  }
  TGeoHMatrix mtmp = *l2g;
  mtmp *= Rxy;
  reco ? setMatrixL2GReco(mtmp) : setMatrixL2G(mtmp);
}
 
//____________________________________________
void AlignableSensorTRD::prepareMatrixT2L()
{
  // extract from geometry T2L matrix
  double alp = math_utils::detail::sector2Angle<float>(mSector);
  mAlp = alp;
  TGeoHMatrix Rs;
  Rs.RotateZ(-alp * TMath::RadToDeg());
  TGeoHMatrix m0 = getMatrixL2GIdeal();
  m0.MultiplyLeft(Rs);
  TGeoHMatrix t2l = m0.Inverse();
  setMatrixT2L(t2l);
  double loc[3] = {0, 0, 0}, glo[3];
  t2l.MasterToLocal(loc, glo);
  mX = glo[0];
  //
}
 
//____________________________________________
void AlignableSensorTRD::dPosTraDParCalib(const AlignmentPoint* pnt, double* deriv, int calibID, const AlignableVolume* parent) const
{
  // calculate point position X,Y,Z derivatives wrt calibration parameter calibID of given parent
  // parent=0 means top detector object calibration
  //
  deriv[0] = deriv[1] = deriv[2] = 0;
  //
  if (!parent) { // TRD detector global calibration
    //
    switch (calibID) {
      case AlignableDetectorTRD::CalibNRCCorrDzDtgl: // correction for Non-Crossing tracklets Z,Y shift: Z -> Z + calib*tgl, Y -> Y + calib*tgl*tilt*sign(tilt);
      {
        double sgYZ = pnt->getYZErrTracking()[1]; // makes sense only for nonRC tracklets
        if (std::abs(sgYZ) > 0.01) {
          const double kTilt = 2. * TMath::DegToRad();
          deriv[2] = pnt->getTrParamWSA()[AlignmentPoint::kParTgl];
          deriv[1] = deriv[2] * Sign(kTilt, sgYZ);
        }
        break;
      }
      case AlignableDetectorTRD::CalibDVT: // correction for bias in VdriftT
      {
        // error in VdriftT equivalent to shift in X at which Y measurement is evaluated
        // Y -> Y + dVdriftT * tg_phi, where tg_phi is the slope of the track in YX plane
        double snp = pnt->getTrParamWSA(AlignmentPoint::kParSnp), slpY = snp / std::sqrt((1.f - snp) * (1.f + snp));
        deriv[1] = slpY;
        break;
      }
      default:
        break;
    }
  }
}
 
} // namespace align
} // namespace o2