Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/solids/CSG/include/G4Cons.icc

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/CSG/include/G4Cons.icc (Version 11.3.0) and /geometry/solids/CSG/include/G4Cons.icc (Version 10.5)


  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 //
                                                   >>  27 //
                                                   >>  28 // --------------------------------------------------------------------
                                                   >>  29 // GEANT 4 inline definitions file
                                                   >>  30 //
                                                   >>  31 // G4Cons.icc
                                                   >>  32 //
 26 // Implementation of inline methods of G4Cons      33 // Implementation of inline methods of G4Cons
 27 // -------------------------------------------     34 // --------------------------------------------------------------------
 28                                                    35 
 29 inline                                             36 inline
 30 G4double G4Cons::GetInnerRadiusMinusZ() const      37 G4double G4Cons::GetInnerRadiusMinusZ() const
 31 {                                                  38 {
 32   return fRmin1 ;                                  39   return fRmin1 ;
 33 }                                                  40 }
 34                                                    41 
 35 inline                                             42 inline
 36 G4double G4Cons::GetOuterRadiusMinusZ() const      43 G4double G4Cons::GetOuterRadiusMinusZ() const
 37 {                                                  44 {
 38   return fRmax1 ;                                  45   return fRmax1 ;
 39 }                                                  46 }
 40                                                    47 
 41 inline                                             48 inline
 42 G4double G4Cons::GetInnerRadiusPlusZ() const       49 G4double G4Cons::GetInnerRadiusPlusZ() const
 43 {                                                  50 {
 44   return fRmin2 ;                                  51   return fRmin2 ;
 45 }                                                  52 }
 46                                                    53 
 47 inline                                             54 inline
 48 G4double G4Cons::GetOuterRadiusPlusZ() const       55 G4double G4Cons::GetOuterRadiusPlusZ() const
 49 {                                                  56 {
 50   return fRmax2 ;                                  57   return fRmax2 ;
 51 }                                                  58 }
 52                                                    59 
 53 inline                                             60 inline
 54 G4double G4Cons::GetZHalfLength() const            61 G4double G4Cons::GetZHalfLength() const
 55 {                                                  62 {
 56   return fDz ;                                     63   return fDz ;
 57 }                                                  64 }
 58                                                    65 
 59 inline                                         <<  66 inline  
 60 G4double G4Cons::GetStartPhiAngle() const          67 G4double G4Cons::GetStartPhiAngle() const
 61 {                                                  68 {
 62   return fSPhi ;                                   69   return fSPhi ;
 63 }                                                  70 }
 64                                                    71 
 65 inline                                             72 inline
 66 G4double G4Cons::GetDeltaPhiAngle() const          73 G4double G4Cons::GetDeltaPhiAngle() const
 67 {                                                  74 {
 68   return fDPhi;                                    75   return fDPhi;
 69 }                                                  76 }
 70                                                    77 
 71 inline                                             78 inline
 72 G4double G4Cons::GetSinStartPhi() const            79 G4double G4Cons::GetSinStartPhi() const
 73 {                                                  80 {
 74   return sinSPhi;                                  81   return sinSPhi;
 75 }                                                  82 }
 76                                                    83 
 77 inline                                             84 inline
 78 G4double G4Cons::GetCosStartPhi() const            85 G4double G4Cons::GetCosStartPhi() const
 79 {                                                  86 {
 80   return cosSPhi;                                  87   return cosSPhi;
 81 }                                                  88 }
 82                                                    89 
 83 inline                                             90 inline
 84 G4double G4Cons::GetSinEndPhi() const              91 G4double G4Cons::GetSinEndPhi() const
 85 {                                                  92 {
 86   return sinEPhi;                                  93   return sinEPhi;
 87 }                                                  94 }
 88                                                    95 
 89 inline                                             96 inline
 90 G4double G4Cons::GetCosEndPhi() const              97 G4double G4Cons::GetCosEndPhi() const
 91 {                                                  98 {
 92   return cosEPhi;                                  99   return cosEPhi;
 93 }                                                 100 }
 94                                                   101 
 95 inline                                            102 inline
 96 void G4Cons::Initialize()                         103 void G4Cons::Initialize()
 97 {                                                 104 {
 98   fCubicVolume = 0.;                              105   fCubicVolume = 0.;
 99   fSurfaceArea = 0.;                              106   fSurfaceArea = 0.;
100   fRebuildPolyhedron = true;                      107   fRebuildPolyhedron = true;
101 }                                                 108 }
102                                                   109 
103 inline                                         << 110 inline 
104 void G4Cons::InitializeTrigonometry()             111 void G4Cons::InitializeTrigonometry()
105 {                                                 112 {
106   G4double hDPhi = 0.5*fDPhi;                     113   G4double hDPhi = 0.5*fDPhi;                       // half delta phi
107   G4double cPhi  = fSPhi + hDPhi;              << 114   G4double cPhi  = fSPhi + hDPhi; 
108   G4double ePhi  = fSPhi + fDPhi;                 115   G4double ePhi  = fSPhi + fDPhi;
109                                                   116 
110   sinCPhi    = std::sin(cPhi);                    117   sinCPhi    = std::sin(cPhi);
111   cosCPhi    = std::cos(cPhi);                    118   cosCPhi    = std::cos(cPhi);
112   cosHDPhi   = std::cos(hDPhi);                << 
113   cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolera    119   cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolerance); // inner/outer tol half dphi
114   cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolera    120   cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolerance);
115   sinSPhi = std::sin(fSPhi);                      121   sinSPhi = std::sin(fSPhi);
116   cosSPhi = std::cos(fSPhi);                      122   cosSPhi = std::cos(fSPhi);
117   sinEPhi = std::sin(ePhi);                       123   sinEPhi = std::sin(ePhi);
118   cosEPhi = std::cos(ePhi);                       124   cosEPhi = std::cos(ePhi);
119 }                                                 125 }
120                                                   126 
121 inline void G4Cons::CheckSPhiAngle(G4double sP    127 inline void G4Cons::CheckSPhiAngle(G4double sPhi)
122 {                                                 128 {
123   // Ensure fSphi in 0-2PI or -2PI-0 range if     129   // Ensure fSphi in 0-2PI or -2PI-0 range if shape crosses 0
124                                                   130 
125   if ( sPhi < 0 )                                 131   if ( sPhi < 0 )
126   {                                               132   {
127     fSPhi = CLHEP::twopi - std::fmod(std::fabs    133     fSPhi = CLHEP::twopi - std::fmod(std::fabs(sPhi),CLHEP::twopi);
128   }                                               134   }
129   else                                            135   else
130   {                                               136   {
131     fSPhi = std::fmod(sPhi,CLHEP::twopi) ;        137     fSPhi = std::fmod(sPhi,CLHEP::twopi) ;
132   }                                               138   }
133   if ( fSPhi+fDPhi > CLHEP::twopi )               139   if ( fSPhi+fDPhi > CLHEP::twopi )
134   {                                               140   {
135     fSPhi -= CLHEP::twopi ;                       141     fSPhi -= CLHEP::twopi ;
136   }                                               142   }
137 }                                                 143 }
138                                                   144 
139 inline void G4Cons::CheckDPhiAngle(G4double dP    145 inline void G4Cons::CheckDPhiAngle(G4double dPhi)
140 {                                                 146 {
141   fPhiFullCone = true;                            147   fPhiFullCone = true;
142   if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5     148   if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5 )
143   {                                               149   {
144     fDPhi=CLHEP::twopi;                           150     fDPhi=CLHEP::twopi;
145     fSPhi=0;                                      151     fSPhi=0;
146   }                                               152   }
147   else                                            153   else
148   {                                               154   {
149     fPhiFullCone = false;                         155     fPhiFullCone = false;
150     if ( dPhi > 0 )                               156     if ( dPhi > 0 )
151     {                                             157     {
152       fDPhi = dPhi;                               158       fDPhi = dPhi;
153     }                                             159     }
154     else                                          160     else
155     {                                             161     {
156       std::ostringstream message;                 162       std::ostringstream message;
157       message << "Invalid dphi." << G4endl        163       message << "Invalid dphi." << G4endl
158               << "Negative or zero delta-Phi (    164               << "Negative or zero delta-Phi (" << dPhi << ") in solid: "
159               << GetName();                       165               << GetName();
160       G4Exception("G4Cons::CheckDPhiAngle()",     166       G4Exception("G4Cons::CheckDPhiAngle()", "GeomSolids0002",
161                   FatalException, message);       167                   FatalException, message);
162     }                                             168     }
163   }                                               169   }
164 }                                                 170 }
165                                                   171 
166 inline void G4Cons::CheckPhiAngles(G4double sP    172 inline void G4Cons::CheckPhiAngles(G4double sPhi, G4double dPhi)
167 {                                                 173 {
168   CheckDPhiAngle(dPhi);                           174   CheckDPhiAngle(dPhi);
169   if ( (fDPhi<CLHEP::twopi) && ((sPhi) != 0.0) << 175   if ( (fDPhi<CLHEP::twopi) && (sPhi) ) { CheckSPhiAngle(sPhi); }
170   InitializeTrigonometry();                       176   InitializeTrigonometry();
171 }                                                 177 }
172                                                   178 
173 inline                                            179 inline
174 void G4Cons::SetInnerRadiusMinusZ( G4double Rm    180 void G4Cons::SetInnerRadiusMinusZ( G4double Rmin1 )
175 {                                                 181 {
176   fRmin1= Rmin1 ;                                 182   fRmin1= Rmin1 ;
177   Initialize();                                   183   Initialize();
178 }                                                 184 }
179                                                   185 
180 inline                                            186 inline
181 void G4Cons::SetOuterRadiusMinusZ( G4double Rm    187 void G4Cons::SetOuterRadiusMinusZ( G4double Rmax1 )
182 {                                                 188 {
183   fRmax1= Rmax1 ;                                 189   fRmax1= Rmax1 ;
184   Initialize();                                   190   Initialize();
185 }                                                 191 }
186                                                   192 
187 inline                                            193 inline
188 void G4Cons::SetInnerRadiusPlusZ ( G4double Rm    194 void G4Cons::SetInnerRadiusPlusZ ( G4double Rmin2 )
189 {                                                 195 {
190   fRmin2= Rmin2 ;                                 196   fRmin2= Rmin2 ;
191   Initialize();                                   197   Initialize();
192 }                                                 198 }
193                                                   199 
194 inline                                            200 inline
195 void G4Cons::SetOuterRadiusPlusZ ( G4double Rm    201 void G4Cons::SetOuterRadiusPlusZ ( G4double Rmax2 )
196 {                                                 202 {
197   fRmax2= Rmax2 ;                                 203   fRmax2= Rmax2 ;
198   Initialize();                                   204   Initialize();
199 }                                                 205 }
200                                                   206 
201 inline                                            207 inline
202 void G4Cons::SetZHalfLength ( G4double newDz )    208 void G4Cons::SetZHalfLength ( G4double newDz )
203 {                                                 209 {
204   fDz= newDz ;                                    210   fDz= newDz ;
205   Initialize();                                   211   Initialize();
206 }                                                 212 }
207                                                   213 
208 inline                                            214 inline
209 void G4Cons::SetStartPhiAngle ( G4double newSP    215 void G4Cons::SetStartPhiAngle ( G4double newSPhi, G4bool compute )
210 {                                                 216 {
211   // Flag 'compute' can be used to explicitely    217   // Flag 'compute' can be used to explicitely avoid recomputation of
212   // trigonometry in case SetDeltaPhiAngle() i    218   // trigonometry in case SetDeltaPhiAngle() is invoked afterwards
213                                                   219 
214   CheckSPhiAngle(newSPhi);                        220   CheckSPhiAngle(newSPhi);
215   fPhiFullCone = false;                           221   fPhiFullCone = false;
216   if (compute)  { InitializeTrigonometry(); }     222   if (compute)  { InitializeTrigonometry(); }
217   Initialize();                                   223   Initialize();
218 }                                                 224 }
219                                                   225 
220 void G4Cons::SetDeltaPhiAngle ( G4double newDP    226 void G4Cons::SetDeltaPhiAngle ( G4double newDPhi )
221 {                                                 227 {
222   CheckPhiAngles(fSPhi, newDPhi);                 228   CheckPhiAngles(fSPhi, newDPhi);
223   Initialize();                                   229   Initialize();
224 }                                                 230 }
225                                                   231 
                                                   >> 232 // Old access methods ...
                                                   >> 233 
                                                   >> 234 inline
                                                   >> 235 G4double G4Cons::GetRmin1() const
                                                   >> 236 {
                                                   >> 237   return GetInnerRadiusMinusZ();
                                                   >> 238 }
                                                   >> 239 
                                                   >> 240 inline
                                                   >> 241 G4double G4Cons::GetRmax1() const
                                                   >> 242 {
                                                   >> 243   return GetOuterRadiusMinusZ();
                                                   >> 244 }
                                                   >> 245 
                                                   >> 246 inline
                                                   >> 247 G4double G4Cons::GetRmin2() const
                                                   >> 248 {
                                                   >> 249   return GetInnerRadiusPlusZ();
                                                   >> 250 }
                                                   >> 251 
                                                   >> 252 inline
                                                   >> 253 G4double G4Cons::GetRmax2() const
                                                   >> 254 {
                                                   >> 255   return GetOuterRadiusPlusZ();
                                                   >> 256 }
                                                   >> 257 
                                                   >> 258 inline  
                                                   >> 259 G4double G4Cons::GetDz() const
                                                   >> 260 {
                                                   >> 261   return GetZHalfLength();
                                                   >> 262 }
                                                   >> 263 
                                                   >> 264 inline
                                                   >> 265 G4double G4Cons::GetSPhi() const
                                                   >> 266 {
                                                   >> 267   return GetStartPhiAngle();
                                                   >> 268 }
                                                   >> 269 
                                                   >> 270 inline
                                                   >> 271 G4double G4Cons::GetDPhi() const
                                                   >> 272 {
                                                   >> 273   return GetDeltaPhiAngle();
                                                   >> 274 }
                                                   >> 275 
226 inline                                            276 inline
227 G4double G4Cons::GetCubicVolume()                 277 G4double G4Cons::GetCubicVolume()
228 {                                                 278 {
229   if(fCubicVolume != 0.) {;}                      279   if(fCubicVolume != 0.) {;}
230   else                                         << 280   else   
231   {                                               281   {
232     G4double Rmean, rMean, deltaR, deltar;        282     G4double Rmean, rMean, deltaR, deltar;
233                                                   283 
234     Rmean  = 0.5*(fRmax1+fRmax2);                 284     Rmean  = 0.5*(fRmax1+fRmax2);
235     deltaR = fRmax1-fRmax2;                       285     deltaR = fRmax1-fRmax2;
236                                                   286 
237     rMean  = 0.5*(fRmin1+fRmin2);                 287     rMean  = 0.5*(fRmin1+fRmin2);
238     deltar = fRmin1-fRmin2;                       288     deltar = fRmin1-fRmin2;
239     fCubicVolume = fDPhi*fDz*(Rmean*Rmean-rMea    289     fCubicVolume = fDPhi*fDz*(Rmean*Rmean-rMean*rMean
240                             +(deltaR*deltaR-de << 290                             +(deltaR*deltaR-deltar*deltar)/12); 
241   }                                               291   }
242   return fCubicVolume;                            292   return fCubicVolume;
243 }                                                 293 }
244                                                   294 
245 inline                                            295 inline
246 G4double G4Cons::GetSurfaceArea()                 296 G4double G4Cons::GetSurfaceArea()
247 {                                                 297 {
248   if(fSurfaceArea != 0.) {;}                      298   if(fSurfaceArea != 0.) {;}
249   else                                         << 299   else   
250   {                                               300   {
251     G4double mmin, mmax, dmin, dmax;              301     G4double mmin, mmax, dmin, dmax;
252                                                   302 
253     mmin= (fRmin1+fRmin2)*0.5;                    303     mmin= (fRmin1+fRmin2)*0.5;
254     mmax= (fRmax1+fRmax2)*0.5;                    304     mmax= (fRmax1+fRmax2)*0.5;
255     dmin= (fRmin2-fRmin1);                        305     dmin= (fRmin2-fRmin1);
256     dmax= (fRmax2-fRmax1);                        306     dmax= (fRmax2-fRmax1);
257                                                << 307     
258     fSurfaceArea = fDPhi*( mmin * std::sqrt(dm    308     fSurfaceArea = fDPhi*( mmin * std::sqrt(dmin*dmin+4*fDz*fDz)
259                          + mmax * std::sqrt(dm    309                          + mmax * std::sqrt(dmax*dmax+4*fDz*fDz)
260                          + 0.5*(fRmax1*fRmax1-    310                          + 0.5*(fRmax1*fRmax1-fRmin1*fRmin1
261                                +fRmax2*fRmax2-    311                                +fRmax2*fRmax2-fRmin2*fRmin2 ));
262     if(!fPhiFullCone)                             312     if(!fPhiFullCone)
263     {                                             313     {
264       fSurfaceArea = fSurfaceArea+4*fDz*(mmax-    314       fSurfaceArea = fSurfaceArea+4*fDz*(mmax-mmin);
265     }                                             315     }
266   }                                               316   }
267   return fSurfaceArea;                            317   return fSurfaceArea;
268 }                                                 318 }
269                                                   319