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