| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file exoticphysics/phonon/src/XAluminumEl
27 /// \brief Implementation of the XAluminumElec
28 //
29 //
30 #include "XAluminumElectrodeSensitivity.hh"
31
32 #include "XAluminumElectrodeHit.hh"
33
34 #include "G4AutoLock.hh"
35 #include "G4HCofThisEvent.hh"
36 #include "G4Navigator.hh"
37 #include "G4SDManager.hh"
38 #include "G4Step.hh"
39 #include "G4SystemOfUnits.hh"
40 #include "G4Threading.hh"
41 #include "G4TouchableHistory.hh"
42 #include "G4Track.hh"
43
44 #include <fstream>
45
46 std::fstream* XAluminumElectrodeSensitivity::f
47 std::fstream* XAluminumElectrodeSensitivity::f
48
49 G4Mutex theMutex = G4MUTEX_INITIALIZER; // Ju
50
51 //....oooOO0OOooo........oooOO0OOooo........oo
52
53 XAluminumElectrodeSensitivity::XAluminumElectr
54 : G4VSensitiveDetector(name)
55 {
56 collectionName.insert("XAluminumElectrodeHit
57 fHCID = -1;
58
59 G4AutoLock lockIt(&theMutex); // Only one t
60 fWriter = new std::fstream("caustic.ssv", st
61 if (!fWriter->is_open()) {
62 G4cerr << "XAluminumElectrodeSensitivity::
63 << "\n\tFailed to open caustic.ssv
64 << "\n\tCreating caustic.ssv" << G4
65 fWriter->open("caustic.ssv");
66 }
67
68 fWriter2 = new std::fstream("timing.ssv", st
69 if (!fWriter2->is_open()) {
70 G4cerr << "XAluminumElectrodeSensitivity::
71 << "\n\tFailed to open timing.ssv f
72 << "\n\tCreating timing.ssv." << G4
73 fWriter2->open("timing.ssv");
74 }
75
76 if (!(fWriter->is_open() && fWriter2->is_ope
77 G4cerr << "XAluminumElectrodeSensitivity::
78 << "\nERROR: COULD NOT CREATE OUTPU
79 }
80 }
81
82 //....oooOO0OOooo........oooOO0OOooo........oo
83
84 XAluminumElectrodeSensitivity::~XAluminumElect
85 {
86 G4AutoLock lockIt(&theMutex); // Only one t
87
88 if (fWriter) {
89 fWriter->close();
90 delete fWriter;
91 fWriter = 0;
92 }
93
94 if (fWriter2) {
95 fWriter2->close();
96 delete fWriter2;
97 fWriter2 = 0;
98 }
99 }
100
101 //....oooOO0OOooo........oooOO0OOooo........oo
102
103 XAluminumElectrodeHitsCollection* XAluminumEle
104 {
105 return fHitsCollection;
106 }
107
108 //....oooOO0OOooo........oooOO0OOooo........oo
109
110 void XAluminumElectrodeSensitivity::Initialize
111 {
112 fHitsCollection = new XAluminumElectrodeHits
113 if (fHCID < 0) {
114 fHCID = G4SDManager::GetSDMpointer()->GetC
115 }
116 HCE->AddHitsCollection(fHCID, fHitsCollectio
117 }
118
119 //....oooOO0OOooo........oooOO0OOooo........oo
120
121 G4bool XAluminumElectrodeSensitivity::ProcessH
122 {
123 // if(aStep->GetTrack()->GetDefinition()!=Ph
124 G4double edp = aStep->GetNonIonizingEnergyDe
125 if (edp == 0.) return true;
126
127 G4StepPoint* preStepPoint = aStep->GetPreSte
128 G4StepPoint* postStepPoint = aStep->GetPostS
129 G4TouchableHistory* theTouchable = (G4Toucha
130 G4ThreeVector fWorldPos = postStepPoint->Get
131 G4ThreeVector fLocalPos = theTouchable->GetH
132
133 XAluminumElectrodeHit* aHit = new XAluminumE
134 aHit->fTime = postStepPoint->GetGlobalTime()
135 aHit->fEdep = edp;
136 aHit->fWorldPos = fWorldPos;
137 aHit->fLocalPos = fLocalPos;
138
139 fHitsCollection->insert(aHit);
140
141 return true;
142 }
143
144 //....oooOO0OOooo........oooOO0OOooo........oo
145
146 void XAluminumElectrodeSensitivity::EndOfEvent
147 {
148 if (!fHitsCollection || fHitsCollection->Get
149
150 for (size_t i = 0; i < fHitsCollection->GetS
151 WriteHitInfo(dynamic_cast<XAluminumElectro
152 }
153 }
154
155 //....oooOO0OOooo........oooOO0OOooo........oo
156
157 void XAluminumElectrodeSensitivity::WriteHitIn
158 {
159 if (!aHit) return;
160
161 G4AutoLock lockIt(&theMutex); // Only one e
162
163 *fWriter << aHit->fWorldPos.getX() / mm << "
164 << aHit->fWorldPos.getZ() / mm << "
165
166 *fWriter2 << aHit->fTime / ns << " " << aHit
167 }
168