#include <MucGeoStrip.h>

Public Member Functions
	MucGeoStrip (const int strip, MucGeoGap *const pGap)
	Default constructor.

	~MucGeoStrip ()
	Desctructor.

int	Part ()
	Get the softID.

int	Seg ()
	Get Seg No.

int	Gap ()
	Get Gap No.

int	Strip ()
	Get Strip No.

MucGeoGap *	GetGap () const
	Get the pointer to the gap that contains the strip.

MucGeoStrip *	GetLeftNeighbor () const
	Get the pointer to the adjacent strip on the -X or -Y side of this one.

MucGeoStrip *	GetRightNeighbor () const
	Get the pointer to the adjacent strip on the +X or +Y side of this one.

float	GetXmin () const
	Get position of low-X edge in the gap coordinate system.

float	GetXmax () const
	Get position of high-X edge in the gap coordinate system.

float	GetYmin () const
	Get position of low-Y edge in the gap coordinate system.

float	GetYmax () const
	Get position of high-Y edge in the gap coordinate system.

float	GetZmin () const
	Get position of low-Z edge in the gap coordinate system.

float	GetZmax () const
	Get position of high-Z edge in the gap coordinate system.

void	GetCenterPos (float &x, float &y, float &z) const
	Get center position of this strip (in the gap coordinate system).

void	GetCenterSigma (float &sx, float &sy, float &sz)
	Get uncertainty in the position of this strip (in the gap coordinate system).

bool	GetHitStatus () const
	Get the strip hit status.

void	SetStrip (const float x1, const float x2, const float y1, const float y2, const float xSize, const float ySize, const float zSize)
	Set the edge, center and sigma of the strip (in the gap coordinate system).

void	SetStrip (TGeoPhysicalNode *stripPhysicalNode)
	Set strip from ROOT.

void	SetLeftNeighbor (MucGeoStrip *p)
	Set pointer to the adjacent strip on the -X or -Y side of this one.

void	SetRightNeighbor (MucGeoStrip *p)
	Set pointer to the adjacent strip on the +X or +Y side of this one.

void	SetHitStatus (const bool hit)
	Set the strip hit status.

bool	IsInStrip (const float x, const float y, const float z) const
	Is the point within the boundaries of this strip?

bool	CrossGasChamber (const HepPoint3D linePoint, const Hep3Vector lineDir) const
	Does the line cross this strip?

Detailed Description

Class MucGeoStrip describes a single strip in a gap of the muon chamber

Author: Zhengyun You \URL{youzy.nosp@m.@pku.nosp@m..edu..nosp@m.cn}

Definition at line 40 of file MucGeoStrip.h.

Constructor & Destructor Documentation

◆ MucGeoStrip()

MucGeoStrip::MucGeoStrip	(	const int	strip,
		MucGeoGap *const	pGap
	)

Default constructor.

Definition at line 20 of file MucGeoStrip.cxx.

  : m_Strip(strip),
    m_pGap(pGap),
    m_pLeftNeighbor(0L),
    m_pRightNeighbor(0L),
    m_CenterPos(0.0,0.0,0.0),
    m_CenterSigma(0.0,0.0,0.0)
{
  // Constructor (dummy strip position and sizes).
}

◆ ~MucGeoStrip()

MucGeoStrip::~MucGeoStrip ( )

Desctructor.

Definition at line 31 of file MucGeoStrip.cxx.

{
  // Destructor.
  m_pLeftNeighbor = 0L;
  m_pRightNeighbor = 0L;
  m_pGap = 0L;
}

Member Function Documentation

◆ CrossGasChamber()

bool MucGeoStrip::CrossGasChamber	(	const HepPoint3D	linePoint,
		const Hep3Vector	lineDir
	)		const

Does the line cross this strip?

Definition at line 166 of file MucGeoStrip.cxx.

{
  // Does the line cross this strip?
 
  // We need to define the 6 planes then loop over to check for intersection.
  // All are defined in the gap coordinate.
  HepPlane3D stripSide[6];
 
  // Xmin.
  HepPoint3D p0Xmin(m_Xmin, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p1Xmin(m_Xmin, m_Ymax, m_dzFarFrontGas);
  HepPoint3D p2Xmin(m_Xmin, m_Ymin, m_dzFarBackGas);
  HepPlane3D plane0(p0Xmin, p1Xmin, p2Xmin);
  stripSide[0] = plane0;
 
  // Xmax.
  HepPoint3D p0Xmax(m_Xmax, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p1Xmax(m_Xmax, m_Ymax, m_dzFarFrontGas);
  HepPoint3D p2Xmax(m_Xmax, m_Ymin, m_dzFarBackGas);
  HepPlane3D plane1(p0Xmax, p1Xmax, p2Xmax);
  stripSide[1] = plane1;
 
  // Ymin.
  HepPoint3D p0Ymin(m_Xmin, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p1Ymin(m_Xmax, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p2Ymin(m_Xmin, m_Ymin, m_dzFarBackGas);
  HepPlane3D plane2(p0Ymin, p1Ymin, p2Ymin);
  stripSide[2] = plane2;
 
  // Ymax.
  HepPoint3D p0Ymax(m_Xmin, m_Ymax, m_dzFarFrontGas);
  HepPoint3D p1Ymax(m_Xmax, m_Ymax, m_dzFarFrontGas);
  HepPoint3D p2Ymax(m_Xmin, m_Ymax, m_dzFarBackGas);
  HepPlane3D plane3(p0Ymax, p1Ymax, p2Ymax);
  stripSide[3] = plane3;
 
  // Zmin.
  HepPoint3D p0Zmin(m_Xmin, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p1Zmin(m_Xmax, m_Ymin, m_dzFarFrontGas);
  HepPoint3D p2Zmin(m_Xmin, m_Ymax, m_dzFarFrontGas);
  HepPlane3D plane4(p0Zmin, p1Zmin, p2Zmin);
  stripSide[4] = plane4;
 
  // Zmax.
  HepPoint3D p0Zmax(m_Xmin, m_Ymin, m_dzFarBackGas);
  HepPoint3D p1Zmax(m_Xmax, m_Ymin, m_dzFarBackGas);
  HepPoint3D p2Zmax(m_Xmin, m_Ymax, m_dzFarBackGas);
  HepPlane3D plane5(p0Zmax, p1Zmax, p2Zmax);
  stripSide[5] = plane5;
 
  //cout << m_Xmin << "  " << m_Xmax << endl
  //   << m_Ymin << "  " << m_Ymax << endl
  //   << m_dzFarFrontGas << "  " << m_dzFarBackGas << endl;
    
  //If particle intersects this strip there should only be two crossings.
  int numIntersections = 0;
  MucGeometron geometron;
  bool intersectFlag = 0;
  HepPoint3D gCross(0,0,0);
 
  for (int sideID = 0; sideID<6; sideID++)
  {
    intersectFlag = geometron.GetIntersectionLinePlane(linePoint, 
                               lineDir,
                               stripSide[sideID],
                               gCross);
    if( intersectFlag && 
  ( gCross.x() >= m_Xmin - 1) &&  // expend strip size, to avoid gap between two strip
  ( gCross.x() <= m_Xmax + 1) &&
  ( gCross.y() >= m_Ymin - 1) &&
  ( gCross.y() <= m_Ymax + 1) &&
  ( ( gCross.z() >= m_dzFarFrontGas && gCross.z() <= m_dzNearFrontGas ) ||
      ( gCross.z() >= m_dzNearBackGas && gCross.z() <= m_dzFarBackGas) ) ) {
      numIntersections++;
      //cout << m_StripID.GetStrip() << gCross << endl;
    }
  }
 
  //if( numIntersections == 1 ) {
  //  cout<<"MucGeoStrip::IntersectStrip(): Warning only one intersection.  Should be two!" << endl;
  //}
  
  if( numIntersections > 2 ) {
//    cout<<"MucGeoStrip::CrossGasChamber(): Warning too many intersections " 
//  << numIntersections 
//  << endl;
  }
 
  if( numIntersections > 0 ) return true;
  
  return false;
}

Referenced by MucGeoGeneral::FindIntersectStrips().

◆ Gap()

int MucGeoStrip::Gap ( )

inline

Get Gap No.

Definition at line 59 of file MucGeoStrip.h.

59{ return m_pGap->Gap(); }

MucGeoGap::Gap

int Gap() const

Get gap identifier (0-8).

Definition: MucGeoGap.h:104

◆ GetCenterPos()

void MucGeoStrip::GetCenterPos	(	float &	x,
		float &	y,
		float &	z
	)		const

Get center position of this strip (in the gap coordinate system).

Definition at line 40 of file MucGeoStrip.cxx.

{
  // Get center position of this strip (in the gap coordinate system).
  x = m_CenterPos.x();
  y = m_CenterPos.y();
  z = m_CenterPos.z();
  
  return;
}

Referenced by MucGeoGeneral::FindIntersectStrips(), MucGeoGap::GuessStrip(), main(), and MucRecHit::MucRecHit().

◆ GetCenterSigma()

void MucGeoStrip::GetCenterSigma	(	float &	sx,
		float &	sy,
		float &	sz
	)

Get uncertainty in the position of this strip (in the gap coordinate system).

Definition at line 51 of file MucGeoStrip.cxx.

{
  // Get uncertainty in the position of this strip (in the gap coordinate system).
  sx = m_CenterSigma.x();
  sy = m_CenterSigma.y();
  sz = m_CenterSigma.z();
 
  return;
}

Referenced by ExtMucKal::ExtMucFilter(), MucRecHit::MucRecHit(), and ExtMucKal::TrackHit().

◆ GetGap()

MucGeoGap * MucGeoStrip::GetGap ( ) const

inline

Get the pointer to the gap that contains the strip.

Definition at line 65 of file MucGeoStrip.h.

65{ return m_pGap; }

Referenced by main(), and SetStrip().

◆ GetHitStatus()

bool MucGeoStrip::GetHitStatus ( ) const

inline

Get the strip hit status.

Definition at line 98 of file MucGeoStrip.h.

98{ return m_HitStatus; }

◆ GetLeftNeighbor()

MucGeoStrip * MucGeoStrip::GetLeftNeighbor ( ) const

inline

Get the pointer to the adjacent strip on the -X or -Y side of this one.

Definition at line 68 of file MucGeoStrip.h.

68{ return m_pLeftNeighbor; }

◆ GetRightNeighbor()

MucGeoStrip * MucGeoStrip::GetRightNeighbor ( ) const

inline

Get the pointer to the adjacent strip on the +X or +Y side of this one.

Definition at line 71 of file MucGeoStrip.h.

71{ return m_pRightNeighbor; }

◆ GetXmax()

float MucGeoStrip::GetXmax ( ) const

inline

Get position of high-X edge in the gap coordinate system.

Definition at line 77 of file MucGeoStrip.h.

77{ return m_Xmax; }

Referenced by MucGeoGeneral::FindIntersectStrips().

◆ GetXmin()

float MucGeoStrip::GetXmin ( ) const

inline

Get position of low-X edge in the gap coordinate system.

Definition at line 74 of file MucGeoStrip.h.

74{ return m_Xmin; }

Referenced by MucGeoGeneral::FindIntersectStrips().

◆ GetYmax()

float MucGeoStrip::GetYmax ( ) const

inline

Get position of high-Y edge in the gap coordinate system.

Definition at line 83 of file MucGeoStrip.h.

83{ return m_Ymax; }

Referenced by MucGeoGeneral::FindIntersectStrips().

◆ GetYmin()

float MucGeoStrip::GetYmin ( ) const

inline

Get position of low-Y edge in the gap coordinate system.

Definition at line 80 of file MucGeoStrip.h.

80{ return m_Ymin; }

Referenced by MucGeoGeneral::FindIntersectStrips().

◆ GetZmax()

float MucGeoStrip::GetZmax ( ) const

inline

Get position of high-Z edge in the gap coordinate system.

Definition at line 89 of file MucGeoStrip.h.

89{ return m_Zmax; }

◆ GetZmin()

float MucGeoStrip::GetZmin ( ) const

inline

Get position of low-Z edge in the gap coordinate system.

Definition at line 86 of file MucGeoStrip.h.

86{ return m_Zmin; }

◆ IsInStrip()

bool MucGeoStrip::IsInStrip	(	const float	x,
		const float	y,
		const float	z
	)		const

Is the point within the boundaries of this strip?

Definition at line 153 of file MucGeoStrip.cxx.

{
  // Is the point within the boundaries of this strip?
  return ( (x > m_Xmin) && (x < m_Xmax) 
       && (y > m_Ymin) && (y < m_Ymax)
       && (z >= m_dzFarFrontGas)
       && (z <= m_dzFarBackGas) );
}

◆ Part()

int MucGeoStrip::Part ( )

inline

Get the softID.

Get Part No.

Definition at line 53 of file MucGeoStrip.h.

53{ return m_pGap->Part(); }

MucGeoGap::Part

int Part() const

Get part identifier (0,2 for cap, 1 for barrel).

Definition: MucGeoGap.h:98

◆ Seg()

int MucGeoStrip::Seg ( )

inline

Get Seg No.

Definition at line 56 of file MucGeoStrip.h.

56{ return m_pGap->Seg(); }

MucGeoGap::Seg

int Seg() const

Get seg identifier (0-7).

Definition: MucGeoGap.h:101

◆ SetHitStatus()

void MucGeoStrip::SetHitStatus ( const bool hit )

inline

Set the strip hit status.

Definition at line 119 of file MucGeoStrip.h.

119{ m_HitStatus = hit; }

◆ SetLeftNeighbor()

void MucGeoStrip::SetLeftNeighbor ( MucGeoStrip * p )

Set pointer to the adjacent strip on the -X or -Y side of this one.

Definition at line 135 of file MucGeoStrip.cxx.

{
  /// Set pointer to the adjacent strip on the -X or -Y side of this one.
  // (0 if no neighbor in this gap).
  if (p)
    m_pLeftNeighbor = p;
}

Referenced by MucGeoGap::AddStrip().

◆ SetRightNeighbor()

void MucGeoStrip::SetRightNeighbor ( MucGeoStrip * p )

Set pointer to the adjacent strip on the +X or +Y side of this one.

Definition at line 144 of file MucGeoStrip.cxx.

{
  /// Set pointer to the adjacent strip on the +X or +Y side of this one.
  // (0 if no neighbor in this gap).
  if (p)
    m_pRightNeighbor = p;
}

Referenced by MucGeoGap::AddStrip().

◆ SetStrip() [1/2]

void MucGeoStrip::SetStrip	(	const float	x1,
		const float	x2,
		const float	y1,
		const float	y2,
		const float	xSize,
		const float	ySize,
		const float	zSize
	)

Set the edge, center and sigma of the strip (in the gap coordinate system).

Definition at line 62 of file MucGeoStrip.cxx.

{
  // Set the edge, center and sigma of the strip (in the gap coordinate system).
  double un = 1.0;  // Uncertainty.
  
  m_Xmin = x1;
  m_Xmax = x2;
  m_Ymin = y1;
  m_Ymax = y2;
  m_Zmin = -0.5*zSize;
  m_Zmax =  0.5*zSize;
 
  m_dzFarFrontGas  = GetGap()->GetdzFarFrontGas(); 
  m_dzNearFrontGas = GetGap()->GetdzNearFrontGas(); 
  m_dzNearBackGas  = GetGap()->GetdzNearBackGas(); 
  m_dzFarBackGas   = GetGap()->GetdzFarBackGas(); 
 
  m_CenterPos.setX(0.5*(m_Xmin + m_Xmax));
  m_CenterPos.setY(0.5*(m_Ymin + m_Ymax));
  m_CenterPos.setZ(0.0);
 
  m_CenterSigma.setX(un*xSize);
  m_CenterSigma.setY(un*ySize);
  m_CenterSigma.setZ(un*zSize);
}

Referenced by MucGeoGeneral::InitFromASCII(), and MucGeoGeneral::InitFromXML().

◆ SetStrip() [2/2]

void MucGeoStrip::SetStrip ( TGeoPhysicalNode * stripPhysicalNode )

Set strip from ROOT.

Definition at line 93 of file MucGeoStrip.cxx.

{
  // Set the edge, center and sigma of the strip (in the gap coordinate system).
  double un = 1.0;  // Uncertainty.
 
  TGeoBBox *stripBox = (TGeoBBox*)stripPhysicalNode->GetShape();
 
  double xSize = 2.0 * stripBox->GetDX();
  double ySize = 2.0 * stripBox->GetDY();
  double zSize = 2.0 * stripBox->GetDZ();
  
  m_dzFarFrontGas  = GetGap()->GetdzFarFrontGas(); 
  m_dzNearFrontGas = GetGap()->GetdzNearFrontGas(); 
  m_dzNearBackGas  = GetGap()->GetdzNearBackGas(); 
  m_dzFarBackGas   = GetGap()->GetdzFarBackGas(); 
 
  double eTrans[3], *pTrans;
  pTrans = eTrans;
  //int level = stripPhysicalNode->GetLevel();
  double eLocal[3], *pLocal;
  pLocal = eLocal;
  for (int i = 0; i < 3; i++) pLocal[i] = 0.0;
  stripPhysicalNode->GetNode()->LocalToMaster(pLocal, pTrans);
  
  m_CenterPos.setX( pTrans[0] );
  m_CenterPos.setY( pTrans[1] );
  m_CenterPos.setZ( pTrans[2] );
 
  m_Xmin = pTrans[0] - 0.5 * xSize;
  m_Xmax = pTrans[0] + 0.5 * xSize;
  m_Ymin = pTrans[1] - 0.5 * ySize;
  m_Ymax = pTrans[1] + 0.5 * ySize;
  m_Zmin = pTrans[2] - 0.5 * zSize;
  m_Zmax = pTrans[2] + 0.5 * zSize;
 
  m_CenterSigma.setX(un*xSize);
  m_CenterSigma.setY(un*ySize);
  m_CenterSigma.setZ(un*zSize);
  //cout << m_CenterPos << endl;
}

◆ Strip()

int MucGeoStrip::Strip ( )

inline

Get Strip No.

Definition at line 62 of file MucGeoStrip.h.

62{ return m_Strip; }

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ MucGeoStrip()

◆ ~MucGeoStrip()

Member Function Documentation

◆ CrossGasChamber()

◆ Gap()

◆ GetCenterPos()

◆ GetCenterSigma()

◆ GetGap()

◆ GetHitStatus()

◆ GetLeftNeighbor()

◆ GetRightNeighbor()

◆ GetXmax()

◆ GetXmin()

◆ GetYmax()

◆ GetYmin()

◆ GetZmax()

◆ GetZmin()

◆ IsInStrip()

◆ Part()

◆ Seg()

◆ SetHitStatus()

◆ SetLeftNeighbor()

◆ SetRightNeighbor()

◆ SetStrip() [1/2]

◆ SetStrip() [2/2]

◆ Strip()