| 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 //
27 /// \file Par02FastSimModelTracker.cc
28 /// \brief Implementation of the Par02FastSimM
29
30 #include "Par02FastSimModelTracker.hh"
31
32 #include "Par02EventInformation.hh"
33 #include "Par02Output.hh"
34 #include "Par02PrimaryParticleInformation.hh"
35 #include "Par02Smearer.hh"
36
37 #include "G4AnalysisManager.hh"
38 #include "G4Electron.hh"
39 #include "G4Event.hh"
40 #include "G4FieldTrack.hh"
41 #include "G4FieldTrackUpdator.hh"
42 #include "G4Gamma.hh"
43 #include "G4PathFinder.hh"
44 #include "G4Positron.hh"
45 #include "G4RunManager.hh"
46 #include "G4SystemOfUnits.hh"
47 #include "G4Track.hh"
48 #include "Randomize.hh"
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51
52 Par02FastSimModelTracker::Par02FastSimModelTra
53 G4String aModelName, G4Region* aEnvelope, Pa
54 : G4VFastSimulationModel(aModelName, aEnvelo
55 fCalculateParametrisation(),
56 fParametrisation(aType)
57 {}
58
59 //....oooOO0OOooo........oooOO0OOooo........oo
60
61 Par02FastSimModelTracker::Par02FastSimModelTra
62 : G4VFastSimulationModel(aModelName, aEnvelo
63 fCalculateParametrisation(),
64 fParametrisation(Par02DetectorParametrisat
65 {}
66
67 //....oooOO0OOooo........oooOO0OOooo........oo
68
69 Par02FastSimModelTracker::Par02FastSimModelTra
70 : G4VFastSimulationModel(aModelName),
71 fCalculateParametrisation(),
72 fParametrisation(Par02DetectorParametrisat
73 {}
74
75 //....oooOO0OOooo........oooOO0OOooo........oo
76
77 Par02FastSimModelTracker::~Par02FastSimModelTr
78
79 //....oooOO0OOooo........oooOO0OOooo........oo
80
81 G4bool Par02FastSimModelTracker::IsApplicable(
82 {
83 return aParticleType.GetPDGCharge() != 0; /
84 }
85
86 //....oooOO0OOooo........oooOO0OOooo........oo
87
88 G4bool Par02FastSimModelTracker::ModelTrigger(
89 {
90 return true; // No kinematical restrictions
91 }
92
93 //....oooOO0OOooo........oooOO0OOooo........oo
94
95 void Par02FastSimModelTracker::DoIt(const G4Fa
96 {
97 G4cout << " ________Tracker model triggered
98
99 // Calculate the final position (at the oute
100 // of the particle with the momentum at the
101
102 G4Track track = *aFastTrack.GetPrimaryTrack(
103 G4FieldTrack aFieldTrack('0');
104 G4FieldTrackUpdator::Update(&aFieldTrack, &t
105
106 G4double retSafety = -1.0;
107 ELimited retStepLimited;
108 G4FieldTrack endTrack('a');
109 G4double currentMinimumStep = 10.0 * m; //
110 //
111 G4PathFinder* fPathFinder = G4PathFinder::Ge
112 /*G4double lengthAlongCurve = */
113 fPathFinder->ComputeStep(aFieldTrack, curren
114 aFastTrack.GetPrima
115 retStepLimited, end
116
117 // Place the particle at the tracking detect
118 // (at the place it would reach without the
119 aFastStep.ProposePrimaryTrackFinalPosition(e
120
121 // Consider only primary tracks (do nothing
122 G4ThreeVector Porg = aFastTrack.GetPrimaryTr
123 if (!aFastTrack.GetPrimaryTrack()->GetParent
124 auto info = (Par02EventInformation*)G4Even
125 if (info->GetDoSmearing()) {
126 // Smearing according to the tracking de
127 G4double res = fCalculateParametrisation
128 Par02DetectorParametrisation::eTRACKER
129 G4double eff = fCalculateParametrisation
130 Par02DetectorParametrisation::eTRACKER
131 G4ThreeVector Psm;
132 Psm = Par02Smearer::Instance()->SmearMom
133 Par02Output::Instance()->FillHistogram(0
134 // Setting the values of Psm, res and ef
135 ((Par02PrimaryParticleInformation*)(cons
136 aF
137
138
139 ->
140 ->SetTrackerMomentum(Psm);
141 ((Par02PrimaryParticleInformation*)(cons
142 aF
143
144
145 ->
146 ->SetTrackerResolution(res);
147 ((Par02PrimaryParticleInformation*)(cons
148 aF
149
150
151 ->
152 ->SetTrackerEfficiency(eff);
153 }
154 else {
155 // No smearing: simply setting the value
156 ((Par02PrimaryParticleInformation*)(cons
157 aF
158
159
160 ->
161 ->SetTrackerMomentum(Porg);
162 }
163 }
164 }
165
166 //....oooOO0OOooo........oooOO0OOooo........oo
167