Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4INCLIntersection.hh
Go to the documentation of this file.
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// INCL++ intra-nuclear cascade model
27// Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
28// Davide Mancusi, CEA
29// Alain Boudard, CEA
30// Sylvie Leray, CEA
31// Joseph Cugnon, University of Liege
32//
33// INCL++ revision: v5.1.8
34//
35#define INCLXX_IN_GEANT4_MODE 1
36
37#include "globals.hh"
38
39/** \file G4INCLIntersection.hh
40 * \brief Simple class for computing intersections between a straight line and a sphere.
41 *
42 * \date 12 December 2011
43 * \author Davide Mancusi
44 */
45
46#ifndef G4INCLINTERSECTION_HH
47#define G4INCLINTERSECTION_HH 1
48
49#include "G4INCLThreeVector.hh"
50#include <utility>
51
52namespace G4INCL {
53 /** \brief Intersection-point structure
54 *
55 * The structure contains the time and position of the intersection point
56 * of a trajectory with a surface, if it exists.
57 */
58 struct Intersection {
59 Intersection(const G4bool e, const G4double t, const ThreeVector &p) :
60 exists(e), time(t), position(p) {}
64 };
65
67 public:
68 /** \brief Compute the first intersection of a straight particle
69 * trajectory with a sphere.
70 *
71 * \param x0 the starting position of the trajectory
72 * \param p the trajectory direction
73 * \param r the radius of the sphere (centred in the origin)
74 * \return an Intersection. The G4bool is true if an intersection exists,
75 * in which case its position is stored in the ThreeVector and
76 * its time in the G4double.
77 */
79 return getTrajectoryIntersection(x0, p, r, true);
80 }
81 static inline Intersection getLaterTrajectoryIntersection(const ThreeVector &x0, const ThreeVector &p, const G4double r) {
82 return getTrajectoryIntersection(x0, p, r, false);
83 }
84 static inline std::pair<Intersection,Intersection> getTrajectoryIntersections(const ThreeVector &x0, const ThreeVector &p, const G4double r) {
85 return std::make_pair(
86 getTrajectoryIntersection(x0, p, r, true),
87 getTrajectoryIntersection(x0, p, r, false)
88 );
89 }
90 private:
91 static Intersection getTrajectoryIntersection(const ThreeVector &x0, const ThreeVector &p, const G4double r, const G4bool earliest);
92 };
93}
94
95#endif /* G4INCLINTERSECTION_HH */
double G4double
Definition: G4Types.hh:64
bool G4bool
Definition: G4Types.hh:67
static Intersection getLaterTrajectoryIntersection(const ThreeVector &x0, const ThreeVector &p, const G4double r)
static Intersection getEarlierTrajectoryIntersection(const ThreeVector &x0, const ThreeVector &p, const G4double r)
Compute the first intersection of a straight particle trajectory with a sphere.
static std::pair< Intersection, Intersection > getTrajectoryIntersections(const ThreeVector &x0, const ThreeVector &p, const G4double r)
Intersection-point structure.
Intersection(const G4bool e, const G4double t, const ThreeVector &p)