[fluka-discuss]: Re: Question about 150MeV electron in Aluminum CSDA

From: Santana, Mario <msantana_at_slac.stanford.edu>
Date: Thu, 9 Nov 2023 02:08:12 +0000

Some comments about your input file:

  * Your aluminum cylinder is 5 cm in radius, but you are scoring 10 cm in radius (?)
  * You did not set transport thresholds for EM
  * Fudgem=1 seems OK, I don’t think you need to go below 100 keV.



From: neu lzf <lzfneu_at_live.com>
Date: Wednesday, November 8, 2023 at 4:47 AM
To: Santana, Mario <msantana_at_slac.stanford.edu>
Cc: fluka-discuss_at_fluka.org <fluka-discuss_at_fluka.org>, neu lzf <lzfneu_at_live.com>
Subject: Re: Question about 150MeV electron in Aluminum CSDA
Thanks for your reply.

I want to consult you that if I make the right settings to use EMFCUT card with sdum=PROD-CUT to suppress photon production.
My settings are: the electron production threshold set 1e-4 GeV and photon threshold set 0.151 GeV, and Fudgem parameter set 1.0.

For Fudgem parameter the manual says for production and transport cutoffs larger than 100 keV Fudgem parameter must be set = 1.0, while in the keV region it must be set = 0.0. Therefore, I'm not sure if I am right to set Fudgem parameter equal to 1.0, because electron production threshold are equal to 100 keV while photon production threshold are larger than 100 keV.

I have attached my input and result files in the attachment. It seems that the peak starting at 11.6cm then sharp drop extending to ~15cm. Therefore, I am not sure if my settings are right. Could you help to check the input file and thanks in advance.

________________________________
From: Santana, Mario <msantana_at_slac.stanford.edu>
Sent: 08 November 2023 3:31
To: neu lzf <lzfneu_at_live.com>; fluka-discuss_at_fluka.org <fluka-discuss_at_fluka.org>
Cc: neu lzf <lzfneu_at_live.com>
Subject: Re: Question about 150MeV electron in Aluminum CSDA


Hi,



Since CSDA is an artificial average construct, what you would want to do is switch radiative fluctuations to the range by setting a high photon production threshold (e.g., 151 MeV). Other than that, you can also switch photonuclear reactions (not sure that the impact will be huge, though). The other options that you discuss relate to analogous calculations for heavier particles such as muons and ions, thus no need to worry about those.

If you then plot the 1-D fluence plot for BEAMPART as you described, then you will see a flat function with a sharp drop starting at around 14 cm depth, which will extend beyond 14.71cm (remember CSDA is an average number not accounting for various fluctuations), but with very small probability (by the way, it might be interesting to play with the EMFCUT thresholds to see if/how they impact the tail).



-M



From: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> on behalf of neu lzf <lzfneu_at_live.com>
Date: Monday, November 6, 2023 at 10:13 PM
To: fluka-discuss_at_fluka.org <fluka-discuss_at_fluka.org>
Cc: neu lzf <lzfneu_at_live.com>
Subject: [fluka-discuss]: Question about 150MeV electron in Aluminum CSDA

Dear All,



As is known that 150MeV electron in Aluminium CSDA is 14.71cm. I want to calculate it by FLUKA.



To do this, I am turning off hadron elastic and inelastic reactions, turning off the MCS of primary and secondary charged hadrons, turning off FLUC and turning off production of secondary delta rays, according to the post:

http://www.fluka.org/web_archive/earchive/new-fluka-discuss/2808.html



One of my questions is where to find the correct result.

I try to use USRBIN energy deposition coupled with BEAMPART to calculate 150MeV electron in Aluminium CSDA range. However, the result is only about 7.5cm and is not 14.71cm as I expected. Is it the plot range setting in FLAIR problem and could you please help to check my input file. I have attached my input file and flair file in the attachment.



Thanks in advance.

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Nov 09 2023 - 05:02:59 CET

This archive was generated by hypermail 2.3.0 : Thu Nov 09 2023 - 05:03:00 CET