source: trunk/examples/extended/electromagnetic/TestEm13/src/RunAction.cc@ 1357

Last change on this file since 1357 was 1337, checked in by garnier, 15 years ago

tag geant4.9.4 beta 1 + modifs locales
File size: 6.1 KB
Line 
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// $Id: RunAction.cc,v 1.8 2010/04/05 13:45:17 maire Exp $
27// GEANT4 tag $Name: geant4-09-04-beta-01 $
28//
29//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
32#include "RunAction.hh"
33
34#include "DetectorConstruction.hh"
35#include "PrimaryGeneratorAction.hh"
36
37#include "G4Run.hh"
38#include "G4RunManager.hh"
39#include "G4UnitsTable.hh"
40#include "G4EmCalculator.hh"
41#include "G4Gamma.hh"
42
43#include "Randomize.hh"
44#include <iomanip>
45
46//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47
48RunAction::RunAction(DetectorConstruction* det, PrimaryGeneratorAction* prim)
49 : detector(det), primary(prim)
50{ }
51
52//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53
54RunAction::~RunAction()
55{ }
56
57//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58
59void RunAction::BeginOfRunAction(const G4Run* aRun)
60{
61 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
62
63 // save Rndm status
64 G4RunManager::GetRunManager()->SetRandomNumberStore(false);
65 CLHEP::HepRandom::showEngineStatus();
66}
67
68//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69
70void RunAction::EndOfRunAction(const G4Run* aRun)
71{
72 G4int NbOfEvents = aRun->GetNumberOfEvent();
73 if (NbOfEvents == 0) return;
74
75 G4int prec = G4cout.precision(5);
76
77 G4Material* material = detector->GetMaterial();
78 G4double density = material->GetDensity();
79 G4double tickness = detector->GetSize();
80
81 G4ParticleDefinition* particle =
82 primary->GetParticleGun()->GetParticleDefinition();
83 G4String Particle = particle->GetParticleName();
84 G4double energy = primary->GetParticleGun()->GetParticleEnergy();
85 G4cout << "\n The run consists of " << NbOfEvents << " "<< Particle << " of "
86 << G4BestUnit(energy,"Energy") << " through "
87 << G4BestUnit(tickness,"Length") << " of "
88 << material->GetName() << " (density: "
89 << G4BestUnit(density,"Volumic Mass") << ")" << G4endl;
90
91 //frequency of processes
92 G4int totalCount = 0;
93 G4int survive = 0;
94 G4cout << "\n Process calls frequency --->";
95 std::map<G4String,G4int>::iterator it;
96 for (it = procCounter.begin(); it != procCounter.end(); it++) {
97 G4String procName = it->first;
98 G4int count = it->second;
99 totalCount += count;
100 G4cout << "\t" << procName << " = " << count;
101 if (procName == "Transportation") survive = count;
102 }
103 G4cout << G4endl;
104 if (totalCount == 0) return;
105
106 G4double ratio = double(survive)/totalCount;
107
108 G4cout << "\n Nb of incident particles unaltered after "
109 << G4BestUnit(tickness,"Length") << " of "
110 << material->GetName() << " : " << survive
111 << " over " << totalCount << " incident particles."
112 << " Ratio = " << 100*ratio << " %" << G4endl;
113
114 if (ratio == 0.) return;
115
116 //compute cross section and related quantities
117 //
118 G4double CrossSection = - std::log(ratio)/tickness;
119 G4double massicCS = CrossSection/density;
120
121 G4cout << " ---> CrossSection per volume:\t" << CrossSection*cm << " cm^-1 "
122 << "\tCrossSection per mass: " << G4BestUnit(massicCS, "Surface/Mass")
123 << G4endl;
124
125 //check cross section from G4EmCalculator
126 //
127 G4cout << "\n Verification from G4EmCalculator: \n";
128 G4EmCalculator emCalculator;
129 G4double sumc = 0.0;
130 for (it = procCounter.begin(); it != procCounter.end(); it++) {
131 G4String procName = it->first;
132 G4double massSigma =
133 emCalculator.GetCrossSectionPerVolume(energy,particle,
134 procName,material)/density;
135 if (particle == G4Gamma::Gamma())
136 massSigma =
137 emCalculator.ComputeCrossSectionPerVolume(energy,particle,
138 procName,material)/density;
139 sumc += massSigma;
140 if (procName != "Transportation")
141 G4cout << "\t" << procName << "= "
142 << G4BestUnit(massSigma, "Surface/Mass");
143 }
144 G4cout << "\ttotal= "
145 << G4BestUnit(sumc, "Surface/Mass") << G4endl;
146
147 //expected ratio of transmitted particles
148 G4double Ratio = std::exp(-sumc*density*tickness);
149 G4cout << "\tExpected ratio of transmitted particles= "
150 << 100*Ratio << " %" << G4endl;
151
152 //restore default format
153 G4cout.precision(prec);
154
155 // remove all contents in procCounter
156 procCounter.clear();
157
158 // show Rndm status
159 CLHEP::HepRandom::showEngineStatus();
160}
161
162//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.