| 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 /// \file runAndEvent/RE01/src/RE01TrackInform 26 /// \file runAndEvent/RE01/src/RE01TrackInformation.cc
27 /// \brief Implementation of the RE01TrackInfo 27 /// \brief Implementation of the RE01TrackInformation class
28 // 28 //
29 // 29 //
>> 30 // $Id: RE01TrackInformation.cc 97671 2016-06-07 08:25:00Z gcosmo $
30 // 31 //
31 32
32 #include "RE01TrackInformation.hh" 33 #include "RE01TrackInformation.hh"
33 <<
34 #include "G4SystemOfUnits.hh" <<
35 #include "G4ios.hh" 34 #include "G4ios.hh"
>> 35 #include "G4SystemOfUnits.hh"
36 36
37 G4ThreadLocal G4Allocator<RE01TrackInformation << 37 G4ThreadLocal G4Allocator<RE01TrackInformation> *
>> 38 aTrackInformationAllocator = 0;
38 39
39 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
40 RE01TrackInformation::RE01TrackInformation() : << 41 RE01TrackInformation::RE01TrackInformation()
>> 42 : G4VUserTrackInformation()
41 { 43 {
42 fOriginalTrackID = 0; << 44 fOriginalTrackID = 0;
43 fParticleDefinition = 0; << 45 fParticleDefinition = 0;
44 fOriginalPosition = G4ThreeVector(0., 0., 0. << 46 fOriginalPosition = G4ThreeVector(0.,0.,0.);
45 fOriginalMomentum = G4ThreeVector(0., 0., 0. << 47 fOriginalMomentum = G4ThreeVector(0.,0.,0.);
46 fOriginalEnergy = 0.; << 48 fOriginalEnergy = 0.;
47 fOriginalTime = 0.; << 49 fOriginalTime = 0.;
48 fTrackingStatus = 1; << 50 fTrackingStatus = 1;
49 fSourceTrackID = -1; << 51 fSourceTrackID = -1;
50 fSourceDefinition = 0; << 52 fSourceDefinition = 0;
51 fSourcePosition = G4ThreeVector(0., 0., 0.); << 53 fSourcePosition = G4ThreeVector(0.,0.,0.);
52 fSourceMomentum = G4ThreeVector(0., 0., 0.); << 54 fSourceMomentum = G4ThreeVector(0.,0.,0.);
53 fSourceEnergy = 0.; << 55 fSourceEnergy = 0.;
54 fSourceTime = 0.; << 56 fSourceTime = 0.;
55 fSuspendedStepID = -1; << 57 fSuspendedStepID = -1;
56 } 58 }
57 59
58 //....oooOO0OOooo........oooOO0OOooo........oo 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59 RE01TrackInformation::RE01TrackInformation(con << 61 RE01TrackInformation::RE01TrackInformation(const G4Track* aTrack)
60 { << 62 : G4VUserTrackInformation()
61 fOriginalTrackID = aTrack->GetTrackID(); << 63 {
62 fParticleDefinition = aTrack->GetDefinition( << 64 fOriginalTrackID = aTrack->GetTrackID();
63 fOriginalPosition = aTrack->GetPosition(); << 65 fParticleDefinition = aTrack->GetDefinition();
64 fOriginalMomentum = aTrack->GetMomentum(); << 66 fOriginalPosition = aTrack->GetPosition();
65 fOriginalEnergy = aTrack->GetTotalEnergy(); << 67 fOriginalMomentum = aTrack->GetMomentum();
66 fOriginalTime = aTrack->GetGlobalTime(); << 68 fOriginalEnergy = aTrack->GetTotalEnergy();
67 fTrackingStatus = 1; << 69 fOriginalTime = aTrack->GetGlobalTime();
68 fSourceTrackID = -1; << 70 fTrackingStatus = 1;
69 fSourceDefinition = 0; << 71 fSourceTrackID = -1;
70 fSourcePosition = G4ThreeVector(0., 0., 0.); << 72 fSourceDefinition = 0;
71 fSourceMomentum = G4ThreeVector(0., 0., 0.); << 73 fSourcePosition = G4ThreeVector(0.,0.,0.);
72 fSourceEnergy = 0.; << 74 fSourceMomentum = G4ThreeVector(0.,0.,0.);
73 fSourceTime = 0.; << 75 fSourceEnergy = 0.;
74 fSuspendedStepID = -1; << 76 fSourceTime = 0.;
>> 77 fSuspendedStepID = -1;
75 } 78 }
76 79
77 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78 RE01TrackInformation ::RE01TrackInformation(co << 81 RE01TrackInformation
>> 82 ::RE01TrackInformation(const RE01TrackInformation* aTrackInfo)
79 : G4VUserTrackInformation() 83 : G4VUserTrackInformation()
80 { 84 {
81 fOriginalTrackID = aTrackInfo->fOriginalTrac << 85 fOriginalTrackID = aTrackInfo->fOriginalTrackID;
82 fParticleDefinition = aTrackInfo->fParticleD << 86 fParticleDefinition = aTrackInfo->fParticleDefinition;
83 fOriginalPosition = aTrackInfo->fOriginalPos << 87 fOriginalPosition = aTrackInfo->fOriginalPosition;
84 fOriginalMomentum = aTrackInfo->fOriginalMom << 88 fOriginalMomentum = aTrackInfo->fOriginalMomentum;
85 fOriginalEnergy = aTrackInfo->fOriginalEnerg << 89 fOriginalEnergy = aTrackInfo->fOriginalEnergy;
86 fOriginalTime = aTrackInfo->fOriginalTime; << 90 fOriginalTime = aTrackInfo->fOriginalTime;
87 fTrackingStatus = aTrackInfo->fTrackingStatu << 91 fTrackingStatus = aTrackInfo->fTrackingStatus;
88 fSourceTrackID = aTrackInfo->fSourceTrackID; << 92 fSourceTrackID = aTrackInfo->fSourceTrackID;
89 fSourceDefinition = aTrackInfo->fSourceDefin << 93 fSourceDefinition = aTrackInfo->fSourceDefinition;
90 fSourcePosition = aTrackInfo->fSourcePositio << 94 fSourcePosition = aTrackInfo->fSourcePosition;
91 fSourceMomentum = aTrackInfo->fSourceMomentu << 95 fSourceMomentum = aTrackInfo->fSourceMomentum;
92 fSourceEnergy = aTrackInfo->fSourceEnergy; << 96 fSourceEnergy = aTrackInfo->fSourceEnergy;
93 fSourceTime = aTrackInfo->fSourceTime; << 97 fSourceTime = aTrackInfo->fSourceTime;
94 fSuspendedStepID = -1; << 98 fSuspendedStepID = -1;
95 } 99 }
96 100
97 //....oooOO0OOooo........oooOO0OOooo........oo 101 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
98 RE01TrackInformation::~RE01TrackInformation() 102 RE01TrackInformation::~RE01TrackInformation()
99 { << 103 {;}
100 ; <<
101 } <<
102 104
103 //....oooOO0OOooo........oooOO0OOooo........oo 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
104 RE01TrackInformation& RE01TrackInformation ::o << 106 RE01TrackInformation& RE01TrackInformation
>> 107 ::operator =(const RE01TrackInformation& aTrackInfo)
105 { 108 {
106 fOriginalTrackID = aTrackInfo.fOriginalTrack << 109 fOriginalTrackID = aTrackInfo.fOriginalTrackID;
107 fParticleDefinition = aTrackInfo.fParticleDe << 110 fParticleDefinition = aTrackInfo.fParticleDefinition;
108 fOriginalPosition = aTrackInfo.fOriginalPosi << 111 fOriginalPosition = aTrackInfo.fOriginalPosition;
109 fOriginalMomentum = aTrackInfo.fOriginalMome << 112 fOriginalMomentum = aTrackInfo.fOriginalMomentum;
110 fOriginalEnergy = aTrackInfo.fOriginalEnergy << 113 fOriginalEnergy = aTrackInfo.fOriginalEnergy;
111 fOriginalTime = aTrackInfo.fOriginalTime; << 114 fOriginalTime = aTrackInfo.fOriginalTime;
112 fTrackingStatus = aTrackInfo.fTrackingStatus << 115 fTrackingStatus = aTrackInfo.fTrackingStatus;
113 fSourceTrackID = aTrackInfo.fSourceTrackID; << 116 fSourceTrackID = aTrackInfo.fSourceTrackID;
114 fSourceDefinition = aTrackInfo.fSourceDefini << 117 fSourceDefinition = aTrackInfo.fSourceDefinition;
115 fSourcePosition = aTrackInfo.fSourcePosition << 118 fSourcePosition = aTrackInfo.fSourcePosition;
116 fSourceMomentum = aTrackInfo.fSourceMomentum << 119 fSourceMomentum = aTrackInfo.fSourceMomentum;
117 fSourceEnergy = aTrackInfo.fSourceEnergy; << 120 fSourceEnergy = aTrackInfo.fSourceEnergy;
118 fSourceTime = aTrackInfo.fSourceTime; << 121 fSourceTime = aTrackInfo.fSourceTime;
119 fSuspendedStepID = -1; << 122 fSuspendedStepID = -1;
120 123
121 return *this; << 124 return *this;
122 } 125 }
123 126
124 //....oooOO0OOooo........oooOO0OOooo........oo 127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125 void RE01TrackInformation::SetSourceTrackInfor 128 void RE01TrackInformation::SetSourceTrackInformation(const G4Track* aTrack)
126 { 129 {
127 fSourceTrackID = aTrack->GetTrackID(); << 130 fSourceTrackID = aTrack->GetTrackID();
128 fSourceDefinition = aTrack->GetDefinition(); << 131 fSourceDefinition = aTrack->GetDefinition();
129 fSourcePosition = aTrack->GetPosition(); << 132 fSourcePosition = aTrack->GetPosition();
130 fSourceMomentum = aTrack->GetMomentum(); << 133 fSourceMomentum = aTrack->GetMomentum();
131 fSourceEnergy = aTrack->GetTotalEnergy(); << 134 fSourceEnergy = aTrack->GetTotalEnergy();
132 fSourceTime = aTrack->GetGlobalTime(); << 135 fSourceTime = aTrack->GetGlobalTime();
133 } 136 }
134 137
135 //....oooOO0OOooo........oooOO0OOooo........oo 138 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
136 void RE01TrackInformation::Print() const 139 void RE01TrackInformation::Print() const
137 { 140 {
138 G4cout << "Source track ID " << fSourceTrack << 141 G4cout
139 << "," << fSourceEnergy / GeV << "[Ge << 142 << "Source track ID " << fSourceTrackID << " ("
140 G4cout << "Original primary track ID " << fO << 143 << fSourceDefinition->GetParticleName() << ","
141 << fParticleDefinition->GetParticleNa << 144 << fSourceEnergy/GeV << "[GeV]) at " << fSourcePosition << G4endl;
142 << G4endl; << 145 G4cout
>> 146 << "Original primary track ID " << fOriginalTrackID << " ("
>> 147 << fParticleDefinition->GetParticleName() << ","
>> 148 << fOriginalEnergy/GeV << "[GeV])" << G4endl;
143 } 149 }
>> 150
144 151