Dear Alex,
I could be wrong, but you might be in a situation where Note 4 of the
{BEAM} card section of the manual applies:
4) All options governed by WHAT(3,4,5) are meaningful only if the beam
direction is along the positive z axis, unless a command BEAMAXES is
issued to establish a beam reference frame different from the
geometry frame (see command BEAMAXES). If the beam is not in the
positive z direction and no BEAMAXES command has been given,
WHAT(3)-WHAT(5) must be set = 0.0 (unpredictable effects would arise
otherwise).
I gather your beam is along x and you wanted to pass the (spatial) beam
width via WHAT(4) and WHAT(5) of the BEAM card, foreseen for beams along
the positive z axis instead, as per the note above in this particular
scenario.
You can try to set the beam direction with the TXFLK,TYFLK,TZFLK
variables in source.f (mind the normalization to ONEONE).
Hope this is not too far off!
Cheers,
Cesc
On Mon, Aug 13 2018, at 20:41 +0000, Ševčik Aleksandras wrote:
>
>Dear experts,
>
>I have this source.f (attached) routine which samples and generates rectangular beam. When I use BEAMPOS for defining various angles and starting points, I can see that it distorts the rectangular shape then. For example if I set BEAMPOS x=-110 cosx=1 to make it parallel to x axis (see attached input file) I can see that starting point of the beam does not look like rectangular one (attached jpg). If I use BEAMPOS without angular adjustment, i.e. without setting cos, and only for initial coordinates of the beam point, it works well.
>Does this mean that angular position of the custom beam must be defined in source.f , not in BEAMPOS ? If yes, what exactly should be added in the source.f to enable the angular adjustment of the custom beam?
>Thank you,
>Alex
--
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?id=acc_info
Received on Tue Aug 14 2018 - 11:20:45 CEST