CbmMustTube.cxx

Go to the documentation of this file.

/* Copyright (C) 2025 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Radoslaw Karabowicz [committer] */

// CbmMustTube header file
//
// Copied from PndFtsTube
//
// authors: Radoslaw Karabowicz, GSI, 2024
 
#include "CbmMustTube.h"
 
#include "CbmMustAddress.h"
 
CbmMustTube::CbmMustTube() {}
 
CbmMustTube::CbmMustTube(CbmMustTube& tube)
  : TObject(tube)
  , fTubeIndex(tube.GetTubeIndex())
  , fCenPosition(tube.GetPosition())
  , fRotationMatrix(TMatrixT<double>(3, 3))
  , fRadIn(tube.GetRadIn())
  , fRadOut(tube.GetRadOut())
  , fHalfLength(tube.GetHalfLength())
{
  fRotationMatrix.ResizeTo(3, 3);
  fRotationMatrix = tube.GetRotationMatrix();
}
 
CbmMustTube::CbmMustTube(uint32_t address, int tubeIndex, Double_t x, Double_t y, Double_t z, Double_t r11,
                         Double_t r12, Double_t r13, Double_t r21, Double_t r22, Double_t r23, Double_t r31,
                         Double_t r32, Double_t r33, Double_t radin, Double_t radout, Double_t hl)
  : fAddress(address)
  , fTubeIndex(tubeIndex)
  , fCenPosition(TVector3(x, y, z))
  , fRotationMatrix(TMatrixT<double>(3, 3))
  , fRadIn(radin)
  , fRadOut(radout)
  , fHalfLength(hl)
{
 
  fRotationMatrix[0][0] = r11;
  fRotationMatrix[0][1] = r12;
  fRotationMatrix[0][2] = r13;
 
  fRotationMatrix[1][0] = r21;
  fRotationMatrix[1][1] = r22;
  fRotationMatrix[1][2] = r23;
 
  fRotationMatrix[2][0] = r31;
  fRotationMatrix[2][1] = r32;
  fRotationMatrix[2][2] = r33;
}
 
CbmMustTube::~CbmMustTube() {}
 
int CbmMustTube::CompareTubeIndex(const CbmMustTube* compTube) const
{
  if (CbmMustAddress::GetElementAddress(compTube->GetAddress(), kMustModule)
      == CbmMustAddress::GetElementAddress(this->GetAddress(), kMustModule)) {
    if (this->GetTubeIndex() == compTube->GetTubeIndex()) return 0;
    if (this->GetTubeIndex() > compTube->GetTubeIndex()) return 1;
    return -1;
  }
  return -2;
}