ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
TrackingAction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file TrackingAction.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collaboration publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // and papers
31 // M. Batmunkh et al. J Radiat Res Appl Sci 8 (2015) 498-507
32 // O. Belov et al. Physica Medica 32 (2016) 1510-1520
33 // The Geant4-DNA web site is available at http://geant4-dna.org
34 //
35 // -------------------------------------------------------------------
36 // November 2016
37 // -------------------------------------------------------------------
38 //
41 
42 #include "TrackingAction.hh"
43 #include "G4Track.hh"
44 #include "G4VSolid.hh"
45 #include "G4Region.hh"
46 #include "G4Electron.hh"
47 #include "G4Gamma.hh"
48 #include "G4RegionStore.hh"
49 //
50 #include "G4UnitsTable.hh"
51 #include "G4SystemOfUnits.hh"
52 #include "G4ios.hh"
53 #include "G4LogicalVolume.hh"
54 #include "G4LogicalVolumeStore.hh"
55 #include "G4VPhysicalVolume.hh"
56 #include "G4PhysicalVolumeStore.hh"
57 //
58 // Bosons
59 #include "G4ChargedGeantino.hh"
60 #include "G4Geantino.hh"
61 #include "G4Gamma.hh"
62 // leptons
63 #include "G4Electron.hh"
64 #include "G4Positron.hh"
65 #include "G4MuonPlus.hh"
66 #include "G4MuonMinus.hh"
67 #include "G4NeutrinoMu.hh"
68 #include "G4AntiNeutrinoMu.hh"
69 #include "G4NeutrinoE.hh"
70 #include "G4AntiNeutrinoE.hh"
71 // Mesons
72 #include "G4PionPlus.hh"
73 #include "G4PionMinus.hh"
74 #include "G4PionZero.hh"
75 #include "G4Eta.hh"
76 #include "G4EtaPrime.hh"
77 
78 #include "G4KaonPlus.hh"
79 #include "G4KaonMinus.hh"
80 #include "G4KaonZero.hh"
81 #include "G4AntiKaonZero.hh"
82 #include "G4KaonZeroLong.hh"
83 #include "G4KaonZeroShort.hh"
84 // Baryons
85 #include "G4Proton.hh"
86 #include "G4AntiProton.hh"
87 #include "G4Neutron.hh"
88 #include "G4AntiNeutron.hh"
89 // Nuclei
90 #include "G4Alpha.hh"
91 #include "G4Deuteron.hh"
92 #include "G4Triton.hh"
93 #include "G4He3.hh"
94 #include "G4GenericIon.hh"
95 #include "G4Step.hh"
96 #include "G4StepStatus.hh"
97 #include "Run.hh"
98 #include "G4RunManager.hh"
99 
100 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
101 
104  RunInitObserver(), fpTargetRegion(0)
105 {}
106 
107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
108 
110 {
111  fpTargetRegion = 0;
112 }
113 
114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
115 
117 {
118 
119 // Target volumes
120  G4VPhysicalVolume* volumeStep = track->GetTouchableHandle()->GetVolume();
121  G4VPhysicalVolume* volumeMedium =
123  G4VPhysicalVolume* volumeSlice =
124  G4PhysicalVolumeStore::GetInstance()->GetVolume("BoundingSlice");
125  Run* run = static_cast<Run*>(
127 
128  //count secondary particles
129  if (track->GetTrackID() == 1) return;
131  G4double energy = track->GetKineticEnergy();
132 
133  // outside bounding slice
134  //if (volumeStep == volumeMedium) run->ParticleCount(name,energy);
135  // inside bounding slice
136  //if (volumeStep != volumeMedium) run->ParticleCountNeuron(name,energy);
137 
138  // particles outside neuron structure
139  if (volumeStep == volumeMedium || volumeStep == volumeSlice)
140  {
141  run->ParticleCount(name,energy);
142  }
143  // count secondary particles in neuron
144  else //if (volumeStep != volumeMedium && volumeStep != volumeSlice)
145  {
146  run->ParticleCountNeuron(name,energy);
147  }
148 
149 }
150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
151 
153 {
154 
155  if(track->GetTrackID() == 1)
156  {
157  G4cout<<"End of tracking primary particle, its final energy is :"
158  << G4BestUnit(track->GetKineticEnergy(), "Energy")<< G4endl;
159  }
160 
161  Run* run = static_cast<Run*>(
163  G4double tracklen = track->GetTrackLength();
164  if (track->GetTrackID() == 1) run->SetTrackLength(tracklen);
165  //G4cout<< " track length = " << tracklen / um << G4endl;
166 }
167 
168 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
169 void
171 {
172  //fpTargetRegion = G4RegionStore::GetInstance()->GetRegion("BoundingSlice");
173 }