DDTrackerPhiAlgo.cc

Go to the documentation of this file.

#include "DD4hep/DetFactoryHelper.h"
#include "DataFormats/Math/interface/CMSUnits.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;
using namespace dd4hep;
using namespace cms;
using namespace cms_units::operators;

static long algorithm(Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
  cms::DDNamespace ns(ctxt, e, true);
  DDAlgoArguments args(ctxt, e);
  Volume mother = ns.volume(args.parentName());
  Volume child = ns.volume(args.childName());
  int startcn = args.find("StartCopyNo") ? args.value<int>("StartCopyNo") : 1;
  int incrcn = args.find("IncrCopyNo") ? args.value<int>("IncrCopyNo") : 1;
  vector<double> phi = args.value<vector<double> >("Phi");    // Phi values
  vector<double> zpos = args.value<vector<double> >("ZPos");  // Z positions
  int numcopies = args.find("NumCopies") ? args.value<int>("NumCopies") : int(phi.size());
  double radius = args.value<double>("Radius");
  double tilt = args.value<double>("Tilt");

  if (numcopies != int(phi.size())) {
    edm::LogVerbatim("TrackerGeom") << "error: Parameter "
                                    << "NumCopies does not agree with the size "
                                    << "of the Phi vector. It was adjusted to "
                                    << "be the size of the Phi vector and may "
                                    << "lead to crashes or errors.";
  }
  edm::LogVerbatim("TrackerGeom") << "debug: Parameters for position"
                                  << "ing:: "
                                  << " Radius " << radius << " Tilt " << tilt << " Copies " << phi.size() << " at";
  for (int i = 0; i < (int)(phi.size()); i++)
    edm::LogVerbatim("TrackerGeom") << "\t[" << i << "] phi = " << phi[i] << " z = " << zpos[i];
  edm::LogVerbatim("TrackerGeom") << "debug: Parent " << mother.name() << "\tChild " << child.name() << " NameSpace "
                                  << ns.name();

  double theta = 90._deg;
  int ci = startcn;
  for (int i = 0; i < numcopies; ++i) {
    double phix = phi[i] + tilt;
    double phiy = phix + 90._deg;
    double xpos = radius * cos(phi[i]);
    double ypos = radius * sin(phi[i]);
    Rotation3D rot = makeRotation3D(theta, phix, theta, phiy, 0., 0.);
    Position tran(xpos, ypos, zpos[i]);
    /* PlacedVolume pv = */ mother.placeVolume(child, ci, Transform3D(rot, tran));
    edm::LogVerbatim("TrackerGeom") << "test: " << child.name() << " number " << ci << " positioned in "
                                    << mother.name() << " at " << tran << " with " << rot;
    ci = ci + incrcn;
  }
  return 1;
}

// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_track_DDTrackerPhiAlgo, algorithm)