From: Ketil Rĝed (ketil.roed@student.uib.no)
Date: Fri Jul 06 2007 - 02:05:39 CEST
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
1 05308 1 0 100 1 0 0.099952 -0.000012 -0.000001 0.004324 1.0
2 05308 1 1 100 1 0 0 0.09978 -0.20106 0.06857 0.97718 1.0
3 05308 1 1 100 0 14 28 0.00017 0.08739 -0.03155 0.01052 1.0
0
1 06861 1 0 100 1 0 0.099938 -0.000005 0.000002 0.006099 1.0
2 06861 1 1 100 1 0 0 0.09677 -0.45276 -0.68072 0.57588 1.0
3 06861 1 1 100 0 14 28 0.00317 0.19744 0.29802 0.19278 1.0
0
1 06984 1 0 100 1 0 0.099973 0.000000 0.000001 0.002450 1.0
2 06984 1 1 100 1 0 0 0.09985 0.08948 0.15282 0.98420 1.0
3 06984 1 1 100 0 14 28 0.00012 -0.03970 -0.06747 0.00730 1.0
0
1 07965 1 0 100 1 0 0.099961 -0.000002 -0.000003 0.003618 1.0
2 07965 1 1 100 1 0 0 0.09982 0.16721 0.08945 0.98185 1.0
3 07965 1 1 100 0 14 28 0.00014 -0.07429 -0.04046 0.00838 1.0
0
1 11208 1 0 100 1 0 0.099810 -0.000002 0.000002 0.005360 1.0
2 11208 1 1 100 1 0 0 0.09955 -0.26159 -0.01203 0.96511 1.0
3 11208 1 1 100 0 14 28 0.00026 0.11562 0.00540 0.01609 1.0
0
1 11554 1 0 100 1 0 0.099967 0.000001 0.000000 0.002003 1.0
2 11554 1 1 100 1 0 0 0.09986 0.17014 -0.02108 0.98519 1.0
3 11554 1 1 100 0 14 29 0.00011 -0.07537 0.00929 0.00683 1.0
0
1 11745 1 0 100 1 0 0.099928 -0.000003 0.000001 0.004439 1.0
2 11745 1 1 100 1 0 0 0.09901 0.24275 0.41020 0.87909 1.0
3 11745 1 1 100 0 14 28 0.00092 -0.10793 -0.18159 0.05562 1.0
0
1 12775 1 0 101 6 1 0.099897 -0.000010 -0.000001 0.008404 1.0
2 12775 1 1 101 -3 0 0 0.04614 0.03127 0.07749 0.99650 1.0
2 12775 1 2 101 1 0 0 0.00537 -0.02550 -0.55178 0.83360 1.0
2 12775 1 3 101 7 0 0 0.00269 -0.60047 -0.06463 0.79703 1.0
2 12775 1 4 101 7 0 0 0.00240 0.76231 0.13393 0.63320 1.0
2 12775 1 5 101 7 0 0 0.00183 0.21218 0.87002 -0.44501 1.0
2 12775 1 6 101 7 0 0 0.00410 0.66122 0.32756 0.67491 1.0
3 12775 1 1 101 0 12 24 0.00047 -0.08696 -0.04796 -0.10468 1.0
4 12775 1 1 101 3 1 2 0.00309 0.67231 0.62590 0.39528 1.0
0
1 17201 1 0 100 1 0 0.099968 -0.000003 -0.000001 0.003796 1.0
2 17201 1 1 100 1 0 0 0.09984 -0.10894 0.14444 0.98350 1.0
3 17201 1 1 100 0 14 28 0.00012 0.04779 -0.06411 0.00762 1.0
0
1 18750 1 0 100 1 0 0.099948 -0.000003 -0.000004 0.004149 1.0
2 18750 1 1 100 1 0 0 0.09994 0.02080 0.02932 0.99935 1.0
3 18750 1 1 100 0 14 28 0.00001 -0.01011 -0.01354 0.00030 1.0
0
1 19591 1 0 100 1 0 0.099975 0.000002 0.000001 0.002725 1.0
2 19591 1 1 100 1 0 0 0.09865 0.25945 0.50377 0.82396 1.0
3 19591 1 1 100 0 14 28 0.00132 -0.11394 -0.22251 0.08081 1.0
0
1 21685 1 0 101 8 0 0.099999 0.000000 0.000000 0.000318 1.0
2 21685 1 1 101 1 0 0 0.03238 -0.10595 -0.84531 0.52367 1.0
2 21685 1 2 101 1 0 0 0.01069 -0.38620 0.57732 0.71941 1.0
2 21685 1 3 101 1 0 0 0.00186 -0.91535 -0.35922 -0.18189 1.0
2 21685 1 4 101 8 0 0 0.00322 -0.20489 0.66600 0.71726 1.0
2 21685 1 5 101 8 0 0 0.00462 -0.24681 -0.42086 0.87290 1.0
2 21685 1 6 101 7 0 0 0.00136 -0.06368 -0.41481 0.90768 1.0
2 21685 1 7 101 7 0 0 0.00495 0.63714 -0.28637 -0.71557 1.0
2 21685 1 8 101 7 0 0 0.00138 0.97895 0.20316 0.01943 1.0
3 21685 1 1 101 0 12 24 0.00125 0.16976 0.13859 0.08792 1.0
0
1 21686 1 0 100 1 0 0.099902 0.000003 -0.000007 0.007449 1.0
2 21686 1 1 100 1 0 0 0.09982 -0.07144 -0.12888 0.98908 1.0
3 21686 1 1 100 0 14 28 0.00008 0.03167 0.05650 0.00504 1.0
0
1 22672 1 0 100 1 0 0.099964 0.000001 0.000000 0.004262 1.0
2 22672 1 1 100 1 0 0 0.09929 0.40733 0.06836 0.91072 1.0
3 22672 1 1 100 0 14 28 0.00067 -0.18030 -0.03040 0.04112 1.0
0
1 22866 1 0 101 5 0 0.099899 0.000020 0.000027 0.009746 1.0
2 22866 1 1 101 1 0 0 0.07723 -0.09978 0.53492 0.83899 1.0
2 22866 1 2 101 1 0 0 0.00342 0.76801 -0.52608 0.36525 1.0
2 22866 1 3 101 7 0 0 0.00307 -0.08663 -0.48928 -0.86781 1.0
2 22866 1 4 101 7 0 0 0.00287 -0.15067 0.68749 0.71040 1.0
2 22866 1 5 101 7 0 0 0.00101 0.32893 0.41908 0.84627 1.0
3 22866 1 1 101 0 13 27 0.00071 -0.02001 -0.16536 0.08893 1.0
0
1 23492 1 0 101 5 0 0.099911 0.000004 0.000004 0.009344 1.0
2 23492 1 1 101 1 0 0 0.05999 -0.44674 0.25119 0.85868 1.0
2 23492 1 2 101 8 0 0 0.00161 -0.24584 0.96521 -0.08908 1.0
2 23492 1 3 101 8 0 0 0.01036 -0.65871 0.43470 0.61412 1.0
2 23492 1 4 101 7 0 0 0.00213 0.69593 -0.51986 0.49541 1.0
2 23492 1 5 101 7 0 0 0.00464 -0.28044 -0.93879 0.20008 1.0
3 23492 1 1 101 0 14 28 0.00209 0.25824 -0.19396 0.06870 1.0
0
1 23927 1 0 101 2 0 0.099878 -0.000003 -0.000002 0.003789 1.0
2 23927 1 1 101 1 0 0 0.09753 0.38290 -0.00180 0.92379 1.0
2 23927 1 2 101 7 0 0 0.00177 0.82467 -0.05988 0.56243 1.0
3 23927 1 1 101 0 14 28 0.00058 -0.16985 0.00020 0.03797 1.0
0
1 26900 1 0 100 1 0 0.099917 0.000003 0.000001 0.004228 1.0
2 26900 1 1 100 1 0 0 0.09960 -0.27770 -0.07060 0.95807 1.0
3 26900 1 1 100 0 14 28 0.00032 0.12361 0.03166 0.01935 1.0
0
1 28194 1 0 100 1 0 0.099950 0.000003 -0.000003 0.004031 1.0
2 28194 1 1 100 1 0 0 0.09988 0.09048 0.10427 0.99042 1.0
3 28194 1 1 100 0 14 28 0.00007 -0.03997 -0.04710 0.00443 1.0
0
1 29179 1 0 101 4 0 0.099918 -0.000005 -0.000004 0.007397 1.0
2 29179 1 1 101 1 0 0 0.08242 -0.34161 -0.09185 0.93534 1.0
2 29179 1 2 101 1 0 0 0.00283 0.07267 -0.74201 -0.66644 1.0
2 29179 1 3 101 7 0 0 0.00128 0.48816 -0.21149 -0.84674 1.0
2 29179 1 4 101 7 0 0 0.00102 0.74887 -0.64392 0.15672 1.0
3 29179 1 1 101 0 13 27 0.00078 0.12993 0.09189 0.11811 1.0
0
1 31866 1 0 100 1 0 0.099897 -0.000002 0.000003 0.007697 1.0
2 31866 1 1 100 1 0 0 0.09845 -0.50813 -0.30025 0.80725 1.0
3 31866 1 1 100 0 14 28 0.00145 0.22406 0.13279 0.08839 1.0
0
1 33261 1 0 101 7 0 0.099903 0.000000 0.000007 0.006632 1.0
2 33261 1 1 101 8 0 0 0.04615 -0.02119 -0.22338 0.97450 1.0
2 33261 1 2 101 8 0 0 0.00220 0.82185 0.11257 0.55847 1.0
2 33261 1 3 101 1 0 0 0.01124 0.65957 0.11245 0.74318 1.0
2 33261 1 4 101 1 0 0 0.00555 -0.43488 0.19101 0.88000 1.0
2 33261 1 5 101 7 0 0 0.00237 -0.78045 0.38111 0.49563 1.0
2 33261 1 6 101 7 0 0 0.00056 -0.13234 0.91013 0.39262 1.0
2 33261 1 7 101 7 0 0 0.00084 -0.47657 -0.87908 -0.01024 1.0
3 33261 1 1 101 0 13 27 0.00032 -0.09550 0.02349 -0.08159 1.0
0
1 35977 1 0 101 3 5 0.099975 -0.000001 -0.000001 0.002201 1.0
2 35977 1 1 101 -5 0 0 0.01057 -0.99860 -0.05292 0.00233 1.0
2 35977 1 2 101 8 0 0 0.00358 -0.04751 -0.99541 -0.08304 1.0
2 35977 1 3 101 7 0 0 0.00440 0.28331 0.95784 0.04772 1.0
3 35977 1 1 101 0 0 0 0.00000 0.00000 0.00000 0.00000 1.0
4 35977 1 1 101 3 1 2 0.00313 -0.36632 0.79534 0.48295 1.0
4 35977 1 2 101 6 2 4 0.00498 -0.14248 0.80914 0.57008 1.0
4 35977 1 3 101 6 2 4 0.00477 0.11065 -0.78856 0.60491 1.0
4 35977 1 4 101 6 2 4 0.00453 0.08498 0.98947 0.11713 1.0
4 35977 1 5 101 7 6 12 0.00599 0.75548 -0.50539 0.41693 1.0
0
1 37186 1 0 100 1 0 0.099920 0.000006 0.000010 0.006952 1.0
2 37186 1 1 100 1 0 0 0.09979 -0.16139 -0.07945 0.98369 1.0
3 37186 1 1 100 0 14 28 0.00013 0.07278 0.03602 0.00754 1.0
0
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
0
1 42882 1 0 100 1 0 0.099910 0.000000 0.000006 0.004968 1.0
2 42882 1 1 100 1 0 0 0.09980 -0.15886 -0.05624 0.98570 1.0
3 42882 1 1 100 0 14 28 0.00011 0.07103 0.02583 0.00661 1.0
0
1 43419 1 0 100 1 0 0.099880 0.000002 0.000000 0.004967 1.0
2 43419 1 1 100 1 0 0 0.09962 -0.24953 -0.06284 0.96633 1.0
3 43419 1 1 100 0 14 28 0.00026 0.11099 0.02772 0.01554 1.0
0
1 44325 1 0 100 1 0 0.099913 -0.000001 0.000001 0.002309 1.0
2 44325 1 1 100 1 0 0 0.09956 0.17189 0.24656 0.95376 1.0
3 44325 1 1 100 0 14 28 0.00035 -0.07660 -0.10922 0.02133 1.0
0
1 48143 1 0 100 1 0 0.099986 0.000000 0.000000 0.001022 1.0
2 48143 1 1 100 1 0 0 0.09989 -0.12038 0.10606 0.98705 1.0
3 48143 1 1 100 0 14 28 0.00010 0.05326 -0.04704 0.00598 1.0
0
1 48720 1 0 101 4 0 0.099988 0.000001 -0.000001 0.001348 1.0
2 48720 1 1 101 8 0 0 0.07389 0.03117 0.43401 0.90037 1.0
2 48720 1 2 101 1 0 0 0.00333 -0.56320 -0.82433 -0.05727 1.0
2 48720 1 3 101 7 0 0 0.00343 -0.19207 0.62396 0.75749 1.0
2 48720 1 4 101 7 0 0 0.00172 0.20338 -0.92753 -0.31356 1.0
3 48720 1 1 101 0 14 27 0.00044 0.03318 -0.10062 0.10500 1.0
0
1 48888 1 0 101 4 0 0.099914 -0.000003 0.000006 0.008955 1.0
2 48888 1 1 101 8 0 0 0.02678 -0.41490 0.14409 0.89838 1.0
2 48888 1 2 101 1 0 0 0.02654 0.28656 -0.83065 0.47739 1.0
2 48888 1 3 101 1 0 0 0.01072 0.00133 -0.87566 -0.48293 1.0
2 48888 1 4 101 1 0 0 0.00177 -0.65257 -0.74970 0.11002 1.0
3 48888 1 1 101 0 12 25 0.00316 0.06637 0.32230 0.19649 1.0
0
1 49270 1 0 100 1 0 0.099844 0.000004 -0.000001 0.009826 1.0
2 49270 1 1 100 1 0 0 0.09912 0.42765 -0.02203 0.90368 1.0
3 49270 1 1 100 0 14 28 0.00072 -0.18891 0.00979 0.04432 1.0
0
1 49648 1 0 101 2 4 0.099953 -0.000002 -0.000002 0.003929 1.0
2 49648 1 1 101 1 0 0 0.04577 0.62500 0.41647 0.66025 1.0
2 49648 1 2 101 7 0 0 0.00622 -0.03825 -0.68959 0.72319 1.0
3 49648 1 1 101 0 0 0 0.00000 0.00000 0.00000 0.00000 1.0
4 49648 1 1 101 6 2 4 0.00611 0.54375 0.40127 0.73710 1.0
4 49648 1 2 101 6 2 4 0.00541 -0.07802 0.16945 0.98245 1.0
4 49648 1 3 101 6 2 4 0.00466 -0.19015 0.29717 0.93570 1.0
4 49648 1 4 101 7 8 16 0.00775 -0.52080 -0.61377 -0.59333 1.0
0
1 52709 1 0 101 2 0 0.099988 0.000000 0.000001 0.001187 1.0
2 52709 1 1 101 -3 0 0 0.07732 0.10695 0.49098 0.86458 1.0
2 52709 1 2 101 7 0 0 0.00625 -0.80566 -0.40084 0.43616 1.0
3 52709 1 1 101 0 14 27 0.00146 -0.05317 -0.26410 -0.02857 1.0
0
1 53576 1 0 100 1 0 0.099989 0.000000 0.000000 0.001073 1.0
2 53576 1 1 100 1 0 0 0.09988 0.07641 0.14952 0.98580 1.0
3 53576 1 1 100 0 14 28 0.00011 -0.03380 -0.06636 0.00656 1.0
0
1 53907 1 0 100 1 0 0.099988 0.000000 0.000001 0.001388 1.0
2 53907 1 1 100 1 0 0 0.09590 -0.52475 -0.72225 0.45055 1.0
3 53907 1 1 100 0 14 28 0.00409 0.22839 0.31453 0.24860 1.0
0
1 54217 1 0 101 4 1 0.099955 0.000005 -0.000001 0.003824 1.0
2 54217 1 1 101 1 0 0 0.06276 -0.52915 0.14999 0.83517 1.0
2 54217 1 2 101 1 0 0 0.00387 0.75007 0.20383 -0.62916 1.0
2 54217 1 3 101 7 0 0 0.00199 -0.21242 -0.55050 0.80736 1.0
2 54217 1 4 101 7 0 0 0.00208 0.08570 0.12889 -0.98795 1.0
3 54217 1 1 101 0 11 23 0.00028 0.02723 -0.10687 -0.00270 1.0
4 54217 1 1 101 6 2 4 0.00730 0.40546 0.16380 0.89932 1.0
0
1 55410 1 0 100 1 0 0.099911 0.000003 -0.000002 0.007257 1.0
2 55410 1 1 100 1 0 0 0.09990 0.00813 0.04975 0.99873 1.0
3 55410 1 1 100 0 14 28 0.00001 -0.00344 -0.02243 0.00059 1.0
0
1 56363 1 0 100 1 0 0.099931 0.000001 -0.000003 0.005400 1.0
2 56363 1 1 100 1 0 0 0.09980 -0.12566 0.13970 0.98219 1.0
3 56363 1 1 100 0 14 28 0.00014 0.05578 -0.06244 0.00823 1.0
0
1 56751 1 0 100 1 0 0.099957 -0.000005 -0.000003 0.004578 1.0
2 56751 1 1 100 1 0 0 0.09679 -0.32743 -0.75073 0.57375 1.0
3 56751 1 1 100 0 14 28 0.00317 0.14261 0.32730 0.19373 1.0
0
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
0
1 59278 1 0 100 1 0 0.099894 -0.000002 -0.000009 0.009071 1.0
2 59278 1 1 100 1 0 0 0.09971 0.08489 0.20220 0.97566 1.0
3 59278 1 1 100 0 14 28 0.00019 -0.03748 -0.09055 0.01124 1.0
0
1 59666 1 0 100 1 0 0.099907 0.000004 -0.000001 0.009320 1.0
2 59666 1 1 100 1 0 0 0.09990 -0.01282 0.03425 0.99933 1.0
3 59666 1 1 100 0 14 28 0.00001 0.00598 -0.01548 0.00031 1.0
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
0
1 60726 1 0 101 3 0 0.099995 0.000000 0.000000 0.000522 1.0
2 60726 1 1 101 -3 0 0 0.06367 -0.50810 -0.09974 0.85550 1.0
2 60726 1 2 101 1 0 0 0.00531 -0.52000 0.23111 -0.82230 1.0
2 60726 1 3 101 1 0 0 0.00192 -0.97110 -0.06568 -0.22947 1.0
3 60726 1 1 101 0 12 26 0.00300 0.36056 0.03017 0.11891 1.0
0
1 60766 1 0 100 1 0 0.099915 -0.000010 0.000012 0.006024 1.0
2 60766 1 1 100 1 0 0 0.09988 0.08780 0.00190 0.99614 1.0
3 60766 1 1 100 0 14 28 0.00003 -0.04090 -0.00010 0.00179 1.0
0
1 61605 1 0 100 1 0 0.099854 -0.000002 -0.000007 0.006803 1.0
2 61605 1 1 100 1 0 0 0.09969 -0.20707 0.00844 0.97829 1.0
3 61605 1 1 100 0 14 28 0.00016 0.09138 -0.00438 0.01002 1.0
0
1 63143 1 0 101 4 0 0.099913 0.000001 0.000002 0.006656 1.0
2 63143 1 1 101 1 0 0 0.08791 0.02987 0.01027 0.99950 1.0
2 63143 1 2 101 7 0 0 0.00328 -0.42572 -0.49406 0.75807 1.0
2 63143 1 3 101 7 0 0 0.00693 -0.55801 0.55492 0.61700 1.0
2 63143 1 4 101 7 0 0 0.00178 0.45568 0.83708 0.30274 1.0
3 63143 1 1 101 0 14 28 0.00001 -0.00769 -0.00781 0.02173 1.0
0
1 63251 1 0 100 1 0 0.099868 0.000000 -0.000012 0.007722 1.0
2 63251 1 1 100 1 0 0 0.09958 0.15566 0.22184 0.96258 1.0
3 63251 1 1 100 0 14 28 0.00029 -0.06964 -0.09950 0.01727 1.0
0
1 64229 1 0 100 1 0 0.099975 0.000001 0.000001 0.001739 1.0
2 64229 1 1 100 1 0 0 0.09645 0.32613 0.78537 0.52615 1.0
3 64229 1 1 100 0 14 28 0.00352 -0.14187 -0.34238 0.21500 1.0
0
1 64325 1 0 100 1 0 0.099914 -0.000020 -0.000003 0.008704 1.0
2 64325 1 1 100 1 0 0 0.09961 -0.01469 0.28293 0.95903 1.0
3 64325 1 1 100 0 14 28 0.00031 0.00478 -0.12494 0.01889 1.0
0
1 64380 1 0 100 1 0 0.099935 -0.000006 0.000006 0.005740 1.0
2 64380 1 1 100 1 0 0 0.09992 -0.01909 -0.06460 0.99773 1.0
3 64380 1 1 100 0 14 28 0.00002 0.00814 0.02983 0.00105 1.0
0
1 66544 1 0 100 1 0 0.099901 -0.000003 -0.000005 0.008805 1.0
2 66544 1 1 100 1 0 0 0.09897 -0.47140 0.09838 0.87642 1.0
3 66544 1 1 100 0 14 28 0.00093 0.20796 -0.04379 0.05682 1.0
0
1 67261 1 0 100 1 0 0.099930 -0.000005 0.000000 0.004625 1.0
2 67261 1 1 100 1 0 0 0.09974 0.15151 -0.16000 0.97542 1.0
3 67261 1 1 100 0 14 28 0.00019 -0.06810 0.07068 0.01135 1.0
0
1 70303 1 0 101 2 2 0.099982 0.000001 -0.000001 0.002026 1.0
2 70303 1 1 101 -3 0 0 0.03017 -0.99472 0.02828 0.09869 1.0
2 70303 1 2 101 8 0 0 0.00090 -0.38714 0.36117 -0.84834 1.0
3 70303 1 1 101 0 10 18 0.00497 0.20517 -0.10310 0.33762 1.0
4 70303 1 1 101 6 2 4 0.00795 -0.39524 0.86504 0.30900 1.0
4 70303 1 2 101 6 2 4 0.01041 0.87200 -0.47481 0.11907 1.0
0
1 70574 1 0 100 1 0 0.099839 -0.000003 0.000010 0.007450 1.0
2 70574 1 1 100 1 0 0 0.09939 0.33580 0.02630 0.94157 1.0
3 70574 1 1 100 0 14 28 0.00045 -0.14958 -0.01085 0.02694 1.0
0
1 70993 1 0 100 1 0 0.099962 0.000001 -0.000004 0.004363 1.0
2 70993 1 1 100 1 0 0 0.09973 0.18082 0.16660 0.96930 1.0
3 70993 1 1 100 0 14 28 0.00023 -0.08045 -0.07462 0.01418 1.0
0
1 79181 1 0 100 1 0 0.099902 -0.000005 0.000005 0.008514 1.0
2 79181 1 1 100 1 0 0 0.09973 0.02630 -0.21025 0.97729 1.0
3 79181 1 1 100 0 14 28 0.00017 -0.01187 0.09395 0.01049 1.0
0
1 79748 1 0 100 1 0 0.099962 0.000000 0.000001 0.004335 1.0
2 79748 1 1 100 1 0 0 0.09989 0.07783 -0.11826 0.98993 1.0
3 79748 1 1 100 0 14 28 0.00008 -0.03446 0.05271 0.00465 1.0
0
1 83282 1 0 101 4 1 0.099975 0.000000 0.000001 0.001346 1.0
2 83282 1 1 101 1 0 0 0.04611 -0.85616 0.24649 0.45414 1.0
2 83282 1 2 101 8 0 0 0.00427 0.02307 0.21617 -0.97608 1.0
2 83282 1 3 101 1 0 0 0.00313 -0.30070 -0.43684 0.84780 1.0
2 83282 1 4 101 7 0 0 0.00423 0.95165 -0.05222 -0.30272 1.0
3 83282 1 1 101 0 12 24 0.00478 0.30882 -0.11985 0.32217 1.0
4 83282 1 1 101 3 1 2 0.00140 -0.51156 0.84596 0.15050 1.0
0
1 83585 1 0 100 1 0 0.099862 -0.000011 -0.000001 0.009366 1.0
2 83585 1 1 100 1 0 0 0.09955 -0.21244 0.18919 0.95868 1.0
3 83585 1 1 100 0 14 28 0.00031 0.09343 -0.08442 0.01905 1.0
0
1 84369 1 0 101 3 0 0.099997 0.000000 0.000000 0.000269 1.0
2 84369 1 1 101 1 0 0 0.02468 -0.17059 0.21431 0.96175 1.0
2 84369 1 2 101 8 0 0 0.04124 0.28468 -0.68584 0.66976 1.0
2 84369 1 3 101 1 0 0 0.00881 0.84188 0.51234 0.16954 1.0
3 84369 1 1 101 0 13 26 0.00062 -0.15171 0.08035 0.02591 1.0
0
1 85076 1 0 101 1 3 0.099995 0.000000 0.000000 0.000274 1.0
2 85076 1 1 101 1 0 0 0.05735 -0.19067 0.17728 0.96551 1.0
3 85076 1 1 101 0 8 16 0.00061 -0.10192 -0.01530 0.08628 1.0
4 85076 1 1 101 6 2 4 0.01107 0.25081 -0.78060 0.57250 1.0
4 85076 1 2 101 6 2 4 0.00287 0.55165 0.04875 -0.83265 1.0
4 85076 1 3 101 6 2 4 0.00407 0.07275 0.99676 -0.03444 1.0
0
1 85570 1 0 100 1 0 0.099967 0.000043 0.000046 0.003172 1.0
2 85570 1 1 100 1 0 0 0.09991 0.08259 0.11756 0.98963 1.0
3 85570 1 1 100 0 14 28 0.00005 -0.02191 -0.04709 0.00446 1.0
0
1 85592 1 0 101 5 0 0.099954 0.000001 0.000007 0.004276 1.0
2 85592 1 1 101 1 0 0 0.05275 -0.30805 0.59923 0.73894 1.0
2 85592 1 2 101 8 0 0 0.01738 0.60411 -0.78893 0.11243 1.0
2 85592 1 3 101 8 0 0 0.00349 0.06374 -0.94513 -0.32042 1.0
2 85592 1 4 101 7 0 0 0.00453 -0.49717 -0.04714 -0.86637 1.0
2 85592 1 5 101 7 0 0 0.00178 0.38515 -0.92244 0.02758 1.0
3 85592 1 1 101 0 14 28 0.00094 -0.01416 0.03156 0.21818 1.0
0
1 87763 1 0 101 5 4 0.099911 -0.000006 0.000007 0.006482 1.0
2 87763 1 1 101 1 0 0 0.01180 -0.81744 -0.06820 0.57196 1.0
2 87763 1 2 101 1 0 0 0.00523 0.41158 0.49594 0.76462 1.0
2 87763 1 3 101 1 0 0 0.00231 -0.40248 -0.74337 -0.53423 1.0
2 87763 1 4 101 8 0 0 0.00257 0.74442 0.04321 0.66632 1.0
2 87763 1 5 101 8 0 0 0.00201 -0.11781 -0.46219 -0.87892 1.0
3 87763 1 1 101 0 0 0 0.00000 0.00000 0.00000 0.00000 1.0
4 87763 1 1 101 6 2 4 0.00814 0.33652 0.60069 0.72521 1.0
4 87763 1 2 101 6 2 4 0.00285 0.86943 0.08458 0.48676 1.0
4 87763 1 3 101 6 2 4 0.00215 -0.28588 0.30992 -0.90676 1.0
4 87763 1 4 101 7 6 12 0.00337 -0.40259 -0.59441 0.69613 1.0
0
1 89101 1 0 101 4 0 0.099924 -0.000003 -0.000005 0.006639 1.0
2 89101 1 1 101 8 0 0 0.07387 -0.19240 -0.26553 0.94471 1.0
2 89101 1 2 101 1 0 0 0.00452 0.27503 -0.56215 0.77996 1.0
2 89101 1 3 101 7 0 0 0.00259 -0.23971 0.23499 0.94198 1.0
2 89101 1 4 101 7 0 0 0.00126 0.51012 -0.83044 0.22395 1.0
3 89101 1 1 101 0 14 27 0.00051 0.04712 0.15273 0.01097 1.0
0
1 91149 1 0 100 1 0 0.099925 -0.000004 0.000005 0.007421 1.0
2 91149 1 1 100 1 0 0 0.09791 -0.50880 -0.45210 0.73262 1.0
3 91149 1 1 100 0 14 28 0.00201 0.22365 0.19946 0.12228 1.0
0
1 91270 1 0 101 6 0 0.099922 0.000001 -0.000006 0.005168 1.0
2 91270 1 1 101 1 0 0 0.00608 -0.23811 -0.16551 0.95703 1.0
2 91270 1 2 101 8 0 0 0.03837 -0.86279 -0.29294 -0.41203 1.0
2 91270 1 3 101 8 0 0 0.00043 -0.48744 -0.01453 -0.87304 1.0
2 91270 1 4 101 1 0 0 0.00352 -0.30922 -0.45890 0.83294 1.0
2 91270 1 5 101 1 0 0 0.00459 -0.17667 -0.98013 -0.09016 1.0
2 91270 1 6 101 7 0 0 0.00135 -0.65787 -0.72098 0.21771 1.0
3 91270 1 1 101 0 12 24 0.00730 0.31544 0.22616 0.41881 1.0
0
1 96245 1 0 100 1 0 0.099949 -0.000002 -0.000001 0.001835 1.0
2 96245 1 1 100 1 0 0 0.09839 -0.24462 0.55672 0.79387 1.0
3 96245 1 1 100 0 14 28 0.00156 0.10717 -0.24638 0.09451 1.0
0
1 96743 1 0 101 4 1 0.099871 0.000038 -0.000004 0.009717 1.0
2 96743 1 1 101 1 0 0 0.04361 0.28910 0.17720 0.94075 1.0
2 96743 1 2 101 -3 0 0 0.01190 0.34579 0.56928 0.74589 1.0
2 96743 1 3 101 7 0 0 0.00102 0.93628 -0.11080 -0.33331 1.0
2 96743 1 4 101 7 0 0 0.00058 -0.92434 0.23326 0.30196 1.0
3 96743 1 1 101 0 11 22 0.00167 0.07220 -0.20794 0.14132 1.0
4 96743 1 1 101 6 2 4 0.00922 -0.86512 0.13150 -0.48402 1.0
0
1 97535 1 0 100 1 0 0.099987 0.000000 0.000001 0.001551 1.0
2 97535 1 1 100 1 0 0 0.09998 0.03868 0.01349 0.99916 1.0
3 97535 1 1 100 0 14 28 0.00001 -0.01758 -0.00570 0.00039 1.0
0
1 98589 1 0 101 5 0 0.099907 0.000003 0.000001 0.008832 1.0
2 98589 1 1 101 8 0 0 0.05633 -0.36591 0.66559 0.65046 1.0
2 98589 1 2 101 1 0 0 0.00178 -0.98536 0.16978 -0.01575 1.0
2 98589 1 3 101 1 0 0 0.00368 0.23982 -0.94263 -0.23224 1.0
2 98589 1 4 101 1 0 0 0.00482 -0.06670 -0.96810 0.24153 1.0
2 98589 1 5 101 7 0 0 0.00059 0.96081 0.22374 -0.16366 1.0
3 98589 1 1 101 0 12 25 0.00176 0.16389 -0.05911 0.22692 1.0
0
*$ CREATE MGDRAW.FOR
*COPY MGDRAW
* *
*=== mgdraw ===========================================================*
* *
SUBROUTINE MGDRAW ( ICODE, MREG )
INCLUDE '(DBLPRC)'
INCLUDE '(DIMPAR)'
INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2006 by Alfredo Ferrari *
* All Rights Reserved. *
* *
* *
* MaGnetic field trajectory DRAWing: actually this entry manages *
* all trajectory dumping for *
* drawing *
* *
* Created on 01 march 1990 by Alfredo Ferrari *
* INFN - Milan *
* Last change 05-may-06 by Alfredo Ferrari *
* INFN - Milan *
* *
*----------------------------------------------------------------------*
*
INCLUDE '(CASLIM)'
INCLUDE '(COMPUT)'
INCLUDE '(SOURCM)'
INCLUDE '(FHEAVY)'
INCLUDE '(FLKSTK)'
INCLUDE '(GENSTK)'
INCLUDE '(MGDDCM)'
INCLUDE '(PAPROP)'
INCLUDE '(QUEMGD)'
INCLUDE '(SUMCOU)'
INCLUDE '(TRACKR)'
INCLUDE '(RESNUC)'
*
DIMENSION DTQUEN ( MXTRCK, MAXQMG )
*
CHARACTER*20 FILNAM
LOGICAL LFCOPE
SAVE LFCOPE
DATA LFCOPE / .FALSE. /
INTEGER ZPART !default Z
INTEGER APART !default A
INTEGER WPART !default weight
*------------------------------------------
*Convert region name to number
INTEGER NTARG1
INTEGER NTARG2
CHARACTER*8 NAMET1
CHARACTER*8 NAMET2
SAVE NTARG1
SAVE NTARG2
PARAMETER (NAMET1 = 'TARG1')
PARAMETER (NAMET2 = 'TARG2')
* Open files for usdraw
LOGICAL UDFIRST
DATA UDFIRST / .TRUE. /
SAVE UDFIRST
LOGICAL LFIRST
DATA LFIRST / .TRUE. /
SAVE LFIRST
IF ( LFIRST ) THEN
ZPART = 0
APART = 0
WPART = 0
* +------------------------------------------------------------*
* |Convert region name to number
CALL GEON2R(NAMET1, NTARG1,IERR)
WRITE(LUNOUT,*) 'NTARG1= ', NTARG1
CALL GEON2R(NAMET2,NTARG2,IERR)
LFIRST = .FALSE.
END IF
*
*----------------------------------------------------------------------*
* *
* Icode = 1: call from Kaskad *
* Icode = 2: call from Emfsco *
* Icode = 3: call from Kasneu *
* Icode = 4: call from Kashea *
* Icode = 5: call from Kasoph *
* *
*----------------------------------------------------------------------*
* *
*
* IF ( .NOT. LFCOPE ) THEN
* LFCOPE = .TRUE.
* IF ( KOMPUT .EQ. 2 ) THEN
* FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
* ELSE
* FILNAM = CFDRAW
* END IF
* OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
* & 'UNFORMATTED' )
* END IF
* WRITE (IODRAW) NTRACK, MTRACK, JTRACK, SNGL (ETRACK),
* & SNGL (WTRACK)
* WRITE (IODRAW) ( SNGL (XTRACK (I)), SNGL (YTRACK (I)),
* & SNGL (ZTRACK (I)), I = 0, NTRACK ),
* & ( SNGL (DTRACK (I)), I = 1, MTRACK ),
* & SNGL (CTRACK)
* +-------------------------------------------------------------------*
* | Quenching is activated
IF ( LQEMGD ) THEN
IF ( MTRACK .GT. 0 ) THEN
RULLL = ZERZER
CALL QUENMG ( ICODE, MREG, RULLL, DTQUEN )
* WRITE (IODRAW) ( ( SNGL (DTQUEN (I,JBK)), I = 1, MTRACK ),
* & JBK = 1, NQEMGD )
END IF
END IF
* | End of quenching
* +-------------------------------------------------------------------*
RETURN
*
*======================================================================*
* *
* Boundary-(X)crossing DRAWing: *
* *
* Icode = 1x: call from Kaskad *
* 19: boundary crossing *
* Icode = 2x: call from Emfsco *
* 29: boundary crossing *
* Icode = 3x: call from Kasneu *
* 39: boundary crossing *
* Icode = 4x: call from Kashea *
* 49: boundary crossing *
* Icode = 5x: call from Kasoph *
* 59: boundary crossing *
* *
*======================================================================*
* *
ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
IF (DEBUG.EQ.1) THEN
* WRITE(LUNOUT,*) '************DEBUG***************'
WRITE(LUNOUT,*) 'NTARG1= ', NTARG1, ' MREG= ', MREG
WRITE(LUNOUT,*) 'NTARG2= ', NTARG2, ' NEWREG= ', NEWREG
* WRITE(LUNOUT,*) '************DEBUG END***************'
END IF
IF(MREG.EQ.NTARG1.AND.NEWREG.EQ.NTARG2) THEN ! select boundary of interest
IF(JTRACK.EQ.1) THEN ! select proton
IF(ETRACK.GT.AM(JTRACK)) THEN ! proton has survived
WRITE(41.0,*) ETRACK-AM(JTRACK),CXTRCK,CYRCK,CZTRCK,WTRACK,NCASE !,XSCO,YSCO,ZSCO,
ENDIF
ENDIF
ENDIF
IF(MREG.EQ.NTARG1.AND.NEWREG.EQ.NTARG2) THEN
IF(JTRACK.EQ.8) THEN ! select neutron
IF(ETRACK.GT.AM(JTRACK)) THEN
WRITE(42.0,*) ETRACK-AM(JTRACK),XSCO,YSCO,ZSCO,NCASE !CXTRCK,CYRCK,CZTRCK,WTRACK,KTRACK,NCASE
ENDIF
ENDIF
ENDIF
RETURN
*
*======================================================================*
* *
* Event End DRAWing: *
* *
*======================================================================*
* *
ENTRY EEDRAW ( ICODE )
RETURN
*
*======================================================================*
* *
* ENergy deposition DRAWing: *
* *
* Icode = 1x: call from Kaskad *
* 10: elastic interaction recoil *
* 11: inelastic interaction recoil *
* 12: stopping particle *
* 13: pseudo-neutron deposition *
* 14: escape *
* 15: time kill *
* Icode = 2x: call from Emfsco *
* 20: local energy deposition (i.e. photoelectric) *
* 21: below threshold, iarg=1 *
* 22: below threshold, iarg=2 *
* 23: escape *
* 24: time kill *
* Icode = 3x: call from Kasneu *
* 30: target recoil *
* 31: below threshold *
* 32: escape *
* 33: time kill *
* Icode = 4x: call from Kashea *
* 40: escape *
* 41: time kill *
* 42: delta ray stack overflow *
* Icode = 5x: call from Kasoph *
* 50: optical photon absorption *
* 51: escape *
* 52: time kill *
* *
*======================================================================*
* *
ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
* IF ( .NOT. LFCOPE ) THEN
* LFCOPE = .TRUE.
* IF ( KOMPUT .EQ. 2 ) THEN
* FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
* ELSE
* FILNAM = CFDRAW
* END IF
* OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
* & 'UNFORMATTED' )
* END IF
* WRITE (IODRAW) 0, ICODE, JTRACK, SNGL (ETRACK), SNGL (WTRACK)
* WRITE (IODRAW) SNGL (XSCO), SNGL (YSCO), SNGL (ZSCO), SNGL (RULL)
* +-------------------------------------------------------------------*
* | Quenching is activated : calculate quenching factor
* | and store quenched energy in DTQUEN(1, jbk)
* IF ( LQEMGD ) THEN
* RULLL = RULL
* CALL QUENMG ( ICODE, MREG, RULLL, DTQUEN )
* WRITE (IODRAW) ( SNGL (DTQUEN(1, JBK)), JBK = 1, NQEMGD )
* END IF
* | end quenching
* +-------------------------------------------------------------------*
RETURN
*
*======================================================================*
* *
* SOurce particle DRAWing: *
* *
*======================================================================*
*
ENTRY SODRAW
* IF ( .NOT. LFCOPE ) THEN
* LFCOPE = .TRUE.
* IF ( KOMPUT .EQ. 2 ) THEN
* FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
* ELSE
* FILNAM = CFDRAW
* END IF
* OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
* & 'UNFORMATTED' )
* END IF
* WRITE (IODRAW) -NCASE, NPFLKA, NSTMAX, SNGL (TKESUM),
* & SNGL (WEIPRI)
* +-------------------------------------------------------------------*
* | (Radioactive) isotope: it works only for 1 source particle on
* | the stack for the time being
* IF ( ILOFLK (NPFLKA) .GE. 100000 .AND. LRADDC (NPFLKA) ) THEN
* IARES = MOD ( ILOFLK (NPFLKA), 100000 ) / 100
* IZRES = MOD ( ILOFLK (NPFLKA), 10000000 ) / 100000
* IISRES = ILOFLK (NPFLKA) / 10000000
* IONID = ILOFLK (NPFLKA)
* WRITE (IODRAW) ( IONID,SNGL(-TKEFLK(I)),
* & SNGL (WTFLK(I)), SNGL (XFLK (I)),
* & SNGL (YFLK (I)), SNGL (ZFLK (I)),
* & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
* & SNGL (TZFLK(I)), I = 1, NPFLKA )
* |
* +-------------------------------------------------------------------*
* | Patch for heavy ions: it works only for 1 source particle on
* | the stack for the time being
* ELSE IF ( ABS (ILOFLK (NPFLKA)) .GE. 10000 ) THEN
* IONID = ILOFLK (NPFLKA)
* CALL DCDION ( IONID )
* WRITE (IODRAW) ( IONID,SNGL(TKEFLK(I)+AMNHEA(-IONID)),
* & SNGL (WTFLK(I)), SNGL (XFLK (I)),
* & SNGL (YFLK (I)), SNGL (ZFLK (I)),
* & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
* & SNGL (TZFLK(I)), I = 1, NPFLKA )
* |
* +-------------------------------------------------------------------*
* | Patch for heavy ions: ???
* ELSE IF ( ILOFLK (NPFLKA) .LT. -6 ) THEN
* WRITE (IODRAW) ( IONID,SNGL(TKEFLK(I)+AMNHEA(-ILOFLK(NPFLKA))),
* & SNGL (WTFLK(I)), SNGL (XFLK (I)),
* & SNGL (YFLK (I)), SNGL (ZFLK (I)),
* & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
* & SNGL (TZFLK(I)), I = 1, NPFLKA )
* |
* +-------------------------------------------------------------------*
* |
* ELSE
* WRITE (IODRAW) ( ILOFLK(I), SNGL (TKEFLK(I)+AM(ILOFLK(I))),
* & SNGL (WTFLK(I)), SNGL (XFLK (I)),
* & SNGL (YFLK (I)), SNGL (ZFLK (I)),
* & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
* & SNGL (TZFLK(I)), I = 1, NPFLKA )
* END IF
* |
* +-------------------------------------------------------------------*
RETURN
*
*======================================================================*
* *
* USer dependent DRAWing: *
* *
* Icode = 10x: call from Kaskad *
* 100: elastic interaction secondaries *
* 101: inelastic interaction secondaries *
* 102: particle decay secondaries *
* 103: delta ray generation secondaries *
* 104: pair production secondaries *
* 105: bremsstrahlung secondaries *
* 110: decay products *
* Icode = 20x: call from Emfsco *
* 208: bremsstrahlung secondaries *
* 210: Moller secondaries *
* 212: Bhabha secondaries *
* 214: in-flight annihilation secondaries *
* 215: annihilation at rest secondaries *
* 217: pair production secondaries *
* 219: Compton scattering secondaries *
* 221: photoelectric secondaries *
* 225: Rayleigh scattering secondaries *
* Icode = 30x: call from Kasneu *
* 300: interaction secondaries *
* Icode = 40x: call from Kashea *
* 400: delta ray generation secondaries *
* For all interactions secondaries are put on GENSTK common (kp=1,np) *
* but for KASHEA delta ray generation where only the secondary elec- *
* tron is present and stacked on FLKSTK common for kp=npflka *
* *
*======================================================================*
*
ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
* IF ( .NOT. LFCOPE ) THEN
* LFCOPE = .TRUE.
* IF ( KOMPUT .EQ. 2 ) THEN
* FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
* ELSE
* FILNAM = CFDRAW
* END IF
* OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
* & 'UNFORMATTED' )
* END IF
IF (UDFIRST) THEN
* OPEN(43,FILE='primary.log',STATUS='UNKNOWN')
OPEN(44,FILE='secondary.log',STATUS='UNKNOWN')
UDFIRST = .FALSE.
ENDIF
* No output by default:
* |
* +-------------------------------------------------------------------*
* |DUMP#1 PRIMARY PARTICLE INFO FROM (TRACK)
* |
* |Dump : |DUMPID|NCASE|JTRACK|NPFLKA|ICODE|ETRACK-AM(JTRACK)|XSCO|-
* | |YSCO|ZSCO|NP|NPHEAVY|WTRACK|
WRITE( 44.0, '(A2,I5.5,1X,I2,1X,I1,1X,I3,1X,I1,1X,I1,1X,
&F8.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F3.1,I1)') '1 ',NCASE,JTRACK,
&NPFLKA,ICODE,NP,NPHEAV,ETRACK-AM(JTRACK),XSCO,YSCO,ZSCO,WTRACK
* WRITE( 43.0, * ) '0 ',NCASE,JTRACK,NPFLKA,
* &ICODE,ETRACK-AM(JTRACK),XSCO,YSCO,ZSCO,NP,NPHEAV,WTRACK
* WRITE( 43.0, * ) '1 ', NCASE,NPFLKA ,ICODE, JTRACK,WTRACK, NP
* & ETRACK-AM(JTRACK),
* |
* +-------------------------------------------------------------------*
* |DUMP#2 LOOP OVER SECONDARIES (GENSTK)
* |For the moment the Z and A are 0. Taken care of in post processing
* |
* |* |Dump : |DUMPID|NCASE|JTRACK|ip|ICODE|KPART(ip)|ZPART|APART|-
* | |TKI(ip)|CXR(ip)|CYR(ip)|CZR(ip)|WEI(ip)|
do 10 ip = 1, NP
WRITE( 44.0, '(A2,I5.5,1X,I2,1X,I1,1X,I3,1X,I2,1X,I2,1X,I2,
&1X,F7.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F3.1)')'2 ',NCASE,JTRACK,
&ip,ICODE,KPART(ip),ZPART,APART,TKI(ip),CXR(ip),CYR(ip),
&CZR(ip),WEI(ip)
10 continue
* |
* +-------------------------------------------------------------------*
* |DUMP#3 RESIDUAL NUCLEUS (RESNUC)
* |
* |Position 4: '1' is replacing the ip counter since only 1 residual
* |nucleus will be produce from RESNUC
* |Position 6: '0' is replacing KPART and KHEAVY
* |Could not find a weighting factor in (RESNUC) so last position is
* |default '1.0'
* |
* IF (ICRES.NE.0.and.IBRES.NE.0) THEN
WRITE( 44.0, '(A2,I5.5,1X,I2,1X,A1,1X,I3,1X,A2,1X,I2,1X,I2,
&1X,F7.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,A3)')'3 ',NCASE,JTRACK,
&'1',ICODE,'0',ICRES,IBRES,EKRES,PXRES,PYRES,PZRES,'1.0'
* END IF
* |
* +-------------------------------------------------------------------*
* |DUMP#4 LOOP OVER HEAVY SECONDARIES (FHEAVY)
* |
do 20 ip = 1, NPHEAV
WRITE( 44.0, '(A2,I5.5,1X,I2,1X,I1,1X,I3,1X,I2,1X,I2,1X,I2,
&1X,F7.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F3.1)')'4 ',NCASE,JTRACK,
&ip,ICODE,KHEAVY(ip),ICHEAV(KHEAVY(ip)),IBHEAV(KHEAVY(ip)),
&TKHEAV(ip),CXHEAV(ip),CYHEAV(ip),CZHEAV(ip),WHEAVY(ip)
20 continue
WRITE(44.0, * ) '0' !output 0 for easy reading
RETURN
*=== End of subrutine Mgdraw ==========================================*
END
TITLE
Proton beam
* FLAIR: Use names
GLOBAL 1.0 1.0
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
DEFAULTS NEW-DEFA
BEAM -0.1 -0.0 0.0 0.0 0.0 -1.0PROTON
BEAMPOS 0.0 0.0 -0.5 0.0 0.0
EVENTYPE 2.0 EVAP
PHYSICS 3.0 EVAPORAT
PHYSICS 1.0 COALESCE
GEOBEGIN COMBNAME
0 0 Cylindrical Target
* *AAA*IIII_________+_________+_________+_________+_________+_________+
SPH BLK 0.0 0.0 0.0 10000.0
* vacuum box
RPP VOI -1000.0 1000.0 -1000.0 1000.0 -1000.0 1000.0
* target
ZCC TARG 0.0 0.0 10.0
* planes limiting the target
XYP START 0.0
XYP END 0.01
YZP TOP 5.0
YZP BOTTOM -5.0
XZP LEFT -5.0
XZP RIGHT 5.0
* Dummy region before target (for boundary scoring)
XYP P1 -0.1
* Dummy region after target (for boundary scoring)
XYP P2 1.1
END
* Regions
* _AAA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....AA.....
* Black Hole
BLKHOLE 5 +BLK -VOI
* Dummy
DUMMY1 5 +TARG -P1 +START -LEFT +RIGHT -BOTTOM +TOP
* Target1
TARG1 5 +TARG -START +END -LEFT +RIGHT -BOTTOM +TOP
* Target2
TARG2 5 +TARG +P2 -END -LEFT +RIGHT -BOTTOM +TOP
* Vac
VAC 5 +VOI -( +TARG -P1 +P2 -LEFT +RIGHT -BOTTOM +TOP )
END
* AAAAAAAxx_________+_________+_________+_________+_________+_________+
GEOEND
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA BLCKHOLE BLKHOLE
ASSIGNMA VACUUM DUMMY1
ASSIGNMA SILICON TARG1
ASSIGNMA VACUUM TARG2
ASSIGNMA VACUUM VAC
PART-THR -0.0001 PROTON PROTON 0.0
USERDUMP 100.0 37.0 0.0 1.0
RESNUCLE 3.0 -26.0 TARG1 1.0Resnuclei
START 100000.0 0.0
OPEN 41.0 NEW
bdx_proton4
OPEN 42.0 NEW
bdx_neutron4
#if 0
OPEN 43.0 NEW
primary
#endif
#if 0
OPEN -44.0 NEW
secondary_bin
#endif
#if 0
OPEN 45.0 NEW
nucleus
#endif
#if 0
OPEN 46.0 NEW
heavy
#endif
#if 0
OPEN 47.0 NEW
summary
#endif
#if 0
OPEN 48.0 NEW
mgdraw_call
#endif
RANDOMIZ 1.0
STOP
This archive was generated by hypermail 2.1.6 : Fri Jul 06 2007 - 10:06:36 CEST