DetectorDescription/Core/plugins/DDAngular.cc - master - githubmirror/cmssw

DetectorDescription/Core/plugins/DDAngular.cc
154 строки6 KB
Cms Build
Clang-Format
29 май 2019, 08:28
29 май 2019, 08:2896591f2
#include "DetectorDescription/Core/interface/DDAlgorithm.h"
#include "DetectorDescription/Core/interface/DDAlgorithmFactory.h"
#include "DetectorDescription/Core/interface/DDCompactView.h"
#include "DetectorDescription/Core/interface/DDCurrentNamespace.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/DDSplit.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
#include "DetectorDescription/Core/interface/DDRotationMatrix.h"
#include "DetectorDescription/Core/interface/DDTranslation.h"
#include "DetectorDescription/Core/interface/DDTypes.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/PluginManager/interface/PluginFactory.h"

#include <cmath>
#include <memory>

using namespace geant_units::operators;

class DDAngular : public DDAlgorithm {
public:
  DDAngular()
      : m_n(1), m_startCopyNo(1), m_incrCopyNo(1), m_startAngle(0.), m_rangeAngle(360._deg), m_radius(0.), m_delta(0.) {}

  void initialize(const DDNumericArguments& nArgs,
                  const DDVectorArguments& vArgs,
                  const DDMapArguments& mArgs,
                  const DDStringArguments& sArgs,
                  const DDStringVectorArguments& vsArgs) override;

  void execute(DDCompactView& cpv) override;

private:
  DD3Vector fUnitVector(double theta, double phi);
  int m_n;                            //Number of copies
  int m_startCopyNo;                  //Start Copy number
  int m_incrCopyNo;                   //Increment in Copy number
  double m_startAngle;                //Start angle
  double m_rangeAngle;                //Range in angle
  double m_radius;                    //Radius
  double m_delta;                     //Increment in phi
  std::vector<double> m_center;       //Phi values
  std::vector<double> m_rotateSolid;  //Rotation of the solid values

  std::string m_idNameSpace;                        //Namespace of this and ALL sub-parts
  std::pair<std::string, std::string> m_childNmNs;  //Child name
                                                    //Namespace of the child

  DDRotationMatrix m_solidRot;  //Rotation of the solid
};

void DDAngular::initialize(const DDNumericArguments& nArgs,
                           const DDVectorArguments& vArgs,
                           const DDMapArguments&,
                           const DDStringArguments& sArgs,
                           const DDStringVectorArguments&) {
  m_n = int(nArgs["N"]);
  m_startCopyNo = int(nArgs["StartCopyNo"]);
  m_incrCopyNo = int(nArgs["IncrCopyNo"]);
  m_rangeAngle = nArgs["RangeAngle"];
  m_startAngle = nArgs["StartAngle"];
  m_radius = nArgs["Radius"];
  m_center = vArgs["Center"];
  m_rotateSolid = vArgs["RotateSolid"];

  m_solidRot = DDRotationMatrix();

  if (std::abs(m_rangeAngle - 360.0_deg) < 0.001_deg) {
    m_delta = m_rangeAngle / double(m_n);
  } else {
    if (m_n > 1) {
      m_delta = m_rangeAngle / double(m_n - 1);
    } else {
      m_delta = 0.;
    }
  }

  LogDebug("DDAlgorithm") << "DDAngular: Parameters for position"
                          << "ing:: n " << m_n << " Start, Range, Delta " << convertRadToDeg(m_startAngle) << " "
                          << convertRadToDeg(m_rangeAngle) << " " << convertRadToDeg(m_delta) << " Radius " << m_radius
                          << " Centre " << m_center[0] << ", " << m_center[1] << ", " << m_center[2];

  //======= collect data concerning the rotation of the solid
  using sz_type = std::map<std::string, std::vector<double> >::mapped_type::size_type;
  sz_type sz = m_rotateSolid.size();
  if (sz % 3) {
    LogDebug("DDAlgorithm") << "\trotateSolid must occur 3*n times (defining n subsequent rotations)\n"
                            << "\t  currently it appears " << sz << " times!\n";
  }
  for (sz_type i = 0; i < sz; i += 3) {
    if ((m_rotateSolid[i] > 180._deg) || (m_rotateSolid[i] < 0._deg)) {
      LogDebug("DDAlgorithm") << "\trotateSolid \'theta\' must be in range [0,180*deg]\n"
                              << "\t  currently it is " << convertRadToDeg(m_rotateSolid[i]) << "*deg in rotateSolid["
                              << double(i) << "]!\n";
    }
    DDAxisAngle temp(fUnitVector(m_rotateSolid[i], m_rotateSolid[i + 1]), m_rotateSolid[i + 2]);
    LogDebug("DDAlgorithm") << "  rotsolid[" << i << "] axis=" << temp.Axis()
                            << " rot.angle=" << convertRadToDeg(temp.Angle());
    m_solidRot = temp * m_solidRot;
  }

  m_idNameSpace = DDCurrentNamespace::ns();
  m_childNmNs = DDSplit(sArgs["ChildName"]);
  if (m_childNmNs.second.empty())
    m_childNmNs.second = DDCurrentNamespace::ns();

  DDName parentName = parent().name();
  LogDebug("DDAlgorithm") << "DDAngular: Parent " << parentName << "\tChild " << m_childNmNs.first << "\tNameSpace "
                          << m_childNmNs.second;
}

void DDAngular::execute(DDCompactView& cpv) {
  DDName mother = parent().name();
  DDName ddname(m_childNmNs.first, m_childNmNs.second);
  double theta = 90._deg;
  int copy = m_startCopyNo;
  double phi = m_startAngle;

  for (int i = 0; i < m_n; ++i) {
    double phix = phi;
    double phiy = phix + 90._deg;
    double phideg = convertRadToDeg(phix);

    std::string rotstr = m_childNmNs.first + "_" + std::to_string(phideg * 10.);
    DDRotation rotation = DDRotation(DDName(rotstr, m_idNameSpace));
    if (!rotation) {
      LogDebug("DDAlgorithm") << "DDAngular: Creating a new "
                              << "rotation: " << rotstr << "\t90., " << convertRadToDeg(phix) << ", 90.,"
                              << convertRadToDeg(phiy) << ", 0, 0";

      rotation = DDrot(
          DDName(rotstr, m_idNameSpace),
          std::make_unique<DDRotationMatrix>(*DDcreateRotationMatrix(theta, phix, theta, phiy, 0., 0.) * m_solidRot));
    }

    double xpos = m_radius * cos(phi) + m_center[0];
    double ypos = m_radius * sin(phi) + m_center[1];
    double zpos = m_center[2];
    DDTranslation tran(xpos, ypos, zpos);

    cpv.position(ddname, mother, copy, tran, rotation);
    LogDebug("DDAlgorithm") << "DDAngular: child " << m_childNmNs.second << ":" << m_childNmNs.first << " number "
                            << copy << " positioned in " << mother << " at " << tran << " with " << rotation << "\n";
    copy += m_incrCopyNo;
    phi += m_delta;
  }
}

DD3Vector DDAngular::fUnitVector(double theta, double phi) {
  return DD3Vector(cos(phi) * sin(theta), sin(phi) * sin(theta), cos(theta));
}

DEFINE_EDM_PLUGIN(DDAlgorithmFactory, DDAngular, "global:DDAngular");