LayerArrayCreator.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/

// LayerArrayCreator.cxx, (c) ATLAS Detector software
 
// Trk include
#include "TrkDetDescrTools/LayerArrayCreator.h"
#include "TrkDetDescrUtils/BinnedArray1D.h"
#include "TrkDetDescrUtils/BinUtility.h"
#include "TrkDetDescrUtils/GeometryStatics.h"
#include "TrkGeometry/Layer.h"
#include "TrkGeometry/NavigationLayer.h"
#include "TrkSurfaces/CylinderSurface.h"
#include "TrkSurfaces/DiscSurface.h"
#include "TrkSurfaces/DiscBounds.h"
#include "TrkSurfaces/PlaneSurface.h"
#include "TrkSurfaces/SurfaceBounds.h"
#include "TrkSurfaces/TrapezoidBounds.h"
#include "TrkSurfaces/RectangleBounds.h"
// Amg
#include "GeoPrimitives/GeoPrimitives.h"
 
 
// constructor
Trk::LayerArrayCreator::LayerArrayCreator(const std::string& t, const std::string& n, const IInterface* p) : AthAlgTool(t,n,p)
{
    declareInterface<ILayerArrayCreator>(this);
}
 
 
std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>> Trk::LayerArrayCreator::cylinderLayerArray(
  const std::vector<Trk::CylinderLayer*>& cylLayersInput,
  double rmin, double rmax, Trk::BinningType btype) const
{
    ATH_MSG_VERBOSE( " build LayerArray with " << cylLayersInput.size() << " cylindrical material layers." );
    ATH_MSG_VERBOSE( "       rmin/rmax provided : " << rmin << " / " << rmax );
 
    //copy so that you can sort
    std::vector<Trk::CylinderLayer*> cylLayers(cylLayersInput);
    // sort the vector
    Trk::CylinderLayerSorterR rSorter;
    std::sort(cylLayers.begin(), cylLayers.end(), rSorter);
 
    // needed for all cases
    Trk::Layer*                         cylinderLayer      = nullptr;
    std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>>    cylinderLayerArray = nullptr;
    std::vector< std::pair<std::shared_ptr<Trk::Layer>, Amg::Vector3D> >   layerOrderVector;
 
    switch (btype) {
 
        // equidistant binning - no navigation layers built
        case Trk::equidistant :
        {
            // count the layers
            unsigned int layers = cylLayers.size();
            // loop over layers and put them in
            for (auto& layIter : cylLayers ) {
                // get the R
                const Trk::CylinderSurface& layerSurface = layIter->surfaceRepresentation();
                double currentR = layerSurface.bounds().r();
                ATH_MSG_VERBOSE( "equidistant : registering cylindrical MaterialLayer at   radius : " << currentR );
                layerOrderVector.emplace_back(std::shared_ptr<Layer>(layIter),
                                              Amg::Vector3D(currentR, 0.,0.));
            }
            // create the binUtility
            auto binUtility = Trk::BinUtility(layers,rmin,rmax,Trk::open, Trk::binR);
            ATH_MSG_VERBOSE( "equidistant : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            cylinderLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
        } break;
 
        // bi-equidistant binning - take care the layers have to be binned equidistant + same thickness
        case Trk::biequidistant :
        {
 
            // count the layers
            unsigned int layers = cylLayers.size();
            // take a reference thinkness
            double layerThickness = cylLayers[0]->thickness();
            // the radialstep
            double radialStep = (rmax-rmin)/(layers-1);
            // the next step
            double navigationR                     = 0.;
            double navLayerHalflengthZ             = 0.;
            const Amg::Transform3D* layerTransform = nullptr;
 
            // loop over layers
            for (auto& layIter : cylLayers ) {
              // get the dimensions
              const Trk::CylinderSurface& layerSurface = layIter->surfaceRepresentation();
              layerTransform = layerSurface.transform().isApprox(Amg::Transform3D::Identity()) ? nullptr : &layerSurface.transform();
 
              double currentR = layerSurface.bounds().r();
              navigationR = currentR - 0.5*radialStep;
 
              navLayerHalflengthZ = layerSurface.bounds().halflengthZ();
              ATH_MSG_VERBOSE( "bi-equidistant : creating cylindrical NavigationLayer at   radius : " << navigationR );
              auto navLayerSurface =
                  layerTransform ? std::make_unique<Trk::CylinderSurface>(
                                       Amg::Transform3D(*layerTransform),navigationR, navLayerHalflengthZ)
                                 : std::make_unique<Trk::CylinderSurface>(
                                       navigationR, navLayerHalflengthZ);
              // the navigation layer
              auto navLayer = std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface));
              // push the navigation layer in
              layerOrderVector.emplace_back(std::move(navLayer),
                                            Amg::Vector3D(navigationR, 0., 0.));
              ATH_MSG_VERBOSE( "bi-equidistant : registering cylindrical MaterialLayer at   radius : " << currentR );
              // push the original layer in
              layerOrderVector.emplace_back(std::shared_ptr<Trk::Layer>(layIter),
                                            Amg::Vector3D(currentR, 0., 0.));
            }
 
            // special treatment for the last one
            ATH_MSG_VERBOSE( "bi-equidistant : creating cylindrical NavigationLayer at   radius : " << navigationR+radialStep);
            auto navLayerSurfacFinal = layerTransform ?
                            std::make_unique<Trk::CylinderSurface>(
                              Amg::Transform3D(*layerTransform), navigationR+radialStep, navLayerHalflengthZ) :
                            std::make_unique<Trk::CylinderSurface>(navigationR+radialStep, navLayerHalflengthZ);
            // the navigation layer
            auto navLayer = std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurfacFinal));
            // push the navigation layer in
            layerOrderVector.emplace_back(std::move(navLayer),
                                          Amg::Vector3D(navigationR+radialStep, 0., 0.));
 
            ATH_MSG_VERBOSE( layerOrderVector.size() << " cylindrical Layers (material + navigation) built. " );
 
            // create the binUtility
            double rMinBoundary = rmin-radialStep+0.5*layerThickness;
            double rMaxBoundary = rmax+radialStep+0.5*layerThickness;
            auto binUtility = Trk::BinUtility(layers, layerThickness, rMinBoundary, rMaxBoundary, Trk::open, Trk::binR);
            ATH_MSG_VERBOSE( "bi-equidistant : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            cylinderLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // arbitrary binning
        case Trk::arbitrary :
        {
            std::vector<float> boundaries;
            // maz z extension
            double halfLengthZ                     = 0;
            const Amg::Transform3D* layerTransform = nullptr;
 
            // initial step
            boundaries.push_back(rmin);
 
            // loop over the provided layers and put Navigation layers in between
            for (auto& layIter : cylLayers) {
                // get the cylinder surface
                const Trk::CylinderSurface& layerSurface = layIter->surfaceRepresentation();
                layerTransform = layerSurface.transform().isApprox(Amg::Transform3D::Identity()) ? nullptr : &layerSurface.transform();
                // and get the halflength
                double currentHalfLengthZ = layerSurface.bounds().halflengthZ();
                takeBigger(halfLengthZ, currentHalfLengthZ);
                double layerRadius = layerSurface.bounds().r();
                double layerThickness = layIter->thickness();
                // navigation layer : previous bin
                double navLayerRadius = 0.5*( (layerRadius-0.5*layerThickness) + boundaries[boundaries.size()-1] );
                auto navLayerSurface = layerTransform ?
                                    std::make_unique<Trk::CylinderSurface>(
                                      Amg::Transform3D(*layerTransform), navLayerRadius, halfLengthZ) :
                                    std::make_unique<Trk::CylinderSurface>(navLayerRadius, halfLengthZ);
                // material layer : current bin
                cylinderLayer = checkAndReplaceEmptyLayer(layIter);
                if (cylinderLayer){
                    ATH_MSG_VERBOSE( "arbitrary : creating cylindrical NavigationLayer at radius : " << navLayerRadius );
                    layerOrderVector.emplace_back(
                                               std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                               Amg::Vector3D(navLayerRadius, 0., 0.));
                    ATH_MSG_VERBOSE( "arbitrary : registering cylindrical MaterialLayer at   radius :" << layerRadius );
                    layerOrderVector.emplace_back(
                                               std::shared_ptr<Trk::Layer>(cylinderLayer),
                                               Amg::Vector3D(layerRadius, 0.,0.));
                    boundaries.push_back(layerRadius-0.5*layerThickness);
                    boundaries.push_back(layerRadius+0.5*layerThickness);
                } else {
                    ATH_MSG_VERBOSE( "arbitrary : empty layer configuration cancelled this building of navigation layer.");
                }
            }
            // close up the array with last bin
            double navLayerRadiusFinal = 0.5*(rmax+boundaries[boundaries.size()-1]);
            auto navLayerSurfaceFinal = layerTransform ?
                                        std::make_unique<Trk::CylinderSurface>(
                                          Amg::Transform3D(*layerTransform), navLayerRadiusFinal, halfLengthZ) :
                                        std::make_unique<Trk::CylinderSurface>(navLayerRadiusFinal, halfLengthZ);
            boundaries.push_back(rmax);
            ATH_MSG_VERBOSE( "arbitrary : creating cylindrical NavigationLayer at radius : " << navLayerRadiusFinal );
            layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurfaceFinal)),
                                          Amg::Vector3D(navLayerRadiusFinal, 0., 0.));
 
            ATH_MSG_VERBOSE( layerOrderVector.size() << " cylindrical Layers (material + navigation) built. " );
            // create the BinUtility
            auto binUtility = Trk::BinUtility(boundaries, Trk::open, Trk::binR);
            ATH_MSG_VERBOSE( "arbitrary : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            cylinderLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // default return 0
        default : { return nullptr; }
    }
 
    return cylinderLayerArray;
}
 
 
std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>> Trk::LayerArrayCreator::discLayerArray(
  const std::vector<Trk::DiscLayer*>& discLayersInput,
  double zmin,
  double zmax,
  Trk::BinningType btype) const
{
 
    ATH_MSG_VERBOSE( " build LayerArray with " << discLayersInput.size() << " disc-like material layers." );
    ATH_MSG_VERBOSE( "       zmin/zmax provided : " << zmin << " / " << zmax );
 
    // needed for all cases
    std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>>  discLayerArray = nullptr;
    std::vector<std::pair<std::shared_ptr<Trk::Layer>, Amg::Vector3D>>   layerOrderVector;
 
    //copy so that you can sort
    std::vector<Trk::DiscLayer*> discLayers(discLayersInput);
    // sort the vector
    Trk::DiscLayerSorterZ zSorter;
    std::sort(discLayers.begin(), discLayers.end(), zSorter);
    // the layer to be pushed back
    Trk::Layer* discLayer      = nullptr;
 
    switch (btype) {
 
        // equidistant binning
        case Trk::equidistant :
        {
            // count the layers
            size_t layers = discLayers.size();
            // loop over layers and put them in
            for (Trk::DiscLayer*  layIter : discLayers ) {
                discLayer = checkAndReplaceEmptyLayer(layIter);
                // get the Z
                const Trk::Surface& layerSurface = discLayer->surfaceRepresentation();
                ATH_MSG_VERBOSE( "equidistant : registering disc-like MaterialLayer   at z-Position : " << layerSurface.center().z() );
                layerOrderVector.emplace_back(std::shared_ptr<Layer>(discLayer),
                                              layerSurface.center());
            }
            // create the binUitlity
            auto binUtility = Trk::BinUtility(layers,zmin,zmax,Trk::open,Trk::binZ);
            ATH_MSG_VERBOSE( "equidistant : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            discLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // bi-equidistant binning
        case Trk::biequidistant :
        {
            // count the layers
            unsigned int layers = discLayers.size();
            // take a reference layer thickness
            double layerThickness = discLayers[0]->thickness();
            ATH_MSG_VERBOSE( "bi-equidistant : zmin / zmax re-evaluated as = " << zmin << " / " << zmax );
            // the radialstep
            double zStep = (zmax-zmin)/(layers-1);
            double minR = 0.;
            double maxR = 0.;
 
            Amg::Transform3D navLayerTransform;
            double navigationZ                  = 0.;
            // loop over layers
            for (auto& layIter : discLayers) {
                // get the dimensions
                const Trk::DiscSurface& layerSurface = layIter->surfaceRepresentation();
                double currentZ = layerSurface.center().z();
                // create the navigation Z from current Z
                navigationZ = currentZ - 0.5*(zStep);
                navLayerTransform = Amg::Transform3D(Amg::Translation3D(0.,0.,navigationZ));
                auto navLayerSurface = std::make_unique<Trk::DiscSurface>(navLayerTransform, minR, maxR);
                // push that layer back
                ATH_MSG_VERBOSE( "bi-equidistant : creating disc-like NavigationLayer at z-Position : " << navigationZ );
                const Amg::Vector3D center = navLayerSurface->center();
                layerOrderVector.emplace_back(
                                            std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                            center);
                // and get dimensions
                const Trk::DiscBounds* dbounds = dynamic_cast<const Trk::DiscBounds*>(&(layerSurface.bounds()));
                if (dbounds) {
                    minR = dbounds->rMin();
                    maxR = dbounds->rMax();
                } else { // protection
                    minR = 0.;
                    maxR = 100000.;
                }
                // get the material layer first
                ATH_MSG_VERBOSE( "bi-equidistant : registering disc-like MaterialLayer at z-Position : " << currentZ );
                layerOrderVector.emplace_back(std::shared_ptr<Trk::Layer>(layIter),
                                              layerSurface.center());
            }
            // special treatment for last bin
            ATH_MSG_VERBOSE( "bi-equidistant : creating disc-like NavigationLayer at z-Position : " << navigationZ + zStep );
            navLayerTransform = Amg::Transform3D(Amg::Translation3D(0.,0.,navigationZ+zStep));
            auto navLayerSurface = std::make_unique<Trk::DiscSurface>(navLayerTransform, minR, maxR);
            const Amg::Vector3D center = navLayerSurface->center();
            layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                          center);
            // verbose output
            ATH_MSG_VERBOSE( layerOrderVector.size() << " disc Layers (material + navigation) built. " );
 
            // create the binUtility
            double zminBoundary = zmin-zStep+0.5*layerThickness;
            double zmaxBoundary = zmax+zStep+0.5*layerThickness;
            auto binUtility = Trk::BinUtility(layers, layerThickness, zminBoundary, zmaxBoundary, Trk::open, Trk::binZ);
            ATH_MSG_VERBOSE( "bi-equidistant : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            discLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // arbitrary binning
        case Trk::arbitrary :
        {
            std::vector<float>  boundaries;
 
            // max disc dimension
            double minR = 10e10;
            double maxR = 0.;
 
            // initial boundary
            boundaries.push_back(zmin);
 
            // loop over the provided layers and put NavigationLayers in between
            for (auto& layIter :  discLayers ) {
                // get the cylinder surface
                const Trk::DiscSurface& layerSurface = layIter->surfaceRepresentation();
                // and get dimensions
                const Trk::DiscBounds* dbounds = dynamic_cast<const Trk::DiscBounds*>(&(layerSurface.bounds()));
                if (dbounds) {
                    double layInnerR = dbounds->rMin();
                    double layOuterR = dbounds->rMax();
                    minR = (layInnerR < minR) ? layInnerR : minR;
                    maxR = (layOuterR > maxR) ? layOuterR : maxR;
                } else { // protection
                    minR = 0.;
                    maxR = 100000.;
                }
                // the radius & position
                double layerPositionZ = layerSurface.center().z();
                double layerThickness = layIter->thickness();
                // navigation Layer
                double navLayerPositionZ = 0.5*((layerPositionZ-0.5*layerThickness)+boundaries[boundaries.size()-1]);
                // now fill the layer post slot after navigation layer has been determined
                // the transform for this
                Amg::Transform3D navLayerTransform = Amg::Transform3D(Amg::Translation3D(0.,0.,navLayerPositionZ));
                auto navLayerSurface = std::make_unique<Trk::DiscSurface>(navLayerTransform, minR, maxR);
 
                // the material layer
                discLayer = checkAndReplaceEmptyLayer(layIter);
                if (discLayer) {
                    ATH_MSG_VERBOSE( "arbitrary : creating disc-like NavigationLayer at z-Position : " << navLayerPositionZ );
                    layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                                  Amg::Vector3D(0., 0., navLayerPositionZ));
                    ATH_MSG_VERBOSE( "arbitrary : registering disc-like MaterialLayer at z-Position : " << layerPositionZ );
                    layerOrderVector.emplace_back(
                                                std::shared_ptr<Trk::Layer>(discLayer),
                                                Amg::Vector3D(0.,0., layerPositionZ));
                    boundaries.push_back(layerPositionZ-0.5*layerThickness);
                    boundaries.push_back(layerPositionZ+0.5*layerThickness);
                } else {
                    ATH_MSG_VERBOSE( "arbitrary : empty layer configuration cancelled this building of navigation layer.");
                }
            }
            // final material layer
            double navLayerPositionZFinal = 0.5*(zmax+boundaries[boundaries.size()-1]);
            Amg::Transform3D navLayerTransformFinal = Amg::Transform3D(
              Amg::Translation3D(0., 0., navLayerPositionZFinal));
            auto navLayerSurfaceFinal = std::make_unique<Trk::DiscSurface>(navLayerTransformFinal, minR, maxR);
            ATH_MSG_VERBOSE( "arbitrary : creating disc-like NavigationLayer at z-Position : " << navLayerPositionZFinal );
            layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurfaceFinal)),
                                          Amg::Vector3D(0., 0., navLayerPositionZFinal));
            ATH_MSG_VERBOSE( layerOrderVector.size() << " disc Layers (material + navigation) built. " );
            // register the last bounday
            boundaries.push_back(zmax);
            // create the BinUtility
            auto binUtility = Trk::BinUtility(boundaries, Trk::open, Trk::binZ);
            ATH_MSG_VERBOSE( "arbitrary : created a BinUtility as " << binUtility );
 
            // create the BinnedArray; BinnedArray now owns the binUtility pointer
            // cppcheck-suppress memleak
            discLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // default return 0
        default : { return nullptr; }
    }
 
    return discLayerArray;
}
 
 
std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>> Trk::LayerArrayCreator::planeLayerArray(
  const std::vector<Trk::PlaneLayer*>& planeLayersInput,
  double posmin,
  double posmax,
  Trk::BinningType btype,
  Trk::BinningValue bv) const
{
    ATH_MSG_VERBOSE( " build LayerArray with " << planeLayersInput.size() << " plane-like material layers." );
 
    // needed for all cases
    std::unique_ptr<Trk::BinnedArray1D<Trk::Layer>>  planeLayerArray = nullptr;
    std::vector< std::pair< std::shared_ptr<Trk::Layer>, Amg::Vector3D> >   layerOrderVector;
    Amg::Vector3D layerCenter(0.,0.,0.);
 
    //copy so that you can sort
    std::vector<Trk::PlaneLayer*> planeLayers(planeLayersInput);
 
    auto sortBegin = planeLayers.begin();
    auto sortEnd   = planeLayers.end();
    switch  (bv) {
        case Trk::binX : { std::sort(sortBegin, sortEnd, Trk::PlaneLayerSorterX()); } break;
        case Trk::binY : { std::sort(sortBegin, sortEnd, Trk::PlaneLayerSorterY()); } break;
        case Trk::binZ : { std::sort(sortBegin, sortEnd, Trk::PlaneLayerSorterZ()); } break;
        default : {
            ATH_MSG_WARNING("Plane Layers can only be sorted in x/y/z. Returning 0.");
            return nullptr;
        }
    }
 
    // the iterator
    auto layIter = planeLayers.begin();
 
    switch (btype) {
 
        // equidistant binning
        case Trk::equidistant :
        {
            // count the layers
            unsigned int layers = planeLayers.size();
            // loop over layers and put them in
            for ( ; layIter != planeLayers.end(); ++layIter) {
                // get the X
                const Trk::PlaneSurface& layerSurface = (*layIter)->surfaceRepresentation();
                ATH_MSG_VERBOSE( "equidistant : registering plane-like MaterialLayer   at position : " << layerSurface.center() );
 
                layerOrderVector.emplace_back(
                                            std::shared_ptr<Layer>(*layIter),
                                            layerSurface.center());
            }
            // create the binUitlity
            auto binUtility = Trk::BinUtility(layers,posmin,posmax, Trk::open, bv);
            // create the BinnedArray
            planeLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // bi-equidistant binning
        case Trk::biequidistant :
        {
            // count the layers
            unsigned int layers = planeLayers.size();
            // the x-step
            double posStep = (posmax-posmin)/(layers+1);
 
            double currentPos = posmin + posStep;
            double lastPos    = posmin;
 
            double minHalfX = 0.;
            double maxHalfX = 0.;
            double halfY    = 0.;
 
            double layerThickness = 0.;
 
            // loop over layers
            for ( ; layIter != planeLayers.end() ; ++layIter) {
 
                // get the dimensions
                const Trk::PlaneSurface& layerSurface = (*layIter)->surfaceRepresentation();
                // get dimensions
                const Trk::RectangleBounds* recbounds = dynamic_cast<const Trk::RectangleBounds*>(&(layerSurface.bounds()));
                // try rectangular hypothesis
                if (recbounds) {
                    maxHalfX = recbounds->halflengthX();
                    halfY    = recbounds->halflengthY();
                } else {
                    // try trapezoidal hypothesis
                    const Trk::TrapezoidBounds* trapbounds = dynamic_cast<const Trk::TrapezoidBounds*>(&(layerSurface.bounds()));
                    if (trapbounds) {
                        minHalfX = trapbounds->minHalflengthX();
                        maxHalfX = trapbounds->maxHalflengthX();
                        halfY    = trapbounds->halflengthY();
                    } else {
                        // protection
                        minHalfX = 0.;
                        maxHalfX = 10e10;
                        halfY    = 10e10;
                    }
                }
 
                layerThickness = ((*layIter)->thickness() > layerThickness ) ? (*layIter)->thickness() : layerThickness;
 
                // the navigation Layer
                double navigationPos = 0.5*(currentPos+lastPos);
                double navigationX = (bv == Trk::binX) ? navigationPos : 0.;
                double navigationY = (bv == Trk::binY) ? navigationPos : 0.;
                double navigationZ = (bv == Trk::binZ) ? navigationPos : 0.;
                Amg::Translation3D(navigationX,navigationY,navigationZ);
 
                std::unique_ptr<Trk::PlaneSurface> navLayerSurface = nullptr;
                Amg::Transform3D navLayerTransform(Amg::Translation3D(navigationX,0.,0.));
 
                if (std::abs(minHalfX)<10e-5) {
                    navLayerSurface = std::make_unique<Trk::PlaneSurface>(navLayerTransform,
                                                                          maxHalfX,
                                                                          halfY);
                } else {
                    navLayerSurface = std::make_unique<Trk::PlaneSurface>(navLayerTransform,
                                                                          minHalfX,
                                                                          maxHalfX,
                                                                          halfY);
                }
 
                ATH_MSG_VERBOSE( "bi-equidistant : creating plane-like NavigationLayer at position : " << navigationX );
 
                layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                            Amg::Vector3D(navigationX, 0.,0.));
                // restore
                lastPos = currentPos;
                // the material Layer
 
                ATH_MSG_VERBOSE( "bi-equidistant : registering plane-like MaterialLayer at position : " << currentPos );
                layerOrderVector.emplace_back(
                                            std::shared_ptr<Trk::Layer>(*layIter),
                                            layerSurface.center());
 
                // increase the Step
                currentPos += posStep;
            }
 
            // the final navigation layer
            double navigationPosFinal = 0.5*(currentPos+lastPos);
            double navigationXFinal   = (bv == Trk::binX) ? navigationPosFinal : 0.;
            double navigationYFinal   = (bv == Trk::binY) ? navigationPosFinal : 0.;
            double navigationZFinal   = (bv == Trk::binZ) ? navigationPosFinal : 0.;
 
            Amg::Transform3D navLayerTransform(Amg::Translation3D(navigationXFinal,navigationYFinal,navigationZFinal));
 
            auto navLayerSurface = (std::abs(minHalfX)<10e-5) ?
                    std::make_unique<Trk::PlaneSurface>(navLayerTransform, maxHalfX,halfY) :
                    std::make_unique<Trk::PlaneSurface>(navLayerTransform, minHalfX, maxHalfX, halfY);
 
            ATH_MSG_VERBOSE( "bi-equidistant : creating plane-like NavigationLayer at position : " << navLayerSurface->center() );
            const Amg::Vector3D center = navLayerSurface->center();
            layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                          center);
 
            // create the binUtility
            auto binUtility = Trk::BinUtility(layers, layerThickness, posmin, posmax, Trk::open, bv);
 
            // create the BinnedArray
            planeLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
 
        // arbitrary binning
        case Trk::arbitrary :
        {
 
            std::vector<float> boundaries;
            boundaries.push_back(posmin);
            // max plane dimension
            double minHalfX = 0.;
            double maxHalfX = 0.;
            double halfY    = 0.;
 
            // loop over the layers and register navigation layers in between
            for ( ; layIter != planeLayers.end(); ++layIter) {
 
                // get the cylinder surface
                const Trk::PlaneSurface& layerSurface = (*layIter)->surfaceRepresentation();
                // get dimensions
                const Trk::RectangleBounds* recbounds = dynamic_cast<const Trk::RectangleBounds*>(&(layerSurface.bounds()));
                // try rectangular hypothesis
                if (recbounds) {
                    maxHalfX = recbounds->halflengthX();
                    halfY    = recbounds->halflengthY();
                } else {
                    // try trapezoidal hypothesis
                    const Trk::TrapezoidBounds* trapbounds = dynamic_cast<const Trk::TrapezoidBounds*>(&(layerSurface.bounds()));
                    if (trapbounds) {
                        minHalfX = trapbounds->minHalflengthX();
                        maxHalfX = trapbounds->maxHalflengthX();
                        halfY    = trapbounds->halflengthY();
                    } else {
                        // protection
                        minHalfX = 0.;
                        maxHalfX = 10e10;
                        halfY    = 10e10;
                    }
                }
                // the x position
                // and get dimensions
                layerCenter = layerSurface.center();
                double layerPosition = layerCenter[bv];
                // get the thickness
                double layerThickness = (*layIter)->thickness();
                // register
                boundaries.push_back(layerPosition-0.5*layerThickness);
                double navLayerPositionX = (bv == Trk::binX) ? 0.5*(layerPosition+boundaries[boundaries.size()-1]) : layerCenter.x();
                double navLayerPositionY = (bv == Trk::binY) ? 0.5*(layerPosition+boundaries[boundaries.size()-1]) : layerCenter.y();
                double navLayerPositionZ = (bv == Trk::binZ) ? 0.5*(layerPosition+boundaries[boundaries.size()-1]) : layerCenter.z();
                Amg::Translation3D navLayerPosition(navLayerPositionX,navLayerPositionY,navLayerPositionZ);
                Amg::Transform3D navLayerTransform(navLayerPosition);
                // create the navigation plane layer
                auto navLayerSurface = (std::abs(minHalfX)<10e-5) ?
                        std::make_unique<Trk::PlaneSurface>( navLayerTransform, maxHalfX, halfY ) :
                        std::make_unique<Trk::PlaneSurface>( navLayerTransform, minHalfX, maxHalfX, halfY );
                ATH_MSG_VERBOSE( "arbitrary : creating plane-like NavigationLayer at position : " << navLayerPositionX );
                layerOrderVector.emplace_back(
                                            std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurface)),
                                            Amg::Vector3D(navLayerPositionX, navLayerPositionY, navLayerPositionZ));
                // register the material layer
                boundaries.push_back(layerPosition+0.5*layerThickness);
                // material layer
                layerOrderVector.emplace_back(
                                            std::shared_ptr<Trk::Layer>(*layIter),
                                            layerSurface.center());
 
            }
            // last NavigationLayer
            double navLayerPositionXFinal = (bv == Trk::binX) ? 0.5*(posmax+boundaries[boundaries.size()-1]) : layerCenter.x();
            double navLayerPositionYFinal = (bv == Trk::binY) ? 0.5*(posmax+boundaries[boundaries.size()-1]) : layerCenter.y();
            double navLayerPositionZFinal = (bv == Trk::binZ) ? 0.5*(posmax+boundaries[boundaries.size()-1]) : layerCenter.z();
            Amg::Translation3D navLayerPositionFinal(navLayerPositionXFinal,navLayerPositionYFinal,navLayerPositionZFinal);
            Amg::Transform3D navLayerTransformFinal(navLayerPositionFinal);
            // create the navigation plane layer
            auto navLayerSurfaceFinal = (std::abs(minHalfX)<10e-5) ?
                        std::make_unique<Trk::PlaneSurface>( navLayerTransformFinal, maxHalfX, halfY ) :
                        std::make_unique<Trk::PlaneSurface>( navLayerTransformFinal, minHalfX, maxHalfX, halfY );
            ATH_MSG_VERBOSE( "arbitrary : creating plane-like NavigationLayer at position : " << 0.5*(posmax+boundaries[boundaries.size()-1]) );
            layerOrderVector.emplace_back(std::make_shared<Trk::NavigationLayer>(std::move(navLayerSurfaceFinal)),
                                          Amg::Vector3D(navLayerPositionXFinal, navLayerPositionYFinal, navLayerPositionZFinal));
 
            ATH_MSG_VERBOSE( layerOrderVector.size() << " plane Layers (material + navigation) built. " );
 
            // create the BinUtility
            auto binUtility = Trk::BinUtility(boundaries, Trk::open, bv);
            // and the BinnedArray
            planeLayerArray = std::make_unique<Trk::BinnedArray1D<Trk::Layer>>(layerOrderVector, binUtility);
 
        } break;
        // default return 0
        default : { return nullptr; }
    }
 
    return planeLayerArray;
    //cppcheck-suppress memleak
}
 
Trk::Layer* Trk::LayerArrayCreator::checkAndReplaceEmptyLayer(Trk::Layer* lay) const {
    // empty layers will be replaced by navigation layers
    if (m_emptyLayerMode){
        if (lay->layerMaterialProperties() || lay->surfaceArray()) return lay;
        ATH_MSG_VERBOSE(
            "         replacing dummyMaterial layer with "
            << (m_emptyLayerMode > 1 ? " nothing" : " NavigationLayer."));
        Trk::NavigationLayer* nLayer =
            m_emptyLayerMode > 1
                ? nullptr
                : new Trk::NavigationLayer(
                      lay->surfaceRepresentation().uniqueClone(), 1.);
        delete lay;
        return nLayer;
    }
    // don't replace - just give back what you had
    return lay;
}