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