*$ CREATE GENSTK2.ADD *COPY GENSTK2 * *=== Genstk2 ==========================================================* * *----------------------------------------------------------------------* * * * Partial (extensions to the original one by the Leipzig group) * * Copyright (C) 2019-2022 by Alberto Fasso`, Alfredo Ferrari, * * All Rights Reserved. their collaborators, and INFN * * (Copyright (C) 1996-2019 by INFN and CERN, for the * * extensions only) * * * * * * 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 !!!! * * !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! * * * * Authors: Alfredo Ferrari, Johannes Ranft, * * & Paola Sala * * * * * * Created on 20 January 1996 by Alfredo Ferrari & Paola Sala * * Infn - Milan * * * * Last change on 16-Mar-22 by Alfredo Ferrari * * Private * * * * /Genstk/ is the storage for secondaries created in hadronic * * events * * Np = total number of secondaries * * Itrn (ip) = (Paprop) id of the ip_th secondary * * If |Itrn| > 100000 it means a heavy ion with * * |Itrn| = Kh + 100 A + 100000 Z + 10000000 m * * + 1000000000 h, where Kh is the correspon-* * ding Kheavy index, m is a (possible) * * isomeric state (m=99 means excited state* * defined by ilvsec), h is the number of * * hyperons for a (possible) hypernucleus * * Ilvsec (ip) = possible excited state index of the ip_th secon- * * dary + 10000, if < 99 is the isomeric state index * * * Infext (k,ip) = possible extra infos for the ip_th secondary * * (k=1,-1000000000 means primary temporary on the * * stack for bookkeeping purposes only, * * k=1,-2000000000 means secondary temporary on the * * stack for bookkeeping purposes only ) * * Cxrn (ip) = x-axis direction cosine of the ip_th secondary * * Cyrn (ip) = y-axis direction cosine of the ip_th secondary * * Czrn (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) * * Elr (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 of the ip_th secondary * * Eexsec (ip) = if positive, it is the excitation energy of a * * possible excited state (0 for isomers!) of the * * ip_th secondary (GeV) * * Tmnsec (ip) = If positive, it is the mean life (s) of a possible* * excited state/isomer of the ip_th secondary * * 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) * * * *----------------------------------------------------------------------* * * Integer/double precision declarations added in order to make easier * the compatibility with the latest Dpmjet-3 versions: INTEGER MXP , ITRN , INFEXT, ILVSEC, NP0 , IRN DOUBLE PRECISION CXRN , CYRN , CZRN , CXRPOL, CYRPOL, CZRPOL, & ELR , PLR , WEI , AGESEC, EEXSEC, TMNSEC, & TV , TVCMS , TVRECL, TVHEAV, TVBIND PARAMETER ( MXP = MXPSCS ) * COMMON / GENSTK / CXRN (MXPSCS), CYRN (MXPSCS), & CZRN (MXPSCS), CXRPOL (MXPSCS), CYRPOL (MXPSCS), & CZRPOL (MXPSCS), ELR (MXPSCS), PLR (MXPSCS), & WEI (MXPSCS), AGESEC (MXPSCS), EEXSEC (MXPSCS), & TMNSEC (MXPSCS), TV , TVCMS , TVRECL, TVHEAV, & TVBIND, & ITRN (MXPSCS), INFEXT (4,MXPSCS), & ILVSEC (MXPSCS), NP0 , IRN SAVE / GENSTK /