Re: FLUKA: Access to models ?
Hi all
the excitation energy CANNOT be accessed at mgdraw level. TVCMS if
everything is going correctly should come out =0 (no residual excitation
energy). The excitation energy after PEANUT cannot be accessed and it
would not be physically meaningful, since it depends a lot on the
termination conditions and thresholds of the preequilibrium part (where
the code assumes equilibrium is reached). The excitation energy after the
cascade part which perhaps was the one Laurent was referring to, is also
of course strongly dependent on the boundary between cascade and
prequilibrium and does not bear any real physics meaning.
Also Ekres,Px,y,zres are not necessarily updated correctly: their status
depends on the physics input cards set by the user and if
fission/fragmentation
occurred etc. DON'T rely on any of these variables.
Ciao
Alfredo
On Fri, 6 Apr 2001, Laurent APHECETCHE wrote:
> 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
> >
>
>
--
+----------------------------------------------------------------------------+
| Alfredo Ferrari || Tel.: +41.22.767.6119 |
| CERN-SL || Fax.: +41.22.767.7555 |
| 1211 Geneva 23 || e-mail: Alfredo.Ferrari@cern.ch |
| Switzerland || Alfredo.Ferrari@mi.infn.it |
+----------------------------------------------------------------------------+
Partial thread listing:
- Re: FLUKA: Access to models ?, (continued)