BTLDetLayerGeometryBuilder Class Reference

#include <BTLDetLayerGeometryBuilder.h>

Public Member Functions
	BTLDetLayerGeometryBuilder ()
	Constructor. More...
virtual	~BTLDetLayerGeometryBuilder ()
	Destructor. More...

Static Public Member Functions
static std::vector< DetLayer * >	buildLayers (const MTDGeometry &geo, const MTDTopology &topo)
	Operations. More...

Detailed Description

Build the BTL DetLayers.

Author

L. Gray - FNAL

Definition at line 17 of file BTLDetLayerGeometryBuilder.h.

Constructor & Destructor Documentation

BTLDetLayerGeometryBuilder()

BTLDetLayerGeometryBuilder::BTLDetLayerGeometryBuilder

(

)

Constructor.

Definition at line 20 of file BTLDetLayerGeometryBuilder.cc.

{}

~BTLDetLayerGeometryBuilder()

BTLDetLayerGeometryBuilder::~BTLDetLayerGeometryBuilder ( )

virtual

Destructor.

Definition at line 22 of file BTLDetLayerGeometryBuilder.cc.

{}

Member Function Documentation

buildLayers()

vector< DetLayer * > BTLDetLayerGeometryBuilder::buildLayers ( const MTDGeometry &  geo, const MTDTopology &  topo  )

static

Operations.

Definition at line 24 of file BTLDetLayerGeometryBuilder.cc.

References MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, MTDGeometry::detsBTL(), MTDTopology::getMTDTopologyMode(), BTLDetId::HALF_ROD, MTDGeometry::idToDet(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, BTLDetId::kCrystalTypes, BTLDetId::kModulesPerRODBarPhiFlat, BTLDetId::kModulesPerRUV2, BTLDetId::kModulesPerTrkV2, BTLDetId::kRUPerTypeV2, LogTrace, callgraph::module, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeomDet::position(), precomputed_value_sort(), DetId::rawId(), mps_fire::result, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), BTLDetId::v2, and BTLDetId::v3.

Referenced by MTDDetLayerGeometry::buildLayers().

                                                                                                         {
  vector<DetLayer*> detlayers;
  vector<MTDTrayBarrelLayer*> result;

  vector<const DetRod*> btlDetTrays;

  const int mtdTopologyMode = topo.getMTDTopologyMode();
  BTLDetId::CrysLayout btlL = MTDTopologyMode::crysLayoutFromTopoMode(mtdTopologyMode);

  if (btlL != BTLDetId::CrysLayout::v2 && btlL != BTLDetId::CrysLayout::v3) {
    for (unsigned tray = 1; tray <= BTLDetId::HALF_ROD; ++tray) {
      vector<const GeomDet*> geomDets;
      for (unsigned module = 1; module <= BTLDetId::kModulesPerRODBarPhiFlat; ++module) {
        for (unsigned side = 0; side <= 1; ++side) {
          const GeomDet* geomDet = geo.idToDet(BTLDetId(side, tray, module, 0, 1));
          if (geomDet != nullptr) {
            geomDets.push_back(geomDet);
            LogTrace("MTDDetLayers") << "get BTL module " << std::hex << BTLDetId(side, tray, module, 0, 1).rawId()
                                     << std::dec << " at R=" << geomDet->position().perp()
                                     << ", phi=" << geomDet->position().phi();
          }
        }
      }

      if (!geomDets.empty()) {
        precomputed_value_sort(geomDets.begin(), geomDets.end(), geomsort::DetZ());
        btlDetTrays.push_back(new MTDDetTray(geomDets));
        LogTrace("MTDDetLayers") << "  New BTLDetTray with " << geomDets.size()
                                 << " modules at R=" << btlDetTrays.back()->position().perp()
                                 << ", phi=" << btlDetTrays.back()->position().phi();
      }
    }

    precomputed_value_sort(btlDetTrays.begin(), btlDetTrays.end(), geomsort::ExtractPhi<GeometricSearchDet, float>());
    result.push_back(new MTDTrayBarrelLayer(btlDetTrays));
    LogTrace("MTDDetLayers") << "BTLDetLayerGeometryBuilder: new MTDTrayBarrelLayer with " << btlDetTrays.size()
                             << " rods, at R " << result.back()->specificSurface().radius();
  } else {
    vector<const GeomDet*> geomDets;

    // logical tracking trays are now rows along z of modules, 3 per each mechanical tray, running from -z to z
    // MTDGeometry is already built with the proper ordering, it is enough to exploit that
    const uint32_t modulesPerDetTray =
        2 * BTLDetId::kRUPerTypeV2 * BTLDetId::kCrystalTypes * BTLDetId::kModulesPerRUV2 / BTLDetId::kModulesPerTrkV2;
    geomDets.reserve(modulesPerDetTray);

    uint32_t index(0);
    for (const auto& det : geo.detsBTL()) {
      index++;
      geomDets.emplace_back(det);
      if (index == modulesPerDetTray) {
        btlDetTrays.emplace_back(new MTDDetTray(geomDets));
        LogTrace("MTDDetLayers") << "  New BTLDetTray with " << geomDets.size()
                                 << " modules at R=" << btlDetTrays.back()->position().perp()
                                 << ", phi=" << btlDetTrays.back()->position().phi();
        index = 0;
        geomDets.clear();
      }
    }

    result.emplace_back(new MTDTrayBarrelLayer(btlDetTrays));
    LogTrace("MTDDetLayers") << "BTLDetLayerGeometryBuilder: new MTDTrayBarrelLayer with " << btlDetTrays.size()
                             << " rods, at R " << result.back()->specificSurface().radius();
  }

  for (vector<MTDTrayBarrelLayer*>::const_iterator it = result.begin(); it != result.end(); it++)
    detlayers.push_back((DetLayer*)(*it));

  return detlayers;
}