| 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 // 26 //
27 /// \file optical/LXe/src/LXePMTSD.cc 27 /// \file optical/LXe/src/LXePMTSD.cc
28 /// \brief Implementation of the LXePMTSD clas 28 /// \brief Implementation of the LXePMTSD class
29 // 29 //
30 // 30 //
31 #include "LXePMTSD.hh" 31 #include "LXePMTSD.hh"
32 32
33 #include "LXeDetectorConstruction.hh" 33 #include "LXeDetectorConstruction.hh"
34 #include "LXePMTHit.hh" 34 #include "LXePMTHit.hh"
35 #include "LXeUserTrackInformation.hh" 35 #include "LXeUserTrackInformation.hh"
36 36
>> 37 #include "G4ios.hh"
37 #include "G4LogicalVolume.hh" 38 #include "G4LogicalVolume.hh"
38 #include "G4ParticleDefinition.hh" 39 #include "G4ParticleDefinition.hh"
39 #include "G4ParticleTypes.hh" 40 #include "G4ParticleTypes.hh"
40 #include "G4SDManager.hh" 41 #include "G4SDManager.hh"
41 #include "G4Step.hh" 42 #include "G4Step.hh"
42 #include "G4TouchableHistory.hh" 43 #include "G4TouchableHistory.hh"
43 #include "G4Track.hh" 44 #include "G4Track.hh"
44 #include "G4VPhysicalVolume.hh" 45 #include "G4VPhysicalVolume.hh"
45 #include "G4VTouchable.hh" 46 #include "G4VTouchable.hh"
46 #include "G4ios.hh" <<
47 47
48 //....oooOO0OOooo........oooOO0OOooo........oo 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 49
50 LXePMTSD::LXePMTSD(G4String name) : G4VSensiti << 50 LXePMTSD::LXePMTSD(G4String name)
>> 51 : G4VSensitiveDetector(name)
>> 52 , fPMTHitCollection(nullptr)
>> 53 , fPMTPositionsX(nullptr)
>> 54 , fPMTPositionsY(nullptr)
>> 55 , fPMTPositionsZ(nullptr)
>> 56 , fHitCID(-1)
51 { 57 {
52 collectionName.insert("pmtHitCollection"); 58 collectionName.insert("pmtHitCollection");
53 } 59 }
54 60
55 //....oooOO0OOooo........oooOO0OOooo........oo 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
56 62
57 LXePMTSD::~LXePMTSD() 63 LXePMTSD::~LXePMTSD()
58 { 64 {
59 delete fPMTPositionsX; 65 delete fPMTPositionsX;
60 delete fPMTPositionsY; 66 delete fPMTPositionsY;
61 delete fPMTPositionsZ; 67 delete fPMTPositionsZ;
62 } 68 }
63 69
64 //....oooOO0OOooo........oooOO0OOooo........oo 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65 71
66 void LXePMTSD::SetPmtPositions(const std::vect 72 void LXePMTSD::SetPmtPositions(const std::vector<G4ThreeVector>& positions)
67 { 73 {
68 for (size_t i = 0; i < positions.size(); ++i << 74 for(size_t i = 0; i < positions.size(); ++i)
69 if (fPMTPositionsX) fPMTPositionsX->push_b << 75 {
70 if (fPMTPositionsY) fPMTPositionsY->push_b << 76 if(fPMTPositionsX)
71 if (fPMTPositionsZ) fPMTPositionsZ->push_b << 77 fPMTPositionsX->push_back(positions[i].x());
>> 78 if(fPMTPositionsY)
>> 79 fPMTPositionsY->push_back(positions[i].y());
>> 80 if(fPMTPositionsZ)
>> 81 fPMTPositionsZ->push_back(positions[i].z());
72 } 82 }
73 } 83 }
74 84
75 //....oooOO0OOooo........oooOO0OOooo........oo 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
76 86
77 void LXePMTSD::Initialize(G4HCofThisEvent* hit 87 void LXePMTSD::Initialize(G4HCofThisEvent* hitsCE)
78 { 88 {
79 fPMTHitCollection = new LXePMTHitsCollection << 89 fPMTHitCollection =
>> 90 new LXePMTHitsCollection(SensitiveDetectorName, collectionName[0]);
80 91
81 if (fHitCID < 0) { << 92 if(fHitCID < 0)
>> 93 {
82 fHitCID = G4SDManager::GetSDMpointer()->Ge 94 fHitCID = G4SDManager::GetSDMpointer()->GetCollectionID(fPMTHitCollection);
83 } 95 }
84 hitsCE->AddHitsCollection(fHitCID, fPMTHitCo 96 hitsCE->AddHitsCollection(fHitCID, fPMTHitCollection);
85 } 97 }
86 98
87 //....oooOO0OOooo........oooOO0OOooo........oo 99 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
88 100
89 G4bool LXePMTSD::ProcessHits(G4Step*, G4Toucha << 101 G4bool LXePMTSD::ProcessHits(G4Step*, G4TouchableHistory*) { return false; }
90 { <<
91 return false; <<
92 } <<
93 102
94 //....oooOO0OOooo........oooOO0OOooo........oo 103 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
95 104
96 // Generates a hit and uses the postStepPoint' 105 // Generates a hit and uses the postStepPoint's mother volume replica number
97 // PostStepPoint because the hit is generated 106 // PostStepPoint because the hit is generated manually when the photon is
98 // absorbed by the photocathode 107 // absorbed by the photocathode
99 108
100 G4bool LXePMTSD::ProcessHits_boundary(const G4 109 G4bool LXePMTSD::ProcessHits_boundary(const G4Step* aStep, G4TouchableHistory*)
101 { 110 {
102 // need to know if this is an optical photon 111 // need to know if this is an optical photon
103 if (aStep->GetTrack()->GetDefinition() != G4 << 112 if(aStep->GetTrack()->GetDefinition() !=
>> 113 G4OpticalPhoton::OpticalPhotonDefinition())
104 return false; 114 return false;
105 115
106 // User replica number 1 since photocathode 116 // User replica number 1 since photocathode is a daughter volume
107 // to the pmt which was replicated 117 // to the pmt which was replicated
108 G4int pmtNumber = aStep->GetPostStepPoint()- << 118 G4int pmtNumber =
109 G4VPhysicalVolume* physVol = aStep->GetPostS << 119 aStep->GetPostStepPoint()->GetTouchable()->GetReplicaNumber(1);
>> 120 G4VPhysicalVolume* physVol =
>> 121 aStep->GetPostStepPoint()->GetTouchable()->GetVolume(1);
110 122
111 // Find the correct hit collection 123 // Find the correct hit collection
112 size_t n = fPMTHitCollection->entries(); << 124 size_t n = fPMTHitCollection->entries();
113 LXePMTHit* hit = nullptr; 125 LXePMTHit* hit = nullptr;
114 for (size_t i = 0; i < n; ++i) { << 126 for(size_t i = 0; i < n; ++i)
115 if ((*fPMTHitCollection)[i]->GetPMTNumber( << 127 {
>> 128 if((*fPMTHitCollection)[i]->GetPMTNumber() == pmtNumber)
>> 129 {
116 hit = (*fPMTHitCollection)[i]; 130 hit = (*fPMTHitCollection)[i];
117 break; 131 break;
118 } 132 }
119 } 133 }
120 134
121 if (hit == nullptr) { // this pmt wasn't pr << 135 if(hit == nullptr)
>> 136 { // this pmt wasn't previously hit in this event
122 hit = new LXePMTHit(); // so create new h 137 hit = new LXePMTHit(); // so create new hit
123 hit->SetPMTNumber(pmtNumber); 138 hit->SetPMTNumber(pmtNumber);
124 hit->SetPMTPhysVol(physVol); 139 hit->SetPMTPhysVol(physVol);
125 fPMTHitCollection->insert(hit); 140 fPMTHitCollection->insert(hit);
126 hit->SetPMTPos((*fPMTPositionsX)[pmtNumber 141 hit->SetPMTPos((*fPMTPositionsX)[pmtNumber], (*fPMTPositionsY)[pmtNumber],
127 (*fPMTPositionsZ)[pmtNumber 142 (*fPMTPositionsZ)[pmtNumber]);
128 } 143 }
129 144
130 hit->IncPhotonCount(); // increment hit for 145 hit->IncPhotonCount(); // increment hit for the selected pmt
131 146
132 if (!LXeDetectorConstruction::GetSphereOn()) << 147 if(!LXeDetectorConstruction::GetSphereOn())
>> 148 {
133 hit->SetDrawit(true); 149 hit->SetDrawit(true);
134 // If the sphere is disabled then this hit 150 // If the sphere is disabled then this hit is automaticaly drawn
135 } 151 }
136 else { // sphere enabled << 152 else
137 auto trackInfo = (LXeUserTrackInformation* << 153 { // sphere enabled
138 if (trackInfo->GetTrackStatus() & hitSpher << 154 LXeUserTrackInformation* trackInfo =
>> 155 (LXeUserTrackInformation*) aStep->GetTrack()->GetUserInformation();
>> 156 if(trackInfo->GetTrackStatus() & hitSphere)
139 // only draw this hit if the photon has 157 // only draw this hit if the photon has hit the sphere first
140 hit->SetDrawit(true); 158 hit->SetDrawit(true);
141 } 159 }
142 160
143 return true; 161 return true;
144 } 162 }
145 163