Geant4
10.03.p03
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
FCALTestbeamSetupSD.cc
Go to the documentation of this file.
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: FCALTestbeamSetupSD.cc 89572 2015-04-20 07:18:15Z gcosmo $
27
//
28
//
29
30
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
31
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
32
33
#include <iostream>
34
35
#include "
FCALTestbeamSetupSD.hh
"
36
37
#include "
FCALCalorHit.hh
"
38
39
#include "
FCALTestbeamSetup.hh
"
40
41
#include "
G4SystemOfUnits.hh
"
42
#include "
G4VPhysicalVolume.hh
"
43
#include "
G4Step.hh
"
44
#include "
G4Track.hh
"
45
#include "
G4VTouchable.hh
"
46
#include "
G4TouchableHistory.hh
"
47
#include "
G4SDManager.hh
"
48
49
#include "
G4ios.hh
"
50
51
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
52
53
FCALTestbeamSetupSD::FCALTestbeamSetupSD
(
G4String
name
) :
G4VSensitiveDetector
(name)
54
{}
55
56
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57
58
FCALTestbeamSetupSD::~FCALTestbeamSetupSD
()
59
{}
60
61
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
62
63
void
FCALTestbeamSetupSD::Initialize
(
G4HCofThisEvent
*)
64
{
65
EBeamS1
=
EBeamS2
=
EBeamS3
= 0.;
66
EHoleScint
=
EBeamHole
= 0.;
67
EBeamDead
= 0;
68
G4int
j;
69
for
(j =0 ; j<
NLENGTH
; j++) {
70
ETailVis
[j] = 0.;
71
ETailDep
[j] = 0.;
72
}
73
TailCatcherID
= 0;
74
}
75
76
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
77
78
G4bool
FCALTestbeamSetupSD::ProcessHits
(
G4Step
* aStep,
G4TouchableHistory
*)
79
{
80
81
G4double
edep = aStep->
GetTotalEnergyDeposit
();
82
if
(edep==0.)
return
true
;
83
84
G4TouchableHistory
* theTouchable
85
= (
G4TouchableHistory
*)(aStep->
GetPreStepPoint
()->
GetTouchable
());
86
G4VPhysicalVolume
* physVol = theTouchable->
GetVolume
();
87
88
G4String
name
= physVol->
GetName
();
89
TailCatcherID
= physVol->
GetCopyNo
();
90
91
if
(name ==
"ScintS1Physical"
) {
92
EBeamS1
=
EBeamS1
+ edep;}
93
94
else
if
(name ==
"ScintS2Physical"
) {
95
EBeamS2
=
EBeamS2
+ edep;}
96
97
else
if
(name ==
"ScintS3Physical"
) {
98
EBeamS3
=
EBeamS3
+ edep;}
99
100
else
if
(name ==
"HoleScintPhysical"
){
EHoleScint
+= edep;}
101
else
if
(name ==
"HoleCntrScintPhysical"
){
102
EBeamHole
=
EBeamHole
+ edep;}
103
104
else
if
(name ==
"MWPCPhysical"
) {
EBeamDead
+= edep;}
105
else
if
(name ==
"HoleCntrPbPhysical"
) {
EBeamDead
+= edep;}
106
else
if
(name ==
"HoleCntrAlPhysical"
) {
EBeamDead
+= edep;}
107
else
if
(name ==
"LeadWallPhysical"
) {
EBeamDead
+= edep;}
108
else
if
(name ==
"IronWallPhysical"
) {
EBeamDead
+= edep;}
109
else
if
(
TailCatcherID
>= 0 &&
TailCatcherID
<
NLENGTH
) {
110
if
(name ==
"BigScintPhysical"
) {
111
ETailVis
[
TailCatcherID
] += edep;
112
}
113
else
if
(name ==
"SmallScintPhysical"
) {
114
if
(
TailCatcherID
+3 <
NLENGTH
)
ETailVis
[
TailCatcherID
+ 3] += edep;
115
}
116
else
if
(name ==
"BigIronPhysical"
) {
117
ETailDep
[
TailCatcherID
] += edep;
118
}
119
else
if
(name ==
"SmallIronPhysical"
) {
120
if
(
TailCatcherID
+2 <
NLENGTH
)
ETailDep
[
TailCatcherID
+2] += edep;
121
}
122
}
123
124
return
true
;
125
}
126
127
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
128
129
void
FCALTestbeamSetupSD::EndOfEvent
(
G4HCofThisEvent
*)
130
{
131
G4cout
<<
" Visisble Energy in S1 , S2 , S3 in (MeV)"
<<
G4endl
;
132
G4cout
<<
EBeamS1
/
MeV
<<
" "
<<
EBeamS2
/
MeV
<<
" "
<<
EBeamS3
/
MeV
<<
" "
<<
G4endl
;
133
134
G4cout
<<
" Visible Energy in Hole Counter (MeV) "
<<
G4endl
;
135
G4cout
<<
EHoleScint
/
MeV
<<
" "
<<
EBeamHole
/
MeV
<<
G4endl
;
136
137
G4cout
<<
" Visible Energy in Upstream Dead Materials "
<<
G4endl
;
138
G4cout
<<
EBeamDead
/
MeV
<<
G4endl
;
139
140
G4cout
<<
" Visible Energy in Tail Catcher Scintillator"
<<
G4endl
;
141
G4int
j;
142
for
(j=1; j<8 ; j++) {
G4cout
<<
ETailVis
[j]/
MeV
<<
" "
;};
G4cout
<<
G4endl
;
143
144
G4cout
<<
" Visible Energy in Tail Catcher Absorber"
<<
G4endl
;
145
for
(j=1; j<7 ; j++) {
G4cout
<<
ETailDep
[j]/
MeV
<<
" "
;};
G4cout
<<
G4endl
;
146
147
}
148
149
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
150
151
void
FCALTestbeamSetupSD::clear
()
152
{}
153
154
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
155
156
void
FCALTestbeamSetupSD::DrawAll
()
157
{}
158
159
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
160
161
void
FCALTestbeamSetupSD::PrintAll
()
162
{}
163
164
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
165
FCALTestbeamSetupSD::ProcessHits
G4bool ProcessHits(G4Step *, G4TouchableHistory *)
Definition:
FCALTestbeamSetupSD.cc:78
FCALTestbeamSetupSD::EHoleScint
G4double EHoleScint
Definition:
FCALTestbeamSetupSD.hh:74
name
const XML_Char * name
Definition:
expat.h:151
G4TouchableHistory.hh
FCALTestbeamSetupSD::clear
void clear()
Definition:
FCALTestbeamSetupSD.cc:151
G4TouchableHistory::GetVolume
G4VPhysicalVolume * GetVolume(G4int depth=0) const
G4Track.hh
G4VPhysicalVolume
Definition:
G4VPhysicalVolume.hh:82
FCALTestbeamSetupSD::EBeamS3
G4double EBeamS3
Definition:
FCALTestbeamSetupSD.hh:73
G4VSensitiveDetector
Definition:
G4VSensitiveDetector.hh:50
G4Step.hh
G4StepPoint::GetTouchable
const G4VTouchable * GetTouchable() const
FCALTestbeamSetupSD::PrintAll
void PrintAll()
Definition:
FCALTestbeamSetupSD.cc:161
G4int
int G4int
Definition:
G4Types.hh:78
FCALTestbeamSetupSD::ETailDep
G4double ETailDep[NLENGTH]
Definition:
FCALTestbeamSetupSD.hh:77
G4Step::GetPreStepPoint
G4StepPoint * GetPreStepPoint() const
FCALTestbeamSetupSD.hh
FCALTestbeamSetupSD::DrawAll
void DrawAll()
Definition:
FCALTestbeamSetupSD.cc:156
G4cout
G4GLOB_DLL std::ostream G4cout
G4VPhysicalVolume::GetName
const G4String & GetName() const
G4SDManager.hh
FCALTestbeamSetupSD::EBeamS1
G4double EBeamS1
Definition:
FCALTestbeamSetupSD.hh:73
G4bool
bool G4bool
Definition:
G4Types.hh:79
FCALCalorHit.hh
FCALTestbeamSetupSD::FCALTestbeamSetupSD
FCALTestbeamSetupSD(G4String)
Definition:
FCALTestbeamSetupSD.cc:53
G4Step
Definition:
G4Step.hh:76
G4Step::GetTotalEnergyDeposit
G4double GetTotalEnergyDeposit() const
FCALTestbeamSetupSD::EndOfEvent
void EndOfEvent(G4HCofThisEvent *)
Definition:
FCALTestbeamSetupSD.cc:129
G4TouchableHistory
Definition:
G4TouchableHistory.hh:53
FCALTestbeamSetupSD::~FCALTestbeamSetupSD
~FCALTestbeamSetupSD()
Definition:
FCALTestbeamSetupSD.cc:58
FCALTestbeamSetupSD::Initialize
void Initialize(G4HCofThisEvent *)
Definition:
FCALTestbeamSetupSD.cc:63
G4ios.hh
G4VPhysicalVolume::GetCopyNo
virtual G4int GetCopyNo() const =0
G4HCofThisEvent
Definition:
G4HCofThisEvent.hh:64
FCALTestbeamSetup.hh
FCALTestbeamSetupSD::ETailVis
G4double ETailVis[NLENGTH]
Definition:
FCALTestbeamSetupSD.hh:77
G4VTouchable.hh
G4endl
#define G4endl
Definition:
G4ios.hh:61
MeV
static constexpr double MeV
Definition:
G4SIunits.hh:214
FCALTestbeamSetupSD::TailCatcherID
G4int TailCatcherID
Definition:
FCALTestbeamSetupSD.hh:76
G4VPhysicalVolume.hh
G4double
double G4double
Definition:
G4Types.hh:76
G4SystemOfUnits.hh
NLENGTH
static const G4int NLENGTH
Definition:
FCALTestbeamSetupSD.hh:40
FCALTestbeamSetupSD::EBeamDead
G4double EBeamDead
Definition:
FCALTestbeamSetupSD.hh:75
FCALTestbeamSetupSD::EBeamS2
G4double EBeamS2
Definition:
FCALTestbeamSetupSD.hh:73
FCALTestbeamSetupSD::EBeamHole
G4double EBeamHole
Definition:
FCALTestbeamSetupSD.hh:74
G4String
Definition:
G4String.hh:45
source
geant4.10.03.p03
examples
advanced
lAr_calorimeter
src
FCALTestbeamSetupSD.cc
Generated on Tue Nov 28 2017 21:43:46 for Geant4 by
1.8.5