source: trunk/source/processes/hadronic/util/include/G4ReactionDynamics.hh@ 1196

Last change on this file since 1196 was 819, checked in by garnier, 17 years ago

import all except CVS
File size: 6.8 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27//
28 // Hadronic Process: Reaction Dynamics
29 // original by H.P. Wellisch
30 // Modified by J.L.Chuma 19-Nov-96
31 // Modified by J.L.Chuma 27-Mar-97
32 // Modified by J.L.Chuma 30-Apr-97
33 // Modified by J.L.Chuma 06-Aug-97 to include the original incident particle
34 // before Fermi motion and evaporation effects
35
36#ifndef G4ReactionDynamics_h
37#define G4ReactionDynamics_h 1
38
39#include "G4ParticleTypes.hh"
40#include "G4DynamicParticle.hh"
41#include "G4ReactionProduct.hh"
42#include "G4Nucleus.hh"
43#include "G4FastVector.hh"
44#include "G4HadProjectile.hh"
45
46enum{ GHADLISTSIZE=256};
47
48 class G4ReactionDynamics
49 {
50 public:
51
52 G4ReactionDynamics() {}
53
54 virtual ~G4ReactionDynamics() {}
55
56 virtual G4double FindInelasticity()
57 { return 0.0; }
58
59 virtual G4double FindTimeDelay()
60 { return 0.0; }
61
62 G4bool GenerateXandPt( // derived from GENXPT
63 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
64 G4int &vecLen,
65 G4ReactionProduct &modifiedOriginal, // Fermi motion & evap. effect included
66 const G4HadProjectile *originalIncident,
67 G4ReactionProduct &currentParticle,
68 G4ReactionProduct &targetParticle,
69 const G4DynamicParticle* originalTarget,
70 const G4Nucleus &targetNucleus,
71 G4bool &incidentHasChanged,
72 G4bool &targetHasChanged,
73 G4bool leadFlag,
74 G4ReactionProduct &leadingStrangeParticle );
75
76 void SuppressChargedPions(
77 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
78 G4int &vecLen,
79 const G4ReactionProduct &modifiedOriginal,
80 G4ReactionProduct &currentParticle,
81 G4ReactionProduct &targetParticle,
82 const G4Nucleus &targetNucleus,
83 G4bool &incidentHasChanged,
84 G4bool &targetHasChanged );
85
86 G4bool TwoCluster( // derived from TWOCLU
87 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
88 G4int &vecLen,
89 G4ReactionProduct &modifiedOriginal, // Fermi motion & evap. effect included
90 const G4HadProjectile *originalIncident,
91 G4ReactionProduct &currentParticle,
92 G4ReactionProduct &targetParticle,
93 const G4DynamicParticle* originalTarget,
94 const G4Nucleus &targetNucleus,
95 G4bool &incidentHasChanged,
96 G4bool &targetHasChanged,
97 G4bool leadFlag,
98 G4ReactionProduct &leadingStrangeParticle );
99
100 void TwoBody( // derived from TWOB
101 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
102 G4int &vecLen,
103 G4ReactionProduct &modifiedOriginal,
104 const G4DynamicParticle *originalTarget,
105 G4ReactionProduct &currentParticle,
106 G4ReactionProduct &targetParticle,
107 const G4Nucleus &targetNucleus,
108 G4bool &targetHasChanged );
109
110 G4int Factorial( G4int n );
111
112 G4double GenerateNBodyEvent( // derived from PHASP
113 const G4double totalEnergy,
114 const G4bool constantCrossSection,
115 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
116 G4int &vecLen );
117
118 void ProduceStrangeParticlePairs(
119 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
120 G4int &vecLen,
121 const G4ReactionProduct &modifiedOriginal,
122 const G4DynamicParticle *originalTarget,
123 G4ReactionProduct &currentParticle,
124 G4ReactionProduct &targetParticle,
125 G4bool &incidentHasChanged,
126 G4bool &targetHasChanged );
127
128 void NuclearReaction( // derived from NUCREC
129 G4FastVector<G4ReactionProduct,4> &vec,
130 G4int &vecLen,
131 const G4HadProjectile *originalIncident,
132 const G4Nucleus &aNucleus,
133 const G4double theAtomicMass,
134 const G4double *massVec );
135
136 private:
137
138 void Rotate(
139 const G4double numberofFinalStateNucleons,
140 const G4ThreeVector &temp,
141 const G4ReactionProduct &modifiedOriginal, // Fermi motion & evap. effect included
142 const G4HadProjectile *originalIncident,
143 const G4Nucleus &targetNucleus,
144 G4ReactionProduct &currentParticle,
145 G4ReactionProduct &targetParticle,
146 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
147 G4int &vecLen );
148
149 void Defs1(
150 const G4ReactionProduct &modifiedOriginal,
151 G4ReactionProduct &currentParticle,
152 G4ReactionProduct &targetParticle,
153 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
154 G4int &vecLen );
155
156 void AddBlackTrackParticles(
157 const G4double epnb,
158 const G4int npnb,
159 const G4double edta,
160 const G4int ndta,
161 const G4double sprob,
162 const G4double kineticMinimum,
163 const G4double kineticFactor,
164 const G4ReactionProduct &modifiedOriginal,
165 G4int PinNucleus,
166 G4int NinNucleus,
167 const G4Nucleus &aNucleus,
168 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
169 G4int &vecLen );
170
171 std::pair<G4int, G4int> GetFinalStateNucleons(
172 const G4DynamicParticle* originalTarget,
173 const G4FastVector<G4ReactionProduct,GHADLISTSIZE>& vec,
174 const G4int& vecLen );
175
176 void MomentumCheck(
177 const G4ReactionProduct &modifiedOriginal,
178 G4ReactionProduct &currentParticle,
179 G4ReactionProduct &targetParticle,
180 G4FastVector<G4ReactionProduct,GHADLISTSIZE> &vec,
181 G4int &vecLen );
182
183 G4double normal();
184
185 G4int Poisson( G4double x );
186
187 };
188
189#endif
190
Note: See TracBrowser for help on using the repository browser.