Geant4  10.00.p02
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 81641 2014-06-04 09:11:38Z 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 
160  mutable G4Polyhedron* fpPolyhedron; // pointer to polyhedron for vis
161 
162  private:
163 
164  void CreateSurfaces();
165 
166  private:
167 
170 
174 
178 
179  G4double fDz; // Half-length along the z axis
180 
181  G4double fDx ; // maximum side in x
182  G4double fDy ; // maximum side in y
183 
185  G4double fTAlph ; // std::tan(fAlph)
186 
189 
190  G4double fPhiTwist; // twist angle ( dphi in surface equation)
191 
192  G4VTwistSurface *fLowerEndcap ; // surface of -ve z
193  G4VTwistSurface *fUpperEndcap ; // surface of +ve z
194 
195  G4VTwistSurface *fSide0 ; // Twisted Side at phi = 0 deg
196  G4VTwistSurface *fSide90 ; // Twisted Side at phi = 90 deg
197  G4VTwistSurface *fSide180 ; // Twisted Side at phi = 180 deg
198  G4VTwistSurface *fSide270 ; // Twisted Side at phi = 270 deg
199 
200  G4double fCubicVolume ; // volume of the solid
201  G4double fSurfaceArea ; // area of the solid
202 
203  class LastState // last Inside result
204  {
205  public:
207  {
209  }
211  LastState(const LastState& r) : p(r.p), inside(r.inside){}
213  {
214  if (this == &r) { return *this; }
215  p = r.p; inside = r.inside;
216  return *this;
217  }
218  public:
221  };
222 
223  class LastVector // last SurfaceNormal result
224  {
225  public:
227  {
230  surface = new G4VTwistSurface*[1];
231  }
233  {
234  delete [] surface;
235  }
236  LastVector(const LastVector& r) : p(r.p), vec(r.vec)
237  {
238  surface = new G4VTwistSurface*[1];
239  surface[0] = r.surface[0];
240  }
242  {
243  if (&r == this) { return *this; }
244  p = r.p; vec = r.vec;
245  delete [] surface; surface = new G4VTwistSurface*[1];
246  surface[0] = r.surface[0];
247  return *this;
248  }
249  public:
253  };
254 
255  class LastValue // last G4double value
256  {
257  public:
259  {
261  value = DBL_MAX;
262  }
264  LastValue(const LastValue& r) : p(r.p), value(r.value){}
266  {
267  if (this == &r) { return *this; }
268  p = r.p; value = r.value;
269  return *this;
270  }
271  public:
274  };
275 
276  class LastValueWithDoubleVector // last G4double value
277  {
278  public:
280  {
283  value = DBL_MAX;
284  }
287  : p(r.p), vec(r.vec), value(r.value){}
289  {
290  if (this == &r) { return *this; }
291  p = r.p; vec = r.vec; value = r.value;
292  return *this;
293  }
294  public:
298  };
299 
306 
307  };
308 
309 //=====================================================================
310 
311 inline
313 {
314  if(fCubicVolume != 0.) ;
315  else fCubicVolume = 2 * fDz
316  * ( ( fDx1 + fDx2 ) * fDy1 + ( fDx3 + fDx4 ) * fDy2 );
317  return fCubicVolume;
318 }
319 
320 inline
322 {
323  if(fSurfaceArea != 0.) ;
325  return fSurfaceArea;
326 }
327 
328 inline
330 {
331  return ( fDx4 + fDx2 + ( fDx4 - fDx2 ) * ( 2 * phi ) / fPhiTwist ) ;
332 }
333 
334 inline
336 {
337  return ( fDx3 + fDx1 + ( fDx3 - fDx1 ) * ( 2 * phi ) / fPhiTwist ) ;
338 }
339 
340 inline
342 {
343  return ( fDy2 + fDy1 + ( fDy2 - fDy1 ) * ( 2 * phi ) / fPhiTwist ) ;
344 }
345 
346 inline
348 {
349  return GetValueA(phi)/2. + (GetValueD(phi)-GetValueA(phi))/4.
350  - u*( ( GetValueD(phi)-GetValueA(phi) ) / ( 2 * GetValueB(phi) ) - ftg );
351 }
352 
353 #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