DDEcalEndcapTrap Class Reference

#include <DDEcalEndcapTrap.h>

Public Member Functions
void	calculateCentres ()
void	calculateCorners ()
DDTranslation	centrePos ()
DDTranslation	cornerPos (const int icorner)
void	cornerPos (const int icorner, const DDTranslation &cc)
	DDEcalEndcapTrap ()=delete
	DDEcalEndcapTrap (const int hand, const double front, const double rear, const double length)
double	elevationAngle ()
double	elevationAngle (const DDTranslation &trans)
DDTranslation	fcentrePos ()
void	moveto (const DDTranslation &frontCentre, const DDTranslation &rearCentre)
double	polarAngle ()
double	polarAngle (const DDTranslation &trans)
void	print ()
DDTranslation	rcentrePos ()
void	rotate (const DDRotationMatrix &rot)
void	rotate (const DDTranslation &frontCentre, const DDTranslation &rearCentre)
void	rotateX (const double angle)
void	rotateY (const double angle)
DDRotationMatrix	rotation ()
void	translate ()
void	translate (const DDTranslation &trans)

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 DDEcalEndcapTrap.h.

Constructor & Destructor Documentation

DDEcalEndcapTrap() [1/2]

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

Definition at line 12 of file DDEcalEndcapTrap.cc.

                                                                                                             {
  //
  //  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();
}

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

DDEcalEndcapTrap() [2/2]

DDEcalEndcapTrap::DDEcalEndcapTrap ( )

delete

Member Function Documentation

calculateCentres()

void DDEcalEndcapTrap::calculateCentres

(

)

Definition at line 176 of file DDEcalEndcapTrap.cc.

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

References m_centre, m_corners, m_fcentre, and m_rcentre.

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

calculateCorners()

void DDEcalEndcapTrap::calculateCorners

(

)

centrePos()

DDTranslation DDEcalEndcapTrap::centrePos

(

)

Definition at line 220 of file DDEcalEndcapTrap.cc.

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

References m_centre.

Referenced by algorithm(), DDEcalEndcapAlgo::EEPositionCRs(), DDEcalEndcapAlgo::execute(), moveto(), print(), and translate().

cornerPos() [1/2]

DDTranslation DDEcalEndcapTrap::cornerPos

(

const int

icorner

)

Definition at line 203 of file DDEcalEndcapTrap.cc.

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

References m_corners.

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

cornerPos() [2/2]

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

Definition at line 210 of file DDEcalEndcapTrap.cc.

                                                                                  {
  //
  //  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()));
    ;
  }
}

References m_corners.

elevationAngle() [1/2]

double DDEcalEndcapTrap::elevationAngle

(

)

Definition at line 249 of file DDEcalEndcapTrap.cc.

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

References fcentrePos(), and rcentrePos().

Referenced by moveto().

elevationAngle() [2/2]

double DDEcalEndcapTrap::elevationAngle

(

const DDTranslation &

trans

)

Definition at line 241 of file DDEcalEndcapTrap.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);
}

fcentrePos()

DDTranslation DDEcalEndcapTrap::fcentrePos

(

)

Definition at line 227 of file DDEcalEndcapTrap.cc.

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

References m_fcentre.

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

moveto()

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

Definition at line 132 of file DDEcalEndcapTrap.cc.

                                                                                               {
  //
  //  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)
  // edm::LogInfo("EcalGeom") << "moveto: frontCentre " << frontCentre << std::endl;
  // edm::LogInfo("EcalGeom") << "moveto: rearCentre  " << rearCentre << std::endl;
  // edm::LogInfo("EcalGeom") << "moveto: X rotation: " << targetTheta << " " << currentTheta << " " << targetTheta-currentTheta << std::endl;
  // edm::LogInfo("EcalGeom") << "moveto: Y rotation: " << targetPhi << " " << currentPhi << " " << " " << targetPhi-currentPhi << std::endl;
  rotateX(targetTheta - currentTheta);
  rotateY(targetPhi - currentPhi);
 
  //  Translate SC to final position
  DDTranslation targetCentre = 0.5 * (frontCentre + rearCentre);
  // edm::LogInfo("EcalGeom") << "moveto: translation " << targetCentre-centrePos() << std::endl;
  translate(targetCentre - centrePos());
}

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

Referenced by algorithm(), DDEcalEndcapAlgo::EEPositionCRs(), and DDEcalEndcapAlgo::execute().

polarAngle() [1/2]

double DDEcalEndcapTrap::polarAngle

(

)

Definition at line 265 of file DDEcalEndcapTrap.cc.

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

References fcentrePos(), and rcentrePos().

Referenced by moveto().

polarAngle() [2/2]

double DDEcalEndcapTrap::polarAngle

(

const DDTranslation &

trans

)

Definition at line 257 of file DDEcalEndcapTrap.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);
}

print()

void DDEcalEndcapTrap::print

(

void

)

Definition at line 273 of file DDEcalEndcapTrap.cc.

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

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

rcentrePos()

DDTranslation DDEcalEndcapTrap::rcentrePos

(

)

Definition at line 234 of file DDEcalEndcapTrap.cc.

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

References m_rcentre.

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

rotate() [1/2]

void DDEcalEndcapTrap::rotate

(

const DDRotationMatrix &

rot

)

Definition at line 94 of file DDEcalEndcapTrap.cc.

                                                         {
  //
  //  Rotate supercrystal by specified rotation about (0,0,0)
  //
 
  int icorner;
  DDTranslation cc;
  //  edm::LogInfo("EcalGeom") << "DDEcalEndcapTrap::rotate - rotation " << rot << std::endl;
  for (icorner = 1; icorner <= 8; icorner++) {
    cc = cornerPos(icorner);
    //     edm::LogInfo("EcalGeom") << "   Corner (orig) " << icorner << cc << std::endl;
    cc = rot * cc;
    //     edm::LogInfo("EcalGeom") << "   Corner (rot)  " << icorner << cc << std::endl;
    cornerPos(icorner, cc);
  }
  m_rotation = rot * m_rotation;
  calculateCentres();
}

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

Referenced by rotateX(), and rotateY().

rotate() [2/2]

void DDEcalEndcapTrap::rotate	(	const DDTranslation &	frontCentre,
		const DDTranslation &	rearCentre
	)

Definition at line 86 of file DDEcalEndcapTrap.cc.

                                                                                               {
  //
  //  Rotate supercrystal to bring front and rear face centres to specified points
  //
  edm::LogInfo("EcalGeom") << "DDEcalEndcapTrap::rotate(DDTranslation,DDTranslation) - not yet implemented"
                           << std::endl;
}

rotateX()

void DDEcalEndcapTrap::rotateX

(

const double

angle

)

Definition at line 157 of file DDEcalEndcapTrap.cc.

                                                 {
  //
  //  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()));
}

References angle(), rotate(), and createJobs::tmp.

Referenced by moveto().

rotateY()

void DDEcalEndcapTrap::rotateY

(

const double

angle

)

Definition at line 167 of file DDEcalEndcapTrap.cc.

                                                 {
  //
  //  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()));
}

References angle(), rotate(), and createJobs::tmp.

Referenced by moveto().

rotation()

DDRotationMatrix DDEcalEndcapTrap::rotation ( void  )

inline

Definition at line 35 of file DDEcalEndcapTrap.h.

{ return m_rotation; }

References m_rotation.

Referenced by algorithm(), DDEcalEndcapAlgo::EEPositionCRs(), and DDEcalEndcapAlgo::execute().

translate() [1/2]

void DDEcalEndcapTrap::translate

(

)

Definition at line 113 of file DDEcalEndcapTrap.cc.

                                 {
  //  edm::LogInfo("EcalGeom") << "DDEcalEndcapTrap::translate() not yet implemented" << std::endl;
  translate(-1. * centrePos());
}

References centrePos().

Referenced by DDEcalEndcapTrap(), and moveto().

translate() [2/2]

void DDEcalEndcapTrap::translate

(

const DDTranslation &

trans

)

Definition at line 118 of file DDEcalEndcapTrap.cc.

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

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

Referenced by algorithm(), and DDEcalEndcapAlgo::execute().

Member Data Documentation

m_centre

double DDEcalEndcapTrap::m_centre[4]

private

Definition at line 42 of file DDEcalEndcapTrap.h.

Referenced by calculateCentres(), and centrePos().

m_corners

double DDEcalEndcapTrap::m_corners[25]

private

Definition at line 45 of file DDEcalEndcapTrap.h.

Referenced by calculateCentres(), cornerPos(), and DDEcalEndcapTrap().

m_fcentre

double DDEcalEndcapTrap::m_fcentre[4]

private

Definition at line 43 of file DDEcalEndcapTrap.h.

Referenced by calculateCentres(), and fcentrePos().

m_front

double DDEcalEndcapTrap::m_front

private

Definition at line 46 of file DDEcalEndcapTrap.h.

Referenced by DDEcalEndcapTrap().

m_hand

int DDEcalEndcapTrap::m_hand

private

Definition at line 50 of file DDEcalEndcapTrap.h.

Referenced by DDEcalEndcapTrap().

m_length

double DDEcalEndcapTrap::m_length

private

Definition at line 48 of file DDEcalEndcapTrap.h.

Referenced by DDEcalEndcapTrap().

m_rcentre

double DDEcalEndcapTrap::m_rcentre[4]

private

Definition at line 44 of file DDEcalEndcapTrap.h.

Referenced by calculateCentres(), and rcentrePos().

m_rear

double DDEcalEndcapTrap::m_rear

private

Definition at line 47 of file DDEcalEndcapTrap.h.

Referenced by DDEcalEndcapTrap().

m_rotation

DDRotationMatrix DDEcalEndcapTrap::m_rotation

private

Definition at line 39 of file DDEcalEndcapTrap.h.

Referenced by rotate(), and rotation().

m_translation

DDTranslation DDEcalEndcapTrap::m_translation

private

Definition at line 40 of file DDEcalEndcapTrap.h.

Referenced by translate().

m_update

int DDEcalEndcapTrap::m_update

private

Definition at line 51 of file DDEcalEndcapTrap.h.