RE: usdraw & mgdraw.f

From: Vasilis Vlachoudis (Vasilis.Vlachoudis@cern.ch)
Date: Sat Jul 07 2007 - 09:47:49 CEST

  • Next message: Raul Hennings-Yeomans: "peanut crash?"

    Dear Ketil,

    for the moment the information you want is not existing in the RESNUC
    check the mail in the fluka-list http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1000.html

    I didn't go through your mgdraw routine to check if it is ok, but I saw in your input you are using natural Silicon with a content of 3% of Si-30 and 4.7% of Si-29.

    Vasilis

    -----Original Message-----
    From: owner-fluka-discuss@mi.infn.it [mailto:owner-fluka-discuss@mi.infn.it] On Behalf Of Ketil Rĝed
    Sent: Friday, July 06, 2007 02:06
    To: fluka-discuss@fluka.org
    Subject: usdraw & mgdraw.f

    Dear all FLUKA users!

    I am running a simple example to get familiar with the mgdraw.f routine
    of FLUKA. I am interested in scoring the inelastic
    cross section for protons on a silicon target. In addition I would like
    to score information on the
    secondary particle production.
    My geometry is very simple as can be seen in the attahced .inp file; to
    reduce the probability of multiple
    scattering I am using a thin (0.01 cm) silicon target.

    I am using the entry USDRAW in mgdraw.f to score particles produced
    through inelasting and elastic
    interaction. Through information available in (TRACKR), (GENSTK),
    (RESNUC) and (FHEAVY) I can score different types
    of particle parameters.
    (TRACKR): TRACKs Recording (properties of the currently transported
    particle and its path)
    (GENSTK): properties of each secondary created in a hadronic event
    (RESNUC): properties of the current residual nucleus
    (FHEAVY): stack of heavy secondaries created in nuclear evaporation

    When an interaction takes place I dump information of interest about the
    primary particle,secondaries,
    heavy secondaries and the residual nucleus. See my attached mgdraw.f for
    detailed information.
    Output is formatted and structured for easy post processing, and an
    example is given here:

    Column 1 indicates source of information:

    1: primary particle (TRACKR) (+NP & NPHEAV)
    '1
    ',NCASE,JTRACK,NPFLKA,ICODE,NP,NPHEAV,ETRACK-AM(JTRACK),XSCO,YSCO,ZSCO,WTRACK

    2: secondaries (GENSTK)
    '2
    ',NCASE,JTRACK,ip,ICODE,KPART(ip),ZPART,APART,TKI(ip),CXR(ip),CYR(ip),CZR(ip),WEI(ip)

    3: residual nucleus (RESNUC)
    '3 ',NCASE,JTRACK,'1',ICODE,'0',ICRES,IBRES,EKRES,PXRES,PYRES,PZRES,'1.0'

    4: heavy secondaries (FHEAVY)
    '4',NCASE,JTRACK,ip,ICODE,KHEAVY(ip),ICHEAV(KHEAVY(ip)),IBHEAV(KHEAVY(ip)),TKHEAV(ip),CXHEAV(ip),CYHEAV(ip),CZHEAV(ip),WHEAVY(ip)

    1 37735 1 0 100 1 0 0.099993 0.000000 0.000000 0.000882 1.0
    2 37735 1 1 100 1 0 0 0.09992 0.08741 -0.10571 0.99055 1.0
    3 37735 1 1 100 0 14 28 0.00007 -0.03891 0.04684 0.00437 1.0
     0
    1 39744 1 0 100 1 0 0.099975 -0.000003 -0.000001 0.002867 1.0
    2 39744 1 1 100 1 0 0 0.09997 0.04786 -0.00442 0.99884 1.0
    3 39744 1 1 100 0 14 28 0.00001 -0.02183 0.00147 0.00053 1.0
     0
    1 42027 1 0 101 4 2 0.099872 0.000000 0.000001 0.006844 1.0
    2 42027 1 1 101 1 0 0 0.03241 0.81883 0.06937 0.56982 1.0
    2 42027 1 2 101 1 0 0 0.01079 -0.49355 -0.86955 0.01712 1.0
    2 42027 1 3 101 7 0 0 0.00906 -0.95249 0.29478 0.07666 1.0
    2 42027 1 4 101 7 0 0 0.00439 0.51260 0.49615 -0.70077 1.0
    3 42027 1 1 101 0 0 0 0.00000 0.00000 0.00000 0.00000 1.0
    4 42027 1 1 101 7 6 12 0.01188 -0.12700 0.67350 0.72819 1.0
    4 42027 1 2 101 8 7 15 0.00247 -0.23353 -0.93209 -0.27690 1.0

    Different information is given for the primaries compared to the
    secondaries
    The main information I want is the particle type,
    energy, direction, and type of interaction. For the secondaries and
    heavy secondaries
    GENSTK and FHEAVY contains information on direction, however I can not
    find the same information
    in RESNUC. Thus, for the moment I give the momentum components instead.
    Can anyone tell me where to find the direction of the residual nucleus?
    I would also like if someone can check if my usdraw/mgdraw.f routine and
    my settings in the inp file
    is correct with regards to my task; scoring the inelastic cross section
    and all particles produced
    in the interactions (at the interaction point).
    The reason for the last question is that in a few cases I see something
    I don't understand for the elastic
    interactions. The two blocks of output data below are taken from the 2nd
    of 5 runs with
    a 100 MeV proton beam (100000 primaries per run). See attached file
    "proton100mev002_secondary.log" for source (produced by mgdraw.f).

    For ncase 57831 and 60628 I have problems balancing numbers of particles
    before and after interaction.
    How can I get a Silicon with nuclear mass numbers of 29 & 30 from a
    proton on Si-28 elastic interaction?

    1 57831 1 0 100 1 0 0.099980 -0.000001 0.000000 0.001614 1.0
    2 57831 1 1 100 1 0 0 0.09998 -0.01587 0.02210 0.99963 1.0
    3 57831 1 1 100 0 14 30 0.00000 0.00676 -0.00985 0.00017 1.0

    1 60628 1 0 100 1 0 0.099947 0.000000 0.000001 0.002527 1.0
    2 60628 1 1 100 1 0 0 0.09820 -0.27552 0.58946 0.75936 1.0
    3 60628 1 1 100 0 14 29 0.00174 0.12146 -0.25924 0.11006 1.0

    Am I using the mgdraw.f routine in an incorrect way? Are any of my
    three settings related to physics
    incorrect? Or is my interpretation of the RESNUC ICRES and IBRES values
    wrong?

    I appreciate any corrections or comments.
    Best regards,
    Ketil Rĝed
    PhD student
    University of Bergen


  • Next message: Raul Hennings-Yeomans: "peanut crash?"

    This archive was generated by hypermail 2.1.6 : Sat Jul 07 2007 - 13:07:20 CEST