#include <OptOMirror.h>

Inheritance diagram for OptOMirror:

Public Member Functions
void	constructSolidShape ()

void	defaultBehaviour (LightRay &lightray, Measurement &meas)

virtual void	detailedDeviatesLightRay (LightRay &lightray)

virtual void	detailedTraversesLightRay (LightRay &lightray)

virtual void	fastDeviatesLightRay (LightRay &lightray)

virtual void	fastTraversesLightRay (LightRay &lightray)

	OptOMirror ()

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

	~OptOMirror ()

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

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

Constructor & Destructor Documentation

OptOMirror::OptOMirror ( )

inline

Definition at line 23 of file OptOMirror.h.

23 { };

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

inline

Definition at line 24 of file OptOMirror.h.

24 :

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

type

Definition: HCALResponse.h:21

OpticalObject::OpticalObject

OpticalObject()

Definition: OpticalObject.h:41

OpticalObject::name

const ALIstring & name() const

Definition: OpticalObject.h:60

OptOMirror::~OptOMirror ( )

inline

Definition at line 26 of file OptOMirror.h.

26 { };

Member Function Documentation

void OptOMirror::constructSolidShape ( )

virtual

Reimplemented from OpticalObject.

Definition at line 213 of file OptOMirror.cc.

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

 {
   ALIdouble go;
   GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
   gomgr->getGlobalOptionValue("VisScale", go );
 
   theSolidShape = new CocoaSolidShapeBox( "Box", go*5.*cm/m, go*5.*cm/m, go*1.*cm/m ); //COCOA internal units are meters
 }

void OptOMirror::defaultBehaviour	(	LightRay &	lightray,
		Measurement &	meas
	)

virtual

Reimplemented from OpticalObject.

Definition at line 24 of file OptOMirror.cc.

 {
   detailedDeviatesLightRay( lightray );
 }

void OptOMirror::detailedDeviatesLightRay ( LightRay & lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 33 of file OptOMirror.cc.

References ALI_DBL_MAX, gather_cfg::cout, ALIUtils::debug, LightRay::direction(), HLT_25ns14e33_v1_cff::distance, ALIUtils::dump3v(), LightRay::dumpData(), LightRay::intersect(), mag(), mergeVDriftHistosByStation::name, LightRay::point(), LightRay::setDirection(), and create_public_lumi_plots::width.

 {
   if (ALIUtils::debug >= 2) std::cout << "LR: DETAILED REFLECTION IN MIRROR " << name() << std::endl;
   if (ALIUtils::debug >= 3) ALIUtils::dump3v( centreGlob(), " centre Global ");
 
   //---------- Get forward plate and intersect lightray with it
   ALIPlane plate = getPlate(1, 0);
   lightray.intersect( plate );
   CLHEP::Hep3Vector inters = lightray.point( );
 
   //---------- Get centre of forward plate
   //----- Get Z axis
   CLHEP::Hep3Vector ZAxis = getZAxis();
   CLHEP::HepRotation rmt = rmGlob();
   ZAxis = rmt * ZAxis;
   //----- Get centre
   ALIdouble width = findExtraEntryValue("width");
   CLHEP::Hep3Vector plate_centre = centreGlob() - 0.5 * width * ZAxis; 
   //-  if(ALIUtils::debug >= 4) std::cout << " mirror width " << width << std::endl;
 
   //---------- Get the distance between the intersection point and the centre of the forward plate
   ALIdouble distance = ( plate_centre - inters ).mag();
 
   //---------- Get normal to mirror at intersection point 
   //------- Get angle of mirror surface (angle between plate centre and intersection)
   ALIdouble flatness = findExtraEntryValue("flatness");
   //-- flatness is defined as number of 632 nm wavelengths
   flatness *= 632.E-9;
   ALIdouble length = 0.;
 
   ALIdouble curvature_radius;
   ALIdouble angFlatness;
   if( flatness != 0) {
     length = findExtraEntryValueMustExist("length");
     curvature_radius = (flatness*flatness + length*length) / (2*flatness);
     angFlatness = asin( distance / curvature_radius);
   } else {
     curvature_radius = ALI_DBL_MAX;
     angFlatness = 0;
   }
 
   if (ALIUtils::debug >= 3) {
     std::cout << " intersection with plate " << inters << std::endl;
     std::cout << " plate_centre " << plate_centre << std::endl;
     std::cout << " distance plate_centre - intersection " << distance << std::endl;
     std::cout << " flatness " << flatness << ", length " << length;
     std::cout << ", curvature radius " << curvature_radius << " angle of flatness " << angFlatness << std::endl;
   }
 
   //----- Axis of rotation is perpendicular to Z Axis and to line plate_centre - intersection
   CLHEP::Hep3Vector ipcV = inters - plate_centre;
   if( ipcV.mag() != 0) ipcV *= 1./ipcV.mag();
   CLHEP::HepRotation rtm = rmGlob();
   ipcV = rtm*ipcV;
   CLHEP::Hep3Vector rotationAxis = ipcV.cross(ZAxis);
   //----- normal is object normal rotated around this axis
   CLHEP::Hep3Vector inters_normal = CLHEP::Hep3Vector(0.,0.,1.);
   inters_normal.rotate( angFlatness, rotationAxis );  
   inters_normal = rmt * inters_normal;
 
   if (ALIUtils::debug >= 2) {
      ALIUtils::dump3v( ipcV, " intersection -  plate_centre std::vector ");
      std::cout << "rotation Axis " << rotationAxis << std::endl;
      std::cout << " plate normal at intersection point " << inters_normal << std::endl;
 }
   //---------- Reflect in plate 
   ALIdouble cosang = -( inters_normal * lightray.direction() ) / 
            inters_normal.mag() / lightray.direction().mag();
   CLHEP::Hep3Vector lrold = lightray.direction();
   lightray.setDirection( lightray.direction() + inters_normal*2*cosang );
 
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Reflected in mirror"); 
   }
 
 }

void OptOMirror::detailedTraversesLightRay ( LightRay & lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 138 of file OptOMirror.cc.

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

 {
   if (ALIUtils::debug >= 2) std::cout << "LR: DETAILED TRAVERSE IN MIRROR " << name() << std::endl;
 
   //---------- Get forward plate
   ALIPlane plate = getPlate(1, 1);
   //---------- If width is 0, just keep the same point 
   ALIdouble width = findExtraEntryValue("width");
   if( width == 0 ) {
     if(ALIUtils::debug >= 3) lightray.dumpData("Traversed with 0 width"); 
     return;
   }
 
   //---------- Refract while entering mirror
   ALIdouble refra_ind1 = 1.;
   ALIdouble refra_ind2 = findExtraEntryValue("refra_ind");
   lightray.refract( plate, refra_ind1, refra_ind2 );
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Refracted in first plate"); 
   }
 
   //---------- Get backward plate
   plate = getPlate(0, 1);
   //---------- Refract while exiting mirror
   lightray.refract( plate, refra_ind2, refra_ind1 );
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Refracted in first plate"); 
   }
 
 }

void OptOMirror::fastDeviatesLightRay ( LightRay & lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 114 of file OptOMirror.cc.

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

 {
   if (ALIUtils::debug >= 2) std::cout << "LR: FAST REFLECTION IN MIRROR " << name() << std::endl;
 
   //---------- Get forward plate
   ALIPlane plate = getPlate(1, 0);
 
   //---------- Reflect in plate (including intersection with it)
   lightray.reflect( plate );
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Reflected in plate"); 
   }
   //---------- Deviate Lightray 
   //  ALIdouble deviX = findExtraEntryValue("deviX");
   // ALIdouble deviY = findExtraEntryValue("deviY");
   //  lightray.shiftAndDeviateWhileTraversing( this, 0., 0., 0., deviX, deviY, 0.);
   lightray.shiftAndDeviateWhileTraversing( this, 'R' );
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Deviated ");
   }
 
 }

void OptOMirror::fastTraversesLightRay ( LightRay & lightray )

virtual

Reimplemented from OpticalObject.

Definition at line 171 of file OptOMirror.cc.

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

 {
   if (ALIUtils::debug >= 2) std::cout << "LR: TRAVERSE MIRROR  " << 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' );
   if (ALIUtils::debug >= 2) {
     lightray.dumpData("Shifted and Deviated");
   }
 
 }

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation