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 // This example is provided by the Geant4-DNA 27 // Any report or published results obtained us 28 // shall cite the following Geant4-DNA collabo 29 // Med. Phys. 45 (2018) e722-e739 30 // Phys. Med. 31 (2015) 861-874 31 // Med. Phys. 37 (2010) 4692-4708 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 33 // 34 // The Geant4-DNA web site is available at htt 35 // 36 /// \file medical/dna/range/src/TrackingAction 37 /// \brief Implementation of the TrackingActio 38 39 #include "TrackingAction.hh" 40 41 #include "HistoManager.hh" 42 #include "PrimaryGeneratorAction.hh" 43 #include "Run.hh" 44 45 #include "G4Positron.hh" 46 #include "G4RunManager.hh" 47 48 //....oooOO0OOooo........oooOO0OOooo........oo 49 50 TrackingAction::TrackingAction(PrimaryGenerato 51 : G4UserTrackingAction(), fPrimary(prim) 52 {} 53 54 //....oooOO0OOooo........oooOO0OOooo........oo 55 56 void TrackingAction::PreUserTrackingAction(con 57 { 58 G4int trackID = track->GetTrackID(); 59 if (trackID == 1) fTrackLength = 0; 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oo 63 64 void TrackingAction::PostUserTrackingAction(co 65 { 66 // Histograms 67 G4AnalysisManager* analysisManager = G4Analy 68 69 G4int trackID = track->GetTrackID(); 70 71 Run* run = static_cast<Run*>(G4RunManager::G 72 73 // 74 75 G4double tracklen = track->GetTrackLength(); 76 77 // *** FOR ELECTRONS 78 79 // Track length of primary particle or charg 80 if (trackID == 1 && track->GetDefinition()-> 81 run->AddTrackLength(tracklen); 82 analysisManager->FillH1(1, tracklen); 83 }; 84 85 // Extract projected range of primary partic 86 if (trackID == 1 && track->GetDefinition()-> 87 G4double pr = 88 (track->GetPosition()) * (fPrimary->GetP 89 run->AddProjRange(pr); 90 analysisManager->FillH1(2, pr); 91 } 92 93 // Extract penetration range of primary part 94 if (trackID == 1 && track->GetDefinition()-> 95 G4double pene = std::sqrt((track->GetPosit 96 run->AddPenetration(pene); 97 analysisManager->FillH1(3, pene); 98 }; 99 100 // *** FOR Geant4-DNA PROTONS, ALPHA AND CHA 101 // track->GetDefinition()->GetPDGMass() > 0 102 103 // Track length of primary particle or charg 104 if (track->GetDefinition()->GetPDGCharge() > 105 && track->GetDefinition() != G4Positron: 106 { 107 // Local cumulation is required before sco 108 fTrackLength = fTrackLength + tracklen; 109 110 if (track->GetKineticEnergy() == 0) { 111 run->AddTrackLength(fTrackLength); 112 analysisManager->FillH1(1, fTrackLength) 113 } 114 } 115 116 // Extract projected range of primary partic 117 if (track->GetDefinition()->GetPDGCharge() > 118 && track->GetDefinition() != G4Positron: 119 { 120 G4double pr = 121 (track->GetPosition()) * (fPrimary->GetP 122 123 if (track->GetKineticEnergy() == 0) { 124 run->AddProjRange(pr); 125 analysisManager->FillH1(2, pr); 126 } 127 } 128 129 // Extract penetration range of primary part 130 if (track->GetDefinition()->GetPDGCharge() > 131 && track->GetDefinition() != G4Positron: 132 { 133 G4double pene = std::sqrt((track->GetPosit 134 135 if (track->GetKineticEnergy() == 0) { 136 run->AddPenetration(pene); 137 analysisManager->FillH1(3, pene); 138 } 139 }; 140 } 141