Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/solids/specific/src/G4EnclosingCylinder.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/specific/src/G4EnclosingCylinder.cc (Version 11.3.0) and /geometry/solids/specific/src/G4EnclosingCylinder.cc (Version 11.0.p4)


  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 G4EnclosingCylinder, a ut     26 // Implementation of G4EnclosingCylinder, a utility class
 27 // for a quick check of geometry.                  27 // for a quick check of geometry.
 28 //                                                 28 //
 29 // Author: David C. Williams (davidw@scipp.ucs     29 // Author: David C. Williams (davidw@scipp.ucsc.edu)
 30 // -------------------------------------------     30 // --------------------------------------------------------------------
 31                                                    31 
 32 #include "G4EnclosingCylinder.hh"                  32 #include "G4EnclosingCylinder.hh"
 33 #include "G4PhysicalConstants.hh"                  33 #include "G4PhysicalConstants.hh"
 34 #include "G4ReduciblePolygon.hh"                   34 #include "G4ReduciblePolygon.hh"
 35 #include "G4GeometryTolerance.hh"                  35 #include "G4GeometryTolerance.hh"
 36                                                    36 
 37 // Constructor                                     37 // Constructor
 38 //                                                 38 //
 39 G4EnclosingCylinder::G4EnclosingCylinder( cons     39 G4EnclosingCylinder::G4EnclosingCylinder( const G4ReduciblePolygon* rz,
 40                                                    40                                                 G4bool thePhiIsOpen, 
 41                                                    41                                                 G4double theStartPhi,
 42                                                    42                                                 G4double theTotalPhi )
 43   : startPhi(theStartPhi), totalPhi(theTotalPh     43   : startPhi(theStartPhi), totalPhi(theTotalPhi),
 44     concave(theTotalPhi > pi)                      44     concave(theTotalPhi > pi)
 45 {                                                  45 {
 46   //                                               46   //
 47   // Obtain largest r and smallest and largest     47   // Obtain largest r and smallest and largest z
 48   //                                               48   //
 49   radius = rz->Amax();                             49   radius = rz->Amax();
 50   zHi = rz->Bmax();                                50   zHi = rz->Bmax();
 51   zLo = rz->Bmin();                                51   zLo = rz->Bmin();
 52                                                    52 
 53   G4double kCarTolerance = G4GeometryTolerance     53   G4double kCarTolerance = G4GeometryTolerance::GetInstance()
 54                            ->GetSurfaceToleran     54                            ->GetSurfaceTolerance();
 55   //                                               55   //
 56   // Save phi info                                 56   // Save phi info
 57   //                                               57   //
 58   phiIsOpen = thePhiIsOpen;                        58   phiIsOpen = thePhiIsOpen;
 59   if ( phiIsOpen )                                 59   if ( phiIsOpen )
 60   {                                                60   {    
 61     rx1 = std::cos(startPhi);                      61     rx1 = std::cos(startPhi);
 62     ry1 = std::sin(startPhi);                      62     ry1 = std::sin(startPhi);
 63     dx1 = +ry1*10*kCarTolerance;                   63     dx1 = +ry1*10*kCarTolerance;
 64     dy1 = -rx1*10*kCarTolerance;                   64     dy1 = -rx1*10*kCarTolerance;
 65                                                    65     
 66     rx2 = std::cos(startPhi+totalPhi);             66     rx2 = std::cos(startPhi+totalPhi);
 67     ry2 = std::sin(startPhi+totalPhi);             67     ry2 = std::sin(startPhi+totalPhi);
 68     dx2 = -ry2*10*kCarTolerance;                   68     dx2 = -ry2*10*kCarTolerance;
 69     dy2 = +rx2*10*kCarTolerance;                   69     dy2 = +rx2*10*kCarTolerance;
 70   }                                                70   }
 71                                                    71   
 72   //                                               72   //
 73   // Add safety                                    73   // Add safety
 74   //                                               74   //
 75   radius += 10*kCarTolerance;                      75   radius += 10*kCarTolerance;
 76   zLo    -= 10*kCarTolerance;                      76   zLo    -= 10*kCarTolerance;
 77   zHi    += 10*kCarTolerance;                      77   zHi    += 10*kCarTolerance;
 78 }                                                  78 }
 79                                                    79 
 80 // Fake default constructor - sets only member     80 // Fake default constructor - sets only member data and allocates memory
 81 //                            for usage restri     81 //                            for usage restricted to object persistency.
 82 //                                                 82 //
 83 G4EnclosingCylinder::G4EnclosingCylinder( __vo     83 G4EnclosingCylinder::G4EnclosingCylinder( __void__& )
 84 : radius(0.), zLo(0.), zHi(0.), phiIsOpen(fals <<  84 : radius(0.), zLo(0.), zHi(0.), phiIsOpen(0.), startPhi(0.), totalPhi(0.),
 85   concave(false)                                   85   concave(false)
                                                   >>  86 {
                                                   >>  87 }
                                                   >>  88 
                                                   >>  89 // Destructor
                                                   >>  90 //
                                                   >>  91 G4EnclosingCylinder::~G4EnclosingCylinder()
 86 {                                                  92 {
 87 }                                                  93 }
 88                                                    94 
 89 // Outside                                         95 // Outside
 90 //                                                 96 //
 91 // Decide very rapidly if the point is outside     97 // Decide very rapidly if the point is outside the cylinder
 92 //                                                 98 //
 93 // If one is not certain, return false             99 // If one is not certain, return false
 94 //                                                100 //
 95 G4bool G4EnclosingCylinder::MustBeOutside( con    101 G4bool G4EnclosingCylinder::MustBeOutside( const G4ThreeVector& p ) const
 96 {                                                 102 {
 97   if (p.perp() > radius) return true;             103   if (p.perp() > radius) return true;
 98   if (p.z() < zLo) return true;                   104   if (p.z() < zLo) return true;
 99   if (p.z() > zHi) return true;                   105   if (p.z() > zHi) return true;
100                                                   106 
101   if (phiIsOpen)                                  107   if (phiIsOpen)
102   {                                               108   {
103     if (concave)                                  109     if (concave)
104     {                                             110     {
105       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1)    111       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1) < 0) return false;
106       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2)    112       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2) > 0) return false;
107     }                                             113     }
108     else                                          114     else
109     {                                             115     {
110       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1)    116       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1) > 0) return true;
111       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2)    117       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2) < 0) return true;
112     }                                             118     }
113   }                                               119   }
114                                                   120   
115   return false;                                   121   return false;
116 }                                                 122 }
117                                                   123 
118 // Misses                                         124 // Misses
119 //                                                125 //
120 // Decide very rapidly if the trajectory is go    126 // Decide very rapidly if the trajectory is going to miss the cylinder
121 //                                                127 //
122 // If one is not sure, return false               128 // If one is not sure, return false
123 //                                                129 //
124 G4bool G4EnclosingCylinder::ShouldMiss( const     130 G4bool G4EnclosingCylinder::ShouldMiss( const G4ThreeVector& p,
125                                         const     131                                         const G4ThreeVector& v ) const
126 {                                                 132 {
127   if (!MustBeOutside(p)) return false;            133   if (!MustBeOutside(p)) return false;
128                                                   134   
129   G4double cross = p.x()*v.y() - p.y()*v.x();     135   G4double cross = p.x()*v.y() - p.y()*v.x();
130   if (cross > radius) return true;                136   if (cross > radius) return true;
131                                                   137   
132   if (p.perp() > radius)                          138   if (p.perp() > radius)
133   {                                               139   {
134     G4double dot = p.x()*v.x() + p.y()*v.y();     140     G4double dot = p.x()*v.x() + p.y()*v.y();
135     if (dot > 0) return true;                     141     if (dot > 0) return true;
136   }                                               142   }
137                                                   143 
138   return false;                                   144   return false;
139 }                                                 145 }
140                                                   146