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 // $Id: GammaRayTelAnticoincidenceSD.cc,v 1.6 2006/06/29 15:56:15 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-08-01-patch-01 $ 27 // ------------------------------------------- 29 // ------------------------------------------------------------ 28 // GEANT 4 class implementation file 30 // GEANT 4 class implementation file 29 // CERN Geneva Switzerland 31 // CERN Geneva Switzerland 30 // 32 // 31 // 33 // 32 // ------------ GammaRayTelAnticoincidenc 34 // ------------ GammaRayTelAnticoincidenceSD ------ 33 // by R.Giannitrapani, F.Longo & G. 35 // by R.Giannitrapani, F.Longo & G.Santin (13 nov 2000) 34 // 36 // 35 // ******************************************* 37 // ************************************************************ 36 << 38 #include "G4RunManager.hh" 37 #include "GammaRayTelAnticoincidenceSD.hh" 39 #include "GammaRayTelAnticoincidenceSD.hh" 38 #include "GammaRayTelAnticoincidenceHit.hh" 40 #include "GammaRayTelAnticoincidenceHit.hh" 39 #include "GammaRayTelDetectorConstruction.hh" 41 #include "GammaRayTelDetectorConstruction.hh" 40 42 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" 43 #include "G4VPhysicalVolume.hh" >> 44 #include "G4Step.hh" 48 #include "G4VTouchable.hh" 45 #include "G4VTouchable.hh" >> 46 #include "G4TouchableHistory.hh" >> 47 #include "G4SDManager.hh" >> 48 #include "G4ios.hh" 49 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 51 51 52 GammaRayTelAnticoincidenceSD::GammaRayTelAntic << 52 GammaRayTelAnticoincidenceSD::GammaRayTelAnticoincidenceSD(G4String name) 53 auto *runManager = G4RunManager::GetRunMan << 53 :G4VSensitiveDetector(name) 54 detector = (GammaRayTelDetectorConstructio << 54 { 55 << 55 G4RunManager* runManager = G4RunManager::GetRunManager(); 56 numberOfACDLateralTiles = detector->GetNbO << 56 Detector = 57 numberOfACDTopTiles = detector->GetNbOfACD << 57 (GammaRayTelDetectorConstruction*)(runManager->GetUserDetectorConstruction()); 58 << 58 59 G4cout << numberOfACDLateralTiles << " LAT << 59 NbOfACDLateralTiles = Detector->GetNbOfACDLateralTiles(); 60 G4cout << numberOfACDTopTiles << " TOP " < << 60 NbOfACDTopTiles = Detector->GetNbOfACDTopTiles(); 61 << 61 62 hitLateralID = new G4int[numberOfACDLatera << 62 G4cout << NbOfACDLateralTiles << " LAT " << G4endl; 63 hitTopID = new G4int[numberOfACDTopTiles]; << 63 G4cout << NbOfACDTopTiles << " TOP " << G4endl; 64 collectionName.insert("AnticoincidenceColl << 64 >> 65 HitLateralID = new G4int[NbOfACDLateralTiles]; >> 66 HitTopID = new G4int[NbOfACDTopTiles]; >> 67 collectionName.insert("AnticoincidenceCollection"); >> 68 } >> 69 >> 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 71 >> 72 GammaRayTelAnticoincidenceSD::~GammaRayTelAnticoincidenceSD() >> 73 { >> 74 delete [] HitLateralID; >> 75 delete [] HitTopID; >> 76 } >> 77 >> 78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 79 >> 80 void GammaRayTelAnticoincidenceSD::Initialize(G4HCofThisEvent*) >> 81 { >> 82 AnticoincidenceCollection = new GammaRayTelAnticoincidenceHitsCollection >> 83 (SensitiveDetectorName,collectionName[0]); >> 84 for (G4int i=0;i<NbOfACDLateralTiles;i++) >> 85 { >> 86 HitLateralID[i] = -1; >> 87 }; >> 88 >> 89 for (G4int j=0;j<NbOfACDTopTiles;j++) >> 90 { >> 91 >> 92 HitTopID[j] = -1; >> 93 }; >> 94 } >> 95 >> 96 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 97 >> 98 G4bool GammaRayTelAnticoincidenceSD::ProcessHits(G4Step* aStep, G4TouchableHistory*) >> 99 { >> 100 >> 101 G4double edep = aStep->GetTotalEnergyDeposit(); >> 102 if ((edep/keV == 0.)) return false; >> 103 >> 104 // This TouchableHistory is used to obtain the physical volume >> 105 // of the hit >> 106 G4TouchableHistory* theTouchable >> 107 = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable()); >> 108 >> 109 G4VPhysicalVolume* acd_tile = theTouchable->GetVolume(); >> 110 >> 111 G4int ACDTileNumber=acd_tile->GetCopyNo(); >> 112 G4String ACDTileName = acd_tile->GetName(); >> 113 >> 114 // G4cout << ACDTileName << " " << edep/keV << G4endl; >> 115 >> 116 if (ACDTileName == "ACT" ) >> 117 // The hit is on an top ACD tile (ACDType 0) >> 118 { >> 119 // This is a new hit >> 120 if (HitTopID[ACDTileNumber]==-1) >> 121 { >> 122 GammaRayTelAnticoincidenceHit* >> 123 AnticoincidenceHit = new GammaRayTelAnticoincidenceHit; >> 124 AnticoincidenceHit->SetACDType(0); >> 125 AnticoincidenceHit->AddEnergy(edep); >> 126 AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition()); >> 127 AnticoincidenceHit->SetACDTileNumber(ACDTileNumber); >> 128 HitTopID[ACDTileNumber] = >> 129 AnticoincidenceCollection->insert(AnticoincidenceHit) -1; >> 130 } >> 131 else // This is not new >> 132 { >> 133 (*AnticoincidenceCollection) >> 134 [HitTopID[ACDTileNumber]]->AddEnergy(edep); >> 135 } >> 136 } >> 137 >> 138 if (ACDTileName == "ACL1") >> 139 // The hit is on an lateral (left-right) ACD tile (ACDType 1) >> 140 { >> 141 // This is a new hit >> 142 if (HitLateralID[ACDTileNumber]==-1) >> 143 { >> 144 GammaRayTelAnticoincidenceHit* >> 145 AnticoincidenceHit = new GammaRayTelAnticoincidenceHit; >> 146 AnticoincidenceHit->SetACDType(1); >> 147 AnticoincidenceHit->AddEnergy(edep); >> 148 AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition()); >> 149 AnticoincidenceHit->SetACDTileNumber(ACDTileNumber); >> 150 HitLateralID[ACDTileNumber] = >> 151 AnticoincidenceCollection->insert(AnticoincidenceHit) -1; >> 152 } >> 153 else // This is not new >> 154 { >> 155 (*AnticoincidenceCollection) >> 156 [HitLateralID[ACDTileNumber]]->AddEnergy(edep); >> 157 } >> 158 } >> 159 >> 160 if (ACDTileName == "ACL2") >> 161 // The hit is on an lateral (rear - front) ACD tile (ACDType 2) >> 162 { >> 163 // This is a new hit >> 164 if (HitLateralID[ACDTileNumber]==-1) >> 165 { >> 166 GammaRayTelAnticoincidenceHit* >> 167 AnticoincidenceHit = new GammaRayTelAnticoincidenceHit; >> 168 AnticoincidenceHit->SetACDType(2); >> 169 AnticoincidenceHit->AddEnergy(edep); >> 170 AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition()); >> 171 AnticoincidenceHit->SetACDTileNumber(ACDTileNumber); >> 172 HitLateralID[ACDTileNumber] = >> 173 AnticoincidenceCollection->insert(AnticoincidenceHit) -1; >> 174 } >> 175 else // This is not new >> 176 { >> 177 (*AnticoincidenceCollection) >> 178 [HitLateralID[ACDTileNumber]]->AddEnergy(edep); >> 179 } >> 180 } >> 181 return true; >> 182 } >> 183 >> 184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 185 >> 186 void GammaRayTelAnticoincidenceSD::EndOfEvent(G4HCofThisEvent* HCE) >> 187 { >> 188 static G4int HCID = -1; >> 189 if(HCID<0) >> 190 { >> 191 HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); >> 192 } >> 193 HCE->AddHitsCollection(HCID,AnticoincidenceCollection); >> 194 >> 195 for (G4int i=0;i<NbOfACDLateralTiles;i++) >> 196 { >> 197 HitLateralID[i] = -1; >> 198 }; >> 199 >> 200 for (G4int j=0;j<NbOfACDTopTiles;j++) >> 201 { >> 202 >> 203 HitTopID[j] = -1; >> 204 }; >> 205 65 } 206 } 66 207 67 //....oooOO0OOooo........oooOO0OOooo........oo 208 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 68 209 69 GammaRayTelAnticoincidenceSD::~GammaRayTelAnti << 210 void GammaRayTelAnticoincidenceSD::clear() 70 delete[] hitLateralID; << 211 {} 71 delete[] hitTopID; << 72 } << 73 212 74 //....oooOO0OOooo........oooOO0OOooo........oo 213 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 75 214 76 void GammaRayTelAnticoincidenceSD::Initialize( << 215 void GammaRayTelAnticoincidenceSD::DrawAll() 77 anticoincidenceCollection = new GammaRayTe << 216 {} 78 for (auto i = 0; i < numberOfACDLateralTil << 79 hitLateralID[i] = -1; << 80 } << 81 << 82 for (auto j = 0; j < numberOfACDTopTiles; << 83 hitTopID[j] = -1; << 84 } << 85 } << 86 217 87 //....oooOO0OOooo........oooOO0OOooo........oo 218 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 88 219 89 auto GammaRayTelAnticoincidenceSD::ProcessHits << 220 void GammaRayTelAnticoincidenceSD::PrintAll() 90 auto depositedEnergy = step->GetTotalEnerg << 221 {} 91 if (depositedEnergy / keV == 0.) { << 92 return false; << 93 } << 94 222 95 // This TouchableHistory is used to obtain << 223 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 96 auto *theTouchable = (G4TouchableHistory*) << 97 224 98 auto *acdTile = theTouchable->GetVolume(); << 99 auto acdTileNumber = acdTile->GetCopyNo(); << 100 auto acdTileName = acdTile->GetName(); << 101 << 102 // G4cout << acdTileName << " " << deposit << 103 << 104 if (acdTileName == "ACT") { // The hit is << 105 if (hitTopID[acdTileNumber] == -1) { / << 106 auto *anticoincidenceHit = new Gam << 107 anticoincidenceHit->SetACDType(0); << 108 anticoincidenceHit->AddEnergy(depo << 109 anticoincidenceHit->SetPosition(st << 110 anticoincidenceHit->SetACDTileNumb << 111 hitTopID[acdTileNumber] = anticoin << 112 } else { // This is not new << 113 (*anticoincidenceCollection)[hitTo << 114 } << 115 } << 116 225 117 if (acdTileName == "ACL1") { // The hit is << 118 if (hitLateralID[acdTileNumber] == -1) << 119 auto *anticoincidenceHit = new Gam << 120 anticoincidenceHit->SetACDType(1); << 121 anticoincidenceHit->AddEnergy(depo << 122 anticoincidenceHit->SetPosition(st << 123 anticoincidenceHit->SetACDTileNumb << 124 hitLateralID[acdTileNumber] = anti << 125 } else { // This is not new << 126 (*anticoincidenceCollection)[hitLa << 127 } << 128 } << 129 226 130 if (acdTileName == "ACL2") { // The hit is << 131 if (hitLateralID[acdTileNumber] == -1) << 132 auto *anticoincidenceHit = new Gam << 133 anticoincidenceHit->SetACDType(2); << 134 anticoincidenceHit->AddEnergy(depo << 135 anticoincidenceHit->SetPosition(st << 136 anticoincidenceHit->SetACDTileNumb << 137 hitLateralID[acdTileNumber] = anti << 138 } else { // This is not new << 139 (*anticoincidenceCollection)[hitLa << 140 } << 141 } << 142 227 143 return true; << 144 } << 145 228 146 //....oooOO0OOooo........oooOO0OOooo........oo << 147 229 148 void GammaRayTelAnticoincidenceSD::EndOfEvent( << 149 static G4int collectionIdentifier = -1; << 150 if (collectionIdentifier < 0) { << 151 collectionIdentifier = G4SDManager::Ge << 152 } << 153 collection->AddHitsCollection(collectionId << 154 230 155 for (auto i = 0; i < numberOfACDLateralTil << 156 hitLateralID[i] = -1; << 157 } << 158 231 159 for (auto j = 0; j < numberOfACDTopTiles; << 160 hitTopID[j] = -1; << 161 } << 162 } << 163 232 164 //....oooOO0OOooo........oooOO0OOooo........oo << 165 233 166 void GammaRayTelAnticoincidenceSD::clear() { << 167 } << 168 234 169 //....oooOO0OOooo........oooOO0OOooo........oo << 170 235 171 void GammaRayTelAnticoincidenceSD::DrawAll() { << 172 } << 173 236 174 //....oooOO0OOooo........oooOO0OOooo........oo << 175 237 176 void GammaRayTelAnticoincidenceSD::PrintAll() << 177 } << 178 238