finuc replaced by genstk
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 1 Sep 2005 13:43:26 +0000 (13:43 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 1 Sep 2005 13:43:26 +0000 (13:43 +0000)
TFluka/Ffinuc.h [deleted file]
TFluka/Fgenstk.h [new file with mode: 0644]

diff --git a/TFluka/Ffinuc.h b/TFluka/Ffinuc.h
deleted file mode 100644 (file)
index 80fe77b..0000000
+++ /dev/null
@@ -1,88 +0,0 @@
-#ifndef FFINUC
-#define FFINUC_H 1
-                                                                                
-#include "Rtypes.h"
-#include "cfortran.h"
-                                                                                
-#include "Fdimpar.h"
-extern "C" {
-//*$ create finuc.add
-//*copy finuc
-//*
-//*=== finuc ============================================================*
-//*
-//*----------------------------------------------------------------------*
-//*                                                                      *
-//*     include file: finuc (new version of old finuc of fluka86)        *
-//*                                                                      *
-//*     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      *
-//*     !!!!     s e e   a l s o   i n c l u d e   f i l e     !!!!      *
-//*     !!!!                 f i n u c 2                       !!!!      *
-//*     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      *
-//*                                                                      *
-//*     created on  20 january 1996  by    alfredo ferrari & paola sala  *
-//*                                                   infn - milan       *
-//*                                                                      *
-//*     last change on 26-jul-97     by    alfredo ferrari               *
-//*                                                                      *
-//*     included in the following subroutines or functions: not updated  *
-//*                                                                      *
-//*     description of the common block(s) and variable(s)               *
-//*                                                                      *
-//*     /finuc/ is the storage for secondaries created in event          *
-//*        np        = number of secondaries                             *
-//*       kpart (ip) = type of the secondary ip                          *
-//*         cxr (ip) = direction cosine of the secondary ip              *
-//*                    with respect to x-axis                            *
-//*         cyr (ip) = direction cosine of the secondary ip              *
-//*                    with respect to y-axis                            *
-//*         czr (ip) = direction cosine of the secondary ip              *
-//*                    with respect to z-axis                            *
-//*      cxrpol (ip) = direction cosine of the secondary ip polarization *
-//*                    with respect to x-axis                            *
-//*      cyrpol (ip) = direction cosine of the secondary ip polarization *
-//*                    with respect to y-axis                            *
-//*      czrpol (ip) = direction cosine of the secondary ip polarization *
-//*                    with respect to z-axis                            *
-//*         tki (ip) = kinetic energy of secondary ip                    *
-//*         plr (ip) = momentum of the secondary ip                      *
-//*         wei (ip) = weight of the secondary ip                        *
-//*      agesec (ip) = "age" of the secondary ip with respect to the     *
-//*                    interaction time                                  *
-//*        tv        = excitation energy                                 *
-//*        tvcms     = actual excitation energy of the residual nucleus  *
-//*        tvrecl    = recoil kinetic energy of the residual nucleus     *
-//*        tvheav    = recoil kinetic energies of heavy (2-h, 3-h, 3-he, *
-//*                    4-he) fragments after evaporation                 *
-//*        tvbind    = approximate energy wasted in nuclear binding      *
-//*                    effects (not yet operational)                     *
-//*                                                                      *
-//*----------------------------------------------------------------------*
-//*
-const Int_t mxp = mxpscs;
-//*
-
-typedef struct {
-   Double_t cxr[mxp];
-   Double_t cyr[mxp];
-   Double_t czr[mxp];
-   Double_t cxrpol[mxp];
-   Double_t cyrpol[mxp];
-   Double_t czrpol[mxp];
-   Double_t tki[mxp];
-   Double_t plr[mxp];
-   Double_t wei[mxp];
-   Double_t agesec[mxp];
-   Double_t tv;
-   Double_t tvcms;
-   Double_t tvrecl;
-   Double_t tvheav;
-   Double_t tvbind;
-   Int_t    np0;
-   Int_t    np;
-   Int_t    kpart[mxp];
-} finucCommon;
-#define FINUC COMMON_BLOCK(FINUC,finuc)
-COMMON_BLOCK_DEF(finucCommon,FINUC);
-}
-#endif
diff --git a/TFluka/Fgenstk.h b/TFluka/Fgenstk.h
new file mode 100644 (file)
index 0000000..888231e
--- /dev/null
@@ -0,0 +1,84 @@
+#ifndef FGENSTK
+#define FGENSTK_H 1
+                                                                                
+#include "Rtypes.h"
+#include "cfortran.h"
+                                                                                
+#include "Fdimpar.h"
+//*$ CREATE GENSTK.ADD
+//*COPY GENSTK
+//*
+//*=== Genstk ===========================================================*
+//*
+//*----------------------------------------------------------------------*
+//*                                                                      *
+//*     hadron GENerator STacK for FLUKA: (new version of old Finuc of   *
+//*     FLUKA86 by J.Ranft)                                              *
+//*                                                                      *
+//*     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      *
+//*     !!!!     S E E   A L S O   I N C L U D E   F I L E     !!!!      *
+//*     !!!!                 G E N S T K 2                     !!!!      *
+//*     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      *
+//*                                                                      *
+//*     Created on  20 january 1996  by    Alfredo Ferrari & Paola Sala  *
+//*                                                   Infn - Milan       *
+//*                                                                      *
+//*     Last change on 15-may-05     by    Alfredo Ferrari               *
+//*                                                                      *
+//*                                                                      *
+//*     /Genstk/ is the storage for secondaries created in hadronic      *
+//*              events                                                  *
+//*        Np        = total number of secondaries                       *
+//*       Kpart (ip) = (Paprop) id of the ip_th secondary                *
+//*         Cxr (ip) = x-axis direction cosine of the ip_th secondary    *
+//*         Cyr (ip) = y-axis direction cosine of the ip_th secondary    *
+//*         Czr (ip) = z-axis direction cosine of the ip_th secondary    *
+//*      Cxrpol (ip) = x-axis direction cosine of the ip_th secondary    *
+//*                    polarization vector (rest frame when applicable)  *
+//*      Cyrpol (ip) = y-axis direction cosine of the ip_th secondary    *
+//*                    polarization vector (rest frame when applicable)  *
+//*      Czrpol (ip) = z-axis direction cosine of the ip_th secondary    *
+//*                    polarization vector (rest frame when applicable)  *
+//*         Tki (ip) = laboratory kinetic energy of ip_th secondary (GeV)*
+//*         Plr (ip) = laboratory momentum of the ip_th secondary (GeV/c)*
+//*         Wei (ip) = statistical weight of the ip_th secondary         *
+//*      Agesec (ip) = "age" of the ip_th secondary with respect to the  *
+//*                    interaction time                                  *
+//*        Tv        = excitation energy (GeV)                           *
+//*        Tvcms     = actual excitation energy of the residual nucleus  *
+//*        Tvrecl    = recoil kinetic energy of the residual nucleus     *
+//*        Tvheav    = recoil kinetic energies of heavy (2-H, 3-H, 3-He, *
+//*                    4-He) fragments after evaporation                 *
+//*        Tvbind    = approximate energy wasted in nuclear binding      *
+//*                    effects (not yet operational)                     *
+//*                                                                      *
+//*----------------------------------------------------------------------*
+
+extern "C" {
+const Int_t mxp = mxpscs;
+//*
+
+typedef struct {
+   Double_t cxr[mxp];
+   Double_t cyr[mxp];
+   Double_t czr[mxp];
+   Double_t cxrpol[mxp];
+   Double_t cyrpol[mxp];
+   Double_t czrpol[mxp];
+   Double_t tki[mxp];
+   Double_t plr[mxp];
+   Double_t wei[mxp];
+   Double_t agesec[mxp];
+   Double_t tv;
+   Double_t tvcms;
+   Double_t tvrecl;
+   Double_t tvheav;
+   Double_t tvbind;
+    Int_t    kpart[mxp];
+   Int_t    np0;
+   Int_t    np;
+} genstkCommon;
+#define GENSTK COMMON_BLOCK(GENSTK,genstk)
+COMMON_BLOCK_DEF(genstkCommon,GENSTK);
+}
+#endif