MuonChamberProjectionHelper::Imp Class Reference

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 58 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 502 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 410 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 143 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 84 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 318 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 344 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 170 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 153 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 302 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 399 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) );
  firstEndWall_point= Amg::Vector3D{ x,y+epsilon,z};
  secondEndWall_point= Amg::Vector3D{x,-y-epsilon,z};
}

Member Data Documentation

detectorStore

StoreGateSvc* MuonChamberProjectionHelper::Imp::detectorStore {}

Definition at line 40 of file MuonChamberProjectionHelper.cxx.

{};

initattempted

bool MuonChamberProjectionHelper::Imp::initattempted {}

Definition at line 41 of file MuonChamberProjectionHelper.cxx.

{};

initsucceeded

bool MuonChamberProjectionHelper::Imp::initsucceeded {}

Definition at line 42 of file MuonChamberProjectionHelper.cxx.

{};

itLastMDTChamberLookedUp

std::map<GeoPVConstLink,MDTChamberInfo>::iterator MuonChamberProjectionHelper::Imp::itLastMDTChamberLookedUp {}

Definition at line 56 of file MuonChamberProjectionHelper.cxx.

{};

mdtchambervolinfo

std::map<GeoPVConstLink,MDTChamberInfo> MuonChamberProjectionHelper::Imp::mdtchambervolinfo

Definition at line 57 of file MuonChamberProjectionHelper.cxx.

theclass

MuonChamberProjectionHelper* MuonChamberProjectionHelper::Imp::theclass {}

Definition at line 39 of file MuonChamberProjectionHelper.cxx.

{};

---

The documentation for this class was generated from the following file:

• MuonChamberProjectionHelper.cxx