Hi again,
Just to confirm that ROTDEFIni works perfectly and there is no need for the angular adjustment of beam itself…sometimes obvious ideas came a lit bit later, sorry and thank you for the help!
Aleks
From: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> On Behalf Of Ševcik Aleksandras
Sent: Tuesday, August 21, 2018 18:29
To: Mikhail Polkovnikov <Михаил Полковников> <pmk_at_ihep.ru>; fluka-discuss_at_fluka.org
Subject: RE: [fluka-discuss]: changing beam direction in source.f
Hi Mikhail,
Yes, exactly. Now I am studying the possibility to use ROTDEFIni card to rotate my voxel geometry – this also would be an option. I hope that using this card is straight-forward and would not require more fortran programming. I will test it. I will appreciate any comments/insights. In general, my situation is like this: I have a complex voxel geometry and want to be able to irradiate it with custom source spectra from any angle.
Rgds
Aleksandras
From: Mikhail Polkovnikov <Михаил Полковников> <pmk_at_ihep.ru<mailto:pmk_at_ihep.ru>>
Sent: Tuesday, August 21, 2018 12:50
To: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>; Ševčik Aleksandras <aleksandras.sevcik_at_ktu.edu<mailto:aleksandras.sevcik_at_ktu.edu>>
Subject: Re: [fluka-discuss]: changing beam direction in source.f
Hi Aleksandras,
Do you need something like a rotating source (gantry) around the Z axis with an isocenter?
Best regards,
Mikhail
On 20.08.2018 13:36, Ševčik Aleksandras wrote:
Dear Francesc, Mikhail,
Thank you very much, it works now very well. Due to my very limited knowledge in fortran I am very dependent of such help. So now I have x-axis aligned rectangular shape customized beam from source.f. One final thing now I need is to be able to move this beam around Z axis. How exactly I can modify source.f to enable this, i.e. one of these variables should be multiplied by cos function?
Rgds
A.
-----Original Message-----
From: Francesc Salvat-Pujol <francesc.salvat.pujol_at_cern.ch><mailto:francesc.salvat.pujol_at_cern.ch>
Sent: Monday, August 20, 2018 11:33
To: Ševčik Aleksandras <aleksandras.sevcik_at_ktu.edu><mailto:aleksandras.sevcik_at_ktu.edu>; fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: changing beam direction in source.f
Should of course be
YSPT = 2.5D+00
ZSPT = 2.5D+00
Cesc
On Mon, Aug 20 2018, at 09:53 +0200, Francesc Salvat-Pujol wrote:
Dear Alex,
If you find it more comfortable, you can consider adapting something
along the following lines in source.f:
XSPT = 2.5D+00
YSPT = 2.5D+00
XFLK(NPFLKA) = XBEAM
YFLK(NPFLKA) = YBEAM + (TWOTWO * FLRNDM(YDUMMY) - ONEONE)
& * YSPT
ZFLK(NPFLKA) = ZBEAM + (TWOTWO * FLRNDM(ZDUMMY) - ONEONE)
& * ZSPT
You can set the beam direction either with the BEAMPOS card or with the
TXFLK,TYFLK,TZFLK variables, i.e., something a la:
TXFLK(NPFLKA) = ONEONE
TYFLK(NPFLKA) = ZERZER
TZFLK(NPFLKA) = ZERZER
To be adapted (and tested/debugged...) as needed.
With kind regards,
Cesc
On Sun, Aug 19 2018, at 20:51 +0000, Ševčik Aleksandras wrote:
Dear experts
I'm still looking for the solution for this problem as it looks that
it's not so straightforward to use BEAMAXES option with source.f, at
least it doesn't work for me. I think maybe it is possible to define
the particular size of rectangular beam shape inside the source.f,
let's say 5x5 cm, so I could avoid using BEAM WHAT(3-5) parameters and
could use only BEAMPOS card to align it with x axis ? How could I do
this?
Problem is that I have a complex voxel geometry (head) and I want to
irradiate the particular spot
Alex
-----Original Message-----
From: owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it> <owner-fluka-discuss_at_mi.infn.it><mailto:owner-fluka-discuss_at_mi.infn.it>
On Behalf Of Ševcik Aleksandras
Sent: Tuesday, August 14, 2018 12:30
To: Francesc Salvat-Pujol <francesc.salvat.pujol_at_cern.ch><mailto:francesc.salvat.pujol_at_cern.ch>
Cc: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: RE: [fluka-discuss]: changing beam direction in source
Dear Francesc
Kindly thanks for the insights. I totally forgot BEAMAXES option, have
never used before. So, I guess it can be applied in this situation
even without modifying the source.f ,
Regards
Alex
-----Original Message-----
From: Francesc Salvat-Pujol <francesc.salvat.pujol_at_cern.ch><mailto:francesc.salvat.pujol_at_cern.ch>
Sent: Tuesday, August 14, 2018 10:52
To: Ševčik Aleksandras <aleksandras.sevcik_at_ktu.edu><mailto:aleksandras.sevcik_at_ktu.edu>
Cc: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: changing beam direction in source.f
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
--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
b��jy���‑�{.n�+��b�w�ka�ɚ�X��X���a��l��0��n��+���F������w
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Aug 22 2018 - 09:44:30 CEST