RecMucTrack Class Reference

#include <RecMucTrack.h>

Inheritance diagram for RecMucTrack:

List of all members.

Public Member Functions
virtual const CLID &	clID () const
	RecMucTrack ()
	Constructor.
RecMucTrack &	operator= (const RecMucTrack &orig)
	Assignment constructor.
	RecMucTrack (const RecMucTrack &source)
	Copy constructor.
RecMucTrack &	operator= (const DstMucTrack &dstTrack)
	Assignment constructor from DstMucTrack.
	RecMucTrack (const DstMucTrack &dstTrack)
	Copy constructor from DstMucTrack.
	~RecMucTrack ()
	Destructor.
void	SetDefault ()
void	setTrackId (const int trackId)
	set the index for this track.
void	SetMdcPos (const float x, const float y, const float z)
	set start position of the track in Mdc.
void	SetMdcMomentum (const float px, const float py, const float pz)
	set start moment of the track in Mdc.
void	SetExtMucPos (const float x, const float y, const float z)
	set start position of ext track in Muc. (compute from MdcPos MdcMomentum or get from ExtTrack)
void	SetExtMucMomentum (const float px, const float py, const float pz)
	set start moment of ext track in Muc.
void	SetMucPos (const float x, const float y, const float z)
	set start position of the track in Muc. (after line fit and correction)
void	SetMucPosSigma (const float sigmax, const float sigmay, const float sigmaz)
void	SetMucMomentum (const float px, const float py, const float pz)
	set start moment of the track in Muc.
void	SetCurrentPos (const float x, const float y, const float z)
	set current position of the trajectory.
void	SetCurrentDir (const float x, const float y, const float z)
	set current direction of the trajectory.
void	SetCurrentInsct (const float x, const float y, const float z)
	set current intersection of the trajectory with strip plane.
void	SetExtTrack (RecExtTrack *extTrack)
	set Ext track point.
void	SetExtTrackID (int id)
	set Ext track id. for compute from mdc myself
void	GetMdcExtTrack (Hep3Vector mdcStartPos, Hep3Vector mdcStartMomentum, int charge, Hep3Vector &mucStartPos, Hep3Vector &mucStartMomentum)
	compute ext track myself from mdc.
bool	IsInsideSuperConductor (Hep3Vector pos)
void	AttachHit (MucRecHit *hit)
	Attach the given hit to this track.
void	Project (const int &part, const int &gap, float &x, float &y, float &z, int &seg)
	Where does the trajectory of this track intersect a specific gap?
float	GetHitDistance (const MucRecHit *hit)
	Calculate the distance of the hit to the intersection in read direction.
float	GetHitDistance2 (const MucRecHit *hit)
	no abs value
Hep3Vector	CalculateInsct (const int part, const int seg, const int gap)
	Calculate intersection from all hits attached on this gap.
void	AttachInsct (Hep3Vector insct)
	Attach the intersection to this trajectory.
void	AttachDirection (Hep3Vector dir)
	Attach the direction to this trajectory.
void	CorrectDir ()
	Correct direction of this trajectory.
void	CorrectPos ()
	Correct current position of this trajectory.
int	getExtTrackID () const
	The identifier of the Ext track as seed.
Hep3Vector	getExtMucPos () const
	start position of the Ext track in Muc.
Hep3Vector	getExtMucMomentum () const
	Start momentum of the Ext track in Muc.
Hep3Vector	getMucPos () const
	start position of this track in Muc.
Hep3Vector	getMucPosSigma () const
Hep3Vector	getMucMomentum () const
	Start momentum of this track in Muc.
Hep3Vector	getMdcMomentum () const
	momentum of this track in Mdc
Hep3Vector	GetCurrentPos () const
	Current position.
Hep3Vector	GetCurrentDir () const
	Current direction.
Hep3Vector	GetCurrentInsct () const
	Current intersection.
void	ComputeLastGap ()
	Comute last gap in barrel and endcap.
int	brFirstLayer () const
	Which gap on Barrel is the first one with hits attached to this track?
int	ecFirstLayer () const
	Which gap on EndCap is the first one with hits attached to this track?
void	LineFit (int fittingMethod)
	Line fit with hits on a seg with max hits.
void	ComputeDepth ()
	Compute depth.
int	ComputeDepth (int method)
	Compute depth.
void	ComputeDistanceMatch ()
	Compute distance match //2006.11.08.
void	Extend ()
	Extend mucpos and extmucpos to first layer of muc.
float	GetDepth3 () const
	Length of the track penetrating in iron absorber.
void	SetRecMode (int recmode)
int	GetRecMode () const
Hep3Vector	GetExtMucDistance () const
	Distance match of the ext track with muc track in first layer.
Hep3Vector	GetMucStripPos () const
void	ComputeNGapsWithHits ()
	ComputeNGapsWithHits;.
int	GetNGapsWithHits () const
	How many gaps provide hits attached to this track?
int	GetTotalHits () const
	How many hits in all does this track contain?
int	GetHitInGap (const int part, const int gap) const
	How many hits per gap does this track contain?
int	GetHitInSeg (const int part, const int seg) const
	How many hits does a segment contains.
int	GetHitInSegOrient (const int part, const int seg, const int orient) const
	How many hits does a segment contains in one orient.
int	FindSegWithMaxHits (int &part, int &seg)
	Find the segment which contains most hits, return max hits number.
void	ComputeMaxHitsInGap ()
	ComputeMaxHitsInGap;.
bool	HasHit (const int part, const int seg, const int gap, const int strip) const
	Does this track contains an assigned hit?
bool	HasHitInGap (const int part, const int gap) const
	Does this track contain any hits in the given gap?
int	GetNSharedHits (const RecMucTrack *track) const
	How many hits do two tracks share?
MucRecHit *	GetHit (const int part, const int gap) const
	Get a pointer to the first hit attached in a particular gap.
vector< MucRecHit * >	GetHits () const
	Get all hits on this track.
vector< int >	getVecHits () const
vector< MucRecHit * >	GetExpectedHits () const
vector< int >	getExpHits () const
vector< float >	getDistHits () const
vector< float >	getQuadDistHits () const
vector< float >	getExtDistHits () const
vector< long >	GetHitIndices () const
	Get indices of all hits in the track.
void	ComputeTrackInfo (int fittingmethod)
	Compute all infomation of this track;.
void	PrintHitsInfo () const
	Print Hits Infomation.
void	OutputUnitChange ()
	change unit
void	setVecHits (vector< MucRecHit * > &pHits)
	reload setVecHits
void	setExpHits (vector< MucRecHit * > &pHits)
void	setVecHits (vector< int > &vechits)
void	setExpHits (vector< int > &exphits)
void	setDistHits (vector< float > &disthits)
void	setQuadDistHits (vector< float > &disthits)
void	setExtDistHits (vector< float > &disthits)
void	pushExtDistHits (float dist)
int	trackId () const
int	id () const
int	status () const
int	type () const
int	startPart () const
int	endPart () const
int	brLastLayer () const
int	ecLastLayer () const
int	numHits () const
int	numLayers () const
int	maxHitsInLayer () const
double	depth () const
double	chi2 () const
int	dof () const
double	rms () const
double	xPos () const
double	yPos () const
double	zPos () const
double	xPosSigma () const
double	yPosSigma () const
double	zPosSigma () const
double	px () const
double	py () const
double	pz () const
double	distance () const
double	deltaPhi () const
vector< int >	vecHits () const
double	kalRechi2 () const
int	kaldof () const
double	kaldepth () const
int	kalbrLastLayer () const
int	kalecLastLayer () const
void	setId (int id)
void	setStatus (int st)
void	setType (int type)
void	setStartPart (int part)
void	setEndPart (int part)
void	setBrLastLayer (int layer)
void	setEcLastLayer (int layer)
void	setNumHits (int numHits)
void	setNumLayers (int numLayers)
void	setMaxHitsInLayer (int maxHits)
void	setDepth (double depth)
void	setChi2 (double chi2)
void	setDof (int dof)
void	setRms (double rms)
void	setXPos (double x)
void	setYPos (double y)
void	setZPos (double z)
void	setXPosSigma (double xsigma)
void	setYPosSigma (double ysigma)
void	setZPosSigma (double zsigma)
void	setPx (double px)
void	setPy (double py)
void	setPz (double pz)
void	setDistance (double dist)
void	setDeltaPhi (double dphi)
void	setkalRechi2 (double ch)
void	setkalDof (int f)
void	setkalDepth (double de)
void	setkalbrLastLayer (int br)
void	setkalecLastLayer (int ec)
Static Public Member Functions
static const CLID &	classID ()
Protected Attributes
int	m_trackId
int	m_id
int	m_status
int	m_type
int	m_startPart
int	m_endPart
int	m_brLastLayer
int	m_ecLastLayer
int	m_numHits
int	m_numLayers
int	m_maxHitsInLayer
double	m_depth
double	m_chi2
int	m_dof
double	m_rms
double	m_xPos
double	m_yPos
double	m_zPos
double	m_xPosSigma
double	m_yPosSigma
double	m_zPosSigma
double	m_px
double	m_py
double	m_pz
double	m_distance
double	m_deltaPhi
double	m_kalrechi2
int	m_kaldof
double	m_kaldepth
int	m_kalbrLastLayer
int	m_kalecLastLayer
Private Attributes
int	m_ExtTrackID
RecExtTrack *	m_ExtTrack
Hep3Vector	m_MdcPos
Hep3Vector	m_MdcMomentum
Hep3Vector	m_ExtMucPos
Hep3Vector	m_ExtMucMomentum
Hep3Vector	m_MucPos
Hep3Vector	m_MucMomentum
Hep3Vector	m_MucPosSigma
Hep3Vector	m_CurrentPos
Hep3Vector	m_CurrentDir
Hep3Vector	m_CurrentInsct
Hep3Vector	m_Dist_muc_ext
Hep3Vector	m_StripPhi
int	m_brFirstLayer
int	m_ecFirstLayer
int	m_ecPart
float	m_depth_3
int	m_recmode
int	m_Good3DLine
int	m_Good3DPart
bool	m_changeUnit
Hep3Vector	m_IntersectionInner [9]
Hep3Vector	m_IntersectionOuter [9]
vector< int >	m_vecHits
vector< int >	m_expHits
vector< float >	m_distHits
vector< float >	m_distHits_quad
vector< float >	m_distHits_ext
vector< MucRecHit * >	m_pHits
vector< MucRecHit * >	m_pExpectedHits
vector< Hep3Vector >	m_Intersections
vector< Hep3Vector >	m_Directions

---

Detailed Description

Definition at line 40 of file RecMucTrack.h.

---

Constructor & Destructor Documentation

RecMucTrack::RecMucTrack

(

)

Constructor.

Definition at line 26 of file RecMucTrack.cxx.

References m_brFirstLayer, DstMucTrack::m_brLastLayer, m_changeUnit, DstMucTrack::m_chi2, DstMucTrack::m_deltaPhi, DstMucTrack::m_depth, DstMucTrack::m_dof, m_ecFirstLayer, DstMucTrack::m_ecLastLayer, m_ecPart, DstMucTrack::m_endPart, DstMucTrack::m_id, m_IntersectionInner, m_IntersectionOuter, DstMucTrack::m_kalbrLastLayer, DstMucTrack::m_kaldepth, DstMucTrack::m_kaldof, DstMucTrack::m_kalecLastLayer, DstMucTrack::m_kalrechi2, DstMucTrack::m_maxHitsInLayer, DstMucTrack::m_numHits, DstMucTrack::m_numLayers, m_recmode, DstMucTrack::m_rms, DstMucTrack::m_startPart, DstMucTrack::m_status, DstMucTrack::m_type, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPosSigma, DstMucTrack::m_zPosSigma, and eformat::write::set().

  : //m_trackId(-1),    //--------------->
  m_ExtTrackID(-1),
  m_MdcPos(0.0, 0.0, 0.0),
  m_MdcMomentum(0.0, 0.0, 0.0),
  m_MucPos(0.0, 0.0, 0.0),
  m_MucPosSigma(0.0, 0.0, 0.0),
  m_MucMomentum(0.0, 0.0, 0.0),
  m_CurrentPos(0.0, 0.0, 0.0),
  m_CurrentDir(0.0, 0.0, 0.0),
  m_CurrentInsct(0.0, 0.0, 0.0),
  m_Good3DLine(0),
  m_pHits(0),
  m_pExpectedHits(0),
  m_Intersections(0),
  m_Directions(0)
{ 
  // initialize m_IntersectionInner/Outer.
  for(int igap = 0; igap < 9; igap++){
    m_IntersectionInner[igap].set(-9999,-9999,-9999);
    m_IntersectionOuter[igap].set(-9999,-9999,-9999);    
  }
  m_id = 0;
  m_status = -1;
  m_type = -1;

  m_numHits = 0;  
  m_startPart = -1;
  m_endPart = -1;
  m_brLastLayer = -1;
  m_ecLastLayer = -1;
  m_brFirstLayer = -1;
  m_ecFirstLayer = -1;
  m_ecPart = -1;
  m_numLayers = 0;
  m_maxHitsInLayer = 0;
  m_depth = -99;
  m_dof = 0;
  m_chi2 = 0.0;
  m_rms = 0.0;
  m_deltaPhi = 0.0; 

  m_xPosSigma = 0.0;
  m_yPosSigma = 0.0;
  m_zPosSigma = 0.0;

  m_changeUnit = false;
  m_recmode = 0;
  //added by LI Chunhua
  m_kalrechi2 =0. ;
  m_kaldof =0;
  m_kaldepth = -99;
  m_kalbrLastLayer = -1;
  m_kalecLastLayer = -1;
}

RecMucTrack::RecMucTrack

(

const RecMucTrack &

source

)

Copy constructor.

Definition at line 134 of file RecMucTrack.cxx.

References DstMucTrack::m_brLastLayer, m_changeUnit, DstMucTrack::m_chi2, DstMucTrack::m_deltaPhi, DstMucTrack::m_depth, DstMucTrack::m_dof, DstMucTrack::m_ecLastLayer, DstMucTrack::m_endPart, DstMucTrack::m_id, DstMucTrack::m_kalbrLastLayer, DstMucTrack::m_kaldepth, DstMucTrack::m_kaldof, DstMucTrack::m_kalecLastLayer, DstMucTrack::m_kalrechi2, DstMucTrack::m_maxHitsInLayer, DstMucTrack::m_numHits, DstMucTrack::m_numLayers, m_recmode, DstMucTrack::m_rms, DstMucTrack::m_startPart, DstMucTrack::m_status, DstMucTrack::m_type, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPosSigma, and DstMucTrack::m_zPosSigma.

  : //m_trackId        (source.m_trackId),    //--------------->
    m_ExtTrackID   (source.m_ExtTrackID),
    m_MdcPos       (source.m_MdcPos),
    m_MdcMomentum  (source.m_MdcMomentum),
    m_MucPos       (source.m_MucPos),
    m_MucPosSigma  (source.m_MucPosSigma),
    m_MucMomentum  (source.m_MucMomentum),
    m_CurrentPos   (source.m_CurrentPos),
    m_CurrentDir   (source.m_CurrentDir),
    m_CurrentInsct (source.m_CurrentInsct),
    m_pHits        (source.m_pHits),
    m_pExpectedHits(source.m_pExpectedHits),
    m_Intersections(source.m_Intersections),
    m_Directions   (source.m_Directions)
{ 
    m_id           = source.m_id;
    m_status       = source.m_status;
    m_type         = source.m_type;
    m_numHits      = source.m_numHits;    //--------------->
    m_startPart    = source.m_startPart;
    m_endPart      = source.m_endPart;
    m_brLastLayer  = source.m_brLastLayer;
    m_ecLastLayer  = source.m_ecLastLayer;
    m_numLayers    = source.m_numLayers;
    m_maxHitsInLayer= source.m_maxHitsInLayer;
    m_depth        = source.m_depth;
    m_dof          = source.m_dof;
    m_chi2         = source.m_chi2;
    m_rms          = source.m_rms;
    m_deltaPhi     = source.m_deltaPhi;
    m_xPosSigma    = source.m_xPosSigma;
    m_yPosSigma    = source.m_yPosSigma;
    m_zPosSigma    = source.m_zPosSigma;
    m_changeUnit   = source.m_changeUnit;
    m_recmode      = source.m_recmode;
    //******
    m_kalrechi2     = source.m_kalrechi2;
    m_kaldof        = source.m_kaldof;
    m_kaldepth      = source.m_kaldepth;
    m_kalbrLastLayer= source.m_kalbrLastLayer;
    m_kalecLastLayer= source.m_kalecLastLayer;
}

RecMucTrack::RecMucTrack

(

const DstMucTrack &

dstTrack

)

Copy constructor from DstMucTrack.

Definition at line 178 of file RecMucTrack.cxx.

References SetDefault().

  :DstMucTrack(dstTrack)
{

  SetDefault();
  
}

RecMucTrack::~RecMucTrack

(

)

Destructor.

Definition at line 201 of file RecMucTrack.cxx.

References genRecEmupikp::i, and m_pExpectedHits.

{ 
for(int i = 0 ; i < m_pExpectedHits.size(); i++)
          delete m_pExpectedHits[i];
}

---

Member Function Documentation

void RecMucTrack::AttachDirection

(

Hep3Vector

dir

)

Attach the direction to this trajectory.

Definition at line 614 of file RecMucTrack.cxx.

References m_Directions.

{
  m_Directions.push_back(dir);
}

void RecMucTrack::AttachHit

(

MucRecHit *

hit

)

Attach the given hit to this track.

Definition at line 392 of file RecMucTrack.cxx.

References MucRecHit::Gap(), MucID::getGapNum(), m_pHits, and MucRecHit::Part().

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  //  cout << "Muc2DRoad::AttachHit-I0  hit = " << hit << endl;

  if (!hit) {
    cout << "RecMucTrack::AttachHit-E1  null hit pointer!" << endl;
    return ;
  }
  
  int part = hit->Part();
  int gap  = hit->Gap();
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    // The gap number of the hit is out of range.
    cout << "Muc2DRoad::AttachHit(MucRecHit*), bad gap number = " << gap
         << endl;
    return;
  }

  // Attach the hit to the road.
  m_pHits.push_back(hit);
  
  //  m_HitDistance[gap] = dX;

  // Now recalculate the total number of hits and the max. number of
  // hits per gap.
  //CountHits();
}

void RecMucTrack::AttachInsct

(

Hep3Vector

insct

)

Attach the intersection to this trajectory.

Definition at line 608 of file RecMucTrack.cxx.

References m_Intersections.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  m_Intersections.push_back(insct);
}

int RecMucTrack::brFirstLayer

(

)

const [inline]

Which gap on Barrel is the first one with hits attached to this track?

Definition at line 221 of file RecMucTrack.h.

References m_brFirstLayer.

Referenced by ComputeDepth().

{ return m_brFirstLayer; }

int DstMucTrack::brLastLayer

(

)

const [inline, inherited]

Definition at line 39 of file DstMucTrack.h.

References DstMucTrack::m_brLastLayer.

Referenced by ComputeDepth().

{ return   m_brLastLayer;  } // last layer with hits in barrel 

Hep3Vector RecMucTrack::CalculateInsct	(	const int	part,
		const int	seg,
		const int	gap
	)

Calculate intersection from all hits attached on this gap.

Definition at line 541 of file RecMucTrack.cxx.

References MucGeoGeneral::GetGap(), genRecEmupikp::i, MucGeoGeneral::Instance(), kInsctWeight, m_CurrentInsct, m_pHits, MucGeoGap::Orient(), MucGeoGap::TransformToGap(), MucGeoGap::TransformToGlobal(), and x.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  MucGeoGap *gapPtr = MucGeoGeneral::Instance()->GetGap(part, seg, gap);
  vector<int> hitSeq;
  for(int i = 0; i < (int)m_pHits.size(); i++) {
    MucRecHit *aHit = m_pHits[i];
    if(aHit->Part() == part &&
       aHit->Seg()  == seg  &&
       aHit->Gap()  == gap) {
      hitSeq.push_back(i);
    }  
  }    
  int nHitInGap = hitSeq.size();
  //cout << "nHitInGap " << nHitInGap << endl;
  
  HepPoint3D insctLocal = gapPtr->TransformToGap(m_CurrentInsct);
  //HepPoint3D newInsct(0.0, 0.0, 0.0);
  HepPoint3D newInsctLocal = insctLocal;

  vector<float> x;
  for(int i = 0; i < nHitInGap; i++) x.push_back(0);
  float xInsct = 0, xNewInsct = 0;


  int orient = gapPtr->Orient();
  if(orient == 1) xInsct = insctLocal.x();
  if(orient == 0) xInsct = insctLocal.y();
  
  for(int i = 0; i < nHitInGap; i++) {
    float xStrip, yStrip, zStrip;
    (m_pHits[hitSeq[i]])->GetStrip()->GetCenterPos(xStrip, yStrip, zStrip);
    if(orient == 1) x[i] = xStrip;
    if(orient == 0) x[i] = yStrip;
  }
  //cout << "local Old" << insctLocal << endl;

  //if == 0, no direction change;
  if(nHitInGap > 0) {
    xNewInsct = xInsct * kInsctWeight;

    //float minDist = kInfinity;
    for(int i = 0; i < nHitInGap; i++) {
      //if(fabs(x[i] - xInsct) < minDist) {
      //xNewInsct = x[i];    //}
      xNewInsct += x[i] * ((1.0 - kInsctWeight) / nHitInGap);
    }
    
    if(orient == 1) {
      newInsctLocal.setX(xNewInsct);
      newInsctLocal.setY(insctLocal.y());
    }
    if(orient == 0) {
      newInsctLocal.setX(insctLocal.x());
      newInsctLocal.setY(xNewInsct);
    }
  }

  m_CurrentInsct = gapPtr->TransformToGlobal(newInsctLocal);
  //cout << "local  New" << newInsctLocal << endl;
  //cout << "global New" << m_CurrentInsct << endl;
  
  return m_CurrentInsct;
}

double DstMucTrack::chi2

(

)

const [inline, inherited]

Definition at line 46 of file DstMucTrack.h.

References DstMucTrack::m_chi2.

Referenced by MucRecTrkExt::execute(), and MucPID::particleIDCalculation().

{ return   m_chi2;   }

static const CLID& RecMucTrack::classID

(

)

[inline, static]

Reimplemented from DstMucTrack.

Definition at line 49 of file RecMucTrack.h.

References CLID_RecMucTrack.

Referenced by clID().

    {
      return CLID_RecMucTrack;
    }

virtual const CLID& RecMucTrack::clID

(

)

const [inline, virtual]

Reimplemented from DstMucTrack.

Definition at line 44 of file RecMucTrack.h.

References classID().

    {
      return RecMucTrack::classID();
    }

int RecMucTrack::ComputeDepth

(

int

method

)

Compute depth.

Definition at line 673 of file RecMucTrack.cxx.

References brFirstLayer(), DstMucTrack::brLastLayer(), cos(), DstMucTrack::deltaPhi(), DstMucTrack::ecLastLayer(), MucGeoGeneral::GetGap(), MucID::getGapMax(), MucGeoGap::GetIronThickness(), MucGeoGeneral::Instance(), kPi, DstMucTrack::m_brLastLayer, DstMucTrack::m_depth, m_depth_3, DstMucTrack::m_ecLastLayer, m_Good3DLine, m_Good3DPart, m_IntersectionInner, m_IntersectionOuter, m_MucMomentum, DstMucTrack::m_numLayers, m_pHits, sin(), and x.

{
  if( m_numLayers == 0 ) {
    m_depth = m_depth_3 = -99; return 0; 
  } 
  else if( m_numLayers == 1 && (m_brLastLayer == 0 || m_ecLastLayer == 0) ) {
    m_depth = m_depth_3 = 0; return 0;
  }
  
  m_depth_3 = 0.0;
  
  float brThickness = 0.0;  float ecThickness = 0.0;  float deltaPhi = 0.0;
  int betweenSeg = 0;

  float phi = m_MucMomentum.phi();
  float theta = m_MucMomentum.theta();

  vector<MucRecHit*>::const_iterator iHit;
  int ngap = MucID::getGapMax();
  //cout<<"ngap:\t"<<ngap<<endl;
  
  int Seg[9]; int part, seg, gap, strip;
  for(int gap = 0; gap < ngap; gap++){Seg[gap] = -1;}
  for(iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if(*iHit) {  // Check for a null pointer.
      part   = (*iHit)->Part();
      seg    = (*iHit)->Seg();
      gap    = (*iHit)->Gap();
      strip  = (*iHit)->Strip();
      if(part==1) Seg[gap] = seg;
    }
  }

  int segTrackBr = -1;
  for(int gap = 0; gap <= brLastLayer(); gap++){
    if(Seg[gap] != -1) segTrackBr = Seg[gap];
  }
  // BR, the innermost layer is RPC module
  for(int gap = 0; gap <= brLastLayer(); gap++){
    float thickness = 0.0;
    if(Seg[gap] != -1 && Seg[brLastLayer()-1] != -1 && Seg[gap] != Seg[brLastLayer()-1]) betweenSeg = 1; 
    thickness = MucGeoGeneral::Instance()->GetGap(1, 0, gap)->GetIronThickness();
    //cout<<"RecMucTrack  gap="<<gap<<" brlastgap="<<brLastLayer()<<" "<<thickness<<endl;
    if(sin(m_MucMomentum.theta()) != 0 ) thickness /= sin(m_MucMomentum.theta());
    else ;//cout<<"RecMucTrack::ComputeDepth, In Barrel,theta=0?"<<endl;
    deltaPhi = m_MucMomentum.phi() - segTrackBr*kPi/4;
    if(Seg[gap] == -1 && betweenSeg == 1) {
      thickness += 40;  // gap width
    }
    if(cos(deltaPhi) != 0 ) thickness /= cos(deltaPhi);
    else ;//cout<<"RecMucTrack::ComputeDepth, In Barrel,Cos(phi)=0?"<<endl;
    if(deltaPhi == 0 && Seg[brLastLayer()-1]==2 ) thickness = 0;
    //cout<<"in muctrack "<<thickness<<" "<<brThickness<<" theta "<<m_MucMomentum.theta()<<" phi="<<deltaPhi<<" "<<m_MucMomentum.phi()<<endl;

    if(brFirstLayer()<brLastLayer())  brThickness += thickness; 
    else if(brFirstLayer()==brLastLayer()) {  //only one gap
      if(m_MucMomentum.mag()>1000 || brFirstLayer()<2)  brThickness += thickness;  //high momentum or only one or two gap
      //cout<<"mom="<<m_MucMomentum.mag()<<" "<<brFirstLayer()<<" "<<brThickness<<" "<<thickness<<endl;
    }
    else cout<<"in RecMucTrack: Wrong Gap Info"<<endl;

    //cout<<brThickness<<endl;
  }

  //cout<<"eclastgap= "<<ecLastLayer()<<"  ecfirstgap= "<<ecFirstLayer()<<endl;

  // EC, the innermost layer is Iron
  //for (int gap = ecFirstLayer(); gap!=-1&&gap <= ecLastLayer(); gap++) {
  for (int gap = 0; gap!=-1&&gap <= ecLastLayer(); gap++) {
    ecThickness += MucGeoGeneral::Instance()->GetGap(0, 0, gap)->GetIronThickness();
  }
  
  if (cos(theta) != 0.0) ecThickness /= cos(theta);
  else ;//cout << "RecMucTrack::ComputeDepth, In EndCap, Track theta = 90.0 ? " << endl;
  ecThickness  = fabs(ecThickness);
  
  //cout<<"eclastgap= "<<ecLastLayer()<<" ecthickness="<<ecThickness<<endl;

  if(method == 2){
    //barrel first
    if((m_Good3DLine == 1) &&(m_Good3DPart == 1)){
      if(m_IntersectionInner[0].x()!=-9999){
        for(int gap = 1;gap <= brLastLayer(); gap++)//Inner[gap1]-Outer[gap0] ...
          {
          if(m_IntersectionInner[gap].x() != -9999 && m_IntersectionInner[gap-1].x() != -9999)
            m_depth_3 += (m_IntersectionInner[gap] - m_IntersectionOuter[gap-1]).mag();
          }
        //may be pass some gap in endcap!
        m_depth_3 += ecThickness;
      }
      else m_depth_3 = brThickness + ecThickness;
    }
    if((m_Good3DLine == 1) &&(m_Good3DPart != 1)){
      for(int gap = 1;gap <= ecLastLayer(); gap++)//Inner[gap1]-Outer[gap0] ...
      {
        if(m_IntersectionInner[gap].x() != -9999 && m_IntersectionInner[gap-1].x() != -9999)
          m_depth_3 += (m_IntersectionInner[gap] - m_IntersectionOuter[gap-1]).mag();
      }
      //may be pass some gap in barrel!
      m_depth_3 += brThickness;
      if (cos(theta) != 0.0) m_depth_3 += 40/cos(theta); //there is 1 absorber before first gap in endcap!
    }
    if(m_Good3DLine == 0) m_depth_3 = brThickness + ecThickness;
    if(m_depth_3>2000) m_depth_3 = brThickness + ecThickness;  //unreasonable depth! so use previous method!
  }
  else  //method == 1 linefit
    {
      m_depth_3 = brThickness + ecThickness;

    }
  
  double offset = 50.0;
  //if(GetTotalHits() > 0) m_depth_3 += offset;

  m_depth = m_depth_3; 

  //if(m_depth<0||m_depth>2000) m_depth = 0;
  if(m_depth>2000) m_depth = -99;
  //cout<<"depth= "<<m_depth<<endl;

}

void RecMucTrack::ComputeDepth

(

)

Compute depth.

Definition at line 796 of file RecMucTrack.cxx.

References DstMucTrack::brLastLayer(), cos(), DstMucTrack::deltaPhi(), DstMucTrack::ecLastLayer(), MucGeoGeneral::GetGap(), MucID::getGapMax(), MucGeoGap::GetIronThickness(), GetTotalHits(), HasHitInGap(), MucGeoGeneral::Instance(), kPi, DstMucTrack::m_depth, m_depth_3, m_MucMomentum, m_pHits, and sin().

Referenced by ComputeTrackInfo().

{
  m_depth = -99.0;
  // Part 1 first.
  float brThickness = 0.0;
  for (int gap = 0; gap <= brLastLayer(); gap++) {
    brThickness += MucGeoGeneral::Instance()->GetGap(1, 0, gap)->GetIronThickness();
  }

  // second alg
  float brThickness_2 = 0.0;
  for (int gap = 0; gap <= brLastLayer(); gap++) {
    if(HasHitInGap(1,gap)){
      brThickness_2 += MucGeoGeneral::Instance()->GetGap(1, 0, gap)->GetIronThickness();
    }
  }
  // third alg
  float brThickness_3 = 0.0;
  vector<MucRecHit*>::const_iterator iHit;
  int ngap = MucID::getGapMax();
  int Seg[9]; int part, seg, gap, strip;
  for(int gap = 0; gap < ngap; gap++){Seg[gap] = -1;}
  for(iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if(*iHit) {  // Check for a null pointer.
      part   = (*iHit)->Part();
      seg    = (*iHit)->Seg();
      gap    = (*iHit)->Gap();
      strip  = (*iHit)->Strip();
      if(part==1) Seg[gap] = seg;
    }
  }
  
  float deltaPhi_3 = 0.0;
  int betweenSeg = 0;

  int segTrackBr = -1;
  for(int gap = 0; gap <= brLastLayer(); gap++){
    if(Seg[gap] != -1) segTrackBr = Seg[gap];
  }

  for(int gap = 0; gap <= brLastLayer(); gap++){
    float thickness = 0.0;
    if(Seg[gap] != -1 && Seg[brLastLayer()-1] != -1 && Seg[gap] != Seg[brLastLayer()-1]) betweenSeg = 1; 
    thickness = MucGeoGeneral::Instance()->GetGap(1, 0, gap)->GetIronThickness();
    if(sin(m_MucMomentum.theta()) != 0 ) thickness /= sin(m_MucMomentum.theta());
    else cout<<"RecMucTrack::ComputeDepth, In Barrel,theta=0?"<<endl;
    //if(Seg[gap] != -1) deltaPhi_3 = m_MucMomentum.phi() - Seg[gap]*kPi/4;
    deltaPhi_3 = m_MucMomentum.phi() - segTrackBr*kPi/4;   //some times, no hits in a gap, but a good track exist!
    if(Seg[gap] == -1 && betweenSeg == 1) {
      cout<<"between segment"<<endl;
      thickness += 40;  // gap width
    }

    if(cos(deltaPhi_3) != 0 ) thickness /= cos(deltaPhi_3);
    else cout<<"RecMucTrack::ComputeDepth, In Barrel,Cos(phi)=0?"<<endl;

    if(deltaPhi_3 == 0 && Seg[brLastLayer()-1]==2 ) thickness = 0;
    //cout<<"in muctrack "<<thickness<<" "<<brThickness_3<<" theta "<<m_MucMomentum.theta()<<" phi="<<deltaPhi_3<<" "<<m_MucMomentum.phi()<<endl;
    brThickness_3 += thickness;
    
  }

  //cout<<"in RecMucTrack: compare thickness "<<brThickness<<" "<<brThickness_2<<" "<<brThickness_3<<endl;

  float phi = m_MucMomentum.phi();
  float deltaPhi = phi - kPi/4*(int)(phi/(kPi/4));
  if (deltaPhi > kPi/8) deltaPhi = kPi/4 - deltaPhi;
  float theta = m_MucMomentum.theta();
//   cout << "br LastGap " << brLastLayer() << " Thick " << brThickness 
//     << " 1/sin(theta) " << 1/sin(theta) << " 1/cos(deltaPhi) " << 1/cos(deltaPhi) << endl;

  if (sin(theta) != 0.0) brThickness /= sin(theta);
  else cout << "RecMucTrack::ComputeDepth, In Barrel, Track theta = 0.0 ? " << endl; 

  brThickness /= cos(deltaPhi);
  brThickness  = fabs(brThickness);
  //cout << "br Depth " << brThickness << endl;
  
  // EC, then
  float ecThickness = 0.0;
  for (int gap = 0; gap <= ecLastLayer(); gap++) {
    ecThickness += MucGeoGeneral::Instance()->GetGap(0, 0, gap)->GetIronThickness();
  }
  //cout << "ec LastGap " << ecLastLayer() << " Thick " << ecThickness
  //   << " 1/cos(theta) " << 1/cos(theta) << endl;
  
  if (cos(theta) != 0.0) ecThickness /= cos(theta);
  else ; //cout << "RecMucTrack::ComputeDepth, In EndCap, Track theta = 90.0 ? " << endl;
  ecThickness  = fabs(ecThickness);
  //cout << "ec Depth " << ecThickness << endl;

  m_depth = brThickness + ecThickness;
  m_depth_3 = brThickness_3 + ecThickness;
  cout << "Depth " << m_depth << " Depth3 = "<<m_depth_3<<endl;
  m_depth = m_depth_3;
  double offset = 50.0;
  if(GetTotalHits() > 0) m_depth += offset; // since brThickness on gap 0 is zero, give an offset for track arriving muc.

}

void RecMucTrack::ComputeDistanceMatch

(

)

Compute distance match //2006.11.08.

Definition at line 897 of file RecMucTrack.cxx.

References DstMucTrack::distance(), FindSegWithMaxHits(), MucGeoGap::Gap(), MucGeoGeneral::GetGap(), MucGeoGeneral::Instance(), m_Dist_muc_ext, DstMucTrack::m_distance, m_ExtMucPos, m_MucPos, m_pHits, m_StripPhi, MucGeoGap::Orient(), MucGeoGap::Part(), MucGeoGap::Seg(), and MucGeoGap::TransformToGap().

Referenced by ComputeTrackInfo().

{
  bool firstHitFound = false;
  MucGeoGap *firstGap = 0; 
  vector<MucRecHit*>::const_iterator iHit;
  vector<MucRecHit*> hitsGap0;
  float stripLocal[3]={0.0, 0.0, 0.0};
  float stripGlobal[3]={0.0, 0.0, 0.0};
  int nStrip = 0;

  int part, seg, gap, strip;
  int barrel_gap0_exist = 0;

  for(iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if(*iHit) {  // Check for a null pointer.
      part   = (*iHit)->Part();
      seg    = (*iHit)->Seg();
      gap    = (*iHit)->Gap();
      strip  = (*iHit)->Strip();
      if(!firstHitFound && gap == 0) {
        firstGap = MucGeoGeneral::Instance()->GetGap(part, seg, gap);
        firstHitFound = true;
      }
      if(firstGap && part == firstGap->Part() && seg == firstGap->Seg() && gap == firstGap->Gap()) {
        //cout<<"in RecMucTrack "<<part<<" "<<seg<<" "<<gap<<" "<<strip<<endl;
        HepPoint3D posHit        = (*iHit)->GetCenterPos();
        HepPoint3D posHitLocal   = (*iHit)->GetGap()->TransformToGap(posHit);
        if(part==1&&gap==0) barrel_gap0_exist = 1;  //exist
        
        stripLocal[0] += posHitLocal.x();
        stripLocal[1] += posHitLocal.y();
        stripLocal[2] += posHitLocal.z(); 

        stripGlobal[0] += posHit.x();   //to calc phi of this strip 
        stripGlobal[1] += posHit.y();
        stripGlobal[2] += posHit.z();
 
        nStrip++;
      }
    }
  }

  //cout<<"in RecMucTrack: extpos "<<m_ExtMucPos<<" mucpos "<< m_MucPos<<endl;
  
  int apart = -1, aseg = -1;
  int nHits = FindSegWithMaxHits(apart, aseg);
  MucGeoGap *fakefirstGap = 0;   // maybe not exist!
  if(apart == -1 && aseg== -1)
    {m_Dist_muc_ext.set(0,0,0);}
  else {
    fakefirstGap = MucGeoGeneral::Instance()->GetGap(apart, aseg, 0);
    HepPoint3D fextLocal = fakefirstGap->TransformToGap(m_ExtMucPos);
    HepPoint3D fmucLocal = fakefirstGap->TransformToGap(m_MucPos);
    float dist_x = fextLocal.x() - fmucLocal.x();
    float dist_y = fextLocal.y() - fmucLocal.y();
    float dist_z = fextLocal.z() - fmucLocal.z();
    m_Dist_muc_ext.set(dist_x,dist_y,dist_z);
    //cout<<"in RecMucTrack dist = "<<dist_x<<" "<<dist_y<<" "<<dist_z<<endl;
    if (fakefirstGap->Orient() == 0) { // part 0,2
      //cout<<"in RecMucTrack "<< extLocal.y()<<" "<<stripLocal[1]<<endl;
    }
    else {
      //cout<<"in RecMucTrack1 "<< extLocal.x()<<" "<<stripLocal[0]<<endl;
    }
  }

  float distance = -9990;

  if (nStrip == 0 || !firstGap) {
    distance = -9990;
    //cout << "no hits on first gap" << endl;
  }  
  else{
    for (int k = 0; k < 3; k++) stripLocal[k] /= nStrip;
    for (int k = 0; k < 3; k++) stripGlobal[k] /= nStrip;
    
    m_StripPhi.set(stripGlobal[0],stripGlobal[1],stripGlobal[2]);

    HepPoint3D extLocal = firstGap->TransformToGap(m_ExtMucPos);
    // extmucpos may be can be used to identify mu/pion

//     cout<<"in RecMucTrack mom "<<m_MucMomentum.x()<<" "<<m_MucMomentum.y()<<" "<<m_MucMomentum.z()<<endl;
//     cout<<"in RecMucTrack extpos "<<m_ExtMucPos.x()<<" "<<m_ExtMucPos.y()<<" "<<m_ExtMucPos.z()<<endl;
//     cout<<"in RecMucTrack extpos2 "<<ExtMucPos_2.x()<<" "<<ExtMucPos_2.y()<<" "<<ExtMucPos_2.z()<<endl;
//     cout<<"in RecMucTrack extloc "<<extLocal.x()<<" "<<extLocal.y()<<" "<<extLocal.z()<<endl;
    if (firstGap->Orient() == 0) { // part 0,2
      distance = extLocal.y() - stripLocal[1];
      //cout<<"in RecMucTrack "<< extLocal.y()<<" "<<stripLocal[1]<<endl;
    }
    else {
      distance = extLocal.x() - stripLocal[0];
      //cout<<"in RecMucTrack1 "<< extLocal.x()<<" "<<stripLocal[0]<<endl;
    }
  }

  m_distance = distance;

  //use m_chi2 temporary 
  //m_chi2 = distance;
  // cout<<"in RecMucTrack  distance= "<<m_distance<<" n= "<<nStrip<<endl;
}

void RecMucTrack::ComputeLastGap

(

)

Comute last gap in barrel and endcap.

Definition at line 635 of file RecMucTrack.cxx.

References m_brFirstLayer, DstMucTrack::m_brLastLayer, m_ecFirstLayer, DstMucTrack::m_ecLastLayer, m_ecPart, DstMucTrack::m_endPart, and m_pHits.

Referenced by ComputeTrackInfo(), and LineFit().

{
  int lastGap1, lastGap2;
  int firstGap1, firstGap2;
  lastGap1 = lastGap2 = -1;
  firstGap1 = firstGap2 = 99;
  vector<MucRecHit*>::const_iterator iHit;
  iHit = m_pHits.begin();
  for( ;iHit != m_pHits.end(); iHit++) 
  {
    if(*iHit) 
    {  // Check for a null pointer.
      int part = (*iHit)->Part();
      int gap  = (*iHit)->Gap();
      
      if(part == 1) {
        if(gap > lastGap1)  lastGap1 = gap;      
        if(gap < firstGap1) firstGap1 = gap;
      }
      else {
        if(gap > lastGap2)  { m_ecPart = part; m_endPart = part; lastGap2 = gap; }
        if(gap < firstGap2) firstGap2 = gap;
      }
    }
  }
  
  m_brLastLayer = lastGap1;
  if(firstGap1 == 99) m_brFirstLayer = -1;
  else if(firstGap1>=0&&firstGap1<9) m_brFirstLayer = firstGap1;

  m_ecLastLayer = lastGap2;
  if(firstGap2 == 99) m_ecFirstLayer = -1;
  else if(firstGap2>=0&&firstGap2<8) m_ecFirstLayer = firstGap2;

  //cout<<"MucTrack, br: "<<m_brFirstLayer<<", "<<m_brLastLayer<<"\tec: "<<m_ecFirstLayer<<", "<<m_ecLastLayer<<endl;
}

void RecMucTrack::ComputeMaxHitsInGap

(

)

ComputeMaxHitsInGap;.

Definition at line 1528 of file RecMucTrack.cxx.

References MucID::getGapNum(), GetHitInGap(), MucID::getPartNum(), and DstMucTrack::m_maxHitsInLayer.

Referenced by ComputeTrackInfo().

{
  int maxHits = 0;
  int hits = 0;
  for(int part = 0; part < (int)MucID::getPartNum(); part++) {
    for(int gap = 0; gap < (int)MucID::getGapNum(part); gap++) {
      hits = GetHitInGap(part, gap);
      if(hits > maxHits) maxHits = hits;
    }
  }

  m_maxHitsInLayer = maxHits;
}

void RecMucTrack::ComputeNGapsWithHits

(

)

ComputeNGapsWithHits;.

Definition at line 1412 of file RecMucTrack.cxx.

References count, MucID::getGapMax(), DstMucTrack::m_numHits, DstMucTrack::m_numLayers, and m_pHits.

Referenced by ComputeTrackInfo().

{
  int ngap = MucID::getGapMax();
  int gap, count = 0;
  vector<int> firedGap;
  for(gap = 0; gap < ngap; gap++) firedGap.push_back(0);

  vector<MucRecHit*>::const_iterator iHit;
  for (iHit=m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      gap = (*iHit)->Gap();
      firedGap[gap] = 1;
    }
  }
  
  for ( gap = 0; gap < ngap; gap++) {
    count += firedGap[gap];
  }

  m_numHits   = m_pHits.size();
  m_numLayers = count;
  //cout<<"nLayer:\t"<<m_numLayers<<endl;
}

void RecMucTrack::ComputeTrackInfo

(

int

fittingmethod

)

Compute all infomation of this track;.

Definition at line 1670 of file RecMucTrack.cxx.

References ComputeDepth(), ComputeDistanceMatch(), ComputeLastGap(), ComputeMaxHitsInGap(), ComputeNGapsWithHits(), m_pExpectedHits, m_pHits, OutputUnitChange(), setExpHits(), and setVecHits().

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  setVecHits(m_pHits);
  setExpHits(m_pExpectedHits);
  ComputeLastGap();
  ComputeNGapsWithHits(); 
  ComputeMaxHitsInGap();
  //ComputeDepth();
  ComputeDepth(fittingmethod);
  ComputeDistanceMatch();
  OutputUnitChange();
}

void RecMucTrack::CorrectDir

(

)

Correct direction of this trajectory.

Definition at line 620 of file RecMucTrack.cxx.

References m_CurrentDir, m_CurrentInsct, and m_CurrentPos.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  //cout << "Before CorrectDir(), fCurrentDir " << m_CurrentDir << endl;
  m_CurrentDir = m_CurrentInsct - m_CurrentPos;
  m_CurrentDir.setMag(1.0);
  //cout << "After CorrectDir(), fCurrentDir " << m_CurrentDir << endl;
}

void RecMucTrack::CorrectPos

(

)

Correct current position of this trajectory.

Definition at line 629 of file RecMucTrack.cxx.

References m_CurrentInsct, and m_CurrentPos.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  m_CurrentPos = m_CurrentInsct;
}

double DstMucTrack::deltaPhi

(

)

const [inline, inherited]

Definition at line 63 of file DstMucTrack.h.

References DstMucTrack::m_deltaPhi.

Referenced by ComputeDepth(), GetMdcExtTrack(), LineFit(), and MucPID::particleIDCalculation().

{ return m_deltaPhi;}

double DstMucTrack::depth

(

)

const [inline, inherited]

Definition at line 45 of file DstMucTrack.h.

References DstMucTrack::m_depth.

Referenced by MucRecTrkExt::execute(), BestDTagSvc::isMuon(), DTagTool::isMuon(), and MucPID::particleIDCalculation().

{ return   m_depth;  }    // depth of the track transport in iron

double DstMucTrack::distance

(

)

const [inline, inherited]

Definition at line 62 of file DstMucTrack.h.

References DstMucTrack::m_distance.

Referenced by ComputeDistanceMatch(), MucRecTrkExt::execute(), GetHitDistance(), GetHitDistance2(), and MucPID::particleIDCalculation().

{ return m_distance;}

int DstMucTrack::dof

(

)

const [inline, inherited]

Definition at line 47 of file DstMucTrack.h.

References DstMucTrack::m_dof.

{ return   m_dof;    }

int RecMucTrack::ecFirstLayer

(

)

const [inline]

Which gap on EndCap is the first one with hits attached to this track?

Definition at line 224 of file RecMucTrack.h.

References m_ecFirstLayer.

{ return m_ecFirstLayer; }

int DstMucTrack::ecLastLayer

(

)

const [inline, inherited]

Definition at line 40 of file DstMucTrack.h.

References DstMucTrack::m_ecLastLayer.

Referenced by ComputeDepth().

{ return   m_ecLastLayer;  } // last layer with hits in end cap

int DstMucTrack::endPart

(

)

const [inline, inherited]

Definition at line 38 of file DstMucTrack.h.

References DstMucTrack::m_endPart.

{ return   m_endPart; }      // end   position of track locates in which part

void RecMucTrack::Extend

(

)

Extend mucpos and extmucpos to first layer of muc.

Definition at line 1000 of file RecMucTrack.cxx.

References count, m_ExtMucPos, m_MucMomentum, m_MucPos, DstMucTrack::m_xPos, DstMucTrack::m_yPos, DstMucTrack::m_zPos, and x.

Referenced by LineFit().

{
  float x = m_ExtMucPos.x(), y = m_ExtMucPos.y(), z = m_ExtMucPos.z();
  //cout<<"in MucTrac, MucPos= "<<x<<" "<<y<<" "<<z<<endl;
  float step = 0.005;
  float myMucR = 1720.0;
  float myMucZ = 2110.0;
   int count_=0; //added by LI Chunhua to avoid loop infinitely
   while ( ( (fabs(x)<=myMucR) && (fabs(y)<=myMucR) && ((fabs(x-y)/sqrt(2.))<=myMucR) && ((fabs(x+y)/sqrt(2.))<=myMucR) && (fabs(z)<=myMucZ) )&&count_<1000 ) {
    x += step * m_MucMomentum.x();
    y += step * m_MucMomentum.y();
    z += step * m_MucMomentum.z();
    count_++;
    //cout << "in RecMucTrack+ x " << x << " y " << y << " z " << z <<" mom "<< m_MucMomentum.x()<< " "<<m_MucMomentum.y()<<" "<<m_MucMomentum.z()<<endl;
  }

   int count = 0; 
  while( !( (fabs(x)<=myMucR) && (fabs(y)<=myMucR) && ((fabs(x-y)/sqrt(2.))<=myMucR) && ((fabs(x+y)/sqrt(2.))<=myMucR) && (fabs(z)<=myMucZ) )&&count<1000 ) {
    x -= step * m_MucMomentum.x();
    y -= step * m_MucMomentum.y();
    z -= step * m_MucMomentum.z();
    count++;
    //cout << "in RecMucTrack- x " << x << " y " << y << " z " << z <<" mom "<< m_MucMomentum.x()<< " "<<m_MucMomentum.y()<<" "<<m_MucMomentum.z()<<endl;
  } 
 
  if(count<1000&&count_<1000){
    if(fabs(x)<2600&&fabs(y)<2600&&fabs(z)<2800){
      m_ExtMucPos.set(x, y, z);
      m_MucPos.set(x,y,z);
      m_xPos = x;
      m_yPos = y;
      m_zPos = z;
    }
  }
 
}

int RecMucTrack::FindSegWithMaxHits	(	int &	part,
		int &	seg
	)

Find the segment which contains most hits, return max hits number.

Definition at line 1509 of file RecMucTrack.cxx.

References GetHitInSeg(), MucID::getPartNum(), and MucID::getSegNum().

Referenced by ComputeDistanceMatch(), and LineFit().

{
  int maxHits = 0;
  int hits = 0;
  for(int part = 0; part < (int)MucID::getPartNum(); part++) {
    for(int seg = 0; seg < (int)MucID::getSegNum(part); seg++) {
      hits = GetHitInSeg(part, seg);
      if (hits > maxHits) {
        maxHits  = hits;
        findPart = part;
        findSeg  = seg;
      }
    }
  }

  return maxHits;
}

Hep3Vector RecMucTrack::GetCurrentDir

(

)

const [inline]

Current direction.

Definition at line 212 of file RecMucTrack.h.

References m_CurrentDir.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{ return m_CurrentDir/1000; }  

Hep3Vector RecMucTrack::GetCurrentInsct

(

)

const [inline]

Current intersection.

Definition at line 215 of file RecMucTrack.h.

References m_CurrentInsct.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{ return m_CurrentInsct; } //internal  

Hep3Vector RecMucTrack::GetCurrentPos

(

)

const [inline]

Current position.

Definition at line 209 of file RecMucTrack.h.

References m_CurrentPos.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{ return m_CurrentPos/10; }  

float RecMucTrack::GetDepth3

(

)

const [inline]

Length of the track penetrating in iron absorber.

Definition at line 252 of file RecMucTrack.h.

References m_depth_3.

{ return m_depth_3/10; }

vector<float> RecMucTrack::getDistHits

(

)

const [inline]

Definition at line 311 of file RecMucTrack.h.

References m_distHits.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{ return  m_distHits ;}

vector< MucRecHit * > RecMucTrack::GetExpectedHits

(

)

const

Definition at line 1643 of file RecMucTrack.cxx.

References m_pExpectedHits.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
    return m_pExpectedHits;
}

vector<int> RecMucTrack::getExpHits

(

)

const [inline]

Definition at line 309 of file RecMucTrack.h.

References m_expHits.

{ return  m_expHits ;}

vector<float> RecMucTrack::getExtDistHits

(

)

const [inline]

Definition at line 315 of file RecMucTrack.h.

References m_distHits_ext.

Referenced by MucRecTrkExt::execute().

{ return  m_distHits_ext ;} 

Hep3Vector RecMucTrack::GetExtMucDistance

(

)

const [inline]

Distance match of the ext track with muc track in first layer.

Definition at line 259 of file RecMucTrack.h.

References m_Dist_muc_ext.

Referenced by MucRecTrkExt::execute().

{ return m_Dist_muc_ext;}

Hep3Vector RecMucTrack::getExtMucMomentum

(

)

const [inline]

Start momentum of the Ext track in Muc.

Definition at line 195 of file RecMucTrack.h.

References m_ExtMucMomentum.

{ return m_ExtMucMomentum/1000; }

Hep3Vector RecMucTrack::getExtMucPos

(

)

const [inline]

start position of the Ext track in Muc.

Definition at line 192 of file RecMucTrack.h.

References m_ExtMucPos.

{ return m_ExtMucPos/10; }

int RecMucTrack::getExtTrackID

(

)

const [inline]

The identifier of the Ext track as seed.

Definition at line 189 of file RecMucTrack.h.

References m_ExtTrackID.

{ return m_ExtTrackID; }

MucRecHit * RecMucTrack::GetHit	(	const int	part,
		const int	gap
	)			const

Get a pointer to the first hit attached in a particular gap.

Definition at line 1615 of file RecMucTrack.cxx.

References MucID::getGapNum(), and m_pHits.

{
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "RecMucTrack::Hit-E1  invalid gap = " << gap << endl;
    return 0;
  }
  
  vector<MucRecHit*>::const_iterator iHit;
  
  for (iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      if ( (*iHit)->Part() == part &&
           (*iHit)->Gap()  == gap) { 
        return (*iHit);
      }
    }
  }
  
  return 0L;
}

float RecMucTrack::GetHitDistance

(

const MucRecHit *

hit

)

Calculate the distance of the hit to the intersection in read direction.

Definition at line 483 of file RecMucTrack.cxx.

References DstMucTrack::distance(), MucRecHit::Gap(), MucRecHit::GetCenterPos(), MucRecHit::GetGap(), MucID::getGapNum(), m_CurrentInsct, MucGeoGap::Orient(), MucRecHit::Part(), and MucGeoGap::TransformToGap().

Referenced by MucRecRoadFinder::TrackFinding().

{
  if (!hit) {
    cout << "RecMucTrack:GetHitDistance-E1  null hit pointer!" << endl;
    return kInfinity;
  }

  int part = hit->Part();
  int gap  = hit->Gap();
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "RecMucTrack::GetHitDistance()  bad gap number = " << gap << endl;
    return kInfinity;
  }

  HepPoint3D posHit        = hit->GetCenterPos();
  HepPoint3D posHitLocal   = hit->GetGap()->TransformToGap(posHit);
  HepPoint3D posInsctLocal = hit->GetGap()->TransformToGap(m_CurrentInsct);
  int orient = hit->GetGap()->Orient();

  float distance = -9990;
  if(orient == 1) distance = fabs(posInsctLocal.x() - posHitLocal.x());
  if(orient == 0) distance = fabs(posInsctLocal.y() - posHitLocal.y());

  return distance;
}

float RecMucTrack::GetHitDistance2

(

const MucRecHit *

hit

)

no abs value

Definition at line 511 of file RecMucTrack.cxx.

References DstMucTrack::distance(), MucRecHit::Gap(), MucRecHit::GetCenterPos(), MucRecHit::GetGap(), MucID::getGapNum(), m_CurrentInsct, MucGeoGap::Orient(), MucRecHit::Part(), and MucGeoGap::TransformToGap().

Referenced by MucRecTrkExt::execute(), LineFit(), and MucRecTrkExt::TrackFinding().

{
  if (!hit) {
    cout << "RecMucTrack:GetHitDistance-E1  null hit pointer!" << endl;
    return kInfinity;
  } 

  int part = hit->Part();
  int gap  = hit->Gap();
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "RecMucTrack::GetHitDistance()  bad gap number = " << gap << endl;
    return kInfinity;
  }
  
  HepPoint3D posHit        = hit->GetCenterPos();
  HepPoint3D posHitLocal   = hit->GetGap()->TransformToGap(posHit);
  HepPoint3D posInsctLocal = hit->GetGap()->TransformToGap(m_CurrentInsct);
  int orient = hit->GetGap()->Orient();
      
  float distance = -9990;
  if(orient == 1) distance = (posInsctLocal.x() - posHitLocal.x());
  if(orient == 0) distance = (posInsctLocal.y() - posHitLocal.y());
 
  //cout<<"========in RecMucTrack: line insct: "<<posInsctLocal.x()<<" "<<posInsctLocal.y()<<"   ->  "<<posHitLocal.x()<<" "<<posHitLocal.y()<<endl;
 
  return distance;
} 

vector< long > RecMucTrack::GetHitIndices

(

)

const

Get indices of all hits in the track.

Definition at line 1649 of file RecMucTrack.cxx.

References m_pHits, and v.

{
  vector<long> v;
  
  vector<MucRecHit*>::const_iterator iHit;
  for ( iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      long index = (*iHit)->GetID();
      v.push_back(index);
      /*
      cout << " Muc2DRoad::HitIndices  gap  orientation  twopack= "
           <<  (*iHit)->ChannelID().Plane() << "   "
           <<  (*iHit)->ChannelID().Orient()  << "   "
           <<  (*iHit)->ChannelID().TwoPack() << endl; 
      */
    }
  }
  return v;
}

int RecMucTrack::GetHitInGap	(	const int	part,
		const int	gap
	)			const

How many hits per gap does this track contain?

Definition at line 1443 of file RecMucTrack.cxx.

References MucID::getGapNum(), and m_pHits.

Referenced by ComputeMaxHitsInGap().

{
  if ( part < 0 || part > 2 ) {
    cout << "RecMucTrack::GetHitInGap(), invalid part " << part << endl;
    return 0;
  }
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "RecMucTrack::GetHitInGap(), invalid gap " << gap << endl;
    return 0;
  } 
  
  vector<MucRecHit*>::const_iterator iHit;
  int hitsInGap = 0;
  
  for (iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    
    if ( !(*iHit) ) {
      cout << "RecMucTrack::GetHitInGap() null hit pointer !"  << endl;
      return 0;
    }
    else {
      if( part == (*iHit)->Part() && gap  == (*iHit)->Gap() ) hitsInGap++;
    }
  }
  
  return hitsInGap;
}

int RecMucTrack::GetHitInSeg	(	const int	part,
		const int	seg
	)			const

How many hits does a segment contains.

Definition at line 1472 of file RecMucTrack.cxx.

References GetHitInSegOrient(), and num.

Referenced by FindSegWithMaxHits().

{
  int num = GetHitInSegOrient(part, seg, 0) + GetHitInSegOrient(part, seg, 1);
  return num;
}

int RecMucTrack::GetHitInSegOrient	(	const int	part,
		const int	seg,
		const int	orient
	)			const

How many hits does a segment contains in one orient.

Definition at line 1479 of file RecMucTrack.cxx.

References MucID::getSegNum(), and m_pHits.

Referenced by GetHitInSeg().

{
  if ( part < 0 || part > 2 ) {
    cout << "RecMucTrack::GetHitInSeg(), invalid part " << part << endl;
    return 0;
  }
  if ( (seg < 0) || (seg >= (int)MucID::getSegNum(part)) ) {
    cout << "RecMucTrack::GetHitInSeg(), invalid seg = " << seg << endl;
    return 0;
  } 
  
  vector<MucRecHit*>::const_iterator iHit;
  int hitsInSeg = 0;
  
  for (iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    
    if ( !(*iHit) ) {
      cout << "RecMucTrack::GetHitInSeg(), null hit pointer !"  << endl;
      return 0;
    }
    else {
      if( part == (*iHit)->Part() && seg == (*iHit)->Seg() && 
          orient == (*iHit)->GetGap()->Orient() ) hitsInSeg++;
    }
  }
  
  return hitsInSeg;
}

vector< MucRecHit * > RecMucTrack::GetHits

(

)

const

Get all hits on this track.

Definition at line 1637 of file RecMucTrack.cxx.

References m_pHits.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  return m_pHits;
}

void RecMucTrack::GetMdcExtTrack	(	Hep3Vector	mdcStartPos,
		Hep3Vector	mdcStartMomentum,
		int	charge,
		Hep3Vector &	mucStartPos,
		Hep3Vector &	mucStartMomentum
	)

compute ext track myself from mdc.

Definition at line 337 of file RecMucTrack.cxx.

References abs, cos(), DstMucTrack::deltaPhi(), IsInsideSuperConductor(), iter(), kDeltaPhi, kMagnetField, kMdcExtIterMax, kvC, and sin().

Referenced by MucRecTrkExt::execute().

{

  //cm->mm; GeV->MeV//
  mdcStartPos*=10;
  mdcStartMomentum*=1000;


  Hep3Vector pt = mdcStartMomentum;
  pt.setZ(0);
  //cout << "pt " << pt.mag() << endl;
  double radius = (pt.mag() * 1e6 / kvC * 1e3) / (fabs(charge * kMagnetField)) ;
  //double startPhi = startP.phi();
  double deltaPhi = -1.0 * (charge / abs(charge)) * kDeltaPhi;
  double deltaZ = (mdcStartMomentum.z() * 1e6 / kvC * 1e3) * kDeltaPhi / (abs(charge) * kMagnetField);

  //cout << "r = " << radius << endl;
  mucStartPos = mdcStartPos;
  mucStartMomentum = mdcStartMomentum;
  double phi;
  int    iter = 0;
  do {
    phi = mucStartMomentum.getPhi() + deltaPhi;
    mucStartPos.setX(mucStartPos.x() + radius * kDeltaPhi * cos(phi));
    mucStartPos.setY(mucStartPos.y() + radius * kDeltaPhi * sin(phi));
    mucStartPos.setZ(mucStartPos.z() + deltaZ);
    
    mucStartMomentum.setPhi(mucStartMomentum.phi() + deltaPhi);
    iter++;
    //cout << endP << "  " << mucStartPos << endl;
  }
  while(IsInsideSuperConductor(mucStartPos) && iter < kMdcExtIterMax);

  //mm->cm; MeV->GeV//
  mucStartPos/=10;
  mucStartMomentum/=1000;
}

Hep3Vector RecMucTrack::getMdcMomentum

(

)

const [inline]

momentum of this track in Mdc

Definition at line 206 of file RecMucTrack.h.

References m_MdcMomentum.

Referenced by MucRecTrkExt::execute(), and MucRecTrkExt::TrackFinding().

{ return m_MdcMomentum/1000; }

Hep3Vector RecMucTrack::getMucMomentum

(

)

const [inline]

Start momentum of this track in Muc.

Definition at line 203 of file RecMucTrack.h.

References m_MucMomentum.

Referenced by MucRecRoadFinder::execute().

{ return m_MucMomentum/1000; }

Hep3Vector RecMucTrack::getMucPos

(

)

const [inline]

start position of this track in Muc.

Definition at line 198 of file RecMucTrack.h.

References m_MucPos.

Referenced by MucRecRoadFinder::execute().

{ return m_MucPos/10; }

Hep3Vector RecMucTrack::getMucPosSigma

(

)

const [inline]

Definition at line 200 of file RecMucTrack.h.

References m_MucPosSigma.

{ return m_MucPosSigma/10; }

Hep3Vector RecMucTrack::GetMucStripPos

(

)

const [inline]

Definition at line 261 of file RecMucTrack.h.

References m_StripPhi.

Referenced by MucRecTrkExt::execute().

{ return m_StripPhi;}

int RecMucTrack::GetNGapsWithHits

(

)

const [inline]

How many gaps provide hits attached to this track?

Definition at line 267 of file RecMucTrack.h.

References DstMucTrack::m_numLayers.

Referenced by LineFit().

{ return m_numLayers; }

int RecMucTrack::GetNSharedHits

(

const RecMucTrack *

track

)

const

How many hits do two tracks share?

Definition at line 1586 of file RecMucTrack.cxx.

References count, MucRecHit::GetID(), and m_pHits.

{
  if (!track2) {
    return 0;
  }

  int count = 0;
  vector<MucRecHit*>::const_iterator iHit1;
  vector<MucRecHit*>::const_iterator iHit2;
  MucRecHit *hit1, *hit2;

  for( iHit1 = m_pHits.begin(); iHit1 != m_pHits.end(); iHit1++){
    for( iHit2 = track2->m_pHits.begin(); 
         iHit2 != track2->m_pHits.end(); iHit2++){
      hit1 = (*iHit1);
      hit2 = (*iHit2);
      
      if ( (hit1 != 0) && (hit2 != 0) ) {
        if (hit1->GetID() == hit2->GetID()) {
          count++;
        }
      }
    }
  }
  
  return count;
}

vector<float> RecMucTrack::getQuadDistHits

(

)

const [inline]

Definition at line 313 of file RecMucTrack.h.

References m_distHits_quad.

Referenced by MucRecTrkExt::execute().

{ return  m_distHits_quad ;}

int RecMucTrack::GetRecMode

(

)

const [inline]

Definition at line 256 of file RecMucTrack.h.

References m_recmode.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{return m_recmode;}

int RecMucTrack::GetTotalHits

(

)

const

How many hits in all does this track contain?

Definition at line 1437 of file RecMucTrack.cxx.

References m_pHits.

Referenced by ComputeDepth(), MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{ 
  return m_pHits.size();
}

vector<int> RecMucTrack::getVecHits

(

void

)

const [inline]

Definition at line 305 of file RecMucTrack.h.

References m_vecHits.

{ return  m_vecHits ;}

bool RecMucTrack::HasHit	(	const int	part,
		const int	seg,
		const int	gap,
		const int	strip
	)			const

Does this track contains an assigned hit?

Definition at line 1543 of file RecMucTrack.cxx.

References m_pHits.

{
  bool found = false;
  vector<MucRecHit*>::const_iterator iHit;

  for (iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      if ( (*iHit)->Part()  == part &&
           (*iHit)->Seg()   == seg &&
           (*iHit)->Gap()   == gap &&
           (*iHit)->Strip() == strip ) {
        found = true;
      }
    }
  }

  return found;
}

bool RecMucTrack::HasHitInGap	(	const int	part,
		const int	gap
	)			const

Does this track contain any hits in the given gap?

Definition at line 1563 of file RecMucTrack.cxx.

References MucID::getGapNum(), and m_pHits.

Referenced by ComputeDepth().

{
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "RecMucTrack::HasHitInGap-E2  invalid gap = " << gap << endl;
    return false;
  }
  
  bool found = false;
  vector<MucRecHit*>::const_iterator iHit;

  for (iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      if ( (*iHit)->Part() == part &&
           (*iHit)->Gap() == gap ) {
        found = true;
      }
    }
  }

  return found;
}

int DstMucTrack::id

(

)

const [inline, inherited]

Definition at line 33 of file DstMucTrack.h.

References DstMucTrack::m_id.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::TrackFinding().

{ return   m_id; }

bool RecMucTrack::IsInsideSuperConductor

(

Hep3Vector

pos

)

Definition at line 379 of file RecMucTrack.cxx.

References kSuperConductorR, and kSuperConductorZ.

Referenced by GetMdcExtTrack().

{
  if(pos.mag() < kSuperConductorR && fabs(pos.z()) < kSuperConductorZ) {
    return true;
  }
  else {
    return false;
  }
}

int DstMucTrack::kalbrLastLayer

(

)

const [inline, inherited]

Definition at line 71 of file DstMucTrack.h.

References DstMucTrack::m_kalbrLastLayer.

{return m_kalbrLastLayer;}

double DstMucTrack::kaldepth

(

)

const [inline, inherited]

Definition at line 70 of file DstMucTrack.h.

References DstMucTrack::m_kaldepth.

{return m_kaldepth;}

int DstMucTrack::kaldof

(

)

const [inline, inherited]

Definition at line 69 of file DstMucTrack.h.

References DstMucTrack::m_kaldof.

{return m_kaldof;}

int DstMucTrack::kalecLastLayer

(

)

const [inline, inherited]

Definition at line 72 of file DstMucTrack.h.

References DstMucTrack::m_kalecLastLayer.

{return m_kalecLastLayer;}

double DstMucTrack::kalRechi2

(

)

const [inline, inherited]

Definition at line 68 of file DstMucTrack.h.

References DstMucTrack::m_kalrechi2.

{return m_kalrechi2;}

void RecMucTrack::LineFit

(

int

fittingMethod

)

Line fit with hits on a seg with max hits.

Definition at line 1038 of file RecMucTrack.cxx.

References MucRec2DRoad::AttachHit(), ComputeLastGap(), cos(), DstMucTrack::deltaPhi(), check_raw_filter::dist, Extend(), FindSegWithMaxHits(), MucID::gap(), MucRec3DRoad::GetDegreesOfFreedom(), MucID::getGapNum(), GetHitDistance2(), MucRec2DRoad::GetIntercept(), GetNGapsWithHits(), MucRec3DRoad::GetPart(), MucRec3DRoad::GetReducedChiSquare(), MucRec2DRoad::GetSlope(), MucRec3DRoad::GetSlopeZX(), MucRec3DRoad::GetSlopeZY(), genRecEmupikp::i, kMinor, kPi, DstMucTrack::m_chi2, DstMucTrack::m_deltaPhi, DstMucTrack::m_depth, m_distHits, m_distHits_quad, DstMucTrack::m_dof, DstMucTrack::m_ecLastLayer, DstMucTrack::m_endPart, m_Good3DLine, m_Good3DPart, m_IntersectionInner, m_IntersectionOuter, m_MdcMomentum, m_MucMomentum, m_MucPos, m_MucPosSigma, m_pExpectedHits, m_pHits, DstMucTrack::m_px, DstMucTrack::m_py, DstMucTrack::m_pz, DstMucTrack::m_rms, DstMucTrack::m_status, DstMucTrack::m_xPos, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPos, DstMucTrack::m_yPosSigma, DstMucTrack::m_zPos, DstMucTrack::m_zPosSigma, momentum, MucID::part(), MucRec3DRoad::Project(), MucRec3DRoad::ProjectNoCurrentGap(), MucRec3DRoad::ProjectToStrip(), MucRec3DRoad::ProjectWithSigma(), DstMucTrack::px(), DstMucTrack::py(), DstMucTrack::pz(), MucRec3DRoad::RefitNoCurrentGap(), MucID::seg(), SetCurrentInsct(), MucRec3DRoad::SetSimpleFitParams(), sin(), DstMucTrack::status(), MucID::strip(), and MucRec3DRoad::TransformPhiRToXY().

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  Extend();

  int part = -1, seg = -1;
  int nHits = FindSegWithMaxHits(part, seg);
//   cout << "RecMucTrack::ComputeDepth(), part " << part << " seg " << seg 
//     << " contains most hits " << nHits << endl;
  if (part < 0 || seg < 0) {
    m_depth = 0;
    //cout << "No hit in this track" << endl;
    return;
  }

  float startPos[3] = {0.0, 0.0, 0.0};
  MucRec2DRoad road2D[2];
  vector<MucRecHit*>::const_iterator iHit;
  for (int orient = 0; orient <= 1; orient++) {
    road2D[orient] = MucRec2DRoad(part, seg, orient, startPos[0], startPos[1], startPos[2],fittingMethod );
    for(iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
      if(*iHit) {  // Check for a null pointer.
        int hitPart   = (*iHit)->Part();
        int hitSeg    = (*iHit)->Seg();
        int hitOrient = (*iHit)->GetGap()->Orient();
        if (hitPart == part && hitSeg == seg && hitOrient == orient) {
          road2D[orient].AttachHit(*iHit);
        }
      }
    }
    //cout << "Orient " << orient
    // << " Slope " << road2D[orient].GetSlope()
    // << " Intercept " << road2D[orient].GetIntercept()
    // << " LastGap " << road2D[orient].GetLastGap()
    // << endl;

  }
  MucRec3DRoad road3D(&road2D[0], &road2D[1]);
  float vx, vy, x0, y0;
  if ( road3D.GetPart() == 1 ){
    road3D.TransformPhiRToXY( road2D[1].GetSlope(),     road2D[0].GetSlope(),     
        road2D[1].GetIntercept(), road2D[0].GetIntercept(),
        vx, vy, x0, y0);
  }
  else {
    vx = road2D[1].GetSlope();
    x0 = road2D[1].GetIntercept();
    vy = road2D[0].GetSlope();
    y0 = road2D[0].GetIntercept();
  }
  
  //cout << "road3D Last Gap " << road3D.GetLastGap() << endl;
  //cout << "vx " << vx << " x0 " << x0 << " vy " << vy << " y0 " << y0 << endl;

  // if number of gaps with hits >= 2 on both orient, change muc pos and momentum
  if (road2D[0].GetNGapsWithHits() >= 2 && road2D[1].GetNGapsWithHits() >= 2) {

    //good 3d road!!!
    m_Good3DLine = 1;
    m_status     = 1;

    m_Good3DPart = road3D.GetPart();
    road3D.SetSimpleFitParams(vx, x0, vy, y0);
    //road3D.PrintHitsInfo();
    
    float startx = 0.0, starty = 0.0, startz = 0.0;
    float startxSigma = 0.0, startySigma = 0.0, startzSigma = 0.0;
    float x1 = 0.0, y1 = 0.0, z1 = 0.0, x2 = 0.0, y2 = 0.0, z2 = 0.0;
    int gap = 0;
    
    if(fittingMethod==2){    //if choose quadratic fitting method!
      for(int igap=0; igap<9;igap++){
        road3D.Project(igap, x1, y1, z1, x2, y2, z2);
        m_IntersectionInner[igap].set(x1,y1,z1);
        m_IntersectionOuter[igap].set(x2,y2,z2);
        //cout<<"3dproject sur "<<x1<<" "<<y1<<" "<<z1<<" "<<x2<<" "<<y2<<" "<<z2<<endl;      
      }
    }

    do {
      //road3D.Project(gap, startx, starty, startz);
      road3D.ProjectWithSigma(gap, startx, starty, startz, startxSigma, startySigma, startzSigma);
      gap++;
    }
    while ( (startx*startx + starty*starty + startz*startz) < kMinor &&
              gap < (int)MucID::getGapNum(part) );
    
    if(fabs(startx)<2600&&fabs(starty)<2600&&fabs(startz)<2800){
      Hep3Vector MucPos_self;
      MucPos_self.set(startx, starty, startz);
      float dist = (MucPos_self - m_MucPos).mag();
    
      if(dist < 1000 ){  // (mm) maybe the fit is bad, if the dist is too big.
        m_MucPos.set(startx, starty, startz);
        m_xPos = startx;
        m_yPos = starty;
        m_zPos = startz;
      }
    }
    m_MucPosSigma.set(startxSigma, startySigma, startzSigma);
    m_xPosSigma = startxSigma;
    m_yPosSigma = startySigma;
    m_zPosSigma = startzSigma;
    
    //cout<<"in RecMucTrack gap= "<<gap<<" start= "<<startx<<" "<<starty<<" "<<startz<<" "<<endl;
    float momentum = m_MucMomentum.mag();
    float zDir = 1.0;
    if (m_MucMomentum.z() != 0.0) zDir  = m_MucMomentum.z() / fabs(m_MucMomentum.z());
    
    float px = road3D.GetSlopeZX()*zDir;
    float py = road3D.GetSlopeZY()*zDir;
    float pz = zDir;
    float segPhi = 0.25*kPi*seg;
    // if vr is opposite to original, but both are very small, e.x. -0.01 -> + 0.01, or you can say, pt~p, change it to guarantee vx, vy is correct.
    if (part == 1 && px*cos(segPhi)+py*sin(segPhi) < 0.0) { // when in barrel, pt dot segPhi < 0, angle between them > 90deg
      px *= -1.0;
      py *= -1.0;
      pz *= -1.0;
    }

    //if(sqrt(px*px+py*py)>0.01)m_MucMomentum.set(px, py, pz);
    m_MucMomentum.set(px, py, pz);
    m_MucMomentum.setMag(momentum);

    m_px = px; m_py = py; m_pz = pz;
    Hep3Vector phi_mdc, phi_muc;
    phi_mdc.set(m_MdcMomentum.x(),m_MdcMomentum.y(),0);
    phi_muc.set(px,py,0);
    double deltaPhi = phi_mdc.angle(phi_muc);

    m_deltaPhi = deltaPhi;
    m_dof  = road3D.GetDegreesOfFreedom();
    m_chi2 = road3D.GetReducedChiSquare();
    if(m_chi2<0||m_chi2>1000)  m_chi2 = 0;
    m_rms  = 0.0; // road3D.GetReducedChiSquare(); // ??? in MucRecRoad3D


   if(road2D[0].GetNGapsWithHits() >= 3 && road2D[1].GetNGapsWithHits() >= 3){ //need 3+ gap to do refit
      //calc distance between each hit with this track.
      //use GetHitDistance(), so we should provide m_CurrentInsct. 
      float xInsct = 0.0, yInsct = 0.0, zInsct = 0.0, sigmax = 0.0, sigmay = 0.0, sigmaz = 0.0;
      float quad_x1 = 0.0, quad_y1 = 0.0, quad_z1 = 0.0, quad_x2 = 0.0, quad_y2 = 0.0, quad_z2 = 0.0;
      for(int ihit = 0; ihit < m_pHits.size(); ihit++){
        int igap = (m_pHits[ihit])->Gap();
        road3D.ProjectNoCurrentGap(igap, xInsct, yInsct, zInsct, sigmax, sigmay, sigmaz, 
                                   quad_x1, quad_y1, quad_z1, quad_x2, quad_y2, quad_z2);

        SetCurrentInsct(xInsct, yInsct, zInsct);

        float dist_hit_track = GetHitDistance2(m_pHits[ihit]);
        m_distHits.push_back(dist_hit_track);

        //cout<<"===========in RecMucTrack: line dist: "<<dist_hit_track<<endl;

        if(fittingMethod==2) // ------calc local insct in a gap with quad line.
          {
           
            Hep3Vector center = m_pHits[ihit]->GetCenterPos();
            int iPart = m_pHits[ihit]->Part();
            int iSeg  = m_pHits[ihit]->Seg();
            int iGap  = m_pHits[ihit]->Gap();
            float xHit, yHit, zHit = 0.;
            if (iPart == 1) {
                    if ( iGap %2  == 1) {
                            xHit = center.z();
                            yHit = sqrt(center.x()*center.x()
                                            + center.y()*center.y());
                            if(iSeg==2) yHit = center.y();  //deal with seg2 seperately! because there is a hole here. 2006.11.9
                    } 
                    else {
                            xHit = center.x();
                            yHit = center.y();
                    }
            }
            else {
                    if ( iGap %2 == 0) {
                            xHit = center.z();
                            yHit = center.y();
                    }
                    else {
                            xHit = center.z();
                            yHit = center.x();
                    }
            }
            
            float distance1 = fabs(xHit - quad_x1)/(xHit - quad_x1) * sqrt((xHit - quad_x1)*(xHit - quad_x1) + (yHit - quad_y1)*(yHit - quad_y1));
            float distance2 = fabs(xHit - quad_x2)/(xHit - quad_x2) * sqrt((xHit - quad_x2)*(xHit - quad_x2) + (yHit - quad_y2)*(yHit - quad_y2));

            float dist_quad = distance1;
            if(fabs(distance1) > fabs(distance2)) dist_quad = distance2;
            
      //cout<<"============in RecMucTrack: quad xyz: "<<quad_x1<<" "<<quad_y1<<" "<<quad_z1<<" "<<quad_x2<<" "<<quad_y2<<endl;
            //cout<<"============in RecMucTrack: hit pos: "<<xHit<<" "<<yHit<<" "<<zHit<<endl;
            //cout<<"============in RecMucTrack: dist1 = "<<distance1<<" dist2 = "<<distance2<<" dist ="<<dist_quad<<endl; 

            if(quad_x1 == -9999) m_distHits_quad.push_back(-99);
            else m_distHits_quad.push_back(dist_quad);

          }
      }
    }
    else{
        for(int ihit = 0; ihit < m_pHits.size(); ihit++){
            m_distHits.push_back(-99);
            m_distHits_quad.push_back(-99);
        }
    }

    ComputeLastGap();
    //HepPoint3D initPos(startx, starty,startz);
    HepPoint3D initPos(startx-m_MucMomentum.x()/momentum, starty-m_MucMomentum.y()/momentum,startz-m_MucMomentum.z()/momentum);

    //for(int igap = 0; igap <= m_brLastLayer; igap++){
    for(int igap = 0; igap < 9 ; igap++){
      vector<int> padID;
      vector<float> intersect_x;
      vector<float> intersect_y;
      vector<float> intersect_z;

      //refit---------------------------------
      float zx, zy, x0, y0;
      int status = road3D.RefitNoCurrentGap(igap,zx,x0,zy,y0);
      Hep3Vector mom_refit;
      if(status == 0){ //refit succeed
              float zDir = 1.0;
              if (m_MucMomentum.z() != 0.0) zDir  = m_MucMomentum.z() / fabs(m_MucMomentum.z());
              float px_refit = zx * zDir;  
              float py_refit = zy * zDir;       
              float pz_refit = zDir;
              float segPhi = 0.25*kPi*seg;
              // if vr is opposite to original, but both are very small, e.x. -0.01 -> + 0.01, or you can say, pt~p, change it to guarantee vx, vy is correct.
              if (part == 1 && px_refit*cos(segPhi)+py_refit*sin(segPhi) < 0.0) { // when in barrel, pt dot segPhi < 0, angle between them > 90deg
                      px_refit *= -1.0;
                      py_refit *= -1.0;
                      pz_refit *= -1.0;
              }

              mom_refit.setX(px_refit);
              mom_refit.setY(py_refit);
              mom_refit.setZ(pz_refit);
              mom_refit.setMag(momentum);
      }
      else mom_refit = m_MucMomentum;  //use the former momentum

      HepPoint3D initPos_refit;
      if(status == 0){
              float initPosx_refit, initPosy_refit, initPosz_refit;
              float sigmax_refit, sigmay_refit, sigmaz_refit;
              road3D.Project(0, zx, x0, zy, y0, initPosx_refit, initPosy_refit, initPosz_refit);
              initPos_refit.setX(initPosx_refit-mom_refit.x()/momentum);
              initPos_refit.setY(initPosy_refit-mom_refit.y()/momentum);
              initPos_refit.setZ(initPosz_refit-mom_refit.z()/momentum);
      }
      else initPos_refit = initPos;

      //cout<<"initPos: "<<initPos<<"  "<<initPos_refit<<endl;
      //if((initPos - initPos_refit).mag()>100) cout<<"--------=====------to far"<<endl;
      //refit---------------------------------
      //cout<<"no gap: "<<igap<<"  mom: "<<m_MucMomentum<<"  refit: "<<mom_refit<<endl; 
      
      //vector<Identifier> MuId = road3D.ProjectToStrip(1,igap,initPos,m_MucMomentum,padID,intersect_x,intersect_y,intersect_z);
      vector<Identifier> MuId = road3D.ProjectToStrip(1,igap,initPos_refit,mom_refit,padID,intersect_x,intersect_y,intersect_z);

      //vector<Identifier> MuId = road3D.ProjectToStrip(1,igap,initPos,mom_refit,padID,intersect_x,intersect_y,intersect_z);

      vector<Identifier>::const_iterator mucid;
      int i = 0;
      for(mucid = MuId.begin(); mucid != MuId.end(); mucid++,i++)
      {
        MucRecHit *pHit = new MucRecHit(MucID::part(*mucid),MucID::seg(*mucid),MucID::gap(*mucid),MucID::strip(*mucid));
        pHit->SetPadID(padID[i]);
        pHit->SetIntersectX(intersect_x[i]);
        pHit->SetIntersectY(intersect_y[i]);
        pHit->SetIntersectZ(intersect_z[i]);

        m_pExpectedHits.push_back(pHit);
      }
    }

    if(m_endPart!=-1&&m_ecLastLayer!=-1)
    {
      //for(int igap = 0; igap <= m_ecLastLayer; igap++){
      for(int igap = 0; igap < 8; igap++){
        vector<int> padID;
        vector<float> intersect_x;
        vector<float> intersect_y;
        vector<float> intersect_z;

        //refit---------------------------------
        float zx, zy, x0, y0;
        //int status = road3D.RefitNoCurrentGap(igap,zy,y0,zx,x0);  //!!!!!different sequence
        int status = road3D.RefitNoCurrentGap(igap,zx,x0,zy,y0); 
        Hep3Vector mom_refit;
        if(status == 0){ //refit succeed
                float zDir = 1.0;
                if (m_MucMomentum.z() != 0.0) zDir  = m_MucMomentum.z() / fabs(m_MucMomentum.z());
                float px_refit = zx * zDir;
                float py_refit = zy * zDir;
                float pz_refit = zDir;
                float segPhi = 0.25*kPi*seg;
                // if vr is opposite to original, but both are very small, e.x. -0.01 -> + 0.01, or you can say, pt~p, change it to guarantee vx, vy is correct.
                if (part == 1 && px_refit*cos(segPhi)+py_refit*sin(segPhi) < 0.0) { // when in barrel, pt dot segPhi < 0, angle between them > 90deg
                        px_refit *= -1.0;
                        py_refit *= -1.0;
                        pz_refit *= -1.0;
                }

                mom_refit.setX(px_refit);
                mom_refit.setY(py_refit);
                mom_refit.setZ(pz_refit);
                mom_refit.setMag(momentum);
        }
        else mom_refit = m_MucMomentum;  //use the former momentum

        HepPoint3D initPos_refit;
        if(status == 0){
                float initPosx_refit, initPosy_refit, initPosz_refit;
                float sigmax_refit, sigmay_refit, sigmaz_refit;
                road3D.Project(0, zx, x0, zy, y0, initPosx_refit, initPosy_refit, initPosz_refit);
                initPos_refit.setX(initPosx_refit-mom_refit.x()/momentum);
                initPos_refit.setY(initPosy_refit-mom_refit.y()/momentum);
                initPos_refit.setZ(initPosz_refit-mom_refit.z()/momentum);
        }
        else initPos_refit = initPos;
        //refit---------------------------------
        //cout<<"mom: "<<m_MucMomentum<<" "<<mom_refit<<" pos: "<<initPos<<"  "<<initPos_refit<<endl;  
        //vector<Identifier> MuId = road3D.ProjectToStrip(m_endPart,igap,initPos,m_MucMomentum,padID,intersect_x,intersect_y,intersect_z);

        vector<Identifier> MuId = road3D.ProjectToStrip(m_endPart,igap,initPos_refit,mom_refit,padID,intersect_x,intersect_y,intersect_z);

        vector<Identifier>::const_iterator mucid;
        int i = 0;
        for(mucid = MuId.begin(); mucid != MuId.end(); mucid++,i++)
        {
          MucRecHit *pHit = new MucRecHit(MucID::part(*mucid),MucID::seg(*mucid),MucID::gap(*mucid),MucID::strip(*mucid));
          pHit->SetPadID(padID[i]);
          pHit->SetIntersectX(intersect_x[i]);
          pHit->SetIntersectY(intersect_y[i]);
          pHit->SetIntersectZ(intersect_z[i]);

          m_pExpectedHits.push_back(pHit);
        }
      }
    }

    //cout<<"in RecMucTrack push back expected hits  size= "<<m_pExpectedHits.size()<<" ? "<<m_pHits.size()<<endl;
    for(int i=0; i< m_pExpectedHits.size(); i++)
    {
      MucRecHit *ihit = m_pExpectedHits[i];
      //cout<<"expected Hits: "<<ihit->Part()<<"  "<<ihit->Seg()<<"  "<<ihit->Gap()<<"  "<<ihit->Strip()<<endl;
      //cout<<"pad: "<<ihit->GetPadID()<<"  pos: "<<ihit->GetIntersectX()<<" "<<ihit->GetIntersectY()<<" "<<ihit->GetIntersectZ()<<" "<<endl;
    }

    for(int i=0; i< m_pHits.size(); i++)
    {
      MucRecHit *ihit = m_pHits[i];
      //cout<<"attachedd Hits: "<<ihit->Part()<<"  "<<ihit->Seg()<<"  "<<ihit->Gap()<<"  "<<ihit->Strip()<<" isseed: "<<ihit->HitIsSeed()<<endl;
    }
    //cout<<"1st: "<<brFirstLayer()<<"  "<<ecFirstLayer()<<"  last: "<<brLastLayer()<<"  "<<ecLastLayer()<<endl;
  
  }
  else {
    //calc distance between each hit with this track.
    //initialize distHits
    for(int ihit = 0; ihit < m_pHits.size(); ihit++){
      m_distHits.push_back(-99);
      m_distHits_quad.push_back(-99);
    }
  }
  //cout << "Muc Pos: x " << m_MucPos.x() << " y " << m_MucPos.y() << " z " << m_MucPos.z() << endl;
}

int DstMucTrack::maxHitsInLayer

(

)

const [inline, inherited]

Definition at line 43 of file DstMucTrack.h.

References DstMucTrack::m_maxHitsInLayer.

Referenced by MucRecTrkExt::execute(), MucRecRoadFinder::execute(), and MucPID::particleIDCalculation().

{ return   m_maxHitsInLayer; }

int DstMucTrack::numHits

(

)

const [inline, inherited]

Definition at line 41 of file DstMucTrack.h.

References DstMucTrack::m_numHits.

Referenced by MucRecRoadFinder::execute(), and DQARhopi::execute().

{ return   m_numHits;      } // total hits on the track   

int DstMucTrack::numLayers

(

)

const [inline, inherited]

Definition at line 42 of file DstMucTrack.h.

References DstMucTrack::m_numLayers.

Referenced by MucRecRoadFinder::execute(), DQARhopi::execute(), and Ppjrhopi::execute().

{ return   m_numLayers;    } // number of layers with hits 

RecMucTrack & RecMucTrack::operator=

(

const DstMucTrack &

dstTrack

)

Assignment constructor from DstMucTrack.

Definition at line 186 of file RecMucTrack.cxx.

References SetDefault().

{
  SetDefault();
  DstMucTrack::operator=(dstTrack);
  return *this;
}

RecMucTrack & RecMucTrack::operator=

(

const RecMucTrack &

orig

)

Assignment constructor.

Definition at line 84 of file RecMucTrack.cxx.

References DstMucTrack::m_brLastLayer, m_changeUnit, DstMucTrack::m_chi2, m_CurrentDir, m_CurrentInsct, m_CurrentPos, DstMucTrack::m_deltaPhi, DstMucTrack::m_depth, m_Directions, DstMucTrack::m_dof, DstMucTrack::m_ecLastLayer, DstMucTrack::m_endPart, m_ExtTrackID, DstMucTrack::m_id, m_Intersections, DstMucTrack::m_kalbrLastLayer, DstMucTrack::m_kaldepth, DstMucTrack::m_kaldof, DstMucTrack::m_kalecLastLayer, DstMucTrack::m_kalrechi2, DstMucTrack::m_maxHitsInLayer, m_MdcMomentum, m_MdcPos, m_MucMomentum, m_MucPos, m_MucPosSigma, DstMucTrack::m_numHits, DstMucTrack::m_numLayers, m_pExpectedHits, m_pHits, m_recmode, DstMucTrack::m_rms, DstMucTrack::m_startPart, DstMucTrack::m_status, DstMucTrack::m_type, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPosSigma, and DstMucTrack::m_zPosSigma.

{
  // Assignment operator.
  if ( this != &orig ) {             // Watch out for self-assignment!
    //m_trackId         = orig.m_trackId;    //--------------->
    m_ExtTrackID    = orig.m_ExtTrackID;
    m_MdcPos        = orig.m_MdcPos;
    m_MdcMomentum   = orig.m_MdcMomentum;
    m_MucPos        = orig.m_MucPos;
    m_MucPosSigma   = orig.m_MucPosSigma;
    m_MucMomentum   = orig.m_MucMomentum;
    m_CurrentPos    = orig.m_CurrentPos;
    m_CurrentDir    = orig.m_CurrentDir;
    m_CurrentInsct  = orig.m_CurrentInsct;
    m_id            = orig.m_id;
    m_status        = orig.m_status;
    m_type          = orig.m_type;
    m_numHits       = orig.m_numHits;    //--------------->
    m_startPart     = orig.m_startPart;
    m_endPart       = orig.m_endPart;
    m_brLastLayer   = orig.m_brLastLayer;
    m_ecLastLayer   = orig.m_ecLastLayer;
    m_numLayers     = orig.m_numLayers;
    m_maxHitsInLayer= orig.m_maxHitsInLayer;
    m_depth         = orig.m_depth;
    m_dof           = orig.m_dof;
    m_chi2          = orig.m_chi2;
    m_rms           = orig.m_rms;
    m_deltaPhi      = orig.m_deltaPhi;
    m_pHits         = orig.m_pHits;
    m_pExpectedHits = orig.m_pExpectedHits;
    m_Intersections = orig.m_Intersections;
    m_Directions    = orig.m_Directions;
    m_xPosSigma     = orig.m_xPosSigma;
    m_yPosSigma     = orig.m_yPosSigma;
    m_zPosSigma     = orig.m_zPosSigma;
    m_changeUnit    = orig.m_changeUnit;
    m_recmode       = orig.m_recmode;
    //*******
    m_kalrechi2     = orig.m_kalrechi2;
    m_kaldof        = orig.m_kaldof;
    m_kaldepth      =orig.m_kaldepth;
    m_kalbrLastLayer=orig.m_kalbrLastLayer;
    m_kalecLastLayer=orig.m_kalecLastLayer;
  }

  return *this;
}

void RecMucTrack::OutputUnitChange

(

)

change unit

Definition at line 1706 of file RecMucTrack.cxx.

References m_changeUnit, DstMucTrack::m_depth, DstMucTrack::m_distance, DstMucTrack::m_px, DstMucTrack::m_py, DstMucTrack::m_pz, DstMucTrack::m_xPos, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPos, DstMucTrack::m_yPosSigma, DstMucTrack::m_zPos, and DstMucTrack::m_zPosSigma.

Referenced by ComputeTrackInfo().

{
  if(m_changeUnit == false){
    m_depth/=10;

    m_xPos /=10;
    m_yPos /=10;
    m_zPos /=10; 

    m_xPosSigma/=10;
    m_yPosSigma/=10;
    m_zPosSigma/=10;

    m_px /= 1000;
    m_py /= 1000;
    m_pz /= 1000;

    m_distance /= 10;
  }

  m_changeUnit = true;

}

void RecMucTrack::PrintHitsInfo

(

)

const

Print Hits Infomation.

Definition at line 1684 of file RecMucTrack.cxx.

References m_pHits.

{
  vector<MucRecHit*>::const_iterator iHit;
  for ( iHit = m_pHits.begin(); iHit != m_pHits.end(); iHit++) {
    if (*iHit) {  // Check for a null pointer.
      float xl, yl, zl;
      (*iHit)->GetStrip()->GetCenterPos(xl, yl, zl);
      HepPoint3D vl(xl, yl, zl);
      HepPoint3D vg = (*iHit)->GetGap()->TransformToGlobal(vl);
      
      cout << " part "   << (*iHit)->Part() 
           << " seg  "   << (*iHit)->Seg()
           << " gap  "   << (*iHit)->Gap()
           << " strip "  << (*iHit)->Strip()
           << " pos ("   << vg.x() << ", " << vg.y() << ", " << vg.z() << ")"  
           << endl;
    }
  }

}

void RecMucTrack::Project	(	const int &	part,
		const int &	gap,
		float &	x,
		float &	y,
		float &	z,
		int &	seg
	)

Where does the trajectory of this track intersect a specific gap?

Definition at line 422 of file RecMucTrack.cxx.

References MucGeoGeneral::FindIntersections(), MucID::getGapNum(), MucGeoGeneral::Instance(), intersection(), kMinor, kPi, m_CurrentDir, m_CurrentPos, and boss::pos.

Referenced by MucRecTrkExt::execute(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  seg = -1;
  x = 0.0; y = 0.0; z = 0.0;
  
  if ( (gap < 0) || (gap >= (int)MucID::getGapNum(part)) ) {
    cout << "Muc2DRoad::Project(), invalid gap = " << gap << endl;
    return;
  }

  HepPoint3D pos = m_CurrentPos;
  Hep3Vector dir = m_CurrentDir;

  vector<HepPoint3D> insctCon = MucGeoGeneral::Instance()->FindIntersections(part, gap, pos, dir);
  if( insctCon.size() == 0 ) {
    return;
  }

  HepPoint3D intersection = insctCon[0];
  
  x = intersection.x();
  y = intersection.y();
  z = intersection.z();

  //cout << " x " << x << " y " << y << " z " << z << endl;

  float phi;
  if( (x*x+y*y+z*z) < kMinor) {
    return;
  }
  else {
    if(part == 1) {
      phi = acos(x/sqrt(x*x+y*y));
      if(y < 0) phi = 2.0*kPi - phi;
      phi = fmod((phi + kPi/8.0), 2.0*kPi);
      seg = int(phi/(kPi/4.0)); 
    }
    else {
      if(x >= 0.0) {
        if(y >= 0.0) {
          seg = 0;
        } 
        else {
          seg = 3;
        }
      }
      else {
        if(y >= 0.0) {
          seg = 1;
        }
        else {
          seg = 2;
        }
      }
    }
  }
}

void RecMucTrack::pushExtDistHits

(

float

dist

)

[inline]

Definition at line 347 of file RecMucTrack.h.

References m_distHits_ext.

Referenced by MucRecTrkExt::execute().

{m_distHits_ext.push_back(dist);}

double DstMucTrack::px

(

)

const [inline, inherited]

Definition at line 58 of file DstMucTrack.h.

References DstMucTrack::m_px.

Referenced by LineFit().

{ return m_px; }   // momentum on the start of track in muc

double DstMucTrack::py

(

)

const [inline, inherited]

Definition at line 59 of file DstMucTrack.h.

References DstMucTrack::m_py.

Referenced by LineFit().

{ return m_py; }

double DstMucTrack::pz

(

)

const [inline, inherited]

Definition at line 60 of file DstMucTrack.h.

References DstMucTrack::m_pz.

Referenced by LineFit().

{ return m_pz; }

double DstMucTrack::rms

(

)

const [inline, inherited]

Definition at line 48 of file DstMucTrack.h.

References DstMucTrack::m_rms.

{ return   m_rms;    }

void DstMucTrack::setBrLastLayer

(

int

layer

)

[inline, inherited]

Definition at line 82 of file DstMucTrack.h.

References DstMucTrack::m_brLastLayer.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_brLastLayer = layer; }

void DstMucTrack::setChi2

(

double

chi2

)

[inline, inherited]

Definition at line 89 of file DstMucTrack.h.

References DstMucTrack::m_chi2.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_chi2   = chi2; }

void RecMucTrack::SetCurrentDir	(	const float	x,
		const float	y,
		const float	z
	)

set current direction of the trajectory.

Definition at line 292 of file RecMucTrack.cxx.

References m_CurrentDir.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_CurrentDir.set(x*1000, y*1000, z*1000);  //unnecessary, because it's dir, not mom
}

void RecMucTrack::SetCurrentInsct	(	const float	x,
		const float	y,
		const float	z
	)

set current intersection of the trajectory with strip plane.

Definition at line 300 of file RecMucTrack.cxx.

References m_CurrentInsct.

Referenced by MucRecTrkExt::execute(), LineFit(), MucRecTrkExt::TrackFinding(), and MucRecRoadFinder::TrackFinding().

{
  m_CurrentInsct.set(x, y, z);  
}

void RecMucTrack::SetCurrentPos	(	const float	x,
		const float	y,
		const float	z
	)

set current position of the trajectory.

Definition at line 284 of file RecMucTrack.cxx.

References m_CurrentPos.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_CurrentPos.set(x*10, y*10, z*10);
}

void RecMucTrack::SetDefault

(

)

Definition at line 193 of file RecMucTrack.cxx.

References m_brFirstLayer, m_ecFirstLayer, and m_Good3DPart.

Referenced by operator=(), and RecMucTrack().

{
  m_brFirstLayer = -99;
  m_ecFirstLayer = -99;
  m_Good3DPart   = -99;
}

void DstMucTrack::setDeltaPhi

(

double

dphi

)

[inline, inherited]

Definition at line 106 of file DstMucTrack.h.

References DstMucTrack::m_deltaPhi.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_deltaPhi = dphi; }

void DstMucTrack::setDepth

(

double

depth

)

[inline, inherited]

Definition at line 88 of file DstMucTrack.h.

References DstMucTrack::m_depth.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_depth  = depth; }

void DstMucTrack::setDistance

(

double

dist

)

[inline, inherited]

Definition at line 105 of file DstMucTrack.h.

References DstMucTrack::m_distance.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_distance = dist; }

void RecMucTrack::setDistHits

(

vector< float > &

disthits

)

[inline]

Definition at line 341 of file RecMucTrack.h.

References m_distHits.

Referenced by RecMucTrackCnv::TObjectToDataObject().

{ m_distHits = disthits;}

void DstMucTrack::setDof

(

int

dof

)

[inline, inherited]

Definition at line 90 of file DstMucTrack.h.

References DstMucTrack::m_dof.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_dof    = dof ; }

void DstMucTrack::setEcLastLayer

(

int

layer

)

[inline, inherited]

Definition at line 83 of file DstMucTrack.h.

References DstMucTrack::m_ecLastLayer.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_ecLastLayer = layer; }

void DstMucTrack::setEndPart

(

int

part

)

[inline, inherited]

Definition at line 81 of file DstMucTrack.h.

References DstMucTrack::m_endPart.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_endPart = part;   }

void RecMucTrack::setExpHits

(

vector< int > &

exphits

)

[inline]

Definition at line 339 of file RecMucTrack.h.

References m_expHits.

{ m_expHits = exphits;}  

void RecMucTrack::setExpHits

(

vector< MucRecHit * > &

pHits

)

Definition at line 330 of file RecMucTrack.cxx.

References genRecEmupikp::i, and m_expHits.

Referenced by ComputeTrackInfo(), and RecMucTrackCnv::TObjectToDataObject().

{ 
    for(int i = 0 ; i < pHits.size(); i++)
          m_expHits.push_back((pHits[i]->GetID()).get_value());
}

void RecMucTrack::setExtDistHits

(

vector< float > &

disthits

)

[inline]

Definition at line 345 of file RecMucTrack.h.

References m_distHits_ext.

{ m_distHits_ext = disthits;}

void RecMucTrack::SetExtMucMomentum	(	const float	px,
		const float	py,
		const float	pz
	)

set start moment of ext track in Muc.

Definition at line 276 of file RecMucTrack.cxx.

References m_ExtMucMomentum.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_ExtMucMomentum.set(px*1000, py*1000, pz*1000);
}

void RecMucTrack::SetExtMucPos	(	const float	x,
		const float	y,
		const float	z
	)

set start position of ext track in Muc. (compute from MdcPos MdcMomentum or get from ExtTrack)

Definition at line 268 of file RecMucTrack.cxx.

References m_ExtMucPos.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_ExtMucPos.set(x*10, y*10, z*10);
}

void RecMucTrack::SetExtTrack

(

RecExtTrack *

extTrack

)

set Ext track point.

Definition at line 314 of file RecMucTrack.cxx.

References DstExtTrack::GetTrackId(), m_ExtTrack, and m_ExtTrackID.

Referenced by MucRecTrkExt::execute().

{
  m_ExtTrack = extTrack;
  if (m_ExtTrack) m_ExtTrackID = extTrack->GetTrackId(); 
}

void RecMucTrack::SetExtTrackID

(

int

id

)

[inline]

set Ext track id. for compute from mdc myself

Definition at line 132 of file RecMucTrack.h.

References m_ExtTrackID.

Referenced by MucRecTrkExt::execute().

{ m_ExtTrackID = id; }

void DstMucTrack::setId

(

int

id

)

[inline, inherited]

Definition at line 76 of file DstMucTrack.h.

References DstMucTrack::m_id.

Referenced by MucRecTrkExt::execute(), MucRecRoadFinder::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_id = id ; }

void DstMucTrack::setkalbrLastLayer

(

int

br

)

[inline, inherited]

Definition at line 113 of file DstMucTrack.h.

References DstMucTrack::m_kalbrLastLayer.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_kalbrLastLayer = br;}

void DstMucTrack::setkalDepth

(

double

de

)

[inline, inherited]

Definition at line 112 of file DstMucTrack.h.

References DstMucTrack::m_kaldepth.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_kaldepth = de;}

void DstMucTrack::setkalDof

(

int

f

)

[inline, inherited]

Definition at line 111 of file DstMucTrack.h.

References DstMucTrack::m_kaldof.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_kaldof = f;}

void DstMucTrack::setkalecLastLayer

(

int

ec

)

[inline, inherited]

Definition at line 114 of file DstMucTrack.h.

References DstMucTrack::m_kalecLastLayer.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_kalecLastLayer = ec;}

void DstMucTrack::setkalRechi2

(

double

ch

)

[inline, inherited]

Definition at line 110 of file DstMucTrack.h.

References DstMucTrack::m_kalrechi2.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_kalrechi2=ch;}

void DstMucTrack::setMaxHitsInLayer

(

int

maxHits

)

[inline, inherited]

Definition at line 86 of file DstMucTrack.h.

References DstMucTrack::m_maxHitsInLayer.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_maxHitsInLayer = maxHits; }

void RecMucTrack::SetMdcMomentum	(	const float	px,
		const float	py,
		const float	pz
	)

set start moment of the track in Mdc.

Definition at line 227 of file RecMucTrack.cxx.

References m_MdcMomentum.

Referenced by MucRecTrkExt::execute().

{
  m_MdcMomentum.set(px*1000, py*1000, pz*1000);
}

void RecMucTrack::SetMdcPos	(	const float	x,
		const float	y,
		const float	z
	)

set start position of the track in Mdc.

Definition at line 219 of file RecMucTrack.cxx.

References m_MdcPos.

Referenced by MucRecTrkExt::execute().

{
  m_MdcPos.set(x*10, y*10, z*10);
}

void RecMucTrack::SetMucMomentum	(	const float	px,
		const float	py,
		const float	pz
	)

set start moment of the track in Muc.

Definition at line 257 of file RecMucTrack.cxx.

References m_MucMomentum, DstMucTrack::m_px, DstMucTrack::m_py, and DstMucTrack::m_pz.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_MucMomentum.set(px*1000, py*1000, pz*1000);
  m_px = px*1000;
  m_py = py*1000;
  m_pz = pz*1000;
}

void RecMucTrack::SetMucPos	(	const float	x,
		const float	y,
		const float	z
	)

set start position of the track in Muc. (after line fit and correction)

Definition at line 235 of file RecMucTrack.cxx.

References m_MucPos, DstMucTrack::m_xPos, DstMucTrack::m_yPos, and DstMucTrack::m_zPos.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{
  m_MucPos.set(x*10, y*10, z*10);
  m_xPos = x*10;
  m_yPos = y*10;
  m_zPos = z*10;
}

void RecMucTrack::SetMucPosSigma	(	const float	sigmax,
		const float	sigmay,
		const float	sigmaz
	)

Definition at line 246 of file RecMucTrack.cxx.

References m_MucPosSigma, DstMucTrack::m_xPosSigma, DstMucTrack::m_yPosSigma, and DstMucTrack::m_zPosSigma.

{
  m_MucPosSigma.set(x*10, y*10, z*10);
  m_xPosSigma = x*10;
  m_yPosSigma = y*10;
  m_zPosSigma = z*10;
}

void DstMucTrack::setNumHits

(

int

numHits

)

[inline, inherited]

Definition at line 84 of file DstMucTrack.h.

References DstMucTrack::m_numHits.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_numHits  = numHits; }

void DstMucTrack::setNumLayers

(

int

numLayers

)

[inline, inherited]

Definition at line 85 of file DstMucTrack.h.

References DstMucTrack::m_numLayers.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_numLayers = numLayers; }

void DstMucTrack::setPx

(

double

px

)

[inline, inherited]

Definition at line 101 of file DstMucTrack.h.

References DstMucTrack::m_px.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_px = px;  }

void DstMucTrack::setPy

(

double

py

)

[inline, inherited]

Definition at line 102 of file DstMucTrack.h.

References DstMucTrack::m_py.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_py = py;  }

void DstMucTrack::setPz

(

double

pz

)

[inline, inherited]

Definition at line 103 of file DstMucTrack.h.

References DstMucTrack::m_pz.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_pz = pz;  }

void RecMucTrack::setQuadDistHits

(

vector< float > &

disthits

)

[inline]

Definition at line 343 of file RecMucTrack.h.

References m_distHits_quad.

{ m_distHits_quad = disthits;}

void RecMucTrack::SetRecMode

(

int

recmode

)

[inline]

Definition at line 254 of file RecMucTrack.h.

References m_recmode.

Referenced by MucRecTrkExt::execute(), and MucRecRoadFinder::execute().

{m_recmode = recmode;}

void DstMucTrack::setRms

(

double

rms

)

[inline, inherited]

Definition at line 91 of file DstMucTrack.h.

References DstMucTrack::m_rms.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_rms    = rms ; }

void DstMucTrack::setStartPart

(

int

part

)

[inline, inherited]

Definition at line 80 of file DstMucTrack.h.

References DstMucTrack::m_startPart.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{m_startPart = part; }

void DstMucTrack::setStatus

(

int

st

)

[inline, inherited]

Definition at line 77 of file DstMucTrack.h.

References DstMucTrack::m_status.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_status = st;  }

void RecMucTrack::setTrackId

(

const int

trackId

)

set the index for this track.

Reimplemented from DstMucTrack.

Definition at line 209 of file RecMucTrack.cxx.

References DstMucTrack::m_trackId.

Referenced by MucRecTrkExt::execute(), MucRecRoadFinder::execute(), and RecMucTrackCnv::TObjectToDataObject().

{
  if(trackId >= 0) {
    m_trackId = trackId;
  }
}

void DstMucTrack::setType

(

int

type

)

[inline, inherited]

Definition at line 78 of file DstMucTrack.h.

References DstMucTrack::m_type.

Referenced by MucRecTrkExt::execute(), RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_type = type; }

void RecMucTrack::setVecHits

(

vector< int > &

vechits

)

[inline]

Reimplemented from DstMucTrack.

Definition at line 337 of file RecMucTrack.h.

References m_vecHits.

{ m_vecHits = vechits;}

void RecMucTrack::setVecHits

(

vector< MucRecHit * > &

pHits

)

reload setVecHits

Definition at line 323 of file RecMucTrack.cxx.

References genRecEmupikp::i, and m_vecHits.

Referenced by ComputeTrackInfo(), and RecMucTrackCnv::TObjectToDataObject().

{
  for(int i = 0 ; i < pHits.size(); i++)
    m_vecHits.push_back((pHits[i]->GetID()).get_value());
}

void DstMucTrack::setXPos

(

double

x

)

[inline, inherited]

Definition at line 93 of file DstMucTrack.h.

References DstMucTrack::m_xPos.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_xPos = x; }

void DstMucTrack::setXPosSigma

(

double

xsigma

)

[inline, inherited]

Definition at line 97 of file DstMucTrack.h.

References DstMucTrack::m_xPosSigma.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_xPosSigma = xsigma; }

void DstMucTrack::setYPos

(

double

y

)

[inline, inherited]

Definition at line 94 of file DstMucTrack.h.

References DstMucTrack::m_yPos.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_yPos = y; }

void DstMucTrack::setYPosSigma

(

double

ysigma

)

[inline, inherited]

Definition at line 98 of file DstMucTrack.h.

References DstMucTrack::m_yPosSigma.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_yPosSigma = ysigma; }

void DstMucTrack::setZPos

(

double

z

)

[inline, inherited]

Definition at line 95 of file DstMucTrack.h.

References DstMucTrack::m_zPos.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_zPos = z; }

void DstMucTrack::setZPosSigma

(

double

zsigma

)

[inline, inherited]

Definition at line 99 of file DstMucTrack.h.

References DstMucTrack::m_zPosSigma.

Referenced by RecMucTrackCnv::TObjectToDataObject(), and MucTrackCnv::TObjectToDataObject().

{ m_zPosSigma = zsigma; }

int DstMucTrack::startPart

(

)

const [inline, inherited]

Definition at line 37 of file DstMucTrack.h.

References DstMucTrack::m_startPart.

{ return   m_startPart; }    // start position of track locates in which part

int DstMucTrack::status

(

)

const [inline, inherited]

Definition at line 34 of file DstMucTrack.h.

References DstMucTrack::m_status.

Referenced by LineFit().

{ return   m_status; } 

int DstMucTrack::trackId

(

)

const [inline, inherited]

Definition at line 32 of file DstMucTrack.h.

References DstMucTrack::m_trackId.

Referenced by MucRecTrkExt::execute(), MucRecRoadFinder::execute(), and MucRecRoadFinder::TrackFinding().

{ return   m_trackId ; }

int DstMucTrack::type

(

)

const [inline, inherited]

Definition at line 35 of file DstMucTrack.h.

References DstMucTrack::m_type.

{ return   m_type;  }

vector<int> DstMucTrack::vecHits

(

)

const [inline, inherited]

Definition at line 65 of file DstMucTrack.h.

References DstMucTrack::m_vecHits.

{ return  m_vecHits ;}

double DstMucTrack::xPos

(

)

const [inline, inherited]

Definition at line 50 of file DstMucTrack.h.

References DstMucTrack::m_xPos.

{ return m_xPos; }   // position on the start of track in muc (vertex) 

double DstMucTrack::xPosSigma

(

)

const [inline, inherited]

Definition at line 54 of file DstMucTrack.h.

References DstMucTrack::m_xPosSigma.

{ return m_xPosSigma; } 

double DstMucTrack::yPos

(

)

const [inline, inherited]

Definition at line 51 of file DstMucTrack.h.

References DstMucTrack::m_yPos.

{ return m_yPos; }

double DstMucTrack::yPosSigma

(

)

const [inline, inherited]

Definition at line 55 of file DstMucTrack.h.

References DstMucTrack::m_yPosSigma.

{ return m_yPosSigma; }

double DstMucTrack::zPos

(

)

const [inline, inherited]

Definition at line 52 of file DstMucTrack.h.

References DstMucTrack::m_zPos.

{ return m_zPos; }

double DstMucTrack::zPosSigma

(

)

const [inline, inherited]

Definition at line 56 of file DstMucTrack.h.

References DstMucTrack::m_zPosSigma.

{ return m_zPosSigma; }

---

Member Data Documentation

int RecMucTrack::m_brFirstLayer [private]

Definition at line 369 of file RecMucTrack.h.

Referenced by brFirstLayer(), ComputeLastGap(), RecMucTrack(), and SetDefault().

int DstMucTrack::m_brLastLayer [protected, inherited]

Definition at line 125 of file DstMucTrack.h.

Referenced by DstMucTrack::brLastLayer(), ComputeDepth(), ComputeLastGap(), operator=(), RecMucTrack(), and DstMucTrack::setBrLastLayer().

bool RecMucTrack::m_changeUnit [private]

Definition at line 377 of file RecMucTrack.h.

Referenced by operator=(), OutputUnitChange(), and RecMucTrack().

double DstMucTrack::m_chi2 [protected, inherited]

Definition at line 132 of file DstMucTrack.h.

Referenced by DstMucTrack::chi2(), LineFit(), operator=(), RecMucTrack(), and DstMucTrack::setChi2().

Hep3Vector RecMucTrack::m_CurrentDir [private]

Definition at line 364 of file RecMucTrack.h.

Referenced by CorrectDir(), GetCurrentDir(), operator=(), Project(), and SetCurrentDir().

Hep3Vector RecMucTrack::m_CurrentInsct [private]

Definition at line 365 of file RecMucTrack.h.

Referenced by CalculateInsct(), CorrectDir(), CorrectPos(), GetCurrentInsct(), GetHitDistance(), GetHitDistance2(), operator=(), and SetCurrentInsct().

Hep3Vector RecMucTrack::m_CurrentPos [private]

Definition at line 363 of file RecMucTrack.h.

Referenced by CorrectDir(), CorrectPos(), GetCurrentPos(), operator=(), Project(), and SetCurrentPos().

double DstMucTrack::m_deltaPhi [protected, inherited]

Definition at line 149 of file DstMucTrack.h.

Referenced by DstMucTrack::deltaPhi(), LineFit(), operator=(), RecMucTrack(), and DstMucTrack::setDeltaPhi().

double DstMucTrack::m_depth [protected, inherited]

Definition at line 131 of file DstMucTrack.h.

Referenced by ComputeDepth(), DstMucTrack::depth(), LineFit(), operator=(), OutputUnitChange(), RecMucTrack(), and DstMucTrack::setDepth().

float RecMucTrack::m_depth_3 [private]

Definition at line 373 of file RecMucTrack.h.

Referenced by ComputeDepth(), and GetDepth3().

vector<Hep3Vector> RecMucTrack::m_Directions [private]

Definition at line 390 of file RecMucTrack.h.

Referenced by AttachDirection(), and operator=().

Hep3Vector RecMucTrack::m_Dist_muc_ext [private]

Definition at line 366 of file RecMucTrack.h.

Referenced by ComputeDistanceMatch(), and GetExtMucDistance().

double DstMucTrack::m_distance [protected, inherited]

Definition at line 148 of file DstMucTrack.h.

Referenced by ComputeDistanceMatch(), DstMucTrack::distance(), OutputUnitChange(), and DstMucTrack::setDistance().

vector<float> RecMucTrack::m_distHits [private]

Definition at line 384 of file RecMucTrack.h.

Referenced by getDistHits(), LineFit(), and setDistHits().

vector<float> RecMucTrack::m_distHits_ext [private]

Definition at line 386 of file RecMucTrack.h.

Referenced by getExtDistHits(), pushExtDistHits(), and setExtDistHits().

vector<float> RecMucTrack::m_distHits_quad [private]

Definition at line 385 of file RecMucTrack.h.

Referenced by getQuadDistHits(), LineFit(), and setQuadDistHits().

int DstMucTrack::m_dof [protected, inherited]

Definition at line 133 of file DstMucTrack.h.

Referenced by DstMucTrack::dof(), LineFit(), operator=(), RecMucTrack(), and DstMucTrack::setDof().

int RecMucTrack::m_ecFirstLayer [private]

Definition at line 370 of file RecMucTrack.h.

Referenced by ComputeLastGap(), ecFirstLayer(), RecMucTrack(), and SetDefault().

int DstMucTrack::m_ecLastLayer [protected, inherited]

Definition at line 126 of file DstMucTrack.h.

Referenced by ComputeDepth(), ComputeLastGap(), DstMucTrack::ecLastLayer(), LineFit(), operator=(), RecMucTrack(), and DstMucTrack::setEcLastLayer().

int RecMucTrack::m_ecPart [private]

Definition at line 371 of file RecMucTrack.h.

Referenced by ComputeLastGap(), and RecMucTrack().

int DstMucTrack::m_endPart [protected, inherited]

Definition at line 124 of file DstMucTrack.h.

Referenced by ComputeLastGap(), DstMucTrack::endPart(), LineFit(), operator=(), RecMucTrack(), and DstMucTrack::setEndPart().

vector<int> RecMucTrack::m_expHits [private]

Definition at line 383 of file RecMucTrack.h.

Referenced by getExpHits(), and setExpHits().

Hep3Vector RecMucTrack::m_ExtMucMomentum [private]

Definition at line 357 of file RecMucTrack.h.

Referenced by getExtMucMomentum(), and SetExtMucMomentum().

Hep3Vector RecMucTrack::m_ExtMucPos [private]

Definition at line 356 of file RecMucTrack.h.

Referenced by ComputeDistanceMatch(), Extend(), getExtMucPos(), and SetExtMucPos().

RecExtTrack* RecMucTrack::m_ExtTrack [private]

Definition at line 351 of file RecMucTrack.h.

Referenced by SetExtTrack().

int RecMucTrack::m_ExtTrackID [private]

Definition at line 350 of file RecMucTrack.h.

Referenced by getExtTrackID(), operator=(), SetExtTrack(), and SetExtTrackID().

int RecMucTrack::m_Good3DLine [private]

Definition at line 375 of file RecMucTrack.h.

Referenced by ComputeDepth(), and LineFit().

int RecMucTrack::m_Good3DPart [private]

Definition at line 376 of file RecMucTrack.h.

Referenced by ComputeDepth(), LineFit(), and SetDefault().

int DstMucTrack::m_id [protected, inherited]

Definition at line 119 of file DstMucTrack.h.

Referenced by DstMucTrack::id(), operator=(), RecMucTrack(), and DstMucTrack::setId().

Hep3Vector RecMucTrack::m_IntersectionInner[9] [private]

Definition at line 378 of file RecMucTrack.h.

Referenced by ComputeDepth(), LineFit(), and RecMucTrack().

Hep3Vector RecMucTrack::m_IntersectionOuter[9] [private]

Definition at line 379 of file RecMucTrack.h.

Referenced by ComputeDepth(), LineFit(), and RecMucTrack().

vector<Hep3Vector> RecMucTrack::m_Intersections [private]

Definition at line 389 of file RecMucTrack.h.

Referenced by AttachInsct(), and operator=().

int DstMucTrack::m_kalbrLastLayer [protected, inherited]

Definition at line 157 of file DstMucTrack.h.

Referenced by DstMucTrack::kalbrLastLayer(), operator=(), RecMucTrack(), and DstMucTrack::setkalbrLastLayer().

double DstMucTrack::m_kaldepth [protected, inherited]

Definition at line 156 of file DstMucTrack.h.

Referenced by DstMucTrack::kaldepth(), operator=(), RecMucTrack(), and DstMucTrack::setkalDepth().

int DstMucTrack::m_kaldof [protected, inherited]

Definition at line 155 of file DstMucTrack.h.

Referenced by DstMucTrack::kaldof(), operator=(), RecMucTrack(), and DstMucTrack::setkalDof().

int DstMucTrack::m_kalecLastLayer [protected, inherited]

Definition at line 158 of file DstMucTrack.h.

Referenced by DstMucTrack::kalecLastLayer(), operator=(), RecMucTrack(), and DstMucTrack::setkalecLastLayer().

double DstMucTrack::m_kalrechi2 [protected, inherited]

Definition at line 154 of file DstMucTrack.h.

Referenced by DstMucTrack::kalRechi2(), operator=(), RecMucTrack(), and DstMucTrack::setkalRechi2().

int DstMucTrack::m_maxHitsInLayer [protected, inherited]

Definition at line 129 of file DstMucTrack.h.

Referenced by ComputeMaxHitsInGap(), DstMucTrack::maxHitsInLayer(), operator=(), RecMucTrack(), and DstMucTrack::setMaxHitsInLayer().

Hep3Vector RecMucTrack::m_MdcMomentum [private]

Definition at line 354 of file RecMucTrack.h.

Referenced by getMdcMomentum(), LineFit(), operator=(), and SetMdcMomentum().

Hep3Vector RecMucTrack::m_MdcPos [private]

Definition at line 353 of file RecMucTrack.h.

Referenced by operator=(), and SetMdcPos().

Hep3Vector RecMucTrack::m_MucMomentum [private]

Definition at line 360 of file RecMucTrack.h.

Referenced by ComputeDepth(), Extend(), getMucMomentum(), LineFit(), operator=(), and SetMucMomentum().

Hep3Vector RecMucTrack::m_MucPos [private]

Definition at line 359 of file RecMucTrack.h.

Referenced by ComputeDistanceMatch(), Extend(), getMucPos(), LineFit(), operator=(), and SetMucPos().

Hep3Vector RecMucTrack::m_MucPosSigma [private]

Definition at line 361 of file RecMucTrack.h.

Referenced by getMucPosSigma(), LineFit(), operator=(), and SetMucPosSigma().

int DstMucTrack::m_numHits [protected, inherited]

Definition at line 127 of file DstMucTrack.h.

Referenced by ComputeNGapsWithHits(), DstMucTrack::numHits(), operator=(), RecMucTrack(), and DstMucTrack::setNumHits().

int DstMucTrack::m_numLayers [protected, inherited]

Definition at line 128 of file DstMucTrack.h.

Referenced by ComputeDepth(), ComputeNGapsWithHits(), GetNGapsWithHits(), DstMucTrack::numLayers(), operator=(), RecMucTrack(), and DstMucTrack::setNumLayers().

vector<MucRecHit*> RecMucTrack::m_pExpectedHits [private]

Definition at line 388 of file RecMucTrack.h.

Referenced by ComputeTrackInfo(), GetExpectedHits(), LineFit(), operator=(), and ~RecMucTrack().

vector<MucRecHit*> RecMucTrack::m_pHits [private]

Definition at line 387 of file RecMucTrack.h.

Referenced by AttachHit(), CalculateInsct(), ComputeDepth(), ComputeDistanceMatch(), ComputeLastGap(), ComputeNGapsWithHits(), ComputeTrackInfo(), GetHit(), GetHitIndices(), GetHitInGap(), GetHitInSegOrient(), GetHits(), GetNSharedHits(), GetTotalHits(), HasHit(), HasHitInGap(), LineFit(), operator=(), and PrintHitsInfo().

double DstMucTrack::m_px [protected, inherited]

Definition at line 144 of file DstMucTrack.h.

Referenced by LineFit(), OutputUnitChange(), DstMucTrack::px(), SetMucMomentum(), and DstMucTrack::setPx().

double DstMucTrack::m_py [protected, inherited]

Definition at line 145 of file DstMucTrack.h.

Referenced by LineFit(), OutputUnitChange(), DstMucTrack::py(), SetMucMomentum(), and DstMucTrack::setPy().

double DstMucTrack::m_pz [protected, inherited]

Definition at line 146 of file DstMucTrack.h.

Referenced by LineFit(), OutputUnitChange(), DstMucTrack::pz(), SetMucMomentum(), and DstMucTrack::setPz().

int RecMucTrack::m_recmode [private]

Definition at line 374 of file RecMucTrack.h.

Referenced by GetRecMode(), operator=(), RecMucTrack(), and SetRecMode().

double DstMucTrack::m_rms [protected, inherited]

Definition at line 134 of file DstMucTrack.h.

Referenced by LineFit(), operator=(), RecMucTrack(), DstMucTrack::rms(), and DstMucTrack::setRms().

int DstMucTrack::m_startPart [protected, inherited]

Definition at line 123 of file DstMucTrack.h.

Referenced by operator=(), RecMucTrack(), DstMucTrack::setStartPart(), and DstMucTrack::startPart().

int DstMucTrack::m_status [protected, inherited]

Definition at line 120 of file DstMucTrack.h.

Referenced by LineFit(), operator=(), RecMucTrack(), DstMucTrack::setStatus(), and DstMucTrack::status().

Hep3Vector RecMucTrack::m_StripPhi [private]

Definition at line 367 of file RecMucTrack.h.

Referenced by ComputeDistanceMatch(), and GetMucStripPos().

int DstMucTrack::m_trackId [protected, inherited]

Definition at line 118 of file DstMucTrack.h.

Referenced by setTrackId(), DstMucTrack::setTrackId(), and DstMucTrack::trackId().

int DstMucTrack::m_type [protected, inherited]

Definition at line 121 of file DstMucTrack.h.

Referenced by operator=(), RecMucTrack(), DstMucTrack::setType(), and DstMucTrack::type().

vector<int> RecMucTrack::m_vecHits [private]

Reimplemented from DstMucTrack.

Definition at line 382 of file RecMucTrack.h.

Referenced by getVecHits(), and setVecHits().

double DstMucTrack::m_xPos [protected, inherited]

Definition at line 136 of file DstMucTrack.h.

Referenced by Extend(), LineFit(), OutputUnitChange(), SetMucPos(), DstMucTrack::setXPos(), and DstMucTrack::xPos().

double DstMucTrack::m_xPosSigma [protected, inherited]

Definition at line 140 of file DstMucTrack.h.

Referenced by LineFit(), operator=(), OutputUnitChange(), RecMucTrack(), SetMucPosSigma(), DstMucTrack::setXPosSigma(), and DstMucTrack::xPosSigma().

double DstMucTrack::m_yPos [protected, inherited]

Definition at line 137 of file DstMucTrack.h.

Referenced by Extend(), LineFit(), OutputUnitChange(), SetMucPos(), DstMucTrack::setYPos(), and DstMucTrack::yPos().

double DstMucTrack::m_yPosSigma [protected, inherited]

Definition at line 141 of file DstMucTrack.h.

Referenced by LineFit(), operator=(), OutputUnitChange(), RecMucTrack(), SetMucPosSigma(), DstMucTrack::setYPosSigma(), and DstMucTrack::yPosSigma().

double DstMucTrack::m_zPos [protected, inherited]

Definition at line 138 of file DstMucTrack.h.

Referenced by Extend(), LineFit(), OutputUnitChange(), SetMucPos(), DstMucTrack::setZPos(), and DstMucTrack::zPos().

double DstMucTrack::m_zPosSigma [protected, inherited]

Definition at line 142 of file DstMucTrack.h.

Referenced by LineFit(), operator=(), OutputUnitChange(), RecMucTrack(), SetMucPosSigma(), DstMucTrack::setZPosSigma(), and DstMucTrack::zPosSigma().

---