Geant4  10.02.p01
G4VTwistedFaceted.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // $Id: G4VTwistedFaceted.hh 83572 2014-09-01 15:23:27Z gcosmo $
28 //
29 // --------------------------------------------------------------------
30 // GEANT 4 class header file
31 //
32 //
33 // G4VTwistedFaceted
34 //
35 // Class description:
36 //
37 // G4VTwistedFaceted is an abstract base class for twisted boxoids:
38 // G4TwistedTrd, G4TwistedTrap and G4TwistedBox
39 
40 // Author:
41 //
42 // 27-Oct-2004 - O.Link (Oliver.Link@cern.ch)
43 //
44 // --------------------------------------------------------------------
45 
46 #ifndef __G4VTWISTEDFACETED__
47 #define __G4VTWISTEDFACETED__
48 
49 #include "G4VSolid.hh"
50 #include "G4TwistTrapAlphaSide.hh"
52 #include "G4TwistBoxSide.hh"
53 #include "G4TwistTrapFlatSide.hh"
54 
55 class G4SolidExtentList;
56 class G4ClippablePolygon;
57 
59 {
60  public: // with description
61 
62  G4VTwistedFaceted(const G4String &pname, // Name of instance
63  G4double PhiTwist, // twist angle
64  G4double pDz, // half z lenght
65  G4double pTheta, // direction between end planes
66  G4double pPhi, // defined by polar & azim. angles
67  G4double pDy1, // half y length at -pDz
68  G4double pDx1, // half x length at -pDz,-pDy
69  G4double pDx2, // half x length at -pDz,+pDy
70  G4double pDy2, // half y length at +pDz
71  G4double pDx3, // half x length at +pDz,-pDy
72  G4double pDx4, // half x length at +pDz,+pDy
73  G4double pAlph // tilt angle at +pDz
74  );
75 
76  virtual ~G4VTwistedFaceted();
77 
79  const G4int,
80  const G4VPhysicalVolume* );
81 
82  virtual G4bool CalculateExtent(const EAxis pAxis,
83  const G4VoxelLimits &pVoxelLimit,
84  const G4AffineTransform &pTransform,
85  G4double &pMin,
86  G4double &pMax ) const;
87 
88  virtual G4double DistanceToIn (const G4ThreeVector &p,
89  const G4ThreeVector &v ) const;
90 
91  virtual G4double DistanceToIn (const G4ThreeVector &p ) const;
92 
93  virtual G4double DistanceToOut(const G4ThreeVector &p,
94  const G4ThreeVector &v,
95  const G4bool calcnorm = false,
96  G4bool *validnorm = 0,
97  G4ThreeVector *n=0 ) const;
98 
99  virtual G4double DistanceToOut(const G4ThreeVector &p) const;
100 
101  virtual EInside Inside (const G4ThreeVector &p) const;
102 
103  virtual G4ThreeVector SurfaceNormal(const G4ThreeVector &p) const;
104 
107 
108  virtual inline G4double GetCubicVolume() ;
109  virtual inline G4double GetSurfaceArea() ;
110 
111  virtual void DescribeYourselfTo (G4VGraphicsScene &scene) const;
112  virtual G4Polyhedron *CreatePolyhedron () const ;
113  virtual G4Polyhedron *GetPolyhedron () const;
114 
115  virtual std::ostream &StreamInfo(std::ostream& os) const;
116 
117  // accessors
118 
119  inline G4double GetTwistAngle () const { return fPhiTwist; }
120 
121  inline G4double GetDx1 () const { return fDx1 ; }
122  inline G4double GetDx2 () const { return fDx2 ; }
123  inline G4double GetDx3 () const { return fDx3 ; }
124  inline G4double GetDx4 () const { return fDx4 ; }
125  inline G4double GetDy1 () const { return fDy1 ; }
126  inline G4double GetDy2 () const { return fDy2 ; }
127  inline G4double GetDz () const { return fDz ; }
128  inline G4double GetPhi () const { return fPhi ; }
129  inline G4double GetTheta () const { return fTheta ; }
130  inline G4double GetAlpha () const { return fAlph ; }
131 
132  inline G4double Xcoef(G4double u,G4double phi, G4double ftg) const ;
133  // For calculating the w(u) function
134 
135  inline G4double GetValueA(G4double phi) const;
136  inline G4double GetValueB(G4double phi) const;
137  inline G4double GetValueD(G4double phi) const;
138 
139  virtual G4VisExtent GetExtent () const;
140  virtual G4GeometryType GetEntityType() const;
141 
142  public: // without description
143 
144  G4VTwistedFaceted(__void__&);
145  // Fake default constructor for usage restricted to direct object
146  // persistency for clients requiring preallocation of memory for
147  // persistifiable objects.
148 
151  // Copy constructor and assignment operator.
152 
153  protected: // with description
154 
156  CreateRotatedVertices(const G4AffineTransform& pTransform) const;
157  // Create the List of transformed vertices in the format required
158  // for G4VSolid:: ClipCrossSection and ClipBetweenSections.
159 
161  mutable G4Polyhedron* fpPolyhedron; // pointer to polyhedron for vis
162 
163  private:
164 
165  void CreateSurfaces();
166 
167  private:
168 
171 
175 
179 
180  G4double fDz; // Half-length along the z axis
181 
182  G4double fDx ; // maximum side in x
183  G4double fDy ; // maximum side in y
184 
186  G4double fTAlph ; // std::tan(fAlph)
187 
190 
191  G4double fPhiTwist; // twist angle ( dphi in surface equation)
192 
193  G4VTwistSurface *fLowerEndcap ; // surface of -ve z
194  G4VTwistSurface *fUpperEndcap ; // surface of +ve z
195 
196  G4VTwistSurface *fSide0 ; // Twisted Side at phi = 0 deg
197  G4VTwistSurface *fSide90 ; // Twisted Side at phi = 90 deg
198  G4VTwistSurface *fSide180 ; // Twisted Side at phi = 180 deg
199  G4VTwistSurface *fSide270 ; // Twisted Side at phi = 270 deg
200 
201  G4double fCubicVolume ; // volume of the solid
202  G4double fSurfaceArea ; // area of the solid
203 
204  class LastState // last Inside result
205  {
206  public:
208  {
210  }
212  LastState(const LastState& r) : p(r.p), inside(r.inside){}
214  {
215  if (this == &r) { return *this; }
216  p = r.p; inside = r.inside;
217  return *this;
218  }
219  public:
222  };
223 
224  class LastVector // last SurfaceNormal result
225  {
226  public:
228  {
231  surface = new G4VTwistSurface*[1];
232  }
234  {
235  delete [] surface;
236  }
237  LastVector(const LastVector& r) : p(r.p), vec(r.vec)
238  {
239  surface = new G4VTwistSurface*[1];
240  surface[0] = r.surface[0];
241  }
243  {
244  if (&r == this) { return *this; }
245  p = r.p; vec = r.vec;
246  delete [] surface; surface = new G4VTwistSurface*[1];
247  surface[0] = r.surface[0];
248  return *this;
249  }
250  public:
254  };
255 
256  class LastValue // last G4double value
257  {
258  public:
260  {
262  value = DBL_MAX;
263  }
265  LastValue(const LastValue& r) : p(r.p), value(r.value){}
267  {
268  if (this == &r) { return *this; }
269  p = r.p; value = r.value;
270  return *this;
271  }
272  public:
275  };
276 
277  class LastValueWithDoubleVector // last G4double value
278  {
279  public:
281  {
284  value = DBL_MAX;
285  }
288  : p(r.p), vec(r.vec), value(r.value){}
290  {
291  if (this == &r) { return *this; }
292  p = r.p; vec = r.vec; value = r.value;
293  return *this;
294  }
295  public:
299  };
300 
307 
308  };
309 
310 //=====================================================================
311 
312 inline
314 {
315  if(fCubicVolume != 0.) ;
316  else fCubicVolume = 2 * fDz
317  * ( ( fDx1 + fDx2 ) * fDy1 + ( fDx3 + fDx4 ) * fDy2 );
318  return fCubicVolume;
319 }
320 
321 inline
323 {
324  if(fSurfaceArea != 0.) ;
326  return fSurfaceArea;
327 }
328 
329 inline
331 {
332  return ( fDx4 + fDx2 + ( fDx4 - fDx2 ) * ( 2 * phi ) / fPhiTwist ) ;
333 }
334 
335 inline
337 {
338  return ( fDx3 + fDx1 + ( fDx3 - fDx1 ) * ( 2 * phi ) / fPhiTwist ) ;
339 }
340 
341 inline
343 {
344  return ( fDy2 + fDy1 + ( fDy2 - fDy1 ) * ( 2 * phi ) / fPhiTwist ) ;
345 }
346 
347 inline
349 {
350  return GetValueA(phi)/2. + (GetValueD(phi)-GetValueA(phi))/4.
351  - u*( ( GetValueD(phi)-GetValueA(phi) ) / ( 2 * GetValueB(phi) ) - ftg );
352 }
353 
354 #endif
virtual void ComputeDimensions(G4VPVParameterisation *, const G4int, const G4VPhysicalVolume *)
LastValueWithDoubleVector(const LastValueWithDoubleVector &r)
virtual G4double GetSurfaceArea()
G4ThreeVectorList * CreateRotatedVertices(const G4AffineTransform &pTransform) const
G4VTwistSurface * fSide90
virtual G4double GetCubicVolume()
static const G4double kInfinity
Definition: geomdefs.hh:42
G4VTwistedFaceted(const G4String &pname, G4double PhiTwist, G4double pDz, G4double pTheta, G4double pPhi, G4double pDy1, G4double pDx1, G4double pDx2, G4double pDy2, G4double pDx3, G4double pDx4, G4double pAlph)
G4VTwistSurface * fUpperEndcap
CLHEP::Hep3Vector G4ThreeVector
LastValueWithDoubleVector fLastDistanceToOutWithV
G4double z
Definition: TRTMaterials.hh:39
virtual G4GeometryType GetEntityType() const
G4ThreeVector GetPointOnSurface() const
G4double GetDy2() const
G4ThreeVector GetPointInSolid(G4double z) const
LastValue & operator=(const LastValue &r)
virtual G4Polyhedron * CreatePolyhedron() const
G4double GetDy1() const
virtual G4double DistanceToIn(const G4ThreeVector &p, const G4ThreeVector &v) const
G4double GetDx3() const
LastValueWithDoubleVector & operator=(const LastValueWithDoubleVector &r)
G4VTwistedFaceted & operator=(const G4VTwistedFaceted &rhs)
G4double GetTwistAngle() const
int G4int
Definition: G4Types.hh:78
virtual G4ThreeVector SurfaceNormal(const G4ThreeVector &p) const
G4VTwistSurface * fSide180
G4double Xcoef(G4double u, G4double phi, G4double ftg) const
G4double GetValueD(G4double phi) const
G4VTwistSurface * fSide270
G4double GetTheta() const
virtual G4bool CalculateExtent(const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pMin, G4double &pMax) const
bool G4bool
Definition: G4Types.hh:79
G4double GetValueB(G4double phi) const
G4VTwistSurface * fLowerEndcap
std::vector< G4ThreeVector > G4ThreeVectorList
Definition: G4VSolid.hh:79
const G4int n
LastValueWithDoubleVector fLastDistanceToInWithV
G4double GetDx2() const
LastVector & operator=(const LastVector &r)
virtual G4double DistanceToOut(const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcnorm=false, G4bool *validnorm=0, G4ThreeVector *n=0) const
G4double GetPhi() const
virtual G4Polyhedron * GetPolyhedron() const
EInside
Definition: geomdefs.hh:58
virtual G4VisExtent GetExtent() const
virtual EInside Inside(const G4ThreeVector &p) const
EAxis
Definition: geomdefs.hh:54
virtual void DescribeYourselfTo(G4VGraphicsScene &scene) const
G4double GetAlpha() const
G4VTwistSurface * fSide0
G4double GetDx1() const
double G4double
Definition: G4Types.hh:76
virtual std::ostream & StreamInfo(std::ostream &os) const
G4double GetDx4() const
#define DBL_MAX
Definition: templates.hh:83
virtual G4double GetSurfaceArea()
Definition: G4VSolid.cc:250
G4Polyhedron * fpPolyhedron
LastState & operator=(const LastState &r)
G4double GetValueA(G4double phi) const
G4double GetDz() const