| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // -------------------------------------------
28 // GEANT 4 class implementation file
29 // -------------------------------------------
30 //
31
32 #include "G4ErrorTrackLengthTarget.hh"
33
34 #include "G4ParticleTable.hh"
35 #include "G4ParticleDefinition.hh"
36 #include "G4VProcess.hh"
37 #include "G4ProcessVector.hh"
38 #include "G4ProcessManager.hh"
39
40 #ifdef G4VERBOSE
41 # include "G4ErrorPropagatorData.hh" //for v
42 #endif
43
44 //--------------------------------------------
45 G4ErrorTrackLengthTarget::G4ErrorTrackLengthTa
46 : G4VDiscreteProcess("G4ErrorTrackLengthTarg
47 , theMaximumTrackLength(maxTrkLength)
48 {
49 theType = G4ErrorTarget_TrkL;
50
51 G4ParticleTable::G4PTblDicIterator* theParti
52 G4ParticleTable::GetParticleTable()->GetIt
53
54 // loop over all particles in G4ParticleTabl
55
56 theParticleIterator->reset();
57 while((*theParticleIterator)()) // Loop che
58 {
59 G4ParticleDefinition* particle = thePartic
60 G4ProcessManager* pmanager = particle-
61 if(!particle->IsShortLived())
62 {
63 // Add transportation process for all pa
64 if(pmanager == 0)
65 {
66 // Error !! no process manager
67 G4String particleName = particle->GetP
68 G4Exception("G4ErrorTrackLengthTarget:
69 "No process manager", RunM
70 }
71 else
72 {
73 G4ProcessVector* procvec = pmanager->G
74 std::size_t isiz = (G4int)procvec->siz
75
76 for(G4int ii = 0; ii < (G4int)isiz; ++
77 {
78 if(((*procvec)[ii])->GetProcessName(
79 {
80 pmanager->RemoveProcess((*procvec)
81 }
82 }
83 pmanager->AddDiscreteProcess(this, 4);
84 isiz = procvec->size();
85 }
86 }
87 else
88 {
89 // shortlived particle case
90 }
91 }
92 }
93
94 //--------------------------------------------
95 G4double G4ErrorTrackLengthTarget::PostStepGet
96 const G4Track& track, G4double, G4ForceCondi
97 {
98 *condition = NotForced;
99 return GetMeanFreePath(track, 0., condition)
100 }
101
102 //--------------------------------------------
103 G4double G4ErrorTrackLengthTarget::GetMeanFree
104
105
106 {
107 #ifdef G4VERBOSE
108 if(G4ErrorPropagatorData::verbose() >= 3)
109 {
110 G4cout << " G4ErrorTrackLengthTarget::GetM
111 << theMaximumTrackLength - track.Ge
112 }
113 #endif
114
115 return theMaximumTrackLength - track.GetTrac
116 }
117
118 G4VParticleChange* G4ErrorTrackLengthTarget::P
119
120 {
121 theParticleChange.Initialize(aTrack);
122 return &theParticleChange;
123 }
124
125 //--------------------------------------------
126 void G4ErrorTrackLengthTarget::Dump(const G4St
127 {
128 G4cout << msg << "G4ErrorTrackLengthTarget:
129 << theMaximumTrackLength << G4endl;
130 }
131