Hi Mikhail,
Your formula for uniform distribution in polar coordinates is not wright,
as it favors probability closer to the axis.
Just replace RADIUS by SQRT(RADIUS) both for XFLK and YFLK and you will
obtain what you seek.
Mario
On 10/17/14 3:12 AM, ""Mikhail Polkovnikov <Михаил Полковников>""
<pmk_at_ihep.ru> wrote:
>Dear FLUKA users,
>
>
>For 70 GeV proton beam with annular form and Rmax = 4.0cm
>
>*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....
>+....
>BEAM -70.0 0.0 0.0 4.0 0.0 -1.0PROTON
>BEAMPOS 0.0 0.0 -79.0 0.0 0.0
>
> i need to add some spectrum in source.f, but in the user routine the
>annular beam always changes its profile even if add uniform distribution
>like
>
> CALL SFECFE( SFE, CFE)
> RADIUS = XSPOT * FLRNDM(XDUMMY)
> XFLK (NPFLKA) = RADIUS * SFE
> YFLK (NPFLKA) = RADIUS * CFE
> ZFLK (NPFLKA) = ZBEAM
>
>Is there any way to peek uniform beam profile in user routine?
>
>Images of beam profiles along x-axis in attachment
>
>
>Best regards,
>Mikhail
>
Received on Fri Oct 17 2014 - 21:52:48 CEST