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
RE05CalorimeterSD.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file RE05CalorimeterSD.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
//
26
//
29
//
30
31
#include "
RE05CalorimeterSD.hh
"
32
#include "
RE05CalorimeterHit.hh
"
33
#include "
G4VPhysicalVolume.hh
"
34
#include "
G4LogicalVolume.hh
"
35
#include "
G4Track.hh
"
36
#include "
G4Step.hh
"
37
#include "
G4ParticleDefinition.hh
"
38
#include "
G4VTouchable.hh
"
39
#include "
G4TouchableHistory.hh
"
40
#include "
G4SystemOfUnits.hh
"
41
#include "
G4ios.hh
"
42
43
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44
45
RE05CalorimeterSD::RE05CalorimeterSD
(
G4String
name
)
46
:
G4VSensitiveDetector
(name),
47
fNumberOfCellsInZ(20),fNumberOfCellsInPhi(48)
48
{
49
G4String
HCname;
50
collectionName
.
insert
(HCname=
"calCollection"
);
51
}
52
53
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54
55
RE05CalorimeterSD::~RE05CalorimeterSD
()
56
{}
57
58
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59
60
void
RE05CalorimeterSD::Initialize
(
G4HCofThisEvent
*)
61
{
62
fCalCollection
=
new
RE05CalorimeterHitsCollection
63
(
SensitiveDetectorName
,
collectionName
[0]);
64
for
(
G4int
j=0;j<
fNumberOfCellsInZ
;j++)
65
for
(
G4int
k
=0;
k
<
fNumberOfCellsInPhi
;
k
++)
66
{
67
fCellID
[j][
k
] = -1;
68
}
69
verboseLevel
= 0;
70
}
71
72
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73
74
G4bool
RE05CalorimeterSD::ProcessHits
(
G4Step
*aStep,
G4TouchableHistory
*)
75
{
76
//***** RE05CalorimeterSD has been migrated to Geant4 version 10 that does not
77
//***** support Readout Geometry in multi-threaded mode. Now RE05CalorimeterSD
78
//***** is assigned to a dedicated parallel world. The pointer "aStep" points to
79
//***** a G4Step object for the parallel world.
80
81
G4double
edep
= aStep->
GetTotalEnergyDeposit
();
82
if
(
verboseLevel
>1)
G4cout
<<
"Next step edep(MeV) = "
<< edep/
MeV
<<
G4endl
;
83
if
(edep==0.)
return
false
;
84
85
const
G4VTouchable
* touchable = aStep->
GetPreStepPoint
()->
GetTouchable
();
86
G4int
copyIDinZ = touchable->
GetReplicaNumber
(0);
87
G4int
copyIDinPhi = touchable->
GetReplicaNumber
(1);
88
89
if
(
fCellID
[copyIDinZ][copyIDinPhi]==-1)
90
{
91
RE05CalorimeterHit
* calHit
92
=
new
RE05CalorimeterHit
93
(touchable->
GetVolume
()->
GetLogicalVolume
(),copyIDinZ,copyIDinPhi);
94
calHit->
SetEdep
( edep );
95
G4AffineTransform
aTrans = touchable->
GetHistory
()->
GetTopTransform
();
96
aTrans.
Invert
();
97
calHit->
SetPos
(aTrans.
NetTranslation
());
98
calHit->
SetRot
(aTrans.
NetRotation
());
99
G4int
icell =
fCalCollection
->
insert
( calHit );
100
fCellID
[copyIDinZ][copyIDinPhi] = icell - 1;
101
if
(
verboseLevel
>0)
102
{
G4cout
<<
" New Calorimeter Hit on fCellID "
103
<< copyIDinZ <<
" "
<< copyIDinPhi <<
G4endl
; }
104
}
105
else
106
{
107
(*fCalCollection)[
fCellID
[copyIDinZ][copyIDinPhi]]->AddEdep(edep);
108
if
(
verboseLevel
>0)
109
{
G4cout
<<
" Energy added to fCellID "
110
<< copyIDinZ <<
" "
<< copyIDinPhi <<
G4endl
; }
111
}
112
113
return
true
;
114
}
115
116
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
117
118
void
RE05CalorimeterSD::EndOfEvent
(
G4HCofThisEvent
*HCE)
119
{
120
static
G4int
HCID = -1;
121
if
(HCID<0)
122
{ HCID =
GetCollectionID
(0); }
123
HCE->
AddHitsCollection
( HCID,
fCalCollection
);
124
}
125
126
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127
128
void
RE05CalorimeterSD::clear
()
129
{}
130
131
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
132
133
void
RE05CalorimeterSD::DrawAll
()
134
{}
135
136
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
137
138
void
RE05CalorimeterSD::PrintAll
()
139
{}
140
141
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
geant4
tree
geant4-10.6-release
examples
extended
runAndEvent
RE05
src
RE05CalorimeterSD.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:10
using
1.8.2 with
ECCE GitHub integration