BTLDetLayerGeometryBuilder.cc

Go to the documentation of this file.

#include <RecoMTD/DetLayers/src/BTLDetLayerGeometryBuilder.h>

#include <DataFormats/ForwardDetId/interface/BTLDetId.h>
#include <Geometry/CommonDetUnit/interface/GeomDet.h>
#include <RecoMTD/DetLayers/interface/MTDTrayBarrelLayer.h>
#include <RecoMTD/DetLayers/interface/MTDDetTray.h>

#include <Utilities/General/interface/precomputed_value_sort.h>
#include <Geometry/CommonDetUnit/interface/DetSorting.h>

#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include <iostream>

using namespace std;

BTLDetLayerGeometryBuilder::BTLDetLayerGeometryBuilder() {
}

BTLDetLayerGeometryBuilder::~BTLDetLayerGeometryBuilder() {
}

vector<DetLayer*> 
BTLDetLayerGeometryBuilder::buildLayers(const MTDGeometry& geo) {
        
  const std::string metname = "MTD|RecoMTD|RecoMTDDetLayers|BTLDetLayerGeometryBuilder";

  vector<DetLayer*> detlayers;
  vector<MTDTrayBarrelLayer*> result;
  
  vector<const DetRod*> btlDetTrays;
      
  for(unsigned tray = BTLDetId::MIN_ROD; tray <= BTLDetId::HALF_ROD; ++tray) {
    
    vector<const GeomDet*> geomDets;
    for(unsigned module = 1; module <= BTLDetId::kModulesPerROD; ++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(metname) 
        << "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(metname) << "  New BTLDetTray with " << geomDets.size()
            << " modules at R=" << btlDetTrays.back()->position().perp()
            << ", phi=" << btlDetTrays.back()->position().phi() << std::endl;
    }    
  }

  precomputed_value_sort(btlDetTrays.begin(), btlDetTrays.end(), geomsort::ExtractPhi<GeometricSearchDet,float>());
  result.push_back(new MTDTrayBarrelLayer(btlDetTrays));  
  LogDebug(metname)
   << "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;
}