Geant4
10.02
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
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
G4TouchableHistory.hh
MeV
static const double MeV
Definition:
G4SIunits.hh:211
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
name
G4String name
Definition:
TRTMaterials.hh:40
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
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
geant4.10.02
examples
advanced
lAr_calorimeter
src
FCALTestbeamSetupSD.cc
Generated on Thu Dec 31 2015 10:43:28 for Geant4 by
1.8.8