| 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 class implementation file
29 // CERN Geneva Switzerland
30 //
31 //
32 // ------------ GammaRayTelAnticoincidenceSD ------
33 // by R.Giannitrapani, F.Longo & G.Santin (13 nov 2000)
34 //
35 // ************************************************************
36
37 #include "GammaRayTelAnticoincidenceSD.hh"
38 #include "GammaRayTelAnticoincidenceHit.hh"
39 #include "GammaRayTelDetectorConstruction.hh"
40
41 #include "G4ios.hh"
42 #include "G4RunManager.hh"
43 #include "G4SDManager.hh"
44 #include "G4Step.hh"
45 #include "G4SystemOfUnits.hh"
46 #include "G4TouchableHistory.hh"
47 #include "G4VPhysicalVolume.hh"
48 #include "G4VTouchable.hh"
49
50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
51
52 GammaRayTelAnticoincidenceSD::GammaRayTelAnticoincidenceSD(G4String name) : G4VSensitiveDetector(name) {
53 auto *runManager = G4RunManager::GetRunManager();
54 detector = (GammaRayTelDetectorConstruction*) (runManager->GetUserDetectorConstruction());
55
56 numberOfACDLateralTiles = detector->GetNbOfACDLateralTiles();
57 numberOfACDTopTiles = detector->GetNbOfACDTopTiles();
58
59 G4cout << numberOfACDLateralTiles << " LATERAL " << G4endl;
60 G4cout << numberOfACDTopTiles << " TOP " << G4endl;
61
62 hitLateralID = new G4int[numberOfACDLateralTiles];
63 hitTopID = new G4int[numberOfACDTopTiles];
64 collectionName.insert("AnticoincidenceCollection");
65 }
66
67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
68
69 GammaRayTelAnticoincidenceSD::~GammaRayTelAnticoincidenceSD() {
70 delete[] hitLateralID;
71 delete[] hitTopID;
72 }
73
74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
75
76 void GammaRayTelAnticoincidenceSD::Initialize(G4HCofThisEvent*) {
77 anticoincidenceCollection = new GammaRayTelAnticoincidenceHitsCollection(SensitiveDetectorName, collectionName[0]);
78 for (auto i = 0; i < numberOfACDLateralTiles; i++) {
79 hitLateralID[i] = -1;
80 }
81
82 for (auto j = 0; j < numberOfACDTopTiles; j++) {
83 hitTopID[j] = -1;
84 }
85 }
86
87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
88
89 auto GammaRayTelAnticoincidenceSD::ProcessHits(G4Step *step, G4TouchableHistory*) -> G4bool {
90 auto depositedEnergy = step->GetTotalEnergyDeposit();
91 if (depositedEnergy / keV == 0.) {
92 return false;
93 }
94
95 // This TouchableHistory is used to obtain the physical volume of the hit
96 auto *theTouchable = (G4TouchableHistory*) (step->GetPreStepPoint()->GetTouchable());
97
98 auto *acdTile = theTouchable->GetVolume();
99 auto acdTileNumber = acdTile->GetCopyNo();
100 auto acdTileName = acdTile->GetName();
101
102 // G4cout << acdTileName << " " << depositedEnergy / keV << G4endl;
103
104 if (acdTileName == "ACT") { // The hit is on a top ACD tile (ACDType: 0)
105 if (hitTopID[acdTileNumber] == -1) { // This is a new hit
106 auto *anticoincidenceHit = new GammaRayTelAnticoincidenceHit;
107 anticoincidenceHit->SetACDType(0);
108 anticoincidenceHit->AddEnergy(depositedEnergy);
109 anticoincidenceHit->SetPosition(step->GetPreStepPoint()->GetPosition());
110 anticoincidenceHit->SetACDTileNumber(acdTileNumber);
111 hitTopID[acdTileNumber] = anticoincidenceCollection->insert(anticoincidenceHit) - 1;
112 } else { // This is not new
113 (*anticoincidenceCollection)[hitTopID[acdTileNumber]]->AddEnergy(depositedEnergy);
114 }
115 }
116
117 if (acdTileName == "ACL1") { // The hit is on a lateral (left-right) ACD tile (ACDType: 1)
118 if (hitLateralID[acdTileNumber] == -1) { // This is a new hit
119 auto *anticoincidenceHit = new GammaRayTelAnticoincidenceHit;
120 anticoincidenceHit->SetACDType(1);
121 anticoincidenceHit->AddEnergy(depositedEnergy);
122 anticoincidenceHit->SetPosition(step->GetPreStepPoint()->GetPosition());
123 anticoincidenceHit->SetACDTileNumber(acdTileNumber);
124 hitLateralID[acdTileNumber] = anticoincidenceCollection->insert(anticoincidenceHit) - 1;
125 } else { // This is not new
126 (*anticoincidenceCollection)[hitLateralID[acdTileNumber]]->AddEnergy(depositedEnergy);
127 }
128 }
129
130 if (acdTileName == "ACL2") { // The hit is on a lateral (rear-front) ACD tile (ACDType: 2)
131 if (hitLateralID[acdTileNumber] == -1) { // This is a new hit
132 auto *anticoincidenceHit = new GammaRayTelAnticoincidenceHit;
133 anticoincidenceHit->SetACDType(2);
134 anticoincidenceHit->AddEnergy(depositedEnergy);
135 anticoincidenceHit->SetPosition(step->GetPreStepPoint()->GetPosition());
136 anticoincidenceHit->SetACDTileNumber(acdTileNumber);
137 hitLateralID[acdTileNumber] = anticoincidenceCollection->insert(anticoincidenceHit) - 1;
138 } else { // This is not new
139 (*anticoincidenceCollection)[hitLateralID[acdTileNumber]]->AddEnergy(depositedEnergy);
140 }
141 }
142
143 return true;
144 }
145
146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
147
148 void GammaRayTelAnticoincidenceSD::EndOfEvent(G4HCofThisEvent *collection) {
149 static G4int collectionIdentifier = -1;
150 if (collectionIdentifier < 0) {
151 collectionIdentifier = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
152 }
153 collection->AddHitsCollection(collectionIdentifier, anticoincidenceCollection);
154
155 for (auto i = 0; i < numberOfACDLateralTiles; i++) {
156 hitLateralID[i] = -1;
157 }
158
159 for (auto j = 0; j < numberOfACDTopTiles; j++) {
160 hitTopID[j] = -1;
161 }
162 }
163
164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
165
166 void GammaRayTelAnticoincidenceSD::clear() {
167 }
168
169 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
170
171 void GammaRayTelAnticoincidenceSD::DrawAll() {
172 }
173
174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
175
176 void GammaRayTelAnticoincidenceSD::PrintAll() {
177 }
178