Update to FLUKA 2006.3b.

diff --git a/TFluka/Fdimpar.h b/TFluka/Fdimpar.h
index 44fdc2b0e9da35b816efbc50f2b7e3f047a9068b..6843be1c5b4b7b1871620eec34f1ed9f6a68f0b8 100644 (file)
--- a/TFluka/Fdimpar.h
+++ b/TFluka/Fdimpar.h
@@ -90,8 +90,8 @@ extern "C" {
 //*     const Int_t mkbmx2 =    3;
     const Int_t mkbmx1 = 11;
     const Int_t mkbmx2 = 11;
-    const Int_t mxirrd = 100;
-    const Int_t mxtrdc = 120;
+    const Int_t mxirrd = 500;
+    const Int_t mxtrdc = 500;
     const Int_t nktl   = 17;
 }
 

diff --git a/TFluka/Ftrackr.h b/TFluka/Ftrackr.h
index 8528cab2f58f249e59ac0f452bef3479eef26c3f..da33d593f2d8717d04d140513f40720f97416faa 100644 (file)
--- a/TFluka/Ftrackr.h
+++ b/TFluka/Ftrackr.h
@@ -41,7 +41,19 @@ extern "C" {
 //*          dtrack = energy deposition of the jth deposition event      *
 //*          dptrck = momentum loss of the jth deposition event          *
 //*                                                                      *
-//*          jtrack = identity number of the particle                    *
+//*          Jtrack = identity number of the particle: for recoils or    *
+//*                   kerma deposition it can be outside the allowed     *
+//*                   particle id range, assuming values like:           *
+//*                     208: "heavy" recoil                              *
+//*                     211: EM below threshold                          *
+//*                     308: low energy neutron kerma                    *
+//*                   in those cases the id of the particle originating  *
+//*                   the interaction is saved inside J0trck (which othe-*
+//*                   rwise is zero)                                     *
+//*          J0trck = see above                                          *
+//*          Ifltrk = flag used for internal debugging (trying to solve  *
+//*                   possible residual issues with Mgdraw driven        *
+//*                   quenching)                                         *
 //*          etrack = total energy of the particle                       *
 //*          ptrack = momentum of the particle (not always defined, if   *
 //*                 < 0 must be obtained from etrack)                    *
@@ -99,46 +111,48 @@ extern "C" {
 const Int_t mxtrck = 2500;
 
 typedef struct {
-   Double_t xtrack[mxtrck+1];
-   Double_t ytrack[mxtrck+1];
-   Double_t ztrack[mxtrck+1];
-   Double_t ttrack[mxtrck];
-   Double_t dtrack[mxtrck];
-   Double_t dptrck[mxtrck][3];
-   Double_t etrack;
-   Double_t ptrack;
-   Double_t cxtrck;
-   Double_t cytrck;
-   Double_t cztrck;
-   Double_t wtrack;
-   Double_t cxtrpl;
-   Double_t cytrpl;
-   Double_t cztrpl;
-   Double_t zfftrk;
-   Double_t zfrttk;
-   Double_t atrack;
-   Double_t ctrack;
-   Double_t cmtrck;
-   Double_t akshrt;
-   Double_t aklong;
-   Double_t wscrng;
-   Double_t wninou;
-   Double_t spausr[mkbmx1];
-   Double_t sttrck;
-   Double_t satrck;
-   Int_t    ntrack;
-   Int_t    mtrack;
-   Int_t    jtrack;
-   Int_t    ktrack;
-   Int_t    mmtrck;
-   Int_t    lt1trk;
-   Int_t    lt2trk;
-   Int_t    ihspnt;
-   Int_t    ltrack;
-   Int_t    llouse;
-   Int_t    ispusr[mkbmx2];
-   Int_t    lfsssc;
-   Int_t    lpkill;
+    Double_t xtrack[mxtrck+1];
+    Double_t ytrack[mxtrck+1];
+    Double_t ztrack[mxtrck+1];
+    Double_t ttrack[mxtrck];
+    Double_t dtrack[mxtrck];
+    Double_t dptrck[mxtrck][3];
+    Double_t etrack;
+    Double_t ptrack;
+    Double_t cxtrck;
+    Double_t cytrck;
+    Double_t cztrck;
+    Double_t wtrack;
+    Double_t cxtrpl;
+    Double_t cytrpl;
+    Double_t cztrpl;
+    Double_t zfftrk;
+    Double_t zfrttk;
+    Double_t atrack;
+    Double_t ctrack;
+    Double_t cmtrck;
+    Double_t akshrt;
+    Double_t aklong;
+    Double_t wscrng;
+    Double_t wninou;
+    Double_t spausr[mkbmx1];
+    Double_t sttrck;
+    Double_t satrck;
+    Int_t    ntrack;
+    Int_t    mtrack;
+    Int_t    ifltrk;
+    Int_t    jtrack;
+    Int_t    j0trck;
+    Int_t    ktrack;
+    Int_t    mmtrck;
+    Int_t    lt1trk;
+    Int_t    lt2trk;
+    Int_t    ihspnt;
+    Int_t    ltrack;
+    Int_t    llouse;
+    Int_t    ispusr[mkbmx2];
+    Int_t    lfsssc;
+    Int_t    lpkill;
 } trackrCommon;
 #define TRACKR COMMON_BLOCK(TRACKR,trackr)
 COMMON_BLOCK_DEF(trackrCommon,TRACKR);