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
//
31
//
32
33
#include "RunAction.hh"
34
#include "RunMessenger.hh"
35
#include "Analysis.hh"
36
37
#include "
G4Run.hh
"
38
#include "
G4ios.hh
"
39
40
#include "
G4PhysicalConstants.hh
"
41
#include "
G4SystemOfUnits.hh
"
42
#include "
Randomize.hh
"
43
44
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
45
46
RunAction::RunAction
()
47
:
G4UserRunAction
(),
48
fRunMessenger(0), fRndmFreq(1)
49
{
50
fRunMessenger
=
new
RunMessenger
(
this
);
51
52
BookHisto
();
53
}
54
55
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
56
57
RunAction::~RunAction
()
58
{
59
delete
fRunMessenger
;
60
}
61
62
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
63
64
void
RunAction::BookHisto
()
65
{
66
// Create or get analysis manager
67
// The choice of analysis technology is done via selection of a namespace
68
// in HistoManager.hh
69
G4AnalysisManager
* analysisManager =
G4AnalysisManager::Instance
();
70
analysisManager->
SetFileName
(
"testem10"
);
71
analysisManager->
SetVerboseLevel
(1);
72
analysisManager->
SetFirstHistoId
(1);
// start histogram numbering from 1
73
analysisManager->
SetActivation
(
true
);
// enable inactivation of histograms
74
75
// Define histograms start values
76
const
G4int
kMaxHisto
= 5;
77
const
G4String
id
[] =
78
{
"1"
,
"2"
,
"3"
,
"4"
,
"5"
};
79
const
G4String
title
[] =
80
{
"Edep"
,
// 1
81
"XTR Gamma spectrum"
,
// 2
82
"Secondary Gamma spectrum"
,
// 3
83
"Secondary e- spectrum"
,
// 4
84
"Edep.old"
};
// 5
85
// Default values (to be reset via /analysis/h1/set command)
86
G4int
nbins = 100;
87
G4double
vmin = 0.;
88
G4double
vmax = 100.;
89
90
// Create all histograms as inactivated
91
// as we have not yet set nbins, vmin, vmax
92
for
(
G4int
k
=0;
k
<
kMaxHisto
;
k
++) {
93
G4int
ih = analysisManager->
CreateH1
(
id
[
k
], title[k], nbins, vmin, vmax);
94
analysisManager->
SetH1Activation
(ih,
false
);
95
}
96
97
}
98
99
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
100
101
void
RunAction::BeginOfRunAction
(
const
G4Run
* aRun)
102
{
103
G4cout
<<
"### Run "
<< aRun->
GetRunID
() <<
" start."
<<
G4endl
;
104
105
// save Rndm status
106
if
(
fRndmFreq
> 0) {
107
CLHEP::HepRandom::showEngineStatus
();
108
CLHEP::HepRandom::saveEngineStatus
(
"beginOfRun.rndm"
);
109
}
110
111
// histograms
112
G4AnalysisManager
* analysisManager =
G4AnalysisManager::Instance
();
113
if
( analysisManager->
IsActive
() ) {
114
analysisManager->
OpenFile
();
115
}
116
}
117
118
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
119
120
void
RunAction::EndOfRunAction
(
const
G4Run
* run)
121
{
122
// print run statisctics
123
G4int
nofEvents = run->
GetNumberOfEvent
();
124
G4cout
<<
" ================== run summary ====================="
<<
G4endl
;
125
G4cout
<<
" End of Run TotNbofEvents = "
<< nofEvents <<
G4endl
;
126
127
G4AnalysisManager
* analysisManager =
G4AnalysisManager::Instance
();
128
if
( analysisManager->
GetH1
(1) ) {
129
G4cout
<<
" Mean energy deposit in absorber = "
130
<< analysisManager->
GetH1
(1)->mean()/
MeV
<<
" +-"
131
<< analysisManager->
GetH1
(1)->rms()/
MeV
<<
" MeV "
<<
G4endl
;
132
}
133
if
( analysisManager->
GetH1
(2) ) {
134
G4cout
<<
" Total number of XTR gammas = "
135
<< analysisManager->
GetH1
(2)->entries() <<
G4endl
;
136
}
137
if
( analysisManager->
GetH1
(3) ) {
138
G4cout
<<
" Total number of all gammas = "
139
<< analysisManager->
GetH1
(3)->entries() <<
G4endl
;
140
}
141
142
// save Rndm status
143
if
(
fRndmFreq
== 1) {
144
CLHEP::HepRandom::showEngineStatus
();
145
CLHEP::HepRandom::saveEngineStatus
(
"endOfRun.rndm"
);
146
}
147
148
//save histograms
149
if
( analysisManager->
IsActive
() ) {
150
analysisManager->
Write
();
151
analysisManager->
CloseFile
();
152
}
153
}
154
155
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
geant4
tree
geant4-10.6-release
examples
extended
electromagnetic
TestEm10
src
RunAction.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:52
using
1.8.2 with
ECCE GitHub integration