DDAngular Class Reference

Inheritance diagram for DDAngular:

Public Member Functions
	DDAngular ()
void	execute (DDCompactView &cpv) override
void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override

Private Member Functions
DD3Vector	fUnitVector (double theta, double phi)

Private Attributes
std::vector< double >	m_center
std::pair< std::string, std::string >	m_childNmNs
double	m_delta
std::string	m_idNameSpace
int	m_incrCopyNo
int	m_n
double	m_radius
double	m_rangeAngle
std::vector< double >	m_rotateSolid
DDRotationMatrix	m_solidRot
double	m_startAngle
int	m_startCopyNo

Detailed Description

Definition at line 21 of file DDAngular.cc.

Constructor & Destructor Documentation

DDAngular::DDAngular ( )

inline

Definition at line 23 of file DDAngular.cc.

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

Member Function Documentation

void DDAngular::execute ( DDCompactView &  cpv )

override

Definition at line 113 of file DDAngular.cc.

References angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, funct::cos(), DDcreateRotationMatrix(), DDrot(), mps_fire::i, LogDebug, SpecificationBuilder_cfi::parent(), DDCompactView::position(), idealTransformation::rotation, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, and theta().

                                          {
  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  )

private

Definition at line 150 of file DDAngular.cc.

References funct::cos(), and funct::sin().

                                                         {
  return DD3Vector(cos(phi) * sin(theta), sin(phi) * sin(theta), cos(theta));
}

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

override

Definition at line 53 of file DDAngular.cc.

References funct::abs(), angle_units::operators::convertRadToDeg(), DDSplit(), mps_fire::i, LogDebug, DDCurrentNamespace::ns(), SpecificationBuilder_cfi::parent(), and groupFilesInBlocks::temp.

                                                           {
  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;
}

Member Data Documentation

std::vector<double> DDAngular::m_center

private

Definition at line 43 of file DDAngular.cc.

std::pair<std::string, std::string> DDAngular::m_childNmNs

private

Definition at line 47 of file DDAngular.cc.

double DDAngular::m_delta

private

Definition at line 42 of file DDAngular.cc.

std::string DDAngular::m_idNameSpace

private

Definition at line 46 of file DDAngular.cc.

int DDAngular::m_incrCopyNo

private

Definition at line 38 of file DDAngular.cc.

int DDAngular::m_n

private

Definition at line 36 of file DDAngular.cc.

double DDAngular::m_radius

private

Definition at line 41 of file DDAngular.cc.

double DDAngular::m_rangeAngle

private

Definition at line 40 of file DDAngular.cc.

std::vector<double> DDAngular::m_rotateSolid

private

Definition at line 44 of file DDAngular.cc.

DDRotationMatrix DDAngular::m_solidRot

private

Definition at line 50 of file DDAngular.cc.

double DDAngular::m_startAngle

private

Definition at line 39 of file DDAngular.cc.

int DDAngular::m_startCopyNo

private

Definition at line 37 of file DDAngular.cc.