Collaboration diagram for MuonChamberProjectionHelper::Imp:

Classes
class	MDTChamberInfo

Public Types
typedef std::map< GeoPVConstLink, MDTChamberInfo >::iterator	ChamberInfoMapItr

Public Member Functions
	Imp (MuonChamberProjectionHelper tc, StoreGateSvc ds)

bool	ensureInited ()

bool	init ()

bool	nameIsMDTChamber (const std::string &n)

bool	getMDTChamberVolInfo (const GeoPVConstLink &mdtChamber, ChamberInfoMapItr &itChamberInfo, bool silent=false)

void	getMDTChamberXAndZ (ChamberInfoMapItr &itChamberInfo, double &trdX, double &trdZ)

double	pointToPlaneDistAlongLine (const Amg::Vector3D &point, const Amg::Vector3D &lineDirection, const Amg::Vector3D &planePoint, const Amg::Vector3D &planeNormal)

void	projectXZPointToTrdAlongYAxis (const double &x, const double &z, const GeoTrd *trd, Amg::Vector3D &firstEndWall_point, Amg::Vector3D &secondEndWall_point)

bool	clip2DLineSegmentToRectangle (const double &rectX, const double &rectY, double &x0, double &y0, double &x1, double &y1)

const GeoTrd *	findTRDInShape (const GeoShape *shape)

Static Public Member Functions
static bool	constrainPointToRectangleAlongLine (const double &trdX, const double &trdZ, const double &x0, const double &z0, double &x1, double &z1)

Public Attributes
MuonChamberProjectionHelper *	theclass

StoreGateSvc *	detectorStore

bool	initattempted

bool	initsucceeded

std::map< GeoPVConstLink, MDTChamberInfo >::iterator	itLastMDTChamberLookedUp

std::map< GeoPVConstLink, MDTChamberInfo >	mdtchambervolinfo

Detailed Description

Definition at line 32 of file MuonChamberProjectionHelper.cxx.

Member Typedef Documentation

◆ ChamberInfoMapItr

typedef std::map<GeoPVConstLink,MDTChamberInfo>::iterator MuonChamberProjectionHelper::Imp::ChamberInfoMapItr

Definition at line 59 of file MuonChamberProjectionHelper.cxx.

Constructor & Destructor Documentation

◆ Imp()

MuonChamberProjectionHelper::Imp::Imp	(	MuonChamberProjectionHelper *	tc,
		StoreGateSvc *	ds
	)

inline

Definition at line 34 of file MuonChamberProjectionHelper.cxx.

                            : theclass(tc),
                  detectorStore(ds),
                  initattempted(false),
                  initsucceeded(false) {}

Member Function Documentation

◆ clip2DLineSegmentToRectangle()

bool MuonChamberProjectionHelper::Imp::clip2DLineSegmentToRectangle	(	const double &	rectX,
		const double &	rectY,
		double &	x0,
		double &	y0,
		double &	x1,
		double &	y1
	)

Definition at line 523 of file MuonChamberProjectionHelper.cxx.

 {
   if ( fabs(x0)<=rectX && fabs(y0)<=rectY  ) {
     if ( fabs(x1)>rectX || fabs(y1)>rectY  ) {
       //P0 inside R, P1 outside R. We must change (x1,y1) so P0-P1 is inside R.
       if (!constrainPointToRectangleAlongLine( rectX, rectY, x0, y0, x1, y1 ))
     theclass->message("MuonChamberProjectionHelper Error: Should never happen (1)");
     }
   } else {
     if ( fabs(x1)<=rectX && fabs(y1)<=rectY  ) {
       //Point P1 inside R, P0 outside R. We must change (x0,y0) so P0-P1 is inside R.
       if (!constrainPointToRectangleAlongLine( rectX, rectY, x1, y1, x0, y0 ))
     theclass->message("MuonChamberProjectionHelper Error: Should never happen (2)");
     } else {
       //Both points outside - but the line might still R!
  
       //First attempt to put (x1,y1) on edge of R, by sliding towards P0
       if (!constrainPointToRectangleAlongLine( rectX, rectY, x0, y0, x1, y1 )) {
     //failed - thus P0-P1 doesn't intersect target rectangle.
     return false;
       }
  
       //Now change (x0,y0) so P0-P1 is inside the target rectangle.
       if (!constrainPointToRectangleAlongLine( rectX, rectY, x1, y1, x0, y0 ))
     theclass->message("MuonChamberProjectionHelper Error: Should never happen (3)");
     }
   }
  
   return true;
  
 }

◆ constrainPointToRectangleAlongLine()

bool MuonChamberProjectionHelper::Imp::constrainPointToRectangleAlongLine	(	const double &	trdX,
		const double &	trdZ,
		const double &	x0,
		const double &	z0,
		double &	x1,
		double &	z1
	)

static

Definition at line 429 of file MuonChamberProjectionHelper.cxx.

 {
   //rectangle R : { (x,z) | abs(x)<=trdX && abs(z)<=trdZ }
   //
   //Returns false if segment from (x0,z0) to (x1,z1) doesn't cross R.
  
   if ( x1 < -trdX ) {
     if ( x0 < -trdX )
       return false;
     //Move (x1,z1) to x == -trdX line:
     z1 += (-trdX-x1)*(z1-z0)/(x1-x0);
     x1 = -trdX;
     if (fabs(z1)>trdZ)
       return constrainPointToRectangleAlongLine(trdX,trdZ,x0,z0,x1,z1);
   }
  
   if ( x1 > trdX ) {
     if ( x0 > trdX )
       return false;
     //Move (x1,z1) to x == trdX line:
     z1 += (trdX-x1)*(z1-z0)/(x1-x0);
     x1 = trdX;
     if (fabs(z1)>trdZ)
       return constrainPointToRectangleAlongLine(trdX,trdZ,x0,z0,x1,z1);
   }
  
   if ( z1 < -trdZ ) {
     if ( z0 < -trdZ )
       return false;
     //Move (x1,z1) to z == -trdZ line:
     x1 += (-trdZ-z1)*(x1-x0)/(z1-z0);
     z1 = -trdZ;
     if (fabs(x1)>trdX)
       return constrainPointToRectangleAlongLine(trdX,trdZ,x0,z0,x1,z1);
   }
  
   if ( z1 > trdZ ) {
     if ( z0 > trdZ )
       return false;
     //Move (x1,z1) to z == trdZ line:
     x1 += (trdZ-z1)*(x1-x0)/(z1-z0);
     z1 = trdZ;
     if (fabs(x1)>trdX)
       return constrainPointToRectangleAlongLine(trdX,trdZ,x0,z0,x1,z1);
   }
  
   //We were actually inside all along:
   return true;
 }

◆ ensureInited()

bool MuonChamberProjectionHelper::Imp::ensureInited ( )

Definition at line 158 of file MuonChamberProjectionHelper.cxx.

 {
   if (initattempted)
     return initsucceeded;
   initattempted = true;
   initsucceeded = init();
   return initsucceeded;
 }

◆ findTRDInShape()

const GeoTrd* MuonChamberProjectionHelper::Imp::findTRDInShape ( const GeoShape * shape )

inline

Definition at line 85 of file MuonChamberProjectionHelper.cxx.

   {
     if (shape->typeID()==GeoTrd::getClassTypeID())
       return static_cast<const GeoTrd*>(shape);
     if (shape->typeID() == GeoShapeShift::getClassTypeID() ) {
       const GeoShapeShift* theShift = static_cast<const GeoShapeShift*>(shape);
       return findTRDInShape(theShift->getOp());
     }
     if (shape->typeID() == GeoShapeSubtraction::getClassTypeID() ) {
       const GeoShapeSubtraction* theSubtraction = static_cast<const GeoShapeSubtraction*>(shape);
       const GeoTrd * trd = findTRDInShape(theSubtraction->getOpA());
       return trd ? trd : findTRDInShape(theSubtraction->getOpB());
     }
     if (shape->typeID() == GeoShapeUnion::getClassTypeID() ) {
       const GeoShapeUnion* theUnion = static_cast<const GeoShapeUnion*>(shape);
       const GeoTrd * trd = findTRDInShape(theUnion->getOpA());
       return trd ? trd : findTRDInShape(theUnion->getOpB());
     }
     if (shape->typeID() == GeoShapeIntersection::getClassTypeID() ) {
       const GeoShapeIntersection* theIntersection = static_cast<const GeoShapeIntersection*>(shape);
       const GeoTrd * trd = findTRDInShape(theIntersection->getOpA());
       return trd ? trd : findTRDInShape(theIntersection->getOpB());
     }
     return 0;
   }

◆ getMDTChamberVolInfo()

bool MuonChamberProjectionHelper::Imp::getMDTChamberVolInfo	(	const GeoPVConstLink &	mdtChamber,
		ChamberInfoMapItr &	itChamberInfo,
		bool	silent = `false`
	)

Definition at line 333 of file MuonChamberProjectionHelper.cxx.

                                                {
  
   if (!ensureInited())
     return false;
  
   if (itLastMDTChamberLookedUp->first == mdtChamber) {
     itChamberInfo = itLastMDTChamberLookedUp;
     return true;
   }
  
   itChamberInfo = mdtchambervolinfo.find(mdtChamber);
   if (itChamberInfo == mdtchambervolinfo.end()) {
     if (!silent)
       theclass->message("MuonChamberProjectionHelper Error: Can't find MDT chamber among the "
             +QString::number(mdtchambervolinfo.size())+" registered");
     return false;
   }
  
   itLastMDTChamberLookedUp = itChamberInfo;
  
   return true;
 }

◆ getMDTChamberXAndZ()

void MuonChamberProjectionHelper::Imp::getMDTChamberXAndZ	(	ChamberInfoMapItr &	itChamberInfo,
		double &	trdX,
		double &	trdZ
	)

Definition at line 359 of file MuonChamberProjectionHelper.cxx.

 {
   const GeoTrd * trd = itChamberInfo->second.trd;
   trdX = trd->getXHalfLength1();
   trdZ = trd->getZHalfLength();
   if ( trdX != trd->getXHalfLength2() ) {
     theclass->message("MuonChamberProjectionHelper Warning: x1!=x2 in GeoTrd shape. Clippings etc. will be to a too large surface.");
     if ( trdX < trd->getXHalfLength2() )
       trdX = trd->getXHalfLength2();
   }
  
 }

◆ init()

bool MuonChamberProjectionHelper::Imp::init ( )

Definition at line 185 of file MuonChamberProjectionHelper.cxx.

 {
   if (!VP1JobConfigInfo::hasMuonGeometry()) {
     theclass->messageDebug("Warning: Can't init since muon geometry information is not present." );
     return false;
   }
  
   VP1SGAccessHelper * sgaccess(0);
   if (theclass->systemBase())
     sgaccess = new VP1SGAccessHelper(theclass->systemBase());
   else
     sgaccess = new VP1SGAccessHelper(detectorStore);
  
  //Locate the world volume if possible:
   const GeoModelExperiment * theExpt;
   if (!sgaccess->retrieve( theExpt, "ATLAS" )) {
     theclass->message("MuonChamberProjectionHelper Error: Can't retrieve"
               " the ATLAS GeoModelExperiment from detector store.");
     delete sgaccess;
     return false;
   }
   delete sgaccess;
  
   PVConstLink world(theExpt->getPhysVol());
   GeoVolumeCursor av(world);
   const GeoLogVol * logvol(0);
   const GeoShape * shape(0);
   while (!av.atEnd()) {
  
     if (av.getName()!="Muon") {
       av.next();
       continue;
     }
     Amg::Transform3D toptransform = av.getTransform();
     GeoVolumeCursor av2(av.getVolume());
     while (!av2.atEnd()) {
       if (nameIsMDTChamber(av2.getName())) {
     logvol = av2.getVolume()->getLogVol();
     if (!logvol) {
       theclass->message("MuonChamberProjectionHelper Error: Chamber has null logvol");
       av2.next(); // increment volume cursor.
       continue;
     }
     shape = logvol->getShape();
     if (!shape) {
       theclass->message("MuonChamberProjectionHelper Error: Chamber has null shape");
       av2.next(); // increment volume cursor.
       continue;
     }
     const GeoTrd * trd = findTRDInShape(shape);
     if (trd) {
       //Fixme: Test for nan's and clearly crazy numbers (if more than 1km for instance).
       if ( trd->getZHalfLength()>0.0
            && trd->getXHalfLength1() > 0.0
            && trd->getXHalfLength2() > 0.0
            && trd->getYHalfLength1() > 0.0
            && trd->getYHalfLength2() > 0.0 ) {
         trd->ref();
         Amg::Transform3D geovolume_transf = av2.getTransform();
         mdtchambervolinfo.insert(std::pair<const GeoPVConstLink,MDTChamberInfo>( av2.getVolume(), Imp::MDTChamberInfo(toptransform * geovolume_transf, trd) ));
       } else {
         theclass->message("MuonChamberProjectionHelper Error: Chamber trd has non-positive shape parameters!");
       }
     } else {
       theclass->message("MuonChamberProjectionHelper Error: Chamber shape is not a GeoTrd, and is not a boolean with a Trd somewhere");
     }
       }
       av2.next(); // increment volume cursor.
     }
     av.next(); // increment volume cursor.
   }
  
   if (mdtchambervolinfo.empty()) {
     theclass->message("MuonChamberProjectionHelper Error: Found no MDT chambers");
     return false;
   }
  
   itLastMDTChamberLookedUp = mdtchambervolinfo.begin();
  
   return true;
 }

◆ nameIsMDTChamber()

bool MuonChamberProjectionHelper::Imp::nameIsMDTChamber ( const std::string & n )

inline

Definition at line 168 of file MuonChamberProjectionHelper.cxx.

 {
   //MDT chamber volumes starts with:
   //    Endcap: EI, EM, EO or EE
   //    Inner barrel: BI or BEE
   //    Middle barrel: BM
   //    Outer barrel:BO
   if (n.size()<3)
     return false;
   if (n[0]=='E')
     return n[1]=='I' || n[1]=='M' || n[1]=='O' || n[1]=='E';
   if (n[0]=='B')
     return n[1]=='I' || n[1]=='M' || n[1]=='O' || (n[1]=='E'&&n[2]=='E');
   return false;
 }

◆ pointToPlaneDistAlongLine()

double MuonChamberProjectionHelper::Imp::pointToPlaneDistAlongLine	(	const Amg::Vector3D &	point,
		const Amg::Vector3D &	lineDirection,
		const Amg::Vector3D &	planePoint,
		const Amg::Vector3D &	planeNormal
	)

inline

Definition at line 317 of file MuonChamberProjectionHelper.cxx.

 {
  
   double denominator(planeNormal.dot(lineDirection)*lineDirection.mag());
   if (denominator==0.0) {
     theclass->message("MuonChamberProjectionHelper Error: pointToPlaneDistAlongLine is undefined!");
     return -1.0;
   }
   double numerator(planeNormal.x() * (planePoint.x() - point.x())
            + planeNormal.y() * (planePoint.y() - point.y())
            + planeNormal.z() * (planePoint.z() - point.z()));
   return fabs(numerator/denominator);
 }

◆ projectXZPointToTrdAlongYAxis()

void MuonChamberProjectionHelper::Imp::projectXZPointToTrdAlongYAxis	(	const double &	x,
		const double &	z,
		const GeoTrd *	trd,
		Amg::Vector3D &	firstEndWall_point,
		Amg::Vector3D &	secondEndWall_point
	)

Definition at line 418 of file MuonChamberProjectionHelper.cxx.

 {
   const double epsilon(0.1);//100micron
   const double trdY1(trd->getYHalfLength1()), trdY2(trd->getYHalfLength2());
   const double y( trdY1 + 0.5*(1.0+z/trd->getZHalfLength())*(trdY2-trdY1) );
   Amg::setVector3DCartesian(firstEndWall_point, x,y+epsilon,z);
   Amg::setVector3DCartesian(secondEndWall_point, x,-y-epsilon,z);
 }