RecMucTrack Class Reference

#include <RecMucTrack.h>

Inheritance diagram for RecMucTrack:

List of all members.

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

---

Constructor & Destructor Documentation

RecMucTrack::RecMucTrack (   )

Constructor. : //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; }

RecMucTrack::RecMucTrack (  const RecMucTrack &  source  )

Copy constructor. : //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; }

RecMucTrack::RecMucTrack (  const DstMucTrack &  dstTrack  )

Copy constructor from DstMucTrack. :DstMucTrack(dstTrack) { SetDefault(); }

RecMucTrack::~RecMucTrack (   )

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

RecMucTrack::RecMucTrack (   )

Constructor.

RecMucTrack::RecMucTrack (  const RecMucTrack &  source  )

Copy constructor.

RecMucTrack::RecMucTrack (  const DstMucTrack &  dstTrack  )

Copy constructor from DstMucTrack.

RecMucTrack::~RecMucTrack (   )

Destructor.

---

Member Function Documentation

void RecMucTrack::AttachDirection (  Hep3Vector  dir  )

Attach the direction to this trajectory.

void RecMucTrack::AttachDirection (  Hep3Vector  dir  )

Attach the direction to this trajectory. { m_Directions.push_back(dir); }

void RecMucTrack::AttachHit (  MucRecHit *  hit  )

Attach the given hit to this track.

void RecMucTrack::AttachHit (  MucRecHit *  hit  )

Attach the given hit to this track. { // 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.

void RecMucTrack::AttachInsct (  Hep3Vector  insct  )

Attach the intersection to this trajectory. { m_Intersections.push_back(insct); }

int RecMucTrack::brFirstLayer (   )  const [inline]

Which gap on Barrel is the first one with hits attached to this track? { return m_brFirstLayer; }

int RecMucTrack::brFirstLayer (   )  const [inline]

Which gap on Barrel is the first one with hits attached to this track? { return m_brFirstLayer; }

int DstMucTrack::brLastLayer (   )  const [inline, inherited]

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

int DstMucTrack::brLastLayer (   )  const [inline, inherited]

{ 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.

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

Calculate intersection from all hits attached on this gap. { 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]

{ return m_chi2; }

double DstMucTrack::chi2 (   )  const [inline, inherited]

{ return m_chi2; }

const CLID& RecMucTrack::classID (   )  [inline, static]

Reimplemented from DstMucTrack. { return CLID_RecMucTrack; }

const CLID& RecMucTrack::classID (   )  [inline, static]

Reimplemented from DstMucTrack. { return CLID_RecMucTrack; }

virtual const CLID& RecMucTrack::clID (   )  const [inline, virtual]

Reimplemented from DstMucTrack. { return RecMucTrack::classID(); }

virtual const CLID& RecMucTrack::clID (   )  const [inline, virtual]

Reimplemented from DstMucTrack. { return RecMucTrack::classID(); }

int RecMucTrack::ComputeDepth (  int  method  )

Compute depth.

void RecMucTrack::ComputeDepth (   )

Compute depth.

int RecMucTrack::ComputeDepth (  int  method  )

Compute depth. { 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. { 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.

void RecMucTrack::ComputeDistanceMatch (   )

Compute distance match //2006.11.08. { 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.

void RecMucTrack::ComputeLastGap (   )

Comute last gap in barrel and endcap. { 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;.

void RecMucTrack::ComputeMaxHitsInGap (   )

ComputeMaxHitsInGap;. { 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;.

void RecMucTrack::ComputeNGapsWithHits (   )

ComputeNGapsWithHits;. { 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;.

void RecMucTrack::ComputeTrackInfo (  int  fittingmethod  )

Compute all infomation of this track;. { setVecHits(m_pHits); setExpHits(m_pExpectedHits); ComputeLastGap(); ComputeNGapsWithHits(); ComputeMaxHitsInGap(); //ComputeDepth(); ComputeDepth(fittingmethod); ComputeDistanceMatch(); OutputUnitChange(); }

void RecMucTrack::CorrectDir (   )

Correct direction of this trajectory.

void RecMucTrack::CorrectDir (   )

Correct direction of this trajectory. { //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.

void RecMucTrack::CorrectPos (   )

Correct current position of this trajectory. { m_CurrentPos = m_CurrentInsct; }

double DstMucTrack::deltaPhi (   )  const [inline, inherited]

{ return m_deltaPhi;}

double DstMucTrack::deltaPhi (   )  const [inline, inherited]

{ return m_deltaPhi;}

double DstMucTrack::depth (   )  const [inline, inherited]

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

double DstMucTrack::depth (   )  const [inline, inherited]

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

double DstMucTrack::distance (   )  const [inline, inherited]

{ return m_distance;}

double DstMucTrack::distance (   )  const [inline, inherited]

{ return m_distance;}

int DstMucTrack::dof (   )  const [inline, inherited]

{ return m_dof; }

int DstMucTrack::dof (   )  const [inline, inherited]

{ return m_dof; }

int RecMucTrack::ecFirstLayer (   )  const [inline]

Which gap on EndCap is the first one with hits attached to this track? { return m_ecFirstLayer; }

int RecMucTrack::ecFirstLayer (   )  const [inline]

Which gap on EndCap is the first one with hits attached to this track? { return m_ecFirstLayer; }

int DstMucTrack::ecLastLayer (   )  const [inline, inherited]

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

int DstMucTrack::ecLastLayer (   )  const [inline, inherited]

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

int DstMucTrack::endPart (   )  const [inline, inherited]

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

int DstMucTrack::endPart (   )  const [inline, inherited]

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

void RecMucTrack::Extend (   )

Extend mucpos and extmucpos to first layer of muc.

void RecMucTrack::Extend (   )

Extend mucpos and extmucpos to first layer of muc. { 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; while ( ( (fabs(x)<=myMucR) && (fabs(y)<=myMucR) && ((fabs(x-y)/sqrt(2.))<=myMucR) && ((fabs(x+y)/sqrt(2.))<=myMucR) && (fabs(z)<=myMucZ) ) ) { x += step * m_MucMomentum.x(); y += step * m_MucMomentum.y(); z += step * m_MucMomentum.z(); //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){ 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.

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

Find the segment which contains most hits, return max hits number. { 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. { return m_CurrentDir/1000; }

Hep3Vector RecMucTrack::GetCurrentDir (   )  const [inline]

Current direction. { return m_CurrentDir/1000; }

Hep3Vector RecMucTrack::GetCurrentInsct (   )  const [inline]

Current intersection. { return m_CurrentInsct; } //internal

Hep3Vector RecMucTrack::GetCurrentInsct (   )  const [inline]

Current intersection. { return m_CurrentInsct; } //internal

Hep3Vector RecMucTrack::GetCurrentPos (   )  const [inline]

Current position. { return m_CurrentPos/10; }

Hep3Vector RecMucTrack::GetCurrentPos (   )  const [inline]

Current position. { return m_CurrentPos/10; }

float RecMucTrack::GetDepth3 (   )  const [inline]

Length of the track penetrating in iron absorber. { return m_depth_3/10; }

float RecMucTrack::GetDepth3 (   )  const [inline]

Length of the track penetrating in iron absorber. { return m_depth_3/10; }

vector<float> RecMucTrack::getDistHits (   )  const [inline]

{ return m_distHits ;}

vector<float> RecMucTrack::getDistHits (   )  const [inline]

{ return m_distHits ;}

vector<MucRecHit*> RecMucTrack::GetExpectedHits (   )  const

vector< MucRecHit * > RecMucTrack::GetExpectedHits (   )  const

{ return m_pExpectedHits; }

vector<int> RecMucTrack::getExpHits (   )  const [inline]

{ return m_expHits ;}

vector<int> RecMucTrack::getExpHits (   )  const [inline]

{ return m_expHits ;}

vector<float> RecMucTrack::getExtDistHits (   )  const [inline]

{ return m_distHits_ext ;}

vector<float> RecMucTrack::getExtDistHits (   )  const [inline]

{ return m_distHits_ext ;}

Hep3Vector RecMucTrack::GetExtMucDistance (   )  const [inline]

Distance match of the ext track with muc track in first layer. { return m_Dist_muc_ext;}

Hep3Vector RecMucTrack::GetExtMucDistance (   )  const [inline]

Distance match of the ext track with muc track in first layer. { return m_Dist_muc_ext;}

Hep3Vector RecMucTrack::getExtMucMomentum (   )  const [inline]

Start momentum of the Ext track in Muc. { return m_ExtMucMomentum/1000; }

Hep3Vector RecMucTrack::getExtMucMomentum (   )  const [inline]

Start momentum of the Ext track in Muc. { return m_ExtMucMomentum/1000; }

Hep3Vector RecMucTrack::getExtMucPos (   )  const [inline]

start position of the Ext track in Muc. { return m_ExtMucPos/10; }

Hep3Vector RecMucTrack::getExtMucPos (   )  const [inline]

start position of the Ext track in Muc. { return m_ExtMucPos/10; }

int RecMucTrack::getExtTrackID (   )  const [inline]

The identifier of the Ext track as seed. { return m_ExtTrackID; }

int RecMucTrack::getExtTrackID (   )  const [inline]

The identifier of the Ext track as seed. { 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.

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

Get a pointer to the first hit attached in a particular gap. { 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.

float RecMucTrack::GetHitDistance (  const MucRecHit *  hit  )

Calculate the distance of the hit to the intersection in read direction. { 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

float RecMucTrack::GetHitDistance2 (  const MucRecHit *  hit  )

no abs value { 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.

vector< long > RecMucTrack::GetHitIndices (   )  const

Get indices of all hits in the track. { 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?

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

How many hits per gap does this track contain? { 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.

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

How many hits does a segment contains. { 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.

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

How many hits does a segment contains in one orient. { 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.

vector< MucRecHit * > RecMucTrack::GetHits (   )  const

Get all hits on this track. { return m_pHits; }

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

compute ext track myself from mdc.

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

compute ext track myself from mdc. { //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 { return m_MdcMomentum/1000; }

Hep3Vector RecMucTrack::getMdcMomentum (   )  const [inline]

momentum of this track in Mdc { return m_MdcMomentum/1000; }

Hep3Vector RecMucTrack::getMucMomentum (   )  const [inline]

Start momentum of this track in Muc. { return m_MucMomentum/1000; }

Hep3Vector RecMucTrack::getMucMomentum (   )  const [inline]

Start momentum of this track in Muc. { return m_MucMomentum/1000; }

Hep3Vector RecMucTrack::getMucPos (   )  const [inline]

start position of this track in Muc. { return m_MucPos/10; }

Hep3Vector RecMucTrack::getMucPos (   )  const [inline]

start position of this track in Muc. { return m_MucPos/10; }

Hep3Vector RecMucTrack::getMucPosSigma (   )  const [inline]

{ return m_MucPosSigma/10; }

Hep3Vector RecMucTrack::getMucPosSigma (   )  const [inline]

{ return m_MucPosSigma/10; }

Hep3Vector RecMucTrack::GetMucStripPos (   )  const [inline]

{ return m_StripPhi;}

Hep3Vector RecMucTrack::GetMucStripPos (   )  const [inline]

{ return m_StripPhi;}

int RecMucTrack::GetNGapsWithHits (   )  const [inline]

How many gaps provide hits attached to this track? { return m_numLayers; }

int RecMucTrack::GetNGapsWithHits (   )  const [inline]

How many gaps provide hits attached to this track? { return m_numLayers; }

int RecMucTrack::GetNSharedHits (  const RecMucTrack *  track  )  const

How many hits do two tracks share?

int RecMucTrack::GetNSharedHits (  const RecMucTrack *  track  )  const

How many hits do two tracks share? { 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]

{ return m_distHits_quad ;}

vector<float> RecMucTrack::getQuadDistHits (   )  const [inline]

{ return m_distHits_quad ;}

int RecMucTrack::GetRecMode (   )  const [inline]

{return m_recmode;}

int RecMucTrack::GetRecMode (   )  const [inline]

{return m_recmode;}

int RecMucTrack::GetTotalHits (   )  const

How many hits in all does this track contain?

int RecMucTrack::GetTotalHits (   )  const

How many hits in all does this track contain? { return m_pHits.size(); }

vector<int> RecMucTrack::getVecHits (  void   )  const [inline]

{ return m_vecHits ;}

vector<int> RecMucTrack::getVecHits (  void   )  const [inline]

{ 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?

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

Does this track contains an assigned hit? { 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?

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

Does this track contain any hits in the given gap? { 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]

{ return m_id; }

int DstMucTrack::id (   )  const [inline, inherited]

{ return m_id; }

bool RecMucTrack::IsInsideSuperConductor (  Hep3Vector  pos  )

bool RecMucTrack::IsInsideSuperConductor (  Hep3Vector  pos  )

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

void RecMucTrack::LineFit (  int  fittingMethod  )

Line fit with hits on a seg with max hits.

void RecMucTrack::LineFit (  int  fittingMethod  )

Line fit with hits on a seg with max hits. { 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]

{ return m_maxHitsInLayer; }

int DstMucTrack::maxHitsInLayer (   )  const [inline, inherited]

{ return m_maxHitsInLayer; }

int DstMucTrack::numHits (   )  const [inline, inherited]

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

int DstMucTrack::numHits (   )  const [inline, inherited]

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

int DstMucTrack::numLayers (   )  const [inline, inherited]

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

int DstMucTrack::numLayers (   )  const [inline, inherited]

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

RecMucTrack& RecMucTrack::operator= (  const DstMucTrack &  dstTrack  )

Assignment constructor from DstMucTrack.

RecMucTrack& RecMucTrack::operator= (  const RecMucTrack &  orig  )

Assignment constructor.

RecMucTrack & RecMucTrack::operator= (  const DstMucTrack &  dstTrack  )

Assignment constructor from DstMucTrack. { SetDefault(); DstMucTrack::operator=(dstTrack); return *this; }

RecMucTrack & RecMucTrack::operator= (  const RecMucTrack &  orig  )

Assignment constructor. { // 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; } return *this; }

void RecMucTrack::OutputUnitChange (   )

change unit

void RecMucTrack::OutputUnitChange (   )

change unit { 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.

void RecMucTrack::PrintHitsInfo (   )  const

Print Hits Infomation. { 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?

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? { 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]

{m_distHits_ext.push_back(dist);}

void RecMucTrack::pushExtDistHits (  float  dist  )  [inline]

{m_distHits_ext.push_back(dist);}

double DstMucTrack::px (   )  const [inline, inherited]

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

double DstMucTrack::px (   )  const [inline, inherited]

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

double DstMucTrack::py (   )  const [inline, inherited]

{ return m_py; }

double DstMucTrack::py (   )  const [inline, inherited]

{ return m_py; }

double DstMucTrack::pz (   )  const [inline, inherited]

{ return m_pz; }

double DstMucTrack::pz (   )  const [inline, inherited]

{ return m_pz; }

double DstMucTrack::rms (  void   )  const [inline, inherited]

{ return m_rms; }

double DstMucTrack::rms (   )  const [inline, inherited]

{ return m_rms; }

void DstMucTrack::setBrLastLayer (  int  layer  )  [inline, inherited]

{m_brLastLayer = layer; }

void DstMucTrack::setBrLastLayer (  int  layer  )  [inline, inherited]

{m_brLastLayer = layer; }

void DstMucTrack::setChi2 (  double  chi2  )  [inline, inherited]

{ m_chi2 = chi2; }

void DstMucTrack::setChi2 (  double  chi2  )  [inline, inherited]

{ m_chi2 = chi2; }

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

set current direction of the trajectory.

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

set current direction of the trajectory. { 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.

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

set current intersection of the trajectory with strip plane. { m_CurrentInsct.set(x, y, z); }

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

set current position of the trajectory.

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

set current position of the trajectory. { m_CurrentPos.set(x*10, y*10, z*10); }

void RecMucTrack::SetDefault (   )

void RecMucTrack::SetDefault (   )

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

void DstMucTrack::setDeltaPhi (  double  dphi  )  [inline, inherited]

{m_deltaPhi = dphi; }

void DstMucTrack::setDeltaPhi (  double  dphi  )  [inline, inherited]

{m_deltaPhi = dphi; }

void DstMucTrack::setDepth (  double  depth  )  [inline, inherited]

{ m_depth = depth; }

void DstMucTrack::setDepth (  double  depth  )  [inline, inherited]

{ m_depth = depth; }

void DstMucTrack::setDistance (  double  dist  )  [inline, inherited]

{m_distance = dist; }

void DstMucTrack::setDistance (  double  dist  )  [inline, inherited]

{m_distance = dist; }

void RecMucTrack::setDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits = disthits;}

void RecMucTrack::setDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits = disthits;}

void DstMucTrack::setDof (  int  dof  )  [inline, inherited]

{ m_dof = dof ; }

void DstMucTrack::setDof (  int  dof  )  [inline, inherited]

{ m_dof = dof ; }

void DstMucTrack::setEcLastLayer (  int  layer  )  [inline, inherited]

{m_ecLastLayer = layer; }

void DstMucTrack::setEcLastLayer (  int  layer  )  [inline, inherited]

{m_ecLastLayer = layer; }

void DstMucTrack::setEndPart (  int  part  )  [inline, inherited]

{m_endPart = part; }

void DstMucTrack::setEndPart (  int  part  )  [inline, inherited]

{m_endPart = part; }

void RecMucTrack::setExpHits (  vector< int > &  exphits  )  [inline]

{ m_expHits = exphits;}

void RecMucTrack::setExpHits (  vector< MucRecHit * > &  pHits  )

void RecMucTrack::setExpHits (  vector< int > &  exphits  )  [inline]

{ m_expHits = exphits;}

void RecMucTrack::setExpHits (  vector< MucRecHit * > &  pHits  )

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

void RecMucTrack::setExtDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits_ext = disthits;}

void RecMucTrack::setExtDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits_ext = disthits;}

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

set start moment of ext track in Muc.

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

set start moment of ext track in Muc. { 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)

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) { m_ExtMucPos.set(x*10, y*10, z*10); }

void RecMucTrack::SetExtTrack (  RecExtTrack *  extTrack  )

set Ext track point.

void RecMucTrack::SetExtTrack (  RecExtTrack *  extTrack  )

set Ext track point. { 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 { m_ExtTrackID = id; }

void RecMucTrack::SetExtTrackID (  int  id  )  [inline]

set Ext track id. for compute from mdc myself { m_ExtTrackID = id; }

void DstMucTrack::setId (  int  id  )  [inline, inherited]

{ m_id = id ; }

void DstMucTrack::setId (  int  id  )  [inline, inherited]

{ m_id = id ; }

void DstMucTrack::setMaxHitsInLayer (  int  maxHits  )  [inline, inherited]

{ m_maxHitsInLayer = maxHits; }

void DstMucTrack::setMaxHitsInLayer (  int  maxHits  )  [inline, inherited]

{ m_maxHitsInLayer = maxHits; }

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

set start moment of the track in Mdc.

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

set start moment of the track in Mdc. { 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.

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

set start position of the track in Mdc. { 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.

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

set start moment of the track in Muc. { 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)

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

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

{ 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]

{ m_numHits = numHits; }

void DstMucTrack::setNumHits (  int  numHits  )  [inline, inherited]

{ m_numHits = numHits; }

void DstMucTrack::setNumLayers (  int  numLayers  )  [inline, inherited]

{ m_numLayers = numLayers; }

void DstMucTrack::setNumLayers (  int  numLayers  )  [inline, inherited]

{ m_numLayers = numLayers; }

void DstMucTrack::setPx (  double  px  )  [inline, inherited]

{ m_px = px; }

void DstMucTrack::setPx (  double  px  )  [inline, inherited]

{ m_px = px; }

void DstMucTrack::setPy (  double  py  )  [inline, inherited]

{ m_py = py; }

void DstMucTrack::setPy (  double  py  )  [inline, inherited]

{ m_py = py; }

void DstMucTrack::setPz (  double  pz  )  [inline, inherited]

{ m_pz = pz; }

void DstMucTrack::setPz (  double  pz  )  [inline, inherited]

{ m_pz = pz; }

void RecMucTrack::setQuadDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits_quad = disthits;}

void RecMucTrack::setQuadDistHits (  vector< float > &  disthits  )  [inline]

{ m_distHits_quad = disthits;}

void RecMucTrack::SetRecMode (  int  recmode  )  [inline]

{m_recmode = recmode;}

void RecMucTrack::SetRecMode (  int  recmode  )  [inline]

{m_recmode = recmode;}

void DstMucTrack::setRms (  double  rms  )  [inline, inherited]

{ m_rms = rms ; }

void DstMucTrack::setRms (  double  rms  )  [inline, inherited]

{ m_rms = rms ; }

void DstMucTrack::setStartPart (  int  part  )  [inline, inherited]

{m_startPart = part; }

void DstMucTrack::setStartPart (  int  part  )  [inline, inherited]

{m_startPart = part; }

void DstMucTrack::setStatus (  int  st  )  [inline, inherited]

{ m_status = st; }

void DstMucTrack::setStatus (  int  st  )  [inline, inherited]

{ m_status = st; }

void RecMucTrack::setTrackId (  const int  trackId  )

set the index for this track. Reimplemented from DstMucTrack.

void RecMucTrack::setTrackId (  const int  trackId  )

set the index for this track. Reimplemented from DstMucTrack. { if(trackId >= 0) { m_trackId = trackId; } }

void DstMucTrack::setType (  int  type  )  [inline, inherited]

{ m_type = type; }

void DstMucTrack::setType (  int  type  )  [inline, inherited]

{ m_type = type; }

void RecMucTrack::setVecHits (  vector< int > &  vechits  )  [inline]

Reimplemented from DstMucTrack. { m_vecHits = vechits;}

void RecMucTrack::setVecHits (  vector< MucRecHit * > &  pHits  )

reload setVecHits

void RecMucTrack::setVecHits (  vector< int > &  vechits  )  [inline]

Reimplemented from DstMucTrack. { m_vecHits = vechits;}

void RecMucTrack::setVecHits (  vector< MucRecHit * > &  pHits  )

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

void DstMucTrack::setXPos (  double  x  )  [inline, inherited]

{ m_xPos = x; }

void DstMucTrack::setXPos (  double  x  )  [inline, inherited]

{ m_xPos = x; }

void DstMucTrack::setXPosSigma (  double  xsigma  )  [inline, inherited]

{ m_xPosSigma = xsigma; }

void DstMucTrack::setXPosSigma (  double  xsigma  )  [inline, inherited]

{ m_xPosSigma = xsigma; }

void DstMucTrack::setYPos (  double  y  )  [inline, inherited]

{ m_yPos = y; }

void DstMucTrack::setYPos (  double  y  )  [inline, inherited]

{ m_yPos = y; }

void DstMucTrack::setYPosSigma (  double  ysigma  )  [inline, inherited]

{ m_yPosSigma = ysigma; }

void DstMucTrack::setYPosSigma (  double  ysigma  )  [inline, inherited]

{ m_yPosSigma = ysigma; }

void DstMucTrack::setZPos (  double  z  )  [inline, inherited]

{ m_zPos = z; }

void DstMucTrack::setZPos (  double  z  )  [inline, inherited]

{ m_zPos = z; }

void DstMucTrack::setZPosSigma (  double  zsigma  )  [inline, inherited]

{ m_zPosSigma = zsigma; }

void DstMucTrack::setZPosSigma (  double  zsigma  )  [inline, inherited]

{ m_zPosSigma = zsigma; }

int DstMucTrack::startPart (   )  const [inline, inherited]

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

int DstMucTrack::startPart (   )  const [inline, inherited]

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

int DstMucTrack::status (  void   )  const [inline, inherited]

{ return m_status; }

int DstMucTrack::status (   )  const [inline, inherited]

{ return m_status; }

int DstMucTrack::trackId (   )  const [inline, inherited]

{ return m_trackId ; }

int DstMucTrack::trackId (   )  const [inline, inherited]

{ return m_trackId ; }

int DstMucTrack::type (  void   )  const [inline, inherited]

{ return m_type; }

int DstMucTrack::type (   )  const [inline, inherited]

{ return m_type; }

vector<int> DstMucTrack::vecHits (   )  const [inline, inherited]

{ return m_vecHits ;}

vector<int> DstMucTrack::vecHits (   )  const [inline, inherited]

{ return m_vecHits ;}

double DstMucTrack::xPos (   )  const [inline, inherited]

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

double DstMucTrack::xPos (   )  const [inline, inherited]

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

double DstMucTrack::xPosSigma (   )  const [inline, inherited]

{ return m_xPosSigma; }

double DstMucTrack::xPosSigma (   )  const [inline, inherited]

{ return m_xPosSigma; }

double DstMucTrack::yPos (   )  const [inline, inherited]

{ return m_yPos; }

double DstMucTrack::yPos (   )  const [inline, inherited]

{ return m_yPos; }

double DstMucTrack::yPosSigma (   )  const [inline, inherited]

{ return m_yPosSigma; }

double DstMucTrack::yPosSigma (   )  const [inline, inherited]

{ return m_yPosSigma; }

double DstMucTrack::zPos (   )  const [inline, inherited]

{ return m_zPos; }

double DstMucTrack::zPos (   )  const [inline, inherited]

{ return m_zPos; }

double DstMucTrack::zPosSigma (   )  const [inline, inherited]

{ return m_zPosSigma; }

double DstMucTrack::zPosSigma (   )  const [inline, inherited]

{ return m_zPosSigma; }

---

Member Data Documentation

int RecMucTrack::m_brFirstLayer [private]

int DstMucTrack::m_brLastLayer [protected, inherited]

bool RecMucTrack::m_changeUnit [private]

double DstMucTrack::m_chi2 [protected, inherited]

Hep3Vector RecMucTrack::m_CurrentDir [private]

Hep3Vector RecMucTrack::m_CurrentInsct [private]

Hep3Vector RecMucTrack::m_CurrentPos [private]

double DstMucTrack::m_deltaPhi [protected, inherited]

double DstMucTrack::m_depth [protected, inherited]

float RecMucTrack::m_depth_3 [private]

vector<Hep3Vector> RecMucTrack::m_Directions [private]

vector<Hep3Vector> RecMucTrack::m_Directions [private]

Hep3Vector RecMucTrack::m_Dist_muc_ext [private]

double DstMucTrack::m_distance [protected, inherited]

vector<float> RecMucTrack::m_distHits [private]

vector<float> RecMucTrack::m_distHits [private]

vector<float> RecMucTrack::m_distHits_ext [private]

vector<float> RecMucTrack::m_distHits_ext [private]

vector<float> RecMucTrack::m_distHits_quad [private]

vector<float> RecMucTrack::m_distHits_quad [private]

int DstMucTrack::m_dof [protected, inherited]

int RecMucTrack::m_ecFirstLayer [private]

int DstMucTrack::m_ecLastLayer [protected, inherited]

int RecMucTrack::m_ecPart [private]

int DstMucTrack::m_endPart [protected, inherited]

vector<int> RecMucTrack::m_expHits [private]

vector<int> RecMucTrack::m_expHits [private]

Hep3Vector RecMucTrack::m_ExtMucMomentum [private]

Hep3Vector RecMucTrack::m_ExtMucPos [private]

RecExtTrack* RecMucTrack::m_ExtTrack [private]

RecExtTrack* RecMucTrack::m_ExtTrack [private]

int RecMucTrack::m_ExtTrackID [private]

int RecMucTrack::m_Good3DLine [private]

int RecMucTrack::m_Good3DPart [private]

int DstMucTrack::m_id [protected, inherited]

Hep3Vector RecMucTrack::m_IntersectionInner [private]

Hep3Vector RecMucTrack::m_IntersectionOuter [private]

vector<Hep3Vector> RecMucTrack::m_Intersections [private]

vector<Hep3Vector> RecMucTrack::m_Intersections [private]

int DstMucTrack::m_maxHitsInLayer [protected, inherited]

Hep3Vector RecMucTrack::m_MdcMomentum [private]

Hep3Vector RecMucTrack::m_MdcPos [private]

Hep3Vector RecMucTrack::m_MucMomentum [private]

Hep3Vector RecMucTrack::m_MucPos [private]

Hep3Vector RecMucTrack::m_MucPosSigma [private]

int DstMucTrack::m_numHits [protected, inherited]

int DstMucTrack::m_numLayers [protected, inherited]

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

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

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

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

double DstMucTrack::m_px [protected, inherited]

double DstMucTrack::m_py [protected, inherited]

double DstMucTrack::m_pz [protected, inherited]

int RecMucTrack::m_recmode [private]

double DstMucTrack::m_rms [protected, inherited]

int DstMucTrack::m_startPart [protected, inherited]

int DstMucTrack::m_status [protected, inherited]

Hep3Vector RecMucTrack::m_StripPhi [private]

int DstMucTrack::m_trackId [protected, inherited]

int DstMucTrack::m_type [protected, inherited]

vector<int> RecMucTrack::m_vecHits [private]

Reimplemented from DstMucTrack.

vector<int> RecMucTrack::m_vecHits [private]

Reimplemented from DstMucTrack.

double DstMucTrack::m_xPos [protected, inherited]

double DstMucTrack::m_xPosSigma [protected, inherited]

double DstMucTrack::m_yPos [protected, inherited]

double DstMucTrack::m_yPosSigma [protected, inherited]

double DstMucTrack::m_zPos [protected, inherited]

double DstMucTrack::m_zPosSigma [protected, inherited]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/MucRecEvent/MucRecEvent/RecMucTrack.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Muc/MucRecEvent/MucRecEvent-00-02-49/MucRecEvent/RecMucTrack.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Muc/MucRecEvent/MucRecEvent-00-02-49/src/RecMucTrack.cxx

---