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 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 // Implementation of G4EnclosingCylinder, a utility class
 27 // for a quick check of geometry.
 28 //
 29 // Author: David C. Williams (davidw@scipp.ucsc.edu)
 30 // --------------------------------------------------------------------
 31 
 32 #include "G4EnclosingCylinder.hh"
 33 #include "G4PhysicalConstants.hh"
 34 #include "G4ReduciblePolygon.hh"
 35 #include "G4GeometryTolerance.hh"
 36 
 37 // Constructor
 38 //
 39 G4EnclosingCylinder::G4EnclosingCylinder( const G4ReduciblePolygon* rz,
 40                                                 G4bool thePhiIsOpen, 
 41                                                 G4double theStartPhi,
 42                                                 G4double theTotalPhi )
 43   : startPhi(theStartPhi), totalPhi(theTotalPhi),
 44     concave(theTotalPhi > pi)
 45 {
 46   //
 47   // Obtain largest r and smallest and largest z
 48   //
 49   radius = rz->Amax();
 50   zHi = rz->Bmax();
 51   zLo = rz->Bmin();
 52 
 53   G4double kCarTolerance = G4GeometryTolerance::GetInstance()
 54                            ->GetSurfaceTolerance();
 55   //
 56   // Save phi info
 57   //
 58   phiIsOpen = thePhiIsOpen;
 59   if ( phiIsOpen )
 60   {    
 61     rx1 = std::cos(startPhi);
 62     ry1 = std::sin(startPhi);
 63     dx1 = +ry1*10*kCarTolerance;
 64     dy1 = -rx1*10*kCarTolerance;
 65     
 66     rx2 = std::cos(startPhi+totalPhi);
 67     ry2 = std::sin(startPhi+totalPhi);
 68     dx2 = -ry2*10*kCarTolerance;
 69     dy2 = +rx2*10*kCarTolerance;
 70   }
 71   
 72   //
 73   // Add safety
 74   //
 75   radius += 10*kCarTolerance;
 76   zLo    -= 10*kCarTolerance;
 77   zHi    += 10*kCarTolerance;
 78 }
 79 
 80 // Fake default constructor - sets only member data and allocates memory
 81 //                            for usage restricted to object persistency.
 82 //
 83 G4EnclosingCylinder::G4EnclosingCylinder( __void__& )
 84 : radius(0.), zLo(0.), zHi(0.), phiIsOpen(false), startPhi(0.), totalPhi(0.),
 85   concave(false)
 86 {
 87 }
 88 
 89 // Outside
 90 //
 91 // Decide very rapidly if the point is outside the cylinder
 92 //
 93 // If one is not certain, return false
 94 //
 95 G4bool G4EnclosingCylinder::MustBeOutside( const G4ThreeVector& p ) const
 96 {
 97   if (p.perp() > radius) return true;
 98   if (p.z() < zLo) return true;
 99   if (p.z() > zHi) return true;
100 
101   if (phiIsOpen)
102   {
103     if (concave)
104     {
105       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1) < 0) return false;
106       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2) > 0) return false;
107     }
108     else
109     {
110       if ( ((p.x()-dx1)*ry1 - (p.y()-dy1)*rx1) > 0) return true;
111       if ( ((p.x()-dx2)*ry2 - (p.y()-dy2)*rx2) < 0) return true;
112     }
113   }
114   
115   return false;
116 }
117 
118 // Misses
119 //
120 // Decide very rapidly if the trajectory is going to miss the cylinder
121 //
122 // If one is not sure, return false
123 //
124 G4bool G4EnclosingCylinder::ShouldMiss( const G4ThreeVector& p,
125                                         const G4ThreeVector& v ) const
126 {
127   if (!MustBeOutside(p)) return false;
128   
129   G4double cross = p.x()*v.y() - p.y()*v.x();
130   if (cross > radius) return true;
131   
132   if (p.perp() > radius)
133   {
134     G4double dot = p.x()*v.x() + p.y()*v.y();
135     if (dot > 0) return true;
136   }
137 
138   return false;
139 }
140