| 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 electromagnetic/TestEm1/src/Tracking << 26 // $Id: TrackingAction.cc,v 1.10 2006-06-29 16:37:36 gunter Exp $
27 /// \brief Implementation of the TrackingActio << 27 // GEANT4 tag $Name: not supported by cvs2svn $
28 // << 28 //
29 // <<
30 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 31
33 #include "TrackingAction.hh" 32 #include "TrackingAction.hh"
34 <<
35 #include "EventAction.hh" <<
36 #include "HistoManager.hh" <<
37 #include "PrimaryGeneratorAction.hh" 33 #include "PrimaryGeneratorAction.hh"
38 #include "Run.hh" << 34 #include "RunAction.hh"
>> 35 #include "HistoManager.hh"
39 36
40 #include "G4Positron.hh" <<
41 #include "G4PhysicalConstants.hh" <<
42 #include "G4RunManager.hh" <<
43 #include "G4StepStatus.hh" <<
44 #include "G4SystemOfUnits.hh" <<
45 #include "G4Track.hh" 37 #include "G4Track.hh"
46 #include "G4UnitsTable.hh" <<
47 38
48 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 40
50 TrackingAction::TrackingAction(PrimaryGenerato << 41 TrackingAction::TrackingAction(PrimaryGeneratorAction* prim, RunAction* run,
51 : fPrimary(prim),fEventAction(evt) {} << 42 HistoManager* histo)
>> 43 :primary(prim), runAction(run), histoManager(histo)
>> 44 { }
52 45
53 //....oooOO0OOooo........oooOO0OOooo........oo 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54 47
55 void TrackingAction::PreUserTrackingAction(con 48 void TrackingAction::PreUserTrackingAction(const G4Track*)
56 { 49 {
57 /* << 50 // G4cout << "ID= " << aTrack->GetTrackID() << " e(MeV)= "
58 //debug << 51 // << aTrack->GetDynamicParticle()->GetKineticEnergy()/MeV << " "
59 const G4DynamicParticle* dynamic = aTrack- << 52 // << aTrack->GetDynamicParticle()->GetDefinition()->GetParticleName()
60 G4double dynamCharge = dynamic->GetCharge( << 53 // << G4endl;
61 G4int occup = dynamic->GetTotalOc <<
62 G4double pdgMass = dynamic->GetParticl <<
63 G4double invarMass = dynamic->Get4Moment <<
64 G4double dynamMass = dynamic->GetMass(); <<
65 G4double dif1 = invarMass - pdgMass; <<
66 G4double dif2 = dynamMass - pdgMass; <<
67 <<
68 G4cout <<
69 << "\n Begin of track :" <<
70 << "\n charge= " << dynamCharge << " <<
71 << "\n pdgMass= " << G4BestUnit (pdgMa <<
72 /// << "\n invarMass= " << G4BestUnit (in <<
73 /// << " invar-pdg= " << G4BestUnit (di <<
74 << "\n dynamMass= " << G4BestUnit (dynam <<
75 << " dynam-pdg= " << G4BestUnit (dif2 <<
76 << G4endl; <<
77 */ <<
78 } 54 }
79 55
80 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
81 57
82 void TrackingAction::PostUserTrackingAction(co 58 void TrackingAction::PostUserTrackingAction(const G4Track* aTrack)
83 { 59 {
84 // increase nb of processed tracks << 60 //increase nb of processed tracks
85 // count nb of steps of this track << 61 //count nb of steps of this track
86 G4int nbSteps = aTrack->GetCurrentStepNumber << 62 G4int nbSteps = aTrack->GetCurrentStepNumber();
87 G4double Trleng = aTrack->GetTrackLength(); 63 G4double Trleng = aTrack->GetTrackLength();
88 << 64
89 Run* run = static_cast<Run*>(G4RunManager::G <<
90 <<
91 if (aTrack->GetDefinition()->GetPDGCharge() 65 if (aTrack->GetDefinition()->GetPDGCharge() == 0.) {
92 run->CountTraks0(1); << 66 runAction->CountTraks0(1);
93 run->CountSteps0(nbSteps); << 67 runAction->CountSteps0(nbSteps);
94 } <<
95 else { <<
96 run->CountTraks1(1); <<
97 run->CountSteps1(nbSteps); <<
98 } <<
99 <<
100 // true and projected ranges for primary par <<
101 if (aTrack->GetTrackID() == 1) { <<
102 run->AddTrueRange(Trleng); <<
103 G4ThreeVector vertex = fPrimary->GetPartic <<
104 G4ThreeVector position = aTrack->GetPositi <<
105 run->AddProjRange(position.x()); <<
106 run->AddTransvDev(position.y()); <<
107 run->AddTransvDev(position.z()); <<
108 <<
109 G4AnalysisManager* analysisManager = G4Ana <<
110 analysisManager->FillH1(1, Trleng); <<
111 analysisManager->FillH1(2, (float)nbSteps) <<
112 } <<
113 /* <<
114 //debug <<
115 const G4DynamicParticle* dynamic = aTrack- <<
116 G4double dynamCharge = dynamic->GetCharge( <<
117 G4int occup = dynamic->GetTotalOc <<
118 G4double pdgMass = dynamic->GetParticl <<
119 G4double invarMass = dynamic->Get4Moment <<
120 G4double dynamMass = dynamic->GetMass(); <<
121 G4double dif1 = invarMass - pdgMass; <<
122 G4double dif2 = dynamMass - pdgMass; <<
123 <<
124 G4cout <<
125 << "\n End of track :" <<
126 << "\n charge= " << dynamCharge << " <<
127 << "\n pdgMass= " << G4BestUnit (pdgMa <<
128 /// << "\n invarMass= " << G4BestUnit (in <<
129 /// << " invar-pdg= " << G4BestUnit (di <<
130 << "\n dynamMass= " << G4BestUnit (dynam <<
131 << " dynam-pdg= " << G4BestUnit (dif2 <<
132 << G4endl; <<
133 */ <<
134 68
135 // energy leakage << 69 } else {
136 G4StepStatus status = aTrack->GetStep()->Get << 70 runAction->CountTraks1(1);
137 if (status == fWorldBoundary) { << 71 runAction->CountSteps1(nbSteps);
138 G4double eleak = aTrack->GetKineticEnergy( <<
139 if (aTrack->GetDefinition() == G4Positron: <<
140 eleak += 2*electron_mass_c2; <<
141 } <<
142 fEventAction->AddEleak(eleak); <<
143 } 72 }
>> 73
>> 74 //true and projected ranges for primary particle
>> 75 if (aTrack->GetTrackID() == 1) {
>> 76 runAction->AddTrueRange(Trleng);
>> 77 G4ThreeVector vertex = primary->GetParticleGun()->GetParticlePosition();
>> 78 G4ThreeVector position = aTrack->GetPosition() - vertex;
>> 79 runAction->AddProjRange(position.x());
>> 80 runAction->AddTransvDev(position.y());
>> 81 runAction->AddTransvDev(position.z());
>> 82
>> 83 histoManager->FillHisto(1,Trleng);
>> 84 histoManager->FillHisto(2,(float)nbSteps);
>> 85 }
144 } 86 }
145 87
146 //....oooOO0OOooo........oooOO0OOooo........oo 88 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 89
147 90