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 /// \file electromagnetic/TestEm11/src/Trackin 27 /// \brief Implementation of the TrackingActio 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "TrackingAction.hh" 34 35 #include "DetectorConstruction.hh" 36 #include "EventAction.hh" 37 #include "HistoManager.hh" 38 #include "Run.hh" 39 40 #include "G4PhysicalConstants.hh" 41 #include "G4Positron.hh" 42 #include "G4RunManager.hh" 43 #include "G4StepStatus.hh" 44 #include "G4Track.hh" 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 48 TrackingAction::TrackingAction(DetectorConstru 49 :fDetector(det),fEventAct(evt) 50 { } 51 52 //....oooOO0OOooo........oooOO0OOooo........oo 53 54 void TrackingAction::PostUserTrackingAction(co 55 { 56 G4int trackID = track->GetTrackID(); 57 58 G4AnalysisManager* analysisManager = G4Analy 59 Run* run 60 = static_cast<Run*>( 61 G4RunManager::GetRunManager()->GetNonC 62 63 //track length of primary particle or charge 64 // 65 G4double tracklen = track->GetTrackLength(); 66 if (trackID == 1) { 67 run->AddTrackLength(tracklen); 68 analysisManager->FillH1(3, tracklen); 69 } else if (track->GetDefinition()->GetPDGCha 70 analysisManager->FillH1(6, tracklen); 71 72 //extract projected range of primary particl 73 // 74 if (trackID == 1) { 75 G4double x = track->GetPosition().x() + 0. 76 run->AddProjRange(x); 77 analysisManager->FillH1(5, x); 78 } 79 80 //mean step size of primary particle 81 // 82 if (trackID == 1) { 83 G4int nbOfSteps = track->GetCurrentStepNum 84 G4double stepSize = tracklen/nbOfSteps; 85 run->AddStepSize(nbOfSteps,stepSize); 86 } 87 88 //status of primary particle : absorbed, tra 89 // 90 if (trackID == 1) { 91 G4int flag = 0; 92 if (!track->GetNextVolume()) { 93 if (track->GetMomentumDirection().x() > 94 else 95 } 96 run->AddTrackStatus(flag); 97 } 98 99 // where are we ? 100 G4StepStatus status = track->GetStep()->GetP 101 if (status == fWorldBoundary) { 102 G4double eleak = track->GetKineticEnergy() 103 if (track->GetDefinition() == G4Positron:: 104 eleak += 2*electron_mass_c2; 105 } 106 fEventAct->AddEleak(eleak); 107 } 108 } 109 110 //....oooOO0OOooo........oooOO0OOooo........oo 111 112