Re: FLUKA: Access to models ?
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
On Thu, 5 Apr 2001, Laurent APHECETCHE wrote:
> Hi,
>
> Is there a way to access "directly" information from the models used in
> FLUKA during the course of a simulation, e.g. to get distribution of
> excitation energies in nuclei (after PEANUT?) or distribution of recoils
> ?
>
> Regards,
>
> --
> 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)