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 // ------------ GammaRayTelTrackerSD --- 32 // ------------ GammaRayTelTrackerSD ------ 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 "GammaRayTelTrackerSD.hh" 37 #include "GammaRayTelTrackerSD.hh" 38 #include "GammaRayTelTrackerHit.hh" 38 #include "GammaRayTelTrackerHit.hh" 39 #include "GammaRayTelDetectorConstruction.hh" 39 #include "GammaRayTelDetectorConstruction.hh" 40 40 41 #include "G4RunManager.hh" 41 #include "G4RunManager.hh" 42 #include "G4SystemOfUnits.hh" 42 #include "G4SystemOfUnits.hh" 43 #include "G4VPhysicalVolume.hh" 43 #include "G4VPhysicalVolume.hh" 44 #include "G4Step.hh" 44 #include "G4Step.hh" 45 #include "G4VTouchable.hh" 45 #include "G4VTouchable.hh" 46 #include "G4TouchableHistory.hh" 46 #include "G4TouchableHistory.hh" 47 #include "G4SDManager.hh" 47 #include "G4SDManager.hh" 48 #include "G4ios.hh" 48 #include "G4ios.hh" 49 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 51 51 52 GammaRayTelTrackerSD::GammaRayTelTrackerSD(G4S 52 GammaRayTelTrackerSD::GammaRayTelTrackerSD(G4String name) : G4VSensitiveDetector(name) { 53 auto *runManager = G4RunManager::GetRunMan 53 auto *runManager = G4RunManager::GetRunManager(); 54 54 55 detector = (GammaRayTelDetectorConstructio 55 detector = (GammaRayTelDetectorConstruction*) (runManager->GetUserDetectorConstruction()); 56 56 57 auto numberOfTKRTiles = detector->GetNbOfT 57 auto numberOfTKRTiles = detector->GetNbOfTKRTiles(); 58 numberOfTKRStrips = detector->GetNbOfTKRSt 58 numberOfTKRStrips = detector->GetNbOfTKRStrips(); 59 numberOfTKRLayers = detector->GetNbOfTKRLa 59 numberOfTKRLayers = detector->GetNbOfTKRLayers(); 60 numberOfTKRStrips = numberOfTKRStrips * nu 60 numberOfTKRStrips = numberOfTKRStrips * numberOfTKRTiles; 61 numberOfTKRChannels = numberOfTKRStrips * 61 numberOfTKRChannels = numberOfTKRStrips * numberOfTKRTiles * numberOfTKRLayers; 62 62 63 tkrHitXID = new G4int[numberOfTKRChannels] 63 tkrHitXID = new G4int[numberOfTKRChannels]; 64 tkrHitYID = new G4int[numberOfTKRChannels] 64 tkrHitYID = new G4int[numberOfTKRChannels]; 65 65 66 constexpr auto TRACKER_COLLECTION_NAME = " 66 constexpr auto TRACKER_COLLECTION_NAME = "TrackerCollection"; 67 collectionName.insert(TRACKER_COLLECTION_N 67 collectionName.insert(TRACKER_COLLECTION_NAME); 68 } 68 } 69 69 70 //....oooOO0OOooo........oooOO0OOooo........oo 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 71 71 72 GammaRayTelTrackerSD::~GammaRayTelTrackerSD() 72 GammaRayTelTrackerSD::~GammaRayTelTrackerSD() { 73 delete[] tkrHitXID; 73 delete[] tkrHitXID; 74 delete[] tkrHitYID; 74 delete[] tkrHitYID; 75 } 75 } 76 76 77 //....oooOO0OOooo........oooOO0OOooo........oo 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 78 78 79 void GammaRayTelTrackerSD::Initialize(G4HCofTh 79 void GammaRayTelTrackerSD::Initialize(G4HCofThisEvent*) { 80 trackerCollection = new GammaRayTelTracker 80 trackerCollection = new GammaRayTelTrackerHitsCollection(SensitiveDetectorName, collectionName[0]); 81 81 82 for (auto i = 0; i < numberOfTKRChannels; 82 for (auto i = 0; i < numberOfTKRChannels; i++) { 83 tkrHitXID[i] = -1; 83 tkrHitXID[i] = -1; 84 tkrHitYID[i] = -1; 84 tkrHitYID[i] = -1; 85 }; 85 }; 86 } 86 } 87 87 88 //....oooOO0OOooo........oooOO0OOooo........oo 88 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 89 89 90 auto GammaRayTelTrackerSD::ProcessHits(G4Step 90 auto GammaRayTelTrackerSD::ProcessHits(G4Step *step, G4TouchableHistory*) -> G4bool { 91 G4double depositedEnergy = 0.; 91 G4double depositedEnergy = 0.; 92 depositedEnergy = step->GetTotalEnergyDepo 92 depositedEnergy = step->GetTotalEnergyDeposit(); 93 if (depositedEnergy == 0.) { 93 if (depositedEnergy == 0.) { 94 return false; 94 return false; 95 } 95 } 96 96 97 auto totalNumberOfStrips = detector->GetNb 97 auto totalNumberOfStrips = detector->GetNbOfTKRStrips(); 98 auto totalnumberOfTiles = detector->GetNbO 98 auto totalnumberOfTiles = detector->GetNbOfTKRTiles(); 99 99 100 // This TouchableHistory is used to obtain 100 // This TouchableHistory is used to obtain the physical volume of the hit 101 auto *touchable = (G4TouchableHistory*) (s 101 auto *touchable = (G4TouchableHistory*) (step->GetPreStepPoint()->GetTouchable()); 102 auto *plane = touchable->GetVolume(2); 102 auto *plane = touchable->GetVolume(2); 103 103 104 G4int planeNumber = 0; 104 G4int planeNumber = 0; 105 planeNumber = plane->GetCopyNo(); 105 planeNumber = plane->GetCopyNo(); 106 auto planeName = plane->GetName(); 106 auto planeName = plane->GetName(); 107 107 108 // The hits sees now the real strip 108 // The hits sees now the real strip 109 109 110 G4int stripNumber = 0; 110 G4int stripNumber = 0; 111 G4VPhysicalVolume *strip{nullptr}; 111 G4VPhysicalVolume *strip{nullptr}; 112 strip = touchable->GetVolume(); 112 strip = touchable->GetVolume(); 113 113 114 G4String stripName = strip->GetName(); 114 G4String stripName = strip->GetName(); 115 stripNumber = strip->GetCopyNo(); 115 stripNumber = strip->GetCopyNo(); 116 116 117 auto *tile = touchable->GetVolume(1); 117 auto *tile = touchable->GetVolume(1); 118 auto tileNumber = tile->GetCopyNo(); 118 auto tileNumber = tile->GetCopyNo(); 119 auto tileName = tile->GetName(); 119 auto tileName = tile->GetName(); 120 auto NTile = (tileNumber % totalnumberOfTi 120 auto NTile = (tileNumber % totalnumberOfTiles); 121 121 122 G4int channelNumber = 0; 122 G4int channelNumber = 0; 123 123 124 for (auto j = 0; j < totalnumberOfTiles; j 124 for (auto j = 0; j < totalnumberOfTiles; j++) { 125 if (NTile == j) { 125 if (NTile == j) { 126 stripNumber += totalNumberOfStrips 126 stripNumber += totalNumberOfStrips * NTile; 127 } 127 } 128 } 128 } 129 129 130 channelNumber = planeNumber * totalnumberO 130 channelNumber = planeNumber * totalnumberOfTiles * totalNumberOfStrips + stripNumber; 131 131 132 /* 132 /* 133 G4cout << " Channel: " << channelNumber << 133 G4cout << " Channel: " << channelNumber << G4endl; 134 G4cout << " Plane: " << planeNumber << " " 134 G4cout << " Plane: " << planeNumber << " " << planeName << G4endl; 135 G4cout << " Strip: " << stripNumber << " " 135 G4cout << " Strip: " << stripNumber << " " << stripName << G4endl; 136 */ 136 */ 137 137 138 // The hit is on an X silicon plane 138 // The hit is on an X silicon plane 139 if (planeName == "TKRDetectorX") { 139 if (planeName == "TKRDetectorX") { 140 if (tkrHitXID[channelNumber] == -1) { 140 if (tkrHitXID[channelNumber] == -1) { // This is a new hit 141 auto *trackerHit = new GammaRayTel 141 auto *trackerHit = new GammaRayTelTrackerHit; 142 trackerHit->SetPlaneType(1); 142 trackerHit->SetPlaneType(1); 143 trackerHit->AddDepositedEnergy(dep 143 trackerHit->AddDepositedEnergy(depositedEnergy); 144 trackerHit->SetPosition(step->GetP 144 trackerHit->SetPosition(step->GetPreStepPoint()->GetPosition()); 145 trackerHit->SetSiliconPlaneNumber( 145 trackerHit->SetSiliconPlaneNumber(planeNumber); 146 trackerHit->SetStripNumber(stripNu 146 trackerHit->SetStripNumber(stripNumber); 147 tkrHitXID[channelNumber] = tracker 147 tkrHitXID[channelNumber] = trackerCollection->insert(trackerHit) - 1; 148 } else { // This is not new 148 } else { // This is not new 149 (*trackerCollection)[tkrHitXID[cha 149 (*trackerCollection)[tkrHitXID[channelNumber]]->AddDepositedEnergy(depositedEnergy); 150 // G4cout << "X" << planeNumber << 150 // G4cout << "X" << planeNumber << " " << stripNumber << G4endl; 151 } 151 } 152 } 152 } 153 153 154 // The hit is on an Y silicon plane 154 // The hit is on an Y silicon plane 155 if (planeName == "TKRDetectorY") { 155 if (planeName == "TKRDetectorY") { 156 if (tkrHitYID[channelNumber] == -1) { 156 if (tkrHitYID[channelNumber] == -1) { // This is a new hit 157 auto *trackerHit = new GammaRayTel 157 auto *trackerHit = new GammaRayTelTrackerHit; 158 trackerHit->SetPlaneType(0); 158 trackerHit->SetPlaneType(0); 159 trackerHit->AddDepositedEnergy(dep 159 trackerHit->AddDepositedEnergy(depositedEnergy); 160 trackerHit->SetPosition(step->GetP 160 trackerHit->SetPosition(step->GetPreStepPoint()->GetPosition()); 161 trackerHit->SetSiliconPlaneNumber( 161 trackerHit->SetSiliconPlaneNumber(planeNumber); 162 trackerHit->SetStripNumber(stripNu 162 trackerHit->SetStripNumber(stripNumber); 163 tkrHitYID[channelNumber] = tracker 163 tkrHitYID[channelNumber] = trackerCollection->insert(trackerHit) - 1; 164 } else { // This is not new 164 } else { // This is not new 165 (*trackerCollection)[tkrHitYID[cha 165 (*trackerCollection)[tkrHitYID[channelNumber]]->AddDepositedEnergy(depositedEnergy); 166 // G4cout << "Y" << planeNumber << 166 // G4cout << "Y" << planeNumber << " " << stripNumber << G4endl; 167 } 167 } 168 } 168 } 169 169 170 return true; 170 return true; 171 } 171 } 172 172 173 //....oooOO0OOooo........oooOO0OOooo........oo 173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 174 174 175 void GammaRayTelTrackerSD::EndOfEvent(G4HCofTh 175 void GammaRayTelTrackerSD::EndOfEvent(G4HCofThisEvent *collection) { 176 static G4int collectionIdentifier = -1; 176 static G4int collectionIdentifier = -1; 177 if (collectionIdentifier < 0) { 177 if (collectionIdentifier < 0) { 178 collectionIdentifier = G4SDManager::Ge 178 collectionIdentifier = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); 179 } 179 } 180 collection->AddHitsCollection(collectionId 180 collection->AddHitsCollection(collectionIdentifier, trackerCollection); 181 181 182 for (auto i = 0; i < numberOfTKRChannels; 182 for (auto i = 0; i < numberOfTKRChannels; i++) { 183 tkrHitXID[i] = -1; 183 tkrHitXID[i] = -1; 184 tkrHitYID[i] = -1; 184 tkrHitYID[i] = -1; 185 } 185 } 186 } 186 } 187 187 188 //....oooOO0OOooo........oooOO0OOooo........oo 188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 189 189 190 void GammaRayTelTrackerSD::clear() { 190 void GammaRayTelTrackerSD::clear() { 191 } 191 } 192 192 193 //....oooOO0OOooo........oooOO0OOooo........oo 193 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 194 194 195 void GammaRayTelTrackerSD::DrawAll() { 195 void GammaRayTelTrackerSD::DrawAll() { 196 } 196 } 197 197 198 //....oooOO0OOooo........oooOO0OOooo........oo 198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 199 199 200 void GammaRayTelTrackerSD::PrintAll() { 200 void GammaRayTelTrackerSD::PrintAll() { 201 } 201 } 202 202