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

From: <baret_at_apc.univ-paris7.fr>
Date: Thu, 20 Sep 2018 10:41:02 +0200

thanks for your help!
Now I now use the LAM-BIAS option and can retrieve secondary ions (but too many it seems to me) thanks to the code snipset suggested in mdstck.f by Francesc, :
> 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

 
 and not tracking secondaries by adding to usrmed.f


    INCLUDE '(TRACKR)'

    IF (LTRACK .NE. 1) THEN
       WEE = 0.0+D0
    END IF

(I checked that all particle that I get have AGESEC=0)
I get now a range of masses of secondary ions but even too much I think. The IAHELP and ZHELP above when computed with the simple .inp file attached, show ion with A>56 and that puzzles me (see attached pdf for 10 000 events). Am I wrong ?





> Le 23 août 2018 à 18:44, Anton Lechner <Anton.Lechner_at_cern.ch> a écrit :
>
> Hi,
>
> One can always increase the interaction rate in a thin target by biasing the interaction length (see LAM-BIAS). This is much more effective than increasing the density.
>
> Concerning the license, as already mentioned by Cesc, one should consider the relevant paragraphs.
>
> Cheers, Anton
>
>
>
> From: owner-fluka-discuss_at_mi.infn.it <mailto:owner-fluka-discuss_at_mi.infn.it> [owner-fluka-discuss_at_mi.infn.it <mailto:owner-fluka-discuss_at_mi.infn.it>] on behalf of baret_at_apc.univ-paris7.fr <mailto:baret_at_apc.univ-paris7.fr> [baret_at_apc.univ-paris7.fr <mailto:baret_at_apc.univ-paris7.fr>]
> Sent: 22 August 2018 10:21
> To: fluka-discuss_at_fluka.org <mailto:fluka-discuss_at_fluka.org>
> Subject: Re: [fluka-discuss]: using fluka as an event generator, where are the ions?
>
> 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 <http://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 Thu Sep 20 2018 - 12:15:15 CEST

This archive was generated by hypermail 2.3.0 : Thu Sep 20 2018 - 12:15:34 CEST