#include <ZHelix.h>

Inheritance diagram for ZHelix:

Public Member Functions
	ZHelix ()

	ZHelix (Double_t Azim, Double_t QovR, Double_t QxDh, Double_t refx, Double_t refy, Double_t refz, Double_t TDip, Double_t phii, Double_t phio, Float_t Chi2, Int_t NDoF, EZHelixRangeType RangeType=kHelixPhi, Double_t RangeMin=0, Double_t RangeMax=2 *TMath::Pi())

virtual	~ZHelix ()

void	Phi2XYZ (Double_t phi, Double_t &x, Double_t &y, Double_t &z)

void	Phi2ZR (Double_t phi, Double_t &z, Double_t &r)

void	SetRange (EZHelixRangeType RangeType, Double_t RangeMin, Double_t RangeMax)

void	SetTrackType (TrackType tt)

virtual void	SetPoints (Option_t *option="3D")

Double_t	X2Phi (Double_t x)

Double_t	Y2Phi (Double_t y)

Double_t	Z2Phi (Double_t z)

Double_t	R2Phi (Double_t r)

Double_t	Phi2S (Double_t phi)

virtual void	Print (Option_t *option="") const

virtual void	Draw (Option_t *option="3D")

virtual void	Paint (Option_t *option="3D")

virtual Int_t	SetNextPoint (Double_t x, Double_t y, Double_t z)

virtual void	SetPoint (Int_t point, Double_t x, Double_t y, Double_t z)

virtual void	Delete (Option_t *option)

virtual void	DrawClass ()

virtual void	DrawClone (Option_t *option)

virtual void	SetDrawOption (Option_t *option)

virtual void	ExecuteEvent (Int_t event, Int_t px, Int_t py)

Protected Attributes
TrackType	fTrackType

Double_t	fAzim

Double_t	fQovR

Double_t	fQxDh

Double_t	fTDip

Double_t	fRefX

Double_t	fRefY

Double_t	fRefZ

Double_t	fRange [2]

EZHelixRangeType	fRType

Float_t	fChi2

Int_t	fNDoF

Double_t	fSinAzim

Double_t	fCosAzim

Int_t	fRSign

Double_t	fPhiI
	sign of radius

Double_t	fPhiO

Bool_t	fEnable

Detailed Description

Definition at line 42 of file ZHelix.h.

Constructor & Destructor Documentation

◆ ZHelix() [1/2]

ZHelix::ZHelix ( )

◆ ZHelix() [2/2]

ZHelix::ZHelix	(	Double_t	Azim,
		Double_t	QovR,
		Double_t	QxDh,
		Double_t	refx,
		Double_t	refy,
		Double_t	refz,
		Double_t	TDip,
		Double_t	phii,
		Double_t	phio,
		Float_t	Chi2,
		Int_t	NDoF,
		EZHelixRangeType	RangeType = `kHelixPhi`,
		Double_t	RangeMin = `0`,
		Double_t	RangeMax = `2*TMath::Pi()`
	)

Definition at line 52 of file ZHelix.cxx.

                                                     :
        TPolyLine3D() {
    //
    // ZHelix normal constructor
    if ( gDebug ) cout << "ZHelix normal ctor called" << endl;
 
    fAzim     = Azim;
    fSinAzim  = TMath::Sin(Azim);
    fCosAzim  = TMath::Cos(Azim);
    fQovR     = QovR;
    fQxDh     = QxDh;
    fRefX     = refx;
    fRefY     = refy;
    fRefZ     = refz;
    fTDip     = TDip;
    fPhiI     = phii;
    fPhiO     = phio;
    fChi2     = Chi2;
    fNDoF     = NDoF;
    fRType    = RangeType;
    fEnable   = kTRUE;
 
    fRSign = 0;
    //  fSegment = 1.;
    SetRange(RangeType, RangeMin, RangeMax);
}

◆ ~ZHelix()

ZHelix::~ZHelix ( )

virtual

Definition at line 87 of file ZHelix.cxx.

                {
    //
    // ZHelix standard destructor
    if ( gDebug ) {
        cout << "ZHelix dtor called" << endl;
    }
}

Member Function Documentation

◆ Delete()

virtual void ZHelix::Delete ( Option_t * option )

inlinevirtual

Definition at line 97 of file ZHelix.h.

97{ TObject::Delete(option); }

◆ Draw()

void ZHelix::Draw ( Option_t * option = "3D" )

virtual

Definition at line 367 of file ZHelix.cxx.

                                  {
    //
    // Draw this helix with its current attributes.
    AppendPad(option);
}

◆ DrawClass()

virtual void ZHelix::DrawClass ( )

inlinevirtual

Definition at line 98 of file ZHelix.h.

98{ TObject::DrawClass(); }

◆ DrawClone()

virtual void ZHelix::DrawClone ( Option_t * option )

inlinevirtual

Definition at line 99 of file ZHelix.h.

99{ TObject::DrawClone(option); }

◆ ExecuteEvent()

virtual void ZHelix::ExecuteEvent	(	Int_t	event,
		Int_t	px,
		Int_t	py
	)

inlinevirtual

Definition at line 101 of file ZHelix.h.

101{;}

◆ Paint()

void ZHelix::Paint ( Option_t * option = "3D" )

virtual

Definition at line 375 of file ZHelix.cxx.

                                   {
    //
    // Draw this helix with its current attributes
    TString opt = option;
    opt.ToUpper();
 
    if ( ( (fRange[1] < fRange[0]) || (fRange[1] > fRange[0]) ) && fEnable ) {
        SetPoints(option);
        TPolyLine3D::Paint("SAME");
    }
}

◆ Phi2S()

Double_t ZHelix::Phi2S ( Double_t phi )

Definition at line 389 of file ZHelix.cxx.

                                   {
    //
    // calculates pathlength for phi
 
    return -(phi-fAzim)/fQovR;
 
}

◆ Phi2XYZ()

void ZHelix::Phi2XYZ	(	Double_t	phi,
		Double_t &	x,
		Double_t &	y,
		Double_t &	z
	)

Definition at line 246 of file ZHelix.cxx.

                                                                        {
    //
    // Calculate XYZ for a given azimuth
    x = fRefX - TMath::Sin(phi) / fQovR + fSinAzim * (1/fQovR + fQxDh);
    y = fRefY + TMath::Cos(phi) / fQovR - fCosAzim * (1/fQovR + fQxDh);
    z = fRefZ + (fAzim - phi) * fTDip / fQovR;
}

Referenced by Phi2ZR(), R2Phi(), and SetPoints().

◆ Phi2ZR()

void ZHelix::Phi2ZR	(	Double_t	phi,
		Double_t &	z,
		Double_t &	r
	)

Definition at line 256 of file ZHelix.cxx.

                                                          {
    //
    // Calculate ZR for a given azimuth
    Double_t x, y;
    Phi2XYZ(phi, x, y, z);
    r = TMath::Sqrt(TMath::Power(x,2) + TMath::Power(y,2));
}

Referenced by SetPoints().

◆ Print()

void ZHelix::Print ( Option_t * option = "" ) const

virtual

Definition at line 342 of file ZHelix.cxx.

                                         {
    //
    // Dump the helix attributes
    TString opt = option;
    opt.ToUpper();
 
    cout << endl
    << "Parameters of helix " << this->GetName() << ":" << endl
    << "Azimuth   = " << fAzim << endl
    << "Q/R       = " << fQovR << endl
    << "Q*D_h     = " << fQxDh << endl
    << "Z_h       = " << fRefZ << endl
    << "cot(Dip)  = " << fTDip << endl
    << "PhiI      = " << fPhiI << endl
    << "PhiO      = " << fPhiO << endl
    << "Chi2      = " << fChi2 << endl
    << "D.o.F.    = " << fNDoF << endl
    << "RangeType = " << fRType << endl
    << "Range: " << fRange[0] << " -> " << fRange[1] << endl;
    if ( opt.Contains("ALL") ) TPolyLine3D::Print("ALL");
    cout << endl;
}

◆ R2Phi()

Double_t ZHelix::R2Phi ( Double_t r )

Definition at line 296 of file ZHelix.cxx.

                                 {
    //
    // Convert r to phi
    Double_t k = 1/fQovR + fQxDh;
    Double_t my_value = fQovR/(2*k) * (-TMath::Power(r,2) + 1/TMath::Power(fQovR,2) + TMath::Power(k,2));
    // check if my_value is ouside acos validity range
    if ( (my_value < -1) || (my_value > 1 ) )
        return 999999999;
    Double_t my_phi = TMath::ACos(my_value);
 
    Double_t phi1 = fAzim - my_phi;
    Double_t phi2 = fAzim + my_phi;
 
    Double_t x1  = 0;
    Double_t y1  = 0;
    Double_t z1  = 0;
    Double_t x2  = 0;
    Double_t y2  = 0;
    Double_t z2  = 0;
    Double_t phi = 0;;
    this->Phi2XYZ(phi1,x1,y1,z1);
    this->Phi2XYZ(phi2,x2,y2,z2);
 
 
    // check if phi1 would result in a z value with the correct sign corrsponding
    // to fTDip
    if ( TMath::Sign(1.,fTDip) == TMath::Sign(1.,(z1-fRefZ)) ) {
        // check if phi2 would result in a z value with the correct sign corrsponding
        // to fTDip
        if ( TMath::Sign(1.,fTDip) == TMath::Sign(1.,(z2-fRefZ)) ) {
            // take phi with smallest resulting z difference to fRefZ (Zh)
            if ( TMath::Abs(z1-fRefZ) <= TMath::Abs(z2-fRefZ) )
                phi = phi1;
            else
                phi = phi2;
        } else
            phi = phi1;
    } else {
        phi = phi2;
    }
 
    return phi;
}

Referenced by SetRange().

◆ SetDrawOption()

virtual void ZHelix::SetDrawOption ( Option_t * option )

inlinevirtual

Definition at line 100 of file ZHelix.h.

100{ TObject::SetDrawOption(option); }

◆ SetNextPoint()

virtual Int_t ZHelix::SetNextPoint	(	Double_t	x,
		Double_t	y,
		Double_t	z
	)

inlinevirtual

Definition at line 91 of file ZHelix.h.

                                                                        {
    return TPolyLine3D::SetNextPoint(x, y, z);
  }

◆ SetPoint()

virtual void ZHelix::SetPoint	(	Int_t	point,
		Double_t	x,
		Double_t	y,
		Double_t	z
	)

inlinevirtual

Definition at line 94 of file ZHelix.h.

                                                                                {
    TPolyLine3D::SetPoint(point, x, y, z);
  }

Referenced by SetPoints().

◆ SetPoints()

void ZHelix::SetPoints ( Option_t * option = "3D" )

virtual

Definition at line 127 of file ZHelix.cxx.

                                       {
    //
    // Set PolyLine3D points
    TString opt = option;
    opt.ToUpper();
 
    // Get view
    BesView *view = dynamic_cast<BesView*>(gPad->GetView());
 
    // check with which sign the segmnents have to added to fRange[0]
    Int_t sign = 0;
    if ( fQovR >= 0 )
        sign = -1;
    else
        sign = +1;
 
    // calculate distance between ranges in rad
    Double_t degrad  = TMath::Pi() / 180.0;
    Double_t segment = 1. * degrad;   // 1 degree segments
    Double_t df = 0;
    if ( sign == +1 ) {
        if ( fRange[1] > fRange[0] )
            df = fRange[1]-fRange[0];
        else
            df = 2*TMath::Pi() + fRange[1]-fRange[0];
    } else {
        if ( fRange[1] < fRange[0] )
            df = fRange[1]-fRange[0];
        else
            df = 2*TMath::Pi() - fRange[1]+fRange[0];
    }
    // number of segments with degree distance
    Int_t nSeg = Int_t(TMath::Abs(df) / segment) + 1;
 
    // solve problem with nseg = 1 , means df too small compared to
    // segment -> new segment
    if ( nSeg < 2 ) {
        segment = 0.01 * degrad;
        nSeg = Int_t(TMath::Abs(df) / segment) + 1;
    }
    //  if ( nSeg < 10 ) nSeg = 10;
 
    // newly calculation of segment
    segment = TMath::Abs(df) / nSeg;
 
    Double_t phi, x, y, z, r;
    TPolyLine3D::SetPolyLine(nSeg+1);
 
    if ( opt.Contains("3D") ) {
        // 3D view
        for ( Int_t i = 0; i <= nSeg; i++ ) {
            phi = fRange[0] + segment * i * sign;
            Phi2XYZ(phi, x, y, z);
            TPolyLine3D::SetPoint(i, x, y, z);
        }
    } else if ( opt.Contains("XY") ) {
        // XY view
        for ( Int_t i = 0; i <= nSeg; i++ ) {
            phi = fRange[0] + segment * i * sign;
            Phi2XYZ(phi, x, y, z);
            TPolyLine3D::SetPoint(i, x, y, 0);
        }
    } else if ( opt.Contains("ZR") ) {
        // ZR view
 
        // The track is drawn either in the upper or lower region of
        // the ZR view depending on how many of its segements lay in
        // each region. If these numbers are unbalanced, the starting
        // point of the track is taken as criterium
        Int_t isgn = 0;
        for ( Int_t i = 0; i <= nSeg; i++ ) {
 
            // Get ZR coordinates
            phi = fRange[0] + segment * i * sign;
            Phi2ZR(phi, z, r);
 
            // Set point
            TPolyLine3D::SetPoint(i, z, r, 0);
 
            // Get sign for R coordinate
            isgn += view->GetRSign(phi);
            if ( i == 0 ) fRSign = view->GetRSign(phi);
        }
        if ( isgn != 0 ) fRSign = TMath::Sign(1, isgn);
 
        // The vertex fitted tracks need to be displaced otherwise they are
        // are not drawn to emanate from the vertex.
 
        Float_t z, r, rref, rdisp;
 
        rref = TMath::Sqrt((fRefX*fRefX)+(fRefY*fRefY));
 
        if ( (fTrackType != kVctpar) || (fTrackType != kVcparsec) || (fTrackType != kZttrprm) || (fTrackType != kZttrsec) ) {
 
            // Set R sign for each point
            for ( Int_t i = 0; i <= nSeg; i++ ) {
                z = GetP()[i*3];
                r = GetP()[i*3+1];
                SetPoint(i, z, r*fRSign, 0);
            }
 
        } else {
 
            for ( Int_t i = 0; i <= nSeg; i++ ) {
                z = GetP()[i*3];
                r = GetP()[i*3+1];
                rdisp = r*fRSign;
                if (fRSign < 0) rdisp += 2*rref;
                SetPoint(i, z, rdisp, 0);
            }
 
        }
 
    }
 
}

Referenced by Paint().

◆ SetRange()

void ZHelix::SetRange	(	EZHelixRangeType	RangeType,
		Double_t	RangeMin,
		Double_t	RangeMax
	)

Definition at line 97 of file ZHelix.cxx.

                                                            {
    //
    // Set range for drawing
    switch (RangeType) {
    case kHelixPhi:
        fRange[0] = RangeMin;
        fRange[1] = RangeMax;
        break;
    case kHelixX:
        fRange[0] = X2Phi(RangeMin);
        fRange[1] = X2Phi(RangeMax);
        break;
    case kHelixY:
        fRange[0] = Y2Phi(RangeMin);
        fRange[1] = Y2Phi(RangeMax);
        break;
    case kHelixZ:
        fRange[0] = Z2Phi(RangeMin);
        fRange[1] = Z2Phi(RangeMax);
        break;
    case kHelixR:
        fRange[0] = R2Phi(RangeMin);
        fRange[1] = R2Phi(RangeMax);
        break;
    }
}

Referenced by ZHelix().

◆ SetTrackType()

void ZHelix::SetTrackType ( TrackType tt )

inline

Definition at line 79 of file ZHelix.h.

79{fTrackType = tt;};

◆ X2Phi()

Double_t ZHelix::X2Phi ( Double_t x )

Definition at line 266 of file ZHelix.cxx.

                                 {
    //
    // Convert x to phi
    Double_t phi = TMath::ASin(fQovR * (fSinAzim*(1/fQovR + fQxDh) - x + fRefX));
 
    return phi;
}

Referenced by SetRange().

◆ Y2Phi()

Double_t ZHelix::Y2Phi ( Double_t y )

Definition at line 276 of file ZHelix.cxx.

                                 {
    //
    // Convert y to phi
    Double_t phi = TMath::ACos(fQovR * (fCosAzim*(1/fQovR + fQxDh) + y - fRefY));
 
    return phi;
}

Referenced by SetRange().

◆ Z2Phi()

Double_t ZHelix::Z2Phi ( Double_t z )

Definition at line 286 of file ZHelix.cxx.

                                 {
    //
    // Convert z to phi
    Double_t phi = fAzim - fQovR * (z - fRefZ) / fTDip;
 
    return phi;
}

Referenced by SetRange().

Member Data Documentation

◆ fAzim

Double_t ZHelix::fAzim

protected

Definition at line 46 of file ZHelix.h.

Referenced by Phi2S(), Phi2XYZ(), Print(), R2Phi(), Z2Phi(), and ZHelix().

◆ fChi2

Float_t ZHelix::fChi2

protected

Definition at line 55 of file ZHelix.h.

Referenced by Print(), and ZHelix().

◆ fCosAzim

Double_t ZHelix::fCosAzim

protected

Definition at line 58 of file ZHelix.h.

Referenced by Phi2XYZ(), Y2Phi(), and ZHelix().

◆ fEnable

Bool_t ZHelix::fEnable

protected

Definition at line 63 of file ZHelix.h.

Referenced by Paint(), and ZHelix().

◆ fNDoF

Int_t ZHelix::fNDoF

protected

Definition at line 56 of file ZHelix.h.

Referenced by Print(), and ZHelix().

◆ fPhiI

Double_t ZHelix::fPhiI

protected

sign of radius

Definition at line 61 of file ZHelix.h.

Referenced by Print(), and ZHelix().

◆ fPhiO

Double_t ZHelix::fPhiO

protected

Definition at line 62 of file ZHelix.h.

Referenced by Print(), and ZHelix().

◆ fQovR

Double_t ZHelix::fQovR

protected

Definition at line 47 of file ZHelix.h.

Referenced by Phi2S(), Phi2XYZ(), Print(), R2Phi(), SetPoints(), X2Phi(), Y2Phi(), Z2Phi(), and ZHelix().

◆ fQxDh

Double_t ZHelix::fQxDh

protected

Definition at line 48 of file ZHelix.h.

Referenced by Phi2XYZ(), Print(), R2Phi(), X2Phi(), Y2Phi(), and ZHelix().

◆ fRange

Double_t ZHelix::fRange[2]

protected

Definition at line 53 of file ZHelix.h.

Referenced by Paint(), Print(), SetPoints(), and SetRange().

◆ fRefX

Double_t ZHelix::fRefX

protected

Definition at line 50 of file ZHelix.h.

Referenced by Phi2XYZ(), SetPoints(), X2Phi(), and ZHelix().

◆ fRefY

Double_t ZHelix::fRefY

protected

Definition at line 51 of file ZHelix.h.

Referenced by Phi2XYZ(), SetPoints(), Y2Phi(), and ZHelix().

◆ fRefZ

Double_t ZHelix::fRefZ

protected

Definition at line 52 of file ZHelix.h.

Referenced by Phi2XYZ(), Print(), R2Phi(), Z2Phi(), and ZHelix().

◆ fRSign

Int_t ZHelix::fRSign

protected

Definition at line 60 of file ZHelix.h.

Referenced by SetPoints(), and ZHelix().

◆ fRType

EZHelixRangeType ZHelix::fRType

protected

Definition at line 54 of file ZHelix.h.

Referenced by Print(), and ZHelix().

◆ fSinAzim

Double_t ZHelix::fSinAzim

protected

Definition at line 57 of file ZHelix.h.

Referenced by Phi2XYZ(), X2Phi(), and ZHelix().

◆ fTDip

Double_t ZHelix::fTDip

protected

Definition at line 49 of file ZHelix.h.

Referenced by Phi2XYZ(), Print(), R2Phi(), Z2Phi(), and ZHelix().

◆ fTrackType

TrackType ZHelix::fTrackType

protected

Definition at line 45 of file ZHelix.h.

Referenced by SetPoints(), and SetTrackType().

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

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ZHelix() [1/2]

◆ ZHelix() [2/2]

◆ ~ZHelix()

Member Function Documentation

◆ Delete()

◆ Draw()

◆ DrawClass()

◆ DrawClone()

◆ ExecuteEvent()

◆ Paint()

◆ Phi2S()

◆ Phi2XYZ()

◆ Phi2ZR()

◆ Print()

◆ R2Phi()

◆ SetDrawOption()

◆ SetNextPoint()

◆ SetPoint()

◆ SetPoints()

◆ SetRange()

◆ SetTrackType()

◆ X2Phi()

◆ Y2Phi()

◆ Z2Phi()

Member Data Documentation

◆ fAzim

◆ fChi2

◆ fCosAzim

◆ fEnable

◆ fNDoF

◆ fPhiI

◆ fPhiO

◆ fQovR

◆ fQxDh

◆ fRange

◆ fRefX

◆ fRefY

◆ fRefZ

◆ fRSign

◆ fRType

◆ fSinAzim

◆ fTDip

◆ fTrackType