DetGeomDescBuilder.cc

Go to the documentation of this file.

#include "Geometry/VeryForwardGeometryBuilder/interface/DetGeomDescBuilder.h"

#include "DataFormats/CTPPSDetId/interface/CTPPSDetId.h"
#include "DetectorDescription/DDCMS/interface/DDDetector.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

/*
 * Generic function to build geo (tree of DetGeomDesc) from old DD compact view.
 */
std::unique_ptr<DetGeomDesc> detgeomdescbuilder::buildDetGeomDescFromCompactView(const DDCompactView& myCompactView,
                                                                                 const bool isRun2) {
  // Create DDFilteredView (no filter!!)
  DDPassAllFilter filter;
  DDFilteredView fv(myCompactView, filter);

  // Geo info: root node.
  auto geoInfoRoot = std::make_unique<DetGeomDesc>(fv, isRun2);

  // Construct the tree of children geo info (DetGeomDesc).
  detgeomdescbuilder::buildDetGeomDescDescendants(fv, geoInfoRoot.get(), isRun2);

  edm::LogInfo("PPSGeometryESProducer") << "Successfully built geometry.";

  return geoInfoRoot;
}

/*
 * Depth-first search recursion.
 * Construct the tree of children geo info (DetGeomDesc) (old DD).
 */
void detgeomdescbuilder::buildDetGeomDescDescendants(DDFilteredView& fv, DetGeomDesc* geoInfo, const bool isRun2) {
  // Leaf
  if (!fv.firstChild())
    return;

  do {
    // Create node
    DetGeomDesc* child = new DetGeomDesc(fv, isRun2);

    // legacy Run2 z sign fix for diamond detectors
    const auto& detId = child->geographicalID();
    if (isRun2 && detId.subdetId() == CTPPSDetId::sdTimingDiamond)
      child->invertZSign();

    // add the to the geoInfoParent's list.
    geoInfo->addComponent(child);

    // Recursion
    buildDetGeomDescDescendants(fv, child, isRun2);
  } while (fv.nextSibling());

  fv.parent();
}

/*
 * Generic function to build geo (tree of DetGeomDesc) from DD4hep compact view.
 */
std::unique_ptr<DetGeomDesc> detgeomdescbuilder::buildDetGeomDescFromCompactView(
    const cms::DDCompactView& myCompactView, const bool isRun2) {
  // create DDFilteredView (no filter!!)
  const cms::DDDetector* mySystem = myCompactView.detector();
  const dd4hep::Volume& worldVolume = mySystem->worldVolume();
  cms::DDFilteredView fv(mySystem, worldVolume);
  if (fv.next(0) == false) {
    edm::LogError("PPSGeometryESProducer") << "Filtered view is empty. Cannot build.";
  }

  // Geo info: root node.
  auto geoInfoRoot = std::make_unique<DetGeomDesc>(fv, isRun2);

  // Construct the tree of children geo info (DetGeomDesc).
  do {
    // Create node
    DetGeomDesc* child = new DetGeomDesc(fv, isRun2);

    // legacy Run2 z sign fix for diamond detectors
    const auto& detId = child->geographicalID();
    if (isRun2 && detId.subdetId() == CTPPSDetId::sdTimingDiamond)
      child->invertZSign();

    // add the node to the geoInfoRoot's list.
    geoInfoRoot->addComponent(child);
  } while (fv.next(0));

  edm::LogInfo("PPSGeometryESProducer") << "Successfully built geometry, it has " << (geoInfoRoot->components()).size()
                                        << " DetGeomDesc nodes.";

  return geoInfoRoot;
}