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

Geometry/TrackerNumberingBuilder/plugins/CmsTrackerOTLayerBuilder.cc
124 строки5 KB
Ianna Osborne
Address Vincenzo's comments
05 сен 2019, 00:07
05 сен 2019, 00:07c6fbce2
#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerOTLayerBuilder.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 "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerOTRingBuilder.h"
#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerLadderBuilder.h"
#include "Geometry/TrackerNumberingBuilder/interface/trackerStablePhiSort.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <vector>
#include <bitset>

template <class FilteredView>
void CmsTrackerOTLayerBuilder<FilteredView>::buildComponent(FilteredView& fv, GeometricDet* g, const std::string& s) {
  LogTrace("DetConstruction") << " CmsTrackerOTLayerBuilder::buildComponent ";
  CmsTrackerLadderBuilder<FilteredView> theCmsTrackerLadderBuilder;
  CmsTrackerOTRingBuilder<FilteredView> theCmsTrackerOTRingBuilder;

  GeometricDet* subdet = new GeometricDet(&fv,
                                          CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
                                              ExtractStringFromDDD<FilteredView>::getString(s, &fv)));
  switch (CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
      ExtractStringFromDDD<FilteredView>::getString(s, &fv))) {
    case GeometricDet::ladder:
      theCmsTrackerLadderBuilder.build(fv, subdet, s);
      break;
    case GeometricDet::panel:
      theCmsTrackerOTRingBuilder.build(fv, subdet, s);
      break;
    default:
      edm::LogError("CmsTrackerOTLayerBuilder") << " ERROR - I was expecting a ladder or a panel, I got a "
                                                << ExtractStringFromDDD<FilteredView>::getString(s, &fv);
  }
  g->addComponent(subdet);
}

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

  //order ladder and rings together
  GeometricDet::GeometricDetContainer rods;
  GeometricDet::GeometricDetContainer ringsNeg;
  GeometricDet::GeometricDetContainer ringsPos;
  rods.clear();
  ringsNeg.clear();
  ringsPos.clear();

  for (uint32_t i = 0; i < comp.size(); i++) {
    auto component = det->component(i);
    if (component->type() == GeometricDet::ladder) {
      rods.emplace_back(component);
    } else if (component->type() == GeometricDet::panel) {
      if (component->translation().z() < 0.) {
        ringsNeg.emplace_back(component);
      } else if (component->translation().z() > 0.) {
        ringsPos.emplace_back(component);
      }
    } else {
      edm::LogError("CmsTrackerOTLayerBuilder")
          << "ERROR - wrong SubDet to sort..... " << det->components().front()->type();
    }
  }

  // negative rings
  if (!ringsNeg.empty()) {
    std::sort(ringsNeg.begin(), ringsNeg.end(), CmsTrackerLevelBuilderHelper::isLessZ);
    uint32_t totalringsNeg = ringsNeg.size();

    LogTrace("DetConstruction") << " Neg rings ordered by z: ";
    for (uint32_t ring = 0; ring < totalringsNeg; ring++) {
      uint32_t temp = ring + 1;
      temp |= (1 << 8);
      ringsNeg[ring]->setGeographicalID(temp);
      LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp
                                       << "(r: " << sqrt(ringsNeg[ring]->translation().Perp2())
                                       << ", phi: " << ringsNeg[ring]->phi()
                                       << ", z: " << ringsNeg[ring]->translation().z() << ")";
    }
  }

  // rods
  if (!rods.empty()) {
    trackerStablePhiSort(rods.begin(), rods.end(), CmsTrackerLevelBuilderHelper::getPhi);
    uint32_t totalrods = rods.size();

    LogTrace("DetConstruction") << " Rods ordered by phi: ";
    for (uint32_t rod = 0; rod < totalrods; rod++) {
      uint32_t temp = rod + 1;
      temp |= (3 << 8);
      rods[rod]->setGeographicalID(DetId(temp));
      LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp << "(r: " << sqrt(rods[rod]->translation().Perp2())
                                       << ", phi: " << rods[rod]->phi() << ", z: " << rods[rod]->translation().z()
                                       << ")";
    }
  }

  // positive rings
  if (!ringsPos.empty()) {
    std::sort(ringsPos.begin(), ringsPos.end(), CmsTrackerLevelBuilderHelper::isLessZ);
    uint32_t totalringsPos = ringsPos.size();

    LogTrace("DetConstruction") << " Pos rings ordered by z: ";
    for (uint32_t ring = 0; ring < totalringsPos; ring++) {
      uint32_t temp = ring + 1;
      temp |= (2 << 8);
      ringsPos[ring]->setGeographicalID(temp);
      LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp
                                       << "(r: " << sqrt(ringsPos[ring]->translation().Perp2())
                                       << ", phi: " << ringsPos[ring]->phi()
                                       << ", z: " << ringsPos[ring]->translation().z() << ")";
    }
  }

  det->clearComponents();
  det->addComponents(ringsNeg);
  det->addComponents(rods);
  det->addComponents(ringsPos);
}

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