Back to home page

EIC code displayed by LXR
 
 

    

File indexing completed on 2026-09-18 09:30:09 

0001 // @(#)root/base:$Id$
0002 // Author: Andrei Gheata   24/10/01
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_TGeoTube
0013 #define ROOT_TGeoTube
0014 
0015 #include "TGeoBBox.h"
0016 
0017 class TGeoTube : public TGeoBBox {
0018 protected:
0019    // data members
0020    Double_t fRmin; // inner radius
0021    Double_t fRmax; // outer radius
0022    Double_t fDz;   // half length
0023                    // methods
0024 
0025    TGeoTube(const TGeoTube &) = delete;
0026    TGeoTube &operator=(const TGeoTube &) = delete;
0027 
0028 public:
0029    // constructors
0030    TGeoTube();
0031    TGeoTube(Double_t rmin, Double_t rmax, Double_t dz);
0032    TGeoTube(const char *name, Double_t rmin, Double_t rmax, Double_t dz);
0033    TGeoTube(Double_t *params);
0034    // destructor
0035    ~TGeoTube() override;
0036    // methods
0037 
0038    Double_t Capacity() const override;
0039    static Double_t Capacity(Double_t rmin, Double_t rmax, Double_t dz);
0040    void ComputeBBox() override;
0041    void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) const override;
0042    void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize) override;
0043    static void ComputeNormalS(const Double_t *point, const Double_t *dir, Double_t *norm, Double_t rmin, Double_t rmax,
0044                               Double_t dz);
0045    Bool_t Contains(const Double_t *point) const override;
0046    void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const override;
0047    static Double_t
0048    DistFromInsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz);
0049    Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact = 1, Double_t step = TGeoShape::Big(),
0050                            Double_t *safe = nullptr) const override;
0051    void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0052                          Double_t *step) const override;
0053    static Double_t
0054    DistFromOutsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax, Double_t dz);
0055    Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0056                             Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const override;
0057    void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0058                           Double_t *step) const override;
0059    static void DistToTube(Double_t rsq, Double_t nsq, Double_t rdotn, Double_t radius, Double_t &b, Double_t &delta);
0060    Int_t DistancetoPrimitive(Int_t px, Int_t py) override;
0061    TGeoVolume *
0062    Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) override;
0063    const char *GetAxisName(Int_t iaxis) const override;
0064    Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const override;
0065    void GetBoundingCylinder(Double_t *param) const override;
0066    const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const override;
0067    Int_t GetByteCount() const override { return 48; }
0068    Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const override;
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    virtual Double_t GetRmin() const { return fRmin; }
0073    virtual Double_t GetRmax() const { return fRmax; }
0074    virtual Double_t GetDz() const { return fDz; }
0075    inline Bool_t HasRmin() const { return (fRmin > 0) ? kTRUE : kFALSE; }
0076    void InspectShape() const override;
0077    Bool_t IsConvex() const override { return !HasRmin(); }
0078    Bool_t IsCylType() const override { return kTRUE; }
0079    Double_t Safety(const Double_t *point, Bool_t in = kTRUE) const override;
0080    void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const override;
0081    static Double_t
0082    SafetyS(const Double_t *point, Bool_t in, Double_t rmin, Double_t rmax, Double_t dz, Int_t skipz = 0);
0083    void SavePrimitive(std::ostream &out, Option_t *option = "") override;
0084    void SetTubeDimensions(Double_t rmin, Double_t rmax, Double_t dz);
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(TGeoTube, 1) // cylindrical tube class
0092 };
0093 
0094 class TGeoTubeSeg : public TGeoTube {
0095 protected:
0096    // data members
0097    Double_t fPhi1; // first phi limit
0098    Double_t fPhi2; // second phi limit
0099    // Transient trigonometric data
0100    Double_t fS1;   // sin(phi1)
0101    Double_t fC1;   // cos(phi1)
0102    Double_t fS2;   // sin(phi2)
0103    Double_t fC2;   // cos(phi2)
0104    Double_t fSm;   // sin(0.5*(phi1+phi2))
0105    Double_t fCm;   // cos(0.5*(phi1+phi2))
0106    Double_t fCdfi; // cos(0.5*(phi1-phi2))
0107 
0108    void InitTrigonometry();
0109 
0110 public:
0111    // constructors
0112    TGeoTubeSeg();
0113    TGeoTubeSeg(Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2);
0114    TGeoTubeSeg(const char *name, Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2);
0115    TGeoTubeSeg(Double_t *params);
0116    // destructor
0117    ~TGeoTubeSeg() override;
0118    // methods
0119    void AfterStreamer() override;
0120    Double_t Capacity() const override;
0121    static Double_t Capacity(Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2);
0122    void ComputeBBox() override;
0123    void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) const override;
0124    void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize) override;
0125    static void ComputeNormalS(const Double_t *point, const Double_t *dir, Double_t *norm, Double_t rmin, Double_t rmax,
0126                               Double_t dz, Double_t c1, Double_t s1, Double_t c2, Double_t s2);
0127    Bool_t Contains(const Double_t *point) const override;
0128    void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const override;
0129    Int_t DistancetoPrimitive(Int_t px, Int_t py) override;
0130    static Double_t DistFromInsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax,
0131                                    Double_t dz, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm,
0132                                    Double_t sm, Double_t cdfi);
0133    Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact = 1, Double_t step = TGeoShape::Big(),
0134                            Double_t *safe = nullptr) const override;
0135    void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0136                          Double_t *step) const override;
0137    static Double_t DistFromOutsideS(const Double_t *point, const Double_t *dir, Double_t rmin, Double_t rmax,
0138                                     Double_t dz, Double_t c1, Double_t s1, Double_t c2, Double_t s2, Double_t cm,
0139                                     Double_t sm, Double_t cdfi);
0140    Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0141                             Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const override;
0142    void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0143                           Double_t *step) const override;
0144    TGeoVolume *
0145    Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) override;
0146    Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const override;
0147    void GetBoundingCylinder(Double_t *param) const override;
0148    const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const override;
0149    Int_t GetByteCount() const override { return 56; }
0150    Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const override;
0151    TGeoShape *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const override;
0152    Int_t GetNmeshVertices() const override;
0153    void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const override;
0154    Double_t GetPhi1() const { return fPhi1; }
0155    Double_t GetPhi2() const { return fPhi2; }
0156    inline Bool_t HasDphi() const { return (fPhi2 - fPhi1) < 360.; }
0157    void InspectShape() const override;
0158    Bool_t IsConvex() const final { return !HasRmin() && !HasDphi(); }
0159    Double_t Safety(const Double_t *point, Bool_t in = kTRUE) const override;
0160    void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const override;
0161    static Double_t SafetyS(const Double_t *point, Bool_t in, Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1,
0162                            Double_t phi2, Int_t skipz = 0);
0163    void SavePrimitive(std::ostream &out, Option_t *option = "") override;
0164    void SetTubsDimensions(Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2);
0165    void SetDimensions(Double_t *param) override;
0166    void SetPoints(Double_t *points) const override;
0167    void SetPoints(Float_t *points) const override;
0168    void SetSegsAndPols(TBuffer3D &buff) const override;
0169    void Sizeof3D() const override;
0170 
0171    ClassDefOverride(TGeoTubeSeg, 2) // cylindrical tube segment class
0172 };
0173 
0174 class TGeoCtub : public TGeoTubeSeg {
0175 protected:
0176    // data members
0177    Double_t fNlow[3];  // normal to lower cut plane
0178    Double_t fNhigh[3]; // normal to higher cut plane
0179 
0180 public:
0181    // constructors
0182    TGeoCtub();
0183    TGeoCtub(Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2, Double_t lx, Double_t ly,
0184             Double_t lz, Double_t tx, Double_t ty, Double_t tz);
0185    TGeoCtub(const char *name, Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2, Double_t lx,
0186             Double_t ly, Double_t lz, Double_t tx, Double_t ty, Double_t tz);
0187    TGeoCtub(Double_t *params);
0188    // destructor
0189    ~TGeoCtub() override;
0190    // methods
0191    Double_t Capacity() const override;
0192    void ComputeBBox() override;
0193    void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) const override;
0194    void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize) override;
0195    Bool_t Contains(const Double_t *point) const override;
0196    void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const override;
0197    Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact = 1, Double_t step = TGeoShape::Big(),
0198                            Double_t *safe = nullptr) const override;
0199    void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0200                          Double_t *step) const override;
0201    Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0202                             Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const override;
0203    void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize,
0204                           Double_t *step) const override;
0205    TGeoVolume *
0206    Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) override;
0207    Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const override;
0208    const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const override;
0209    Int_t GetByteCount() const override { return 98; }
0210    Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const override;
0211    TGeoShape *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const override;
0212    const Double_t *GetNlow() const { return &fNlow[0]; }
0213    const Double_t *GetNhigh() const { return &fNhigh[0]; }
0214    Double_t GetZcoord(Double_t xc, Double_t yc, Double_t zc) const;
0215    void InspectShape() const override;
0216    Double_t Safety(const Double_t *point, Bool_t in = kTRUE) const override;
0217    void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const override;
0218    void SavePrimitive(std::ostream &out, Option_t *option = "") override;
0219    void SetCtubDimensions(Double_t rmin, Double_t rmax, Double_t dz, Double_t phi1, Double_t phi2, Double_t lx,
0220                           Double_t ly, Double_t lz, Double_t tx, Double_t ty, Double_t tz);
0221    void SetDimensions(Double_t *param) override;
0222    void SetPoints(Double_t *points) const override;
0223    void SetPoints(Float_t *points) const override;
0224 
0225    ClassDefOverride(TGeoCtub, 1) // cut tube segment class
0226 };
0227 
0228 #endif
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!