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 // 27 // -------------------------------------------------------------- 28 // GEANT 4 - ULTRA experiment example 29 // -------------------------------------------------------------- 30 // 31 // Code developed by: 32 // B. Tome, M.C. Espirito-Santo, A. Trindade, P. Rodrigues 33 // 34 // **************************************************** 35 // * UltraPMTSD.cc 36 // **************************************************** 37 // 38 // Class used to define the Ultra photomultiplier as a sensitive detector. 39 // Hits in this sensitive detector are defined in the UltraOpticalHit class 40 // 41 #include "UltraPMTSD.hh" 42 43 #include "G4Material.hh" 44 #include "G4Step.hh" 45 #include "G4VTouchable.hh" 46 #include "G4TouchableHistory.hh" 47 #include "G4SDManager.hh" 48 #include "G4OpticalPhoton.hh" 49 #include "G4ParticleDefinition.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 52 53 UltraPMTSD::UltraPMTSD(G4String name):G4VSensitiveDetector(name) 54 { 55 56 collectionName.insert("OpticalHitsCollection"); 57 58 } 59 60 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 62 63 UltraPMTSD::~UltraPMTSD() 64 {;} 65 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 67 68 void UltraPMTSD::Initialize(G4HCofThisEvent* HCE) 69 { 70 71 72 static int HCID1 = -1; 73 74 75 // SensitiveDetectorName and collectionName are data members of G4VSensitiveDetector 76 77 OpticalHitsCollection = 78 new UltraOpticalHitsCollection(SensitiveDetectorName,collectionName[0]); 79 80 if(HCID1<0) 81 { HCID1 = GetCollectionID(0); } 82 HCE->AddHitsCollection(HCID1,OpticalHitsCollection); 83 84 } 85 86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 87 88 G4bool UltraPMTSD::ProcessHits(G4Step* aStep,G4TouchableHistory*) 89 { 90 91 // Get Material 92 93 auto thisVolume = aStep->GetTrack()->GetVolume()->GetName() ; 94 const auto particle = aStep->GetTrack()->GetDefinition(); 95 96 97 if (thisVolume != "PMT1" && thisVolume != "PMT2") 98 return false; 99 if (particle != G4OpticalPhoton::Definition() ) 100 return false; 101 102 if(particle == G4OpticalPhoton::Definition()) 103 aStep->GetTrack()->SetTrackStatus(fStopAndKill); 104 105 auto kineticEnergy = aStep->GetTrack()->GetKineticEnergy(); 106 auto HitPosition = aStep->GetPreStepPoint()->GetPosition() ; 107 108 auto OpticalHit = new UltraOpticalHit ; 109 OpticalHit->SetEnergy(kineticEnergy); 110 OpticalHit->SetPosition(HitPosition); 111 112 113 OpticalHitsCollection->insert(OpticalHit); 114 115 116 #ifdef ULTRA_VERBOSE 117 G4cout << "*******************************" << G4endl; 118 G4cout << " PMT HIT " << G4endl; 119 G4cout << " Volume: " << thisVolume << G4endl; 120 G4cout << " Photon energy (eV) : " << kineticEnergy/CLHEP::eV << G4endl; 121 G4cout << " POSITION (mm) : " 122 << HitPosition.x()/CLHEP::mm << " " << HitPosition.y()/CLHEP::mm << " " << HitPosition.z()/CLHEP::mm << G4endl; 123 G4cout << "*******************************" << G4endl; 124 #endif 125 126 127 return true; 128 } 129 130 131 void UltraPMTSD::EndOfEvent(G4HCofThisEvent* HCE) 132 { 133 static G4int HCID = -1; 134 if(HCID<0) 135 { 136 HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); 137 } 138 HCE->AddHitsCollection(HCID,OpticalHitsCollection); 139 140 } 141 142 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 143 144 void UltraPMTSD::clear() 145 {;} 146 147 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 148 149 void UltraPMTSD::DrawAll() 150 {;} 151 152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 153 154 void UltraPMTSD::PrintAll() 155 {;} 156 157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 158