PndFtsHoughTrackCand Class Reference

Class for saving a FTS track cand. for Hough transform based FTS PR. More...

#include <PndFtsHoughTrackCand.h>

Inheritance diagram for PndFtsHoughTrackCand:

Public Member Functions
	PndFtsHoughTrackCand (PndFtsHoughTrackerTask *trackerTask=0)
	Set pointer to tracker task (super important as it provides an I/O interface to PandaRoot) More...
	~PndFtsHoughTrackCand ()
void	Print ()
Bool_t	isComplete () const
PndTrackCand	getPndTrackCand ()
PndTrack	getPndTrack ()
FairTrackParP	getTrackParPForHit (const UInt_t i)
Int_t	getCharge () const
TVector3	getP (const Double_t zLabSys) const
TVector3	getPos (const Double_t zLabSys) const
Double_t	getXLabSys (const Double_t zLabSys) const
Double_t	getThetaZxRad (const Double_t zLabSys) const
Double_t	getThetaZyRad (const Double_t) const
Double_t	getZLineParabola ()
Double_t	getZParabolaLine ()
void	SetZxLineBeforeDipole (const PndFtsHoughTracklet zxLineParabola)
void	SetZxParabola (const PndFtsHoughTracklet zxParabola)
void	SetZxLineBehindDipole (const PndFtsHoughTracklet zxParabolaLine)
void	SetZyLine (const PndFtsHoughTracklet zyLine)
	ClassDef (PndFtsHoughTrackCand, 1)

Private Member Functions
void	throwError (const TString s) const
	For error reporting. More...
void	throwIfZOutOfRange (const Double_t &zLabSys) const
void	addUniqueTrackletHits (const PndFtsHoughTracklet inTracklet)
Double_t	getQdivPzx () const
Double_t	getPYLab () const
std::pair< Double_t, Double_t >	getPZPXLabLine (const Double_t &zLabSys, const PndFtsHoughTracklet *const lineTracklet) const
std::pair< Double_t, Double_t >	getPZPXLabParabola (const Double_t &zLabSys) const
Double_t	getXOrYLabForLine (const Double_t &zLabSys, const PndFtsHoughTracklet *const lineTracklet) const
Double_t	getXLabForParabola (const Double_t &zLabSys) const

Private Attributes
PndFtsHoughTrackerTask *	fTrackerTask
Int_t	fVerbose
PndFtsHoughTracklet	fZxLineBeforeDipole
PndFtsHoughTracklet	fZxParabola
PndFtsHoughTracklet	fZxLineBehindDipole
PndFtsHoughTracklet	fZyLine
PndTrackCand	fIntTrackCand
Double_t	fZCoordLineParabola
Double_t	fZCoordParabolaLine

Detailed Description

Class for saving a FTS track cand. for Hough transform based FTS PR.

Author

Martin J. Galuska <martin [dot] j [dot] galuska (at) physik [dot] uni (minus) giessen [dot] de>

The track cand. consists of several tracklets: In zx-plane line + parabola + line In zy-plane line

This class takes care of geometric calculations for tracklets.

Loosely modeled after tools/riemannfit/PndRiemannTrack.h

TODO Maybe it would be better not to derive from PndTrackCand, but to contain an object of PndTrackCand

Created: 24.01.2014

Definition at line 47 of file PndFtsHoughTrackCand.h.

Constructor & Destructor Documentation

PndFtsHoughTrackCand::PndFtsHoughTrackCand

(

PndFtsHoughTrackerTask *

trackerTask = 0

)

Set pointer to tracker task (super important as it provides an I/O interface to PandaRoot)

Definition at line 18 of file PndFtsHoughTrackCand.cxx.

References fTrackerTask, fVerbose, and PndFtsHoughTrackerTask::GetVerbose().

                                                                              :
                                fTrackerTask(trackerTask),

                                fVerbose(0),

                                fZxLineBeforeDipole(0., trackerTask), // TODO: It could be a problem here that I set the z reference value to 0.

                                fZxParabola(0., trackerTask),

                                fZxLineBehindDipole(0., trackerTask),

                                fZyLine(0., trackerTask),

                                fIntTrackCand(),

                                fZCoordLineParabola(0.),
                                fZCoordParabolaLine(0.)
{
        if (0==fTrackerTask){
                std::cout << "PndFtsHoughTrackCand FATAL ERROR Tracker task pointer not set.\n";
        } else {
                fVerbose = fTrackerTask->GetVerbose();
                if(3<fVerbose) std::cout << "PndFtsHoughTrackCand called with tracker ptr " << fTrackerTask << '\n';
        }
}

PndFtsHoughTrackCand::~PndFtsHoughTrackCand

(

)

Definition at line 44 of file PndFtsHoughTrackCand.cxx.

{

}

Member Function Documentation

void PndFtsHoughTrackCand::addUniqueTrackletHits ( const PndFtsHoughTracklet  inTracklet )

private

Definition at line 119 of file PndFtsHoughTrackCand.cxx.

References PndTrackCand::AddHit(), Double_t, fIntTrackCand, PndTrackCandHit::GetDetId(), PndTrackCandHit::GetHitId(), PndTrackCand::GetNHits(), PndTrackCandHit::GetRho(), PndTrackCand::GetSortedHit(), and PndTrackCand::HitInTrack().

Referenced by SetZxLineBeforeDipole(), SetZxLineBehindDipole(), SetZxParabola(), and SetZyLine().

{
        // add all the hits from the tracklet which are not already in the hitId vector
        for (UInt_t iHit = 0; iHit < inTracklet.GetNHits(); ++iHit)
        {
                PndTrackCandHit inHit = inTracklet.GetSortedHit(iHit);
                const Int_t inHitId = inHit.GetHitId();
                const Int_t inDetId = inHit.GetDetId();
                // if hit is NOT in track -1 is returned by HitInTrack, otherwise the index (>=0) in the HitId vector is returned
                if (-1==fIntTrackCand.HitInTrack(inDetId,inHitId)){
                        // add hit to track cand
                        const Double_t inRho = inHit.GetRho();
                        fIntTrackCand.AddHit(inDetId, inHitId, inRho);
                }
        }
}

PndFtsHoughTrackCand::ClassDef	(	PndFtsHoughTrackCand	,
		1
	)

Int_t PndFtsHoughTrackCand::getCharge ( ) const

inline

Definition at line 163 of file PndFtsHoughTrackCand.h.

References getQdivPzx().

Referenced by getPZPXLabLine(), getPZPXLabParabola(), and getTrackParPForHit().

                                            {
        if (0 < getQdivPzx()){
                return 1;
        }
        else{
                return -1;
        }
}

TVector3 PndFtsHoughTrackCand::getP

(

const Double_t

zLabSys

)

const

Definition at line 265 of file PndFtsHoughTrackCand.cxx.

References Double_t, fVerbose, fZCoordLineParabola, fZCoordParabolaLine, fZxLineBeforeDipole, fZxLineBehindDipole, getPYLab(), getPZPXLabLine(), getPZPXLabParabola(), isComplete(), mom, and throwIfZOutOfRange().

Referenced by getTrackParPForHit().

                                                               {
        throwIfZOutOfRange(zLabSys);

        TVector3 mom;

        if (kFALSE == isComplete())
        {
                std::cout << "getPforHit: Track cand. is not complete yet. Momentum will be calculated for incomplete track cand.\n";
        }

        Double_t pZLabSys;
        Double_t pXLabSys;
        std::pair<Double_t, Double_t> pZPXLabSys;

        // track model is assumed to be line+parabola+line in zx and line in zy
        if ( zLabSys <= fZCoordLineParabola ){
                // use 1st line in zx plane
                pZPXLabSys = getPZPXLabLine(zLabSys, &fZxLineBeforeDipole);
        } else if ( zLabSys < fZCoordParabolaLine ){
                // use "tangent" to parabola in zx plane
                pZPXLabSys = getPZPXLabParabola(zLabSys);
        } else {
                // use 2nd line in zx plane
                pZPXLabSys = getPZPXLabLine(zLabSys, &fZxLineBehindDipole);
        }

        pZLabSys = pZPXLabSys.first;
        pXLabSys = pZPXLabSys.second;
        const Double_t pYLabSys = getPYLab();
        mom.SetXYZ(pXLabSys, pYLabSys, pZLabSys);
        if (fVerbose > 10) std::cout << "P-Vector for z=" << zLabSys << " : " << mom.X() << " " << mom.Y() << " " << mom.Z() << '\n';
        return mom;
}

PndTrack PndFtsHoughTrackCand::getPndTrack

(

)

Definition at line 208 of file PndFtsHoughTrackCand.cxx.

References fIntTrackCand, PndTrackCand::GetNHits(), and getTrackParPForHit().

                                           {
        // calculates first and last parameter, but uses an empty PndTrackCand which has to be set lateron using SetTrackCandRef
        FairTrackParP firstPar, lastPar;
        PndTrackCand myCand;
        if (fIntTrackCand.GetNHits() > 0){
                firstPar = getTrackParPForHit(0);
                lastPar = getTrackParPForHit(fIntTrackCand.GetNHits()-1);
        }
        return PndTrack(firstPar, lastPar, myCand);
}

PndTrackCand PndFtsHoughTrackCand::getPndTrackCand

(

)

Definition at line 184 of file PndFtsHoughTrackCand.cxx.

References fIntTrackCand.

                                                   {
        return fIntTrackCand;

        //      // The following is what I used when I derived this class from PndTrackCand. Instead this class now contains a member of PndTrackCand
        //      // Maybe this is not necessary because PndFtsHoughTrackCand is derived from PndTrackCand
        //      PndTrackCand myCand;
        //      // copy all the hits from PndFtsHoughTrackCand *this to PndTrackCand myCand
        //      for (UInt_t iHit = 0; iHit < intTrackCand.GetNHits(); ++iHit)
        //      {
        //              PndTrackCandHit inHit = intTrackCand.GetSortedHit(iHit);
        //              const Int_t inHitId = inHit.GetHitId();
        //              const Int_t inDetId = inHit.GetDetId();
        //              const Double_t inRho = inHit.GetRho();
        //              myCand.AddHit(inDetId, inHitId, inRho);
        //      }
        //      return myCand;
}

TVector3 PndFtsHoughTrackCand::getPos

(

const Double_t

zLabSys

)

const

Definition at line 331 of file PndFtsHoughTrackCand.cxx.

References Double_t, fVerbose, fZyLine, getXLabSys(), getXOrYLabForLine(), isComplete(), and throwIfZOutOfRange().

Referenced by getTrackParPForHit(), and plotHoughTrackCand().

                                                                 {
        // calculates the point on the track based on the results from the Hough transforms using my track model for the given z

        throwIfZOutOfRange(zLabSys);

        if (kFALSE == isComplete()) std::cout << "getPositionForHit: Track cand. is not complete yet. Position will be calculated for incomplete track cand.\n";

        TVector3 position;
        // track model is assumed to be line+parabola+line in zx and line in zy
        const Double_t yLabSys = getXOrYLabForLine(zLabSys, &fZyLine);

        Double_t xLabSys = getXLabSys(zLabSys);

        position.SetXYZ(xLabSys, yLabSys, zLabSys);
        if (fVerbose > 10) std::cout << "Pos-Vector for z=" << zLabSys << " : " << position.X() << " " << position.Y() << " " << position.Z() << '\n';
        return position;
}

Double_t PndFtsHoughTrackCand::getPYLab ( ) const

inlineprivate

Definition at line 217 of file PndFtsHoughTrackCand.h.

References Double_t, fZxLineBeforeDipole, fZxParabola, fZyLine, getPZPXLabLine(), PndFtsHoughTracklet::getThetaRadVal(), and PndFtsHoughTracklet::getZRefLabSys().

Referenced by getP().

                                              {
                // theta in radian in zy plane given at z where first zx line meets the parabola
                const Double_t zRefLabSys = fZxParabola.getZRefLabSys();
                const Double_t thetaZyRad = fZyLine.getThetaRadVal();
                std::pair<Double_t, Double_t> pZPXLabSys = getPZPXLabLine(zRefLabSys, &fZxLineBeforeDipole);
                const Double_t pZLabSys = pZPXLabSys.first;
                Double_t pYLabSys = tan(thetaZyRad)*pZLabSys;
                return pYLabSys;
        }

std::pair< Double_t, Double_t > PndFtsHoughTrackCand::getPZPXLabLine ( const Double_t &  zLabSys, const PndFtsHoughTracklet *const  lineTracklet  ) const

inlineprivate

Definition at line 203 of file PndFtsHoughTrackCand.h.

References cos(), Double_t, fZxParabola, getCharge(), PndFtsHoughTracklet::getSecondVal(), PndFtsHoughTracklet::getThetaRadVal(), and sin().

Referenced by getP(), and getPYLab().

                                                                                                                                       { // zLabSys //[R.K.03/2017] unused variable(s)
                // theta in radian in zx plane given at z = zRefLabSys
                const Double_t thetaRad = lineTracklet->getThetaRadVal();
                const Double_t qDivPzx = fZxParabola.getSecondVal(); // Q/pzx
                const Double_t pZx = getCharge() / qDivPzx;

                const Double_t pZLabSys = pZx*cos(thetaRad);
                const Double_t pXLabSys = pZx*sin(thetaRad);
                std::pair<Double_t, Double_t> pZPXLabSys(pZLabSys, pXLabSys);

                return pZPXLabSys;
        }

std::pair< Double_t, Double_t > PndFtsHoughTrackCand::getPZPXLabParabola ( const Double_t &  zLabSys ) const

inlineprivate

Definition at line 188 of file PndFtsHoughTrackCand.h.

References cos(), Double_t, getCharge(), getQdivPzx(), getThetaZxRad(), and sin().

Referenced by getP().

                                                                                                  {

                const Double_t currentThetaRad = getThetaZxRad(zLabSys);

                const Double_t qDivPzx = getQdivPzx(); // Q/pzx
                const Double_t pZx = getCharge() / qDivPzx;

                const Double_t pZLabSys = pZx*cos(currentThetaRad);
                const Double_t pXLabSys = pZx*sin(currentThetaRad);
                std::pair<Double_t, Double_t> pZPXLabSys(pZLabSys, pXLabSys);

                return pZPXLabSys;
}

Double_t PndFtsHoughTrackCand::getQdivPzx ( ) const

inlineprivate

Definition at line 103 of file PndFtsHoughTrackCand.h.

References fZxParabola, and PndFtsHoughTracklet::getSecondVal().

Referenced by getCharge(), getPZPXLabParabola(), and getXLabForParabola().

{ return fZxParabola.getSecondVal(); };

Double_t PndFtsHoughTrackCand::getThetaZxRad

(

const Double_t

zLabSys

)

const

Definition at line 300 of file PndFtsHoughTrackCand.cxx.

References atan2(), Double_t, getXLabSys(), and throwIfZOutOfRange().

Referenced by getPZPXLabParabola().

                                                                        {
        // estimate angle by calculating 2 points (zLabSys1, xLabSys1) and (zLabSys2, xLabSys2) which are close together

        throwIfZOutOfRange(zLabSys);

        const Double_t epsilon = 0.5; // small positive number

        const Double_t xLabSys1 = getXLabSys(zLabSys-epsilon);
        const Double_t xLabSys2 = getXLabSys(zLabSys+epsilon);

        return atan2(xLabSys2-xLabSys1, 2*epsilon);
}

Double_t PndFtsHoughTrackCand::getThetaZyRad ( const Double_t  ) const

inline

Definition at line 75 of file PndFtsHoughTrackCand.h.

References fZyLine, and PndFtsHoughTracklet::getThetaRadVal().

Referenced by PndFtsHoughTrackFinder::LineBehindDipoleMatchesToLinePlusParabola().

{ return fZyLine.getThetaRadVal(); }; // zLabSys //[R.K.03/2017] unused variable(s)

FairTrackParP PndFtsHoughTrackCand::getTrackParPForHit

(

const UInt_t

i

)

Definition at line 220 of file PndFtsHoughTrackCand.cxx.

References Double_t, fIntTrackCand, fTrackerTask, getCharge(), PndFtsHoughTrackerTask::GetFtsHit(), PndFtsHoughTrackerTask::GetFtsHitPosErrors(), PndTrackCandHit::GetHitId(), getP(), getPos(), PndTrackCand::GetSortedHit(), isComplete(), mom, and throwError().

Referenced by getPndTrack().

                                                                         {
        // index is the index of the hit in intTrackCand, not in the FTS hit array
        // this method will sort the hitId vector of intTrackCand and therefore cannot be used for const track candidates

        // TODO: Check if all arguments are correct

        // Warn if we do not have a complete track candidate
        if (!isComplete()) std::cout << "getHit: You try to access hits before we have a complete track candidate.\n";

        // not necessary, I check later if I get a hit or not
        //if (index >= intTrackCand.GetNHits()) throwError("index in PndFtsHoughTrackCand::getTrackParPForHit is too large.");

        // Translate index in track candidate to hitId in FTS hit array
        UInt_t hitId = fIntTrackCand.GetSortedHit(index).GetHitId();

        // get position of hit
        const PndFtsHit *myHit = fTrackerTask->GetFtsHit(hitId);
        if (0 == myHit) throwError("Cannot get hit, probably the tracking has not finished or the index is too large.");

        // Take z of hit and calculate the position and momentum for that z using the results from the Hough transforms and my track model
        Double_t zLabSys = myHit->GetZ();

        // position should NOT come from hit, it should come from the pattern recognition track model
        // position error can be large (like 1* or 2* tube size) as the Kalman filter will adjust it.
        TVector3 hitPos = getPos(zLabSys);


        TVector3 hitPosError = fTrackerTask->GetFtsHitPosErrors(myHit);

        // momentum comes from the pattern recognition track model
        TVector3 mom = getP(zLabSys);
        TVector3 momError = 0.1*mom; // TODO: add correct values here

        // set plane as detector plane in which the hit is, FTS planes are parallel to xy plane
        TVector3 origin; // set origin of plane = position of hit
        myHit->Position(origin);
        TVector3 dj(1,0,0); // unit vector along x // TODO: Check if that is set correctly for FTS
        TVector3 dk(0,1,0); // unit vector along y // TODO: Check if that is set correctly for FTS

        // -----   Constructor with track parameters in LAB -----------------------------------
        // FairTrackParP::FairTrackParP(TVector3 pos, TVector3 Mom, TVector3 posErr, TVector3 MomErr, Int_t Q, TVector3 o, TVector3 dj, TVector3 dk)
        FairTrackParP result(hitPos, mom, hitPosError, momError, getCharge(), origin, dj, dk);
        return result;
}

Double_t PndFtsHoughTrackCand::getXLabForParabola ( const Double_t &  zLabSys ) const

private

Definition at line 138 of file PndFtsHoughTrackCand.cxx.

References a, Double_t, fabs(), fZxLineBeforeDipole, fZxParabola, getQdivPzx(), PndFtsHoughTracklet::getSecondVal(), PndFtsHoughTracklet::getThetaRadVal(), PndFtsHoughTracklet::getZRefLabSys(), sin(), and sqrt().

Referenced by getXLabSys().

                                                                               {
        // does not take option of varying B field into account
        // calculate x in lab sys for a given z position in lab sys for which the parabola assumption holds
        // theta in radian in zx plane given at z = zRefLabSys
        const Double_t thetaRad = fZxParabola.getThetaRadVal();
        const Double_t qDivPzx = getQdivPzx(); // Q/pzx
        const Double_t zRefLabSys = fZxParabola.getZRefLabSys();
        const Double_t zshifted = zLabSys-zRefLabSys;

        const Double_t interceptLine = fZxLineBeforeDipole.getSecondVal();

        const Double_t By = 1.;

        if (0==thetaRad){
                Double_t xParabolaSys = qDivPzx /2. * By * zshifted;
                // go from local parabola system to the lab system
                Double_t xLabSys = xParabolaSys + interceptLine;
                return xLabSys;
        }


        const Double_t tantheta = tan(thetaRad);
        const Double_t pzstuff = 1. / qDivPzx / By / sin(thetaRad);
        const Double_t ztantheta = zshifted / tantheta;
        const Double_t a = -ztantheta + pzstuff / tantheta;

        const Double_t wurzel = sqrt(a * a - 2. * pzstuff * zshifted - ztantheta * ztantheta);
        const Double_t x1ParabolaSys = a - wurzel;
        const Double_t x2ParabolaSys = a + wurzel;

        // go from local parabola system to the lab system
        Double_t x1LabSys = x1ParabolaSys + interceptLine;
        Double_t x2LabSys = x2ParabolaSys + interceptLine;

        // there are two analytical solutions, take the one closer to 0
        if (fabs(x1LabSys)<fabs(x2LabSys)) {
                return x1LabSys;
        } else {
                return x2LabSys;
        }
}

Double_t PndFtsHoughTrackCand::getXLabSys

(

const Double_t

zLabSys

)

const

Definition at line 313 of file PndFtsHoughTrackCand.cxx.

References fZCoordLineParabola, fZCoordParabolaLine, fZxLineBeforeDipole, fZxLineBehindDipole, getXLabForParabola(), getXOrYLabForLine(), throwError(), and throwIfZOutOfRange().

Referenced by PndFtsHoughSpace::CalculateHitPosFromIntersectionsWithZxTrackModel(), getPos(), and getThetaZxRad().

                                                                      {
        // calculates the corresponding x coordinate to a given z coordinate according to the track model

        throwIfZOutOfRange(zLabSys);

        if (zLabSys <= fZCoordLineParabola) {
                // use 1st line in zx plane
                return getXOrYLabForLine(zLabSys, &fZxLineBeforeDipole);
        } else if (zLabSys < fZCoordParabolaLine) {
                // use parabola in zx plane
                return getXLabForParabola(zLabSys);
        } else {
                // use 2nd line in zx plane
                return getXOrYLabForLine(zLabSys, &fZxLineBehindDipole);
        }
        throwError("Argument was not handled!");
}

Double_t PndFtsHoughTrackCand::getXOrYLabForLine ( const Double_t &  zLabSys, const PndFtsHoughTracklet *const  lineTracklet  ) const

inlineprivate

Definition at line 174 of file PndFtsHoughTrackCand.h.

References Double_t, PndFtsHoughTracklet::getSecondVal(), PndFtsHoughTracklet::getThetaRadVal(), and PndFtsHoughTracklet::getZRefLabSys().

Referenced by getPos(), and getXLabSys().

                                                                                                                              {
        // calculate x or y in lab sys for a given z position in lab sys for which the line assumption holds
        // theta in radian in zx plane given at z = zRefLabSys
        const Double_t thetaRad = lineTracklet->getThetaRadVal();
        const Double_t intercept = lineTracklet->getSecondVal();
        const Double_t zRefLabSys = lineTracklet->getZRefLabSys();
        const Double_t zshifted = zLabSys-zRefLabSys;

        Double_t xOrYLabSys = tan(thetaRad)*zshifted+intercept;

        return xOrYLabSys;
}

Double_t PndFtsHoughTrackCand::getZLineParabola ( )

inline

Definition at line 77 of file PndFtsHoughTrackCand.h.

References fZCoordLineParabola.

Referenced by plotHoughTrackCand().

{ return fZCoordLineParabola; }; // gets z coordinate in laboratory system where I switch from line before dipole to parabola within dipole (in zx plane)

Double_t PndFtsHoughTrackCand::getZParabolaLine ( )

inline

Definition at line 78 of file PndFtsHoughTrackCand.h.

References fZCoordParabolaLine.

Referenced by plotHoughTrackCand().

{ return fZCoordParabolaLine; }; // gets z coordinate in laboratory system where I switch from parabola within dipole to line behind dipole (in zx plane)

Bool_t PndFtsHoughTrackCand::isComplete ( ) const

inline

Definition at line 60 of file PndFtsHoughTrackCand.h.

References fZxLineBeforeDipole, fZxLineBehindDipole, fZxParabola, fZyLine, and PndFtsHoughTracklet::isSet().

Referenced by getP(), getPos(), and getTrackParPForHit().

{ return (fZxLineBeforeDipole.isSet() && fZxParabola.isSet() && fZxLineBehindDipole.isSet() && fZyLine.isSet()); };

void PndFtsHoughTrackCand::Print

(

)

Definition at line 94 of file PndFtsHoughTrackCand.cxx.

References fIntTrackCand, fZxLineBeforeDipole, fZxLineBehindDipole, fZxParabola, fZyLine, PndFtsHoughTracklet::isSet(), PndFtsHoughTracklet::Print(), and PndTrackCand::Print().

Referenced by PndFtsHoughTrackerTask::Exec().

                                 {
        std::cout << "=========== PndFtsHoughTrackCand::Print() ==========\n";
        fIntTrackCand.Print();
        if (kTRUE==fZxLineBeforeDipole.isSet()){
                std::cout << "zx plane\n\n";
                std::cout << "1st line: ";
                fZxLineBeforeDipole.Print();
        }
        if (kTRUE==fZxParabola.isSet()){
                std::cout << "parabola: ";
                fZxParabola.Print();
        }
        if (kTRUE==fZxLineBehindDipole.isSet()){
                std::cout << "2nd line: ";
                fZxLineBehindDipole.Print();
        }
        if (kTRUE==fZyLine.isSet()){
                std::cout << "zy plane\n\n";
                std::cout << "Line: ";
                fZyLine.Print();
        }
}

void PndFtsHoughTrackCand::SetZxLineBeforeDipole

(

const PndFtsHoughTracklet

zxLineParabola

)

Definition at line 50 of file PndFtsHoughTrackCand.cxx.

References addUniqueTrackletHits(), fZCoordLineParabola, fZxLineBeforeDipole, PndFtsHoughTracklet::getZRefLabSys(), and PndFtsHoughTracklet::isSet().

Referenced by PndFtsHoughTrackFinder::FindMatchingParabolaToLineBeforeDipoleZxAndAddLineBehindDipole().

                                                                                        {
        if (fZxLineBeforeDipole.isSet()) {
                std::cerr << "WARNING from PndFtsHoughTrackCand: Line before dipole in zx plane is set more than once! Ignore new values! Potentially FATAL ERROR!\n";
                return;
        }
        fZxLineBeforeDipole = zxLineParabola;
        fZCoordLineParabola = fZxLineBeforeDipole.getZRefLabSys();
        addUniqueTrackletHits(fZxLineBeforeDipole);
}

void PndFtsHoughTrackCand::SetZxLineBehindDipole

(

const PndFtsHoughTracklet

zxParabolaLine

)

Definition at line 71 of file PndFtsHoughTrackCand.cxx.

References addUniqueTrackletHits(), fZCoordParabolaLine, fZxLineBehindDipole, PndFtsHoughTracklet::getZRefLabSys(), and PndFtsHoughTracklet::isSet().

Referenced by PndFtsHoughTrackFinder::FindMatchingParabolaToLineBeforeDipoleZxAndAddLineBehindDipole().

                                                                                        {
        if (fZxLineBehindDipole.isSet()) {
                std::cerr << "WARNING from PndFtsHoughTrackCand: Line behind dipole in zx plane is set more than once! Ignore new values! Potentially FATAL ERROR!\n";
                return;
        }
        fZxLineBehindDipole = zxParabolaLine;
        fZCoordParabolaLine = fZxLineBehindDipole.getZRefLabSys();
        addUniqueTrackletHits(fZxLineBehindDipole);
}

void PndFtsHoughTrackCand::SetZxParabola

(

const PndFtsHoughTracklet

zxParabola

)

Definition at line 59 of file PndFtsHoughTrackCand.cxx.

References addUniqueTrackletHits(), fZCoordLineParabola, fZxParabola, PndFtsHoughTracklet::getZRefLabSys(), and PndFtsHoughTracklet::isSet().

Referenced by PndFtsHoughTrackFinder::FindMatchingParabolaToLineBeforeDipoleZxAndAddLineBehindDipole().

                                                                            {
        if (fZxParabola.isSet()) {
                std::cerr << "WARNING from PndFtsHoughTrackCand: Parabola inside dipole in zx plane is set more than once! Ignore new values! Potentially FATAL ERROR!\n";
                return;
        }
        // warn if line before dipole field and parabola within are not calculated wrt the same z reference value
        if (zxParabola.getZRefLabSys() != fZCoordLineParabola){
                std::cout << "WARNING from PndFtsHoughTrackCand: First line and parabola were not calculated wrt the same z position! Potentially FATAL ERROR!\n";
        }
        fZxParabola = zxParabola;
        addUniqueTrackletHits(fZxParabola);
}

void PndFtsHoughTrackCand::SetZyLine

(

const PndFtsHoughTracklet

zyLine

)

Definition at line 80 of file PndFtsHoughTrackCand.cxx.

References addUniqueTrackletHits(), fZyLine, and PndFtsHoughTracklet::isSet().

Referenced by PndFtsHoughTrackFinder::FindZyLineMatchingToLineParabolaLineInZx().

                                                                    {
        if (fZyLine.isSet()) {
                std::cerr << "WARNING from PndFtsHoughTrackCand: Line in zy plane is set more than once! Ignore new values! Potentially FATAL ERROR!\n";
                return;
        }
        fZyLine=zyLine;
        addUniqueTrackletHits(fZyLine);
}

void PndFtsHoughTrackCand::throwError ( const TString  s ) const

inlineprivate

For error reporting.

Definition at line 95 of file PndFtsHoughTrackCand.h.

Referenced by getTrackParPForHit(), getXLabSys(), and throwIfZOutOfRange().

{ throw std::runtime_error(s.Data()); };

void PndFtsHoughTrackCand::throwIfZOutOfRange ( const Double_t &  zLabSys ) const

inlineprivate

Definition at line 96 of file PndFtsHoughTrackCand.h.

References throwError().

Referenced by getP(), getPos(), getThetaZxRad(), and getXLabSys().

                                                               {
                if ( zLabSys <= 100 ) throwError("zLabSys is too small, track model is not valid in that region.");
                if ( zLabSys >= 1000 ) throwError("zLabSys is too big, track model is not valid in that region.");
        };

Member Data Documentation

PndTrackCand PndFtsHoughTrackCand::fIntTrackCand

private

Definition at line 144 of file PndFtsHoughTrackCand.h.

Referenced by addUniqueTrackletHits(), getPndTrack(), getPndTrackCand(), getTrackParPForHit(), and Print().

PndFtsHoughTrackerTask* PndFtsHoughTrackCand::fTrackerTask

private

Definition at line 92 of file PndFtsHoughTrackCand.h.

Referenced by getTrackParPForHit(), and PndFtsHoughTrackCand().

Int_t PndFtsHoughTrackCand::fVerbose

private

Definition at line 115 of file PndFtsHoughTrackCand.h.

Referenced by getP(), getPos(), and PndFtsHoughTrackCand().

Double_t PndFtsHoughTrackCand::fZCoordLineParabola

private

Definition at line 148 of file PndFtsHoughTrackCand.h.

Referenced by getP(), getXLabSys(), getZLineParabola(), SetZxLineBeforeDipole(), and SetZxParabola().

Double_t PndFtsHoughTrackCand::fZCoordParabolaLine

private

Definition at line 150 of file PndFtsHoughTrackCand.h.

Referenced by getP(), getXLabSys(), getZParabolaLine(), and SetZxLineBehindDipole().

PndFtsHoughTracklet PndFtsHoughTrackCand::fZxLineBeforeDipole

private

Definition at line 119 of file PndFtsHoughTrackCand.h.

Referenced by getP(), getPYLab(), getXLabForParabola(), getXLabSys(), isComplete(), Print(), and SetZxLineBeforeDipole().

PndFtsHoughTracklet PndFtsHoughTrackCand::fZxLineBehindDipole

private

Definition at line 131 of file PndFtsHoughTrackCand.h.

Referenced by getP(), getXLabSys(), isComplete(), Print(), and SetZxLineBehindDipole().

PndFtsHoughTracklet PndFtsHoughTrackCand::fZxParabola

private

Definition at line 125 of file PndFtsHoughTrackCand.h.

Referenced by getPYLab(), getPZPXLabLine(), getQdivPzx(), getXLabForParabola(), isComplete(), Print(), and SetZxParabola().

PndFtsHoughTracklet PndFtsHoughTrackCand::fZyLine

private

Definition at line 138 of file PndFtsHoughTrackCand.h.

Referenced by getPos(), getPYLab(), getThetaZyRad(), isComplete(), Print(), and SetZyLine().

---

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

• PndFtsHoughTrackCand.h
• PndFtsHoughTrackCand.cxx