]>
Commit | Line | Data |
---|---|---|
1 | #ifndef FFLKSTK_H | |
2 | #define FFLKSTK_H 1 | |
3 | ||
4 | #include "cfortran.h" | |
5 | #include "Rtypes.h" | |
6 | ||
7 | #include "Fdimpar.h" | |
8 | ||
9 | extern "C" { | |
10 | //*$ create flkstk.add | |
11 | //*copy flkstk | |
12 | //* | |
13 | //*=== flkstk ============================================================* | |
14 | //* | |
15 | //*----------------------------------------------------------------------* | |
16 | //* * | |
17 | //* FLUKA90-200x particle stack: * | |
18 | //* * | |
19 | //* Changes: last change on 15-may-2005 by Alfredo Ferrari * | |
20 | //* INFN, Milan * | |
21 | //* * | |
22 | //* * | |
23 | //* description of the common block(s) and variable(s) * | |
24 | //* * | |
25 | //* /Flkstk/ stack for the primaries * | |
26 | //* Wtflk = particle statistical weight * | |
27 | //* Pmoflk = particle (laboratory) momentum (GeV/c) * | |
28 | //* Tkeflk = particle (laboratory) kinetic energy (GeV) * | |
29 | //* Xflk = particle position x-coordinate * | |
30 | //* Yflk = particle position y-coordinate * | |
31 | //* Zflk = particle position z-coordinate * | |
32 | //* Txflk = particle direction x-coordinate * | |
33 | //* Tyflk = particle direction y-coordinate * | |
34 | //* Tzflk = particle direction z-coordinate * | |
35 | //* Txpol = x direction cosine of the particle polarization * | |
36 | //* Typol = y direction cosine of the particle polarization * | |
37 | //* Tzpol = z direction cosine of the particle polarization * | |
38 | //* Txnor = x direction cosine of a (possible) surface normal * | |
39 | //* Tynor = y direction cosine of a (possible) surface normal * | |
40 | //* Tznor = z direction cosine of a (possible) surface normal * | |
41 | //* Dfnear = distance to the nearest boundary * | |
42 | //* Agestk = age of the particle (seconds) * | |
43 | //* Aknshr = Kshort component of K0/K0bar * | |
44 | //* Frcphn = cross section for force photonuclear interaction (if * | |
45 | //* < 0), distance to a forced photonuclear interaction * | |
46 | //* (if > 0) * | |
47 | //* Lfrphn = flag for forced photonuclear interaction * | |
48 | //* Raddly = delay (s) in production wrt the nominal primary "0" * | |
49 | //* time for particle produced in radioactive decays * | |
50 | //* (i.e. those coming from decays of daughter isotopes), * | |
51 | //* when in analogue mode, flag for position in the * | |
52 | //* activr array when in non-analogue mode * | |
53 | //* Cmpath = cumulative path travelled by the particle since it * | |
54 | //* was produced (cm) * | |
55 | //* Sparek = spare real variables available for K.W.Burn * | |
56 | //* Ispark = spare integer variables available for K.W.Burn * | |
57 | //* Iloflk = particle identity (Paprop numbering) * | |
58 | //* Igroup = energy group for low energy neutrons * | |
59 | //* Loflk = particle generation * | |
60 | //* Louse = user flag * | |
61 | //* Nrgflk = particle region number * | |
62 | //* Nlattc = particle lattice cell number * | |
63 | //* Nhspnt = pointer to the history object (Geant4 geometry) * | |
64 | //* Nevent = number of the event which created the particle * | |
65 | //* Numpar = particle number * | |
66 | //* Lraddc = flag for particles generated in radioactive decays * | |
67 | //* Nparma = largest particle number ever reached * | |
68 | //* Nstmax = highest value of the stack pointer ever reached * | |
69 | //* in the run * | |
70 | //* Npflka = Fluka stack pointer * | |
71 | //* Nstaol = stack pointer of the last processed particle * | |
72 | //* Igroun = energy group number of the last processed particle * | |
73 | //* if it is a low energy neutron * | |
74 | //*----------------------------------------------------------------------* | |
75 | //* | |
76 | ||
77 | typedef struct { | |
78 | Double_t xflk[mfstck+1]; //(0:MFSTCK) | |
79 | Double_t yflk[mfstck+1]; //(0:MFSTCK) | |
80 | Double_t zflk[mfstck+1]; //(0:MFSTCK) | |
81 | Double_t txflk[mfstck+1]; //(0:MFSTCK) | |
82 | Double_t tyflk[mfstck+1]; //(0:MFSTCK) | |
83 | Double_t tzflk[mfstck+1]; //(0:MFSTCK) | |
84 | Double_t txpol[mfstck+1]; //(0:MFSTCK) | |
85 | Double_t typol[mfstck+1]; //(0:MFSTCK) | |
86 | Double_t tzpol[mfstck+1]; //(0:MFSTCK) | |
87 | Double_t txnor[mfstck+1]; //(0:MFSTCK) | |
88 | Double_t tynor[mfstck+1]; //(0:MFSTCK) | |
89 | Double_t tznor[mfstck+1]; //(0:MFSTCK) | |
90 | Double_t wtflk[mfstck+1]; //(0:MFSTCK) | |
91 | Double_t pmoflk[mfstck+1]; //(0:MFSTCK) | |
92 | Double_t tkeflk[mfstck+1]; //(0:MFSTCK) | |
93 | Double_t dfnear[mfstck+1]; //(0:MFSTCK) | |
94 | Double_t agestk[mfstck+1]; //(0:MFSTCK) | |
95 | Double_t aknshr[mfstck+1]; //(0:MFSTCK) | |
96 | Double_t raddly[mfstck+1]; //(0:MFSTCK) | |
97 | Double_t cmpath[mfstck+1]; //(0:MFSTCK) | |
98 | Double_t frcphn[mfstck+1]; //(0:MFSTCK) | |
99 | Double_t sparek[mfstck+1][mkbmx1]; //(MKBMX1,0:MFSTCK) | |
100 | Int_t ispark[mfstck+1][mkbmx2]; //(MKBMX2,0:MFSTCK) | |
101 | Int_t iloflk[mfstck+1]; //(0:MFSTCK) | |
102 | Int_t igroup[mfstck+1]; //(0:MFSTCK) | |
103 | Int_t loflk[mfstck+1]; //(0:MFSTCK) | |
104 | Int_t louse[mfstck+1]; //(0:MFSTCK) | |
105 | Int_t nrgflk[mfstck+1]; //(0:MFSTCK) | |
106 | Int_t nlattc[mfstck+1]; //(0:MFSTCK) | |
107 | Int_t nhspnt[mfstck+1]; //(0:MFSTCK) | |
108 | Int_t nevent[mfstck+1]; //(0:MFSTCK) | |
109 | Int_t numpar[mfstck+1]; //(0:MFSTCK) | |
110 | Int_t lraddc[mfstck+1]; //(0:MFSTCK) | |
111 | Int_t lfrphn[mfstck+1]; //(0:MFSTCK) | |
112 | Int_t nparma; | |
113 | Int_t nstmax; | |
114 | Int_t npflka; | |
115 | Int_t nstaol; | |
116 | Int_t igroun; | |
117 | } flkstkCommon; | |
118 | #define FLKSTK COMMON_BLOCK(FLKSTK,flkstk) | |
119 | COMMON_BLOCK_DEF(flkstkCommon,FLKSTK); | |
120 | } | |
121 | ||
122 | #endif |