RE:[fluka-discuss]: Double differential source in FLUKA

From: Trinh Ngoc-Duy <Ngoc-Duy.Trinh_at_ganil.fr>
Date: Fri, 29 Jan 2016 08:45:49 +0000

Dear Paola Sala.

Thank you so much for your answer.

I will simulate my source with your advices.

For the party of beam direction in source.f, i have just found that i committed some problems of normalisation (cos(a)^2+cos(b)^2+cos(c)^2=1). It works ok now.

Sincerely.

ND Trinh, GANIL.

________________________________________
De : Paola Sala [paola.sala_at_mi.infn.it]
Envoyé : vendredi 29 janvier 2016 08:51
À : Trinh Ngoc-Duy
Cc : fluka-discuss_at_fluka.org
Objet : Re: [fluka-discuss]: Double differential source in FLUKA

Dear Trinh Ngoc-Duy,
In order to simulate the doubly-differential distribution, you should
- read in the angular distribution, exactly like you are doing for the
energy distribution, and build the cumulative distribution. Pay attention:
if you arngular distribution is given in unit of particles/steradian, the
value in each angular bin must be multiplied by the correspoonding solid
angle.
- sample from this cumulative as you do for the energy
-
-For the energy sampling: transform your ENECUM(I) and ENEPRO(I) vectrors
in two-dimensional matirices ENECUM(I,J) and ENEPRO(I,J), where J
corresponds to the angular bin.
For each angular bin J read the ENECUM(I=1,nene,J) values and calculate
the cumulative.
-Select your angle,, keep track of the angular interval you choose, then
select your energy from the J-th energy distribution

About the problem you have with selectring direction, I still do not
understand, could you please send the inpput and exactly what do you look
at to say that it does not work, please?
Regards
Paola



> Dear FLUKA experts,
>
> I'm using source.f routine to simulate some simple source like a square
> source with an energy distribution read from an external file.
>
> Now i want to extend it to simulate a double differential source (angular,
> energy distribution). But when i try for example to change the beam
> direction by changing the UBEAM, VBEAM and WBEAM value, it does not work.
> Even when i just put WBEAM=1 for a positive Z axis beam direction.
>
> TXFLK (NPFLKA) = UBEAM
> TYFLK (NPFLKA) = VBEAM
> TZFLK (NPFLKA) = WBEAM
>
> Otherwise, would you please show me how can i sample energy distribution
> for the corresponding angle (e.g the energy distribution in 0deg is
> different from the one of the 90deg.) ?
>
> I put my source.f file in input so you can investigate my problem.
>
> Thank you very much for your answer.
>
> Sincerely.
>
> ND Trinh, GANIL, France.
>
> ________________________________
> Préservons notre environnement, n’imprimez ce mail que si nécessaire.
> Preserve our environment, print this email only if necessary.
>


Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148

________________________________
Préservons notre environnement, n’imprimez ce mail que si nécessaire.
Preserve our environment, print this email only if necessary.

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id¬c_info
Received on Fri Jan 29 2016 - 11:00:40 CET

This archive was generated by hypermail 2.3.0 : Fri Jan 29 2016 - 11:00:44 CET