Contains the straight part of a track in one region only. More...

#include <QwPartialTrack.h>

Inheritance diagram for QwPartialTrack:

Collaboration diagram for QwPartialTrack:

Public Member Functions
	QwPartialTrack ()
	Default constructor. More...

	QwPartialTrack (const QwPartialTrack &that)
	Copy constructor by reference. More...

	QwPartialTrack (const QwPartialTrack *that)
	Copy constructor from pointer. More...

	QwPartialTrack (const TVector3 &position, const TVector3 &momentum)
	Constructor with track position and direction. More...

virtual	~QwPartialTrack ()
	Destructor. More...

QwPartialTrack &	operator= (const QwPartialTrack &that)
	Assignment operator. More...

bool	IsVoid () const

bool	IsValid () const

bool	IsGood () const

bool	IsNotValid () const

bool	IsUsed () const

bool	IsNotUsed () const

void	Clear (Option_t *option="")

void	Reset (Option_t *option="")

double	GetChiWeight () const

void	Print (const Option_t *options=0) const

void	PrintValid ()

const TVector3	GetPosition (const double z) const
	Return the vertex at position z. More...

const TVector3	GetMomentumDirection () const
	Return the direction. More...

Double_t	GetMomentumDirectionPhi () const
	Return the phi angle. More...

Double_t	GetMomentumDirectionTheta () const
	Return the theta angle. More...

Double_t	GetVertexZ () const
	Get the vertex position. More...

QwPartialTrack &	SmearPosition (const double sigma_x, const double sigma_y)
	Smear the position. More...

QwPartialTrack &	SmearAngleTheta (const double sigma)
	Smear the theta angle. More...

QwPartialTrack &	SmearAnglePhi (const double sigma)
	Smear the phi angle. More...

const QwVertex *	DeterminePositionInDetector (const QwDetectorInfo *geometry)
	Determine vertex in detector. More...

const QwVertex *	DeterminePositionInTarget (const QwGeometry &geometry)
	Determine vertex in the target. More...

const QwVertex *	DeterminePositionInTriggerScintillators (const QwGeometry &geometry)
	Determine intersection with trigger scintillators. More...

const QwVertex *	DeterminePositionInCerenkovBars (const QwGeometry &geometry)
	Determine intersection with cerenkov bars. More...

const QwVertex *	DeterminePositionInHDC (const QwGeometry &geometry)
	Determine position in first horizontal drift chamber. More...

double	GetAverageResidual () const

void	SetAverageResidual (const double residual)

double	CalculateAverageResidual ()

void	SetAverageResidual ()

void	SetAlone (const int alone)

int	GetAlone () const
	Only 2 Treelines in Plane 0. More...

void	RotateCoordinates ()
	Rotate coordinates to right octant. More...

void	RotateRotator (const QwDetectorInfo *geometry)

Creating, adding, and getting tree lines
QwTreeLine *	CreateNewTreeLine ()
	Create a new tree line. More...

void	AddTreeLine (const QwTreeLine *treeline)
	Add an existing tree line as a copy. More...

void	AddTreeLineList (const QwTreeLine *treelinelist)
	Add a linked list of existing tree lines as a copy. More...

void	AddTreeLineList (const std::vector< QwTreeLine * > &treelinelist)
	Add a standard vector of existing tree lines as a copy. More...

void	ClearTreeLines (Option_t *option="")
	Clear the list of tree lines. More...

void	ResetTreeLines (Option_t *option="")
	Reset the list of tree lines. More...

Int_t	GetNumberOfTreeLines () const
	Get the number of tree lines. More...

const std::vector< QwTreeLine * > &	GetListOfTreeLines () const
	Get the list of tree lines. More...

const QwTreeLine *	GetTreeLine (const int tl) const
	Get the specified tree line. More...

void	PrintTreeLines (Option_t *option="") const
	Print the list of tree lines. More...

Public Member Functions inherited from VQwTrackingElement
	VQwTrackingElement ()
	Default constructor. More...

	VQwTrackingElement (const VQwTrackingElement &that)

virtual	~VQwTrackingElement ()
	Virtual destructor. More...

VQwTrackingElement &	operator= (const VQwTrackingElement &that)
	Assignment operator. More...

const QwDetectorInfo *	GetDetectorInfo () const
	Get the detector info pointer. More...

void	SetDetectorInfo (const QwDetectorInfo *detectorinfo)
	Set the detector info pointer. More...

EQwRegionID	GetRegion () const
	Get the region. More...

void	SetRegion (EQwRegionID region)
	Set the region. More...

EQwDetectorPackage	GetPackage () const
	Get the package. More...

void	SetPackage (EQwDetectorPackage package)
	Set the package. More...

int	GetOctant () const
	Get the octant number. More...

void	SetOctant (int octant)
	Set the octant number. More...

EQwDirectionID	GetDirection () const
	Get the direction. More...

void	SetDirection (EQwDirectionID direction)
	Set the direction. More...

int	GetPlane () const
	Get the plane number. More...

void	SetPlane (int plane)
	Set the plane number. More...

int	GetElement () const
	Get the element number. More...

void	SetElement (int element)
	Set the element number. More...

void	SetGeometryTo (const VQwTrackingElement &e)
	Copy the geometry info from another object. More...

Public Member Functions inherited from QwObjectCounter< QwPartialTrack >
	QwObjectCounter ()
	Default constructor. More...

	QwObjectCounter (const QwObjectCounter &)
	Copy constructor. More...

virtual	~QwObjectCounter ()
	Destructor. More...

Data Fields
Double_t	fOffsetX
	x coordinate (at MAGNET_CENTER) More...

Double_t	fOffsetY
	y coordinate (at MAGNET_CENTER) More...

Double_t	fSlopeX
	x slope More...

Double_t	fSlopeY
	y slope More...

Double_t	fChi
	combined chi square More...

Double_t	fCov [4][4]
	covariance matrix More...

Double_t	TResidual [kNumDirections]

Double_t	fSignedResidual [12]

Bool_t	fIsUsed
	used (part of a track) More...

Bool_t	fIsVoid
	marked as being void More...

Bool_t	fIsGood

Int_t	fNumMiss
	missing hits More...

Int_t	fNumHits
	used hits More...

double	pR2hit [3]
	x-y-z position at R2 More...

double	uvR2hit [3]
	direction at R2 More...

double	pR3hit [3]
	x-y-z position at R3 More...

double	uvR3hit [3]
	direction at R3 More...

int	fAlone

Private Member Functions
void	Initialize ()
	Initialization. More...

	ClassDef (QwPartialTrack, 1)

Private Attributes
Int_t	fNQwTreeLines
	Number of tree lines in this partial track. More...

std::vector< QwTreeLine * >	fQwTreeLines
	List of tree lines in this partial track. More...

double	fAverageResidual
	number of Plane 0 Treelines More...

Friends
std::ostream &	operator<< (std::ostream &stream, const QwPartialTrack &pt)

Additional Inherited Members
Static Public Member Functions inherited from QwObjectCounter< QwPartialTrack >
static size_t	GetObjectsCreated ()
	Get number of objects ever created. More...

static size_t	GetObjectsAlive ()
	Get number of objects still alive. More...

Protected Member Functions inherited from VQwTrackingElement
	ClassDef (VQwTrackingElement, 1)

Protected Attributes inherited from VQwTrackingElement
const QwDetectorInfo *	fDetectorInfo

EQwRegionID	fRegion
	///< Detector info pointer More...

EQwDetectorPackage	fPackage
	Package. More...

int	fOctant
	Octant number. More...

EQwDirectionID	fDirection
	Direction. More...

int	fPlane
	Plane number. More...

int	fElement
	Element number. More...

Detailed Description

Contains the straight part of a track in one region only.

A QwPartialTrack contains tracking information for one region only. It is constructed by combining multiple tree lines together, and in turn used to construct the complete track.

Definition at line 38 of file QwPartialTrack.h.

Constructor & Destructor Documentation

QwPartialTrack::QwPartialTrack ( )

Default constructor.

Default constructor

Definition at line 15 of file QwPartialTrack.cc.

References Initialize().

 {
   // Initialize
   Initialize();
 }

Here is the call graph for this function:

QwPartialTrack::QwPartialTrack ( const QwPartialTrack & that )

Copy constructor by reference.

Copy constructor by reference

Definition at line 43 of file QwPartialTrack.cc.

References Initialize().

 : VQwTrackingElement(that)
 {
   // Initialize
   Initialize();
 
   // Assignment
   *this = that;
 }

Here is the call graph for this function:

QwPartialTrack::QwPartialTrack ( const QwPartialTrack * that )

Copy constructor from pointer.

Copy constructor with pointer

Definition at line 56 of file QwPartialTrack.cc.

References Initialize().

 : VQwTrackingElement(*that)
 {
   // Initialize
   Initialize();
 
   // Null pointer          
   if (that == 0) return;
 
   // Assignment
   *this = *that;
 }

Here is the call graph for this function:

QwPartialTrack::QwPartialTrack	(	const TVector3 &	position,
		const TVector3 &	direction
	)

Constructor with track position and direction.

Constructor with position and direction vectors

Parameters

position	Position of the partial track
direction	Direction of the partial track

Definition at line 26 of file QwPartialTrack.cc.

References fOffsetX, fOffsetY, fSlopeX, fSlopeY, and Initialize().

 {
   // Initialize
   Initialize();
 
   // Calculate slopes
   fSlopeX = direction.X() / direction.Z();
   fSlopeY = direction.Y() / direction.Z();
 
   // Calculate offset
   fOffsetX = position.X() - fSlopeX * position.Z();
   fOffsetY = position.Y() - fSlopeY * position.Z();
 }

Here is the call graph for this function:

QwPartialTrack::~QwPartialTrack ( )

virtual

Destructor.

Destructor

Definition at line 72 of file QwPartialTrack.cc.

References Clear().

 {
   // Clear the partial track
   Clear();
 
   /* for (int i = 0; i < kNumDirections; i++){
     if(fTreeLine[i])
       delete fTreeLine[i];
   }
   */
   // nothing  
 }

Here is the call graph for this function:

Member Function Documentation

void QwPartialTrack::AddTreeLine ( const QwTreeLine * treeline )

Add an existing tree line as a copy.

Definition at line 225 of file QwPartialTrack.cc.

References fNQwTreeLines, and fQwTreeLines.

Referenced by AddTreeLineList(), CreateNewTreeLine(), QwTrackingTreeCombine::r3_PartialTrackFit(), and QwTrackingTreeCombine::TcTreeLineCombine().

 {
   fQwTreeLines.push_back(new QwTreeLine(treeline));
   ++fNQwTreeLines;
 }

Here is the caller graph for this function:

void QwPartialTrack::AddTreeLineList ( const QwTreeLine * treelinelist )

Add a linked list of existing tree lines as a copy.

Definition at line 232 of file QwPartialTrack.cc.

References AddTreeLine(), and QwTreeLine::next.

Referenced by operator=().

 {
   for (const QwTreeLine *treeline = treelinelist;
          treeline; treeline = treeline->next){
     if (treeline->IsValid()){
        AddTreeLine(treeline);
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwPartialTrack::AddTreeLineList ( const std::vector< QwTreeLine * > & treelinelist )

Add a standard vector of existing tree lines as a copy.

Definition at line 243 of file QwPartialTrack.cc.

References AddTreeLine().

 { 
   for (std::vector<QwTreeLine*>::const_iterator treeline = treelinelist.begin();
        treeline != treelinelist.end(); treeline++)
     AddTreeLine(*treeline);
 }

Here is the call graph for this function:

double QwPartialTrack::CalculateAverageResidual ( )

Calculate the average residual of this partial track over all treelines

Definition at line 188 of file QwPartialTrack.cc.

References fAverageResidual, QwTreeLine::GetAverageResidual(), GetNumberOfTreeLines(), GetTreeLine(), and QwTreeLine::IsUsed().

Referenced by SetAverageResidual(), and QwTrackingTreeCombine::TlTreeCombine().

 {
   int numTreeLines = 0;
   double sumResiduals = 0.0;
   for (Int_t i = 0; i < GetNumberOfTreeLines(); i++) {
     const QwTreeLine* treeline = GetTreeLine(i);
     if (treeline->IsUsed()) {
       numTreeLines++;
       sumResiduals += treeline->GetAverageResidual();
     }
   } // end of loop over treelines
   fAverageResidual = sumResiduals / numTreeLines;
   return fAverageResidual;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

QwPartialTrack::ClassDef	(	QwPartialTrack	,
		1
	)

private

void QwPartialTrack::Clear ( Option_t * option = "" )

Definition at line 205 of file QwPartialTrack.cc.

References ClearTreeLines().

Referenced by ~QwPartialTrack().

 {
   ClearTreeLines(option);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwPartialTrack::ClearTreeLines ( Option_t * option = "" )

Clear the list of tree lines.

Definition at line 251 of file QwPartialTrack.cc.

References fNQwTreeLines, and fQwTreeLines.

Referenced by Clear(), operator=(), and ResetTreeLines().

 {
   for (size_t i = 0; i < fQwTreeLines.size(); i++) {
     QwTreeLine* tl = fQwTreeLines.at(i);
     delete tl;
   }
   fQwTreeLines.clear();
   fNQwTreeLines = 0;
 }

Here is the caller graph for this function:

QwTreeLine * QwPartialTrack::CreateNewTreeLine ( )

Create a new tree line.

Definition at line 217 of file QwPartialTrack.cc.

References AddTreeLine().

 {
   QwTreeLine* treeline = new QwTreeLine();
   AddTreeLine(treeline);
   return treeline;
 }

Here is the call graph for this function:

const QwVertex * QwPartialTrack::DeterminePositionInCerenkovBars ( const QwGeometry & geometry )

Determine intersection with cerenkov bars.

This method checks whether the region 3 partial track has an intersection with the trigger scintillators in a detector package.

The cerenkov bars are assumed to be the only (or at least the first) detectors in region 3.

Definition at line 458 of file QwPartialTrack.cc.

References DeterminePositionInDetector(), VQwTrackingElement::GetPackage(), QwGeometry::in(), and kRegionIDCer.

 {
   // Get the Cherenkov detectors in this package
   QwGeometry cerenkov(geometry.in(kRegionIDCer).in(GetPackage()));
 
   // Find the first hit
   for (size_t i = 0; i < cerenkov.size(); i++) {
     const QwVertex* vertex = DeterminePositionInDetector(cerenkov.at(i));
     if (vertex) return vertex;
   }
   return 0;
 }

Here is the call graph for this function:

const QwVertex * QwPartialTrack::DeterminePositionInDetector ( const QwDetectorInfo * detector )

Determine vertex in detector.

This method checks determines the intersection of the partial track with a planar detector.

Definition at line 395 of file QwPartialTrack.cc.

References QwLog::endl(), QwDetectorInfo::GetActiveWidthX(), QwDetectorInfo::GetActiveWidthY(), QwDetectorInfo::GetDetectorPitch(), QwDetectorInfo::GetPosition(), GetPosition(), QwDetectorInfo::GetZPosition(), and QwMessage.

Referenced by DeterminePositionInCerenkovBars(), DeterminePositionInHDC(), and DeterminePositionInTriggerScintillators().

 {
   QwVertex* vertex = 0;
   if (detector->GetDetectorPitch() != 0) {
     // Get intersection with z-plane
     double z = detector->GetZPosition();
     TVector3 position = GetPosition(z);
     // Check acceptance with active width
     if (fabs(position.X() - detector->GetPosition().X()) < detector->GetActiveWidthX()/2
      && fabs(position.Y() - detector->GetPosition().Y()) < detector->GetActiveWidthY()/2) {
       vertex = new QwVertex(position);
     } else {
       QwMessage << "No support for intersections of partial tracks "
           << "with tilted detectors yet..." << QwLog::endl;
     }
   }
   return vertex;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const QwVertex * QwPartialTrack::DeterminePositionInHDC ( const QwGeometry & geometry )

Determine position in first horizontal drift chamber.

This method checks whether the region 2 partial track has an intersection with the horizontal drift chambers in a detector package.

The drift chambers are assumed to be the only (or at least the first) detectors in region 2.

Definition at line 478 of file QwPartialTrack.cc.

References DeterminePositionInDetector(), VQwTrackingElement::GetPackage(), QwGeometry::in(), and kRegionID2.

 {
   // Get the HDC detectors in this package
   QwGeometry hdc(geometry.in(kRegionID2).in(GetPackage()));
   // Find the first hit
   for (size_t i = 0; i < hdc.size(); i++) {
     const QwVertex* vertex = DeterminePositionInDetector(hdc.at(i));
     if (vertex) return vertex;
   }
   return 0;
 }

Here is the call graph for this function:

const QwVertex * QwPartialTrack::DeterminePositionInTarget ( const QwGeometry & geometry )

Determine vertex in the target.

This method checks determines the intersection of the partial track with the target.

Note: At this point the target defined at one z value only. A more careful determination of the primary vertex will be implemented using the general QwVertex class, e.g. QwVertex* primary_vertex = QwVertex (partial_track, beam)

Definition at line 423 of file QwPartialTrack.cc.

References QwLog::endl(), GetPosition(), and QwVerbose.

 {
   // TODO target should be put in in the geometry file
   TVector3 primary = GetPosition(0.0);
   QwVerbose << "Target vertex at : (" << primary.X() << "," << primary.Y() << "," << primary.Z() << ")" << QwLog::endl;
   return 0;
 }

Here is the call graph for this function:

const QwVertex * QwPartialTrack::DeterminePositionInTriggerScintillators ( const QwGeometry & geometry )

Determine intersection with trigger scintillators.

This method checks whether the region 3 partial track has an intersection with the trigger scintillators in a detector package.

The trigger scintillators are assumed to be the only (or at least the first) detectors in region 3.

Definition at line 438 of file QwPartialTrack.cc.

References DeterminePositionInDetector(), VQwTrackingElement::GetPackage(), QwGeometry::in(), and kRegionIDTrig.

 {
   // Get the trigger scintillator detectors in this package
   QwGeometry trigscint(geometry.in(kRegionIDTrig).in(GetPackage()));
 
   // Find the first hit
   for (size_t i = 0; i < trigscint.size(); i++) {
     const QwVertex* vertex = DeterminePositionInDetector(trigscint.at(i));
     if (vertex) return vertex;
   }
   return 0;
 }

Here is the call graph for this function:

int QwPartialTrack::GetAlone ( ) const

inline

Only 2 Treelines in Plane 0.

Definition at line 161 of file QwPartialTrack.h.

References fAlone.

161 { return fAlone; }

QwPartialTrack::fAlone

int fAlone

Definition: QwPartialTrack.h:197

double QwPartialTrack::GetAverageResidual ( ) const

inline

Definition at line 153 of file QwPartialTrack.h.

References fAverageResidual.

153 { return fAverageResidual; };

QwPartialTrack::fAverageResidual

double fAverageResidual

number of Plane 0 Treelines

Definition: QwPartialTrack.h:202

double QwPartialTrack::GetChiWeight ( ) const

Determine the chi^2 for a partial track, weighted by the number of hits

Returns: Weighted chi^2

Definition at line 167 of file QwPartialTrack.cc.

References QwLog::endl(), fChi, fNumHits, fNumMiss, and QwDebug.

Referenced by QwTrackingTreeSort::rcPartConnSort().

 {
   // Determine the weight if there are enough hits
   if (fNumHits >= fNumMiss) {
     // NOTE Added +1 to get this to work if fNumHits == fNumMiss (wdc)
     double weight = (double) (fNumHits + fNumMiss + 1)
                            / (fNumHits - fNumMiss + 1);
     return weight * weight * fChi;
     // TODO Why is the weight squared here, but not in the weighted chi^2 for treelines?
 
   } else {
 
     QwDebug << "miss = " << fNumMiss << ", hit = " << fNumHits << QwLog::endl;
     return 100.0; // This is bad...
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const std::vector<QwTreeLine*>& QwPartialTrack::GetListOfTreeLines ( ) const

inline

Get the list of tree lines.

Definition at line 99 of file QwPartialTrack.h.

References fQwTreeLines.

Referenced by QwTrackingTreeSort::rcPTCommonWires().

99 { return fQwTreeLines; };

QwPartialTrack::fQwTreeLines

std::vector< QwTreeLine * > fQwTreeLines

List of tree lines in this partial track.

Definition: QwPartialTrack.h:45

Here is the caller graph for this function:

const TVector3 QwPartialTrack::GetMomentumDirection ( ) const

Return the direction.

Determines the direction of the track at the given z position

Definition at line 340 of file QwPartialTrack.cc.

References fSlopeX, and fSlopeY.

Referenced by QwMatrixLookup::Bridge(), QwRayTracer::Bridge(), QwEvent::CalculateKinematics(), QwBridgingTrackFilter::Filter(), GetMomentumDirectionPhi(), and GetMomentumDirectionTheta().

 {
   TVector3 direction;
   double kz = sqrt(fSlopeX * fSlopeX + fSlopeY * fSlopeY + 1);
   direction.SetX(fSlopeX / kz);
   direction.SetY(fSlopeY / kz);
   direction.SetZ(1 / kz);
   return direction;
 }

Here is the caller graph for this function:

Double_t QwPartialTrack::GetMomentumDirectionPhi ( ) const

inline

Return the phi angle.

Definition at line 118 of file QwPartialTrack.h.

References GetMomentumDirection().

Referenced by QwBridgingTrackFilter::Filter(), SmearAnglePhi(), and SmearAngleTheta().

                                              {
       return GetMomentumDirection().Phi();
     };

Here is the call graph for this function:

Here is the caller graph for this function:

Double_t QwPartialTrack::GetMomentumDirectionTheta ( ) const

inline

Return the theta angle.

Definition at line 122 of file QwPartialTrack.h.

References GetMomentumDirection().

Referenced by QwEvent::CalculateKinematics(), QwBridgingTrackFilter::Filter(), SmearAnglePhi(), and SmearAngleTheta().

                                                {
       return GetMomentumDirection().Theta();
     };

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwPartialTrack::GetNumberOfTreeLines ( ) const

inline

Get the number of tree lines.

Definition at line 97 of file QwPartialTrack.h.

References fNQwTreeLines.

Referenced by CalculateAverageResidual().

97 { return fNQwTreeLines; };

QwPartialTrack::fNQwTreeLines

Int_t fNQwTreeLines

Number of tree lines in this partial track.

Definition: QwPartialTrack.h:43

Here is the caller graph for this function:

const TVector3 QwPartialTrack::GetPosition ( const double z ) const

Return the vertex at position z.

Determines the position of the track at the given z position

Definition at line 328 of file QwPartialTrack.cc.

References fOffsetX, fOffsetY, fSlopeX, and fSlopeY.

Referenced by QwMatrixLookup::Bridge(), QwRayTracer::Bridge(), QwEvent::CalculateKinematics(), DeterminePositionInDetector(), DeterminePositionInTarget(), and QwBridgingTrackFilter::Filter().

 {
   TVector3 position;
   position.SetX(fOffsetX + fSlopeX * z);
   position.SetY(fOffsetY + fSlopeY * z);
   position.SetZ(z);
   return position;
 }

Here is the caller graph for this function:

const QwTreeLine* QwPartialTrack::GetTreeLine ( const int tl ) const

inline

Get the specified tree line.

Definition at line 101 of file QwPartialTrack.h.

References fQwTreeLines.

Referenced by CalculateAverageResidual(), and DrawTreeLineVDC().

101 { return fQwTreeLines.at(tl); };

QwPartialTrack::fQwTreeLines

std::vector< QwTreeLine * > fQwTreeLines

List of tree lines in this partial track.

Definition: QwPartialTrack.h:45

Here is the caller graph for this function:

Double_t QwPartialTrack::GetVertexZ ( ) const

inline

Get the vertex position.

Definition at line 128 of file QwPartialTrack.h.

References fOffsetX, fOffsetY, fSlopeX, and fSlopeY.

Referenced by QwEvent::CalculateKinematics().

                                 {
       return
          -(fSlopeX * fOffsetX + fSlopeY * fOffsetY) /
           (fSlopeX * fSlopeX  + fSlopeY * fSlopeY);
     }

Here is the caller graph for this function:

void QwPartialTrack::Initialize ( )

private

Initialization.

Perform object initialization

Definition at line 142 of file QwPartialTrack.cc.

References fAlone, fIsGood, fIsUsed, fIsVoid, fOffsetX, fOffsetY, fSignedResidual, fSlopeX, fSlopeY, kNumDirections, and TResidual.

Referenced by QwPartialTrack().

 {
   // Initialize the member fields
   fOffsetX = 0.0; fOffsetY = 0.0;
   fSlopeX = 0.0;  fSlopeY = 0.0;
   fIsVoid = false; fIsUsed = false; fIsGood = false;
 
   // Initialize pointers
   for (int i = 0; i < kNumDirections; ++i){
     TResidual[i]=0.0;
   }
 
   for(int i=0; i< 12;++i)
     fSignedResidual[i]=-10;
 
   //Only 2 Plane 0 treelines
   fAlone=0;
 }

Here is the caller graph for this function:

bool QwPartialTrack::IsGood ( ) const

inline

Definition at line 73 of file QwPartialTrack.h.

References fIsGood.

73 { return fIsGood; };

QwPartialTrack::fIsGood

Bool_t fIsGood

Definition: QwPartialTrack.h:186

bool QwPartialTrack::IsNotUsed ( ) const

inline

Definition at line 76 of file QwPartialTrack.h.

References fIsUsed.

76 { return ! fIsUsed; };

QwPartialTrack::fIsUsed

Bool_t fIsUsed

used (part of a track)

Definition: QwPartialTrack.h:184

bool QwPartialTrack::IsNotValid ( ) const

inline

Definition at line 74 of file QwPartialTrack.h.

References fIsVoid.

74 { return ! fIsVoid; };

QwPartialTrack::fIsVoid

Bool_t fIsVoid

marked as being void

Definition: QwPartialTrack.h:185

bool QwPartialTrack::IsUsed ( ) const

inline

Definition at line 75 of file QwPartialTrack.h.

References fIsUsed.

75 { return fIsUsed; };

QwPartialTrack::fIsUsed

Bool_t fIsUsed

used (part of a track)

Definition: QwPartialTrack.h:184

bool QwPartialTrack::IsValid ( ) const

inline

Definition at line 72 of file QwPartialTrack.h.

References fIsVoid.

Referenced by PrintValid(), and QwTrackingTreeCombine::TlTreeCombine().

72 { return ! fIsVoid; };

QwPartialTrack::fIsVoid

Bool_t fIsVoid

marked as being void

Definition: QwPartialTrack.h:185

Here is the caller graph for this function:

bool QwPartialTrack::IsVoid ( ) const

inline

Definition at line 71 of file QwPartialTrack.h.

References fIsVoid.

Referenced by operator<<().

71 { return fIsVoid; };

QwPartialTrack::fIsVoid

Bool_t fIsVoid

marked as being void

Definition: QwPartialTrack.h:185

Here is the caller graph for this function:

QwPartialTrack & QwPartialTrack::operator= ( const QwPartialTrack & that )

Assignment operator.

Assignment operator

Definition at line 88 of file QwPartialTrack.cc.

References AddTreeLineList(), ClearTreeLines(), fAlone, fAverageResidual, fChi, fCov, fIsGood, fIsUsed, fIsVoid, fNumHits, fNumMiss, fOffsetX, fOffsetY, fQwTreeLines, fSignedResidual, fSlopeX, fSlopeY, kNumDirections, VQwTrackingElement::operator=(), pR2hit, pR3hit, TResidual, uvR2hit, and uvR3hit.

 {
   if (this == &that) return *this;
 
   VQwTrackingElement::operator=(that);
   
   fAlone = that.fAlone;
   
   fOffsetX = that.fOffsetX;
   fOffsetY = that.fOffsetY;
   fSlopeX = that.fSlopeX;
   fSlopeY = that.fSlopeY;
 
   fChi = that.fChi;
   fAverageResidual=that.fAverageResidual;
 
   for (size_t i = 0; i < 4; ++i)
     for (size_t j = 0; j < 4; ++j)
       fCov[i][j] = that.fCov[i][j];
 
   // this is a pointer
   for (size_t i = 0; i < kNumDirections; ++i){
     TResidual[i]=that.TResidual[i];
   }
     
   for (size_t i=0;i< 12;++i)
     fSignedResidual[i]=that.fSignedResidual[i];
   
   fIsUsed = that.fIsUsed;
   fIsVoid = that.fIsVoid;
   fIsGood = that.fIsGood;
 
   fNumMiss = that.fNumMiss;
   fNumHits = that.fNumHits;
 
   for (size_t i = 0; i < 3; ++i) {
     pR2hit[i] = that.pR2hit[i];
     uvR2hit[i] = that.uvR2hit[i];
     pR3hit[i] = that.pR3hit[i];
   
   uvR3hit[i] = that.uvR3hit[i];
   }
 
   // Copy tree lines
   ClearTreeLines();
   AddTreeLineList(that.fQwTreeLines);
 
   return *this;
 }

Here is the call graph for this function:

void QwPartialTrack::Print ( const Option_t * options = 0 ) const

Print some debugging information

Definition at line 285 of file QwPartialTrack.cc.

References QwLog::endl(), and QwVerbose.

 {
   if (!this) return;
   QwVerbose << *this << QwLog::endl;
 }

Here is the call graph for this function:

void QwPartialTrack::PrintTreeLines ( Option_t * option = "" ) const

Print the list of tree lines.

Definition at line 268 of file QwPartialTrack.cc.

References fQwTreeLines, and QwTreeLine::next.

 {
   for (std::vector<QwTreeLine*>::const_iterator treeline = fQwTreeLines.begin();
        treeline != fQwTreeLines.end(); treeline++) {
     std::cout << **treeline << std::endl;
     QwTreeLine* tl = (*treeline)->next;
     while (tl) {
       std::cout << *tl << std::endl;
       tl = tl->next;
     }
   }
 }

void QwPartialTrack::PrintValid ( )

Print some debugging information for valid partial tracks

Definition at line 294 of file QwPartialTrack.cc.

References QwLog::endl(), IsValid(), and QwVerbose.

 {
   if (!this) return;
   if (this->IsValid()) QwVerbose << *this << QwLog::endl;
 }

Here is the call graph for this function:

void QwPartialTrack::Reset ( Option_t * option = "" )

Definition at line 211 of file QwPartialTrack.cc.

References ResetTreeLines().

 {
   ResetTreeLines(option);
 }

Here is the call graph for this function:

void QwPartialTrack::ResetTreeLines ( Option_t * option = "" )

Reset the list of tree lines.

Definition at line 262 of file QwPartialTrack.cc.

References ClearTreeLines().

Referenced by Reset().

 {
   ClearTreeLines();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwPartialTrack::RotateCoordinates ( )

Rotate coordinates to right octant.

This method rotate track geometry information to the right postion according to the octant number. Remember the initial position is at octant 3.

Definition at line 496 of file QwPartialTrack.cc.

References fOffsetX, fOffsetY, fSlopeX, fSlopeY, and VQwTrackingElement::GetOctant().

Referenced by QwTrackingTreeCombine::TlTreeCombine().

                                       {
   //pick two points
   double z1=0,z2=-250;
   double x1=0,y1=0,x2=0,y2=0;
   double x3=0,y3=0,x4=0,y4=0;
   double PI=3.1415926;
   x1=fOffsetX+fSlopeX*z1;
   y1=fOffsetY+fSlopeY*z1;
   x2=fOffsetX+fSlopeX*z2;
   y2=fOffsetY+fSlopeY*z2;
   int oct=GetOctant();
   double rotate=oct==8? 135: (3-oct)*45;
   double Cos=cos(rotate*PI/180);
   double Sin=sin(rotate*PI/180);
   x3=Cos*x1+Sin*y1;
   y3=-Sin*x1+Cos*y1;
   x4=Cos*x2+Sin*y2;
   y4=-Sin*x2+Cos*y2;
 
   fOffsetX=x3;
   fOffsetY=y3;
 
   fSlopeX=(x4-x3)/(z2-z1);
   fSlopeY=(y4-y3)/(z2-z1);
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwPartialTrack::RotateRotator ( const QwDetectorInfo * geometry )

Definition at line 522 of file QwPartialTrack.cc.

References fOffsetX, fOffsetY, fSlopeX, fSlopeY, QwDetectorInfo::GetRotatorPitch(), QwDetectorInfo::GetRotatorRoll(), and QwDetectorInfo::GetRotatorYaw().

Referenced by QwTrackingTreeCombine::TlTreeCombine().

                                                                 {
 
   double z1=476.7,z2=500.0;
   double x1=0,y1=0,x2=0,y2=0;
   double pitch=0,yaw=0,roll=0;
   double rotatorZpos = 476.7;
 
   x1=fOffsetX+fSlopeX*z1;
   y1=fOffsetY+fSlopeY*z1;
   x2=fOffsetX+fSlopeX*z2;
   y2=fOffsetY+fSlopeY*z2;
 
   z1 = z1 - rotatorZpos;
   z2 = z2 - rotatorZpos;
 
   TVector3 v1(x1,y1,z1);
   TVector3 v2(x2,y2,z2);
 
   TVector3 xaxis(1,0,0);   //x-axis of rotator
   TVector3 yaxis(0,1,0);  //y-axis of rotator
   TVector3 zaxis(0,0,1);  //z-axis of rotator
 
    pitch = detector->GetRotatorPitch();
    yaw = detector->GetRotatorYaw();
    roll = detector->GetRotatorRoll();
 
   v1.Rotate(pitch,xaxis);
   v2.Rotate(pitch,xaxis);
   yaxis.Rotate(pitch,xaxis);  //rotating y-axis for yaw next
   zaxis.Rotate(pitch,xaxis);  //rotating z-axis for roll later
 
   v1.Rotate(yaw,yaxis);
   v2.Rotate(yaw,yaxis);
   zaxis.Rotate(yaw,yaxis);  //rotating z-axis for roll next
 
   v1.Rotate(roll,zaxis);
   v2.Rotate(roll,zaxis);
 
   v1.SetZ( v1.Z() + rotatorZpos );
   v2.SetZ( v2.Z() + rotatorZpos );
 
   fSlopeX = (v2.X()-v1.X())/(v2.Z()-v1.Z());
   fSlopeY = (v2.Y()-v1.Y())/(v2.Z()-v1.Z());
 
   fOffsetX = v1.X()-v1.Z()*fSlopeX;
   fOffsetY = v1.Y()-v1.Z()*fSlopeY;
 
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwPartialTrack::SetAlone ( const int alone )

inline

Definition at line 159 of file QwPartialTrack.h.

References fAlone.

Referenced by QwTrackingTreeCombine::TlTreeCombine().

159 { fAlone = alone; }

QwPartialTrack::fAlone

int fAlone

Definition: QwPartialTrack.h:197

Here is the caller graph for this function:

void QwPartialTrack::SetAverageResidual ( const double residual )

inline

Definition at line 154 of file QwPartialTrack.h.

References fAverageResidual.

Referenced by QwTrackingTreeCombine::TcTreeLineCombine().

154 { fAverageResidual = residual; }

QwPartialTrack::fAverageResidual

double fAverageResidual

number of Plane 0 Treelines

Definition: QwPartialTrack.h:202

Here is the caller graph for this function:

void QwPartialTrack::SetAverageResidual ( )

inline

Definition at line 156 of file QwPartialTrack.h.

References CalculateAverageResidual(), and fAverageResidual.

156 { fAverageResidual = CalculateAverageResidual(); }

QwPartialTrack::fAverageResidual

double fAverageResidual

number of Plane 0 Treelines

Definition: QwPartialTrack.h:202

QwPartialTrack::CalculateAverageResidual

double CalculateAverageResidual()

Definition: QwPartialTrack.cc:188

Here is the call graph for this function:

QwPartialTrack & QwPartialTrack::SmearAnglePhi ( const double sigma )

Smear the phi angle.

Smear the phi angle of the partial track

Parameters

sigma Standard deviation in phi

Definition at line 380 of file QwPartialTrack.cc.

References fSlopeX, fSlopeY, GetMomentumDirectionPhi(), and GetMomentumDirectionTheta().

 {
   double theta = GetMomentumDirectionTheta();
   double phi = GetMomentumDirectionPhi();
   phi = gRandom->Gaus(phi, sigma);
   fSlopeX = sin(theta) * cos(phi);
   fSlopeY = sin(theta) * sin(phi);
   return *this;
 }

Here is the call graph for this function:

QwPartialTrack & QwPartialTrack::SmearAngleTheta ( const double sigma )

Smear the theta angle.

Smear the theta angle of the partial track

Parameters

sigma Standard deviation in theta

Definition at line 366 of file QwPartialTrack.cc.

References fSlopeX, fSlopeY, GetMomentumDirectionPhi(), and GetMomentumDirectionTheta().

 {
   double theta = GetMomentumDirectionTheta();
   double phi = GetMomentumDirectionPhi();
   theta = gRandom->Gaus(theta, sigma);
   fSlopeX = sin(theta) * cos(phi);
   fSlopeY = sin(theta) * sin(phi);
   return *this;
 }

Here is the call graph for this function:

QwPartialTrack & QwPartialTrack::SmearPosition	(	const double	sigma_x,
		const double	sigma_y
	)

Smear the position.

Smear the position of the partial track

Parameters

sigma_x	Standard deviation in x
sigma_y	Standard deviation in y

Definition at line 355 of file QwPartialTrack.cc.

References fOffsetX, and fOffsetY.

 {
   fOffsetX = gRandom->Gaus(fOffsetX, sigma_x);
   fOffsetY = gRandom->Gaus(fOffsetY, sigma_y);
   return *this;
 }

Friends And Related Function Documentation

std::ostream& operator<<	(	std::ostream &	stream,
		const QwPartialTrack &	pt
	)

friend

Output stream operator overloading

Definition at line 303 of file QwPartialTrack.cc.

 {
   stream << "pt: ";
   if (pt.GetRegion() != kRegionIDNull)
     stream << "(" << pt.GetRegion() << "/" << "?UD"[pt.GetPackage()] << "); ";
   stream << "x,y(z=0) = (" << pt.fOffsetX/Qw::cm << " cm, " << pt.fOffsetY/Qw::cm << " cm), ";
   stream << "d(x,y)/dz = (" << pt.fSlopeX << ", " << pt.fSlopeY << ")";
 
   double x=pt.fOffsetX/Qw::cm-663*pt.fSlopeX;
   double y=pt.fOffsetY/Qw::cm-663*pt.fSlopeY;
   stream << " downstream target: " << x << ", " << y; 
   if (pt.fChi > 0.0) { // parttrack has been fitted
     stream << ", chi = " << pt.fChi;
   }
   if (pt.fAverageResidual > 0.0) {
     stream << ", res = " << pt.fAverageResidual;
   }
   if (pt.IsVoid()) stream << " (void)";
   return stream;
 }