OptOModifiedRhomboidPrism Class Reference

#include <OptOModifiedRhomboidPrism.h>

Inheritance diagram for OptOModifiedRhomboidPrism:

Public Member Functions
void	constructSolidShape ()
virtual void	detailedDeviatesLightRay (LightRay &lightray)
virtual void	detailedTraversesLightRay (LightRay &lightray)
virtual void	fastDeviatesLightRay (LightRay &lightray)
virtual void	fastTraversesLightRay (LightRay &lightray)
ALIPlane	getRotatedPlate (const ALIbool forwardPlate)
	OptOModifiedRhomboidPrism ()
	OptOModifiedRhomboidPrism (OpticalObject *parent, const ALIstring &type, const ALIstring &name, const ALIbool copy_data)
	~OptOModifiedRhomboidPrism ()
Public Member Functions inherited from OpticalObject
void	addCoordinateEntryToList (Entry *entry)
void	addExtraEntryToList (Entry *entry)
void	addExtraEntryValueOriginalOriginalToList (ALIdouble entry_value)
void	addExtraEntryValueOriginalToList (ALIdouble entry_value)
void	addExtraEntryValueToList (ALIdouble entry_value)
double	addPii (double val)
double	approxTo0 (double val)
CLHEP::HepRotation	buildRmFromEntryValuesOriginalOriginal ()
const CLHEP::Hep3Vector &	centreGlob () const
const CLHEP::Hep3Vector &	centreGlobal () const
const CLHEP::Hep3Vector &	centreGlobOriginal () const
const CLHEP::Hep3Vector &	centreGlobOriginalOriginal () const
const CLHEP::Hep3Vector	centreLocal () const
int	checkMatrixEquations (double angleX, double angleY, double angleZ, CLHEP::HepRotation *rot=0)
void	construct ()
void	constructFromOptAligInfo (const OpticalAlignInfo &oaInfo)
virtual void	constructMaterial ()
const std::vector< Entry * > &	CoordinateEntryList () const
void	createComponentOptOsFromOptAlignInfo ()
virtual void	defaultBehaviour (LightRay &lightray, Measurement &meas)
double	diff2pi (double ang1, double ang2)
void	displaceCentreGlob (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlob (const CLHEP::Hep3Vector &dispVec)
void	displaceCentreGlobOriginal (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlobOriginal (const CLHEP::Hep3Vector &dispVec)
void	displaceCentreGlobOriginalOriginal (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlobOriginalOriginal (const CLHEP::Hep3Vector &dispVec)
void	displaceExtraEntry (const ALIuint entryNo, const ALIdouble disp)
void	displaceExtraEntryOriginal (const ALIuint entryNo, const ALIdouble disp)
void	displaceExtraEntryOriginalOriginal (const ALIuint entryNo, const ALIdouble disp)
void	displaceRmGlobAroundGlobal (OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobAroundLocal (OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobOriginal (const OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobOriginalOriginal (const OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
bool	eq2ang (double ang1, double ang2)
const std::vector< Entry * > &	ExtraEntryList () const
const ALIint	extraEntryNo (const ALIstring &entry_name) const
std::vector< ALIdouble > &	ExtraEntryValueList ()
const std::vector< ALIdouble > &	ExtraEntryValueOriginalList ()
const std::vector< ALIdouble > &	ExtraEntryValueOriginalOriginalList ()
virtual void	fillIguana ()
virtual void	fillVRML ()
const ALIdouble	findExtraEntryValue (const ALIstring &eename) const
const ALIbool	findExtraEntryValueIfExists (const ALIstring &eename, ALIdouble &value) const
const ALIdouble	findExtraEntryValueMustExist (const ALIstring &eename) const
const ALIuint	getCmsswID () const
std::vector< ALIstring >	getCoordinateFromOptAlignParam (const OpticalAlignParam &oaParam)
CLHEP::Hep3Vector	getDisplacementInLocalCoordinates (const XYZcoor coor, const ALIdouble disp)
const double	getEntryCentre (const XYZcoor coor) const
const double	getEntryCentre (const ALIstring &coor) const
const double	getEntryRMangle (const XYZcoor coor) const
const double	getEntryRMangle (const ALIstring &coor) const
std::vector< double >	getLocalRotationAngles (const std::vector< Entry * > &entries) const
CocoaMaterialElementary *	getMaterial () const
ALIPlane	getPlate (const ALIbool forwardPlate, const ALIbool applyWedge)
std::vector< double >	getRotationAnglesFromMatrix (CLHEP::HepRotation &rmLocal, const std::vector< Entry * > &entries) const
std::vector< double >	getRotationAnglesInOptOFrame (const OpticalObject *optoAncestor, const std::vector< Entry * > &entries) const
CocoaSolidShape *	getSolidShape () const
CLHEP::Hep3Vector	getZAxis ()
const ALIuint	ID () const
const ALIstring	longName () const
virtual void	makeMeasurement (LightRay &lightray, Measurement &meas)
Measurement *	meas ()
const ALIstring &	name () const
	OpticalObject ()
	OpticalObject (OpticalObject *parent, const ALIstring &type, const ALIstring &name, const ALIbool copy_data)
const OpticalObject *	parent () const
virtual void	participateInMeasurement (LightRay &lightray, Measurement &meas, const ALIstring &behav)
void	propagateGlobalRMOriginalOriginalChangeToChildren (const CLHEP::HepRotation &rmorioriold, const CLHEP::HepRotation &rmoriorinew)
void	resetGlobalCoordinates ()
void	resetOriginalOriginalCoordinates ()
const CLHEP::HepRotation &	rmGlob () const
const CLHEP::HepRotation &	rmGlobOriginal () const
const CLHEP::HepRotation &	rmGlobOriginalOriginal () const
const CLHEP::HepRotation	rmLocal () const
void	setCmsswID (ALIuint id)
void	setExtraEntryValue (const ALIuint entryNo, const ALIdouble disp)
void	setGlobalCoordinates ()
void	setGlobalRMOriginalOriginal (const CLHEP::HepRotation &rmoriori)
void	setID (ALIuint id)
void	setMeas (Measurement *meas)
	set current measurement More...
void	setOriginalEntryValues ()
void	setRmGlobal (const CLHEP::HepRotation &rm)
void	setRmGlobalOriginal (const CLHEP::HepRotation &rm)
void	setType (const ALIstring &type)
const ALIstring	shortName () const
const ALIstring &	type () const
virtual void	userDefinedBehaviour (LightRay &lightray, Measurement &meas, const ALIstring &behav)
virtual	~OpticalObject ()

Additional Inherited Members
Protected Member Functions inherited from OpticalObject
virtual void	fillExtraEntry (std::vector< ALIstring > &wordlist)
Protected Attributes inherited from OpticalObject
CocoaMaterialElementary *	theMaterial
CocoaSolidShape *	theSolidShape
ALIint	verbose

Detailed Description

Definition at line 18 of file OptOModifiedRhomboidPrism.h.

Constructor & Destructor Documentation

OptOModifiedRhomboidPrism::OptOModifiedRhomboidPrism ( )

inline

Definition at line 23 of file OptOModifiedRhomboidPrism.h.

{ };

OptOModifiedRhomboidPrism::OptOModifiedRhomboidPrism ( OpticalObject *  parent, const ALIstring &  type, const ALIstring &  name, const ALIbool  copy_data  )

inline

Definition at line 24 of file OptOModifiedRhomboidPrism.h.

                                                                                                                          : 
  OpticalObject( parent, type, name, copy_data){ };

OptOModifiedRhomboidPrism::~OptOModifiedRhomboidPrism ( )

inline

Definition at line 26 of file OptOModifiedRhomboidPrism.h.

{ };

Member Function Documentation

void OptOModifiedRhomboidPrism::constructSolidShape ( )

virtual

Reimplemented from OpticalObject.

Definition at line 281 of file OptOModifiedRhomboidPrism.cc.

References GlobalOptionMgr::getGlobalOptionValue(), GlobalOptionMgr::getInstance(), and m.

{
  ALIdouble go;
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  gomgr->getGlobalOptionValue("VisScale", go );

  theSolidShape = new CocoaSolidShapeBox( "Box", go*5.*cm/m, go*5.*cm/m, go*5.*cm/m ); //COCOA internal units are meters
}

void OptOModifiedRhomboidPrism::detailedDeviatesLightRay ( LightRay &  lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 28 of file OptOModifiedRhomboidPrism.cc.

References gather_cfg::cout, ALIUtils::debug, ALIUtils::dump3v(), LightRay::dumpData(), mergeVDriftHistosByStation::name, LightRay::reflect(), and LightRay::refract().

{
  if (ALIUtils::debug >= 2) std::cout << "$$$$$ LR: DETAILED DEVIATION IN MODIFIED RHOMBOID PRISM " << name() << std::endl;

  CLHEP::Hep3Vector XAxis(1.,0.,0.);
  CLHEP::HepRotation rmt = rmGlob();
  XAxis = rmt*XAxis;
  CLHEP::Hep3Vector YAxis(0.,1.,0.);
  YAxis = rmt*YAxis;
  CLHEP::Hep3Vector ZAxis(0.,0.,1.);
  ZAxis = rmt*ZAxis;

  ALIUtils::dump3v( XAxis , " x axis ");
  ALIUtils::dump3v( YAxis , " y axis ");
  ALIUtils::dump3v( ZAxis , " z axis ");
  if (ALIUtils::debug >= 5) {
    ALIUtils::dump3v( centreGlob(), " centre ");
  }

  if (ALIUtils::debug >= 2) std::cout << "$$$ LR: REFRACTION IN FORWARD PLATE " << std::endl;
  //---------- Get forward plate
  ALIPlane plate = getPlate( 1, 1 );
  //---------- Refract in plate while entering
  ALIdouble refra_ind1 = 1.;
  ALIdouble refra_ind2 = findExtraEntryValueMustExist("refra_ind");
  lightray.refract( plate, refra_ind1, refra_ind2 );
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("LightRay after Refraction at entering: "); 
  }

  if (ALIUtils::debug >= 2) std::cout << std::endl << "$$$ LR: REFLECTION IN FIRST PLATE " << std::endl;
  //---------- Get up plate rotated
  plate = getRotatedPlate( 1 );
  //---------- Reflect in plate
  lightray.reflect( plate );

  if (ALIUtils::debug >= 2) std::cout << std::endl << "$$$ LR: REFLECTION IN SECOND PLATE " << std::endl;
  //---------- Get up plate rotated
  plate = getRotatedPlate( 0 );
  //---------- Reflect in plate
  lightray.reflect( plate );

  if (ALIUtils::debug >= 2) std::cout << std::endl << "$$$ LR: REFRACTION IN BACKWARD PLATE " << std::endl;
  //---------- Get backward plate
  plate = getPlate( 0, 1 );
  //---------- Refract in plate while exiting
  lightray.refract( plate, refra_ind2, refra_ind1 );

}

void OptOModifiedRhomboidPrism::detailedTraversesLightRay ( LightRay &  lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 87 of file OptOModifiedRhomboidPrism.cc.

References gather_cfg::cout, ALIUtils::debug, LightRay::dumpData(), mergeVDriftHistosByStation::name, and LightRay::refract().

{
  if (ALIUtils::debug >= 2) std::cout << "LR: DETAILED TRAVERSE MODIFIED RHOMBOID PRISM " << name() << std::endl;

  //---------- Get forward plate
  ALIPlane plate = getPlate(1, 1);
  //---------- Refract while entering splitter
  ALIdouble refra_ind1 = 1.;
  ALIdouble refra_ind2 = findExtraEntryValueMustExist("refra_ind");
  lightray.refract( plate, refra_ind1, refra_ind2 );
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Refracted in first plate"); 
  }

  //---------- Get back ward plate (of triangular piiece)
  plate = getPlate(1, 0);
  //---------- Refract while exiting prism
  lightray.refract( plate, refra_ind2, refra_ind1 );
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Refracted in first plate"); 
  }

}

void OptOModifiedRhomboidPrism::fastDeviatesLightRay ( LightRay &  lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 119 of file OptOModifiedRhomboidPrism.cc.

References gather_cfg::cout, ALIUtils::debug, LightRay::dumpData(), LightRay::intersect(), mergeVDriftHistosByStation::name, and LightRay::shiftAndDeviateWhileTraversing().

{
  if (ALIUtils::debug >= 2) std::cout << "LR: FAST REFLECTION IN MODIFIED RHOMBOID PRISM " << name() << std::endl;

  //---------- Get backward plate
  ALIPlane plate = getPlate(0, 0);
  //---------- Intersect with plate
  lightray.intersect( plate );
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Intersected in plate"); 
  }
  //---------- Deviate Lightray 
  lightray.shiftAndDeviateWhileTraversing( this, 'R');
  /*  ALIdouble deviRX = findExtraEntryValue("deviRX");
  ALIdouble deviRY = findExtraEntryValue("deviRY");
  ALIdouble shiftRX = findExtraEntryValue("shiftRX");
  ALIdouble shiftRY = findExtraEntryValue("shiftRY");
  lightray.shiftAndDeviateWhileTraversing( this, shiftRX, shiftRY, deviRX, deviRY);
  */

  if (ALIUtils::debug >= 2) {
    //    std::cout << " shiftRX " << shiftRX << " shiftRY " << shiftRY << std::endl;
    //   std::cout << " deviRX " << deviRX << " deviRY " << deviRY << std::endl;
    lightray.dumpData("Deviated ");
  }
  
}

void OptOModifiedRhomboidPrism::fastTraversesLightRay ( LightRay &  lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 157 of file OptOModifiedRhomboidPrism.cc.

References gather_cfg::cout, ALIUtils::debug, LightRay::dumpData(), LightRay::intersect(), mergeVDriftHistosByStation::name, and LightRay::shiftAndDeviateWhileTraversing().

{
  if (ALIUtils::debug >= 2) std::cout << "LR: FAST TRAVERSE MODIFIED RHOMBOID PRISM " << name() << std::endl;
  
  //---------- Get backward plate
  ALIPlane plate = getPlate(0, 0);
  lightray.intersect( plate );
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Intersected with plate"); 
  }
  //---------- Shift and Deviate
  lightray.shiftAndDeviateWhileTraversing( this, 'T');
  /*  ALIdouble shiftTX = findExtraEntryValue("shiftTX");
  ALIdouble shiftTY = findExtraEntryValue("shiftTY");
  ALIdouble deviTX = findExtraEntryValue("deviTX");
  ALIdouble deviTY = findExtraEntryValue("deviTY");
  lightray.shiftAndDeviateWhileTraversing( this, shiftTX, shiftTY, deviTX, deviTY);*/

  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Shifted and Deviated");
  }

}

ALIPlane OptOModifiedRhomboidPrism::getRotatedPlate

(

const ALIbool

forwardPlate

)

Definition at line 192 of file OptOModifiedRhomboidPrism.cc.

References gather_cfg::cout, ALIUtils::debug, ALIUtils::deg, ALIUtils::dump3v(), ALIUtils::dumprm(), mergeVDriftHistosByStation::name, and edm::shift.

{
  if (ALIUtils::debug >= 4) std::cout << "% LR: GET ROTATED PLATE " << name() << std::endl;
  //---------- Get OptO variables
  const ALIdouble shift = (findExtraEntryValue("shiftRY"));
  ALIdouble wedgeR = findExtraEntryValue("wedgeR");

  //---------- Get centre of plate
  //----- plate centre = OptO centre +/- 1/2 shift
  CLHEP::Hep3Vector plate_point = centreGlob();
  //--- Add to it half of the shift following the direction of the prism Y. -1/2 if it is forward plate, +1/2 if it is backward plate
  ALIdouble normal_sign = -forwardPlate*2 + 1;
  CLHEP::Hep3Vector YAxis(0.,1.,0.);
  CLHEP::HepRotation rmt = rmGlob();
  YAxis = rmt*YAxis;
  plate_point += normal_sign * shift/2. * YAxis;

  //---------- Get normal of plate
  //----- Plate normal before wedgeR (Z axis of OptO rotated 'angle' around X)
  CLHEP::Hep3Vector ZAxis(0.,0.,1.);
  ALIdouble anglePlanes;
  ALIbool we = findExtraEntryValueIfExists("anglePlanes", anglePlanes);
  if( !we ) { 
    anglePlanes = 45.*ALIUtils::deg;
  }
  ZAxis.rotateX( anglePlanes );

  //----- Rotate with global rotation matrix
  CLHEP::Hep3Vector plate_normal = rmt*ZAxis;
  if (ALIUtils::debug >= 3) { 
    ALIUtils::dump3v( plate_point, "plate_point");
    ALIUtils::dump3v( plate_normal, "plate_normal before wedge");
    ALIUtils::dumprm( rmt, "rmt before wedge angles" );
  }

  //----- Rotate plate normal by 1/2 wedgeR angles
  //--- Around X axis
  CLHEP::Hep3Vector XAxis(0.,0.,1.);
  XAxis = rmt*XAxis;
  plate_normal.rotate( normal_sign * wedgeR/2., XAxis );
  if (ALIUtils::debug >= 3) ALIUtils::dump3v( plate_normal, "plate_normal after wedgeR around X ");
  //--- Around the axis obtained rotating the prism Y axis by 'anglePlanes' around the prism X axis
  YAxis = CLHEP::Hep3Vector(0.,1.,0.);
  YAxis.rotateX( anglePlanes );
  YAxis = rmt*YAxis;
  plate_normal.rotate( normal_sign * wedgeR/2., YAxis );
  if (ALIUtils::debug >= 3) ALIUtils::dump3v( plate_normal, "plate_normal after wedgeR around Y ");

  //---------- Return plate plane
  return ALIPlane(plate_point, plate_normal);

}