| 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 Par02FastSimModelEMCal.cc
28 /// \brief Implementation of the Par02FastSimM
29
30 #include "Par02FastSimModelEMCal.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 "G4Gamma.hh"
41 #include "G4Positron.hh"
42 #include "G4RunManager.hh"
43 #include "G4SystemOfUnits.hh"
44 #include "G4Track.hh"
45 #include "Randomize.hh"
46
47 //....oooOO0OOooo........oooOO0OOooo........oo
48
49 Par02FastSimModelEMCal::Par02FastSimModelEMCal
50
51 : G4VFastSimulationModel(aModelName, aEnvelo
52 fCalculateParametrisation(),
53 fParametrisation(aType)
54 {}
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 Par02FastSimModelEMCal::Par02FastSimModelEMCal
59 : G4VFastSimulationModel(aModelName, aEnvelo
60 fCalculateParametrisation(),
61 fParametrisation(Par02DetectorParametrisat
62 {}
63
64 //....oooOO0OOooo........oooOO0OOooo........oo
65
66 Par02FastSimModelEMCal::Par02FastSimModelEMCal
67 : G4VFastSimulationModel(aModelName),
68 fCalculateParametrisation(),
69 fParametrisation(Par02DetectorParametrisat
70 {}
71
72 //....oooOO0OOooo........oooOO0OOooo........oo
73
74 Par02FastSimModelEMCal::~Par02FastSimModelEMCa
75
76 //....oooOO0OOooo........oooOO0OOooo........oo
77
78 G4bool Par02FastSimModelEMCal::IsApplicable(co
79 {
80 // Applicable for electrons, positrons, and
81 return &aParticleType == G4Electron::Definit
82 || &aParticleType == G4Gamma::Definit
83 }
84
85 //....oooOO0OOooo........oooOO0OOooo........oo
86
87 G4bool Par02FastSimModelEMCal::ModelTrigger(co
88 {
89 return true; // No kinematical restrictions
90 }
91
92 //....oooOO0OOooo........oooOO0OOooo........oo
93
94 void Par02FastSimModelEMCal::DoIt(const G4Fast
95 {
96 // G4cout << " ________EMCal model triggered
97
98 // Kill the parameterised particle at the en
99 aFastStep.KillPrimaryTrack();
100 aFastStep.ProposePrimaryTrackPathLength(0.0)
101 G4double Edep = aFastTrack.GetPrimaryTrack()
102
103 // Consider only primary tracks (do nothing
104 G4ThreeVector Pos = aFastTrack.GetPrimaryTra
105 if (!aFastTrack.GetPrimaryTrack()->GetParent
106 auto info = (Par02EventInformation*)G4Even
107 if (info->GetDoSmearing()) {
108 // Smearing according to the electromagn
109 G4ThreeVector Porg = aFastTrack.GetPrima
110 G4double res = fCalculateParametrisation
111
112 G4double eff = fCalculateParametrisation
113
114 G4double Esm;
115 Esm = std::abs(Par02Smearer::Instance()-
116 Par02Output::Instance()->FillHistogram(1
117 // Setting the values of Pos, Esm, res a
118 ((Par02PrimaryParticleInformation*)(cons
119 aF
120
121
122 ->
123 ->SetEMCalPosition(Pos);
124 ((Par02PrimaryParticleInformation*)(cons
125 aF
126
127
128 ->
129 ->SetEMCalEnergy(Esm);
130 ((Par02PrimaryParticleInformation*)(cons
131 aF
132
133
134 ->
135 ->SetEMCalResolution(res);
136 ((Par02PrimaryParticleInformation*)(cons
137 aF
138
139
140 ->
141 ->SetEMCalEfficiency(eff);
142 // The (smeared) energy of the particle
143 // (which corresponds to the entrance of
144 aFastStep.ProposeTotalEnergyDeposited(Es
145 }
146 else {
147 // No smearing: simply setting the value
148 ((Par02PrimaryParticleInformation*)(cons
149 aF
150
151
152 ->
153 ->SetEMCalEnergy(Edep);
154 // The (initial) energy of the particle
155 // (which corresponds to the entrance of
156 aFastStep.ProposeTotalEnergyDeposited(Ed
157 }
158 }
159 }
160
161 //....oooOO0OOooo........oooOO0OOooo........oo
162