PrimitivesFactory.cpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PrimitivesFactory.cpp

// This file is part of the Acts project.
//
// Copyright (C) 2023 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#include "Acts/Plugins/DD4hep/DD4hepConversionHelpers.hpp"

#include <DD4hep/DetFactoryHelper.h>
#include <DD4hep/Objects.h>
#include <XML/Utilities.h>

#include "DD4hepTestsHelper.hpp"

using namespace std;
using namespace dd4hep;

static Ref_t create_cylinder(Detector &dd, xml_h xml,
                             SensitiveDetector /*sens*/) {
  xml_det_t x_det = xml;
  string detName = x_det.nameStr();

  // Make Volume
  xml_comp_t x_det_tubs = x_det.child(_U(tubs));

  // Make DetElement
  DetElement cylinderElement(detName, x_det.id());
  dd4hep::xml::setDetectorTypeFlag(xml, cylinderElement);

  string shapeName = x_det_tubs.nameStr();
  double phiMin = Acts::getAttrValueOr<double>(x_det_tubs, "phimin", 0.);
  double phiMax = Acts::getAttrValueOr<double>(x_det_tubs, "phimax", 360.);
  Tube tubeShape(shapeName, x_det_tubs.rmin(), x_det_tubs.rmax(),
                 x_det_tubs.dz(), phiMin, phiMax);
  Volume tubeVolume(detName, tubeShape, dd.material(x_det_tubs.materialStr()));
  tubeVolume.setVisAttributes(dd, x_det.visStr());

  // Place it in the mother
  Volume motherVolume = dd.pickMotherVolume(cylinderElement);
  PlacedVolume placedTube = motherVolume.placeVolume(
      tubeVolume, DD4hepTestsHelper::createTransform(x_det_tubs));
  placedTube.addPhysVolID(detName, cylinderElement.id());
  cylinderElement.setPlacement(placedTube);

  // And return the element for further parsing
  return cylinderElement;
}

DECLARE_DETELEMENT(Cylinder, create_cylinder)








static Ref_t create_disc(Detector &dd, xml_h xml, SensitiveDetector /*sens*/) {
  xml_det_t x_det = xml;
  string detName = x_det.nameStr();

  // Make Volume
  xml_comp_t x_det_tubs = x_det.child(_U(tubs));

  // Make DetElement
  DetElement discElement(detName, x_det.id());
  dd4hep::xml::setDetectorTypeFlag(xml, discElement);

  string shapeName = x_det_tubs.nameStr();
  double phiMin = Acts::getAttrValueOr<double>(x_det_tubs, "phimin", 0.);
  double phiMax = Acts::getAttrValueOr<double>(x_det_tubs, "phimax", 360.);

  Tube discShape = Tube(shapeName, x_det_tubs.rmin(), x_det_tubs.rmax(),
                        x_det_tubs.dz(), phiMin, phiMax);

  Volume discVolume(detName, discShape, dd.material(x_det_tubs.materialStr()));
  discVolume.setVisAttributes(dd, x_det.visStr());

  // Place it in the mother
  Volume motherVolume = dd.pickMotherVolume(discElement);
  PlacedVolume placedTube = motherVolume.placeVolume(
      discVolume, DD4hepTestsHelper::createTransform(x_det_tubs));
  placedTube.addPhysVolID(detName, discElement.id());
  discElement.setPlacement(placedTube);

  // And return the element for further parsing
  return discElement;
}

DECLARE_DETELEMENT(Disc, create_disc)








static Ref_t create_rectangle(Detector &dd, xml_h xml,
                              SensitiveDetector /*sens*/) {
  xml_det_t x_det = xml;
  string detName = x_det.nameStr();

  // Make Volume
  xml_comp_t x_det_box = x_det.child(_U(box));

  // Make DetElement
  DetElement rectElement(detName, x_det.id());
  dd4hep::xml::setDetectorTypeFlag(xml, rectElement);

  string shapeName = x_det_box.nameStr();
  Box rectShape(shapeName, 0.5 * x_det_box.dx(), 0.5 * x_det_box.dy(),
                0.5 * x_det_box.dz());

  Volume rectVolume(detName, rectShape, dd.material(x_det_box.materialStr()));
  rectVolume.setVisAttributes(dd, x_det.visStr());

  // Place it in the mother
  Volume motherVolume = dd.pickMotherVolume(rectElement);
  PlacedVolume placedRect = motherVolume.placeVolume(
      rectVolume, DD4hepTestsHelper::createTransform(x_det_box));
  placedRect.addPhysVolID(detName, rectElement.id());
  rectElement.setPlacement(placedRect);

  // And return the element for further parsing
  return rectElement;
}

DECLARE_DETELEMENT(Rectangle, create_rectangle)








static Ref_t create_trapezoid(Detector &dd, xml_h xml,
                              SensitiveDetector /*sens*/) {
  xml_det_t x_det = xml;
  string detName = x_det.nameStr();

  // Make Volume
  xml_comp_t x_det_trap = x_det.child(_U(trap));

  // Make DetElement
  DetElement trapElement(detName, x_det.id());
  dd4hep::xml::setDetectorTypeFlag(xml, trapElement);

  string shapeName = x_det_trap.nameStr();

  // Due to convention this causes an axis flip on x
  Trapezoid trapShape(x_det_trap.x1(), x_det_trap.x2(), 0.5 * x_det_trap.dz(),
                      0.5 * x_det_trap.dz(), 0.5 * x_det_trap.dy());

  Volume trapVolume(detName, trapShape, dd.material(x_det_trap.materialStr()));
  trapVolume.setVisAttributes(dd, x_det.visStr());

  // Place it in the mother
  Volume motherVolume = dd.pickMotherVolume(trapElement);
  PlacedVolume placedTrap = motherVolume.placeVolume(
      trapVolume, DD4hepTestsHelper::createTransform(x_det_trap));
  placedTrap.addPhysVolID(detName, trapElement.id());
  trapElement.setPlacement(placedTrap);

  // And return the element for further parsing
  return trapElement;
}

DECLARE_DETELEMENT(Trapezoid, create_trapezoid)