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

From: George Kharashvili <georgek_at_jlab.org>
Date: Sat, 22 Sep 2018 05:42:39 -0400 (EDT)

Hello Bruny,

A=59 can be a product of protons incident on Fe-58.
I suppose you are correct about the AGESEC=0 criteria for your problem, but I'm not absolutely sure. May be someone else can comment on this.

George


----- Original Message -----
From: baret_at_apc.univ-paris7.fr
To: "George Kharashvili" <georgek_at_jlab.org>
Cc: "FLUKA Discussion List" <fluka-discuss_at_fluka.org>
Sent: Friday, September 21, 2018 7:35:47 PM
Subject: Re: [fluka-discuss]: using fluka as an event generator, where are the ions?

Dear George,
thanks for the input. Indeed, I didn’t put the right inp file (the right one is attached)… Sorry for that.
Indeed if there are other isotopes that might explain some heavier stuff than Fe56 that could explain most of them. However I’m still intrigued by the A=59 since I’ve normally set things so that secondaries don’t interact after creation with the usrmed.f modification, unless I didn’t do it properly.
Particles that can interact after creation should at a point have a AGESEC variable with non zero value right? So far, all mine have AGESEC=0.


Cheers,
Bruny







> Le 21 sept. 2018 à 06:52, George Kharashvili <georgek_at_jlab.org> a écrit :
>
> Dear Bruny,
>
> The input file does not contain the LAM-BIAS card, so I wonder if you meant to send a different input.
> In any case, 300 GeV photon beam will certainly produce a large number of fragments, some of which may reinteract and produce nuclei with higher mass number than the target. Also, keep in mind that you are using natural iron as your target, which contains some Fe-57 and Fe-58, in addition to Fe-56 and Fe-54.
>
> Cheers,
> George
>
> ----- Original Message -----
> From: baret_at_apc.univ-paris7.fr
> To: "Anton Lechner" <Anton.Lechner_at_cern.ch>
> Cc: "FLUKA Discussion List" <fluka-discuss_at_fluka.org>
> Sent: Thursday, September 20, 2018 12:41:02 PM
> Subject: Re: [fluka-discuss]: using fluka as an event generator, where are the ions?
>
> 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 <https://urldefense.proofpoint.com/v2/url?u=http-3A__baret-5Fat-5Fapc.univ-2Dparis7.fr_&d=DwIFaQ&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=VEsjWcjHR_RwJLerQJyK4g&m=ZVvO6vStNXpBSohLqAcCrjzQnMkPhQgVHacnEQsVmmA&s=9zglRe7QEN9lgIa4wBLDc97lvliBzxH5lHUvYxZRYac&e= > 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
>
>
> 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 <https://urldefense.proofpoint.com/v2/url?u=http-3A__baret-5Fat-5Fapc.univ-2Dparis7.fr_&d=DwIFaQ&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=VEsjWcjHR_RwJLerQJyK4g&m=ZVvO6vStNXpBSohLqAcCrjzQnMkPhQgVHacnEQsVmmA&s=9zglRe7QEN9lgIa4wBLDc97lvliBzxH5lHUvYxZRYac&e= > 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
>
>
> 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 <https://urldefense.proofpoint.com/v2/url?u=http-3A__baret-5Fat-5Fapc.univ-2Dparis7.fr_&d=DwIFaQ&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=VEsjWcjHR_RwJLerQJyK4g&m=ZVvO6vStNXpBSohLqAcCrjzQnMkPhQgVHacnEQsVmmA&s=9zglRe7QEN9lgIa4wBLDc97lvliBzxH5lHUvYxZRYac&e= > 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?idc_info
Received on Sat Sep 22 2018 - 13:14:47 CEST

This archive was generated by hypermail 2.3.0 : Sat Sep 22 2018 - 13:14:57 CEST