Warning, file /include/root/TGeoShape.h was not indexed
or was modified since last indexation (in which case cross-reference links may be missing, inaccurate or erroneous).
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012 #ifndef ROOT_TGeoShape
0013 #define ROOT_TGeoShape
0014
0015 #include "TNamed.h"
0016
0017
0018 class TGeoBoolCombinator;
0019 class TGeoBBox;
0020 class TGeoMatrix;
0021 class TGeoHMatrix;
0022 class TGeoVolume;
0023 class TBuffer3D;
0024
0025 class TGeoShape : public TNamed {
0026 private:
0027 static TGeoMatrix *fgTransform;
0028 static Double_t fgEpsMch;
0029 public:
0030 enum EShapeType {
0031 kBitMask32 = 0xffffffff,
0032 kGeoNoShape = 0,
0033 kGeoBad = BIT(0),
0034 kGeoRSeg = BIT(1),
0035 kGeoPhiSeg = BIT(2),
0036 kGeoThetaSeg = BIT(3),
0037 kGeoVisX = BIT(4),
0038 kGeoVisY = BIT(5),
0039 kGeoVisZ = BIT(6),
0040 kGeoRunTimeShape = BIT(7),
0041 kGeoInvalidShape = BIT(8),
0042 kGeoTorus = BIT(9),
0043 kGeoBox = BIT(10),
0044 kGeoPara = BIT(11),
0045 kGeoSph = BIT(12),
0046 kGeoTube = BIT(13),
0047 kGeoTubeSeg = BIT(14),
0048 kGeoCone = BIT(15),
0049 kGeoConeSeg = BIT(16),
0050 kGeoPcon = BIT(17),
0051 kGeoPgon = BIT(18),
0052 kGeoArb8 = BIT(19),
0053 kGeoEltu = BIT(20),
0054 kGeoTrap = BIT(21),
0055 kGeoCtub = BIT(22),
0056 kGeoTrd1 = BIT(23),
0057 kGeoTrd2 = BIT(24),
0058 kGeoComb = BIT(25),
0059 kGeoClosedShape = BIT(26),
0060 kGeoXtru = BIT(27),
0061 kGeoParaboloid = BIT(28),
0062 kGeoHalfSpace = BIT(29),
0063 kGeoHype = BIT(30),
0064 kGeoSavePrimitive = BIT(20)
0065 };
0066 virtual void ClearThreadData() const {}
0067 virtual void CreateThreadData(Int_t) {}
0068
0069 protected:
0070
0071 Int_t fShapeId;
0072 UInt_t fShapeBits;
0073
0074 virtual void FillBuffer3D(TBuffer3D &buffer, Int_t reqSections, Bool_t localFrame) const;
0075 Int_t GetBasicColor() const;
0076 void SetOnBoundary(Bool_t ) {}
0077 void TransformPoints(Double_t *points, UInt_t NbPoints) const;
0078
0079 public:
0080
0081 TGeoShape();
0082 TGeoShape(const char *name);
0083
0084 ~TGeoShape() override;
0085
0086
0087 static Double_t Big() { return 1.E30; }
0088 static TGeoMatrix *GetTransform();
0089 static void SetTransform(TGeoMatrix *matrix);
0090 static Double_t Tolerance() { return 1.E-10; }
0091 static Double_t ComputeEpsMch();
0092 static Double_t EpsMch();
0093 virtual void AfterStreamer(){};
0094 virtual Double_t Capacity() const = 0;
0095 void CheckShape(Int_t testNo, Int_t nsamples = 10000, Option_t *option = "");
0096 virtual void ComputeBBox() = 0;
0097 virtual void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm) = 0;
0098 virtual void ComputeNormal_v(const Double_t *, const Double_t *, Double_t *, Int_t) {}
0099 virtual Bool_t Contains(const Double_t *point) const = 0;
0100 virtual void Contains_v(const Double_t *, Bool_t *, Int_t) const {}
0101 virtual Bool_t CouldBeCrossed(const Double_t *point, const Double_t *dir) const = 0;
0102 Int_t DistancetoPrimitive(Int_t px, Int_t py) override = 0;
0103 virtual Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0104 Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const = 0;
0105 virtual void DistFromInside_v(const Double_t *, const Double_t *, Double_t *, Int_t, Double_t *) const {}
0106 virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact = 1,
0107 Double_t step = TGeoShape::Big(), Double_t *safe = nullptr) const = 0;
0108 virtual void DistFromOutside_v(const Double_t *, const Double_t *, Double_t *, Int_t, Double_t *) const {}
0109 static Double_t DistToPhiMin(const Double_t *point, const Double_t *dir, Double_t s1, Double_t c1, Double_t s2,
0110 Double_t c2, Double_t sm, Double_t cm, Bool_t in = kTRUE);
0111 virtual TGeoVolume *
0112 Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step) = 0;
0113 void Draw(Option_t *option = "") override;
0114 void ExecuteEvent(Int_t event, Int_t px, Int_t py) override;
0115 virtual const char *GetAxisName(Int_t iaxis) const = 0;
0116 virtual Double_t GetAxisRange(Int_t iaxis, Double_t &xlo, Double_t &xhi) const = 0;
0117 virtual void GetBoundingCylinder(Double_t *param) const = 0;
0118 virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const;
0119 virtual Int_t GetByteCount() const = 0;
0120 virtual Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const = 0;
0121 virtual Int_t
0122 GetFittingBox(const TGeoBBox *parambox, TGeoMatrix *mat, Double_t &dx, Double_t &dy, Double_t &dz) const = 0;
0123 Int_t GetId() const { return fShapeId; }
0124 virtual TGeoShape *GetMakeRuntimeShape(TGeoShape *mother, TGeoMatrix *mat) const = 0;
0125 virtual void GetMeshNumbers(Int_t & , Int_t & , Int_t & ) const {}
0126 const char *GetName() const override;
0127 virtual Int_t GetNmeshVertices() const { return 0; }
0128 const char *GetPointerName() const;
0129 virtual Bool_t IsAssembly() const { return kFALSE; }
0130 virtual Bool_t IsComposite() const { return kFALSE; }
0131 virtual Bool_t IsCylType() const = 0;
0132 static Bool_t
0133 IsCloseToPhi(Double_t epsil, const Double_t *point, Double_t c1, Double_t s1, Double_t c2, Double_t s2);
0134 static Bool_t IsCrossingSemiplane(const Double_t *point, const Double_t *dir, Double_t cphi, Double_t sphi,
0135 Double_t &snext, Double_t &rxy);
0136 static Bool_t IsSameWithinTolerance(Double_t a, Double_t b);
0137 static Bool_t IsSegCrossing(Double_t x1, Double_t y1, Double_t x2, Double_t y2, Double_t x3, Double_t y3,
0138 Double_t x4, Double_t y4);
0139 static Bool_t IsInPhiRange(const Double_t *point, Double_t phi1, Double_t phi2);
0140 virtual Bool_t IsReflected() const { return kFALSE; }
0141 virtual Bool_t IsVecGeom() const { return kFALSE; }
0142 Bool_t IsRunTimeShape() const { return TestShapeBit(kGeoRunTimeShape); }
0143 Bool_t IsValid() const { return !TestShapeBit(kGeoInvalidShape); }
0144 virtual Bool_t IsValidBox() const = 0;
0145 virtual void InspectShape() const = 0;
0146 virtual TBuffer3D *MakeBuffer3D() const { return nullptr; }
0147 static void NormalPhi(const Double_t *point, const Double_t *dir, Double_t *norm, Double_t c1, Double_t s1,
0148 Double_t c2, Double_t s2);
0149 void Paint(Option_t *option = "") override;
0150 virtual Double_t Safety(const Double_t *point, Bool_t in = kTRUE) const = 0;
0151 virtual void Safety_v(const Double_t *, const Bool_t *, Double_t *, Int_t) const {}
0152 static Double_t SafetyPhi(const Double_t *point, Bool_t in, Double_t phi1, Double_t phi2);
0153 static Double_t SafetySeg(Double_t r, Double_t z, Double_t r1, Double_t z1, Double_t r2, Double_t z2, Bool_t outer);
0154 virtual void SetDimensions(Double_t *param) = 0;
0155 void SetId(Int_t id) { fShapeId = id; }
0156 virtual void SetPoints(Double_t *points) const = 0;
0157 virtual void SetPoints(Float_t *points) const = 0;
0158 virtual void SetSegsAndPols(TBuffer3D &buff) const = 0;
0159 void SetRuntime(Bool_t flag = kTRUE) { SetShapeBit(kGeoRunTimeShape, flag); }
0160 Int_t ShapeDistancetoPrimitive(Int_t numpoints, Int_t px, Int_t py) const;
0161 virtual void Sizeof3D() const = 0;
0162
0163
0164 void SetShapeBit(UInt_t f, Bool_t set);
0165 void SetShapeBit(UInt_t f) { fShapeBits |= f & kBitMask32; }
0166 void ResetShapeBit(UInt_t f) { fShapeBits &= ~(f & kBitMask32); }
0167 Bool_t TestShapeBit(UInt_t f) const { return (Bool_t)((fShapeBits & f) != 0); }
0168 Int_t TestShapeBits(UInt_t f) const { return (Int_t)(fShapeBits & f); }
0169 void InvertShapeBit(UInt_t f) { fShapeBits ^= f & kBitMask32; }
0170
0171 ClassDefOverride(TGeoShape, 2)
0172 };
0173
0174 #endif