#include <DDEcalEndcapTrapX.h>

Public Member Functions
void	calculateCentres ()

void	calculateCorners ()

DDTranslation	centrePos ()

DDTranslation	cornerPos (const int icorner)

void	cornerPos (const int icorner, const DDTranslation &cc)

	DDEcalEndcapTrapX (const int hand, const double front, const double rear, const double length)

	DDEcalEndcapTrapX ()=delete

double	elevationAngle (const DDTranslation &trans)

double	elevationAngle ()

DDTranslation	fcentrePos ()

void	moveto (const DDTranslation &frontCentre, const DDTranslation &rearCentre)

double	polarAngle (const DDTranslation &trans)

double	polarAngle ()

void	print ()

DDTranslation	rcentrePos ()

void	rotate (const DDRotationMatrix &rot)

void	rotateX (const double angle)

void	rotateY (const double angle)

DDRotationMatrix	rotation ()

void	translate (const DDTranslation &trans)

void	translate ()

Private Attributes
double	m_centre [4]

double	m_corners [25]

double	m_fcentre [4]

double	m_front

int	m_hand

double	m_length

double	m_rcentre [4]

double	m_rear

DDRotationMatrix	m_rotation

DDTranslation	m_translation

int	m_update

Detailed Description

Definition at line 11 of file DDEcalEndcapTrapX.h.

Constructor & Destructor Documentation

DDEcalEndcapTrapX::DDEcalEndcapTrapX	(	const int	hand,
		const double	front,
		const double	rear,
		const double	length
	)

Definition at line 7 of file DDEcalEndcapTrapX.cc.

References calculateCentres(), m_corners, m_front, m_hand, m_length, m_rear, and translate().

                                                                                                                {
   //
   //  Initialise corners of supercrystal.
 
   // Start out with bottom surface on (x,z) plane, front face in (x,y) plane.
 
   double xsign;
 
   if (hand == 2) {
     xsign = -1.;
   } else {
     xsign = 1.;
   }
 
   m_hand = hand;
   m_front = front;
   m_rear = rear;
   m_length = length;
 
   int icorner;
   icorner = 1;
   m_corners[3 * icorner - 3] = xsign * front;
   m_corners[3 * icorner - 2] = front;
   m_corners[3 * icorner - 1] = 0.;
   icorner = 2;
   m_corners[3 * icorner - 3] = xsign * front;
   m_corners[3 * icorner - 2] = 0.;
   m_corners[3 * icorner - 1] = 0.;
   icorner = 3;
   m_corners[3 * icorner - 3] = 0.;
   m_corners[3 * icorner - 2] = 0.;
   m_corners[3 * icorner - 1] = 0.;
   icorner = 4;
   m_corners[3 * icorner - 3] = 0.;
   m_corners[3 * icorner - 2] = front;
   m_corners[3 * icorner - 1] = 0.;
 
   icorner = 5;
   m_corners[3 * icorner - 3] = xsign * rear;
   m_corners[3 * icorner - 2] = rear;
   m_corners[3 * icorner - 1] = length;
   icorner = 6;
   m_corners[3 * icorner - 3] = xsign * rear;
   m_corners[3 * icorner - 2] = 0.;
   m_corners[3 * icorner - 1] = length;
   icorner = 7;
   m_corners[3 * icorner - 3] = 0.;
   m_corners[3 * icorner - 2] = 0.;
   m_corners[3 * icorner - 1] = length;
   icorner = 8;
   m_corners[3 * icorner - 3] = 0.;
   m_corners[3 * icorner - 2] = rear;
   m_corners[3 * icorner - 1] = length;
 
   calculateCentres();
 
   // Move centre of SC to (0,0,0)
 
   translate();
 
   // Rotate into standard position (face centres on z axis)
 
   //  this->rotate();
 
   calculateCentres();
 }

DDEcalEndcapTrapX::DDEcalEndcapTrapX ( )

delete

Member Function Documentation

void DDEcalEndcapTrapX::calculateCentres ( )

Definition at line 170 of file DDEcalEndcapTrapX.cc.

References m_centre, m_corners, m_fcentre, and m_rcentre.

Referenced by DDEcalEndcapTrapX(), rotate(), and translate().

                                          {
   //
   //  Calculate crystal centre and front & rear face centres
   //
 
   int ixyz, icorner;
 
   for (ixyz = 0; ixyz < 3; ixyz++) {
     m_centre[ixyz] = 0;
     m_fcentre[ixyz] = 0;
     m_rcentre[ixyz] = 0;
   }
 
   for (icorner = 1; icorner <= 4; icorner++) {
     for (ixyz = 0; ixyz < 3; ixyz++) {
       m_centre[ixyz] = m_centre[ixyz] + 0.125 * m_corners[3 * icorner - 3 + ixyz];
       m_fcentre[ixyz] = m_fcentre[ixyz] + 0.25 * m_corners[3 * icorner - 3 + ixyz];
     }
   }
   for (icorner = 5; icorner <= 8; icorner++) {
     for (ixyz = 0; ixyz < 3; ixyz++) {
       m_centre[ixyz] = m_centre[ixyz] + 0.125 * m_corners[3 * icorner - 3 + ixyz];
       m_rcentre[ixyz] = m_rcentre[ixyz] + 0.25 * m_corners[3 * icorner - 3 + ixyz];
     }
   }
 }

void DDEcalEndcapTrapX::calculateCorners ( )

DDTranslation DDEcalEndcapTrapX::centrePos ( )

Definition at line 214 of file DDEcalEndcapTrapX.cc.

References m_centre.

Referenced by algorithm(), moveto(), print(), and translate().

                                            {
   //
   //  Return SC centre as a Translation
   //
   return DDTranslation(m_centre[0], m_centre[1], m_centre[2]);
 }

DDTranslation DDEcalEndcapTrapX::cornerPos ( const int icorner )

Definition at line 197 of file DDEcalEndcapTrapX.cc.

References m_corners.

Referenced by print(), rotate(), and translate().

                                                             {
   //
   //  Return specified corner as a Translation
   //
   return DDTranslation(m_corners[3 * icorner - 3], m_corners[3 * icorner - 2], m_corners[3 * icorner - 1]);
 }

void DDEcalEndcapTrapX::cornerPos	(	const int	icorner,
		const DDTranslation &	cc
	)

Definition at line 204 of file DDEcalEndcapTrapX.cc.

References m_corners.

                                                                                    {
   //
   //  Save position of specified corner.
   //
   for (int ixyz = 0; ixyz < 3; ixyz++) {
     m_corners[3 * icorner - 3 + ixyz] = (0 == ixyz ? cornerxyz.x() : (1 == ixyz ? cornerxyz.y() : cornerxyz.z()));
     ;
   }
 }

double DDEcalEndcapTrapX::elevationAngle ( const DDTranslation & trans )

Definition at line 235 of file DDEcalEndcapTrapX.cc.

                                                                    {
   //
   //  Return elevation angle (out of x-z plane) of a given translation (seen as a vector from the origin).
   //
   double sintheta = trans.y() / trans.r();
   return asin(sintheta);
 }

double DDEcalEndcapTrapX::elevationAngle ( )

Definition at line 243 of file DDEcalEndcapTrapX.cc.

References fcentrePos(), and rcentrePos().

Referenced by moveto().

                                          {
   //
   //  Return elevation angle (out of x-z plane) of SC in current position.
   //
   DDTranslation current = fcentrePos() - rcentrePos();
   return elevationAngle(current);
 }

DDTranslation DDEcalEndcapTrapX::fcentrePos ( )

Definition at line 221 of file DDEcalEndcapTrapX.cc.

References m_fcentre.

Referenced by elevationAngle(), polarAngle(), and print().

                                             {
   //
   //  Return SC front face centre as a Translation
   //
   return DDTranslation(m_fcentre[0], m_fcentre[1], m_fcentre[2]);
 }

void DDEcalEndcapTrapX::moveto	(	const DDTranslation &	frontCentre,
		const DDTranslation &	rearCentre
	)

Definition at line 120 of file DDEcalEndcapTrapX.cc.

References centrePos(), elevationAngle(), polarAngle(), rotateX(), rotateY(), and translate().

Referenced by algorithm().

                                                                                                 {
   //
   //  Rotate (about X then about Y) and translate supercrystal to bring axis joining front and rear face centres parallel to line connecting specified points
   //
 
   //  Get azimuthal and polar angles of current axis and target axis
   double currentTheta = elevationAngle();
   double currentPhi = polarAngle();
   double targetTheta = elevationAngle(frontCentre - rearCentre);
   double targetPhi = polarAngle(frontCentre - rearCentre);
 
   //  Rotate to correct angle (X then Y)
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalGeom") << "moveto: frontCentre " << frontCentre << std::endl
                                << "moveto: rearCentre  " << rearCentre << std::endl
                                << "moveto: X rotation: " << targetTheta << " " << currentTheta << " "
                                << targetTheta - currentTheta << std::endl
                                << "moveto: Y rotation: " << targetPhi << " " << currentPhi << " "
                                << " " << targetPhi - currentPhi << std::endl;
 #endif
   rotateX(targetTheta - currentTheta);
   rotateY(targetPhi - currentPhi);
 
   //  Translate SC to final position
   DDTranslation targetCentre = 0.5 * (frontCentre + rearCentre);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalGeom") << "moveto: translation " << targetCentre - centrePos();
 #endif
   translate(targetCentre - centrePos());
 }

double DDEcalEndcapTrapX::polarAngle ( const DDTranslation & trans )

Definition at line 251 of file DDEcalEndcapTrapX.cc.

                                                                {
   //
   //  Return polar angle (from x to z) of a given translation (seen as a vector from the origin).
   //
   double tanphi = trans.x() / trans.z();
   return atan(tanphi);
 }

double DDEcalEndcapTrapX::polarAngle ( )

Definition at line 259 of file DDEcalEndcapTrapX.cc.

References fcentrePos(), and rcentrePos().

Referenced by moveto().

                                      {
   //
   //  Return elevation angle (out of x-z plane) of SC in current position.
   //
   DDTranslation current = fcentrePos() - rcentrePos();
   return polarAngle(current);
 }

void DDEcalEndcapTrapX::print ( void )

Definition at line 267 of file DDEcalEndcapTrapX.cc.

References centrePos(), cornerPos(), fcentrePos(), and rcentrePos().

                               {
   //
   //  Print SC coordinates for debugging
   //
   edm::LogVerbatim("EcalGeom") << "Endcap supercrystal";
   for (int ic = 1; ic <= 8; ic++) {
     DDTranslation cc = cornerPos(ic);
     edm::LogVerbatim("EcalGeom") << "Corner " << ic << " " << cc << std::endl;
   }
   edm::LogVerbatim("EcalGeom") << "    Centre " << centrePos() << std::endl
                                << "   fCentre " << fcentrePos() << std::endl
                                << "   rCentre " << rcentrePos();
 }

DDTranslation DDEcalEndcapTrapX::rcentrePos ( )

Definition at line 228 of file DDEcalEndcapTrapX.cc.

References m_rcentre.

Referenced by elevationAngle(), polarAngle(), and print().

                                             {
   //
   //  Return SC rear face centre as a Translation
   //
   return DDTranslation(m_rcentre[0], m_rcentre[1], m_rcentre[2]);
 }

void DDEcalEndcapTrapX::rotate ( const DDRotationMatrix & rot )

Definition at line 74 of file DDEcalEndcapTrapX.cc.

References calculateCentres(), cornerPos(), m_rotation, and makeMuonMisalignmentScenario::rot.

Referenced by rotateX(), and rotateY().

                                                           {
   //
   //  Rotate supercrystal by specified rotation about (0,0,0)
   //
 
   int icorner;
   DDTranslation cc;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalGeom") << "DDEcalEndcapTrapX::rotate - rotation " << rot;
 #endif
   for (icorner = 1; icorner <= 8; icorner++) {
     cc = cornerPos(icorner);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("EcalGeom") << "   Corner (orig) " << icorner << cc;
 #endif
     cc = rot * cc;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("EcalGeom") << "   Corner (rot)  " << icorner << cc;
 #endif
     cornerPos(icorner, cc);
   }
   m_rotation = rot * m_rotation;
   calculateCentres();
 }

void DDEcalEndcapTrapX::rotateX ( const double angle )

Definition at line 151 of file DDEcalEndcapTrapX.cc.

References rotate(), and createJobs::tmp.

Referenced by moveto().

                                                   {
   //
   //  Rotate SC through given angle about X axis
   //
 
   const CLHEP::HepRotation tmp(CLHEP::Hep3Vector(1., 0., 0.), angle);
 
   rotate(DDRotationMatrix(tmp.xx(), tmp.xy(), tmp.xz(), tmp.yx(), tmp.yy(), tmp.yz(), tmp.zx(), tmp.zy(), tmp.zz()));
 }

void DDEcalEndcapTrapX::rotateY ( const double angle )

Definition at line 161 of file DDEcalEndcapTrapX.cc.

References rotate(), and createJobs::tmp.

Referenced by moveto().

                                                   {
   //
   //  Rotate SC through given angle about Y axis
   //
   const CLHEP::HepRotation tmp(CLHEP::Hep3Vector(0., 1., 0.), angle);
 
   rotate(DDRotationMatrix(tmp.xx(), tmp.xy(), tmp.xz(), tmp.yx(), tmp.yy(), tmp.yz(), tmp.zx(), tmp.zy(), tmp.zz()));
 }

DDRotationMatrix DDEcalEndcapTrapX::rotation ( void )

inline

Definition at line 34 of file DDEcalEndcapTrapX.h.

References m_rotation.

Referenced by algorithm().

34 { return m_rotation; }

DDEcalEndcapTrapX::m_rotation

DDRotationMatrix m_rotation

Definition: DDEcalEndcapTrapX.h:38

void DDEcalEndcapTrapX::translate ( const DDTranslation & trans )

Definition at line 106 of file DDEcalEndcapTrapX.cc.

References calculateCentres(), cornerPos(), and m_translation.

Referenced by algorithm().

                                                             {
   //
   //  Translate supercrystal by specified amount
   //
 
   DDTranslation tcorner;
   for (int icorner = 1; icorner <= 8; icorner++) {
     tcorner = cornerPos(icorner) + trans;
     cornerPos(icorner, tcorner);
   }
   calculateCentres();
   m_translation = trans + m_translation;
 }

void DDEcalEndcapTrapX::translate ( )

Definition at line 99 of file DDEcalEndcapTrapX.cc.

References centrePos().

Referenced by DDEcalEndcapTrapX(), and moveto().

                                   {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalGeom") << "DDEcalEndcapTrapX::translate() not yet implemented";
 #endif
   translate(-1. * centrePos());
 }