ETLDetLayerGeometryBuilder Class Reference

#include <ETLDetLayerGeometryBuilder.h>

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

Private Member Functions
	ETLDetLayerGeometryBuilder ()

Static Private Member Functions
static MTDRingForwardDoubleLayer *	buildLayer (int endcap, int layer, std::vector< unsigned > &rings, const MTDGeometry &geo)
static MTDSectorForwardDoubleLayer *	buildLayerNew (int endcap, int layer, std::vector< unsigned > &sectors, const MTDGeometry &geo, const MTDTopology &topo)
static bool	isFront (int layer, int ring, int module)
static MTDDetRing *	makeDetRing (std::vector< const GeomDet * > &geomDets)
static MTDDetSector *	makeDetSector (std::vector< const GeomDet * > &geomDets, const MTDTopology &topo)

Detailed Description

Build the ETL DetLayers.

Author

L. Gray - FNAL

Definition at line 21 of file ETLDetLayerGeometryBuilder.h.

Constructor & Destructor Documentation

ETLDetLayerGeometryBuilder::ETLDetLayerGeometryBuilder ( )

inlineprivate

Definition at line 30 of file ETLDetLayerGeometryBuilder.h.

{}

Member Function Documentation

MTDRingForwardDoubleLayer * ETLDetLayerGeometryBuilder::buildLayer ( int  endcap, int  layer, std::vector< unsigned > &  rings, const MTDGeometry &  geo  )

staticprivate

Definition at line 77 of file ETLDetLayerGeometryBuilder.cc.

References funct::abs(), cms::cuda::assert(), TauDecayModes::dec, MTDGeometry::idToDet(), ETLDetId::kETLmoduleMask, LogTrace, callgraph::module, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeomDet::position(), DetId::rawId(), mps_fire::result, relativeConstraints::ring, and PV3DBase< T, PVType, FrameType >::z().

                                                                                          {
  MTDRingForwardDoubleLayer* result = nullptr;

  vector<const ForwardDetRing*> frontRings, backRings;

  for (unsigned ring : rings) {
    vector<const GeomDet*> frontGeomDets, backGeomDets;
    for (unsigned module = 1; module <= ETLDetId::kETLmoduleMask; ++module) {
      ETLDetId detId(endcap, ring, module, 0);
      const GeomDet* geomDet = geo.idToDet(detId);
      // we sometimes loop over more chambers than there are in ring
      bool isInFront = isFront(layer, ring, module);
      if (geomDet != nullptr) {
        if (isInFront) {
          frontGeomDets.push_back(geomDet);
        } else {
          backGeomDets.push_back(geomDet);
        }
        LogTrace("MTDDetLayers") << "get ETL module " << std::hex << ETLDetId(endcap, layer, ring, module).rawId()
                                 << std::dec << " at R=" << geomDet->position().perp()
                                 << ", phi=" << geomDet->position().phi() << ", z= " << geomDet->position().z()
                                 << " isFront? " << isInFront;
      }
    }

    if (!backGeomDets.empty()) {
      backRings.push_back(makeDetRing(backGeomDets));
    }

    if (!frontGeomDets.empty()) {
      frontRings.push_back(makeDetRing(frontGeomDets));
      assert(!backGeomDets.empty());
      float frontz = frontRings[0]->position().z();
      float backz = backRings[0]->position().z();
      assert(std::abs(frontz) < std::abs(backz));
    }
  }

  // How should they be sorted?
  //    precomputed_value_sort(muDetRods.begin(), muDetRods.end(), geomsort::ExtractZ<GeometricSearchDet,float>());
  result = new MTDRingForwardDoubleLayer(frontRings, backRings);
  LogTrace("MTDDetLayers") << "New MTDRingForwardLayer with " << frontRings.size() << " and " << backRings.size()
                           << " rings, at Z " << result->position().z()
                           << " R1: " << result->specificSurface().innerRadius()
                           << " R2: " << result->specificSurface().outerRadius();
  return result;
}

MTDSectorForwardDoubleLayer * ETLDetLayerGeometryBuilder::buildLayerNew ( int  endcap, int  layer, std::vector< unsigned > &  sectors, const MTDGeometry &  geo, const MTDTopology &  topo  )

staticprivate

Definition at line 140 of file ETLDetLayerGeometryBuilder.cc.

References funct::abs(), cms::cuda::assert(), MTDGeometry::detsETL(), LogDebug, LogTrace, MTDTopology::orderETLSector(), mps_fire::result, and ForwardDetLayer::specificSurface().

                                                                                                       {
  MTDSectorForwardDoubleLayer* result = nullptr;

  std::vector<const MTDDetSector*> frontSectors, backSectors;

  LogDebug("MTDDetLayers") << "ETL dets array size = " << geo.detsETL().size();

  for (unsigned sector : sectors) {
    std::vector<const GeomDet*> frontGeomDets, backGeomDets;
    LogDebug("MTDDetLayers") << "endcap = " << endcap << " layer = " << layer << " sector = " << sector;
#ifdef EDM_ML_DEBUG
    unsigned int nfront(0), nback(0);
#endif
    for (auto det : geo.detsETL()) {
      ETLDetId theMod(det->geographicalId().rawId());
      if (theMod.mtdSide() == endcap && theMod.nDisc() == layer && theMod.sector() == static_cast<int>(sector)) {
        LogTrace("MTDLayerDump") << std::fixed << "ETLDetId " << theMod.rawId() << " side = " << std::setw(4)
                                 << theMod.mtdSide() << " Disc/Side/Sector = " << std::setw(4) << theMod.nDisc() << " "
                                 << std::setw(4) << theMod.discSide() << " " << std::setw(4) << theMod.sector()
                                 << " mod/type = " << std::setw(4) << theMod.module() << " " << std::setw(4)
                                 << theMod.modType() << " pos = " << det->position();
        // front layer face
        if (theMod.discSide() == 0) {
#ifdef EDM_ML_DEBUG
          nfront++;
          LogTrace("MTDDetLayers") << "Front " << theMod.discSide() << " " << nfront;
#endif
          frontGeomDets.emplace_back(det);
          // back layer face
        } else if (theMod.discSide() == 1) {
#ifdef EDM_ML_DEBUG
          nback++;
          LogTrace("MTDDetLayers") << "Back " << theMod.discSide() << " " << nback;
#endif
          backGeomDets.emplace_back(det);
        }
      }
    }

    if (!backGeomDets.empty()) {
      std::sort(backGeomDets.begin(), backGeomDets.end(), topo.orderETLSector);
      LogDebug("MTDDetLayers") << "backGeomDets size = " << backGeomDets.size();
      backSectors.emplace_back(makeDetSector(backGeomDets, topo));
    }

    if (!frontGeomDets.empty()) {
      std::sort(frontGeomDets.begin(), frontGeomDets.end(), topo.orderETLSector);
      LogDebug("MTDDetLayers") << "frontGeomDets size = " << frontGeomDets.size();
      frontSectors.emplace_back(makeDetSector(frontGeomDets, topo));
      assert(!backGeomDets.empty());
      float frontz = frontSectors.back()->position().z();
      float backz = backSectors.back()->position().z();
      assert(std::abs(frontz) < std::abs(backz));
    }
  }

  result = new MTDSectorForwardDoubleLayer(frontSectors, backSectors);
  LogTrace("MTDDetLayers") << "New MTDSectorForwardDoubleLayer with " << std::fixed << std::setw(14)
                           << frontSectors.size() << " and " << std::setw(14) << backSectors.size() << " rings, at Z "
                           << std::setw(14) << result->specificSurface().position().z() << " R1: " << std::setw(14)
                           << result->specificSurface().innerRadius() << " R2: " << std::setw(14)
                           << result->specificSurface().outerRadius();

  return result;
}

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

static

return.first=forward (+Z), return.second=backward (-Z) both vectors are sorted inside-out

Definition at line 21 of file ETLDetLayerGeometryBuilder.cc.

References MTDTopologyMode::barphiflat, MTDTopologyMode::btlv1etlv4, MTDTopologyMode::btlv1etlv5, Reference_intrackfit_cff::endcap, Exception, MTDTopology::getMTDTopologyMode(), ETLDetId::kETLv1maxRing, ETLDetId::kETLv1nDisc, ETLDetId::kETLv4maxSector, ETLDetId::kETLv4nDisc, ETLDetId::kETLv5maxSector, phase1PixelTopology::layer, mps_fire::result, relativeConstraints::ring, and mergeVDriftHistosByStation::sectors.

Referenced by MTDDetLayerGeometry::buildLayers().

                                                                                                             {
  vector<DetLayer*> result[2];  // one for each endcap

  const int mtdTopologyMode = topo.getMTDTopologyMode();
  if (mtdTopologyMode <= static_cast<int>(MTDTopologyMode::Mode::barphiflat)) {
    for (unsigned endcap = 0; endcap < 2; ++endcap) {
      // there is only one layer for ETL right now, maybe more later
      for (unsigned layer = 0; layer < ETLDetId::kETLv1nDisc; ++layer) {
        vector<unsigned> rings;
        rings.reserve(ETLDetId::kETLv1maxRing + 1);
        for (unsigned ring = 1; ring <= ETLDetId::kETLv1maxRing; ++ring) {
          rings.push_back(ring);
        }
        MTDRingForwardDoubleLayer* thelayer = buildLayer(endcap, layer, rings, geo);
        if (thelayer)
          result[endcap].push_back(thelayer);
      }
    }
  } else {
    // number of layers is identical for post TDR scenarios, pick v4
    // loop on number of sectors per face, two faces per disc (i.e. layer) taken into account in layer building (front/back)
    unsigned int nSector(1);
    switch (mtdTopologyMode) {
      case static_cast<int>(MTDTopologyMode::Mode::btlv1etlv4):
        nSector *= ETLDetId::kETLv4maxSector;
        break;
      case static_cast<int>(MTDTopologyMode::Mode::btlv1etlv5):
        nSector *= ETLDetId::kETLv5maxSector;
        break;
      default:
        throw cms::Exception("MTDDetLayers") << "Not implemented scenario " << mtdTopologyMode;
        break;
    }

    for (unsigned endcap = 0; endcap < 2; ++endcap) {
      // number of layers is two, identical for post TDR scenarios, pick v4
      for (unsigned layer = 1; layer <= ETLDetId::kETLv4nDisc; ++layer) {
        vector<unsigned> sectors;
        sectors.reserve(nSector + 1);
        for (unsigned sector = 1; sector <= nSector; ++sector) {
          sectors.push_back(sector);
        }
        MTDSectorForwardDoubleLayer* thelayer = buildLayerNew(endcap, layer, sectors, geo, topo);
        if (thelayer)
          result[endcap].push_back(thelayer);
      }
    }
  }
  //
  // the first entry is Z+ ( MTD side 1), the second is Z- (MTD side 0)
  //
  pair<vector<DetLayer*>, vector<DetLayer*> > res_pair(result[1], result[0]);
  return res_pair;
}

bool ETLDetLayerGeometryBuilder::isFront ( int  layer, int  ring, int  module  )

staticprivate

Definition at line 128 of file ETLDetLayerGeometryBuilder.cc.

{ return (module + 1) % 2; }

MTDDetRing * ETLDetLayerGeometryBuilder::makeDetRing ( std::vector< const GeomDet * > &  geomDets )

staticprivate

Definition at line 130 of file ETLDetLayerGeometryBuilder.cc.

References LogTrace, precomputed_value_sort(), and mps_fire::result.

                                                                                    {
  precomputed_value_sort(geomDets.begin(), geomDets.end(), geomsort::DetPhi());
  MTDDetRing* result = new MTDDetRing(geomDets);
  LogTrace("MTDDetLayers") << "ETLDetLayerGeometryBuilder: new MTDDetRing with " << geomDets.size()
                           << " chambers at z=" << result->position().z()
                           << " R1: " << result->specificSurface().innerRadius()
                           << " R2: " << result->specificSurface().outerRadius();
  return result;
}

MTDDetSector * ETLDetLayerGeometryBuilder::makeDetSector ( std::vector< const GeomDet * > &  geomDets, const MTDTopology &  topo  )

staticprivate

Definition at line 207 of file ETLDetLayerGeometryBuilder.cc.

References LogTrace, and mps_fire::result.

                                                                                                                 {
  MTDDetSector* result = new MTDDetSector(geomDets, topo);
  LogTrace("MTDDetLayers") << "ETLDetLayerGeometryBuilder::makeDetSector new MTDDetSector with " << std::fixed
                           << std::setw(14) << geomDets.size() << " modules \n"
                           << (*result);

  return result;
}