| 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 #include "SiPMHit.hh"
27
28 #include <cstdlib>
29
30 //....oooOO0OOooo........oooOO0OOooo........oo
31
32 SiPMHit::SiPMHit(G4String aVolumeName, G4int a
33 G4int aCopyNumCell)
34 : fVolumeName(aVolumeName), fCopyNumCell(a
35 fCopyNumSensor(aCopyNumSensor) {}
36
37 //....oooOO0OOooo........oooOO0OOooo........oo
38
39 void SiPMHit::Digitise(const G4double aTimeWin
40 const G4double aToaThre
41
42 // process energy deposits
43 if (fEdep.empty()) {
44 fIsValidHit = false;
45 return;
46 }
47
48 std::sort(fEdep.begin(), fEdep.end(),
49 [](const std::pair<G4double, G4dou
50 const std::pair<G4double, G4dou
51 return left.second < right.secon
52 });
53
54 G4double firstHitTime = fEdep[0].second;
55 fEdepDigi = 0;
56 for (size_t i = 0; i < fEdep.size(); ++i) {
57 if (aTimeWindow < 0 || fEdep[i].second < f
58 fEdepDigi += fEdep[i].first;
59 #ifdef DEBUG
60 else
61 std::cout << "Rejecting hit (" << fEdep[
62 << fEdep[i].second << "(start:
63 << std::endl;
64 #endif
65 if ((fTimeOfArrival == -1) && (fEdepDigi >
66 fTimeOfArrival = fEdep[i].second;
67 }
68 fIsValidHit = (fEdepDigi > 0);
69
70 // non ionizing part, does not contribute to
71 if (fEdepNonIonizing.empty())
72 return;
73
74 std::sort(fEdepNonIonizing.begin(), fEdepNon
75 [](const std::pair<G4double, G4dou
76 const std::pair<G4double, G4dou
77 return left.second < right.secon
78 });
79
80 fEdepNonIonizingDigi = 0;
81 for (size_t i = 0; i < fEdepNonIonizing.size
82 if (aTimeWindow == -1 ||
83 fEdepNonIonizing[i].second < firstHitT
84 fEdepNonIonizingDigi += fEdepNonIonizing
85 }
86 }
87