| 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 /// \file Par01/src/Par01PiModel.cc
27 /// \brief Implementation of the Par01PiModel
28 //
29 //
30 //
31 #include "Par01PiModel.hh"
32
33 #include "G4Gamma.hh"
34 #include "G4PionMinus.hh"
35 #include "G4PionPlus.hh"
36
37 //....oooOO0OOooo........oooOO0OOooo........oo
38
39 Par01PiModel::Par01PiModel(G4Region* anEnvelop
40 : G4VFastSimulationModel("Par01PiModel", anE
41 {
42 ;
43 }
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 Par01PiModel::~Par01PiModel()
48 {
49 ;
50 }
51
52 //....oooOO0OOooo........oooOO0OOooo........oo
53
54 G4bool Par01PiModel::IsApplicable(const G4Part
55 {
56 return &particleType == G4PionMinus::PionMin
57 || &particleType == G4PionPlus::PionP
58 }
59
60 //....oooOO0OOooo........oooOO0OOooo........oo
61
62 G4bool Par01PiModel::ModelTrigger(const G4Fast
63 {
64 //------------------------------------------
65 // UserTrigger() method: method which has to
66 // the parameterisation has to be applied.
67 // Here ModelTrigger() asks the user (ie you
68 //
69 // Note that quantities like the local/globa
70 // are available at this level via the fastT
71 // to check distance from boundaries, see be
72 //------------------------------------------
73 G4cout << "\nPar01PiModel::ModelTrigger() ca
74 G4cout << "--------------------------------"
75 G4cout << "(particle is a " << fastTrack.Get
76 << " )\n"
77 << G4endl;
78
79 // -- Examples of available informations:
80
81 // -- position:
82 G4cout << " Track position: " << fas
83 << "(global coord.)" << fastTrack.Get
84 << G4endl;
85
86 // -- direction:
87 G4cout << " Track direction:" << fas
88 << "(global coord.)" << fastTrack.Get
89 << G4endl;
90
91 return true;
92 }
93
94 //....oooOO0OOooo........oooOO0OOooo........oo
95
96 void Par01PiModel::DoIt(const G4FastTrack& fas
97 //--------------------------------------------
98 //
99 // User method to code the parameterisation pr
100 //
101 //--------------------------------------------
102 {
103 //------------------------------------------
104 // The primary track continues along its dir
105 // One secondary (a photon) is added:
106 //------------------------------------------
107 G4cout << " Pion `model' applied\n" <<
108
109 //------------------------------
110 // Primary:
111 // idem as in "DefaultModel":
112 //
113 //------------------------------
114 G4ThreeVector position;
115 G4double distance;
116 distance = fastTrack.GetEnvelopeSolid()->Dis
117
118 position =
119 fastTrack.GetPrimaryTrackLocalPosition() +
120
121 // -- set final position:
122 fastStep.ProposePrimaryTrackFinalPosition(po
123
124 //---------------------------
125 // Secondary:
126 // Adds one "secondary":
127 //---------------------------
128 // -- First, user has to say how many second
129 fastStep.SetNumberOfSecondaryTracks(1);
130
131 //------------------------
132 // -- Build the secondary:
133 //------------------------
134 // -- direction:
135 G4ParticleMomentum direction(fastTrack.GetPr
136 direction.setZ(direction.z() * 0.5);
137 direction.setY(direction.y() + direction.z()
138 direction = direction.unit(); // necessary
139
140 // -- dynamics (Note that many constructors
141 // -- see prototype/particle+matter/particle
142 G4DynamicParticle dynamique(G4Gamma::GammaDe
143 fastTrack.GetPri
144 // -- position:
145 G4double Dist;
146 Dist = fastTrack.GetEnvelopeSolid()->Distanc
147
148 G4ThreeVector posi;
149 posi = fastTrack.GetPrimaryTrackLocalPositio
150
151 //------------------------------------
152 //-- Creation of the secondary Track:
153 //------------------------------------
154 fastStep.CreateSecondaryTrack(dynamique, pos
155 }
156