MuonCSCDetLayerGeometryBuilder Class Reference

#include <MuonCSCDetLayerGeometryBuilder.h>

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

Private Member Functions
	MuonCSCDetLayerGeometryBuilder ()

Static Private Member Functions
static MuRingForwardDoubleLayer *	buildLayer (int endcap, int station, std::vector< int > &rings, const CSCGeometry &geo)
static bool	isFront (int station, int ring, int chamber)
static MuDetRing *	makeDetRing (std::vector< const GeomDet *> &geomDets)

Detailed Description

Build the CSC DetLayers.

Author

N. Amapane - CERN

Definition at line 18 of file MuonCSCDetLayerGeometryBuilder.h.

Constructor & Destructor Documentation

MuonCSCDetLayerGeometryBuilder()

MuonCSCDetLayerGeometryBuilder::MuonCSCDetLayerGeometryBuilder ( )

inlineprivate

Definition at line 26 of file MuonCSCDetLayerGeometryBuilder.h.

{}

Member Function Documentation

buildLayer()

MuRingForwardDoubleLayer * MuonCSCDetLayerGeometryBuilder::buildLayer ( int  endcap, int  station, std::vector< int > &  rings, const CSCGeometry &  geo  )

staticprivate

Definition at line 60 of file MuonCSCDetLayerGeometryBuilder.cc.

References cms::cuda::assert(), relativeConstraints::chamber, makeMuonMisalignmentScenario::endcap, CSCGeometry::idToDet(), LogTrace, CSCDetId::maxChamberId(), metname, CSCDetId::minChamberId(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeomDet::position(), mps_fire::result, relativeConstraints::ring, relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, and PV3DBase< T, PVType, FrameType >::z().

                                                                                             {
  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|MuonCSCDetLayerGeometryBuilder";
  MuRingForwardDoubleLayer* result = nullptr;

  vector<const ForwardDetRing*> frontRings, backRings;

  for (vector<int>::iterator ring = rings.begin(); ring != rings.end(); ring++) {
    vector<const GeomDet*> frontGeomDets, backGeomDets;
    for (int chamber = CSCDetId::minChamberId(); chamber <= CSCDetId::maxChamberId(); chamber++) {
      CSCDetId detId(endcap, station, (*ring), chamber, 0);
      const GeomDet* geomDet = geo.idToDet(detId);
      // we sometimes loop over more chambers than there are in ring
      bool isInFront = isFront(station, *ring, chamber);
      if (geomDet != nullptr) {
        if (isInFront) {
          frontGeomDets.push_back(geomDet);
        } else {
          backGeomDets.push_back(geomDet);
        }
        LogTrace(metname) << "get CSC chamber " << CSCDetId(endcap, station, (*ring), chamber, 0)
                          << " 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(fabs(frontz) < fabs(backz));
    }
  }

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

buildLayers()

pair< vector< DetLayer * >, vector< DetLayer * > > MuonCSCDetLayerGeometryBuilder::buildLayers ( const CSCGeometry &  geo )

static

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

Definition at line 17 of file MuonCSCDetLayerGeometryBuilder.cc.

References makeMuonMisalignmentScenario::endcap, mps_fire::i, phase1PixelTopology::layer, CSCDetId::maxRingId(), CSCDetId::maxStationId(), CSCDetId::minRingId(), mps_fire::result, relativeConstraints::ring, and relativeConstraints::station.

Referenced by MuonDetLayerGeometryESProducer::produce().

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

  for (int i = 0; i < 2; i++) {
    int endcap = i + 1;

    // ME/1/1a (= station 1, ring 4) and ME/1/1b (= station 1, ring 1)
    {
      vector<int> rings;
      rings.push_back(4);
      rings.push_back(1);

      MuRingForwardDoubleLayer* layer = buildLayer(endcap, 1, rings, geo);
      if (layer)
        result[i].push_back(layer);
    }

    // ME/1/2 and 1/3 (= station 1, ring 2 and 3)
    {
      vector<int> rings;
      rings.push_back(2);
      rings.push_back(3);

      MuRingForwardDoubleLayer* layer = buildLayer(endcap, 1, rings, geo);
      if (layer)
        result[i].push_back(layer);
    }

    // Stations 2,3,4
    for (int station = 2; station <= CSCDetId::maxStationId(); station++) {
      vector<int> rings;
      for (int ring = CSCDetId::minRingId(); ring <= CSCDetId::maxRingId(); ring++) {
        rings.push_back(ring);
      }
      MuRingForwardDoubleLayer* layer = buildLayer(endcap, station, rings, geo);
      if (layer)
        result[i].push_back(layer);
    }
  }
  pair<vector<DetLayer*>, vector<DetLayer*> > res_pair(result[0], result[1]);
  return res_pair;
}

isFront()

bool MuonCSCDetLayerGeometryBuilder::isFront ( int  station, int  ring, int  chamber  )

staticprivate

Definition at line 110 of file MuonCSCDetLayerGeometryBuilder.cc.

References relativeConstraints::chamber, mps_fire::result, relativeConstraints::ring, and relativeConstraints::station.

                                                                               {
  bool result = false;

  bool isOverlapping = !(station == 1 && ring == 3);
  // not overlapping means back
  if (isOverlapping) {
    bool isEven = (chamber % 2 == 0);
    // odd chambers are bolted to the iron, which faces
    // forward in 1&2, backward in 3&4, so...
    result = (station < 3) ? isEven : !isEven;
  }
  return result;
}

makeDetRing()

MuDetRing * MuonCSCDetLayerGeometryBuilder::makeDetRing ( std::vector< const GeomDet *> &  geomDets )

staticprivate

Definition at line 124 of file MuonCSCDetLayerGeometryBuilder.cc.

References LogTrace, metname, precomputed_value_sort(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|MuonCSCDetLayerGeometryBuilder";

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