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