Dear fluka experts and users,
now that experts have unstuck me on the hadronic version linking of user routines I can ask another question.
I want to use fluka as event generator for photon +hadron and hadron +hadron above 100GeV.
Actually I would like to retrieve all the interactions and their secondaries.
To do so I modified mdstk.f to print the variables of GENSTK at each interaction:
WRITE (IO,*) NP,NPHEAV, IFLAG, NPSECN, NTRACK, MTRACK,
& JTRACK, SNGL (ETRACK), SNGL (WTRACK)
WRITE (IO,*) ( SNGL (XTRACK (I)), SNGL (YTRACK (I)),
& SNGL (ZTRACK (I)), I = 0, NTRACK ),
& ( SNGL (DTRACK (I)), I = 1, MTRACK ),
& SNGL (CTRACK)
I am not sure at all it is the best way to do it but it seem to work except that I don’t find any ion in the output file I create even if the .out file shows that some are produced:
Number of secondaries generated in inelastic interactions per beam particle:
Prompt radiation Radioactive decays
8.4490E+02 (100.%) 0.0000E+00 (100.%)
3.3900E+01 ( 4.0%) 0.0000E+00 ( 0.0%) 4-HELIUM
3.6000E+00 ( 0.4%) 0.0000E+00 ( 0.0%) 3-HELIUM
5.9000E+00 ( 0.7%) 0.0000E+00 ( 0.0%) TRITON
1.4600E+01 ( 1.7%) 0.0000E+00 ( 0.0%) DEUTERON
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) HEAVYION
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) OPTIPHOT
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) RAY
1.5820E+02 (18.7%) 0.0000E+00 ( 0.0%) PROTON
1.0000E-01 ( 0.0%) 0.0000E+00 ( 0.0%) APROTON
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) ELECTRON
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) POSITRON
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) NEUTRIE
0.0000E+00 ( 0.0%) 0.0000E+00 ( 0.0%) ANEUTRIE
2.4200E+02 (28.6%) 0.0000E+00 ( 0.0%) PHOTON
3.0730E+02 (36.4%) 0.0000E+00 ( 0.0%) NEUTRON
and so on…
I think I activated all necessay options but can’t be 100% sure
The input file is:
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
Charged pion fluence inside and around a proton-irradiated Fe target
PHOTONUC 1111.
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
BEAM 3.E+02 PROTON
BEAMPOS 0.0 0.0 -50.
*
GEOBEGIN COMBNAME
0 0 A simple Fe target inside vacuum
RPP body1 -5000000. 5000000. -5000000. 5000000. -5000000. 5000000.
RPP body2 -1000000. 1000000. -1000000. 1000000. -100. 1000000.
RPP body3 -10. 10. -10. 10. 0.0 5.
* plane to separate the upstream and downstream part of the target
XYP body4 2.5
END
* black hole
regBH1 11 +body1 -body2
* vacuum around
regVA2 11 +body2 -body3
* Be target 1st half
regBE3 11 +body3 +body4
* Be target 2nd half
regBE4 11 +body3 -body4
END
GEOEND
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
MATERIAL 26. 55.845 52.1 11. IRON
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Be target, 1st and 2nd half
ASSIGNMA IRON regBE3 regBE4
* External Black Hole
ASSIGNMA BLCKHOLE regBH1
* Vacuum
ASSIGNMA VACUUM regVA2
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* e+e- and gamma production threshold set at 10 MeV
EMFCUT -0.010 0.010 1. IRON PROD-CUT
* score in each region energy deposition and stars produced by primaries
SCORE ENERGY BEAMPART
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
IONTRANS -2.0
PHYSICS 2.0 0.0 0.0 0.0 0.0 0.0 EM-DISSO
PHYSICS 3.0 0.0 0.0 0.0 0.0 0.0 EVAPORAT
PHYSICS 1.0 0.0 0.0 0.0 0.0 0.0 COALESC
* opening the file with the random number seeds for the next run
*OPEN 37.
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
USERDUMP 100. 37. 1. MGDRAW
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
RANDOMIZ 1.
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
START 10.0
STOP
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Aug 16 2018 - 17:52:27 CEST