CaloSurfaceBuilder.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
//****************************************************************************
// CaloSurfacaBuilder.cxx, (c) ATLAS detector software
//****************************************************************************
 
#include "CaloTrackingGeometry/CaloSurfaceBuilder.h"
 
#include "Gaudi/Property.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/IService.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/IToolSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SystemOfUnits.h"
 
#include "TrkGeometry/CylinderLayer.h"
#include "TrkSurfaces/CylinderBounds.h"
#include "TrkSurfaces/CylinderSurface.h"
 
#include "TrkGeometry/DiscLayer.h"
#include "TrkSurfaces/DiscBounds.h"
#include "TrkSurfaces/DiscSurface.h"
 
#include "TrkGeometry/HomogeneousLayerMaterial.h"
#include "TrkGeometry/MaterialProperties.h"
 
// Amg
#include "GeoPrimitives/GeoPrimitives.h"
 
// Calo specific stuff :
#include "CaloDetDescr/CaloDetDescrElement.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
#include "CaloDetDescr/CaloDetDescriptor.h"
#include "CaloDetDescr/ICaloCoordinateTool.h"
#include "CaloDetDescr/ICaloRecoMaterialTool.h"
#include "CaloDetDescr/ICaloRecoSimpleGeomTool.h"
 
// Tile specific stuff :
#include "CaloIdentifier/TileID.h"
#include "TileDetDescr/TileDetDescrManager.h"
#include "TileDetDescr/TileDetDescriptor.h"
 
#include "TrkDetDescrInterfaces/ITrackingVolumeBuilder.h"
#include "TrkGeometry/TrackingVolume.h"
#include "TrkGeometrySurfaces/SlidingCylinderSurface.h"
#include "TrkGeometrySurfaces/SlidingDiscSurface.h"
#include "TrkVolumes/CylinderVolumeBounds.h"
 
#include <algorithm> //for std::min, std::max
#include <cmath>
 
CaloSurfaceBuilder::CaloSurfaceBuilder(const std::string& type,
                                       const std::string& name,
                                       const IInterface* parent)
  : AthAlgTool(type, name, parent)
  , m_tile_dd(nullptr)
  , m_lar_mat("LArRecoMaterialTool")
  , m_lar_simplegeom("LArRecoSimpleGeomTool")
{
  declareInterface<ICaloSurfaceBuilder>(this);
  declareProperty("LArRecoMaterialTool", m_lar_mat);
  declareProperty("LarRecoSimpleGeometryTool", m_lar_simplegeom);
}
 
CaloSurfaceBuilder::~CaloSurfaceBuilder() = default;
 
StatusCode
CaloSurfaceBuilder::initialize()
{
  if (m_lar_mat.retrieve().isFailure()) {
    ATH_MSG_FATAL("Failed to retrieve " << m_lar_mat << ". Aborting.");
    return StatusCode::FAILURE;
  }
 
  if (m_lar_simplegeom.retrieve().isFailure()) {
    ATH_MSG_FATAL("Failed to retrieve " << m_lar_simplegeom << ". Aborting.");
    return StatusCode::FAILURE;
  }
 
  if (m_calodepth.retrieve().isFailure()) {
    ATH_MSG_FATAL("Failed to retrieve " << m_calodepth << ". Aborting.");
    return StatusCode::FAILURE;
  }
 
  if (detStore()->retrieve(m_tile_dd).isFailure()) {
    ATH_MSG_FATAL("Could not find TileDetDescrManager in DetStore");
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode
CaloSurfaceBuilder::finalize()
{
  // clear surface cash
  for (auto & layerEntry : m_layerEntries) {
    delete layerEntry.first;
    delete layerEntry.second;
  }
 
  if (!m_layerExits.empty()) {
    delete m_layerExits[CaloCell_ID::TileBar2].first;
    delete m_layerExits[CaloCell_ID::TileBar2].second;
    delete m_layerExits[CaloCell_ID::TileExt2].first;
    delete m_layerExits[CaloCell_ID::TileExt2].second;
    delete m_layerExits[CaloCell_ID::TileGap2].first;
    delete m_layerExits[CaloCell_ID::TileGap2].second;
    delete m_layerExits[CaloCell_ID::TileGap3].first;
    delete m_layerExits[CaloCell_ID::TileGap3].second;
    delete m_layerExits[CaloCell_ID::HEC3].first;
    delete m_layerExits[CaloCell_ID::HEC3].second;
  }
 
  StatusCode sc = StatusCode::SUCCESS;
  return sc;
}
 
Trk::Surface*
CaloSurfaceBuilder::CreateUserSurface(const CaloCell_ID::CaloSample sample,
                                      const double offset,
                                      const double etaCaloLocal,
                                      const CaloDetDescrManager* calo_dd) const
{
 
  if (sample == CaloCell_ID::Unknown) {
    return nullptr;
  }
 
  auto pos = Amg::Transform3D(Trk::s_idTransform);
  Trk::Surface* surf = nullptr;
 
  bool result = false;
  double rmin = 0.;
  double z = 0.;
  double rmax = 0.;
  double radius = 0.;
  double hphi = 0.;
  double hlen = 0.;
  double depth = 0.;
  int side = 1;
 
  bool tile = false;
  if (sample == CaloCell_ID::TileBar0 || sample == CaloCell_ID::TileBar1 ||
      sample == CaloCell_ID::TileBar2 || sample == CaloCell_ID::TileExt0 ||
      sample == CaloCell_ID::TileExt1 || sample == CaloCell_ID::TileExt2 ||
      sample == CaloCell_ID::TileGap1 || sample == CaloCell_ID::TileGap2 ||
      sample == CaloCell_ID::TileGap3)
    tile = true;
 
  bool HEC3 = false;
  if (sample == CaloCell_ID::HEC3)
    HEC3 = true;
 
  if (etaCaloLocal < 0)
    side = -1;
 
  // Specific hack for the testbeam : even if tracks have eta<0, ask for the
  // positive side
  if (calo_dd->lar_geometry() == "H8")
    side = 1;
 
  // too bad, I need this ugly if to decide which type of surface should be
  // built cylinder of disk ?
  // Note: TileGap3 is in the barrel but is a disk!!
  bool barrel = true;
  if (sample == CaloCell_ID::PreSamplerE || sample == CaloCell_ID::EME1 ||
      sample == CaloCell_ID::EME2 || sample == CaloCell_ID::EME3 ||
      sample == CaloCell_ID::HEC0 || sample == CaloCell_ID::HEC1 ||
      sample == CaloCell_ID::HEC2 || sample == CaloCell_ID::HEC3 ||
      sample == CaloCell_ID::FCAL0 || sample == CaloCell_ID::FCAL1 ||
      sample == CaloCell_ID::FCAL2 || sample == CaloCell_ID::TileGap3) {
    barrel = false;
  }
 
  if (barrel) {
 
    /*To avoid the funny numbers given by the CaloDescriptors to  TileBar2, (due
      to the cell D type)
      we'll take only the radius from them, the rest we copy from TileBar1 */
    if (sample == CaloCell_ID::TileBar2) {
      double radius2 = 0;
      auto  pos2 = Amg::Transform3D(Trk::s_idTransform);
      result = this->get_cylinder_surface(
        sample, side, pos2, radius, hphi, hlen, depth, calo_dd);
      result = this->get_cylinder_surface(CaloCell_ID::TileBar1,
                                          side,
                                          pos,
                                          radius2,
                                          hphi,
                                          hlen,
                                          depth,
                                          calo_dd);
    } else {
      result = this->get_cylinder_surface(
        sample, side, pos, radius, hphi, hlen, depth, calo_dd);
    }
    if (!result) {
      return nullptr;
    }
 
    // this correction will only work for LAr
    double betterrad = 0;
    if (!tile) {
      betterrad = m_calodepth->radius(sample, etaCaloLocal, 0., calo_dd);
      radius = offset + betterrad;
    }
    if (tile)
      hlen = hlen / 2;
    // for Atlas configuration, avoid to create a hole at phi ~ pi by hphi
    // rounding effects  (the phi range was introduced for the testbeam)
    if (hphi < 0.9 * M_PI) {
      surf = new Trk::CylinderSurface(pos, radius, hphi, hlen);
    } else {
      surf = new Trk::CylinderSurface(pos, radius, hlen);
    }
  } else {
    bool result = this->get_disk_surface(
      sample, side, pos, z, rmin, rmax, hphi, depth, calo_dd);
    if (!result) {
      return nullptr;
    }
 
    double betterz = 0;
 
    if (!tile && !HEC3) {
      betterz = m_calodepth->radius(sample, etaCaloLocal, 0., calo_dd);
      z = offset + betterz - (pos)(2, 3);
      Amg::Translation3D vec(0., 0., z);
      pos = (pos) * vec;
    }
    // for Atlas configuration, avoid to create a hole a phi ~ pi by hphi
    // rounding effects. the phi range was introduced for the testbeam
    if (hphi < 0.9 * M_PI) {
      surf = new Trk::DiscSurface(pos, rmin, rmax, hphi);
    } else {
      surf = new Trk::DiscSurface(pos, rmin, rmax);
    }
  }
 
  return surf;
}
 
Trk::Surface*
CaloSurfaceBuilder::CreateLastSurface(const CaloCell_ID::CaloSample sample,
                                      const double offset,
                                      const double etaCaloLocal,
                                      const CaloDetDescrManager* calo_dd) const
{
 
  Trk::Surface* surf = nullptr;
 
  bool tile = false;
  if (sample == CaloCell_ID::TileBar0 || sample == CaloCell_ID::TileBar1 ||
      sample == CaloCell_ID::TileBar2 || sample == CaloCell_ID::TileExt0 ||
      sample == CaloCell_ID::TileExt1 || sample == CaloCell_ID::TileExt2 ||
      sample == CaloCell_ID::TileGap1 || sample == CaloCell_ID::TileGap2 ||
      sample == CaloCell_ID::TileGap3)
    tile = true;
 
  auto pos = Amg::Transform3D(Trk::s_idTransform);
 
  bool result = false;
  double rmin = 0.;
  double z = 0.;
  double rmax = 0.;
  double radius = 0.;
  double hphi = 0.;
  double hlen = 0.;
  double depth = 0.;
  int side = 1;
 
  if (etaCaloLocal < 0)
    side = -1;
 
  if (sample == CaloCell_ID::TileBar2 || sample == CaloCell_ID::TileExt2 ||
      sample == CaloCell_ID::TileGap2) {
 
    /*To avoid the funny numbers given by the CaloDescriptors to  TileBar2, (due
      to the cell D type)
      we'll take only the radius from them, the rest we copy from TileBar1 */
    if (sample == CaloCell_ID::TileBar2) {
      double radius2 = 0;
      auto pos2 = Amg::Transform3D(Trk::s_idTransform);
      result = this->get_cylinder_surface(
        sample, side, pos2, radius, hphi, hlen, depth, calo_dd);
      result = this->get_cylinder_surface(CaloCell_ID::TileBar1,
                                          side,
                                          pos,
                                          radius2,
                                          hphi,
                                          hlen,
                                          depth,
                                          calo_dd);
    } else {
 
      result = this->get_cylinder_surface(
        sample, side, pos, radius, hphi, hlen, depth, calo_dd);
    }
    if (!result) {
      return nullptr;
    }
 
    // for Atlas configuration, avoid to create a hole at phi ~ pi by hphi
    // rounding effects  (the phi range was introduced for the testbeam)
    else {
      std::vector<TileDetDescriptor*>::const_iterator Tile_first =
        m_tile_dd->tile_descriptors_begin();
      TileDetDescriptor* TileBar_descr = *(Tile_first + 1);
      TileDetDescriptor* TileExt_descr = *(Tile_first + 3);
      TileDetDescriptor* TileGap_descr = *(Tile_first + 5);
      if (sample == CaloCell_ID::TileBar2) {
        radius = TileBar_descr->rcenter(2) + .5 * TileBar_descr->dr(2);
      } else if (sample == CaloCell_ID::TileExt2) {
        radius = TileExt_descr->rcenter(2) + .5 * TileExt_descr->dr(2);
      } else if (sample == CaloCell_ID::TileGap2) {
        radius = TileGap_descr->rcenter(2) + .5 * TileGap_descr->dr(2);
      } else {
        ATH_MSG_DEBUG("Not a considered cylinder tile");
      }
 
      if (tile) {
        hlen = hlen / 2;
      }
 
      if (hphi < 0.9 * M_PI) {
        return surf = new Trk::CylinderSurface(pos, radius, hphi, hlen);
      } else {
        return surf = new Trk::CylinderSurface(pos, radius, hlen);
      }
    }
  } else if (sample == CaloCell_ID::TileGap3 || sample == CaloCell_ID::HEC3) {
 
    bool result = this->get_disk_surface(
      sample, side, pos, z, rmin, rmax, hphi, depth, calo_dd);
    if (!result) {
      return nullptr;
    }
    double zend = 0;
    if (sample == CaloCell_ID::TileGap3) {
      std::vector<TileDetDescriptor*>::const_iterator Tile_first =
        m_tile_dd->tile_descriptors_begin();
      TileDetDescriptor* TileGap_descr = *(Tile_first + 5);
      zend = TileGap_descr->zcenter(2) + .5 * TileGap_descr->dz(2);
      zend = zend * side;
    }
    if (sample == CaloCell_ID::HEC3) {
      for (const CaloDetDescriptor* descr : calo_dd->calo_descriptors_range()) {
        if (descr) {
          if (descr->getSampling(0) == CaloCell_ID::HEC3) {
            zend = descr->calo_z_max();
            zend = zend * side;
          }
        }
      }
    }
 
    z = offset + zend - (pos)(2, 3);
 
    Amg::Translation3D vec(0., 0., z);
    pos = (pos) * vec;
 
    // for Atlas configuration, avoid to create a hole a phi ~ pi by hphi
    // rounding effects. the phi range was introduced for the testbeam
    if (hphi < 0.9 * M_PI) {
      return new Trk::DiscSurface(pos, rmin, rmax, hphi);
    } else {
      return new Trk::DiscSurface(pos, rmin, rmax);
    }
  } else {
    ATH_MSG_ERROR(
      "This sample is not a 'LAST' surface... returning dummy surface");
    return surf;
  }
 
  return nullptr;
}
 
CaloDepthTool*
CaloSurfaceBuilder::getCaloDepth()
{
  return &(*m_calodepth); // from ToolHandle to pointer
}
 
// End of the general user methods :
 
bool
CaloSurfaceBuilder::get_cylinder_surface(
  CaloCell_ID::CaloSample sample,
  int side,
  Amg::Transform3D& htrans,
  double& radius,
  double& hphi,
  double& hlength,
  double& depth,
  const CaloDetDescrManager* calo_dd) const
{
  bool result = false;
 
  // strips are spread on several descriptor, which all have the same
  // htrans, radius, hphi, but not the same hlength
  double hl{};
  hlength = 0;
 
  // these ones are definitely not cylinders :
  if (sample == CaloCell_ID::PreSamplerE || sample == CaloCell_ID::EME1 ||
      sample == CaloCell_ID::EME2 || sample == CaloCell_ID::EME3 ||
      sample == CaloCell_ID::HEC0 || sample == CaloCell_ID::HEC1 ||
      sample == CaloCell_ID::HEC2 || sample == CaloCell_ID::HEC3 ||
      sample == CaloCell_ID::FCAL0 || sample == CaloCell_ID::FCAL1 ||
      sample == CaloCell_ID::FCAL2 || sample == CaloCell_ID::Unknown ||
      sample == CaloCell_ID::TileGap3) {
    return false;
  }
  if (sample == CaloCell_ID::PreSamplerB || sample == CaloCell_ID::EMB1 ||
      sample == CaloCell_ID::EMB2 || sample == CaloCell_ID::EMB3) {
 
    for (const CaloDetDescriptor* reg : calo_dd->calo_descriptors_range()) {
      if (reg) {
        if (reg->getSampling(0) == sample && reg->calo_sign() * side > 0) {
          result = reg->get_cylinder_surface(htrans, radius, hphi, hl, depth);
          hlength = std::max(hl, hlength);
        }
      }
    }
  }
 
  else {
    // these are not cylinders
    switch (sample) {
      case CaloCell_ID::TileBar0:
      case CaloCell_ID::TileBar1:
      case CaloCell_ID::TileBar2:
      case CaloCell_ID::TileExt0:
      case CaloCell_ID::TileExt1:
      case CaloCell_ID::TileExt2:
      case CaloCell_ID::TileGap1:
      case CaloCell_ID::TileGap2:
        break;
      default:
        ATH_MSG_ERROR(" get_cylinder_surface Not Cylindrical Tile found!!!   ");
        return 0;
    }
 
    for (const CaloDetDescriptor* reg : calo_dd->tile_descriptors_range()) {
      if (reg) {
        if (reg->getSampling(0) == sample && reg->calo_sign() * side > 0) {
          result = reg->get_cylinder_surface(htrans, radius, hphi, hl, depth);
          if (hl > hlength)
            hlength = hl;
          hphi = .5 * reg->dphi() * reg->n_phi();
        }
      }
    }
    result = true;
  }
  return result;
}
 
bool
CaloSurfaceBuilder::get_disk_surface(CaloCell_ID::CaloSample sample,
                                     int side,
                                     Amg::Transform3D& htrans,
                                     double& z,
                                     double& rmin,
                                     double& rmax,
                                     double& hphisec,
                                     double& depth,
                                     const CaloDetDescrManager* calo_dd) const
{
 
  // set defaults :
  bool result = false;
  z = 0.;
  rmin = 999999.;
  rmax = 0.;
 
  // strips are spread on several descriptor, which all have the same
  // htrans, hphisec, but not the same rmin and rmax
  // + for EMEC2 and EMEC3, the InnerWeel entrance z is different :-(
  // ==> as CaloCell_ID::CaloSample does not allow to distinguish,
  //     take OW z and fix the HepTransform accordingly
 
  // these ones are definitely not disks :
  if (sample == CaloCell_ID::PreSamplerB || sample == CaloCell_ID::EMB1 ||
      sample == CaloCell_ID::EMB2 || sample == CaloCell_ID::EMB3 ||
      sample == CaloCell_ID::TileBar0 || sample == CaloCell_ID::TileBar1 ||
      sample == CaloCell_ID::TileBar2 || sample == CaloCell_ID::TileGap1 ||
      sample == CaloCell_ID::TileGap2 || sample == CaloCell_ID::TileExt0 ||
      sample == CaloCell_ID::TileExt1 || sample == CaloCell_ID::TileExt2) {
    return false;
  }
 
  double ri, ra, zow;
  for (const CaloDetDescriptor* reg : calo_dd->calo_descriptors_range()) {
    if (reg) {
      if (reg->getSampling(0) == sample && reg->calo_sign() * side > 0) {
        result = reg->get_disk_surface(htrans, zow, ri, ra, hphisec, depth);
        if (!reg->is_lar_em_endcap_inner())
          z = zow;
        rmin = std::min(ri, rmin);
        rmax = std::max(rmax, ra);
      }
    }
  }
 
  if (sample == CaloCell_ID::TileGap3) {
    for (const CaloDetDescriptor* reg : calo_dd->tile_descriptors_range()) {
      if (reg) {
        if (reg->getSampling(0) == sample && reg->calo_sign() * side > 0) {
          result = reg->get_disk_surface(htrans, zow, ri, ra, hphisec, depth);
          if (!reg->is_lar_em_endcap_inner())
            z = zow;
          rmin = std::min(ri, rmin);
          rmax = std::max(rmax, ra);
        }
      }
    }
  }
 
  if (rmin > rmax) {
    rmin = 0.;
  }
  // Fix the z ( convention is different in Calo's and Tracking )
  Amg::Translation3D vec(0., 0., z - (htrans)(2, 3));
  htrans = (htrans) * vec;
  return result;
}
 
bool
CaloSurfaceBuilder::get_cylinder_surface(CaloSubdetNames::ALIGNVOL alvol,
                                         Amg::Transform3D& htrans,
                                         double& hphi,
                                         std::vector<double>& radius,
                                         std::vector<double>& depth,
                                         std::vector<double>& hlength) const
{
  bool result = m_lar_simplegeom->get_cylinder_surface(
    alvol, htrans, hphi, radius, depth, hlength);
 
  // Shift the POS/NEG sides  ( htrans convention is different in Calo's and
  // Tracking )
  //    FIXME : for this, assume that the 1rst hlength is the good one
 
  if (alvol == CaloSubdetNames::EMB_POS ||
      alvol == CaloSubdetNames::PRESAMPLER_B_POS) {
    result = m_lar_simplegeom->get_cylinder_surface(
      alvol, htrans, hphi, radius, depth, hlength);
    Amg::Translation3D vec(0., 0., hlength[0]);
    htrans = (htrans) * vec;
    return result;
  } else if (alvol == CaloSubdetNames::EMB_NEG ||
             alvol == CaloSubdetNames::PRESAMPLER_B_NEG) {
    result = m_lar_simplegeom->get_cylinder_surface(
      alvol, htrans, hphi, radius, depth, hlength);
    Amg::Translation3D vec(0., 0., -1. * hlength[0]);
    htrans = (htrans) * vec;
    return result;
  }
 
  return result;
}
 
bool
CaloSurfaceBuilder::get_disk_surface(CaloSubdetNames::ALIGNVOL alvol,
                                     Amg::Transform3D& htrans,
                                     double& hphi,
                                     std::vector<double>& z,
                                     std::vector<double>& depth,
                                     std::vector<double>& rmin,
                                     std::vector<double>& rmax) const
{
  bool result = m_lar_simplegeom->get_disk_surface(
    alvol, htrans, hphi, z, depth, rmin, rmax);
 
  // Fix the z ( htrans convention is different in Calo's and Tracking ) ?
  //   = left to the CreateDDEClayer method, because z is a vector => 1 fix per
  //   disk element
 
  return result;
}
 
// store all the surfaces into a vector
void
CaloSurfaceBuilder::fill_tg_surfaces(const CaloDetDescrManager* calo_dd) const
{
 
  // entry surfaces ( id<24 to avoid error messages )
  for (CaloCell_ID::CaloSample sample = CaloCell_ID::PreSamplerB; sample < 24;
       sample = CaloCell_ID::CaloSample(sample + 1)) {
    float etaRef = (sample < 4 || std::abs(sample - 13) < 2) ? 1. : 2.;
    Trk::Surface* spos = CreateUserSurface(sample, 0., etaRef, calo_dd);
    Trk::Surface* sneg = CreateUserSurface(sample, 0., -etaRef, calo_dd);
    if (spos)
      spos->setOwner(Trk::TGOwn);
    if (sneg)
      sneg->setOwner(Trk::TGOwn);
    m_layerEntries.emplace_back(spos, sneg);
  }
  // turn into SlidingSurfaces
  // EMB2 neg
  float etaMin = -1.5;
  float etaMax = 0.;
  float deta = 0.05;
  unsigned int neta = (etaMax - etaMin) / deta;
  std::vector<float> offset(neta);
  float rbase = m_layerEntries[CaloCell_ID::EMB2].second->bounds().r();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    const Trk::Surface* surf =
      CreateUserSurface(CaloCell_ID::EMB2, 0., eta, calo_dd);
    offset[ie] = surf->bounds().r() - rbase;
    delete surf;
  }
  const Trk::CylinderSurface* cyl = dynamic_cast<const Trk::CylinderSurface*>(
    m_layerEntries[CaloCell_ID::EMB2].second);
  if (!cyl) {
    std::abort();
  }
  Trk::SlidingCylinderSurface* tmpSlidCyl = new Trk::SlidingCylinderSurface(
    *cyl,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidCyl->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EMB2].second = tmpSlidCyl;
  delete cyl;
  // EMB2 pos
  etaMin = 0.;
  etaMax = 1.5;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  rbase = m_layerEntries[CaloCell_ID::EMB2].first->bounds().r();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    Trk::Surface* surf = CreateUserSurface(CaloCell_ID::EMB2, 0., eta, calo_dd);
    offset[ie] = surf->bounds().r() - rbase;
    delete surf;
  }
  cyl = dynamic_cast<const Trk::CylinderSurface*>(
    m_layerEntries[CaloCell_ID::EMB2].first);
  if (!cyl) {
    std::abort();
  }
  tmpSlidCyl = new Trk::SlidingCylinderSurface(
    *cyl,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidCyl->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EMB2].first = tmpSlidCyl;
  delete cyl;
  // EMB3 neg
  etaMin = -1.5;
  etaMax = 0.;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  rbase = m_layerEntries[CaloCell_ID::EMB3].second->bounds().r();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    Trk::Surface* surf = CreateUserSurface(CaloCell_ID::EMB3, 0., eta, calo_dd);
    offset[ie] = surf->bounds().r() - rbase;
    delete surf;
  }
  cyl = dynamic_cast<const Trk::CylinderSurface*>(
    m_layerEntries[CaloCell_ID::EMB3].second);
  if (!cyl) {
    std::abort();
  }
  tmpSlidCyl = new Trk::SlidingCylinderSurface(
    *cyl,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidCyl->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EMB3].second = tmpSlidCyl;
  delete cyl;
 
  // EMB3 pos
  etaMin = 0.;
  etaMax = 1.5;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  rbase = m_layerEntries[CaloCell_ID::EMB3].first->bounds().r();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    const Trk::Surface* surf =
      CreateUserSurface(CaloCell_ID::EMB3, 0., eta, calo_dd);
    offset[ie] = surf->bounds().r() - rbase;
    delete surf;
  }
  cyl = dynamic_cast<const Trk::CylinderSurface*>(
    m_layerEntries[CaloCell_ID::EMB3].first);
  if (!cyl) {
    std::abort();
  }
  tmpSlidCyl = new Trk::SlidingCylinderSurface(
    *cyl,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidCyl->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EMB3].first = tmpSlidCyl;
  delete cyl;
 
  // EME2 neg
  etaMin = -3.2;
  etaMax = -1.35;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  float zbase = m_layerEntries[CaloCell_ID::EME2].second->center().z();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    const Trk::Surface* surf =
      CreateUserSurface(CaloCell_ID::EME2, 0., eta, calo_dd);
    offset[ie] = surf->center().z() - zbase;
    delete surf;
  }
  const Trk::DiscSurface* disc = dynamic_cast<const Trk::DiscSurface*>(
    m_layerEntries[CaloCell_ID::EME2].second);
  if (!disc) {
    std::abort();
  }
  Trk::SlidingDiscSurface* tmpSlidDisc = new Trk::SlidingDiscSurface(
    *disc,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidDisc->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EME2].second = tmpSlidDisc;
  delete disc;
 
  // EME2 pos
  etaMin = 1.35;
  etaMax = 3.2;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  zbase = m_layerEntries[CaloCell_ID::EME2].first->center().z();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    Trk::Surface* surf = CreateUserSurface(CaloCell_ID::EME2, 0., eta, calo_dd);
    offset[ie] = surf->center().z() - zbase;
    delete surf;
  }
  disc = dynamic_cast<const Trk::DiscSurface*>(
    m_layerEntries[CaloCell_ID::EME2].first);
  if (!disc) {
    std::abort();
  }
  tmpSlidDisc = new Trk::SlidingDiscSurface(
    *disc,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidDisc->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EME2].first = tmpSlidDisc;
  delete disc;
 
  // EME3 neg
  etaMin = -3.2;
  etaMax = -1.35;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  zbase = m_layerEntries[CaloCell_ID::EME3].second->center().z();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    const Trk::Surface* surf =
      CreateUserSurface(CaloCell_ID::EME3, 0., eta, calo_dd);
    offset[ie] = surf->center().z() - zbase;
    delete surf;
  }
  disc = dynamic_cast<const Trk::DiscSurface*>(
    m_layerEntries[CaloCell_ID::EME3].second);
  if (!disc) {
    std::abort();
  }
  tmpSlidDisc = new Trk::SlidingDiscSurface(
    *disc,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidDisc->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EME3].second = tmpSlidDisc;
  delete disc;
 
  // EME3 pos
  etaMin = 1.35;
  etaMax = 3.2;
  deta = 0.05;
  neta = (etaMax - etaMin) / deta;
  offset.resize(neta);
  zbase = m_layerEntries[CaloCell_ID::EME3].first->center().z();
  for (unsigned int ie = 0; ie < neta; ie++) {
    float eta = etaMin + (ie + 0.5) * deta;
    Trk::Surface* surf = CreateUserSurface(CaloCell_ID::EME3, 0., eta, calo_dd);
    offset[ie] = surf->center().z() - zbase;
    delete surf;
  }
  disc = dynamic_cast<const Trk::DiscSurface*>(
    m_layerEntries[CaloCell_ID::EME3].first);
  if (!disc) {
    std::abort();
  }
  tmpSlidDisc = new Trk::SlidingDiscSurface(
    *disc,
    Trk::BinUtility(neta, etaMin, etaMax, Trk::open, Trk::binEta),
    offset);
  tmpSlidDisc->setOwner(Trk::TGOwn);
  m_layerEntries[CaloCell_ID::EME3].first = tmpSlidDisc;
  delete disc;
 
  // exit surfaces
  for (CaloCell_ID::CaloSample sample = CaloCell_ID::PreSamplerB;
       sample + 1 < CaloCell_ID::Unknown;
       sample = CaloCell_ID::CaloSample(sample + 1)) {
    m_layerExits.push_back(m_layerEntries[sample + 1]);
  }
 
  // fix exit for outer layers
  Trk::Surface* lpos =
    CreateLastSurface(CaloCell_ID::TileBar2, 0., 1., calo_dd);
  Trk::Surface* lneg =
    CreateLastSurface(CaloCell_ID::TileBar2, 0., -1., calo_dd);
  if (lpos) {
    lpos->setOwner(Trk::TGOwn);
  }
  if (lneg) {
    lneg->setOwner(Trk::TGOwn);
  }
  m_layerExits[CaloCell_ID::TileBar2] =
    std::pair<const Trk::Surface*, const Trk::Surface*>(lpos, lneg);
 
  lpos = CreateLastSurface(CaloCell_ID::TileExt2, 0., 1., calo_dd);
  lneg = CreateLastSurface(CaloCell_ID::TileExt2, 0., -1., calo_dd);
  if (lpos) {
    lpos->setOwner(Trk::TGOwn);
  }
  if (lneg) {
    lneg->setOwner(Trk::TGOwn);
  }
  m_layerExits[CaloCell_ID::TileExt2] =
    std::pair<const Trk::Surface*, const Trk::Surface*>(lpos, lneg);
 
  lpos = CreateLastSurface(CaloCell_ID::TileGap2, 0., 1., calo_dd);
  lneg = CreateLastSurface(CaloCell_ID::TileGap2, 0., -1., calo_dd);
  if (lpos) {
    lpos->setOwner(Trk::TGOwn);
  }
  if (lneg) {
    lneg->setOwner(Trk::TGOwn);
  }
  m_layerExits[CaloCell_ID::TileGap2] =
    std::pair<const Trk::Surface*, const Trk::Surface*>(lpos, lneg);
 
  lpos = CreateLastSurface(CaloCell_ID::TileGap3, 0., 1., calo_dd);
  lneg = CreateLastSurface(CaloCell_ID::TileGap3, 0., -1., calo_dd);
  if (lpos) {
    lpos->setOwner(Trk::TGOwn);
  }
  if (lneg) {
    lneg->setOwner(Trk::TGOwn);
  }
  m_layerExits[CaloCell_ID::TileGap3] =
    std::pair<const Trk::Surface*, const Trk::Surface*>(lpos, lneg);
 
  lpos = CreateLastSurface(CaloCell_ID::HEC3, 0., 1., calo_dd);
  lneg = CreateLastSurface(CaloCell_ID::HEC3, 0., -1., calo_dd);
  if (lpos) {
    lpos->setOwner(Trk::TGOwn);
  }
  if (lneg) {
    lneg->setOwner(Trk::TGOwn);
  }
  m_layerExits[CaloCell_ID::HEC3] =
    std::pair<const Trk::Surface*, const Trk::Surface*>(lpos, lneg);
}