ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
RunAction.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file RunAction.cc
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
//
28
//
29
//
30
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
33
#include "RunAction.hh"
34
#include "DetectorConstruction.hh"
35
#include "PhysicsList.hh"
36
#include "StepMax.hh"
37
#include "PrimaryGeneratorAction.hh"
38
#include "Run.hh"
39
#include "HistoManager.hh"
40
41
#include "
G4EmCalculator.hh
"
42
#include "
G4EmParameters.hh
"
43
44
#include "
Randomize.hh
"
45
46
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47
48
RunAction::RunAction
(
DetectorConstruction
* det,
PhysicsList
* phys,
49
PrimaryGeneratorAction
* kin)
50
:
G4UserRunAction
(),
51
fDetector(det),fPhysics(phys),fPrimary(kin),fRun(0),fHistoManager(0)
52
{
53
// Book predefined histograms
54
fHistoManager
=
new
HistoManager
();
55
}
56
57
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58
59
RunAction::~RunAction
()
60
{
61
delete
fHistoManager
;
62
}
63
64
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65
66
G4Run
*
RunAction::GenerateRun
()
67
{
68
fRun
=
new
Run
(
fDetector
);
69
return
fRun
;
70
}
71
72
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73
74
void
RunAction::BeginOfRunAction
(
const
G4Run
*)
75
{
76
// show Rndm status
77
if
(
isMaster
) {
78
G4Random::showEngineStatus();
79
G4EmParameters::Instance
()->
Dump
();
80
}
81
//histograms
82
//
83
G4AnalysisManager
* analysisManager =
G4AnalysisManager::Instance
();
84
if
( analysisManager->
IsActive
() ) {
85
analysisManager->
OpenFile
();
86
}
87
88
if
(!
fPrimary
)
return
;
89
90
// keep run condition
91
G4ParticleDefinition
*
particle
92
=
fPrimary
->
GetParticleGun
()->
GetParticleDefinition
();
93
G4double
energy
=
fPrimary
->
GetParticleGun
()->
GetParticleEnergy
();
94
fRun
->
SetPrimary
(particle, energy);
95
96
//get CsdaRange from EmCalculator
97
//
98
G4EmCalculator
emCalculator;
99
G4Material
*
material
=
fDetector
->
GetAbsorMaterial
();
100
101
G4double
csdaRange =
DBL_MAX
;
102
if
(particle->
GetPDGCharge
() != 0.) {
103
csdaRange = emCalculator.
GetCSDARange
(energy,particle,material);
104
fRun
->
SetCsdaRange
(csdaRange);
105
}
106
107
//set StepMax from histos 1 and 8
108
//
109
G4double
stepMax = csdaRange;
110
G4int
ih = 1;
111
if
(analysisManager->
GetH1Activation
(ih))
112
stepMax = analysisManager->
GetH1Width
(ih)*analysisManager->
GetH1Unit
(ih);
113
ih = 8;
114
if
(analysisManager->
GetH1Activation
(ih)) {
115
G4double
width
= analysisManager->
GetH1Width
(ih);
116
stepMax =
std::min
(stepMax, width*csdaRange);
117
}
118
fPhysics
->
GetStepMaxProcess
()->SetMaxStep2(stepMax);
119
122
}
123
124
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125
126
void
RunAction::EndOfRunAction
(
const
G4Run
*)
127
{
128
if
(
isMaster
)
fRun
->
EndOfRun
();
129
130
// save histograms
131
G4AnalysisManager
* analysisManager =
G4AnalysisManager::Instance
();
132
if
( analysisManager->
IsActive
() ) {
133
analysisManager->
Write
();
134
analysisManager->
CloseFile
();
135
}
136
137
// show Rndm status
138
if
(
isMaster
) G4Random::showEngineStatus();
139
}
140
141
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
geant4
tree
geant4-10.6-release
examples
extended
electromagnetic
TestEm12
src
RunAction.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:52
using
1.8.2 with
ECCE GitHub integration