ROOT  6.06/09
Reference Guide
TGeoPgon.h
Go to the documentation of this file.
1 // @(#)root/geom:$Id$
2 // Author: Andrei Gheata 31/01/02
3 
4 /*************************************************************************
5  * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6  * All rights reserved. *
7  * *
8  * For the licensing terms see $ROOTSYS/LICENSE. *
9  * For the list of contributors see $ROOTSYS/README/CREDITS. *
10  *************************************************************************/
11 
12 #ifndef ROOT_TGeoPgon
13 #define ROOT_TGeoPgon
14 
15 #ifndef ROOT_TGeoPcon
16 #include "TGeoPcon.h"
17 #endif
18 
19 
20 ////////////////////////////////////////////////////////////////////////////
21 // //
22 // TGeoPgon - a polygone. It has at least 10 parameters : //
23 // - the lower phi limit; //
24 // - the range in phi; //
25 // - the number of edges on each z plane; //
26 // - the number of z planes (at least two) where the inner/outer //
27 // radii are changing; //
28 // - z coordinate, inner and outer radius for each z plane //
29 // //
30 ////////////////////////////////////////////////////////////////////////////
31 
32 class TGeoPgon : public TGeoPcon
33 {
34 public:
35  struct ThreadData_t
36  {
37  Int_t *fIntBuffer; //![fNedges+4] temporary int buffer array
38  Double_t *fDblBuffer; //![fNedges+4] temporary double buffer array
39 
40  ThreadData_t();
41  ~ThreadData_t();
42  };
43  ThreadData_t& GetThreadData() const;
44  void ClearThreadData() const;
45  void CreateThreadData(Int_t nthreads);
46 
47 protected:
48  // data members
49  Int_t fNedges; // number of edges (at least one)
50  mutable std::vector<ThreadData_t*> fThreadData; //! Navigation data per thread
51  mutable Int_t fThreadSize; //! Size for the navigation data array
52 
53  // internal utility methods
54  Int_t GetPhiCrossList(const Double_t *point, const Double_t *dir, Int_t istart, Double_t *sphi, Int_t *iphi, Double_t stepmax=TGeoShape::Big()) const;
55  Bool_t IsCrossingSlice(const Double_t *point, const Double_t *dir, Int_t iphi, Double_t sstart, Int_t &ipl, Double_t &snext, Double_t stepmax) const;
56  void LocatePhi(const Double_t *point, Int_t &ipsec) const;
57  Double_t Rpg(Double_t z, Int_t ipl, Bool_t inner, Double_t &a, Double_t &b) const;
58  Double_t Rproj(Double_t z,const Double_t *point, const Double_t *dir, Double_t cphi, Double_t sphi, Double_t &a, Double_t &b) const;
59  Bool_t SliceCrossing(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const;
60  Bool_t SliceCrossingIn(const Double_t *point, const Double_t *dir, Int_t ipl, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const;
61  Bool_t SliceCrossingZ(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const;
62  Bool_t SliceCrossingInZ(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const;
63 
64 public:
65  // constructors
66  TGeoPgon();
67  TGeoPgon(Double_t phi, Double_t dphi, Int_t nedges, Int_t nz);
68  TGeoPgon(const char *name, Double_t phi, Double_t dphi, Int_t nedges, Int_t nz);
69  TGeoPgon(Double_t *params);
70  // destructor
71  virtual ~TGeoPgon();
72  // methods
73  virtual Double_t Capacity() const;
74  virtual void ComputeBBox();
75  virtual void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm);
76  virtual void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize);
77  virtual Bool_t Contains(const Double_t *point) const;
78  virtual void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const;
79  virtual Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact=1,
80  Double_t step=TGeoShape::Big(), Double_t *safe=0) const;
81  virtual void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const;
82  virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1,
83  Double_t step=TGeoShape::Big(), Double_t *safe=0) const;
84  virtual void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const;
85  virtual Int_t DistancetoPrimitive(Int_t px, Int_t py);
86  virtual TGeoVolume *Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv,
87  Double_t start, Double_t step);
88  virtual void GetBoundingCylinder(Double_t *param) const;
89  virtual const TBuffer3D &GetBuffer3D(Int_t reqSections, Bool_t localFrame) const;
90  virtual Int_t GetByteCount() const {return 64+12*fNz;}
91  virtual TGeoShape *GetMakeRuntimeShape(TGeoShape * /*mother*/, TGeoMatrix * /*mat*/) const {return 0;}
92  virtual void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const;
93  Int_t GetNedges() const {return fNedges;}
94  virtual Int_t GetNmeshVertices() const;
95  virtual Int_t GetNsegments() const {return fNedges;}
96  virtual Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const {return TGeoBBox::GetPointsOnSegments(npoints,array);}
97  virtual void InspectShape() const;
98  virtual TBuffer3D *MakeBuffer3D() const;
99  virtual Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const;
100  virtual void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const;
101  Double_t SafetyToSegment(const Double_t *point, Int_t ipl, Int_t iphi, Bool_t in, Double_t safphi, Double_t safmin=TGeoShape::Big()) const;
102  virtual void SavePrimitive(std::ostream &out, Option_t *option = "");
103  virtual void SetDimensions(Double_t *param);
104  void SetNedges(Int_t ne) {if (ne>2) fNedges=ne;}
105  virtual void SetPoints(Double_t *points) const;
106  virtual void SetPoints(Float_t *points) const;
107  virtual void SetSegsAndPols(TBuffer3D &buff) const;
108  virtual void Sizeof3D() const;
109 
110  ClassDef(TGeoPgon, 1) // polygone class
111 };
112 
113 #endif
void ClearThreadData() const
Definition: TGeoPgon.cxx:78
Bool_t SliceCrossing(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const
Check boundary crossing inside phi slices.
Definition: TGeoPgon.cxx:880
#define snext(osub1, osub2)
Definition: triangle.c:1167
virtual Double_t DistFromInside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const
compute distance from inside point to surface of the polygone first find out in which Z section the p...
Definition: TGeoPgon.cxx:410
virtual Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const
Fills array with n random points located on the line segments of the shape mesh.
Definition: TGeoBBox.cxx:647
virtual void ComputeNormal_v(const Double_t *points, const Double_t *dirs, Double_t *norms, Int_t vecsize)
Compute the normal for an array o points so that norm.dot.dir is positive Input: Arrays of point coor...
Definition: TGeoPgon.cxx:1989
virtual void SavePrimitive(std::ostream &out, Option_t *option="")
Save a primitive as a C++ statement(s) on output stream "out".
Definition: TGeoPgon.cxx:1792
virtual Bool_t Contains(const Double_t *point) const
test if point is inside this shape check total z range
Definition: TGeoPgon.cxx:363
virtual Int_t GetNmeshVertices() const
Return number of vertices of the mesh representation.
Definition: TGeoPgon.cxx:1913
float Float_t
Definition: RtypesCore.h:53
const char Option_t
Definition: RtypesCore.h:62
virtual ~TGeoPgon()
destructor
Definition: TGeoPgon.cxx:170
TGeoPgon()
dummy ctor
Definition: TGeoPgon.cxx:114
Double_t * fDblBuffer
[fNedges+4] temporary int buffer array
Definition: TGeoPgon.h:38
int Int_t
Definition: RtypesCore.h:41
bool Bool_t
Definition: RtypesCore.h:59
TArc * a
Definition: textangle.C:12
Int_t fThreadSize
Navigation data per thread.
Definition: TGeoPgon.h:51
~ThreadData_t()
Destructor.
Definition: TGeoPgon.cxx:62
Double_t SafetyToSegment(const Double_t *point, Int_t ipl, Int_t iphi, Bool_t in, Double_t safphi, Double_t safmin=TGeoShape::Big()) const
Compute safety from POINT to segment between planes ipl, ipl+1 within safmin.
Definition: TGeoPgon.cxx:1648
std::vector< ThreadData_t * > fThreadData
Definition: TGeoPgon.h:50
Bool_t SliceCrossingInZ(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const
Performs ray propagation between Z segments.
Definition: TGeoPgon.cxx:593
#define ClassDef(name, id)
Definition: Rtypes.h:254
Bool_t IsCrossingSlice(const Double_t *point, const Double_t *dir, Int_t iphi, Double_t sstart, Int_t &ipl, Double_t &snext, Double_t stepmax) const
Check crossing of a given pgon slice, from a starting point inside the slice.
Definition: TGeoPgon.cxx:968
Int_t GetNedges() const
Definition: TGeoPgon.h:93
virtual TGeoVolume * Divide(TGeoVolume *voldiv, const char *divname, Int_t iaxis, Int_t ndiv, Double_t start, Double_t step)
— Divide this polygone shape belonging to volume "voldiv" into ndiv volumes called divname...
Definition: TGeoPgon.cxx:1270
Bool_t SliceCrossingZ(const Double_t *point, const Double_t *dir, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const
Performs ray propagation between Z segments.
Definition: TGeoPgon.cxx:665
virtual void ComputeBBox()
compute bounding box for a polygone Check if the sections are in increasing Z order ...
Definition: TGeoPgon.cxx:202
virtual void DistFromOutside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const
Compute distance from array of input points having directions specisied by dirs. Store output in dist...
Definition: TGeoPgon.cxx:2005
virtual void ComputeNormal(const Double_t *point, const Double_t *dir, Double_t *norm)
Compute normal to closest surface from POINT.
Definition: TGeoPgon.cxx:269
virtual void Safety_v(const Double_t *points, const Bool_t *inside, Double_t *safe, Int_t vecsize) const
Compute safe distance from each of the points in the input array.
Definition: TGeoPgon.cxx:2015
char * out
Definition: TBase64.cxx:29
ThreadData_t & GetThreadData() const
Definition: TGeoPgon.cxx:70
virtual Double_t Safety(const Double_t *point, Bool_t in=kTRUE) const
computes the closest distance from given point to this shape, according to option.
Definition: TGeoPgon.cxx:1717
point * points
Definition: X3DBuffer.c:20
virtual void InspectShape() const
Inspect the PGON parameters.
Definition: TGeoPgon.cxx:1381
virtual TBuffer3D * MakeBuffer3D() const
Creates a TBuffer3D describing this shape.
Definition: TGeoPgon.cxx:1392
virtual void GetMeshNumbers(Int_t &nvert, Int_t &nsegs, Int_t &npols) const
Returns numbers of vertices, segments and polygons composing the shape mesh.
Definition: TGeoPgon.cxx:1900
Int_t fNz
Definition: TGeoPcon.h:35
Generic 3D primitive description class.
Definition: TBuffer3D.h:19
Int_t GetPhiCrossList(const Double_t *point, const Double_t *dir, Int_t istart, Double_t *sphi, Int_t *iphi, Double_t stepmax=TGeoShape::Big()) const
Size for the navigation data array.
Definition: TGeoPgon.cxx:517
Double_t Rpg(Double_t z, Int_t ipl, Bool_t inner, Double_t &a, Double_t &b) const
Computes projected pgon radius (inner or outer) corresponding to a given Z value. ...
Definition: TGeoPgon.cxx:1599
void CreateThreadData(Int_t nthreads)
Create thread data for n threads max.
Definition: TGeoPgon.cxx:95
ThreadData_t()
[fNedges+4] temporary double buffer array
Double_t Rproj(Double_t z, const Double_t *point, const Double_t *dir, Double_t cphi, Double_t sphi, Double_t &a, Double_t &b) const
Computes projected distance at a given Z for a given ray inside a given sector and fills coefficients...
Definition: TGeoPgon.cxx:1633
virtual Double_t Capacity() const
Computes capacity of the shape in [length^3].
Definition: TGeoPgon.cxx:178
Int_t fNedges
Definition: TGeoPgon.h:49
double Double_t
Definition: RtypesCore.h:55
virtual void SetPoints(Double_t *points) const
create polygone mesh points
Definition: TGeoPgon.cxx:1840
virtual Int_t GetByteCount() const
Definition: TGeoPgon.h:90
virtual void SetSegsAndPols(TBuffer3D &buff) const
Fill TBuffer3D structure for segments and polygons.
Definition: TGeoPgon.cxx:1419
static Double_t Big()
Definition: TGeoShape.h:98
#define name(a, b)
Definition: linkTestLib0.cpp:5
virtual void GetBoundingCylinder(Double_t *param) const
— Fill vector param[4] with the bounding cylinder parameters.
Definition: TGeoPgon.cxx:1357
virtual TGeoShape * GetMakeRuntimeShape(TGeoShape *, TGeoMatrix *) const
In case shape has some negative parameters, these has to be computed in order to fit the mother...
Definition: TGeoPgon.h:91
virtual void SetDimensions(Double_t *param)
Set PGON dimensions starting from an array.
Definition: TGeoPgon.cxx:1814
virtual Int_t DistancetoPrimitive(Int_t px, Int_t py)
compute closest distance from point px,py to each corner
Definition: TGeoPgon.cxx:1255
virtual void DistFromInside_v(const Double_t *points, const Double_t *dirs, Double_t *dists, Int_t vecsize, Double_t *step) const
Compute distance from array of input points having directions specisied by dirs. Store output in dist...
Definition: TGeoPgon.cxx:1997
Bool_t SliceCrossingIn(const Double_t *point, const Double_t *dir, Int_t ipl, Int_t nphi, Int_t *iphi, Double_t *sphi, Double_t &snext, Double_t stepmax) const
Check boundary crossing inside phi slices.
Definition: TGeoPgon.cxx:745
void SetNedges(Int_t ne)
Definition: TGeoPgon.h:104
virtual void Contains_v(const Double_t *points, Bool_t *inside, Int_t vecsize) const
Check the inside status for each of the points in the array.
Definition: TGeoPgon.cxx:1979
virtual Double_t DistFromOutside(const Double_t *point, const Double_t *dir, Int_t iact=1, Double_t step=TGeoShape::Big(), Double_t *safe=0) const
Compute distance from outside point to surface of the polygone.
Definition: TGeoPgon.cxx:1067
const Bool_t kTRUE
Definition: Rtypes.h:91
virtual const TBuffer3D & GetBuffer3D(Int_t reqSections, Bool_t localFrame) const
Fills a static 3D buffer and returns a reference.
Definition: TGeoPgon.cxx:1940
virtual Bool_t GetPointsOnSegments(Int_t npoints, Double_t *array) const
Fills array with n random points located on the line segments of the shape mesh.
Definition: TGeoPgon.h:96
double norm(double *x, double *p)
Definition: unuranDistr.cxx:40
virtual void Sizeof3D() const
Definition: TGeoPgon.cxx:1933
void LocatePhi(const Double_t *point, Int_t &ipsec) const
Locates index IPSEC of the phi sector containing POINT.
Definition: TGeoPgon.cxx:506
virtual Int_t GetNsegments() const
Returns number of segments on each mesh circle segment.
Definition: TGeoPgon.h:95