| Geant4 Cross Reference |
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 // G4TwistedBox implementation
27 //
28 // Author: 10/11/2004 - O.Link (Oliver.Link@cern.ch)
29 // --------------------------------------------------------------------
30
31 #include "G4TwistedBox.hh"
32 #include "G4SystemOfUnits.hh"
33 #include "G4Polyhedron.hh"
34
35 //=====================================================================
36 //* Constructor -------------------------------------------------------
37
38 G4TwistedBox::G4TwistedBox( const G4String& pName,
39 G4double pPhiTwist,
40 G4double pDx,
41 G4double pDy,
42 G4double pDz )
43 : G4VTwistedFaceted( pName, pPhiTwist,pDz,0.,0.,
44 pDy, pDx, pDx, pDy, pDx, pDx,0. )
45 {
46 }
47
48 //=====================================================================
49 // Fake default constructor - sets only member data and allocates memory
50 // for usage restricted to object persistency.
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& rhs)
66 : G4VTwistedFaceted(rhs)
67 {
68 fpPolyhedron = GetPolyhedron();
69 }
70
71 //=====================================================================
72 //* Assignment operator -----------------------------------------------
73
74 G4TwistedBox& G4TwistedBox::operator = (const G4TwistedBox& rhs)
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::ostream& os) const
92 {
93 //
94 // Stream object contents to an output stream
95 //
96 os << "-----------------------------------------------------------\n"
97 << " *** Dump for solid - " << GetName() << " ***\n"
98 << " ===================================================\n"
99 << " Solid type: G4TwistedBox\n"
100 << " Parameters: \n"
101 << " pDx = " << GetXHalfLength()/cm << " cm" << G4endl
102 << " pDy = " << GetYHalfLength()/cm << " cm" << G4endl
103 << " pDz = " << GetZHalfLength()/cm << " cm" << G4endl
104 << " pPhiTwist = " << GetPhiTwist()/degree << " deg" << G4endl
105 << "-----------------------------------------------------------\n";
106
107 return os;
108 }
109
110 //=====================================================================
111 //* GetEntityType -----------------------------------------------------
112
113 G4GeometryType G4TwistedBox::GetEntityType() const
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()*GetYHalfLength()*GetZHalfLength();
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) + hh*std::asinh(dxang/h)/ang) +
161 2.*(dy*std::sqrt(hh + dyang*dyang) + hh*std::asinh(dyang/h)/ang);
162 }
163 }
164 return fSurfaceArea;
165 }
166