DDRotation.cc File Reference

#include <cstdio>
#include <atomic>
#include <cmath>
#include <sstream>
#include <string>
#include "DetectorDescription/Core/interface/DDRotationMatrix.h"
#include "DetectorDescription/Core/interface/DDTranslation.h"
#include "DetectorDescription/Core/interface/DDBase.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "Math/GenVector/AxisAngle.h"
#include "Math/GenVector/Cartesian3D.h"
#include "Math/GenVector/DisplacementVector3D.h"

Go to the source code of this file.

Functions
DDRotation	DDanonymousRot (std::unique_ptr< DDRotationMatrix > rot)
	Defines a anonymous rotation or rotation-reflection matrix. More...
std::unique_ptr< DDRotationMatrix >	DDcreateRotationMatrix (double thetaX, double phiX, double thetaY, double phiY, double thetaZ, double phiZ)
	create a new DDRotationMatrix in the GEANT3 style. More...
DDRotation	DDrot (const DDName &ddname, std::unique_ptr< DDRotationMatrix > rot)
	Definition of a uniquely identifiable rotation matrix named by DDName name. More...
DDRotation	DDrot (const DDName &ddname, double thetaX, double phiX, double thetaY, double phiY, double thetaZ, double phiZ)
	Definition of a uniquely identifiable rotation matrix named by DDName name in the GEANT3 style. More...
std::unique_ptr< DDRotation >	DDrotPtr (const DDName &ddname, std::unique_ptr< DDRotationMatrix > rot)
DDRotation	DDrotReflect (const DDName &ddname, std::unique_ptr< DDRotationMatrix > rot)
DDRotation	DDrotReflect (const DDName &ddname, double thetaX, double phiX, double thetaY, double phiY, double thetaZ, double phiZ)
	Defines a rotation-reflection in the Geant3 way. More...
std::ostream &	operator<< (std::ostream &os, const DDRotation &r)

Function Documentation

DDanonymousRot()

DDRotation DDanonymousRot

(

std::unique_ptr< DDRotationMatrix >

rot

)

Defines a anonymous rotation or rotation-reflection matrix.

It can't be addressed by a unique DDName. Once created, it's the users responsibility to keep the reference object DDRotation! Will be mostly used by algorithmic positioning.

Definition at line 143 of file DDRotation.cc.

References eostools::move(), and makeMuonMisalignmentScenario::rot.

Referenced by DDTIBLayerAlgo::execute().

{ return DDRotation(std::move(rot)); }

DDcreateRotationMatrix()

std::unique_ptr<DDRotationMatrix> DDcreateRotationMatrix	(	double	thetaX,
		double	phiX,
		double	thetaY,
		double	phiY,
		double	thetaZ,
		double	phiZ
	)

create a new DDRotationMatrix in the GEANT3 style.

The Matrix must be orthonormal - left or right handed - otherwise a DDException is thrown; memory of the returned pointer belongs to the caller

Definition at line 120 of file DDRotation.cc.

References RPCNoise_example::check, angle_units::operators::convertRadToDeg(), funct::cos(), Exception, EcalTangentSkim_cfg::o, funct::sin(), DOFs::thetaX, DOFs::thetaY, DOFs::thetaZ, x, relval_data_highstats::y, and ctppsOpticalFunctions_non_DB_cff::z.

Referenced by DDHGCalNoTaperEndcap::createQuarter(), DDAngular::execute(), DDMTDLinear::execute(), and DDTIBLayerAlgo::execute().

                                                                                        {
  // define 3 unit std::vectors forming the new left-handed axes
  DD3Vector x(cos(phiX) * sin(thetaX), sin(phiX) * sin(thetaX), cos(thetaX));
  DD3Vector y(cos(phiY) * sin(thetaY), sin(phiY) * sin(thetaY), cos(thetaY));
  DD3Vector z(cos(phiZ) * sin(thetaZ), sin(phiZ) * sin(thetaZ), cos(thetaZ));

  double tol = 1.0e-3;                 // Geant4 compatible
  double check = (x.Cross(y)).Dot(z);  // in case of a LEFT-handed orthogonal system this must be -1, RIGHT-handed: +1
  if ((1. - fabs(check)) > tol) {
    std::ostringstream o;
    o << "matrix is not an (left or right handed) orthonormal matrix! (in deg)" << std::endl
      << " thetaX=" << convertRadToDeg(thetaX) << " phiX=" << convertRadToDeg(phiX) << std::endl
      << " thetaY=" << convertRadToDeg(thetaY) << " phiY=" << convertRadToDeg(phiY) << std::endl
      << " thetaZ=" << convertRadToDeg(thetaZ) << " phiZ=" << convertRadToDeg(phiZ) << std::endl;
    edm::LogError("DDRotation") << o.str() << std::endl;

    throw cms::Exception("DDException") << o.str();
  }

  return std::make_unique<DDRotationMatrix>(x.x(), y.x(), z.x(), x.y(), y.y(), z.y(), x.z(), y.z(), z.z());
}

DDrot() [1/2]

DDRotation DDrot	(	const DDName &	ddname,
		std::unique_ptr< DDRotationMatrix >	rot
	)

Definition of a uniquely identifiable rotation matrix named by DDName name.

DDrot() returns a reference-object DDRotation representing the rotation matrix rot.

The user must not free memory allocated for rot!

Definition at line 67 of file DDRotation.cc.

References eostools::move(), and makeMuonMisalignmentScenario::rot.

Referenced by DDPixFwdBlades::computeNippleParameters(), DDHCalBarrelAlgo::constructGeneralVolume(), DDHCalEndcapAlgo::constructGeneralVolume(), DDHGCalNoTaperEndcap::createQuarter(), DDLinear::execute(), DDTIDModulePosAlgo::execute(), DDAngular::execute(), DDTotemAngular::execute(), DDTrackerAngular::execute(), DDTrackerPhiAlgo::execute(), DDTrackerPhiAltAlgo::execute(), DDHGCalWafer::execute(), DDMTDLinear::execute(), DDTIDRingAlgo::execute(), DDPixFwdDiskAlgo::execute(), DDTECOptoHybAlgo::execute(), DDTECPhiAlgo::execute(), DDTECPhiAltAlgo::execute(), DDTIDAxialCableAlgo::execute(), DDHCalTestBeamAlgo::execute(), DDHCalAngular::execute(), DDPixPhase1FwdDiskAlgo::execute(), DDHCalXtalAlgo::execute(), DDHGCalWaferAlgo::execute(), DDHCalTBZposAlgo::execute(), DDPixFwdRotation::execute(), DDGEMAngular::execute(), DDMuonAngular::execute(), DDTIBLayerAlgo::execute(), DDPixBarLayerAlgo::execute(), DDPixBarLayerUpgradeAlgo::execute(), DDHCalTBCableAlgo::execute(), DDHCalFibreBundle::execute(), DDPixFwdBladesNew::execute(), DDTrackerRingAlgo::execute(), DDTrackerIrregularRingAlgo::execute(), DDPixFwdBlades::execute(), DDDividedTubsPhi::makeDDRotation(), DDDividedConsPhi::makeDDRotation(), DDDividedPolyconePhi::makeDDRotation(), DDDividedPolyhedraPhi::makeDDRotation(), DDTBH4Algo::myrot(), DDEcalBarrelNewAlgo::myrot(), DDEcalBarrelAlgo::myrot(), DDEcalEndcapAlgo::myrot(), DDLRotationSequence::processElement(), DDLRotationByAxis::processElement(), DDLRotationAndReflection::processElement(), regressionTest_first(), and regressionTest_setup().

                                                                            {
  // memory of rot goes sto DDRotationImpl!!
  return DDRotation(ddname, std::move(rot));
}

DDrot() [2/2]

DDRotation DDrot	(	const DDName &	name,
		double	thetaX,
		double	phiX,
		double	thetaY,
		double	phiY,
		double	thetaZ,
		double	phiZ
	)

Definition of a uniquely identifiable rotation matrix named by DDName name in the GEANT3 style.

DDrot() returns a reference-object DDRotation representing the rotation matrix.

Definition at line 78 of file DDRotation.cc.

References RPCNoise_example::check, funct::cos(), Exception, DDName::name(), funct::sin(), DOFs::thetaX, DOFs::thetaY, DOFs::thetaZ, x, relval_data_highstats::y, and ctppsOpticalFunctions_non_DB_cff::z.

                                                                                                              {
  // define 3 unit std::vectors
  DD3Vector x(cos(phiX) * sin(thetaX), sin(phiX) * sin(thetaX), cos(thetaX));
  DD3Vector y(cos(phiY) * sin(thetaY), sin(phiY) * sin(thetaY), cos(thetaY));
  DD3Vector z(cos(phiZ) * sin(thetaZ), sin(phiZ) * sin(thetaZ), cos(thetaZ));

  double tol = 1.0e-3;                 // Geant4 compatible
  double check = (x.Cross(y)).Dot(z);  // in case of a LEFT-handed orthogonal system this must be -1
  if (fabs(1. - check) > tol) {
    edm::LogError("DDRotation") << ddname << " is not a RIGHT-handed orthonormal matrix!" << std::endl;
    throw cms::Exception("DDException") << ddname.name() << " is not RIGHT-handed!";
  }

  return DDRotation(ddname,
                    std::make_unique<DDRotationMatrix>(x.x(), y.x(), z.x(), x.y(), y.y(), z.y(), x.z(), y.z(), z.z()));
}

DDrotPtr()

std::unique_ptr<DDRotation> DDrotPtr	(	const DDName &	ddname,
		std::unique_ptr< DDRotationMatrix >	rot
	)

Definition at line 72 of file DDRotation.cc.

References eostools::move(), and makeMuonMisalignmentScenario::rot.

Referenced by DDLPosPart::processElement().

                                                                                              {
  // memory of rot goes sto DDRotationImpl!!
  return std::make_unique<DDRotation>(ddname, std::move(rot));
}

DDrotReflect() [1/2]

DDRotation DDrotReflect	(	const DDName &	ddname,
		std::unique_ptr< DDRotationMatrix >	rot
	)

Definition at line 96 of file DDRotation.cc.

References eostools::move(), and makeMuonMisalignmentScenario::rot.

                                                                                   {
  return DDRotation(ddname, std::move(rot));
}

DDrotReflect() [2/2]

DDRotation DDrotReflect	(	const DDName &	ddname,
		double	thetaX,
		double	phiX,
		double	thetaY,
		double	phiY,
		double	thetaZ,
		double	phiZ
	)

Defines a rotation-reflection in the Geant3 way.

The resulting matrix MUST be a LEFThanded orthonormal system, otherwise a DDException will be thrown!

Definition at line 101 of file DDRotation.cc.

References RPCNoise_example::check, funct::cos(), Exception, DDName::name(), funct::sin(), DOFs::thetaX, DOFs::thetaY, DOFs::thetaZ, x, relval_data_highstats::y, and ctppsOpticalFunctions_non_DB_cff::z.

Referenced by DDLRotationAndReflection::processElement().

                                                                                                              {
  // define 3 unit std::vectors forming the new left-handed axes
  DD3Vector x(cos(phiX) * sin(thetaX), sin(phiX) * sin(thetaX), cos(thetaX));
  DD3Vector y(cos(phiY) * sin(thetaY), sin(phiY) * sin(thetaY), cos(thetaY));
  DD3Vector z(cos(phiZ) * sin(thetaZ), sin(phiZ) * sin(thetaZ), cos(thetaZ));

  double tol = 1.0e-3;                 // Geant4 compatible
  double check = (x.Cross(y)).Dot(z);  // in case of a LEFT-handed orthogonal system this must be -1
  if (fabs(1. + check) > tol) {
    edm::LogError("DDRotation") << ddname << " is not a LEFT-handed orthonormal matrix!" << std::endl;
    throw cms::Exception("DDException") << ddname.name() << " is not LEFT-handed!";
  }

  return DDRotation(ddname,
                    std::make_unique<DDRotationMatrix>(x.x(), y.x(), z.x(), x.y(), y.y(), z.y(), x.z(), y.z(), z.z()));
}

operator<<()

std::ostream& operator<<	(	std::ostream &	os,
		const DDRotation &	r
	)

Definition at line 21 of file DDRotation.cc.

References angle_units::operators::convertRadToDeg(), l1tGTMenu_BTagSeeds_cff::os, and eostools::rm().

                                                            {
  DDBase<DDName, DDRotationMatrix*>::def_type defined(r.isDefined());
  if (defined.first) {
    os << *(defined.first) << " ";
    if (defined.second) {
      const DDRotationMatrix& rm = r.rotation();
      DDAxisAngle ra(rm);
      os << "t=" << convertRadToDeg(ra.Axis().Theta()) << "deg "
         << "p=" << convertRadToDeg(ra.Axis().Phi()) << "deg "
         << "a=" << convertRadToDeg(ra.Angle()) << "deg";
    } else {
      os << "* rotation not defined * ";
    }
  } else {
    os << "* rotation not declared * ";
  }
  return os;
}