MuonGEMDetLayerGeometryBuilder Class Reference

#include <MuonGEMDetLayerGeometryBuilder.h>

Public Member Functions
	MuonGEMDetLayerGeometryBuilder ()
	Constructor (disabled, only static access is allowed) More...
virtual	~MuonGEMDetLayerGeometryBuilder ()
	Destructor. More...

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

Static Private Member Functions
static MuRingForwardDoubleLayer *	buildLayer (int endcap, std::vector< int > &rings, int station, int layer, std::vector< int > &chambers, std::vector< int > &rolls, const GEMGeometry &geo)
static bool	isFront (const GEMDetId &gemId)
static MuDetRing *	makeDetRing (std::vector< const GeomDet * > &geomDets)

Detailed Description

Build the GEM DetLayers.

Author

R. Radogna

Definition at line 20 of file MuonGEMDetLayerGeometryBuilder.h.

Constructor & Destructor Documentation

MuonGEMDetLayerGeometryBuilder::MuonGEMDetLayerGeometryBuilder ( )

inline

Constructor (disabled, only static access is allowed)

Definition at line 23 of file MuonGEMDetLayerGeometryBuilder.h.

{}

virtual MuonGEMDetLayerGeometryBuilder::~MuonGEMDetLayerGeometryBuilder ( )

virtual

Destructor.

Member Function Documentation

pair< vector< DetLayer * >, vector< DetLayer * > > MuonGEMDetLayerGeometryBuilder::buildEndcapLayers ( const GEMGeometry &  geo )

static

Builds the forward (+Z, return.first) and backward (-Z, return.second) layers. Both vectors are sorted inside-out

Definition at line 26 of file MuonGEMDetLayerGeometryBuilder.cc.

References chambers, Reference_intrackfit_cff::endcap, GEMDetId::maxChamberId, GEMDetId::maxLayerId, GEMDetId::maxRingId, GEMDetId::maxRollId, GEMDetId::maxStationId, GEMDetId::minChamberId, GEMDetId::minLayerId, GEMDetId::minRingId, GEMDetId::minRollId, GEMDetId::minStationId, mps_fire::result, relativeConstraints::ring, and relativeConstraints::station.

Referenced by MuonDetLayerGeometryESProducer::produce().

                                                                        {
  
  vector<DetLayer*> result[2];

  for (int endcap = -1; endcap<=1; endcap+=2) {
    int iendcap = (endcap==1) ? 0 : 1; // +1: forward, -1: backward

    for(int station = GEMDetId::minStationId; station < GEMDetId::maxStationId; ++station) {
      for(int layer = GEMDetId::minLayerId; layer <= GEMDetId::maxLayerId; ++layer) { 
    vector<int> rolls;      
    std::vector<int> rings;
    std::vector<int> chambers;
    for(int ring = GEMDetId::minRingId; ring <= GEMDetId::maxRingId; ++ring) {
      rings.push_back(ring);
    }
    for(int roll = GEMDetId::minRollId+1; roll <= GEMDetId::maxRollId; ++roll) {
      rolls.push_back(roll);
    }
    for(int chamber = GEMDetId::minChamberId; chamber <= GEMDetId::maxChamberId; chamber++ ){
      chambers.push_back(chamber);
    }

    MuRingForwardDoubleLayer* ringLayer = buildLayer(endcap, rings, station, layer, chambers, rolls, geo);          

    if (ringLayer) result[iendcap].push_back(ringLayer);

      }

    }

    
  }
  pair<vector<DetLayer*>, vector<DetLayer*> > res_pair(result[0], result[1]); 

  return res_pair;

}

MuRingForwardDoubleLayer * MuonGEMDetLayerGeometryBuilder::buildLayer ( int  endcap, std::vector< int > &  rings, int  station, int  layer, std::vector< int > &  chambers, std::vector< int > &  rolls, const GEMGeometry &  geo  )

staticprivate

Definition at line 67 of file MuonGEMDetLayerGeometryBuilder.cc.

References GEMGeometry::idToDet(), LogTrace, metname, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GeomDet::position(), precomputed_value_sort(), mps_fire::result, relativeConstraints::ring, AlCaHLTBitMon_QueryRunRegistry::string, and PV3DBase< T, PVType, FrameType >::z().

                                               {

  const std::string metname = "Muon|RecoMuon|RecoMuonDetLayers|MuonGEMDetLayerGeometryBuilder";
  MuRingForwardDoubleLayer * result = 0;
  vector<const ForwardDetRing*> frontRings, backRings;


  for (std::vector<int>::iterator ring=rings.begin(); ring!=rings.end()-2;ring++){ 

    for (vector<int>::iterator roll = rolls.begin(); roll!=rolls.end(); roll++) {    

      vector<const GeomDet*> frontDets, backDets;

      for(std::vector<int>::iterator chamber=chambers.begin()+1; chamber<chambers.end(); chamber++) {
          GEMDetId gemId(endcap,(*ring), station,layer,(*chamber), (*roll));

      const GeomDet* geomDet = geo.idToDet(gemId);
      
      if (geomDet !=0) {
        bool isInFront = isFront(gemId);
        if(isInFront)
            {
              frontDets.push_back(geomDet);
            }
            else 
            {
              backDets.push_back(geomDet);
            }
        LogTrace(metname) << "get GEM Endcap roll "
                  << gemId
                              << (isInFront ? "front" : "back ")
                  << " at R=" << geomDet->position().perp()
                  << ", phi=" << geomDet->position().phi()
                              << ", Z=" << geomDet->position().z();
      }
      }

      if (frontDets.size()!=0) {
    precomputed_value_sort(frontDets.begin(), frontDets.end(), geomsort::DetPhi());
    frontRings.push_back(new MuDetRing(frontDets));
    LogTrace(metname) << "New front ring with " << frontDets.size()
              << " chambers at z="<< frontRings.back()->position().z();
      }
      if (backDets.size()!=0) {
        precomputed_value_sort(backDets.begin(), backDets.end(), geomsort::DetPhi());
        backRings.push_back(new MuDetRing(backDets));
        LogTrace(metname) << "New back ring with " << backDets.size()
                          << " chambers at z="<< backRings.back()->position().z();
      }

    }

  }

  // How should they be sorted?
  //    precomputed_value_sort(muDetRods.begin(), muDetRods.end(), geomsort::ExtractZ<GeometricSearchDet,float>());                                   
  if(backRings.size()!=0 && frontRings.size()!=0) result = new MuRingForwardDoubleLayer(frontRings, backRings);
    else result = 0;
  if(result != 0){
    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;

}

bool MuonGEMDetLayerGeometryBuilder::isFront ( const GEMDetId &  gemId )

staticprivate

Definition at line 141 of file MuonGEMDetLayerGeometryBuilder.cc.

References GEMDetId::chamber(), and mps_fire::result.

{

  bool result = false;
  int chamber = gemId.chamber();

    if(chamber%2 == 0) result = !result;

    return result;
}

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

staticprivate

Definition at line 152 of file MuonGEMDetLayerGeometryBuilder.cc.

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

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


    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;
}