include/root/TGeoTorus.h

0001 // @(#)root/base:$Id$
0002 // Author: Andrei Gheata   28/07/03
0003
0004 /*************************************************************************
0005  * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers.               *
0006  * All rights reserved.                                                  *
0007  *                                                                       *
0008  * For the licensing terms see $ROOTSYS/LICENSE.                         *
0009  * For the list of contributors see $ROOTSYS/README/CREDITS.             *
0010  *************************************************************************/
0011
0012 #ifndef ROOT_TGeoTorus
0013 #define ROOT_TGeoTorus
0014
0015 #include "TGeoBBox.h"
0016
0017 class TGeoTorus : public TGeoBBox {
0018 protected:
0019    // data members
0020    Double_t fR;    // axial radius
0021    Double_t fRmin; // inner radius
0022    Double_t fRmax; // outer radius
0023    Double_t fPhi1; // starting phi
0024    Double_t fDphi; // phi extent
0025                    // methods
0026    TGeoTorus(const TGeoTorus &) = delete;
0027    TGeoTorus &operator=(const TGeoTorus &) = delete;
0028
0029 public:
0030    Double_t Capacity() const override;
0031    Double_t Daxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
0032    Double_t DDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
0033    Double_t DDDaxis(const Double_t *pt, const Double_t *dir, Double_t t) const;
0034    Double_t ToBoundary(const Double_t *pt, const Double_t *dir, Double_t r, Bool_t in) const;
0035    Int_t SolveCubic(Double_t a, Double_t b, Double_t c, Double_t *x) const;
0036    Int_t SolveQuartic(Double_t a, Double_t b, Double_t c, Double_t d, Double_t *x) const;
0037
0038 public:
0039    // constructors
0040    TGeoTorus();
0041    TGeoTorus(Double_t r, Double_t rmin, Double_t rmax, Double_t phi1 = 0, Double_t dphi = 360);
0042    TGeoTorus(const char *name, Double_t r, Double_t rmin, Double_t rmax, Double_t phi1 = 0, Double_t dphi = 360);
0043    TGeoTorus(Double_t *params);
0044    // destructor
0045    ~TGeoTorus() override {}
0046    // methods
0047
0048    void ComputeBBox() override;
0049    void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) override;
0050    void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize) override;
0051    Bool_t Contains(const Double_t *point) const override;
0052    void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const override;
0053    Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact = 1, Double_t step = TGeoShape::Big(),
0054                            Double_t *safe = nullptr) const override;
0055    void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0056                          Double_t *step) const override;
0057    Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0058                             Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const override;
0059    void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0060                           Double_t *step) const override;
0061    Int_t DistancetoPrimitive(Int_t px, Int_t py) override;
0062    TGeoVolume *
0063    Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) override;
0064    const char *GetAxisName(Int_t iaxis) const override;
0065    Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const override;
0066    void GetBoundingCylinder(Double_t *param) const override;
0067    const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const override;
0068    Int_t GetByteCount() const override { return 56; }
0069    TGeoShape *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const override;
0070    void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const override;
0071    Int_t GetNmeshVertices() const override;
0072    Bool_t GetPointsOnSegments(Int_t /*npoints*/, Double_t * /*array*/) const override { return kFALSE; }
0073    Double_t GetR() const { return fR; }
0074    Double_t GetRmin() const { return fRmin; }
0075    Double_t GetRmax() const { return fRmax; }
0076    Double_t GetPhi1() const { return fPhi1; }
0077    Double_t GetDphi() const { return fDphi; }
0078    void InspectShape() const override;
0079    Bool_t IsCylType() const override { return kTRUE; }
0080    TBuffer3D *MakeBuffer3D() const override;
0081    Double_t Safety(const Double_t *point, Bool_t in = kTRUE) const override;
0082    void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const override;
0083    void SavePrimitive(std::ostream &out, Option_t *option = "") override;
0084    void SetTorusDimensions(Double_t r, Double_t rmin, Double_t rmax, Double_t phi1, Double_t dphi);
0085    void SetDimensions(Double_t *param) override;
0086    void SetPoints(Double_t *points) const override;
0087    void SetPoints(Float_t *points) const override;
0088    void SetSegsAndPols(TBuffer3D &buff) const override;
0089    void Sizeof3D() const override;
0090
0091    ClassDefOverride(TGeoTorus, 1) // torus class
0092 };
0093
0094 #endif