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