Continuing the discussion from Souce for Simulation of beam in a Cyclotron:
Dear FLUKA experts,
as discussed in the above linked post, I tried to implement a circular source using the source_newgen.f subroutine.
The cyclotron center is located at x0, y0, z0 as stated in the coordinate section.
To generate the particles uniformly around the circumference I used:
coordinate_x = x0 + cyclo_radius * COS(phi)
coordinate_y = y0 + cyclo_radius * SIN(phi)
For the direction of the particles I entered:
direction_cosx = SIN(phi)
direction_cosy = COS(phi)
I also include a SOURCE card in the source code and it is compiling and running.
The output of the simulation when scoring BEAMPART in the outer and inner part of the cyclotron, using USRBIN-R-phi-Z binning, looks like this:
My question now is, why the primary particles are distributed in such a weird way instead of around the whole circumference.
Thank you in advance
source_newgen.f (19.0 KB)