Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // 27 // ------------------------------------------- 24 // -------------------------------------------------------------- 28 // GEANT 4 - Underground Dark Matter Detecto 25 // GEANT 4 - Underground Dark Matter Detector Advanced Example 29 // 26 // 30 // For information related to this code c 27 // For information related to this code contact: Alex Howard 31 // e-mail: alexander.howard@cern.ch << 28 // e-mail: a.s.howard@ic.ac.uk 32 // ------------------------------------------- 29 // -------------------------------------------------------------- 33 // Comments 30 // Comments 34 // 31 // 35 // Underground Advanced 32 // Underground Advanced 36 // by A. Howard and H. Araujo 33 // by A. Howard and H. Araujo 37 // (27th November 2001) 34 // (27th November 2001) 38 // 35 // 39 // ScintSD (scintillator sensitive detector de 36 // ScintSD (scintillator sensitive detector definition) program 40 // ------------------------------------------- 37 // -------------------------------------------------------------- 41 38 42 #include "DMXScintSD.hh" 39 #include "DMXScintSD.hh" 43 40 44 #include "DMXScintHit.hh" 41 #include "DMXScintHit.hh" 45 #include "DMXDetectorConstruction.hh" 42 #include "DMXDetectorConstruction.hh" 46 43 47 #include "G4VPhysicalVolume.hh" 44 #include "G4VPhysicalVolume.hh" 48 #include "G4HCofThisEvent.hh" 45 #include "G4HCofThisEvent.hh" 49 #include "G4Step.hh" 46 #include "G4Step.hh" 50 #include "G4SDManager.hh" 47 #include "G4SDManager.hh" 51 #include "G4ParticleDefinition.hh" 48 #include "G4ParticleDefinition.hh" 52 #include "G4ParticleTypes.hh" 49 #include "G4ParticleTypes.hh" 53 #include "G4Ions.hh" 50 #include "G4Ions.hh" 54 #include "G4ios.hh" 51 #include "G4ios.hh" 55 52 56 #include "G4OpticalPhoton.hh" << 57 << 58 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 59 54 60 DMXScintSD::DMXScintSD(G4String name) << 55 DMXScintSD::DMXScintSD(G4String name,DMXDetectorConstruction* DMXSD) 61 :G4VSensitiveDetector(name) << 56 :G4VSensitiveDetector(name),DMXDetector(DMXSD) 62 { 57 { 63 G4String HCname="scintillatorCollection"; 58 G4String HCname="scintillatorCollection"; 64 collectionName.insert(HCname); 59 collectionName.insert(HCname); 65 } 60 } 66 61 67 //....oooOO0OOooo........oooOO0OOooo........oo 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 68 63 69 DMXScintSD::~DMXScintSD(){ } 64 DMXScintSD::~DMXScintSD(){ } 70 65 71 //....oooOO0OOooo........oooOO0OOooo........oo 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 72 67 73 void DMXScintSD::Initialize(G4HCofThisEvent*) 68 void DMXScintSD::Initialize(G4HCofThisEvent*) 74 { 69 { 75 scintillatorCollection = new DMXScintHitsCol 70 scintillatorCollection = new DMXScintHitsCollection 76 (SensitiveDetectorName,collectionName[0]); 71 (SensitiveDetectorName,collectionName[0]); 77 72 78 HitID = -1; 73 HitID = -1; 79 } 74 } 80 75 81 //....oooOO0OOooo........oooOO0OOooo........oo 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 82 77 83 G4bool DMXScintSD::ProcessHits(G4Step* aStep, 78 G4bool DMXScintSD::ProcessHits(G4Step* aStep, G4TouchableHistory*) 84 { 79 { 85 << 86 //need to know if this is an optical photon << 87 if(aStep->GetTrack()->GetDefinition() << 88 == G4OpticalPhoton::OpticalPhotonDefiniti << 89 << 90 << 91 G4double edep = aStep->GetTotalEnergyDeposit 80 G4double edep = aStep->GetTotalEnergyDeposit(); 92 G4ParticleDefinition* particleType = aStep-> 81 G4ParticleDefinition* particleType = aStep->GetTrack()->GetDefinition(); 93 G4String particleName = particleType->GetPar 82 G4String particleName = particleType->GetParticleName(); 94 83 95 G4double stepl = 0.; 84 G4double stepl = 0.; 96 if (particleType->GetPDGCharge() != 0.) 85 if (particleType->GetPDGCharge() != 0.) 97 stepl = aStep->GetStepLength(); 86 stepl = aStep->GetStepLength(); 98 87 99 if ((edep==0.)&&(stepl==0.)) return false; 88 if ((edep==0.)&&(stepl==0.)) return false; 100 89 101 90 102 // fill in hit 91 // fill in hit 103 DMXScintHit* newHit = new DMXScintHit(); 92 DMXScintHit* newHit = new DMXScintHit(); 104 newHit->SetEdep(edep); 93 newHit->SetEdep(edep); 105 newHit->SetPos(aStep->GetPostStepPoint()->Ge 94 newHit->SetPos(aStep->GetPostStepPoint()->GetPosition()); 106 newHit->SetTime(aStep->GetPreStepPoint()->Ge 95 newHit->SetTime(aStep->GetPreStepPoint()->GetGlobalTime()); 107 newHit->SetParticle(particleName); 96 newHit->SetParticle(particleName); 108 newHit->SetParticleEnergy(aStep->GetPreStepP 97 newHit->SetParticleEnergy(aStep->GetPreStepPoint()->GetKineticEnergy() ); 109 98 110 HitID = scintillatorCollection->insert(newHi 99 HitID = scintillatorCollection->insert(newHit); 111 100 112 return true; 101 return true; 113 } 102 } 114 103 115 104 116 //....oooOO0OOooo........oooOO0OOooo........oo 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 117 106 118 void DMXScintSD::EndOfEvent(G4HCofThisEvent* H 107 void DMXScintSD::EndOfEvent(G4HCofThisEvent* HCE) 119 { 108 { 120 109 121 G4String HCname = collectionName[0]; 110 G4String HCname = collectionName[0]; 122 static G4int HCID = -1; 111 static G4int HCID = -1; 123 if(HCID<0) 112 if(HCID<0) 124 HCID = G4SDManager::GetSDMpointer()->GetCo 113 HCID = G4SDManager::GetSDMpointer()->GetCollectionID(HCname); 125 HCE->AddHitsCollection(HCID,scintillatorColl 114 HCE->AddHitsCollection(HCID,scintillatorCollection); 126 115 127 G4int nHits = scintillatorCollection->entrie 116 G4int nHits = scintillatorCollection->entries(); 128 if (verboseLevel>=1) 117 if (verboseLevel>=1) 129 G4cout << " LXe collection: " << nHit 118 G4cout << " LXe collection: " << nHits << " hits" << G4endl; 130 if (verboseLevel>=2) 119 if (verboseLevel>=2) 131 scintillatorCollection->PrintAllHits(); 120 scintillatorCollection->PrintAllHits(); 132 121 133 } 122 } 134 123 135 //....oooOO0OOooo........oooOO0OOooo........oo 124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 136 125 137 void DMXScintSD::clear() 126 void DMXScintSD::clear() 138 {} 127 {} 139 128 140 //....oooOO0OOooo........oooOO0OOooo........oo 129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 141 130 142 void DMXScintSD::DrawAll() 131 void DMXScintSD::DrawAll() 143 {} 132 {} 144 133 145 //....oooOO0OOooo........oooOO0OOooo........oo 134 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 146 135 147 void DMXScintSD::PrintAll() 136 void DMXScintSD::PrintAll() 148 {} 137 {} 149 138 150 //....oooOO0OOooo........oooOO0OOooo........oo 139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 151 140 152 141