ETLDetLayerGeometryBuilder Class Reference

#include <ETLDetLayerGeometryBuilder.h>

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

Private Member Functions
	ETLDetLayerGeometryBuilder ()

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

Detailed Description

Build the ETL DetLayers.

Author

L. Gray - FNAL

Definition at line 18 of file ETLDetLayerGeometryBuilder.h.

Constructor & Destructor Documentation

ETLDetLayerGeometryBuilder::ETLDetLayerGeometryBuilder ( )

inlineprivate

Definition at line 26 of file ETLDetLayerGeometryBuilder.h.

References buildLayer(), makeMuonMisalignmentScenario::endcap, isFront(), makeDetRing(), and relativeConstraints::ring.

{}

Member Function Documentation

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

staticprivate

Definition at line 37 of file ETLDetLayerGeometryBuilder.cc.

References TauDecayModes::dec, MTDGeometry::idToDet(), ETLDetId::kETLmoduleMask, LogTrace, metname, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeomDet::position(), DetId::rawId(), mps_fire::result, relativeConstraints::ring, AlCaHLTBitMon_QueryRunRegistry::string, and PV3DBase< T, PVType, FrameType >::z().

Referenced by ETLDetLayerGeometryBuilder().

                                                                                          {
  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|ETLDetLayerGeometryBuilder";
  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(metname) 
      << "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 << std::endl;
      }
    }

    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 MTDRingForwardDoubleLayer(frontRings, backRings);  
  LogTrace(metname) 
    << "New MTDRingForwardLayer with " << frontRings.size() 
    << " and " << backRings.size()
    << " rings, at Z " << result->position().z()
    << " R1: " << result->specificSurface().innerRadius()
    << " R2: " << result->specificSurface().outerRadius() << std::endl; 
  return result;
}

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

static

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

Definition at line 18 of file ETLDetLayerGeometryBuilder.cc.

References makeMuonMisalignmentScenario::endcap, mps_fire::result, and relativeConstraints::ring.

Referenced by MTDDetLayerGeometryESProducer::produce().

                                                              {

  vector<DetLayer*> result[2]; // one for each endcap
  
  for(unsigned endcap=0; endcap<2; ++endcap) {    
    // there is only one layer for ETL right now, maybe more later
    for(unsigned layer = 0; layer <= 0; ++layer) {
      vector<unsigned> rings;      
      for(unsigned ring = 1; ring <= 12; ++ring) {
        rings.push_back(ring);
      }
      MTDRingForwardDoubleLayer* thelayer = buildLayer(endcap, layer, rings, geo);          
      if (thelayer) result[endcap].push_back(thelayer);
    }
  }
  pair<vector<DetLayer*>, vector<DetLayer*> > res_pair(result[0], result[1]); 
  return res_pair;
}

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

staticprivate

Definition at line 101 of file ETLDetLayerGeometryBuilder.cc.

Referenced by ETLDetLayerGeometryBuilder().

{  
  return (module+1)%2;
}

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

staticprivate

Definition at line 108 of file ETLDetLayerGeometryBuilder.cc.

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

Referenced by ETLDetLayerGeometryBuilder().

{
    const std::string metname = "MTD|RecoMTD|RecoMTDDetLayers|ETLDetLayerGeometryBuilder";


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