RecMucTrack Class Reference

#include <RecMucTrack.h>

Inheritance diagram for RecMucTrack:

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)
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)
	set corresponding monte carlo track pointer.
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
	Where does the trajectory of this track intersect the reference plane?
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 ()
	chi2 of line fit
void	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
	How many hits were attached in the gap with the most attached hits?
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)
	Get corresponding monte carlo track pointer.
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)
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)
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)
	set corresponding monte carlo track pointer.
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
	Where does the trajectory of this track intersect the reference plane?
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 ()
	chi2 of line fit
void	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
	How many hits were attached in the gap with the most attached hits?
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)
	Get corresponding monte carlo track pointer.
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)
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)
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)
	set corresponding monte carlo track pointer.
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
	Where does the trajectory of this track intersect the reference plane?
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 ()
	chi2 of line fit
void	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
	How many hits were attached in the gap with the most attached hits?
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)
	Get corresponding monte carlo track pointer.
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)
Public Member Functions inherited from DstMucTrack
	DstMucTrack ()
virtual	~DstMucTrack ()
	Assignment constructor.
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	setTrackId (int trackId)
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	setVecHits (vector< int > &vecHits)
void	setkalRechi2 (double ch)
void	setkalDof (int f)
void	setkalDepth (double de)
void	setkalbrLastLayer (int br)
void	setkalecLastLayer (int ec)
	DstMucTrack ()
virtual	~DstMucTrack ()
	Assignment constructor.
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	setTrackId (int trackId)
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	setVecHits (vector< int > &vecHits)
void	setkalRechi2 (double ch)
void	setkalDof (int f)
void	setkalDepth (double de)
void	setkalbrLastLayer (int br)
void	setkalecLastLayer (int ec)
	DstMucTrack ()
virtual	~DstMucTrack ()
	Assignment constructor.
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	setTrackId (int trackId)
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	setVecHits (vector< int > &vecHits)
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 ()
static const CLID &	classID ()
static const CLID &	classID ()
Static Public Member Functions inherited from DstMucTrack
static const CLID &	classID ()
static const CLID &	classID ()
static const CLID &	classID ()

Additional Inherited Members
Protected Attributes inherited from DstMucTrack
int	m_trackId { 0 }
int	m_id { 0 }
int	m_status { 0 }
int	m_type { 0 }
int	m_startPart { 0 }
int	m_endPart { 0 }
int	m_brLastLayer { 0 }
int	m_ecLastLayer { 0 }
int	m_numHits { 0 }
int	m_numLayers { 0 }
int	m_maxHitsInLayer { 0 }
double	m_depth { 0.0 }
double	m_chi2 { 0.0 }
int	m_dof { 0 }
double	m_rms { 0.0 }
double	m_xPos { 0.0 }
double	m_yPos { 0.0 }
double	m_zPos { 0.0 }
double	m_xPosSigma { 0.0 }
double	m_yPosSigma { 0.0 }
double	m_zPosSigma { 0.0 }
double	m_px { 0.0 }
double	m_py { 0.0 }
double	m_pz { 0.0 }
double	m_distance { 0.0 }
double	m_deltaPhi { 0.0 }
vector< int >	m_vecHits
double	m_kalrechi2 { 0.0 }
int	m_kaldof { 0 }
double	m_kaldepth { 0.0 }
int	m_kalbrLastLayer { 0 }
int	m_kalecLastLayer { 0 }

Detailed Description

Definition at line 40 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

Constructor & Destructor Documentation

RecMucTrack() [1/9]

RecMucTrack::RecMucTrack

(

)

Constructor.

Definition at line 27 of file RecMucTrack.cxx.

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

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

RecMucTrack() [2/9]

RecMucTrack::RecMucTrack

(

const RecMucTrack &

source

)

Copy constructor.

Definition at line 134 of file RecMucTrack.cxx.

    : // 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() [3/9]

RecMucTrack::RecMucTrack

(

const DstMucTrack &

dstTrack

)

Copy constructor from DstMucTrack.

Definition at line 177 of file RecMucTrack.cxx.

                                                      : DstMucTrack( dstTrack ) {
  SetDefault();
}

~RecMucTrack() [1/3]

RecMucTrack::~RecMucTrack

(

)

Destructor.

Definition at line 194 of file RecMucTrack.cxx.

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

RecMucTrack() [4/9]

RecMucTrack::RecMucTrack

(

)

Constructor.

RecMucTrack() [5/9]

RecMucTrack::RecMucTrack

(

const RecMucTrack &

source

)

Copy constructor.

RecMucTrack() [6/9]

RecMucTrack::RecMucTrack

(

const DstMucTrack &

dstTrack

)

Copy constructor from DstMucTrack.

~RecMucTrack() [2/3]

RecMucTrack::~RecMucTrack

(

)

Destructor.

RecMucTrack() [7/9]

RecMucTrack::RecMucTrack

(

)

Constructor.

RecMucTrack() [8/9]

RecMucTrack::RecMucTrack

(

const RecMucTrack &

source

)

Copy constructor.

RecMucTrack() [9/9]

RecMucTrack::RecMucTrack

(

const DstMucTrack &

dstTrack

)

Copy constructor from DstMucTrack.

~RecMucTrack() [3/3]

RecMucTrack::~RecMucTrack

(

)

Destructor.

Member Function Documentation

AttachDirection() [1/3]

void RecMucTrack::AttachDirection

(

Hep3Vector

dir

)

Attach the direction to this trajectory.

Definition at line 536 of file RecMucTrack.cxx.

{ m_Directions.push_back( dir ); }

AttachDirection() [2/3]

void RecMucTrack::AttachDirection

(

Hep3Vector

dir

)

Attach the direction to this trajectory.

AttachDirection() [3/3]

void RecMucTrack::AttachDirection

(

Hep3Vector

dir

)

Attach the direction to this trajectory.

AttachHit() [1/3]

void RecMucTrack::AttachHit

(

MucRecHit *

hit

)

set corresponding monte carlo track pointer.

Attach the given hit to this track.

Definition at line 325 of file RecMucTrack.cxx.

                                            {
  //  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();
}

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

AttachHit() [2/3]

void RecMucTrack::AttachHit

(

MucRecHit *

hit

)

set corresponding monte carlo track pointer.

Attach the given hit to this track.

AttachHit() [3/3]

void RecMucTrack::AttachHit

(

MucRecHit *

hit

)

set corresponding monte carlo track pointer.

Attach the given hit to this track.

AttachInsct() [1/3]

void RecMucTrack::AttachInsct

(

Hep3Vector

insct

)

Attach the intersection to this trajectory.

Definition at line 534 of file RecMucTrack.cxx.

{ m_Intersections.push_back( insct ); }

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

AttachInsct() [2/3]

void RecMucTrack::AttachInsct

(

Hep3Vector

insct

)

Attach the intersection to this trajectory.

AttachInsct() [3/3]

void RecMucTrack::AttachInsct

(

Hep3Vector

insct

)

Attach the intersection to this trajectory.

brFirstLayer() [1/3]

int RecMucTrack::brFirstLayer ( ) const

inline

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

Definition at line 186 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_brFirstLayer; }

Referenced by ComputeDepth().

brFirstLayer() [2/3]

int RecMucTrack::brFirstLayer ( ) const

inline

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

Definition at line 186 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_brFirstLayer; }

brFirstLayer() [3/3]

int RecMucTrack::brFirstLayer ( ) const

inline

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

Definition at line 186 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_brFirstLayer; }

CalculateInsct() [1/3]

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

Calculate intersection from all hits attached on this gap.

Definition at line 469 of file RecMucTrack.cxx.

                                                                                     {
  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;
}

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

CalculateInsct() [2/3]

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

Calculate intersection from all hits attached on this gap.

CalculateInsct() [3/3]

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

Calculate intersection from all hits attached on this gap.

classID() [1/3]

const CLID & RecMucTrack::classID ( )

inlinestatic

Definition at line 44 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return CLID_RecMucTrack; }

Referenced by clID().

classID() [2/3]

const CLID & RecMucTrack::classID ( )

inlinestatic

Definition at line 44 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return CLID_RecMucTrack; }

classID() [3/3]

const CLID & RecMucTrack::classID ( )

inlinestatic

Definition at line 44 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return CLID_RecMucTrack; }

clID() [1/3]

virtual const CLID & RecMucTrack::clID ( ) const

inlinevirtual

Reimplemented from DstMucTrack.

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return RecMucTrack::classID(); }

clID() [2/3]

virtual const CLID & RecMucTrack::clID ( ) const

inlinevirtual

Reimplemented from DstMucTrack.

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return RecMucTrack::classID(); }

clID() [3/3]

virtual const CLID & RecMucTrack::clID ( ) const

inlinevirtual

Reimplemented from DstMucTrack.

Definition at line 42 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return RecMucTrack::classID(); }

ComputeDepth() [1/6]

void RecMucTrack::ComputeDepth

(

)

chi2 of line fit

degree of freedom in line fit rms of distance from hit to track Compute depth.

Definition at line 736 of file RecMucTrack.cxx.

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

Referenced by ComputeTrackInfo().

ComputeDepth() [2/6]

void RecMucTrack::ComputeDepth

(

)

chi2 of line fit

degree of freedom in line fit rms of distance from hit to track Compute depth.

ComputeDepth() [3/6]

void RecMucTrack::ComputeDepth

(

)

chi2 of line fit

degree of freedom in line fit rms of distance from hit to track Compute depth.

ComputeDepth() [4/6]

void RecMucTrack::ComputeDepth

(

int

method

)

Compute depth.

Definition at line 591 of file RecMucTrack.cxx.

                                           {
  if ( m_numLayers == 0 )
  {
    m_depth = m_depth_3 = -99;
    return;
  }
  else if ( m_numLayers == 1 && ( m_brLastLayer == 0 || m_ecLastLayer == 0 ) )
  {
    m_depth = m_depth_3 = 0;
    return;
  }
 
  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;
}

ComputeDepth() [5/6]

void RecMucTrack::ComputeDepth

(

int

method

)

Compute depth.

ComputeDepth() [6/6]

void RecMucTrack::ComputeDepth

(

int

method

)

Compute depth.

ComputeDistanceMatch() [1/3]

void RecMucTrack::ComputeDistanceMatch

(

)

Compute distance match //2006.11.08.

Definition at line 845 of file RecMucTrack.cxx.

                                       {
  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;
}

Referenced by ComputeTrackInfo().

ComputeDistanceMatch() [2/3]

void RecMucTrack::ComputeDistanceMatch

(

)

Compute distance match //2006.11.08.

ComputeDistanceMatch() [3/3]

void RecMucTrack::ComputeDistanceMatch

(

)

Compute distance match //2006.11.08.

ComputeLastGap() [1/3]

void RecMucTrack::ComputeLastGap

(

)

Comute last gap in barrel and endcap.

Definition at line 547 of file RecMucTrack.cxx.

                                 {
  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;
}

Referenced by ComputeTrackInfo(), and LineFit().

ComputeLastGap() [2/3]

void RecMucTrack::ComputeLastGap

(

)

Comute last gap in barrel and endcap.

ComputeLastGap() [3/3]

void RecMucTrack::ComputeLastGap

(

)

Comute last gap in barrel and endcap.

ComputeMaxHitsInGap() [1/3]

void RecMucTrack::ComputeMaxHitsInGap

(

)

ComputeMaxHitsInGap;.

Definition at line 1555 of file RecMucTrack.cxx.

                                      {
  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;
}

Referenced by ComputeTrackInfo().

ComputeMaxHitsInGap() [2/3]

void RecMucTrack::ComputeMaxHitsInGap

(

)

ComputeMaxHitsInGap;.

ComputeMaxHitsInGap() [3/3]

void RecMucTrack::ComputeMaxHitsInGap

(

)

ComputeMaxHitsInGap;.

ComputeNGapsWithHits() [1/3]

void RecMucTrack::ComputeNGapsWithHits

(

)

ComputeNGapsWithHits;.

Definition at line 1439 of file RecMucTrack.cxx.

                                       {
  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;
}

Referenced by ComputeTrackInfo().

ComputeNGapsWithHits() [2/3]

void RecMucTrack::ComputeNGapsWithHits

(

)

ComputeNGapsWithHits;.

ComputeNGapsWithHits() [3/3]

void RecMucTrack::ComputeNGapsWithHits

(

)

ComputeNGapsWithHits;.

ComputeTrackInfo() [1/3]

void RecMucTrack::ComputeTrackInfo

(

int

fittingmethod

)

Get corresponding monte carlo track pointer.

Compute all infomation of this track;

Definition at line 1679 of file RecMucTrack.cxx.

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

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

ComputeTrackInfo() [2/3]

void RecMucTrack::ComputeTrackInfo

(

int

fittingmethod

)

Get corresponding monte carlo track pointer.

Compute all infomation of this track;

ComputeTrackInfo() [3/3]

void RecMucTrack::ComputeTrackInfo

(

int

fittingmethod

)

Get corresponding monte carlo track pointer.

Compute all infomation of this track;

CorrectDir() [1/3]

void RecMucTrack::CorrectDir

(

)

Correct direction of this trajectory.

Definition at line 538 of file RecMucTrack.cxx.

                             {
  // 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;
}

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

CorrectDir() [2/3]

void RecMucTrack::CorrectDir

(

)

Correct direction of this trajectory.

CorrectDir() [3/3]

void RecMucTrack::CorrectDir

(

)

Correct direction of this trajectory.

CorrectPos() [1/3]

void RecMucTrack::CorrectPos

(

)

Correct current position of this trajectory.

Definition at line 545 of file RecMucTrack.cxx.

{ m_CurrentPos = m_CurrentInsct; }

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

CorrectPos() [2/3]

void RecMucTrack::CorrectPos

(

)

Correct current position of this trajectory.

CorrectPos() [3/3]

void RecMucTrack::CorrectPos

(

)

Correct current position of this trajectory.

ecFirstLayer() [1/3]

int RecMucTrack::ecFirstLayer ( ) const

inline

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

Definition at line 189 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_ecFirstLayer; }

ecFirstLayer() [2/3]

int RecMucTrack::ecFirstLayer ( ) const

inline

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

Definition at line 189 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_ecFirstLayer; }

ecFirstLayer() [3/3]

int RecMucTrack::ecFirstLayer ( ) const

inline

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

Definition at line 189 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_ecFirstLayer; }

Extend() [1/3]

void RecMucTrack::Extend

(

)

Extend mucpos and extmucpos to first layer of muc.

Definition at line 958 of file RecMucTrack.cxx.

                         {
  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;
    }
  }
}

Referenced by LineFit().

Extend() [2/3]

void RecMucTrack::Extend

(

)

Extend mucpos and extmucpos to first layer of muc.

Extend() [3/3]

void RecMucTrack::Extend

(

)

Extend mucpos and extmucpos to first layer of muc.

FindSegWithMaxHits() [1/3]

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

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

Definition at line 1535 of file RecMucTrack.cxx.

                                                                 {
  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;
}

Referenced by ComputeDistanceMatch(), and LineFit().

FindSegWithMaxHits() [2/3]

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

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

FindSegWithMaxHits() [3/3]

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

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

GetCurrentDir() [1/3]

Hep3Vector RecMucTrack::GetCurrentDir ( ) const

inline

Current direction.

Definition at line 177 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentDir / 1000; }

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

GetCurrentDir() [2/3]

Hep3Vector RecMucTrack::GetCurrentDir ( ) const

inline

Current direction.

Definition at line 177 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentDir / 1000; }

GetCurrentDir() [3/3]

Hep3Vector RecMucTrack::GetCurrentDir ( ) const

inline

Current direction.

Definition at line 177 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentDir / 1000; }

GetCurrentInsct() [1/3]

Hep3Vector RecMucTrack::GetCurrentInsct ( ) const

inline

Current intersection.

Definition at line 180 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentInsct; } // internal

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

GetCurrentInsct() [2/3]

Hep3Vector RecMucTrack::GetCurrentInsct ( ) const

inline

Current intersection.

Definition at line 180 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentInsct; } // internal

GetCurrentInsct() [3/3]

Hep3Vector RecMucTrack::GetCurrentInsct ( ) const

inline

Current intersection.

Definition at line 180 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentInsct; } // internal

GetCurrentPos() [1/3]

Hep3Vector RecMucTrack::GetCurrentPos ( ) const

inline

Current position.

Definition at line 174 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentPos / 10; }

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

GetCurrentPos() [2/3]

Hep3Vector RecMucTrack::GetCurrentPos ( ) const

inline

Current position.

Definition at line 174 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentPos / 10; }

GetCurrentPos() [3/3]

Hep3Vector RecMucTrack::GetCurrentPos ( ) const

inline

Current position.

Definition at line 174 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_CurrentPos / 10; }

GetDepth3() [1/3]

float RecMucTrack::GetDepth3 ( ) const

inline

Length of the track penetrating in iron absorber.

Definition at line 216 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_depth_3 / 10; }

GetDepth3() [2/3]

float RecMucTrack::GetDepth3 ( ) const

inline

Length of the track penetrating in iron absorber.

Definition at line 216 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_depth_3 / 10; }

GetDepth3() [3/3]

float RecMucTrack::GetDepth3 ( ) const

inline

Length of the track penetrating in iron absorber.

Definition at line 216 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_depth_3 / 10; }

getDistHits() [1/3]

vector< float > RecMucTrack::getDistHits ( ) const

inline

Definition at line 275 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits; }

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

getDistHits() [2/3]

vector< float > RecMucTrack::getDistHits ( ) const

inline

Definition at line 275 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits; }

getDistHits() [3/3]

vector< float > RecMucTrack::getDistHits ( ) const

inline

Definition at line 275 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits; }

GetExpectedHits() [1/3]

vector< MucRecHit * > RecMucTrack::GetExpectedHits

(

)

const

Definition at line 1656 of file RecMucTrack.cxx.

{ return m_pExpectedHits; }

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

GetExpectedHits() [2/3]

vector< MucRecHit * > RecMucTrack::GetExpectedHits

(

)

const

GetExpectedHits() [3/3]

vector< MucRecHit * > RecMucTrack::GetExpectedHits

(

)

const

getExpHits() [1/3]

vector< int > RecMucTrack::getExpHits ( ) const

inline

Definition at line 273 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_expHits; }

getExpHits() [2/3]

vector< int > RecMucTrack::getExpHits ( ) const

inline

Definition at line 273 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_expHits; }

getExpHits() [3/3]

vector< int > RecMucTrack::getExpHits ( ) const

inline

Definition at line 273 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_expHits; }

getExtDistHits() [1/3]

vector< float > RecMucTrack::getExtDistHits ( ) const

inline

Definition at line 279 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_ext; }

Referenced by MucRecTrkExt::execute().

getExtDistHits() [2/3]

vector< float > RecMucTrack::getExtDistHits ( ) const

inline

Definition at line 279 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_ext; }

getExtDistHits() [3/3]

vector< float > RecMucTrack::getExtDistHits ( ) const

inline

Definition at line 279 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_ext; }

GetExtMucDistance() [1/3]

Hep3Vector RecMucTrack::GetExtMucDistance ( ) const

inline

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

Definition at line 223 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_Dist_muc_ext; }

Referenced by MucRecTrkExt::execute().

GetExtMucDistance() [2/3]

Hep3Vector RecMucTrack::GetExtMucDistance ( ) const

inline

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

Definition at line 223 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_Dist_muc_ext; }

GetExtMucDistance() [3/3]

Hep3Vector RecMucTrack::GetExtMucDistance ( ) const

inline

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

Definition at line 223 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_Dist_muc_ext; }

getExtMucMomentum() [1/3]

Hep3Vector RecMucTrack::getExtMucMomentum ( ) const

inline

Start momentum of the Ext track in Muc.

Definition at line 160 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucMomentum / 1000; }

getExtMucMomentum() [2/3]

Hep3Vector RecMucTrack::getExtMucMomentum ( ) const

inline

Start momentum of the Ext track in Muc.

Definition at line 160 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucMomentum / 1000; }

getExtMucMomentum() [3/3]

Hep3Vector RecMucTrack::getExtMucMomentum ( ) const

inline

Start momentum of the Ext track in Muc.

Definition at line 160 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucMomentum / 1000; }

getExtMucPos() [1/3]

Hep3Vector RecMucTrack::getExtMucPos ( ) const

inline

start position of the Ext track in Muc.

Definition at line 157 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucPos / 10; }

getExtMucPos() [2/3]

Hep3Vector RecMucTrack::getExtMucPos ( ) const

inline

start position of the Ext track in Muc.

Definition at line 157 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucPos / 10; }

getExtMucPos() [3/3]

Hep3Vector RecMucTrack::getExtMucPos ( ) const

inline

start position of the Ext track in Muc.

Definition at line 157 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtMucPos / 10; }

getExtTrackID() [1/3]

int RecMucTrack::getExtTrackID ( ) const

inline

Where does the trajectory of this track intersect the reference plane?

Where does the trajectory of this track intersect a z=const plane? Transform the Phi, ZR cord. to ZX, ZY cord. The identifier of the Ext track as seed.

Definition at line 154 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtTrackID; }

getExtTrackID() [2/3]

int RecMucTrack::getExtTrackID ( ) const

inline

Where does the trajectory of this track intersect the reference plane?

Where does the trajectory of this track intersect a z=const plane? Transform the Phi, ZR cord. to ZX, ZY cord. The identifier of the Ext track as seed.

Definition at line 154 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtTrackID; }

getExtTrackID() [3/3]

int RecMucTrack::getExtTrackID ( ) const

inline

Where does the trajectory of this track intersect the reference plane?

Where does the trajectory of this track intersect a z=const plane? Transform the Phi, ZR cord. to ZX, ZY cord. The identifier of the Ext track as seed.

Definition at line 154 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_ExtTrackID; }

GetHit() [1/3]

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 1634 of file RecMucTrack.cxx.

                                                                    {
  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;
}

GetHit() [2/3]

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

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

GetHit() [3/3]

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

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

GetHitDistance() [1/3]

float RecMucTrack::GetHitDistance

(

const MucRecHit *

hit

)

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

Definition at line 409 of file RecMucTrack.cxx.

                                                        {
  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;
}

Referenced by MucRecRoadFinder::TrackFinding().

GetHitDistance() [2/3]

float RecMucTrack::GetHitDistance

(

const MucRecHit *

hit

)

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

GetHitDistance() [3/3]

float RecMucTrack::GetHitDistance

(

const MucRecHit *

hit

)

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

GetHitDistance2() [1/3]

float RecMucTrack::GetHitDistance2

(

const MucRecHit *

hit

)

no abs value

Definition at line 437 of file RecMucTrack.cxx.

                                                         {
  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;
}

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

GetHitDistance2() [2/3]

float RecMucTrack::GetHitDistance2

(

const MucRecHit *

hit

)

no abs value

GetHitDistance2() [3/3]

float RecMucTrack::GetHitDistance2

(

const MucRecHit *

hit

)

no abs value

GetHitIndices() [1/3]

vector< long > RecMucTrack::GetHitIndices

(

)

const

Get indices of all hits in the track.

Definition at line 1658 of file RecMucTrack.cxx.

                                              {
  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;
}

GetHitIndices() [2/3]

vector< long > RecMucTrack::GetHitIndices

(

)

const

Get indices of all hits in the track.

GetHitIndices() [3/3]

vector< long > RecMucTrack::GetHitIndices

(

)

const

Get indices of all hits in the track.

GetHitInGap() [1/3]

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

How many hits per gap does this track contain?

Definition at line 1464 of file RecMucTrack.cxx.

                                                                  {
  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;
}

Referenced by ComputeMaxHitsInGap().

GetHitInGap() [2/3]

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

How many hits per gap does this track contain?

GetHitInGap() [3/3]

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

How many hits per gap does this track contain?

GetHitInSeg() [1/3]

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

How many hits does a segment contains.

Definition at line 1496 of file RecMucTrack.cxx.

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

Referenced by FindSegWithMaxHits().

GetHitInSeg() [2/3]

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

How many hits does a segment contains.

GetHitInSeg() [3/3]

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

How many hits does a segment contains.

GetHitInSegOrient() [1/3]

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 1501 of file RecMucTrack.cxx.

                                                                                          {
  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;
}

Referenced by GetHitInSeg().

GetHitInSegOrient() [2/3]

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

How many hits does a segment contains in one orient.

GetHitInSegOrient() [3/3]

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

How many hits does a segment contains in one orient.

GetHits() [1/3]

vector< MucRecHit * > RecMucTrack::GetHits

(

)

const

Get all hits on this track.

Definition at line 1654 of file RecMucTrack.cxx.

{ return m_pHits; }

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

GetHits() [2/3]

vector< MucRecHit * > RecMucTrack::GetHits

(

)

const

Get all hits on this track.

GetHits() [3/3]

vector< MucRecHit * > RecMucTrack::GetHits

(

)

const

Get all hits on this track.

GetMdcExtTrack() [1/3]

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

compute ext track myself from mdc.

Definition at line 278 of file RecMucTrack.cxx.

                                                                 {
 
  // 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;
  ////////////////////
}

Referenced by MucRecTrkExt::execute().

GetMdcExtTrack() [2/3]

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

compute ext track myself from mdc.

GetMdcExtTrack() [3/3]

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

compute ext track myself from mdc.

getMdcMomentum() [1/3]

Hep3Vector RecMucTrack::getMdcMomentum ( ) const

inline

momentum of this track in Mdc

Definition at line 171 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_MdcMomentum / 1000; }

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

getMdcMomentum() [2/3]

Hep3Vector RecMucTrack::getMdcMomentum ( ) const

inline

momentum of this track in Mdc

Definition at line 171 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_MdcMomentum / 1000; }

getMdcMomentum() [3/3]

Hep3Vector RecMucTrack::getMdcMomentum ( ) const

inline

momentum of this track in Mdc

Definition at line 171 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_MdcMomentum / 1000; }

getMucMomentum() [1/3]

Hep3Vector RecMucTrack::getMucMomentum ( ) const

inline

Start momentum of this track in Muc.

Definition at line 168 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_MucMomentum / 1000; }

Referenced by MucRecRoadFinder::execute().

getMucMomentum() [2/3]

Hep3Vector RecMucTrack::getMucMomentum ( ) const

inline

Start momentum of this track in Muc.

Definition at line 168 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_MucMomentum / 1000; }

getMucMomentum() [3/3]

Hep3Vector RecMucTrack::getMucMomentum ( ) const

inline

Start momentum of this track in Muc.

Definition at line 168 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_MucMomentum / 1000; }

getMucPos() [1/3]

Hep3Vector RecMucTrack::getMucPos ( ) const

inline

start position of this track in Muc.

Definition at line 163 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPos / 10; }

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

getMucPos() [2/3]

Hep3Vector RecMucTrack::getMucPos ( ) const

inline

start position of this track in Muc.

Definition at line 163 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPos / 10; }

getMucPos() [3/3]

Hep3Vector RecMucTrack::getMucPos ( ) const

inline

start position of this track in Muc.

Definition at line 163 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPos / 10; }

getMucPosSigma() [1/3]

Hep3Vector RecMucTrack::getMucPosSigma ( ) const

inline

Definition at line 165 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPosSigma / 10; }

Referenced by MucRecTrkExt::execute().

getMucPosSigma() [2/3]

Hep3Vector RecMucTrack::getMucPosSigma ( ) const

inline

Definition at line 165 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPosSigma / 10; }

getMucPosSigma() [3/3]

Hep3Vector RecMucTrack::getMucPosSigma ( ) const

inline

Definition at line 165 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_MucPosSigma / 10; }

GetMucStripPos() [1/3]

Hep3Vector RecMucTrack::GetMucStripPos ( ) const

inline

Definition at line 225 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_StripPhi; }

Referenced by MucRecTrkExt::execute().

GetMucStripPos() [2/3]

Hep3Vector RecMucTrack::GetMucStripPos ( ) const

inline

Definition at line 225 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_StripPhi; }

GetMucStripPos() [3/3]

Hep3Vector RecMucTrack::GetMucStripPos ( ) const

inline

Definition at line 225 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_StripPhi; }

GetNGapsWithHits() [1/3]

int RecMucTrack::GetNGapsWithHits ( ) const

inline

How many gaps provide hits attached to this track?

Definition at line 231 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_numLayers; }

Referenced by LineFit().

GetNGapsWithHits() [2/3]

int RecMucTrack::GetNGapsWithHits ( ) const

inline

How many gaps provide hits attached to this track?

Definition at line 231 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_numLayers; }

GetNGapsWithHits() [3/3]

int RecMucTrack::GetNGapsWithHits ( ) const

inline

How many gaps provide hits attached to this track?

Definition at line 231 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_numLayers; }

GetNSharedHits() [1/3]

int RecMucTrack::GetNSharedHits

(

const RecMucTrack *

track

)

const

How many hits do two tracks share?

Definition at line 1609 of file RecMucTrack.cxx.

                                                                 {
  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;
}

GetNSharedHits() [2/3]

int RecMucTrack::GetNSharedHits

(

const RecMucTrack *

track

)

const

How many hits do two tracks share?

GetNSharedHits() [3/3]

int RecMucTrack::GetNSharedHits

(

const RecMucTrack *

track

)

const

How many hits do two tracks share?

getQuadDistHits() [1/3]

vector< float > RecMucTrack::getQuadDistHits ( ) const

inline

Definition at line 277 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_quad; }

Referenced by MucRecTrkExt::execute().

getQuadDistHits() [2/3]

vector< float > RecMucTrack::getQuadDistHits ( ) const

inline

Definition at line 277 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_quad; }

getQuadDistHits() [3/3]

vector< float > RecMucTrack::getQuadDistHits ( ) const

inline

Definition at line 277 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_distHits_quad; }

GetRecMode() [1/3]

int RecMucTrack::GetRecMode ( ) const

inline

Definition at line 220 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_recmode; }

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

GetRecMode() [2/3]

int RecMucTrack::GetRecMode ( ) const

inline

Definition at line 220 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_recmode; }

GetRecMode() [3/3]

int RecMucTrack::GetRecMode ( ) const

inline

Definition at line 220 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_recmode; }

GetTotalHits() [1/3]

int RecMucTrack::GetTotalHits

(

)

const

How many hits in all does this track contain?

Definition at line 1462 of file RecMucTrack.cxx.

{ return m_pHits.size(); }

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

GetTotalHits() [2/3]

int RecMucTrack::GetTotalHits

(

)

const

How many hits in all does this track contain?

GetTotalHits() [3/3]

int RecMucTrack::GetTotalHits

(

)

const

How many hits in all does this track contain?

getVecHits() [1/3]

vector< int > RecMucTrack::getVecHits ( ) const

inline

Definition at line 269 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ return m_vecHits; }

getVecHits() [2/3]

vector< int > RecMucTrack::getVecHits ( ) const

inline

Definition at line 269 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ return m_vecHits; }

getVecHits() [3/3]

vector< int > RecMucTrack::getVecHits ( ) const

inline

Definition at line 269 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ return m_vecHits; }

HasHit() [1/3]

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

How many hits were attached in the gap with the most attached hits?

Does this track contains an assigned hit?

Definition at line 1570 of file RecMucTrack.cxx.

                                                  {
  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;
}

HasHit() [2/3]

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

How many hits were attached in the gap with the most attached hits?

Does this track contains an assigned hit?

HasHit() [3/3]

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

How many hits were attached in the gap with the most attached hits?

Does this track contains an assigned hit?

HasHitInGap() [1/3]

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

Does this track contain any hits in the given gap?

Definition at line 1588 of file RecMucTrack.cxx.

                                                                   {
  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;
}

Referenced by ComputeDepth().

HasHitInGap() [2/3]

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

Does this track contain any hits in the given gap?

HasHitInGap() [3/3]

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

Does this track contain any hits in the given gap?

IsInsideSuperConductor() [1/3]

bool RecMucTrack::IsInsideSuperConductor

(

Hep3Vector

pos

)

Definition at line 318 of file RecMucTrack.cxx.

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

Referenced by GetMdcExtTrack().

IsInsideSuperConductor() [2/3]

bool RecMucTrack::IsInsideSuperConductor

(

Hep3Vector

pos

)

IsInsideSuperConductor() [3/3]

bool RecMucTrack::IsInsideSuperConductor

(

Hep3Vector

pos

)

LineFit() [1/3]

void RecMucTrack::LineFit

(

int

fittingMethod

)

Line fit with hits on a seg with max hits.

Definition at line 1007 of file RecMucTrack.cxx.

                                             {
  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 );
      }
    }
    //////////////////////////////////////
 
    ///////////////find expected strips for this 3D road
    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;
}

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

LineFit() [2/3]

void RecMucTrack::LineFit

(

int

fittingMethod

)

Line fit with hits on a seg with max hits.

LineFit() [3/3]

void RecMucTrack::LineFit

(

int

fittingMethod

)

Line fit with hits on a seg with max hits.

operator=() [1/6]

RecMucTrack & RecMucTrack::operator=

(

const DstMucTrack &

dstTrack

)

Assignment constructor from DstMucTrack.

Definition at line 181 of file RecMucTrack.cxx.

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

operator=() [2/6]

RecMucTrack & RecMucTrack::operator=

(

const DstMucTrack &

dstTrack

)

Assignment constructor from DstMucTrack.

operator=() [3/6]

RecMucTrack & RecMucTrack::operator=

(

const DstMucTrack &

dstTrack

)

Assignment constructor from DstMucTrack.

operator=() [4/6]

RecMucTrack & RecMucTrack::operator=

(

const RecMucTrack &

orig

)

Assignment constructor.

Definition at line 84 of file RecMucTrack.cxx.

                                                           {
  // 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;
}

operator=() [5/6]

RecMucTrack & RecMucTrack::operator=

(

const RecMucTrack &

orig

)

Assignment constructor.

operator=() [6/6]

RecMucTrack & RecMucTrack::operator=

(

const RecMucTrack &

orig

)

Assignment constructor.

OutputUnitChange() [1/3]

void RecMucTrack::OutputUnitChange

(

)

change unit

Definition at line 1709 of file RecMucTrack.cxx.

                                   {
  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;
}

Referenced by ComputeTrackInfo().

OutputUnitChange() [2/3]

void RecMucTrack::OutputUnitChange

(

)

change unit

OutputUnitChange() [3/3]

void RecMucTrack::OutputUnitChange

(

)

change unit

PrintHitsInfo() [1/3]

void RecMucTrack::PrintHitsInfo

(

)

const

Print Hits Infomation.

Definition at line 1691 of file RecMucTrack.cxx.

                                      {
  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;
    }
  }
}

PrintHitsInfo() [2/3]

void RecMucTrack::PrintHitsInfo

(

)

const

Print Hits Infomation.

PrintHitsInfo() [3/3]

void RecMucTrack::PrintHitsInfo

(

)

const

Print Hits Infomation.

Project() [1/3]

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 354 of file RecMucTrack.cxx.

                                      {
  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; }
      }
    }
  }
}

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

Project() [2/3]

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?

Project() [3/3]

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?

pushExtDistHits() [1/3]

void RecMucTrack::pushExtDistHits ( float dist )

inline

Definition at line 311 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext.push_back( dist ); }

Referenced by MucRecTrkExt::execute().

pushExtDistHits() [2/3]

void RecMucTrack::pushExtDistHits ( float dist )

inline

Definition at line 311 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext.push_back( dist ); }

pushExtDistHits() [3/3]

void RecMucTrack::pushExtDistHits ( float dist )

inline

Definition at line 311 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext.push_back( dist ); }

SetCurrentDir() [1/3]

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

set current direction of the trajectory.

Definition at line 246 of file RecMucTrack.cxx.

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

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

SetCurrentDir() [2/3]

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

set current direction of the trajectory.

SetCurrentDir() [3/3]

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

set current direction of the trajectory.

SetCurrentInsct() [1/3]

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

set current intersection of the trajectory with strip plane.

Definition at line 250 of file RecMucTrack.cxx.

                                                                               {
  m_CurrentInsct.set( x, y,
                      z ); /////this is intenal function, it receive "mm", so need not convert.
}

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

SetCurrentInsct() [2/3]

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

set current intersection of the trajectory with strip plane.

SetCurrentInsct() [3/3]

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

set current intersection of the trajectory with strip plane.

SetCurrentPos() [1/3]

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

set current position of the trajectory.

Definition at line 242 of file RecMucTrack.cxx.

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

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

SetCurrentPos() [2/3]

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

set current position of the trajectory.

SetCurrentPos() [3/3]

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

set current position of the trajectory.

SetDefault() [1/3]

void RecMucTrack::SetDefault

(

)

Definition at line 187 of file RecMucTrack.cxx.

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

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

SetDefault() [2/3]

void RecMucTrack::SetDefault

(

)

SetDefault() [3/3]

void RecMucTrack::SetDefault

(

)

setDistHits() [1/3]

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

inline

Definition at line 305 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_distHits = disthits; }

Referenced by RecMucTrackCnv::TObjectToDataObject().

setDistHits() [2/3]

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

inline

Definition at line 305 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_distHits = disthits; }

setDistHits() [3/3]

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

inline

Definition at line 305 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_distHits = disthits; }

setExpHits() [1/6]

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

inline

Definition at line 303 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_expHits = exphits; }

setExpHits() [2/6]

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

inline

Definition at line 303 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_expHits = exphits; }

setExpHits() [3/6]

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

inline

Definition at line 303 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_expHits = exphits; }

setExpHits() [4/6]

void RecMucTrack::setExpHits

(

vector< MucRecHit * > &

pHits

)

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

setExpHits() [5/6]

void RecMucTrack::setExpHits

(

vector< MucRecHit * > &

pHits

)

setExpHits() [6/6]

void RecMucTrack::setExpHits

(

vector< MucRecHit * > &

pHits

)

Definition at line 273 of file RecMucTrack.cxx.

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

setExtDistHits() [1/3]

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

inline

Definition at line 309 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext = disthits; }

setExtDistHits() [2/3]

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

inline

Definition at line 309 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext = disthits; }

setExtDistHits() [3/3]

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

inline

Definition at line 309 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_ext = disthits; }

SetExtMucMomentum() [1/3]

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

set start moment of ext track in Muc.

Definition at line 238 of file RecMucTrack.cxx.

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

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

SetExtMucMomentum() [2/3]

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

set start moment of ext track in Muc.

SetExtMucMomentum() [3/3]

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

set start moment of ext track in Muc.

SetExtMucPos() [1/3]

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 234 of file RecMucTrack.cxx.

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

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

SetExtMucPos() [2/3]

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)

SetExtMucPos() [3/3]

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)

SetExtTrack() [1/3]

void RecMucTrack::SetExtTrack

(

RecExtTrack *

extTrack

)

set Ext track point.

Definition at line 261 of file RecMucTrack.cxx.

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

Referenced by MucRecTrkExt::execute().

SetExtTrack() [2/3]

void RecMucTrack::SetExtTrack

(

RecExtTrack *

extTrack

)

set Ext track point.

SetExtTrack() [3/3]

void RecMucTrack::SetExtTrack

(

RecExtTrack *

extTrack

)

set Ext track point.

SetExtTrackID() [1/3]

void RecMucTrack::SetExtTrackID ( int id )

inline

set Ext track id. for compute from mdc myself

Definition at line 103 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_ExtTrackID = id; }

Referenced by MucRecTrkExt::execute().

SetExtTrackID() [2/3]

void RecMucTrack::SetExtTrackID ( int id )

inline

set Ext track id. for compute from mdc myself

Definition at line 103 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_ExtTrackID = id; }

SetExtTrackID() [3/3]

void RecMucTrack::SetExtTrackID ( int id )

inline

set Ext track id. for compute from mdc myself

Definition at line 103 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_ExtTrackID = id; }

SetMdcMomentum() [1/3]

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

set start moment of the track in Mdc.

Definition at line 209 of file RecMucTrack.cxx.

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

Referenced by MucRecTrkExt::execute().

SetMdcMomentum() [2/3]

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

set start moment of the track in Mdc.

SetMdcMomentum() [3/3]

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

set start moment of the track in Mdc.

SetMdcPos() [1/3]

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

set start position of the track in Mdc.

Definition at line 205 of file RecMucTrack.cxx.

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

Referenced by MucRecTrkExt::execute().

SetMdcPos() [2/3]

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

set start position of the track in Mdc.

SetMdcPos() [3/3]

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

set start position of the track in Mdc.

SetMucMomentum() [1/3]

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

set start moment of the track in Muc.

Definition at line 227 of file RecMucTrack.cxx.

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

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

SetMucMomentum() [2/3]

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

set start moment of the track in Muc.

SetMucMomentum() [3/3]

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

set start moment of the track in Muc.

SetMucPos() [1/3]

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 213 of file RecMucTrack.cxx.

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

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

SetMucPos() [2/3]

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)

SetMucPos() [3/3]

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)

SetMucPosSigma() [1/3]

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

Definition at line 220 of file RecMucTrack.cxx.

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

SetMucPosSigma() [2/3]

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

SetMucPosSigma() [3/3]

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

setQuadDistHits() [1/3]

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

inline

Definition at line 307 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_quad = disthits; }

setQuadDistHits() [2/3]

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

inline

Definition at line 307 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_quad = disthits; }

setQuadDistHits() [3/3]

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

inline

Definition at line 307 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_distHits_quad = disthits; }

SetRecMode() [1/3]

void RecMucTrack::SetRecMode ( int recmode )

inline

Definition at line 218 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_recmode = recmode; }

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

SetRecMode() [2/3]

void RecMucTrack::SetRecMode ( int recmode )

inline

Definition at line 218 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_recmode = recmode; }

SetRecMode() [3/3]

void RecMucTrack::SetRecMode ( int recmode )

inline

Definition at line 218 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_recmode = recmode; }

setTrackId() [1/3]

void RecMucTrack::setTrackId

(

const int

trackId

)

set the index for this track.

Definition at line 199 of file RecMucTrack.cxx.

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

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

setTrackId() [2/3]

void RecMucTrack::setTrackId

(

const int

trackId

)

set the index for this track.

setTrackId() [3/3]

void RecMucTrack::setTrackId

(

const int

trackId

)

set the index for this track.

setVecHits() [1/6]

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

inline

Definition at line 301 of file InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h.

{ m_vecHits = vechits; }

setVecHits() [2/6]

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

inline

Definition at line 301 of file InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h.

{ m_vecHits = vechits; }

setVecHits() [3/6]

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

inline

Definition at line 301 of file Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h.

{ m_vecHits = vechits; }

setVecHits() [4/6]

void RecMucTrack::setVecHits

(

vector< MucRecHit * > &

pHits

)

reload setVecHits

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

setVecHits() [5/6]

void RecMucTrack::setVecHits

(

vector< MucRecHit * > &

pHits

)

reload setVecHits

setVecHits() [6/6]

void RecMucTrack::setVecHits

(

vector< MucRecHit * > &

pHits

)

reload setVecHits

Definition at line 268 of file RecMucTrack.cxx.

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

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/MucRecEvent/RecMucTrack.h
• InstallArea/x86_64-el9-gcc13-opt/include/MucRecEvent/RecMucTrack.h
• Muc/MucRecEvent/include/MucRecEvent/RecMucTrack.h
• RecMucTrack.cxx