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 //GEANT4 - Depth-of-Interaction enabled Positr 27 28 //Authors and contributors 29 30 // Author list to be updated, with names of co 31 32 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 33 // Hideaki Tashima (3), Go Akamatsu (3), Akram 34 // Susanna Guatelli (2), and Mitra Safavi-Naei 35 36 // (1) Australian Nuclear Science and Technolo 37 // (2) University of Wollongong, Australia 38 // (3) National Institute of Radiological Scie 39 40 41 #include "doiPETSteppingAction.hh" 42 #include "doiPETAnalysis.hh" 43 #include "doiPETRun.hh" 44 #include "doiPETEventAction.hh" 45 46 #include "G4Step.hh" 47 #include "G4RunManager.hh" 48 #include "G4SteppingManager.hh" 49 50 #include "G4SystemOfUnits.hh" 51 #include "G4PhysicalConstants.hh" 52 #include <limits> 53 54 ////////// Constructor /////////////////////// 55 /*doiPETSteppingAction::doiPETSteppingAction() 56 {;}*/ 57 doiPETSteppingAction::doiPETSteppingAction() 58 : G4UserSteppingAction() 59 { } 60 61 ///////// Destructor ///////////////////////// 62 doiPETSteppingAction::~doiPETSteppingAction() 63 {;} 64 65 ///////// UserSteppingAction ///////////////// 66 void doiPETSteppingAction::UserSteppingAction( 67 { 68 doiPETRun* run = static_cast<doiPETRun*>(G4R 69 70 G4Track* track = aStep->GetTrack(); 71 G4String volumeName = aStep->GetPreStepPoint 72 G4double edep = aStep->GetTotalEnergyDeposit 73 G4String processName =aStep->GetPostStepPoin 74 G4String particleName = track->GetDynamicPar 75 //G4int eventID = G4RunManager::GetRunManage 76 G4ThreeVector pos = track->GetPosition(); 77 78 79 //G4StepPoint* p2 = aStep->GetPostStepPoint( 80 G4StepPoint* p1 = aStep->GetPreStepPoint(); 81 G4ThreeVector coord1 = p1->GetPosition(); 82 83 //The following is to get the local position 84 const G4AffineTransform transformation = p1 85 G4ThreeVector localPosition = transformation 86 87 G4int blockID; 88 G4int crystalID; 89 90 G4int scatterIndex = 0;//For checking 91 92 //If an event is created in a cold region or 93 if((volumeName != "phantom_physicalV" || vol 94 track->SetTrackStatus(fStopAndKill); 95 return; 96 } 97 98 //Get the source position if the process nam 99 if(processName == "annihil" && volumeName==" 100 doiPETAnalysis::GetInstance()->SetSourcePo 101 run->SetAnnihilationTime((track->GetGlobal 102 } 103 104 //get event ID 105 //doiPETAnalysis::GetInstance()->SetEventID 106 107 //Get scatter information in the phantom by 108 //If there energy deposition in the phantom 109 if(edep > 0 && particleName == "gamma" && v 110 scatterIndex = 1; 111 doiPETAnalysis::GetInstance()->SetScatterI 112 } 113 114 /////////////////// Retrive (Extract) inf 115 if(edep>0. && volumeName=="Crystal_physicalV 116 117 //get the copy number of the block 118 blockID = aStep->GetPreStepPoint()->GetTou 119 120 //get the crystal copy number 121 crystalID = aStep->GetPreStepPoint()->GetT 122 123 //G4cout<<localPosition.x()<<" "<<localPos 124 125 //get event ID 126 //doiPETAnalysis::GetInstance()->GetEventI 127 128 //define a pointer for extracting informat 129 InteractionInformation* ExtractIntInfo = n 130 131 //get the energy deposition of each intera 132 ExtractIntInfo->SetEdep( edep ); 133 134 135 //get the blockID 136 ExtractIntInfo->SetBlockNo( blockID ); 137 138 139 //get the crystal ID 140 ExtractIntInfo->SetCrystalNo( crystalID ); 141 142 143 //get local position interaction in the cr 144 ExtractIntInfo->SetInteractionPositionInCr 145 146 //get the global time of the interaction w 147 ExtractIntInfo->SetGlobalTime( track->GetG 148 149 //pass all the obtained information to the 150 //doiPETAnalysis::GetInstance()->GetIntrac 151 run->GetIntractionInfomation(ExtractIntInf 152 153 } 154 } 155