| 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 #include "F01RunAction.hh"
27
28 #include "G4CoupledTransportation.hh"
29 #include "G4Electron.hh"
30 #include "G4ParticleDefinition.hh"
31 #include "G4ProcessManager.hh"
32 #include "G4Run.hh"
33 #include "G4SystemOfUnits.hh"
34 #include "G4Transportation.hh"
35 #include "globals.hh"
36
37 //....oooOO0OOooo........oooOO0OOooo........oo
38
39 void F01RunAction::BeginOfRunAction(const G4Ru
40 {
41 G4cout << "### Run " << aRun->GetRunID() <<
42
43 G4cout << " Calling F01RunAction::ChangeLoop
44 ChangeLooperParameters(G4Electron::Definitio
45 }
46
47 //....oooOO0OOooo........oooOO0OOooo........oo
48
49 void F01RunAction::ChangeLooperParameters(cons
50 {
51 if (particleDef == nullptr) particleDef = G4
52 auto transport = FindTransportation(particle
53
54 // Note that all particles that share the sa
55 // (or Coupled Transportation) process (wit
56 // will be affected by this change
57 // -- to change these parameters by partic
58 // must instead be changed/reset in a (
59
60 G4cout << " ChangeLooperParameters called wi
61 << " transport process: ";
62 if (transport != nullptr) {
63 G4cout << transport->GetProcessName();
64 }
65 else {
66 G4cout << " UNKNOWN -- it is neither G4Tra
67 }
68 G4cout << G4endl;
69
70 if (transport != nullptr) {
71 if (fWarningEnergy >= 0.0) {
72 transport->SetThresholdWarningEnergy(fWa
73 G4cout << "-- Changed Threshold Warning
74 << fWarningEnergy / CLHEP::MeV <<
75 }
76 if (fImportantEnergy >= 0.0) {
77 transport->SetThresholdImportantEnergy(f
78
79 G4cout << "-- Changed Threshold Importan
80 << fImportantEnergy / CLHEP::MeV
81 }
82
83 if (fNumberOfTrials > 0) {
84 transport->SetThresholdTrials(fNumberOfT
85
86 G4cout << "-- Changed number of Trials (
87 }
88 }
89
90 if (transport == nullptr) {
91 if (fWarningEnergy >= 0.0)
92 G4cerr << " Unknown transport process> C
93 if (fImportantEnergy >= 0.0)
94 G4cerr << " Unknown transport process> C
95 if (fNumberOfTrials > 0)
96 G4cerr << " Unknown transport process> C
97 }
98 }
99
100 //....oooOO0OOooo........oooOO0OOooo........oo
101
102 void F01RunAction::EndOfRunAction(const G4Run*
103 {
104 if (fVerboseLevel > 1)
105 G4cout << G4endl << G4endl << " ##########
106 << " ########### " << G4endl << " #
107 << " ####################### " << G
108
109 auto transport = FindTransportation(G4Electr
110 if (transport) {
111 transport->PrintStatistics(G4cout);
112 }
113 }
114
115 //....oooOO0OOooo........oooOO0OOooo........oo
116
117 G4Transportation* F01RunAction::FindTransporta
118
119 {
120 const auto* partPM = particleDef->GetProcess
121
122 G4VProcess* partTransport = partPM->GetProce
123 auto transport = dynamic_cast<G4Transportati
124
125 if (!transport) {
126 partTransport = partPM->GetProcess("Couple
127 auto coupledTransport = dynamic_cast<G4Cou
128 if (coupledTransport) {
129 transport = coupledTransport;
130 }
131 else {
132 partTransport = partPM->GetProcess("Tran
133 auto transportWithMsc = dynamic_cast<G4C
134 if (transportWithMsc) {
135 transport = transportWithMsc;
136 }
137 }
138 }
139
140 if (reportError && !transport) {
141 G4cerr << "Unable to find Transportation p
142 << particleDef->GetParticleName() <
143 << " ) " << G4endl;
144 }
145
146 return transport;
147 }
148