GeometryBase.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 "TRDBase/GeometryBase.h"
 
using namespace o2::trd;
using namespace o2::trd::constants;
 
//_____________________________________________________________________________
GPUd() int GeometryBase::getStack(float z, int layer) const
{
  //
  // Reconstruct the chamber number from the z position and layer number
  //
  // The return function has to be protected for positiveness !!
  //
 
  if ((layer < 0) || (layer >= NLAYER)) {
    return -1;
  }
 
  int istck = NSTACK;
  float zmin = 0.0;
  float zmax = 0.0;
 
  do {
    istck--;
    if (istck < 0) {
      break;
    }
    const PadPlane* pp = getPadPlane(layer, istck);
    zmax = pp->getRow0();
    int nrows = pp->getNrows();
    zmin = zmax - 2 * pp->getLengthOPad() - (nrows - 2) * pp->getLengthIPad() - (nrows - 1) * pp->getRowSpacing();
  } while ((z < zmin) || (z > zmax));
 
  return istck;
}
 
//_____________________________________________________________________________
GPUd() bool GeometryBase::isOnBoundary(int det, float y, float z, float eps) const
{
  //
  // Checks whether position is at the boundary of the sensitive volume
  //
 
  int ly = getLayer(det);
  if ((ly < 0) || (ly >= NLAYER)) {
    return true;
  }
 
  int stk = getStack(det);
  if ((stk < 0) || (stk >= NSTACK)) {
    return true;
  }
 
  const PadPlane* pp = &mPadPlanes[getDetectorSec(ly, stk)];
  if (!pp) {
    return true;
  }
 
  float max = pp->getRow0();
  int n = pp->getNrows();
  float min = max - 2 * pp->getLengthOPad() - (n - 2) * pp->getLengthIPad() - (n - 1) * pp->getRowSpacing();
  if (z < min + eps || z > max - eps) {
    // printf("z : min[%7.2f (%7.2f)] %7.2f max[(%7.2f) %7.2f]\n", min, min+eps, z, max-eps, max);
    return true;
  }
  min = pp->getCol0();
  n = pp->getNcols();
  max = min + 2 * pp->getWidthOPad() + (n - 2) * pp->getWidthIPad() + (n - 1) * pp->getColSpacing();
  if (y < min + eps || y > max - eps) {
    // printf("y : min[%7.2f (%7.2f)] %7.2f max[(%7.2f) %7.2f]\n", min, min+eps, y, max-eps, max);
    return true;
  }
 
  return false;
}