Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/solids/Boolean/src/G4DisplacedSolid.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /geometry/solids/Boolean/src/G4DisplacedSolid.cc (Version 11.3.0) and /geometry/solids/Boolean/src/G4DisplacedSolid.cc (Version 8.1)


  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 // Implementation of G4DisplacedSolid class fo <<  26 //
                                                   >>  27 // $Id: G4DisplacedSolid.cc,v 1.27 2006/06/29 18:43:41 gunter Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-08-01 $
                                                   >>  29 //
                                                   >>  30 // Implementation for G4DisplacedSolid class for boolean 
 27 // operations between other solids                 31 // operations between other solids
 28 //                                                 32 //
                                                   >>  33 // History:
                                                   >>  34 //
 29 // 28.10.98 V.Grichine: created                    35 // 28.10.98 V.Grichine: created
 30 // 28.02.18 E.Tcherniaev: improved contruction <<  36 // 14.11.99 V.Grichine: modifications in CalculateExtent(...) method
                                                   >>  37 // 22.11.00 V.Grichine: new set methods for matrix/vectors
                                                   >>  38 //
 31 // -------------------------------------------     39 // --------------------------------------------------------------------
 32                                                    40 
 33 #include "G4DisplacedSolid.hh"                     41 #include "G4DisplacedSolid.hh"
 34                                                    42 
 35 #include "G4VoxelLimits.hh"                        43 #include "G4VoxelLimits.hh"
 36                                                    44 
 37 #include "G4VPVParameterisation.hh"                45 #include "G4VPVParameterisation.hh"
 38                                                    46 
 39 #include "G4VGraphicsScene.hh"                     47 #include "G4VGraphicsScene.hh"
 40 #include "G4Polyhedron.hh"                         48 #include "G4Polyhedron.hh"
                                                   >>  49 #include "G4NURBS.hh"
                                                   >>  50 // #include "G4NURBSbox.hh"
 41                                                    51 
 42 //////////////////////////////////////////////     52 ////////////////////////////////////////////////////////////////
 43 //                                                 53 //
 44 // Constructor for transformation like rotatio     54 // Constructor for transformation like rotation of frame then translation 
 45 // in new frame. It is similar to 1st constrac     55 // in new frame. It is similar to 1st constractor in G4PVPlacement
 46                                                    56 
 47 G4DisplacedSolid::G4DisplacedSolid( const G4St     57 G4DisplacedSolid::G4DisplacedSolid( const G4String& pName,
 48                                           G4VS     58                                           G4VSolid* pSolid ,
 49                                           G4Ro     59                                           G4RotationMatrix* rotMatrix,
 50                                     const G4Th     60                                     const G4ThreeVector& transVector    )
 51   : G4VSolid(pName)                            <<  61   : G4VSolid(pName), fpPolyhedron(0)
 52 {                                                  62 {
 53   if (pSolid->GetEntityType() == "G4DisplacedS <<  63   fPtrSolid = pSolid ;
 54   {                                            <<  64   fPtrTransform = new G4AffineTransform(rotMatrix,transVector) ;
 55     fPtrSolid = ((G4DisplacedSolid*)pSolid)->G <<  65   fPtrTransform->Invert() ;
 56     G4AffineTransform t1 = ((G4DisplacedSolid* <<  66   fDirectTransform = new G4AffineTransform(rotMatrix,transVector) ;
 57     G4AffineTransform t2 = G4AffineTransform(r << 
 58     fDirectTransform = new G4AffineTransform(t << 
 59   }                                            << 
 60   else                                         << 
 61   {                                            << 
 62     fPtrSolid = pSolid;                        << 
 63     fDirectTransform = new G4AffineTransform(r << 
 64   }                                            << 
 65   fPtrTransform = new G4AffineTransform(fDirec << 
 66 }                                                  67 }
 67                                                    68 
 68 //////////////////////////////////////////////     69 /////////////////////////////////////////////////////////////////////////////////
 69 //                                                 70 //
 70 // Constructor                                     71 // Constructor
 71                                                    72 
 72 G4DisplacedSolid::G4DisplacedSolid( const G4St     73 G4DisplacedSolid::G4DisplacedSolid( const G4String& pName,
 73                                           G4VS     74                                           G4VSolid* pSolid ,
 74                                     const G4Tr     75                                     const G4Transform3D& transform  )
 75   : G4VSolid(pName)                            <<  76   : G4VSolid(pName), fpPolyhedron(0)
 76 {                                                  77 {
 77   if (pSolid->GetEntityType() == "G4DisplacedS <<  78   fPtrSolid = pSolid ;
 78   {                                            <<  79   fDirectTransform = new G4AffineTransform(transform.getRotation().inverse(),
 79     fPtrSolid = ((G4DisplacedSolid*)pSolid)->G <<  80                                            transform.getTranslation()) ;
 80     G4AffineTransform t1 = ((G4DisplacedSolid* <<  81 
 81     G4AffineTransform t2 = G4AffineTransform(t <<  82   fPtrTransform    = new G4AffineTransform(transform.getRotation().inverse(),
 82                                              t <<  83                                            transform.getTranslation()) ;
 83     fDirectTransform = new G4AffineTransform(t <<  84   fPtrTransform->Invert() ;
 84   }                                            << 
 85   else                                         << 
 86   {                                            << 
 87     fPtrSolid = pSolid;                        << 
 88     fDirectTransform = new G4AffineTransform(t << 
 89                                              t << 
 90   }                                            << 
 91   fPtrTransform = new G4AffineTransform(fDirec << 
 92 }                                                  85 }
 93                                                    86 
 94 //////////////////////////////////////////////     87 ///////////////////////////////////////////////////////////////////
 95 //                                                 88 //
 96 // Constructor for use with creation of Transi     89 // Constructor for use with creation of Transient object
 97 // from Persistent object                          90 // from Persistent object
 98                                                    91 
 99 G4DisplacedSolid::G4DisplacedSolid( const G4St     92 G4DisplacedSolid::G4DisplacedSolid( const G4String& pName,
100                                           G4VS     93                                           G4VSolid* pSolid ,
101                                     const G4Af     94                                     const G4AffineTransform directTransform )
102   : G4VSolid(pName)                            <<  95   : G4VSolid(pName), fpPolyhedron(0)
103 {                                                  96 {
104   if (pSolid->GetEntityType() == "G4DisplacedS <<  97   fPtrSolid = pSolid ;
105   {                                            <<  98   fDirectTransform = new G4AffineTransform( directTransform );
106     fPtrSolid = ((G4DisplacedSolid*)pSolid)->G <<  99   fPtrTransform    = new G4AffineTransform( directTransform.Inverse() ) ; 
107     G4AffineTransform t1 = ((G4DisplacedSolid* << 
108     auto  t2 = G4AffineTransform(directTransfo << 
109     fDirectTransform = new G4AffineTransform(t << 
110   }                                            << 
111   else                                         << 
112   {                                            << 
113     fPtrSolid = pSolid;                        << 
114     fDirectTransform = new G4AffineTransform(d << 
115   }                                            << 
116   fPtrTransform = new G4AffineTransform(fDirec << 
117 }                                                 100 }
118                                                   101 
119 //////////////////////////////////////////////    102 ///////////////////////////////////////////////////////////////////
120 //                                                103 //
121 // Fake default constructor - sets only member    104 // Fake default constructor - sets only member data and allocates memory
122 //                            for usage restri    105 //                            for usage restricted to object persistency.
123                                                   106 
124 G4DisplacedSolid::G4DisplacedSolid( __void__&     107 G4DisplacedSolid::G4DisplacedSolid( __void__& a )
125   : G4VSolid(a)                                << 108   : G4VSolid(a), fPtrSolid(0), fPtrTransform(0),
                                                   >> 109     fDirectTransform(0), fpPolyhedron(0)
126 {                                                 110 {
127 }                                                 111 }
128                                                   112 
129 //////////////////////////////////////////////    113 ///////////////////////////////////////////////////////////////////
130 //                                                114 //
131 // Destructor                                     115 // Destructor
132                                                   116 
133 G4DisplacedSolid::~G4DisplacedSolid()             117 G4DisplacedSolid::~G4DisplacedSolid() 
134 {                                                 118 {
135   CleanTransformations();                         119   CleanTransformations();
136   delete fpPolyhedron; fpPolyhedron = nullptr; << 120   delete fpPolyhedron;
137 }                                                 121 }
138                                                   122 
139 ////////////////////////////////////////////// << 123 G4GeometryType G4DisplacedSolid::GetEntityType() const 
140 //                                             << 
141 // Copy constructor                            << 
142                                                << 
143 G4DisplacedSolid::G4DisplacedSolid(const G4Dis << 
144   : G4VSolid (rhs), fPtrSolid(rhs.fPtrSolid)   << 
145 {                                                 124 {
146   fPtrTransform = new G4AffineTransform(*(rhs. << 125   return G4String("G4DisplacedSolid");
147   fDirectTransform = new G4AffineTransform(*(r << 
148 }                                                 126 }
149                                                   127 
150 ////////////////////////////////////////////// << 
151 //                                             << 
152 // Assignment operator                         << 
153                                                << 
154 G4DisplacedSolid& G4DisplacedSolid::operator = << 
155 {                                              << 
156   // Check assignment to self                  << 
157   //                                           << 
158   if (this == &rhs)  { return *this; }         << 
159                                                << 
160   // Copy base class data                      << 
161   //                                           << 
162   G4VSolid::operator=(rhs);                    << 
163                                                << 
164   // Copy data                                 << 
165   //                                           << 
166   fPtrSolid = rhs.fPtrSolid;                   << 
167   delete fPtrTransform; delete fDirectTransfor << 
168   fPtrTransform = new G4AffineTransform(*(rhs. << 
169   fDirectTransform = new G4AffineTransform(*(r << 
170   fRebuildPolyhedron = false;                  << 
171   delete fpPolyhedron; fpPolyhedron = nullptr; << 
172                                                << 
173   return *this;                                << 
174 }                                              << 
175                                                << 
176 void G4DisplacedSolid::CleanTransformations()     128 void G4DisplacedSolid::CleanTransformations()
177 {                                                 129 {
178   if(fPtrTransform != nullptr)                 << 130   if(fPtrTransform)
179   {                                               131   {
180     delete fPtrTransform; fPtrTransform = null << 132     delete fPtrTransform;  fPtrTransform=0;
181     delete fDirectTransform; fDirectTransform  << 133     delete fDirectTransform;  fDirectTransform=0;
182   }                                               134   }
183 }                                                 135 }
184                                                   136 
185 const G4DisplacedSolid* G4DisplacedSolid::GetD    137 const G4DisplacedSolid* G4DisplacedSolid::GetDisplacedSolidPtr() const   
186 {                                                 138 {
187   return this;                                    139   return this;
188 }                                                 140 }
189                                                   141 
190 G4DisplacedSolid* G4DisplacedSolid::GetDisplac    142 G4DisplacedSolid* G4DisplacedSolid::GetDisplacedSolidPtr() 
191 {                                                 143 {
192   return this;                                    144   return this;
193 }                                                 145 }
194                                                   146 
195 G4VSolid* G4DisplacedSolid::GetConstituentMove    147 G4VSolid* G4DisplacedSolid::GetConstituentMovedSolid() const
196 {                                                 148 { 
197   return fPtrSolid;                               149   return fPtrSolid; 
198 }                                                 150 } 
199                                                   151 
200 //////////////////////////////////////////////    152 /////////////////////////////////////////////////////////////////////////////
201                                                   153 
202 G4AffineTransform  G4DisplacedSolid::GetTransf    154 G4AffineTransform  G4DisplacedSolid::GetTransform() const
203 {                                                 155 {
204   G4AffineTransform aTransform = *fPtrTransfor    156   G4AffineTransform aTransform = *fPtrTransform;
205   return aTransform;                              157   return aTransform;
206 }                                                 158 }
207                                                   159 
208 void G4DisplacedSolid::SetTransform(G4AffineTr    160 void G4DisplacedSolid::SetTransform(G4AffineTransform& transform) 
209 {                                                 161 {
210   fPtrTransform = &transform ;                    162   fPtrTransform = &transform ;
211   fRebuildPolyhedron = true;                   << 163   fpPolyhedron = 0;
212 }                                                 164 }
213                                                   165 
214 //////////////////////////////////////////////    166 //////////////////////////////////////////////////////////////////////////////
215                                                   167 
216 G4AffineTransform  G4DisplacedSolid::GetDirect    168 G4AffineTransform  G4DisplacedSolid::GetDirectTransform() const
217 {                                                 169 {
218   G4AffineTransform aTransform= *fDirectTransf    170   G4AffineTransform aTransform= *fDirectTransform;
219   return aTransform;                              171   return aTransform;
220 }                                                 172 }
221                                                   173 
222 void G4DisplacedSolid::SetDirectTransform(G4Af    174 void G4DisplacedSolid::SetDirectTransform(G4AffineTransform& transform) 
223 {                                                 175 {
224   fDirectTransform = &transform ;                 176   fDirectTransform = &transform ;
225   fRebuildPolyhedron = true;                   << 177   fpPolyhedron = 0;
226 }                                                 178 }
227                                                   179 
228 //////////////////////////////////////////////    180 /////////////////////////////////////////////////////////////////////////////
229                                                   181 
230 G4RotationMatrix G4DisplacedSolid::GetFrameRot    182 G4RotationMatrix G4DisplacedSolid::GetFrameRotation() const
231 {                                                 183 {
232   G4RotationMatrix InvRotation = fDirectTransf << 184   G4RotationMatrix InvRotation= fDirectTransform->NetRotation();
233   return InvRotation;                             185   return InvRotation;
234 }                                                 186 }
235                                                   187 
236 void G4DisplacedSolid::SetFrameRotation(const     188 void G4DisplacedSolid::SetFrameRotation(const G4RotationMatrix& matrix)
237 {                                                 189 {
238   fDirectTransform->SetNetRotation(matrix);       190   fDirectTransform->SetNetRotation(matrix);
239   fRebuildPolyhedron = true;                   << 191   fpPolyhedron = 0;
240 }                                                 192 }
241                                                   193 
242 //////////////////////////////////////////////    194 /////////////////////////////////////////////////////////////////////////////
243                                                   195 
244 G4ThreeVector  G4DisplacedSolid::GetFrameTrans    196 G4ThreeVector  G4DisplacedSolid::GetFrameTranslation() const
245 {                                                 197 {
246   return fPtrTransform->NetTranslation();         198   return fPtrTransform->NetTranslation();
247 }                                                 199 }
248                                                   200 
249 void G4DisplacedSolid::SetFrameTranslation(con    201 void G4DisplacedSolid::SetFrameTranslation(const G4ThreeVector& vector)
250 {                                                 202 {
251   fPtrTransform->SetNetTranslation(vector);       203   fPtrTransform->SetNetTranslation(vector);
252   fRebuildPolyhedron = true;                   << 204   fpPolyhedron = 0;
253 }                                                 205 }
254                                                   206 
255 //////////////////////////////////////////////    207 ///////////////////////////////////////////////////////////////
256                                                   208 
257 G4RotationMatrix G4DisplacedSolid::GetObjectRo    209 G4RotationMatrix G4DisplacedSolid::GetObjectRotation() const
258 {                                                 210 {
259   G4RotationMatrix Rotation = fPtrTransform->N << 211   G4RotationMatrix Rotation= fPtrTransform->NetRotation();
260   return Rotation;                                212   return Rotation;
261 }                                                 213 }
262                                                   214 
263 void G4DisplacedSolid::SetObjectRotation(const    215 void G4DisplacedSolid::SetObjectRotation(const G4RotationMatrix& matrix)
264 {                                                 216 {
265   fPtrTransform->SetNetRotation(matrix);          217   fPtrTransform->SetNetRotation(matrix);
266   fRebuildPolyhedron = true;                   << 218   fpPolyhedron = 0;
267 }                                                 219 }
268                                                   220 
269 //////////////////////////////////////////////    221 ///////////////////////////////////////////////////////////////////////
270                                                   222 
271 G4ThreeVector  G4DisplacedSolid::GetObjectTran    223 G4ThreeVector  G4DisplacedSolid::GetObjectTranslation() const
272 {                                                 224 {
273   return fDirectTransform->NetTranslation();      225   return fDirectTransform->NetTranslation();
274 }                                                 226 }
275                                                   227 
276 void G4DisplacedSolid::SetObjectTranslation(co    228 void G4DisplacedSolid::SetObjectTranslation(const G4ThreeVector& vector)
277 {                                                 229 {
278   fDirectTransform->SetNetTranslation(vector);    230   fDirectTransform->SetNetTranslation(vector);
279   fRebuildPolyhedron = true;                   << 231   fpPolyhedron = 0;
280 }                                                 232 }
281                                                   233 
282 ////////////////////////////////////////////// << 234 ///////////////////////////////////////////////////////////////
283 //                                                235 //
284 // Get bounding box                            << 
285                                                << 
286 void G4DisplacedSolid::BoundingLimits(G4ThreeV << 
287                                       G4ThreeV << 
288 {                                              << 
289   if (!fDirectTransform->IsRotated())          << 
290   {                                            << 
291     // Special case of pure translation        << 
292     //                                         << 
293     fPtrSolid->BoundingLimits(pMin,pMax);      << 
294     G4ThreeVector offset = fDirectTransform->N << 
295     pMin += offset;                            << 
296     pMax += offset;                            << 
297   }                                            << 
298   else                                         << 
299   {                                            << 
300     // General case, use CalculateExtent() to  << 
301     //                                         << 
302     G4VoxelLimits unLimit;                     << 
303     G4double xmin,xmax,ymin,ymax,zmin,zmax;    << 
304     fPtrSolid->CalculateExtent(kXAxis,unLimit, << 
305     fPtrSolid->CalculateExtent(kYAxis,unLimit, << 
306     fPtrSolid->CalculateExtent(kZAxis,unLimit, << 
307     pMin.set(xmin,ymin,zmin);                  << 
308     pMax.set(xmax,ymax,zmax);                  << 
309   }                                            << 
310                                                << 
311   // Check correctness of the bounding box     << 
312   //                                           << 
313   if (pMin.x() >= pMax.x() || pMin.y() >= pMax << 
314   {                                            << 
315     std::ostringstream message;                << 
316     message << "Bad bounding box (min >= max)  << 
317             << GetName() << " !"               << 
318             << "\npMin = " << pMin             << 
319             << "\npMax = " << pMax;            << 
320     G4Exception("G4DisplacedSolid::BoundingLim << 
321                JustWarning, message);          << 
322     DumpInfo();                                << 
323   }                                            << 
324 }                                              << 
325                                                << 
326 ////////////////////////////////////////////// << 
327 //                                                236 //
328 // Calculate extent under transform and specif << 
329                                                   237      
330 G4bool                                            238 G4bool 
331 G4DisplacedSolid::CalculateExtent( const EAxis    239 G4DisplacedSolid::CalculateExtent( const EAxis pAxis,
332                                    const G4Vox    240                                    const G4VoxelLimits& pVoxelLimit,
333                                    const G4Aff    241                                    const G4AffineTransform& pTransform,
334                                          G4dou    242                                          G4double& pMin, 
335                                          G4dou    243                                          G4double& pMax           ) const 
336 {                                                 244 {
337   G4AffineTransform sumTransform ;                245   G4AffineTransform sumTransform ;
338   sumTransform.Product(*fDirectTransform,pTran    246   sumTransform.Product(*fDirectTransform,pTransform) ;
339   return fPtrSolid->CalculateExtent(pAxis,pVox    247   return fPtrSolid->CalculateExtent(pAxis,pVoxelLimit,sumTransform,pMin,pMax) ;
340 }                                                 248 }
341                                                   249  
342 //////////////////////////////////////////////    250 /////////////////////////////////////////////////////
343 //                                                251 //
344 // SurfaceNormal                               << 252 // 
345                                                   253 
346 EInside G4DisplacedSolid::Inside(const G4Three    254 EInside G4DisplacedSolid::Inside(const G4ThreeVector& p) const
347 {                                                 255 {
348   G4ThreeVector newPoint = fPtrTransform->Tran    256   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
349   return fPtrSolid->Inside(newPoint) ;            257   return fPtrSolid->Inside(newPoint) ; 
350 }                                                 258 }
351                                                   259 
352 //////////////////////////////////////////////    260 //////////////////////////////////////////////////////////////
353 //                                                261 //
354 //                                                262 //
355                                                   263 
356 G4ThreeVector                                     264 G4ThreeVector 
357 G4DisplacedSolid::SurfaceNormal( const G4Three    265 G4DisplacedSolid::SurfaceNormal( const G4ThreeVector& p ) const 
358 {                                                 266 {
359   G4ThreeVector newPoint = fPtrTransform->Tran    267   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
360   G4ThreeVector normal = fPtrSolid->SurfaceNor    268   G4ThreeVector normal = fPtrSolid->SurfaceNormal(newPoint) ; 
361   return fDirectTransform->TransformAxis(norma    269   return fDirectTransform->TransformAxis(normal) ;
362 }                                                 270 }
363                                                   271 
364 //////////////////////////////////////////////    272 /////////////////////////////////////////////////////////////
365 //                                                273 //
366 // The same algorithm as in DistanceToIn(p)       274 // The same algorithm as in DistanceToIn(p)
367                                                   275 
368 G4double                                          276 G4double 
369 G4DisplacedSolid::DistanceToIn( const G4ThreeV    277 G4DisplacedSolid::DistanceToIn( const G4ThreeVector& p,
370                                 const G4ThreeV    278                                 const G4ThreeVector& v  ) const 
371 {                                                 279 {    
372   G4ThreeVector newPoint = fPtrTransform->Tran    280   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
373   G4ThreeVector newDirection = fPtrTransform->    281   G4ThreeVector newDirection = fPtrTransform->TransformAxis(v) ;
374   return fPtrSolid->DistanceToIn(newPoint,newD    282   return fPtrSolid->DistanceToIn(newPoint,newDirection) ;   
375 }                                                 283 }
376                                                   284 
377 //////////////////////////////////////////////    285 ////////////////////////////////////////////////////////
378 //                                                286 //
379 // Approximate nearest distance from the point    287 // Approximate nearest distance from the point p to the intersection of
380 // two solids                                     288 // two solids
381                                                   289 
382 G4double                                          290 G4double 
383 G4DisplacedSolid::DistanceToIn( const G4ThreeV    291 G4DisplacedSolid::DistanceToIn( const G4ThreeVector& p ) const 
384 {                                                 292 {
385   G4ThreeVector newPoint = fPtrTransform->Tran    293   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
386   return fPtrSolid->DistanceToIn(newPoint) ;      294   return fPtrSolid->DistanceToIn(newPoint) ;   
387 }                                                 295 }
388                                                   296 
389 //////////////////////////////////////////////    297 //////////////////////////////////////////////////////////
390 //                                                298 //
391 // The same algorithm as DistanceToOut(p)         299 // The same algorithm as DistanceToOut(p)
392                                                   300 
393 G4double                                          301 G4double 
394 G4DisplacedSolid::DistanceToOut( const G4Three    302 G4DisplacedSolid::DistanceToOut( const G4ThreeVector& p,
395                                  const G4Three    303                                  const G4ThreeVector& v,
396                                  const G4bool     304                                  const G4bool calcNorm,
397                                        G4bool     305                                        G4bool *validNorm,
398                                        G4Three    306                                        G4ThreeVector *n   ) const 
399 {                                                 307 {
400   G4ThreeVector solNorm ;                         308   G4ThreeVector solNorm ; 
401   G4ThreeVector newPoint = fPtrTransform->Tran    309   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
402   G4ThreeVector newDirection = fPtrTransform->    310   G4ThreeVector newDirection = fPtrTransform->TransformAxis(v) ;
403   G4double dist = fPtrSolid->DistanceToOut(new    311   G4double dist = fPtrSolid->DistanceToOut(newPoint,newDirection,
404                                            cal    312                                            calcNorm,validNorm,&solNorm) ;
405   if(calcNorm)                                    313   if(calcNorm)
406   {                                               314   { 
407     *n = fDirectTransform->TransformAxis(solNo    315     *n = fDirectTransform->TransformAxis(solNorm) ;
408   }                                               316   }
409   return dist ;                                   317   return dist ;  
410 }                                                 318 }
411                                                   319 
412 //////////////////////////////////////////////    320 //////////////////////////////////////////////////////////////
413 //                                                321 //
414 // Inverted algorithm of DistanceToIn(p)          322 // Inverted algorithm of DistanceToIn(p)
415                                                   323 
416 G4double                                          324 G4double 
417 G4DisplacedSolid::DistanceToOut( const G4Three    325 G4DisplacedSolid::DistanceToOut( const G4ThreeVector& p ) const 
418 {                                                 326 {
419   G4ThreeVector newPoint = fPtrTransform->Tran    327   G4ThreeVector newPoint = fPtrTransform->TransformPoint(p) ;
420   return fPtrSolid->DistanceToOut(newPoint) ;     328   return fPtrSolid->DistanceToOut(newPoint) ;   
421 }                                                 329 }
422                                                   330 
423 //////////////////////////////////////////////    331 //////////////////////////////////////////////////////////////
424 //                                                332 //
425 // ComputeDimensions                           << 333 //
426                                                   334 
427 void                                              335 void 
428 G4DisplacedSolid::ComputeDimensions(       G4V    336 G4DisplacedSolid::ComputeDimensions(       G4VPVParameterisation*,
429                                      const G4i    337                                      const G4int,
430                                      const G4V    338                                      const G4VPhysicalVolume* ) 
431 {                                                 339 {
432   DumpInfo();                                     340   DumpInfo();
433   G4Exception("G4DisplacedSolid::ComputeDimens    341   G4Exception("G4DisplacedSolid::ComputeDimensions()",
434               "GeomSolids0001", FatalException << 342                 "NotApplicable", FatalException,
435               "Method not applicable in this c << 343                 "Method not applicable in this context!");
436 }                                              << 
437                                                << 
438 ////////////////////////////////////////////// << 
439 //                                             << 
440 // Return volume                               << 
441                                                << 
442 G4double G4DisplacedSolid::GetCubicVolume()    << 
443 {                                              << 
444   return fPtrSolid->GetCubicVolume();          << 
445 }                                              << 
446                                                << 
447 ////////////////////////////////////////////// << 
448 //                                             << 
449 // Return surface area                         << 
450                                                << 
451 G4double G4DisplacedSolid::GetSurfaceArea()    << 
452 {                                              << 
453   return fPtrSolid->GetSurfaceArea();          << 
454 }                                                 344 }
455                                                   345 
456 //////////////////////////////////////////////    346 //////////////////////////////////////////////////////////////////////////
457 //                                                347 //
458 // Returns a point (G4ThreeVector) randomly an    348 // Returns a point (G4ThreeVector) randomly and uniformly selected
459 // on the solid surface                           349 // on the solid surface
460 //                                                350 //
461                                                   351 
462 G4ThreeVector G4DisplacedSolid::GetPointOnSurf    352 G4ThreeVector G4DisplacedSolid::GetPointOnSurface() const
463 {                                                 353 {
464   G4ThreeVector p = fPtrSolid->GetPointOnSurfa << 354   G4ThreeVector p =  fPtrSolid->GetPointOnSurface();
465   return fDirectTransform->TransformPoint(p);     355   return fDirectTransform->TransformPoint(p);
466 }                                                 356 }
467                                                   357 
468 //////////////////////////////////////////////    358 //////////////////////////////////////////////////////////////////////////
469 //                                                359 //
470 // Return the number of constituents used for  << 
471                                                << 
472 G4int G4DisplacedSolid::GetNumOfConstituents() << 
473 {                                              << 
474   return fPtrSolid->GetNumOfConstituents();    << 
475 }                                              << 
476                                                << 
477 ////////////////////////////////////////////// << 
478 //                                             << 
479 // Return true if the solid has only planar fa << 
480                                                << 
481 G4bool G4DisplacedSolid::IsFaceted() const     << 
482 {                                              << 
483   return fPtrSolid->IsFaceted();               << 
484 }                                              << 
485                                                << 
486 ////////////////////////////////////////////// << 
487 //                                             << 
488 // Return object type name                     << 
489                                                << 
490 G4GeometryType G4DisplacedSolid::GetEntityType << 
491 {                                              << 
492   return {"G4DisplacedSolid"};                 << 
493 }                                              << 
494                                                << 
495 ////////////////////////////////////////////// << 
496 //                                             << 
497 // Make a clone of the object                  << 
498 //                                             << 
499 G4VSolid* G4DisplacedSolid::Clone() const      << 
500 {                                              << 
501   return new G4DisplacedSolid(*this);          << 
502 }                                              << 
503                                                << 
504 ////////////////////////////////////////////// << 
505 //                                             << 
506 // Stream object contents to an output stream     360 // Stream object contents to an output stream
507                                                   361 
508 std::ostream& G4DisplacedSolid::StreamInfo(std    362 std::ostream& G4DisplacedSolid::StreamInfo(std::ostream& os) const
509 {                                                 363 {
510   os << "-------------------------------------    364   os << "-----------------------------------------------------------\n"
511      << "    *** Dump for Displaced solid - "     365      << "    *** Dump for Displaced solid - " << GetName() << " ***\n"
512      << "    =================================    366      << "    ===================================================\n"
513      << " Solid type: " << GetEntityType() <<     367      << " Solid type: " << GetEntityType() << "\n"
514      << " Parameters of constituent solid: \n"    368      << " Parameters of constituent solid: \n"
515      << "=====================================    369      << "===========================================================\n";
516   fPtrSolid->StreamInfo(os);                      370   fPtrSolid->StreamInfo(os);
517   os << "=====================================    371   os << "===========================================================\n"
518      << " Transformations: \n"                    372      << " Transformations: \n"
519      << "    Direct transformation - translati    373      << "    Direct transformation - translation : \n"
520      << "           " << fDirectTransform->Net    374      << "           " << fDirectTransform->NetTranslation() << "\n"
521      << "                          - rotation     375      << "                          - rotation    : \n"
522      << "           ";                            376      << "           ";
523   fDirectTransform->NetRotation().print(os);      377   fDirectTransform->NetRotation().print(os);
524   os << "\n"                                      378   os << "\n"
525      << "=====================================    379      << "===========================================================\n";
526                                                   380 
527   return os;                                      381   return os;
528 }                                                 382 }
529                                                   383 
530 //////////////////////////////////////////////    384 //////////////////////////////////////////////////////////////////////////
531 //                                                385 //
532 // DescribeYourselfTo                          << 386 //                    
533                                                   387 
534 void                                              388 void 
535 G4DisplacedSolid::DescribeYourselfTo ( G4VGrap    389 G4DisplacedSolid::DescribeYourselfTo ( G4VGraphicsScene& scene ) const 
536 {                                                 390 {
537   scene.AddSolid (*this);                         391   scene.AddSolid (*this);
538 }                                                 392 }
539                                                   393 
540 //////////////////////////////////////////////    394 //////////////////////////////////////////////////////////////////////////
541 //                                                395 //
542 // CreatePolyhedron                            << 396 //
543                                                   397 
544 G4Polyhedron*                                     398 G4Polyhedron* 
545 G4DisplacedSolid::CreatePolyhedron () const       399 G4DisplacedSolid::CreatePolyhedron () const 
546 {                                                 400 {
547   G4Polyhedron* polyhedron = fPtrSolid->Create    401   G4Polyhedron* polyhedron = fPtrSolid->CreatePolyhedron();
548   if (polyhedron != nullptr)                   << 402   polyhedron
549   {                                            << 
550     polyhedron                                 << 
551     ->Transform(G4Transform3D(GetObjectRotatio    403     ->Transform(G4Transform3D(GetObjectRotation(),GetObjectTranslation()));
552   }                                            << 
553   else                                         << 
554   {                                            << 
555     DumpInfo();                                << 
556     G4Exception("G4DisplacedSolid::CreatePolyh << 
557                 "GeomSolids2002", JustWarning, << 
558                 "No G4Polyhedron for displaced << 
559   }                                            << 
560   return polyhedron;                              404   return polyhedron;
561 }                                                 405 }
562                                                   406 
563 //////////////////////////////////////////////    407 //////////////////////////////////////////////////////////////////////////
564 //                                                408 //
565 // GetPolyhedron                               << 409 //
                                                   >> 410 
                                                   >> 411 G4NURBS*      
                                                   >> 412 G4DisplacedSolid::CreateNURBS () const 
                                                   >> 413 {
                                                   >> 414   // Take into account local transformation - see CreatePolyhedron.
                                                   >> 415   // return fPtrSolid->CreateNURBS() ;
                                                   >> 416   return 0;
                                                   >> 417 }
                                                   >> 418 
                                                   >> 419 //////////////////////////////////////////////////////////////////////////
                                                   >> 420 //
                                                   >> 421 //
566                                                   422 
567 G4Polyhedron* G4DisplacedSolid::GetPolyhedron     423 G4Polyhedron* G4DisplacedSolid::GetPolyhedron () const
568 {                                                 424 {
569   if (fpPolyhedron == nullptr ||               << 425   if (!fpPolyhedron ||
570       fRebuildPolyhedron ||                    << 
571       fpPolyhedron->GetNumberOfRotationStepsAt    426       fpPolyhedron->GetNumberOfRotationStepsAtTimeOfCreation() !=
572       fpPolyhedron->GetNumberOfRotationSteps()    427       fpPolyhedron->GetNumberOfRotationSteps())
573     {                                             428     {
                                                   >> 429       delete fpPolyhedron;
574       fpPolyhedron = CreatePolyhedron();          430       fpPolyhedron = CreatePolyhedron();
575       fRebuildPolyhedron = false;              << 
576     }                                             431     }
577   return fpPolyhedron;                            432   return fpPolyhedron;
578 }                                                 433 }
579                                                   434