| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 23
27 #ifndef G4RKPropagation_h 24 #ifndef G4RKPropagation_h
28 #define G4RKPropagation_h 1 25 #define G4RKPropagation_h 1
29 26
30 #include "G4VFieldPropagation.hh" 27 #include "G4VFieldPropagation.hh"
31 #include "G4VNuclearField.hh" 28 #include "G4VNuclearField.hh"
32 #include "G4V3DNucleus.hh" 29 #include "G4V3DNucleus.hh"
33 #include "G4KM_DummyField.hh" 30 #include "G4KM_DummyField.hh"
34 #include "G4Mag_EqRhs.hh" 31 #include "G4Mag_EqRhs.hh"
35 #include <map> 32 #include <map>
36 33
37 34
38 class G4RKPropagation: public G4VFieldPropagat 35 class G4RKPropagation: public G4VFieldPropagation
39 { 36 {
40 37
41 public: 38 public:
42 G4RKPropagation(); 39 G4RKPropagation();
43 virtual ~G4RKPropagation(); 40 virtual ~G4RKPropagation();
44 41
45 private: 42 private:
46 G4RKPropagation(const G4RKPropagation &righ 43 G4RKPropagation(const G4RKPropagation &right);
47 const G4RKPropagation & operator=(const G4RK 44 const G4RKPropagation & operator=(const G4RKPropagation & right);
48 G4bool operator==(const G4RKPropagation & ri << 45 G4int operator==(const G4RKPropagation & right) const;
49 G4bool operator!=(const G4RKPropagation & ri << 46 G4int operator!=(const G4RKPropagation & right) const;
50 47
51 public: 48 public:
52 49
53 virtual void Init(G4V3DNucleus * nucleus); 50 virtual void Init(G4V3DNucleus * nucleus);
54 virtual void Transport(G4KineticTrackVector 51 virtual void Transport(G4KineticTrackVector &theActive,
55 const G4KineticTrackVector &theSpectato 52 const G4KineticTrackVector &theSpectators,
56 G4double theTimeStep); 53 G4double theTimeStep);
57 G4bool GetSphereIntersectionTimes(const G4Ki 54 G4bool GetSphereIntersectionTimes(const G4KineticTrack * track,
58 G4double & t1, G4double & t2); 55 G4double & t1, G4double & t2);
59 G4ThreeVector GetMomentumTransfer() const; 56 G4ThreeVector GetMomentumTransfer() const;
60 private: 57 private:
61 G4double theOuterRadius; 58 G4double theOuterRadius;
62 G4V3DNucleus * theNucleus; 59 G4V3DNucleus * theNucleus;
63 std::map <G4int, G4VNuclearField *, std::les 60 std::map <G4int, G4VNuclearField *, std::less<G4int> > * theFieldMap;
64 std::map <G4int, G4Mag_EqRhs *, std::less<G4 61 std::map <G4int, G4Mag_EqRhs *, std::less<G4int> > * theEquationMap;
65 G4KM_DummyField * theField; 62 G4KM_DummyField * theField;
66 63
67 G4ThreeVector theMomentumTranfer; 64 G4ThreeVector theMomentumTranfer;
68 65
69 G4bool GetSphereIntersectionTimes(const G4do 66 G4bool GetSphereIntersectionTimes(const G4double radius,
70 const G4ThreeVector & currentPos, 67 const G4ThreeVector & currentPos,
71 const G4LorentzVector & momentum, 68 const G4LorentzVector & momentum,
72 G4double & t1, G4double & t2); 69 G4double & t1, G4double & t2);
73 // implementation 70 // implementation
74 71
75 G4bool FieldTransport(G4KineticTrack * track 72 G4bool FieldTransport(G4KineticTrack * track, const G4double timestep);
76 G4bool FreeTransport(G4KineticTrack * track, 73 G4bool FreeTransport(G4KineticTrack * track, const G4double timestep);
77 74
78 void delete_FieldsAndMap( 75 void delete_FieldsAndMap(
79 std::map <G4int, G4VNuclearField *, std::les 76 std::map <G4int, G4VNuclearField *, std::less<G4int> > * aMap);
80 void delete_EquationsAndMap( 77 void delete_EquationsAndMap(
81 std::map <G4int, G4Mag_EqRhs *, std::less<G4 78 std::map <G4int, G4Mag_EqRhs *, std::less<G4int> > * aMap);
82 79
83 public: 80 public:
84 inline G4double GetBarrier(G4int encoding) 81 inline G4double GetBarrier(G4int encoding)
85 { 82 {
86 std::map <G4int, G4VNuclearField *, std::les 83 std::map <G4int, G4VNuclearField *, std::less<G4int> >::iterator iter;
87 iter = theFieldMap->find(encoding); 84 iter = theFieldMap->find(encoding);
88 if(iter == theFieldMap->end()) return 0; 85 if(iter == theFieldMap->end()) return 0;
89 return (*theFieldMap)[encoding]->GetBarrier( 86 return (*theFieldMap)[encoding]->GetBarrier();
90 } 87 }
91 88
92 inline G4double GetField(G4int encoding,G4Th 89 inline G4double GetField(G4int encoding,G4ThreeVector pos)
93 { 90 {
94 std::map <G4int, G4VNuclearField *, std::les 91 std::map <G4int, G4VNuclearField *, std::less<G4int> >::iterator iter;
95 iter = theFieldMap->find(encoding); 92 iter = theFieldMap->find(encoding);
96 if(iter == theFieldMap->end()) return 0; 93 if(iter == theFieldMap->end()) return 0;
97 return (*theFieldMap)[encoding]->GetField(po 94 return (*theFieldMap)[encoding]->GetField(pos);
98 } 95 }
99 96
100 }; 97 };
101 98
102 #endif 99 #endif
103 100
104 101
105 102
106 103
107 104
108 105
109 106
110 107
111 108
112 109
113 110
114 111
115 112
116 113
117 114
118 115