Re: FLUKA: Access to models ?
Hi Stefan,
thanks a lot for those informations.
by the way, where did you get the common comments ? In the version of
FLUKA I have (nov.28), there is no such information...
Regards,
Stefan Roesler wrote:
>
> Laurent:
>
> You have access to the information on a certain interaction using the
> MGDRAW user-routine (ENTRY USDRAW) which is activated by the DUMPTHEM
> card. The variable ICODE allows to trigger on certain classes of
> interactions (see comments in MGDRAW). For example, after an inelastic
> collision (ICODE=101) a lot of useful information is stored in COMMONs
> FINUC, FHEAVY and RESNUC (see decriptions below). For example, the
> 4-momentum of the recoil nucleus is stored in EKRES,PXRES,PYRES,PZRES, the
> excitation energy in (I believe) TVCMS, etc.
>
> However, as this is a quite non-standard application, use the information
> in the above COMMONs with care and always cross check if the variables
> really contain what you think they contain (e.g. energy-momentum
> conservation checks of in- and outgoing particles/fragments,
> 4-momentum-mass consistency checks, etc.)
>
> Stefan
>
> ___________________________________
>
> Stefan Roesler
> SLAC
> P.O. Box 4349, MS 48
> Stanford, CA 94309
>
> Tel: 650-926-2048
> Fax: 650-926-3569
> E-mail: sroesler@slac.stanford.edu
> Stefan.Roesler@cern.ch
>
> *----------------------------------------------------------------------*
> * *
> * 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) *
> * *
> *----------------------------------------------------------------------*
> *
> PARAMETER (MXP=999)
> *
> COMMON / FINUC / CXR (MXP), CYR (MXP), CZR (MXP),
> & CXRPOL (MXP), CYRPOL (MXP), CZRPOL (MXP),
> & TKI (MXP), PLR (MXP), WEI (MXP),
> & AGESEC (MXP), TV, TVCMS, TVRECL, TVHEAV, TVBIND,
> & NP0, NP, KPART (MXP)
>
> *----------------------------------------------------------------------*
> * *
> * include file: fheavy *
> * *
> * Created on 5 april 1990 by Alfredo Ferrari, INFN Milan *
> * *
> * Last change on 26-jul-97 by Alfredo Ferrari, INFN Milan *
> * *
> * included in the following subroutines or functions: not updated *
> * *
> * description of the common block(s) and variable(s) *
> * *
> * /fheavy/ is the storage for heavy secondaries created in the *
> * nuclear evaporation *
> * npheav = number of secondaries *
> * kheavy(ip) = type of the secondary ip *
> * ( 3 = deuteron, 4 = 3-H, 5 = 3-He, 6 = 4-He, *
> * 7-12 = "Heavy" fragment specified by Ibheav and *
> * Icheav ) *
> * cxheav(ip) = direction cosine of the secondary ip *
> * with respect to x-axis *
> * cyheav(ip) = direction cosine of the secondary ip *
> * with respect to y-axis *
> * czheav(ip) = direction cosine of the secondary ip *
> * with respect to z-axis *
> * tkheav(ip) = kinetic energy of secondary ip *
> * pheavy(ip) = momentum of the secondary ip *
> * wheavy(ip) = weight of the secondary ip *
> * agheav(ip) = "age" of the secondary ip with respect to the *
> * interaction time *
> * amheav(kp) = atomic masses of the twelve types of evaporated *
> * or fragmented or fissioned particles *
> * amnhea(kp) = nuclear masses of the twelve types of evaporated *
> * or fragmented or fissioned particles *
> * bhheav(jp,kp) = (nuclear) binding energy of the jp_th hyperon of *
> * the kp-type heavy particle *
> * anheav(kp) = name of the kp-type heavy particle *
> * icheav(kp) = charge of the kp-type heavy particle *
> * ibheav(kp) = mass number of the kp-type heavy particle *
> * imheav(kp) = isomeric state of the kp-type heavy particle *
> * ihheav(kp) = number of hyperons of the kp-type heavy particle *
> * khheav(jp,kp) = id of the jp_th hyperon of the kp-type heavy *
> * particle *
> * !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
> * !!! There is now the possibility to produce up to 6 "heavy" !!!! *
> * !!! fragments besides the residual nucleus recorded in !!!! *
> * !!! Resnuc: they are identified by indeces 7-12, of course !!!! *
> * !!! the corresponding physical properties (Z,A,m..) must be !!!! *
> * !!! updated every time they are produced !!!! *
> * !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
> *----------------------------------------------------------------------*
> *
> PARAMETER ( MXHEAV = 100 )
> PARAMETER ( KXHEAV = 12 )
> CHARACTER*8 ANHEAV
> COMMON / FHEAVY / CXHEAV (MXHEAV), CYHEAV (MXHEAV),
> & CZHEAV (MXHEAV), TKHEAV (MXHEAV),
> & PHEAVY (MXHEAV), WHEAVY (MXHEAV),
> & AGHEAV (MXHEAV), BHHEAV (IHYPMX,KXHEAV),
> & AMHEAV (KXHEAV), AMNHEA (KXHEAV),
> & KHEAVY (MXHEAV), ICHEAV (KXHEAV),
> & IBHEAV (KXHEAV), IMHEAV (KXHEAV),
> & IHHEAV (KXHEAV), KHHEAV (IHYPMX,KXHEAV), NPHEAV
> COMMON / FHEAVC / ANHEAV (KXHEAV)
>
> *---------------------------------------------------------------------*
> * *
> * Include file Resnuc *
> * *
> * Created on 20 april 1990 by Alfredo Ferrari *
> * INFN Milan *
> * *
> * Last change on 09-nov-99 by Alfredo Ferrari, INFN - Milan *
> * *
> * *
> * Description of variables (incomplete): *
> * *
> * Icres = residual nucleus atomic number *
> * Ibres = residual nucleus mass number *
> * Istres = residual nucleus stable level index *
> * Ismres = residual nucleus isomeric state index *
> * Ihyres = residual nucleus hyperon number *
> * Amnres = residual nucleus nuclear mass *
> * Ammres = residual nucleus atomic mass *
> * Eres = residual nucleus total energy *
> * Ekres = residual nucleus kinetic energy *
> * Px,y,zres = residual nucleus momentum components *
> * Ptres2 = residual nucleus squared momentum *
> * Angres = residual nucleus angular momentum (GeV/c fm) *
> * Anx,y,zres = residual nucleus angular momentum components *
> * Khyres(jp) = id of the jp_th hyperon in the residual nucleus *
> * Bhyres(jp) = (nuclear) binding energy of the jp_th hyperon *
> * in the residual nucleus *
> * *
> * *
> * *
> *---------------------------------------------------------------------*
> *
> LOGICAL LRNFSS, LFRAGM
> COMMON /RESNUC/ AMNTAR, AMMTAR, AMNZM1, AMMZM1, AMNNM1, AMMNM1,
> & ANOW, ZNOW, ANCOLL, ZNCOLL, AMMLFT, AMNLFT,
> & ERES, EKRES, AMNRES, AMMRES, PTRES, PXRES,
> & PYRES, PZRES, PTRES2, ANGRES, ANXRES, ANYRES,
> & ANZRES, BHYRES (IHYPMX),
> & KTARP, KTARN, IGREYP, IGREYN, IPDPM0, IPDPM1,
> & IPREEH, IPRDEU, IPRTRI, IPR3HE, IPR4HE, ICRES,
> & IBRES, ISTRES, ISMRES, IHYRES, IEVAPL, IEVAPH,
> & IEVPHO, IEVNEU, IEVPRO, IEVDEU, IEVTRI, IEV3HE,
> & IEV4HE, IDEEXG, IBTAR, ICHTAR, IBLEFT, ICLEFT,
> & IOTHER, KHYRES (IHYPMX), LRNFSS, LFRAGM
>
--
Dr. Laurent APHECETCHE (mailto:aphecetc@in2p3.fr) (IN2P3-CNRS)
SUBATECH-EMN-4 rue Alfred Kastler-BP 20722-44307 NANTES cedex 03
TEL (+33/0) 2 51 85 84 17 - FAX (+33/0) 2 51 85 84 24 (France)
Collaborations PHENIX http://www.phenix.bnl.gov/~aphecetc et MEGAPIE.
Partial thread listing:
- Re: FLUKA: Access to models ?, (continued)