| Geant4 Cross Reference |
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 /// \file exoticphysics/phonon/src/XAluminumElectrodeSensitivity.cc
27 /// \brief Implementation of the XAluminumElectrodeSensitivity class
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::fWriter = 0;
47 std::fstream* XAluminumElectrodeSensitivity::fWriter2 = 0;
48
49 G4Mutex theMutex = G4MUTEX_INITIALIZER; // Just need one
50
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
52
53 XAluminumElectrodeSensitivity::XAluminumElectrodeSensitivity(const G4String& name)
54 : G4VSensitiveDetector(name)
55 {
56 collectionName.insert("XAluminumElectrodeHit");
57 fHCID = -1;
58
59 G4AutoLock lockIt(&theMutex); // Only one thread opens files!
60 fWriter = new std::fstream("caustic.ssv", std::fstream::out | std::fstream::ate);
61 if (!fWriter->is_open()) {
62 G4cerr << "XAluminumElectrodeSensitivity::Constructor:"
63 << "\n\tFailed to open caustic.ssv for appending data."
64 << "\n\tCreating caustic.ssv" << G4endl;
65 fWriter->open("caustic.ssv");
66 }
67
68 fWriter2 = new std::fstream("timing.ssv", std::fstream::out | std::fstream::ate);
69 if (!fWriter2->is_open()) {
70 G4cerr << "XAluminumElectrodeSensitivity::Constructor: "
71 << "\n\tFailed to open timing.ssv for appending data."
72 << "\n\tCreating timing.ssv." << G4endl;
73 fWriter2->open("timing.ssv");
74 }
75
76 if (!(fWriter->is_open() && fWriter2->is_open())) {
77 G4cerr << "XAluminumElectrodeSensitivity::Constructor: "
78 << "\nERROR: COULD NOT CREATE OUTPUT FILES FOR WRITING" << G4endl;
79 }
80 }
81
82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
83
84 XAluminumElectrodeSensitivity::~XAluminumElectrodeSensitivity()
85 {
86 G4AutoLock lockIt(&theMutex); // Only one thread deletes!
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........oooOO0OOooo........oooOO0OOooo....
102
103 XAluminumElectrodeHitsCollection* XAluminumElectrodeSensitivity::GetHitsCollection()
104 {
105 return fHitsCollection;
106 }
107
108 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
109
110 void XAluminumElectrodeSensitivity::Initialize(G4HCofThisEvent* HCE)
111 {
112 fHitsCollection = new XAluminumElectrodeHitsCollection(SensitiveDetectorName, collectionName[0]);
113 if (fHCID < 0) {
114 fHCID = G4SDManager::GetSDMpointer()->GetCollectionID(fHitsCollection);
115 }
116 HCE->AddHitsCollection(fHCID, fHitsCollection);
117 }
118
119 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
120
121 G4bool XAluminumElectrodeSensitivity::ProcessHits(G4Step* aStep, G4TouchableHistory* /*ROhist*/)
122 {
123 // if(aStep->GetTrack()->GetDefinition()!=Phonon::PhononDefinition()) return true;
124 G4double edp = aStep->GetNonIonizingEnergyDeposit();
125 if (edp == 0.) return true;
126
127 G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
128 G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
129 G4TouchableHistory* theTouchable = (G4TouchableHistory*)(preStepPoint->GetTouchable());
130 G4ThreeVector fWorldPos = postStepPoint->GetPosition();
131 G4ThreeVector fLocalPos = theTouchable->GetHistory()->GetTopTransform().TransformPoint(fWorldPos);
132
133 XAluminumElectrodeHit* aHit = new XAluminumElectrodeHit();
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........oooOO0OOooo........oooOO0OOooo....
145
146 void XAluminumElectrodeSensitivity::EndOfEvent(G4HCofThisEvent* /*HCE*/)
147 {
148 if (!fHitsCollection || fHitsCollection->GetSize() == 0) return;
149
150 for (size_t i = 0; i < fHitsCollection->GetSize(); i++) {
151 WriteHitInfo(dynamic_cast<XAluminumElectrodeHit*>(fHitsCollection->GetHit(i)));
152 }
153 }
154
155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
156
157 void XAluminumElectrodeSensitivity::WriteHitInfo(const XAluminumElectrodeHit* aHit)
158 {
159 if (!aHit) return;
160
161 G4AutoLock lockIt(&theMutex); // Only one event can write at a time
162
163 *fWriter << aHit->fWorldPos.getX() / mm << "," << aHit->fWorldPos.getY() / mm << ","
164 << aHit->fWorldPos.getZ() / mm << "\n";
165
166 *fWriter2 << aHit->fTime / ns << " " << aHit->fEdep / eV << "\n";
167 }
168