HelixBase.h

/*
 * HalHelixBase.h
 *
 *  Created on: 30-04-2022
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALHELIXBASE_H_
#define HALHELIXBASE_H_
 
#include <TObject.h>
#include <TVector3.h>
namespace Hal {
  class HelixBase : public TObject {
  private:
    static Double_t fgHelixBz;
    Double_t fX0, fY0, fZ0;
    Double_t fPhi0, fPhiCos, fPhiSin;
    Double_t fLambda, fLambdaCos, fLambdaSin;
    Double_t fCurv;
    Double_t fXc, fYc;
    Double_t fH;
    static Double_t GetBz();
    Double_t FudgePathLenght(const TVector3& vec) const;
    Double_t Distance(const TVector3& p, Bool_t scanPeriods = kTRUE) const;
 
  protected:
    HelixBase(Double_t x,
              Double_t y,
              Double_t z,
              Double_t px,
              Double_t py,
              Double_t pz,
              Double_t charge,
              Double_t conversion = 0.);
    void BasePathLength(Double_t r, Double_t x, Double_t y, Double_t& s1, Double_t& s2) const;
    void BaseSetParams(const TVector3& x, const TVector3& mom, Double_t charge, Double_t conversion = 0.);
    void BaseSetParams(Double_t x,
                       Double_t y,
                       Double_t z,
                       Double_t px,
                       Double_t py,
                       Double_t pz,
                       Double_t charge,
                       Double_t conversion = 0);
    inline Double_t BaseGetXcenter() const { return fXc; };
    inline Double_t BaseGetYcenter() const { return fYc; };
    inline Double_t BaseGetStartX() const { return fX0; };
    inline Double_t BaseGetStartY() const { return fY0; };
    inline Double_t BaseGetStartZ() const { return fZ0; };
    Double_t BasePathLength(const TVector3& point, Bool_t scanPeriods = kTRUE) const;
    Double_t BasePathLength(Double_t x, Double_t y) const;
    TVector3 BasePosition(Double_t s) const;
    TVector3 BaseGetStartPoint() const;
    TVector3 BaseMomentum(Double_t s) const;
    Int_t BaseIntersection(Double_t R, TVector3& x1, TVector3& x2) const;
    void BasePrint() const;
    void BasePathLength(Double_t r, Double_t& s1, Double_t& s2) const;
    void BaseShift(Double_t x, Double_t y, Double_t z);
    void BaseFullEval(Double_t s, TVector3& mom, TVector3& pos) const;
 
  public:
    HelixBase();
    HelixBase(const HelixBase& t);
    HelixBase(const TVector3& x, const TVector3& mom, Double_t charge, Double_t conversion = 0.);
    inline Double_t GetCurv() const { return fCurv; };
    inline Double_t GetPhi0() const { return fPhi0; };
    inline Double_t GetDipAngle() const { return fLambda; };
    inline Double_t GetH() const { return fH; };
    HelixBase& operator=(const HelixBase& helix);
    void PathLengths(const HelixBase& h, Double_t& s1, Double_t& s2) const;
    static void SetMagField(Double_t Bz);
    Double_t GetPathAbsMin(Double_t s1, Double_t s2) const;
    Double_t GetPeriod() const;
    Double_t GetRotationDirection() const;
    static Double_t MaxPath() { return 1E+9; };
    virtual ~HelixBase();
    ClassDef(HelixBase, 1)
  };
}  // namespace Hal
#endif /* HALHELIXBASE_H_ */