Geometry/TrackerNumberingBuilder/plugins/CmsTrackerRingBuilder.cc - master - githubmirror/cmssw

Geometry/TrackerNumberingBuilder/plugins/CmsTrackerRingBuilder.cc

90 строк3 KB

Gabrielle Hugo

Elegant fix of aweful custom-made go down / look at child / go back up in CmsTrackerRingBuilder. This does the same as before for old DD, but allows to be compatible with DDDhep navigation!!

22 сен 2020, 19:01

22 сен 2020, 19:01368d5aa

#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerRingBuilder.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "Geometry/TrackerNumberingBuilder/plugins/ExtractStringFromDDD.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "Geometry/TrackerNumberingBuilder/plugins/CmsDetConstruction.h"
#include "Geometry/TrackerNumberingBuilder/interface/trackerStablePhiSort.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <vector>
#include <bitset>

template <class FilteredView>
void CmsTrackerRingBuilder<FilteredView>::buildComponent(FilteredView& fv, GeometricDet* g, const std::string& s) {
  moduleName_ = ExtractStringFromDDD<FilteredView>::getString("TkDDDStructure", &fv);

  CmsDetConstruction<FilteredView> theCmsDetConstruction;
  theCmsDetConstruction.buildComponent(fv, g, s);
}

template <class FilteredView>
void CmsTrackerRingBuilder<FilteredView>::sortNS(FilteredView& fv, GeometricDet* det) {
  GeometricDet::ConstGeometricDetContainer& comp = det->components();
  GeometricDet::GeometricDetContainer compfw;
  GeometricDet::GeometricDetContainer compbw;

  switch (comp.front()->type()) {
    case GeometricDet::mergedDet:
      trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiGluedModule);
      break;
    case GeometricDet::DetUnit:
      trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhi);
      break;
    default:
      edm::LogError("CmsTrackerRingBuilder")
          << "ERROR - wrong SubDet to sort..... " << det->components().front()->type();
  }

  static std::string const TECGluedDet("TECGluedDet");
  static std::string const TECDet("TECDet");

  // TEC
  // Module Number: 3 bits [1,...,5 at most]
  if (moduleName_ == TECGluedDet || moduleName_ == TECDet) {
    // TEC-
    if (det->translation().z() < 0 && moduleName_ == TECDet) {
      trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiMirror);
    }

    if (det->translation().z() < 0 && moduleName_ == TECGluedDet) {
      trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiGluedModuleMirror);
    }

    for (uint32_t i = 0; i < comp.size(); i++)
      det->component(i)->setGeographicalID(i + 1);

  } else {
    // TID
    // Ring Side: 2 bits [back:1 front:2]
    // Module Number: 5 bits [1,...,20 at most]
    //
    for (uint32_t i = 0; i < comp.size(); i++) {
      if (std::abs(comp[i]->translation().z()) < std::abs(det->translation().z())) {
        compfw.emplace_back(det->component(i));
      } else {
        compbw.emplace_back(det->component(i));
      }
    }

    for (uint32_t i = 0; i < compbw.size(); i++) {
      uint32_t temp = i + 1;
      temp |= (1 << 5);
      compbw[i]->setGeographicalID(temp);
    }

    for (uint32_t i = 0; i < compfw.size(); i++) {
      uint32_t temp = i + 1;
      temp |= (2 << 5);
      compfw[i]->setGeographicalID(temp);
    }

    det->clearComponents();
    det->addComponents(compfw);
    det->addComponents(compbw);
  }
}

template class CmsTrackerRingBuilder<DDFilteredView>;
template class CmsTrackerRingBuilder<cms::DDFilteredView>;