DDTECOptoHybAlgo.cc File Reference

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

Go to the source code of this file.

Functions
static long	algorithm (Detector &, cms::DDParsingContext &ctxt, xml_h e)
	DD4HEP_OPEN_PLUGIN (dd4hep, ddcms_det_element_DDCMS_track_DDTECOptoHybAlgo)

Function Documentation

algorithm()

static long algorithm ( Detector &  , cms::DDParsingContext &  ctxt, xml_h  e  )

static

Definition at line 11 of file DDTECOptoHybAlgo.cc.

References angle(), particleFlowDisplacedVertex_cfi::angles, writedatasetfile::args, angle_units::operators::convertRadToDeg(), funct::cos(), MillePedeFileConverter_cfg::e, cms::makeRotation3D(), cms::DDNamespace::name(), cms::DDNamespace::nsName(), PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), idealTransformation::rotation, cms::DDNamespace::rotation(), cms::DDParsingContext::rotations, funct::sin(), theta(), to_string(), and cms::DDNamespace::volume().

                                                                       {
  cms::DDNamespace ns(ctxt, e, true);
  DDAlgoArguments args(ctxt, e);
  int startCopyNo = args.value<int>("StartCopyNo");
  double rpos = args.value<double>("Rpos");
  double zpos = args.value<double>("Zpos");
  double optoHeight = args.value<double>("OptoHeight");
  double optoWidth = args.value<double>("OptoWidth");
  vector<double> angles = args.value<vector<double> >("Angles");
  Volume child = ns.volume(args.value<string>("ChildName"));
  Volume mother = ns.volume(args.parentName());

  edm::LogVerbatim("TECGeom") << "Parent " << mother.name() << " Child " << child.name() << " NameSpace " << ns.name();
  edm::LogVerbatim("TECGeom") << "Height of the Hybrid " << optoHeight << " and Width " << optoWidth << "Rpos " << rpos
                              << " Zpos " << zpos << " StartCopyNo " << startCopyNo << " Number " << angles.size();

  // given r positions are for the lower left corner
  rpos += optoHeight / 2;
  int copyNo = startCopyNo;
  for (double angle : angles) {
    double phix = -angle;
    // given phi positions are for the lower left corner
    phix += asin(optoWidth / 2 / rpos);
    double xpos = rpos * cos(phix);
    double ypos = rpos * sin(phix);
    Position tran(xpos, ypos, zpos);

    Rotation3D rotation;
    double phiy = phix + 90._deg;
    double phideg = convertRadToDeg(phix);
    if (phideg != 0) {
      string rotstr = ns.nsName(child.name()) + std::to_string(phideg * 1000.);
      auto irot = ctxt.rotations.find(ns.prepend(rotstr));
      if (irot != ctxt.rotations.end()) {
        rotation = ns.rotation(ns.prepend(rotstr));
      } else {
        double theta = 90._deg;
        edm::LogVerbatim("TECGeom") << "test: Creating a new "
                                    << "rotation: " << rotstr << "\t90., " << convertRadToDeg(phix) << ", 90.,"
                                    << convertRadToDeg(phiy) << ", 0, 0";
        rotation = makeRotation3D(theta, phix, theta, phiy, 0., 0.);
      }
    }
    mother.placeVolume(child, copyNo, Transform3D(rotation, tran));
    edm::LogVerbatim("TECGeom") << "test " << child.name() << " number " << copyNo << " positioned in " << mother.name()
                                << " at " << tran << " with " << rotation;
    copyNo++;
  }
  edm::LogVerbatim("TECGeom") << "<<== End of DDTECOptoHybAlgo construction ...";
  return 1;
}

DD4HEP_OPEN_PLUGIN()

DD4HEP_OPEN_PLUGIN	(	dd4hep	,
		ddcms_det_element_DDCMS_track_DDTECOptoHybAlgo
	)

Definition at line 64 of file DDTECOptoHybAlgo.cc.