HGTD_LayerBuilderCond Class Reference

#include <HGTD_LayerBuilderCond.h>

Inheritance diagram for HGTD_LayerBuilderCond:

Collaboration diagram for HGTD_LayerBuilderCond:

Public Member Functions
	HGTD_LayerBuilderCond (const std::string &, const std::string &, const IInterface *)
	AlgTool style constructor. More...
virtual	~HGTD_LayerBuilderCond ()
	Destructor. More...
virtual StatusCode	initialize () override
	AlgTool initialize method. More...
virtual std::unique_ptr< const std::vector< Trk::CylinderLayer * > >	cylindricalLayers (const EventContext &ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &whandle) const override final
	LayerBuilder interface method - returning Barrel-like layers. More...
virtual std::unique_ptr< const std::vector< Trk::DiscLayer * > >	discLayers (const EventContext &ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &whandle) const override final
	LayerBuilder interface method - returning Endcap-like layers. More...
virtual std::unique_ptr< const std::vector< Trk::PlaneLayer * > >	planarLayers (const EventContext &ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &whandle) const override final
	LayerBuilder interface method - returning Planar-like layers. More...
virtual const std::string &	identification () const override final
	Name identification. More...

Private Member Functions
SG::ReadCondHandle< InDetDD::HGTD_DetectorElementCollection >	retrieveHGTDdetElements (const EventContext &ctx) const
	helper method to construct HGTD materia More...
const Trk::BinnedLayerMaterial	discLayerMaterial (double rMin, double rMax) const
	layer association More...
void	registerSurfacesToLayer (std::span< Trk::Surface *const > &surfaces, const Trk::Layer &layer) const

Static Private Member Functions
static void	evaluateBestBinning (std::vector< Trk::SurfaceOrderPosition > &surfaces, std::vector< float > &rBins, float &maxRadius, std::vector< std::vector< float >> &phiBins)

Private Attributes
const HGTD_DetectorManager *	m_hgtdMgr
	the HGTD Detector Manager More...
const HGTD_ID *	m_hgtdHelper
	HGTD Id Helper. More...
bool	m_setLayerAssociation
	Set Layer Association. More...
std::string	m_identification
	string identification More...
int	m_rBins
	set the number of bins More...
int	m_phiBins
	set the number of bins More...
float	m_discEnvelopeR
	set disc envelope More...
float	m_discThickness
	set disc thickness More...
bool	m_runGeometryValidation
	run geometry validation More...
SG::ReadCondHandleKey< InDetDD::HGTD_DetectorElementCollection >	m_HGTD_ReadKey

Detailed Description

The HGTD_LayerBuilderCond parses the senstive detector elments and orders them onto a Disc surface; no cylindrical layers are expected. This implementation is based on what done in the SiLayerBuilderCond, adapted to the HGTD use case.

Definition at line 51 of file HGTD_LayerBuilderCond.h.

Constructor & Destructor Documentation

HGTD_LayerBuilderCond()

HGTD_LayerBuilderCond::HGTD_LayerBuilderCond	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

AlgTool style constructor.

Definition at line 33 of file HGTD_LayerBuilderCond.cxx.

                                                                                                        :
  base_class(t,n,p),
  m_hgtdMgr(nullptr),
  m_hgtdHelper(nullptr),
  m_setLayerAssociation(true),
  m_identification("HGTD"),
  m_rBins(50),
  m_phiBins(100),
  m_discEnvelopeR(50.),
  m_discThickness(0.2),
  m_runGeometryValidation(true)
{
  // general steering
  declareProperty("SetLayerAssociation"           , m_setLayerAssociation);
  // identification
  declareProperty("Identification"                , m_identification);
  // set some parameters
  declareProperty("DiscMaterialBinsR"             , m_rBins);
  declareProperty("DiscMaterialBinsPhi"           , m_phiBins);
  declareProperty("DiscEnvelope"                  , m_discEnvelopeR);
  declareProperty("DiscThickness"                 , m_discThickness);
  // validation
  declareProperty("RunValidation"                 , m_runGeometryValidation);
}

~HGTD_LayerBuilderCond()

HGTD_LayerBuilderCond::~HGTD_LayerBuilderCond ( )

virtualdefault

Destructor.

Member Function Documentation

cylindricalLayers()

std::unique_ptr< const std::vector< Trk::CylinderLayer * > > HGTD_LayerBuilderCond::cylindricalLayers ( const EventContext &  ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &  whandle  ) const

inlinefinaloverridevirtual

LayerBuilder interface method - returning Barrel-like layers.

Definition at line 120 of file HGTD_LayerBuilderCond.h.

{
  return nullptr;
}

discLayerMaterial()

const Trk::BinnedLayerMaterial HGTD_LayerBuilderCond::discLayerMaterial ( double  rMin, double  rMax  ) const

private

layer association

Definition at line 327 of file HGTD_LayerBuilderCond.cxx.

{
  Trk::BinUtility layerBinUtilityR(m_rBins, rMin, rMax, Trk::open, Trk::binR);
  Trk::BinUtility layerBinUtilityPhi(m_phiBins, -M_PI, M_PI, Trk::closed, Trk::binPhi);
  layerBinUtilityR += layerBinUtilityPhi;
  return Trk::BinnedLayerMaterial(layerBinUtilityR);
}

discLayers()

std::unique_ptr< const std::vector< Trk::DiscLayer * > > HGTD_LayerBuilderCond::discLayers ( const EventContext &  ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &  whandle  ) const

finaloverridevirtual

LayerBuilder interface method - returning Endcap-like layers.

position & bounds of the active Layer

Definition at line 90 of file HGTD_LayerBuilderCond.cxx.

{
  ATH_MSG_DEBUG( "calling HGTD_LayerBuilderCond::discLayers()" );
 
  // sanity check for HGTD Helper
  if (!m_hgtdHelper){
    ATH_MSG_ERROR("HGTD Detector Manager or ID Helper could not be retrieved - giving up.");
    return nullptr;
  }
 
  // get general layout
  SG::ReadCondHandle<InDetDD::HGTD_DetectorElementCollection> readHandle = retrieveHGTDdetElements(ctx);
  if(*readHandle == nullptr){
    return nullptr;
  }
  whandle.addDependency (readHandle);
 
  const InDetDD::HGTD_DetectorElementCollection* readCdo{*readHandle};
  InDetDD::HGTD_DetectorElementCollection::const_iterator hgtdDetIter = readCdo->begin();
 
  // loop on all modules (selecting only one endcap side)
  // and evaluates the number of discs
  // assuming you have the same number of modules on both sides
  int nlayers = 0;
  for (; hgtdDetIter != readCdo->end(); ++hgtdDetIter){
    Identifier    currentId((*hgtdDetIter)->identify());
    // skipping negative side
    if (m_hgtdHelper->endcap(currentId)<0) continue;
    if (m_hgtdHelper->layer(currentId)>nlayers)
      nlayers++;
  }
  // adding one layer offset
  nlayers+=1;
 
  ATH_MSG_DEBUG( "Configured to build " << nlayers << " *2 disc-like layers  (+ additional support layers)." );
 
  // prepare the vectors
  std::vector<float>                                     discZpos(2*nlayers,0.);
  std::vector< std::vector<Trk::SurfaceOrderPosition> >  discSurfaces(2*nlayers, std::vector<Trk::SurfaceOrderPosition>());
 
  int hgtdModules = 0;
  int sumCheckhgtdModules = 0;
  unsigned int currentlayer = 0;
  float maxRmax = -std::numeric_limits<float>::max();
  float minRmin =  std::numeric_limits<float>::max();
 
  // get the missing dimensions by loop over DetElements
  hgtdDetIter = readCdo->begin();
  for (; hgtdDetIter != readCdo->end(); ++hgtdDetIter){
  // take it - if
  // a) you have a detector element ... protection
  if ( (*hgtdDetIter) ) {
     // get the identifier
     Identifier    currentId((*hgtdDetIter)->identify());
 
     ATH_MSG_DEBUG("Element : " << m_hgtdHelper->endcap(currentId) << "/"
                                << m_hgtdHelper->layer(currentId) << "/"
                                << m_hgtdHelper->eta_module(currentId) << "/"
                                << m_hgtdHelper->phi_module(currentId));
 
     // increase the counter of HGTD modules
     hgtdModules++;
 
     // parse all z positions for the mean value of the discs
     float currentZ = (*hgtdDetIter)->center().z();
     //calculate current layer and current disk from it
     currentlayer  = m_hgtdHelper->layer(currentId);
     // adding the numbe of layers per side as offset
     currentlayer += currentZ > 0. ? nlayers : 0;
     ATH_MSG_DEBUG( " ----  layer: " << currentlayer );
 
     // evaluate the z-position per layer
     // all modules on the same HGTD layer have the same z
     discZpos[currentlayer] = currentZ;
 
     // evaluate the r-extension per layer
     float currentRmin = (*hgtdDetIter)->rMin();
     float currentRmax = (*hgtdDetIter)->rMax();
     ATH_MSG_DEBUG( " ----  rmin/rmax: " << currentRmin << "/" << currentRmax );
     if (maxRmax<currentRmax)
       maxRmax = currentRmax;
     if (minRmin>currentRmin)
       minRmin = currentRmin;
 
     // fill the elements for the layers into the surface arrays
     // get the center position
     const Amg::Vector3D& orderPosition = (*hgtdDetIter)->center();
 
 
     // Register the chosen side in the object array
     //
     // Passing a no-op deleter (no delete happens)
     // This line can be problematic
     // 1. We couple the DetElement owned surface to the Tracking Geometry.
     // The lifetime is not controlled by the geometry.
     // For now we need to be careful on how we schedule these as we do not
     // want to end up with dangling ptr.
     // 2. We modiy the payload (see const_cast).
     // NOTE !!! We can avoid these if we can clone and not "share".
     // Or if we have also created the layer and attached them
     // to the surfaces when the Det Elements are created.
     // Aka run this Layer Builder as part of the DetElement creation.
     Trk::Surface* mutableSurace = const_cast<Trk::Surface*>(&((*hgtdDetIter)->surface()));
     std::shared_ptr<Trk::Surface> sharedSurface(mutableSurace,[](Trk::Surface*) {});
     //
     Trk::SurfaceOrderPosition surfaceOrder(sharedSurface, orderPosition);
 
     discSurfaces[currentlayer].push_back(surfaceOrder);
 
   } else if (!(*hgtdDetIter))
       ATH_MSG_WARNING("Not valid pointer to HGTD Detector element... something wrong with the Id dictionary?");
  }
 
  // adding some envelope
  maxRmax += m_discEnvelopeR;
  minRmin -= m_discEnvelopeR;
 
  // construct the layers
  auto discLayers = std::make_unique<std::vector<Trk::DiscLayer*> >();
 
  double thickness = m_discThickness;
 
  int discCounter = 0;
  for (auto& thisDiscZpos : discZpos) {
    // screen output
    ATH_MSG_DEBUG( "Building a DiscLayer: " );
    ATH_MSG_DEBUG( "  -> At Z - Position      :  " << thisDiscZpos );
    ATH_MSG_DEBUG( "  -> With Thickness       :  " << thickness   << " i- ncludes envelope tolerance : " << m_discEnvelopeR );
    ATH_MSG_DEBUG( "  -> With Rmin/Rmax (est) :  " << minRmin << " / " << maxRmax );
 
    ATH_MSG_DEBUG( "... creating binned array ... ");
 
    std::vector<float> rBins = {minRmin};
    std::vector<std::vector<float>> phiBins = {{}};
 
    evaluateBestBinning(discSurfaces[discCounter], rBins, maxRmax, phiBins);
 
    // Build the BinUtilities using the bins defined at construction
    // the extension is provided in the previous loop
    auto BinUtilityR = Trk::BinUtility(rBins, Trk::open, Trk::binR);
    auto subBinUtilitiesPhi =  std::vector<Trk::BinUtility>();
    ATH_MSG_DEBUG("BinUtilityR --> " << BinUtilityR );
 
    for (unsigned int bin = 0; bin < rBins.size()-1; bin++) {
      auto BinUtilityY = Trk::BinUtility(phiBins.at(bin), Trk::closed, Trk::binPhi);
      subBinUtilitiesPhi.push_back(BinUtilityY);
      ATH_MSG_DEBUG(bin << ") BinUtilityPhi --> " << BinUtilityY );
    }
 
    // prepare the binned array, it can be with one to several rings
    auto currentBinnedArray =
        std::make_unique<Trk::BinnedArray1D1D<Trk::Surface>>(
            discSurfaces[discCounter], BinUtilityR, subBinUtilitiesPhi);
 
    ATH_MSG_DEBUG( "... done!" );
 
    int discSurfacesNum = (discSurfaces[discCounter]).size();
 
    ATH_MSG_DEBUG( "Constructed BinnedArray for DiscLayer with "<< discSurfacesNum << " SubSurfaces." );
 
    // always run the geometry validation to catch flaws
    if (m_runGeometryValidation) {
      // checking for :
      //   - empty surface bins
      //   - doubly filled bins
      std::map< const Trk::Surface*,Amg::Vector3D > uniqueSurfaceMap;
      std::map< const Trk::Surface*,Amg::Vector3D >::iterator usmIter = uniqueSurfaceMap.end();
      // check the registered surfaces in the binned array
      std::span<Trk::Surface * const> arraySurfaces = currentBinnedArray->arrayObjects();
      size_t dsumCheckSurfaces = 0;
      double lastPhi = 0.;
      for (const auto & asurfIter : arraySurfaces){
        if ( asurfIter ) {
          ++dsumCheckSurfaces;
          usmIter = uniqueSurfaceMap.find(asurfIter);
          lastPhi = asurfIter->center().phi();
          if (usmIter != uniqueSurfaceMap.end()) {
            ATH_MSG_WARNING("Non-unique surface found with eta/phi = "
                            << asurfIter->center().eta() << " / "
                            << asurfIter->center().phi());
          } else {
            uniqueSurfaceMap[asurfIter] = asurfIter->center();
          }
        } else {
          ATH_MSG_WARNING("Zero-pointer in array detected in this ring, last "
                          "valid phi value was = "
                          << lastPhi << " --> discCounter: " << discCounter);
        }
      }
      sumCheckhgtdModules +=  dsumCheckSurfaces;
    }
 
    // get the layer material from the helper method
    const Trk::LayerMaterialProperties& layerMaterial = discLayerMaterial(minRmin,maxRmax);
 
    Amg::Transform3D activeLayerTransform ;
    activeLayerTransform = Amg::Translation3D(0.,0.,thisDiscZpos);
 
    auto activeLayerBounds = std::make_shared<Trk::DiscBounds>(minRmin, maxRmax);
 
    auto olDescriptor = std::make_unique<HGTD_OverlapDescriptor>(
        currentBinnedArray.get(), rBins, phiBins);
 
    // register the layer to the surfaces
    std::span<Trk::Surface * const> layerSurfaces = currentBinnedArray->arrayObjects();
    // layer creation; deletes currentBinnedArray in baseclass 'Layer' upon destruction
    // activeLayerTransform deleted in 'Surface' baseclass
    Trk::DiscLayer* activeLayer = new Trk::DiscLayer(activeLayerTransform,
                                                     activeLayerBounds,
                                                     std::move(currentBinnedArray),
                                                     layerMaterial,
                                                     thickness,
                                                     std::move(olDescriptor));
 
    registerSurfacesToLayer(layerSurfaces,*activeLayer);
    discLayers->push_back(activeLayer);
    // increase the disc counter by one
    ++discCounter;
  }
 
  //
  ATH_MSG_DEBUG( hgtdModules << " HGTD Modules parsed for Disc Layer dimensions." );
  if (m_runGeometryValidation) {
    ATH_MSG_DEBUG( sumCheckhgtdModules << " HGTD Modules filled in Disc Layer Arrays." );
    if ( hgtdModules-sumCheckhgtdModules )
      ATH_MSG_WARNING( hgtdModules-sumCheckhgtdModules << " Modules not registered properly in binned array." );
  }
 
  // sort the vector
  Trk::DiscLayerSorterZ zSorter;
  std::sort(discLayers->begin(), discLayers->end(), zSorter);
 
  return discLayers;
}

evaluateBestBinning()

void HGTD_LayerBuilderCond::evaluateBestBinning ( std::vector< Trk::SurfaceOrderPosition > &  surfaces, std::vector< float > &  rBins, float &  maxRadius, std::vector< std::vector< float >> &  phiBins  )

staticprivate

Definition at line 347 of file HGTD_LayerBuilderCond.cxx.

{
  // get all the centers (r,phi), as you want to play with them
  std::vector < std::pair< float, float> > centers = {};
  centers.reserve(surfaces.size());
  for ( auto& orderedSurface : surfaces) {
    centers.emplace_back(orderedSurface.second.perp(), orderedSurface.second.phi());
  }
 
  // sorting the centers accordingly to r
  std::sort(centers.begin(), centers.end(),
            [](const std::pair< float, float>& a, const std::pair< float, float>& b) -> bool {
              return a.first < b.first;
            });
 
  // at the beginning use a fine binning in phi
  // it is updated later to fit the amount of surfaces
  // once you have defined a bin in radius
  int bins = 100;
  float step = 2*M_PI/float(bins);
  std::vector<float> finerBinning = {};
  finerBinning.reserve(bins);
 
  for (int bin = 0; bin<=bins; bin++) {
    finerBinning.push_back(-M_PI+step*bin);
  }
 
  // use this vector to save the indices and
  // guess when you have to add
  // an additional bin in r
  std::vector<int> phiIndices = {};
  std::vector<float> tmpRadii = {};
 
  for (auto& center : centers) {
    float phi = center.second;
    const auto boundVal =  std::lower_bound(finerBinning.begin(), finerBinning.end(), phi);
    int phiIndex = std::distance(finerBinning.begin(), boundVal);
    // if the element fits in the given r bin, add it,
    // otherwise reset the indices and start a new r bin
    if (std::find(phiIndices.begin(), phiIndices.end(), phiIndex)==phiIndices.end()) {
      phiIndices.push_back(phiIndex);
      tmpRadii.push_back(center.first);
    } else {
      phiIndices.clear();
      for (unsigned int index = (tmpRadii.size()-1); index>0; index--) {
        auto& prevRadius = tmpRadii.at(index);
        if ( std::abs(prevRadius - center.first)<1e-5 ) {
          const auto boundVal =  std::lower_bound(finerBinning.begin(), finerBinning.end(), phi);
          int phiIndex = std::distance(finerBinning.begin(), boundVal);
          phiIndices.push_back(phiIndex);
          continue;
        } else {
          float r = 0.5*(prevRadius+center.first);
          rBins.push_back(r);
          tmpRadii = {prevRadius};
          break;
        }
      }
    }
  }
 
  rBins.push_back(maxRadius);
 
  // now we have the best binning in r and want to
  // map the centers accordingly to this
  std::vector< std::vector < float > > binnedCenters = {{}};
 
  for (auto& center : centers) {
    float r = center.first;
    float phi = center.second;
    const auto boundVal =  std::lower_bound(rBins.begin(), rBins.end(), r);
    int rIndex = std::distance(rBins.begin(), boundVal);
    if (int(binnedCenters.size())<rIndex)
      binnedCenters.push_back({phi});
    else
      binnedCenters.back().push_back(phi);
  }
 
  // now that we have the centers binned in r, we evaluate the best
  // bin in phi for each of those bins
  bool isFirst = true;
  for (auto& centersInBin : binnedCenters) {
    // sorting the centers accordingly to phi_bins
    std::sort(centersInBin.begin(), centersInBin.end());
    if (isFirst) {
      phiBins.back().push_back(-M_PI);
      isFirst=false;
    } else phiBins.push_back({-M_PI});
    for (unsigned int index = 0; index<(centersInBin.size()-1); index++) {
      float phi = 0.5*(centersInBin.at(index)+centersInBin.at(index+1));
      phiBins.back().push_back(phi);
    }
  }
 
  }

identification()

const std::string & HGTD_LayerBuilderCond::identification ( ) const

inlinefinaloverridevirtual

Name identification.

Definition at line 133 of file HGTD_LayerBuilderCond.h.

{ return m_identification; }

initialize()

StatusCode HGTD_LayerBuilderCond::initialize ( )

overridevirtual

AlgTool initialize method.

Definition at line 63 of file HGTD_LayerBuilderCond.cxx.

{
 
    ATH_MSG_DEBUG( "initialize()" );
    // get HGTD Detector Description Manager and HGTD Helper
    ATH_CHECK(detStore()->retrieve(m_hgtdMgr, "HGTD"));
    ATH_CHECK(detStore()->retrieve(m_hgtdHelper, "HGTD_ID"));
 
    // get HGTD detector element collection
    ATH_CHECK(m_HGTD_ReadKey.initialize());
 
    return StatusCode::SUCCESS;
}

planarLayers()

std::unique_ptr< const std::vector< Trk::PlaneLayer * > > HGTD_LayerBuilderCond::planarLayers ( const EventContext &  ctx, SG::WriteCondHandle< Trk::TrackingGeometry > &  whandle  ) const

inlinefinaloverridevirtual

LayerBuilder interface method - returning Planar-like layers.

Definition at line 127 of file HGTD_LayerBuilderCond.h.

{
  return nullptr;
}

registerSurfacesToLayer()

void HGTD_LayerBuilderCond::registerSurfacesToLayer ( std::span< Trk::Surface *const > &  surfaces, const Trk::Layer &  layer  ) const

private

Definition at line 335 of file HGTD_LayerBuilderCond.cxx.

{
   if (!m_setLayerAssociation) return;
   // register the surfaces to the layer
   for (const auto & surfaces : layerSurfaces) {
     if (surfaces) {
       // register the current surfaces
       (*surfaces).associateLayer(lay);
     }
   }
}

retrieveHGTDdetElements()

SG::ReadCondHandle< InDetDD::HGTD_DetectorElementCollection > HGTD_LayerBuilderCond::retrieveHGTDdetElements ( const EventContext &  ctx ) const

private

helper method to construct HGTD materia

Definition at line 78 of file HGTD_LayerBuilderCond.cxx.

{
  auto readHandle = SG::ReadCondHandle<InDetDD::HGTD_DetectorElementCollection> (m_HGTD_ReadKey, ctx);
  if (*readHandle==nullptr) {
    ATH_MSG_ERROR("Null pointer to the read conditions object of " << m_HGTD_ReadKey.key());
  }
  return readHandle;
}

Member Data Documentation

m_discEnvelopeR

float HGTD_LayerBuilderCond::m_discEnvelopeR

private

set disc envelope

Definition at line 105 of file HGTD_LayerBuilderCond.h.

m_discThickness

float HGTD_LayerBuilderCond::m_discThickness

private

set disc thickness

Definition at line 106 of file HGTD_LayerBuilderCond.h.

m_HGTD_ReadKey

SG::ReadCondHandleKey<InDetDD::HGTD_DetectorElementCollection> HGTD_LayerBuilderCond::m_HGTD_ReadKey

private

Initial value:

{
        this,
        "HGTD_ReadKey",
        "HGTD_DetectorElementCollection",
        "Key of output HGTD_DetectorElementCollection for HGTD"
      }

Definition at line 111 of file HGTD_LayerBuilderCond.h.

m_hgtdHelper

const HGTD_ID* HGTD_LayerBuilderCond::m_hgtdHelper

private

HGTD Id Helper.

Definition at line 96 of file HGTD_LayerBuilderCond.h.

m_hgtdMgr

const HGTD_DetectorManager* HGTD_LayerBuilderCond::m_hgtdMgr

private

the HGTD Detector Manager

Definition at line 95 of file HGTD_LayerBuilderCond.h.

m_identification

std::string HGTD_LayerBuilderCond::m_identification

private

string identification

Definition at line 100 of file HGTD_LayerBuilderCond.h.

m_phiBins

int HGTD_LayerBuilderCond::m_phiBins

private

set the number of bins

Definition at line 103 of file HGTD_LayerBuilderCond.h.

m_rBins

int HGTD_LayerBuilderCond::m_rBins

private

set the number of bins

Definition at line 102 of file HGTD_LayerBuilderCond.h.

m_runGeometryValidation

bool HGTD_LayerBuilderCond::m_runGeometryValidation

private

run geometry validation

Definition at line 108 of file HGTD_LayerBuilderCond.h.

m_setLayerAssociation

bool HGTD_LayerBuilderCond::m_setLayerAssociation

private

Set Layer Association.

Definition at line 98 of file HGTD_LayerBuilderCond.h.

---

The documentation for this class was generated from the following files:

• HGTD_LayerBuilderCond.h
• HGTD_LayerBuilderCond.cxx