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 | // $Id: G4Incl.hh,v 1.10 2007/12/03 19:36:06 miheikki Exp $ |
---|
27 | // Translation of INCL4.2/ABLA V3 |
---|
28 | // Pekka Kaitaniemi, HIP (translation) |
---|
29 | // Christelle Schmidt, IPNL (fission code) |
---|
30 | // Alain Boudard, CEA (contact person INCL/ABLA) |
---|
31 | // Aatos Heikkinen, HIP (project coordination) |
---|
32 | |
---|
33 | #ifndef G4Incl_hh |
---|
34 | #define G4Incl_hh 1 |
---|
35 | |
---|
36 | #include "globals.hh" |
---|
37 | #include "G4InclDataDefs.hh" |
---|
38 | #include "G4Abla.hh" |
---|
39 | #include <fstream> |
---|
40 | |
---|
41 | using namespace std; |
---|
42 | /** |
---|
43 | * Class containing INCL4 hadronic cascade algorithm. |
---|
44 | */ |
---|
45 | |
---|
46 | class G4Incl { |
---|
47 | |
---|
48 | public: |
---|
49 | /** |
---|
50 | * |
---|
51 | * Support for Doxygen JavaDoc style. |
---|
52 | * |
---|
53 | * \author{pekka.kaitaniemi@helsinki.fi} |
---|
54 | */ |
---|
55 | |
---|
56 | /** |
---|
57 | * Constructor to be used with Geant4. |
---|
58 | * @param hazard a pointer to G4Hazard structure. |
---|
59 | * @param calincl a pointer to G4Calincl structure. |
---|
60 | * @param ws a pointer to G4Ws structure. |
---|
61 | * @param mat a pointer to G4Mat structure. |
---|
62 | * @param varntp a pointer to G4VarNtp structure. |
---|
63 | */ |
---|
64 | G4Incl(G4Hazard *hazard, G4Calincl *calincl, G4Ws *ws, G4Mat *mat, G4VarNtp *varntp); |
---|
65 | |
---|
66 | /** |
---|
67 | * Constructor for private unit testing purposes. |
---|
68 | * @param hazard a pointer to G4Hazard structure. |
---|
69 | * @param dton a pointer to G4Dton structure. |
---|
70 | * @param saxw a pointer to G4Saxw structure. |
---|
71 | * @param ws a pointer to G4Ws structure. |
---|
72 | */ |
---|
73 | G4Incl(G4Hazard *hazard, G4Dton *dton, G4Saxw *saxw, G4Ws *ws); |
---|
74 | |
---|
75 | G4Incl(); |
---|
76 | |
---|
77 | ~G4Incl(); // Destructor |
---|
78 | |
---|
79 | G4double energyTest(G4int i); // Test for NaN energy of particle i. |
---|
80 | void dumpBl5(std::ofstream& dumpOut); // Dump the contents of G4Bl5. |
---|
81 | void dumpSaxw(std::ofstream& dumpOut); // Dump the contents of G4Saxw. |
---|
82 | void dumpBl1(std::ofstream& dumpOut); // Dump the contents of G4Bl1. |
---|
83 | void dumpBl2(std::ofstream& dumpOut); // Dump the contents of G4Bl2. |
---|
84 | void dumpBl3(std::ofstream& dumpOut); // Dump the contents of G4Bl3. |
---|
85 | |
---|
86 | /** |
---|
87 | * Set verbosity level. |
---|
88 | */ |
---|
89 | void setVerboseLevel(G4int level); |
---|
90 | /** |
---|
91 | * Get verbosity level. |
---|
92 | */ |
---|
93 | G4int getVerboseLevel(); |
---|
94 | |
---|
95 | void setDtonData(G4Dton *newDton); // Set internal data. |
---|
96 | void setWsData(G4Ws *newWs); |
---|
97 | void setHazardData(G4Hazard *newHazard); |
---|
98 | void setSaxwData(G4Saxw *newSaxw); |
---|
99 | void setSpl2Data(G4Spl2 *newSpl2); |
---|
100 | void setMatData(G4Mat *newMat); |
---|
101 | void setCalinclData(G4Calincl *newCalincl); |
---|
102 | void setLightNucData(G4LightNuc *newLightNuc); |
---|
103 | void setLightGausNucData(G4LightGausNuc *newLightGausNuc); |
---|
104 | void setBl1Data(G4Bl1 *newBl1); |
---|
105 | void setBl2Data(G4Bl2 *newBl2); |
---|
106 | void setBl3Data(G4Bl3 *newBl3); |
---|
107 | void setBl4Data(G4Bl4 *newBl4); |
---|
108 | void setBl5Data(G4Bl5 *newBl5); |
---|
109 | void setBl6Data(G4Bl6 *newBl6); |
---|
110 | void setBl8Data(G4Bl8 *newBl8); |
---|
111 | void setBl9Data(G4Bl9 *newBl9); |
---|
112 | void setBl10Data(G4Bl10 *newBl10); |
---|
113 | void setKindData(G4Kind *newKind); |
---|
114 | |
---|
115 | public: |
---|
116 | /** |
---|
117 | * Process one event with INCL4 only. |
---|
118 | */ |
---|
119 | void processEventIncl(); |
---|
120 | |
---|
121 | /** |
---|
122 | * Process one event with INCL4 and built-in ABLA evaporation and fission. |
---|
123 | */ |
---|
124 | void processEventInclAbla(G4int eventnumber); |
---|
125 | |
---|
126 | public: // Methods used to initialize INCL |
---|
127 | /** |
---|
128 | * Initialize INCL. |
---|
129 | * |
---|
130 | * @param initRandomSeed choose whether INCL should initialize random seeds. |
---|
131 | */ |
---|
132 | void initIncl(G4bool initRandomSeed); |
---|
133 | |
---|
134 | /** |
---|
135 | * Initialize target materials. |
---|
136 | * |
---|
137 | * @param izmat charge number |
---|
138 | * @param iamat mass number |
---|
139 | * @param imat material number (array index) |
---|
140 | */ |
---|
141 | void initMaterial(G4int izmat, G4int iamat, G4int imat); |
---|
142 | |
---|
143 | /** |
---|
144 | * A normal member taking two arguments and returning an integer value. |
---|
145 | * |
---|
146 | * @param l an integer argument. |
---|
147 | * @param q a constant character pointer. |
---|
148 | * @return The test results |
---|
149 | */ |
---|
150 | G4double deutv(G4int l, G4double q); |
---|
151 | |
---|
152 | /** |
---|
153 | * Returns the values of the function: |
---|
154 | * \f[ |
---|
155 | * (0.23162461 + (j - 1))^2 |
---|
156 | * \f] |
---|
157 | * |
---|
158 | * @param j an integer parameter |
---|
159 | * @return a double value |
---|
160 | */ |
---|
161 | G4double fm2(G4int j); |
---|
162 | |
---|
163 | /** |
---|
164 | * Interpolates function described by class G4Saxw around a point. |
---|
165 | * |
---|
166 | * @param xv interpolation point |
---|
167 | * @return a double value |
---|
168 | */ |
---|
169 | G4double interpolateFunction(G4double xv); |
---|
170 | |
---|
171 | /** |
---|
172 | * Calculates the first derivative of the function stored in class G4Saxw. |
---|
173 | * |
---|
174 | * @param xv an integer parameter |
---|
175 | */ |
---|
176 | void firstDerivative(G4int k); |
---|
177 | |
---|
178 | /** |
---|
179 | * Returns the values of the function: |
---|
180 | * \f[ |
---|
181 | * \frac{r^2}{1 + e^{\frac{r - r_{0}}{A_{dif}}}} |
---|
182 | * \f] |
---|
183 | * |
---|
184 | * @param r a G4double argument |
---|
185 | * @return a double value |
---|
186 | */ |
---|
187 | G4double wsax(G4double r); |
---|
188 | |
---|
189 | /** |
---|
190 | * Returns the values of the function: |
---|
191 | * \f[ |
---|
192 | * \frac{1}{A_{dif}^{2}} \frac{r^3 e^{r - r_0}}{1 + e^{\frac{r - r_{0}}{A_{dif}}}} |
---|
193 | * \f] |
---|
194 | * |
---|
195 | * @param r a G4double argument |
---|
196 | * @return a double value |
---|
197 | */ |
---|
198 | G4double derivWsax(G4double r); |
---|
199 | |
---|
200 | /** |
---|
201 | * Returns the values of the function: |
---|
202 | * \f[ |
---|
203 | * r^2 (1.0 + r_0 \frac{r^2}{A_{dif}}) e^{-\frac{r^2}{A_{dif}^2}} |
---|
204 | * \f] |
---|
205 | * @param r a G4double argument |
---|
206 | * @return a double value |
---|
207 | */ |
---|
208 | G4double dmho(G4double r); |
---|
209 | |
---|
210 | /** |
---|
211 | * Returns the values of the function: |
---|
212 | * \f[ |
---|
213 | * -\frac{2r^4}{A_{dif}^2} (r_0(1.0 - \frac{r^2}{A_{dif}^2}) - 1.0) e^{-\frac{r^2}{A_{dif}^2}} |
---|
214 | * \f] |
---|
215 | * |
---|
216 | * @param r a G4double argument |
---|
217 | * @return a double value |
---|
218 | */ |
---|
219 | G4double derivMho(G4double r); |
---|
220 | |
---|
221 | /** |
---|
222 | * Returns the values of the function: |
---|
223 | * \f[ |
---|
224 | * \frac{r^4}{A_{dif}^2} e^{-\frac{1}{2} \frac{r^2}{A_{dif}^2}} |
---|
225 | * \f] |
---|
226 | * |
---|
227 | * @param r a G4double argument |
---|
228 | * @return a double value |
---|
229 | */ |
---|
230 | G4double derivGaus(G4double r); |
---|
231 | |
---|
232 | /** |
---|
233 | * Ce subroutine appele sur le premier tir va calculer la densite du deuton |
---|
234 | * dans l'espace des impulsions et preparer l'interpolation permettant ensuite |
---|
235 | * le tir au hasard d'un module de l'impulsion (q). |
---|
236 | * Ce subroutine remplit le G4Spl2: |
---|
237 | * xsp, ysp integrale normalisee de la densite de 0 a q. |
---|
238 | * a(),b(),c() coefs des nsp points pour une interpolation du second degre. |
---|
239 | * q est en fm-1. |
---|
240 | */ |
---|
241 | void densDeut(); |
---|
242 | |
---|
243 | /** |
---|
244 | * Integrate using Alkazhov's method. |
---|
245 | * |
---|
246 | * @param ami a double parameter |
---|
247 | * @param ama a double parameter |
---|
248 | * @param dr a double parameter |
---|
249 | * @param functionChoice an integer parameter |
---|
250 | * @return a double value |
---|
251 | */ |
---|
252 | G4double integrate(G4double ami, G4double ama, G4double step, G4int functionChoice); |
---|
253 | |
---|
254 | /** |
---|
255 | * Deuteron density |
---|
256 | * |
---|
257 | * @param q a double parameter |
---|
258 | * @return a double value |
---|
259 | */ |
---|
260 | G4double dens(G4double q); |
---|
261 | |
---|
262 | /** |
---|
263 | * |
---|
264 | */ |
---|
265 | void spl2ab(); |
---|
266 | |
---|
267 | /** |
---|
268 | * |
---|
269 | * @param xv a double parameter |
---|
270 | * @return a double value |
---|
271 | */ |
---|
272 | G4double splineab(G4double xv); |
---|
273 | |
---|
274 | public: // Main INCL routines |
---|
275 | /** |
---|
276 | * INCL model as a function. |
---|
277 | */ |
---|
278 | void pnu(G4int *ibert_p, G4int *nopart_p, G4int *izrem_p, G4int *iarem_p, G4double *esrem_p, |
---|
279 | G4double *erecrem_p, G4double *alrem_p, G4double *berem_p, G4double *garem_p, |
---|
280 | G4double *bimpact_p, G4int *l_p); |
---|
281 | |
---|
282 | /** |
---|
283 | * Single nucleon-nucleon collision. |
---|
284 | */ |
---|
285 | void collis(G4double *p1_p, G4double *p2_p, G4double *p3_p, G4double *e1_p, G4double *pout11_p, G4double *pout12_p, |
---|
286 | G4double *pout13_p, G4double *eout1_p, G4double *q1_p, G4double *q2_p, G4double *q3_p, |
---|
287 | G4double *q4_p, G4int *np_p, G4int *ip_p, G4int *k2_p, G4int *k3_p, G4int *k4_p, |
---|
288 | G4int *k5_p, G4int *m1_p, G4int *m2_p, G4int *is1_p, G4int *is2_p); |
---|
289 | |
---|
290 | /** |
---|
291 | * This routine describes the anisotropic decay of a particle of |
---|
292 | * mass xi into 2 particles of masses x1,x2. |
---|
293 | * The anisotropy is supposed to follow a 1+3*hel*(std::cos(theta))**2 law |
---|
294 | * with respect to the direction of the incoming particle. |
---|
295 | * |
---|
296 | * In the input, p1,p2,p3 is the momentum of particle xi. |
---|
297 | * |
---|
298 | * In the output, p1,p2,p3 is the momentum of particle x1, |
---|
299 | * while q1,q2,q3 is the momentum of particle x2. |
---|
300 | * |
---|
301 | * @param p1_p pointer to momentum component 1 |
---|
302 | * @param p2_p pointer to momentum component 2 |
---|
303 | * @param p3_p pointer to momentum component 3 |
---|
304 | * @param wp_p pointer to a double parameter |
---|
305 | * @param q1_p pointer to momentum component 1 |
---|
306 | * @param q2_p pointer to momentum component 2 |
---|
307 | * @param q3_p pointer to momentum component 3 |
---|
308 | * @param wq_p pointer to a double parameter |
---|
309 | * @param xi_p pointer to a double parameter |
---|
310 | * @param x1_p pointer to momentum component 1 |
---|
311 | * @param x2_p pointer to momentum component 2 |
---|
312 | * @param x3_p pointer to momentum compone n3 |
---|
313 | * @param hel_p pointer to a double parameter |
---|
314 | */ |
---|
315 | void decay2(G4double *p1_p, G4double *p2_p, G4double *p3_p, G4double *wp_p, G4double *q1_p, |
---|
316 | G4double *q2_p, G4double *q3_p, G4double *wq_p, G4double *xi_p, G4double *x1_p, |
---|
317 | G4double *x2_p, G4double *hel_p); |
---|
318 | |
---|
319 | /** |
---|
320 | * Time calculation. |
---|
321 | * |
---|
322 | * @param i an index of particle 1 |
---|
323 | * @param j an index of particle 2 |
---|
324 | */ |
---|
325 | void time(G4int i, G4int j); |
---|
326 | |
---|
327 | /** |
---|
328 | * New time. |
---|
329 | * |
---|
330 | * @param l1 an integer parameter |
---|
331 | * @param l2 an integer parameter |
---|
332 | */ |
---|
333 | void newt(G4int l1, G4int l2); |
---|
334 | |
---|
335 | /** |
---|
336 | * |
---|
337 | * |
---|
338 | * @param l1 an integer parameter |
---|
339 | */ |
---|
340 | void new1(G4int l1); |
---|
341 | |
---|
342 | /** |
---|
343 | * |
---|
344 | * |
---|
345 | * @param y1 a double parameter |
---|
346 | * @param y2 a double parameter |
---|
347 | * @param y3 a double parameter |
---|
348 | * @param q1 a double parameter |
---|
349 | * @param q2 a double parameter |
---|
350 | * @param q3 a double parameter |
---|
351 | * @param q4 a double parameter |
---|
352 | * @param npion number of pions |
---|
353 | * @param l1 an integer parameter |
---|
354 | */ |
---|
355 | void new2(G4double y1, G4double y2, G4double y3, G4double q1, G4double q2, G4double q3, |
---|
356 | G4double q4, G4int npion, G4int l1); |
---|
357 | |
---|
358 | /** |
---|
359 | * |
---|
360 | * |
---|
361 | * @param y1 a double parameter |
---|
362 | * @param y2 a double parameter |
---|
363 | * @param y3 a double parameter |
---|
364 | * @param q1 a double parameter |
---|
365 | * @param q2 a double parameter |
---|
366 | * @param q3 a double parameter |
---|
367 | * @param q4 a double parameter |
---|
368 | * @param npion number of pions |
---|
369 | * @param l1 an integer parameter |
---|
370 | */ |
---|
371 | void new3(G4double y1, G4double y2, G4double y3, G4double q1, G4double q2, G4double q3, |
---|
372 | G4double q4, G4int npion, G4int l1); |
---|
373 | |
---|
374 | /** |
---|
375 | * Lorentz transformation. |
---|
376 | * |
---|
377 | * @param q1 a double parameter |
---|
378 | * @param q2 a double parameter |
---|
379 | * @param q3 a double parameter |
---|
380 | * @param b1 a double parameter |
---|
381 | * @param b2 a double parameter |
---|
382 | * @param b3 a double parameter |
---|
383 | * @param E energy |
---|
384 | */ |
---|
385 | void loren(G4double *q1, G4double *q2, G4double *q3, G4double *b1, G4double *b2, G4double *b3, G4double *E); |
---|
386 | |
---|
387 | /** |
---|
388 | * Pauli blocking. |
---|
389 | * |
---|
390 | * @param l an integer parameter |
---|
391 | * @param xr a double parameter |
---|
392 | * @param pr a double parameter |
---|
393 | * @return a double value |
---|
394 | */ |
---|
395 | G4double pauliBlocking(G4int l, G4double xr, G4double pr); |
---|
396 | |
---|
397 | /** |
---|
398 | * Fit by J. Vandermeulen. |
---|
399 | * Low energy fit from reference J.Cugnon, D. L'hote and J. Vandermeulen, NIM B111 (1996) 215. |
---|
400 | * |
---|
401 | * @param E energy |
---|
402 | * @param m a double parameter m = 0, 1, 2 for nucleon-nucleon, nucleon-delta, delta-delta |
---|
403 | * @param i a double parameter i = 2, 0, -2 for pp, pn, nn |
---|
404 | * @return a double value |
---|
405 | */ |
---|
406 | G4double lowEnergy(G4double E, G4double m, G4double i); |
---|
407 | |
---|
408 | /** |
---|
409 | * Total cross-sections. |
---|
410 | * |
---|
411 | * @param E energy |
---|
412 | * @param m an integer parameter m=0,1,2 for nucleon-nucleon, nucleon-delta, delta-delta |
---|
413 | * @param i an integer parameter i = 2, 0, -2 for pp, pn, nn |
---|
414 | * @return a double value |
---|
415 | */ |
---|
416 | G4double totalCrossSection(G4double E, G4int m, G4int i); |
---|
417 | |
---|
418 | /** |
---|
419 | * |
---|
420 | * |
---|
421 | * @param Ein energy |
---|
422 | * @param d a double parameter |
---|
423 | * @param i an integer parameter |
---|
424 | * @param isa an integer parameter |
---|
425 | * @return a double value |
---|
426 | */ |
---|
427 | G4double srec(G4double Ein, G4double d, G4int i, G4int isa); |
---|
428 | |
---|
429 | /** |
---|
430 | * Delta production cross section. |
---|
431 | * |
---|
432 | * @param E energy |
---|
433 | * @param i an integer parameter |
---|
434 | * @return a double value |
---|
435 | */ |
---|
436 | G4double deltaProductionCrossSection(G4double E, G4int i); |
---|
437 | |
---|
438 | /** |
---|
439 | * Sigma(pi+ + p) in the (3,3) region. |
---|
440 | * New fit by J. Vandermeulen and constant value above the (3,3) resonance. |
---|
441 | * |
---|
442 | * @param x a double parameter |
---|
443 | * @return a double value |
---|
444 | */ |
---|
445 | G4double pionNucleonCrossSection(G4double x); |
---|
446 | |
---|
447 | /** |
---|
448 | * Transmission probability for a nucleon of kinetic energy |
---|
449 | * E on the edge of the well of depth v0 (nr approximation). |
---|
450 | * , |
---|
451 | * of the nucleus and r is the target radius |
---|
452 | * |
---|
453 | * @param E kinetic energy |
---|
454 | * @param iz the isospin of the nucleon |
---|
455 | * @param izn the instanteneous charge |
---|
456 | * @param r |
---|
457 | * @param v0 |
---|
458 | * @return a double value |
---|
459 | */ |
---|
460 | G4double transmissionProb(G4double E, G4double iz, G4double izn, G4double r, G4double v0); |
---|
461 | |
---|
462 | /** |
---|
463 | * |
---|
464 | * |
---|
465 | * @param x1 a double parameter |
---|
466 | * @param x2 a double parameter |
---|
467 | * @param x3 a double parameter |
---|
468 | * @param p1 a double parameter |
---|
469 | * @param p2 a double parameter |
---|
470 | * @param p3 a double parameter |
---|
471 | * @param E a double parameter |
---|
472 | * @param r2 a double parameter |
---|
473 | * @return a double value |
---|
474 | */ |
---|
475 | G4double ref(G4double x1, G4double x2, G4double x3, G4double p1, |
---|
476 | G4double p2, G4double p3, G4double E, G4double r2); |
---|
477 | |
---|
478 | /** |
---|
479 | * ForceAbsor |
---|
480 | */ |
---|
481 | void forceAbsor(G4int *nopart, G4int *iarem, G4int *izrem, G4double *esrem, G4double *erecrem, |
---|
482 | G4double *alrem, G4double *berem, G4double *garem, G4int *jrem); |
---|
483 | |
---|
484 | /** |
---|
485 | * ForceAbs |
---|
486 | * @param iprojo projectile |
---|
487 | * @param at target mass number |
---|
488 | * @param zt target charge number |
---|
489 | * @param ep projectile energy |
---|
490 | * @param bmax a double parameter |
---|
491 | * @param pt a double parameter |
---|
492 | * @return absorption probability |
---|
493 | */ |
---|
494 | G4double forceAbs(G4double iprojo, G4double at, G4double zt, G4double ep, G4double bmax, G4double pt); |
---|
495 | |
---|
496 | /** |
---|
497 | * absoprption xsec revised version rkt-97/5 |
---|
498 | * neutron data from barashenkov |
---|
499 | * this gives absorption xsec for given zp,ap,zt,at,e (mev/nucleon) |
---|
500 | * arguement changed to mev; then e=ep/ap mev/nucleon |
---|
501 | * can be used for neutrons also. |
---|
502 | * this has coulomb as ours |
---|
503 | * @param zp projectile charge number |
---|
504 | * @param zp projectile mass number |
---|
505 | * @param zt a double parameter |
---|
506 | * @param zt target charge number |
---|
507 | * @param at target mass number |
---|
508 | * @param ep projectile energy |
---|
509 | */ |
---|
510 | G4double xabs2(G4double zp, G4double ap, G4double zt, G4double at, G4double ep); |
---|
511 | |
---|
512 | /** |
---|
513 | * Standard random number generator. |
---|
514 | * @param *rndm pointer to the variable reserved for random number |
---|
515 | * @param *seed pointer to the random seed |
---|
516 | */ |
---|
517 | void standardRandom(G4double *rndm, G4long *seed); |
---|
518 | |
---|
519 | /** |
---|
520 | * First derivative of a gaussian potential. |
---|
521 | * @param *rndm pointer to the variable reserved for random number |
---|
522 | */ |
---|
523 | void gaussianRandom(G4double *rndm); |
---|
524 | |
---|
525 | /** |
---|
526 | * Safe exponential function which eliminates the CPU under and overflows. |
---|
527 | * @param x a double parameter |
---|
528 | * @return a double value |
---|
529 | */ |
---|
530 | G4double safeExp(G4double x); |
---|
531 | |
---|
532 | /** |
---|
533 | * Nuclear radius |
---|
534 | * @param A mass number (double parameter) |
---|
535 | */ |
---|
536 | G4double radius(G4double A); |
---|
537 | |
---|
538 | /** Parametrisation de la section efficace de réaction calculée par incl4.1 |
---|
539 | * iprojo=1 proton incident, iprojo=2, neutron incident). |
---|
540 | * entre al et u, entre 10 et 100 mev protons, 20 et 100 mev neutrons. |
---|
541 | * bon ordre de grandeur pour les noyaux légers (c, o ...), trés faux |
---|
542 | * a energie sup a 100 mev. |
---|
543 | * (Comment needs to be translated) |
---|
544 | * @param projectile an integer parameter (1 = proton, 2 = neutron) |
---|
545 | * @param E energy of the projectile (double parameter) |
---|
546 | * @param A target mass number (double parameter) |
---|
547 | * @return cross section (double value) |
---|
548 | */ |
---|
549 | G4double crossSection(G4int projectile, G4double E, G4double A); |
---|
550 | |
---|
551 | /** |
---|
552 | * coulombTransm |
---|
553 | * subroutine coulomb_transm(e,fm1,z1,fm2,z2,proba) |
---|
554 | * calcul du coulombien dans lahet (proba de transmission ou |
---|
555 | * d'entree dans le potentiel nucleaire). |
---|
556 | * @param E energy (a double parameter) |
---|
557 | * @param fm1 a double parameter |
---|
558 | * @param z1 a double parameter |
---|
559 | * @param fm2 a double parameter |
---|
560 | * @param z2 a double parameter |
---|
561 | * @return a double value |
---|
562 | */ |
---|
563 | G4double coulombTransm(G4double E, G4double fm1, G4double z1, G4double fm2, G4double z2); |
---|
564 | |
---|
565 | /** |
---|
566 | * Clmb1 |
---|
567 | * @param eta a double parameter \f$\eta = c_2*z_1*z_2*\sqrt{m/E}\f$ |
---|
568 | * @param rho a double parameter \f$\rho = c_3*(r_1+r_2)*\sqrt{mE}\f$ |
---|
569 | * @return a double value |
---|
570 | */ |
---|
571 | G4double clmb1(G4double rho, G4double eta, G4double *ml); |
---|
572 | |
---|
573 | /** |
---|
574 | * First derivative of a gaussian potential. |
---|
575 | * @param eta a double parameter \f$\eta = c_2*z_1*z_2*\sqrt{m/E}\f$ |
---|
576 | * @param rho a double parameter \f$\rho = c_3*(r_1+r_2)*\sqrt{mE}\f$ |
---|
577 | * @return a double value |
---|
578 | */ |
---|
579 | G4double clmb2(G4double rho, G4double eta, G4double *t1); |
---|
580 | |
---|
581 | public: // Utilities |
---|
582 | /** |
---|
583 | * Returns the smaller of two numbers. |
---|
584 | * @param a a double value |
---|
585 | * @param b a double value |
---|
586 | * @return a double value |
---|
587 | */ |
---|
588 | G4double min(G4double a, G4double b); |
---|
589 | |
---|
590 | /** |
---|
591 | * Returns the smaller of two numbers. |
---|
592 | * @param a an integer value |
---|
593 | * @param b an integer value |
---|
594 | * @return an integer value |
---|
595 | */ |
---|
596 | G4int min(G4int a, G4int b); |
---|
597 | |
---|
598 | /** |
---|
599 | * Returns the greater of two numbers. |
---|
600 | * @param a a double value |
---|
601 | * @param b a double value |
---|
602 | * @return a double value |
---|
603 | */ |
---|
604 | G4double max(G4double a, G4double b); |
---|
605 | |
---|
606 | /** |
---|
607 | * Returns the greater of two numbers. |
---|
608 | * @param a an integer value |
---|
609 | * @param b an integer value |
---|
610 | * @return an integer value |
---|
611 | */ |
---|
612 | G4int max(G4int a, G4int b); |
---|
613 | |
---|
614 | /** |
---|
615 | * Rounds a double to the nearest int |
---|
616 | * @param a double parameter |
---|
617 | * @return an integer value |
---|
618 | */ |
---|
619 | G4int nint(G4double number); |
---|
620 | |
---|
621 | /** |
---|
622 | * Calls a function |
---|
623 | * @param functionChoice an integer value representing the choice of |
---|
624 | * function (0 = wsax, 1 = derivWsax, 2 = dmho, 3 = derivMho, 4 = derivGaus) |
---|
625 | * @param r a double parameter |
---|
626 | * @return a double value |
---|
627 | */ |
---|
628 | G4double callFunction(G4int functionChoice, G4double r); |
---|
629 | |
---|
630 | G4double am(G4double a, G4double b, G4double c, G4double d); |
---|
631 | G4double pcm(G4double e, G4double a, G4double c); |
---|
632 | G4double sign(G4double a, G4double b); |
---|
633 | G4double utilabs(G4double a); |
---|
634 | G4double amax1(G4double a, G4double b); |
---|
635 | G4double w(G4double a, G4double b, G4double c, G4double d); |
---|
636 | G4int idnint(G4double a); |
---|
637 | private: |
---|
638 | |
---|
639 | /** |
---|
640 | * Random seeds for INCL4 internal random number generators. |
---|
641 | */ |
---|
642 | G4Hazard *hazard; |
---|
643 | |
---|
644 | /** |
---|
645 | * Data structure for INCL4. |
---|
646 | */ |
---|
647 | G4Dton *dton; |
---|
648 | |
---|
649 | /** |
---|
650 | * Data structure for INCL4. Contains the Woods-Saxon potential |
---|
651 | * functions for target nuclei. |
---|
652 | */ |
---|
653 | G4Saxw *saxw; |
---|
654 | |
---|
655 | /** |
---|
656 | * Data structure for INCL4. |
---|
657 | */ |
---|
658 | G4Ws *ws; |
---|
659 | |
---|
660 | /** |
---|
661 | * G4Spl2 |
---|
662 | */ |
---|
663 | G4Spl2 *spl2; |
---|
664 | |
---|
665 | /** |
---|
666 | * G4LightGausNuc |
---|
667 | */ |
---|
668 | G4LightGausNuc *light_gaus_nuc; |
---|
669 | |
---|
670 | /** |
---|
671 | * G4LightNuc |
---|
672 | */ |
---|
673 | G4LightNuc *light_nuc; |
---|
674 | |
---|
675 | /** |
---|
676 | * G4Calincl |
---|
677 | */ |
---|
678 | G4Calincl *calincl; |
---|
679 | |
---|
680 | /** |
---|
681 | * G4Mat |
---|
682 | */ |
---|
683 | G4Mat *mat; |
---|
684 | |
---|
685 | /** |
---|
686 | * |
---|
687 | */ |
---|
688 | G4Bl1 *bl1; |
---|
689 | |
---|
690 | /** |
---|
691 | * |
---|
692 | */ |
---|
693 | G4Bl2 *bl2; |
---|
694 | |
---|
695 | /** |
---|
696 | * |
---|
697 | */ |
---|
698 | G4Bl3 *bl3; |
---|
699 | |
---|
700 | /** |
---|
701 | * |
---|
702 | */ |
---|
703 | G4Bl4 *bl4; |
---|
704 | |
---|
705 | /** |
---|
706 | * |
---|
707 | */ |
---|
708 | G4Bl5 *bl5; |
---|
709 | |
---|
710 | /** |
---|
711 | * |
---|
712 | */ |
---|
713 | G4Bl6 *bl6; |
---|
714 | |
---|
715 | /** |
---|
716 | * |
---|
717 | */ |
---|
718 | G4Bl8 *bl8; |
---|
719 | |
---|
720 | /** |
---|
721 | * |
---|
722 | */ |
---|
723 | G4Bl9 *bl9; |
---|
724 | |
---|
725 | /** |
---|
726 | * |
---|
727 | */ |
---|
728 | G4Bl10 *bl10; |
---|
729 | |
---|
730 | /** |
---|
731 | * |
---|
732 | */ |
---|
733 | G4Kind *kindstruct; |
---|
734 | |
---|
735 | /** |
---|
736 | * |
---|
737 | */ |
---|
738 | G4Paul *paul; |
---|
739 | |
---|
740 | /** |
---|
741 | * Detailed information of the cascade |
---|
742 | */ |
---|
743 | G4VarAvat *varavat; |
---|
744 | |
---|
745 | /** |
---|
746 | * Cascade output. |
---|
747 | */ |
---|
748 | G4VarNtp *varntp; |
---|
749 | |
---|
750 | /** |
---|
751 | * For storing the results of the evaporation. |
---|
752 | */ |
---|
753 | G4VarNtp *evaporationResult; |
---|
754 | |
---|
755 | /** |
---|
756 | * Defines the verbosity of console output. Values can be between 0 |
---|
757 | * and 4 where 0 means silent and 4 the most verbose possible |
---|
758 | * output. |
---|
759 | */ |
---|
760 | G4int verboseLevel; |
---|
761 | |
---|
762 | /** |
---|
763 | * Function ID for wsax. |
---|
764 | * @see wsax |
---|
765 | * @see integrate |
---|
766 | * @see callFunction |
---|
767 | */ |
---|
768 | G4int wsaxFunction; |
---|
769 | |
---|
770 | /** |
---|
771 | * Function ID for derivWsax. |
---|
772 | * @see derivWsax |
---|
773 | * @see integrate |
---|
774 | * @see callFunction |
---|
775 | */ |
---|
776 | G4int derivWsaxFunction; |
---|
777 | |
---|
778 | /** |
---|
779 | * Function ID for dmho. |
---|
780 | * @see derivWsax |
---|
781 | * @see integrate |
---|
782 | * @see callFunction |
---|
783 | */ |
---|
784 | G4int dmhoFunction; |
---|
785 | |
---|
786 | /** |
---|
787 | * Function ID for derivMho. |
---|
788 | * @see derivMho |
---|
789 | * @see integrate |
---|
790 | * @see callFunction |
---|
791 | */ |
---|
792 | G4int derivMhoFunction; |
---|
793 | |
---|
794 | /** |
---|
795 | * Function ID for derivGaus. |
---|
796 | * @see derivGaus |
---|
797 | * @see integrate |
---|
798 | * @see callFunction |
---|
799 | */ |
---|
800 | G4int derivGausFunction; |
---|
801 | |
---|
802 | /** |
---|
803 | * Function ID for dens. |
---|
804 | * @see dens |
---|
805 | * @see integrate |
---|
806 | * @see callFunction |
---|
807 | */ |
---|
808 | G4int densFunction; |
---|
809 | |
---|
810 | G4int kind[300]; //= (*kind_p); |
---|
811 | G4double ep[300]; // = (*ep_p); |
---|
812 | G4double alpha[300]; // = (*alpha_p); |
---|
813 | G4double beta[300]; // = (*beta_p); |
---|
814 | G4double gam[300]; // = (*gam_p); |
---|
815 | |
---|
816 | G4Volant *volant; |
---|
817 | G4Abla *abla; |
---|
818 | }; |
---|
819 | |
---|
820 | #endif |
---|