Re: [fluka-discuss]: using fluka as an event generator, where are the ions?

From: <baret_at_apc.univ-paris7.fr>
Date: Wed, 22 Aug 2018 10:21:40 +0200

Hi and thanks!
concerning the User Liscence , I want to generate events p/gamma+Nuclei over a subGeV-TeV range for Cosmic Ray acceleration and propagation studies. It indeed uses a home made code but no piece of Fluka code and would not make accessible the details of the underlying codes (I think!) and has already been used with outputs of some of the underlying hadronic codes within Fluka. I don’t think this goes out of the scope of the licence but I would of course ask written permission if you think it’s necessary. After some search and discussions with colleagues, it appeared to me that fluka is the most (only?) reliable software mixing all the interaction models for this kind of purpose...
So comes my next question, is there a way to produce forced interaction without propagation of secondaries but with storage of their characteristics?
So far I’ve been trying to get those by playing with found simple examples of beam+target.
Thanks to your post I can now retrieve secondary nuclei with fluka ids=-9…-1 (by the way I’m not sure to have understood from the manual Chap5 the -9…-7 ones ) .
Given the unsatifactory way I do this, I naively at some point have put a high density of Fe to boost interaction rate. Of course it also make interactions of secondaries and so on (useless for me) more likely. I have retrieved some wanted interactions by selecting LTRACK=2 but that’s very inefficient way to say the least.
Would you have a hint?
Cheers,
Bruny


From: Francesc Salvat-Pujol <francesc.salvat.pujol_at_cern.ch <mailto:francesc.salvat.pujol_at_cern.ch?Subject=Re%3A%20%5Bfluka-discuss%5D%3A%20using%20fluka%20as%20an%20event%20generator%2C%20where%20are%20the%20ions%3F>>
Date: Fri, 17 Aug 2018 09:48:58 +0200

Hi,

Regarding your intended use, please observe paragraphs 1 and 6 of the
FLUKA User License.

Regarding the input, it could be a display or copy-paste artifact, but
the SDUMs of the three PHYSICS cards appear misaligned. Ensure column
alignment.

In your pasted output section one sees light ions (4-HELIUM, 3-HELIUM,
TRITON, DEUTERON), but no HEAVYION. It could be that the latter are
generated below transport threshold. In any event, once you are sure of
the input, you can examine the heavy ion secondaries with e.g. mdstck.f
using the following snippet:

      DO I = 1 , NP
         IONID = KPART(I)
         IF (KPART(I) .LT. -6) THEN
            IAHELP = MOD ( ABS(IONID), 100000 ) / 100
            IZHELP = MOD ( ABS(IONID), 10000000 ) / 100000
* And then write out to whichever unit you defined.
         END IF
      ENDDO

With kind regards,

Cesc

PS: WHAT(3) (density in g/cm**3) of the MATERIAL card for IRON appears rather high.
Note that IRON is a pre-defined material.

On Thu, Aug 16 2018, at 16:29 +0200, baret_at_apc.univ-paris7.fr wrote:
>
>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
>
>
>

--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Aug 22 2018 - 11:51:39 CEST

This archive was generated by hypermail 2.3.0 : Wed Aug 22 2018 - 11:51:51 CEST