| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // G4TwistedBox implementation
27 //
28 // Author: 10/11/2004 - O.Link (Oliver.Link@c
29 // -------------------------------------------
30
31 #include "G4TwistedBox.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4Polyhedron.hh"
34
35 //============================================
36 //* Constructor ------------------------------
37
38 G4TwistedBox::G4TwistedBox( const G4String& pN
39 G4double pP
40 G4double pD
41 G4double pD
42 G4double pD
43 : G4VTwistedFaceted( pName, pPhiTwist,pDz,0.
44 pDy, pDx, pDx, pDy, pDx
45 {
46 }
47
48 //============================================
49 // Fake default constructor - sets only member
50 // for usage restri
51
52 G4TwistedBox::G4TwistedBox( __void__& a )
53 : G4VTwistedFaceted(a)
54 {
55 }
56
57 //============================================
58 //* Destructor -------------------------------
59
60 G4TwistedBox::~G4TwistedBox() = default;
61
62 //============================================
63 //* Copy constructor -------------------------
64
65 G4TwistedBox::G4TwistedBox(const G4TwistedBox&
66 : G4VTwistedFaceted(rhs)
67 {
68 fpPolyhedron = GetPolyhedron();
69 }
70
71 //============================================
72 //* Assignment operator ----------------------
73
74 G4TwistedBox& G4TwistedBox::operator = (const
75 {
76 // Check assignment to self
77 //
78 if (this == &rhs) { return *this; }
79
80 // Copy base class data
81 //
82 G4VTwistedFaceted::operator=(rhs);
83 fpPolyhedron = GetPolyhedron();
84
85 return *this;
86 }
87
88 //============================================
89 //* StreamInfo -------------------------------
90
91 std::ostream& G4TwistedBox::StreamInfo(std::os
92 {
93 //
94 // Stream object contents to an output strea
95 //
96 os << "-------------------------------------
97 << " *** Dump for solid - " << GetName
98 << " =================================
99 << " Solid type: G4TwistedBox\n"
100 << " Parameters: \n"
101 << " pDx = " << GetXHalfLength()/cm <
102 << " pDy = " << GetYHalfLength()/cm <
103 << " pDz = " << GetZHalfLength()/cm <
104 << " pPhiTwist = " << GetPhiTwist()/de
105 << "-------------------------------------
106
107 return os;
108 }
109
110 //============================================
111 //* GetEntityType ----------------------------
112
113 G4GeometryType G4TwistedBox::GetEntityType() c
114 {
115 return {"G4TwistedBox"};
116 }
117
118 //============================================
119 //* Clone ------------------------------------
120
121 G4VSolid* G4TwistedBox::Clone() const
122 {
123 return new G4TwistedBox(*this);
124 }
125
126 //============================================
127 //* GetCubicVolume ---------------------------
128
129 double G4TwistedBox::GetCubicVolume()
130 {
131 if (fCubicVolume == 0.)
132 {
133 fCubicVolume = 8.*GetXHalfLength()*GetYHal
134 }
135 return fCubicVolume;
136 }
137
138 //============================================
139 //* GetSurfaceArea ---------------------------
140
141 double G4TwistedBox::GetSurfaceArea()
142 {
143 if (fSurfaceArea == 0.)
144 {
145 G4double ang = GetPhiTwist();
146 G4double dx = GetXHalfLength();
147 G4double dy = GetYHalfLength();
148 G4double dz = GetZHalfLength();
149 if (ang == 0.)
150 {
151 fSurfaceArea = 8.*(dx*dy + dx*dz + dy*dz
152 }
153 else
154 {
155 G4double h = 2.*dz;
156 G4double hh = h*h;
157 G4double dxang = dx*ang;
158 G4double dyang = dy*ang;
159 fSurfaceArea = 8.*dx*dy +
160 2.*(dx*std::sqrt(hh + dxang*dxang) + h
161 2.*(dy*std::sqrt(hh + dyang*dyang) + h
162 }
163 }
164 return fSurfaceArea;
165 }
166