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