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