Re: [fluka-discuss]: Re: regarding x-ray beam simulation of dental machines

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Fri, 28 Jul 2017 13:45:31 +0200

Dear Alex,

> Do I understand correctly that I get results in pSv per photon particle?

yes.

> So the spectrum should be measured or generated via
> special software generators and should already contain the desirable
> intensity?

Not at all.
But, as you say, you still need to know your photon intensity in order to
normalize the above results.

Kind regards

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Thu, 27 Jul 2017, Ševčik Aleksandras wrote:

>
> Dear experts,
>
>
> Still want to ask now you about the scoring in such case when custom
> source.f is used. As in this case I don't need any decays/nucleis/etc, I can
> just use usrbin with associated auxscore for calculation of effective dose.
> Do I understand correctly that I get results in pSv per photon particle? So
> to know the actual dose I need to convert the x-ray tube's voltage kV,
> current mA and time t to nr of photons? If yes, is there any advisable
> convenient method to do so?
>
>
> Or I am on false way and the whole information should be delivered with
> spectrum.dat table? So the spectrum should be measured or generated via
> special software generators and should already contain the desirable
> intensity?
>
>
> Thanks, Alex
>
> ____________________________________________________________________________
> From: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> on
> behalf of Ševčik Aleksandras <aleksandras.sevcik_at_ktu.edu>
> Sent: Wednesday, July 26, 2017 10:30:04 PM
> To: Eleftherios Skordis; Mikhail Polkovnikov <Михаил Полковников>;
> fluka-discuss_at_fluka.org
> Subject: Re: [fluka-discuss]: Re: regarding x-ray beam simulation of dental
> machines  
>
> Dear Mikhail, Lefteris,
>
>
> Thanks again, as always, for your support, it seems it works nicely,
> attached pictures. 
>
>
> Alex
>
> ____________________________________________________________________________
> From: Eleftherios Skordis <eleftherios.skordis_at_cern.ch>
> Sent: Wednesday, July 26, 2017 3:20:10 PM
> To: Mikhail Polkovnikov <Михаил Полковников>; Ševčik Aleksandras;
> fluka-discuss_at_fluka.org
> Subject: RE: [fluka-discuss]: Re: regarding x-ray beam simulation of dental
> machines  
> Indeed Mikhail is correct.
> I misinterpreted your question.
>
> sorry for the confusion
>
> Regards
>
> Lefteris
>
>
> -----------------------
>
> Eleftherios Skordis
> Dep. EN/STI, CERN
> CH-1211 GENEVA 23
> SWITZERLAND
>
> OFFICE: +41-22-7679541
>
> ____________________________________________________________________________
> From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on
> behalf of Mikhail Polkovnikov <Михаил Полковников> [pmk_at_ihep.ru]
> Sent: 26 July 2017 13:00
> To: Ševčik Aleksandras; fluka-discuss_at_fluka.org
> Subject: Re: [fluka-discuss]: Re: regarding x-ray beam simulation of dental
> machines
>
> On 26.07.2017 01:48, Ševčik Aleksandras wrote:
>
> Dear Lefteris, Mikhail,
>
>
> Thank you very much, it works perfectly. May I ask one more
> question:
>
> 1) the shape of the beam is annular. If I want to make it
> rectangular, is it enough to change this in BEAM card
> characteristics, or should I modify the source code too?
>
>
> I have a voxelised head phantom and want to score effective dose
> from the simulated spectra in the special regions.  I will
> proceed further with help from Activation/AdvancedActivation
> slides, and will try to compare the results from simulation with
> mathematical calculations used by our medical physicists. In
> case of strong deviation, I may ask for your advice once again,
>
>
> Alex
>
> ____________________________________________________________________________
> From: Eleftherios Skordis <eleftherios.skordis_at_cern.ch>
> Sent: Friday, July 21, 2017 12:15:48 PM
> To: Ševčik Aleksandras; Mikhail Polkovnikov <Михаил Полковников>;
> fluka-discuss_at_fluka.org
> Subject: RE: [fluka-discuss]: Re: regarding x-ray beam simulation of
> dental machines  
> Hi Alex,
> It could be that there is something missing from the flupix version.
> However, until this is properly fixed you can try running the
> following command in the flupix terminal.
>
> sudo yum install glibc-devel.i686 libgcc.i686 libstdc++-devel.i686 ncurses-d
> evel.i686
> In my case it solved the problem.
> Hope this helps.
>
> Kind regards
>
> Lefteris
>
>
> -----------------------
>
> Eleftherios Skordis
> Dep. EN/STI, CERN
> CH-1211 GENEVA 23
> SWITZERLAND
>
> OFFICE: +41-22-7679541
>
> ____________________________________________________________________________
> From: Ševčik Aleksandras [aleksandras.sevcik_at_ktu.edu]
> Sent: 21 July 2017 09:49
> To: Mikhail Polkovnikov <Михаил Полковников>; fluka-discuss_at_fluka.org;
> Eleftherios Skordis
> Subject: Re: [fluka-discuss]: Re: regarding x-ray beam simulation of
> dental machines
>
> Thanks Mikhail, now it is clear that I am doing smth wrong when
> compiling/running the custom executable, as Eleftherios guessed.
>
>
> When I try to compile it, I get the error, see attached. What it could
> be, looks like cannot find smth ? I am running 64bit flupix, and the
> folder is /home/flupix/New
>
>
> Alex
>
> ____________________________________________________________________________
> From: Mikhail Polkovnikov <Михаил Полковников> <pmk_at_ihep.ru>
> Sent: Friday, July 21, 2017 7:04:32 AM
> To: Ševčik Aleksandras; fluka-discuss_at_fluka.org
> Subject: Re: [fluka-discuss]: Re: regarding x-ray beam simulation of
> dental machines  
> On 20.07.2017 20:13, Ševčik Aleksandras wrote:
>
> Thank you both for your attention, and special thanks
> to Mikhail, for your created file which I am using.
>
>
> I added the proposed line, however, my little experiment
> still shows strong deviation between exemplary ref
> spectrum and scored spectrum, see attached files.
>
>
> I also attach all the input files. Where the problem could
> be?
>
>
> A.
>
> ____________________________________________________________________________
> From: Eleftherios Skordis <eleftherios.skordis_at_cern.ch>
> Sent: Thursday, July 20, 2017 10:43:48 AM
> To: Ševčik Aleksandras; fluka-discuss_at_fluka.org
> Subject: RE: regarding x-ray beam simulation of dental machines
>  
> Hi Alex,
> Your points 1 and 3 are correct. My guess is you are not using
> the customised executable when you are running and therefore not
> generating a multi-energetic beam.
>
> Could you please send us your input and output files so I can
> help you further.
>
> Kind regards
>
> Lefteris
>
>
> -----------------------
>
> Eleftherios Skordis
> Dep. EN/STI, CERN
> CH-1211 GENEVA 23
> SWITZERLAND
>
> OFFICE: +41-22-7679541
>
> ____________________________________________________________________________
> From: owner-fluka-discuss_at_mi.infn.it
> [owner-fluka-discuss_at_mi.infn.it] on behalf of Ševčik Aleksandras
> [aleksandras.sevcik_at_ktu.edu]
> Sent: 19 July 2017 22:26
> To: fluka-discuss_at_fluka.org
> Subject: [fluka-discuss]: Re: regarding x-ray beam simulation of
> dental machines
>
> and hello again, sorry for multiple post, but I thought that
> it's always better to try a practical example, so here it is:
> using the attached files and artificially made spectrum.dat for
> example purposes, I simulated 5x1E6 histories and scored the
> spectra which does not look similar, maybe anybody could comment
> what I misunderstand?
>
>
> Thank you very much,
>
>
> Alex
>
> ____________________________________________________________________________
> From: owner-fluka-discuss_at_mi.infn.it
> <owner-fluka-discuss_at_mi.infn.it> on behalf of Ševčik Aleksandras
> <aleksandras.sevcik_at_ktu.edu>
> Sent: Wednesday, July 19, 2017 9:02:49 PM
> To: fluka-discuss_at_fluka.org
> Subject: [fluka-discuss]: Re: regarding x-ray beam simulation of
> dental machines  
>
> Addition: I was able to find the answer to 2nd question,
> it contains pairs of values of energy in keV and flux in photons
> per second, so it can be ignored.
>
>
> A.
>
> ____________________________________________________________________________
> From: Ševčik Aleksandras
> Sent: Wednesday, July 19, 2017 7:59:02 PM
> To: Vasilis Vlachoudis; fluka-discuss_at_fluka.org
> Subject: Re: regarding x-ray beam simulation of dental machines
>  
>
> Thank you very much. I looked all through it, however, it is too
> complicated for me to pursue it in right way from very
> beginning without detail step-by-step guidance and clear
> examples. Therefore  I am happy that fluka archive helped me  to
> find the similar case with example I could probably use for my
> needs. Yet few clarifications are needed:
>
>
> Attached modified source.f and spectrum source file. Also
> attached dental x-ray spectra for reference.
>
>
> 1) Do I correctly understand that this source.f file generates
> polyenergetic photon spectrum, .i.e. such as shown in the
> attached picture?
>
> 2) Could you please help me to understand what exactly is given
> in spectrum.dat? The first number is just a number of bin, and
> second - ?
>
> 3) Do I understand correctly that in addition I must use BEAM
> with energy set equal to x-ray used kVp, i.e. maximum energy ?
>
>
> Thank you very much,
>
> Alex
>
> ____________________________________________________________________________
> From: Vasilis Vlachoudis <Vasilis.Vlachoudis_at_cern.ch>
> Sent: Monday, July 3, 2017 10:08:55 AM
> To: Ševčik Aleksandras; fluka-discuss_at_fluka.org
> Subject: RE: regarding x-ray beam simulation of dental machines
>  
> Hi Aleksandras,
>
> my guess is that the best way will be to create a custom
> source.f routine
> that you sample directly from the spectrum of your x-ray
> equipment.
>
> Look the advanced course on how to make a custom source
>
> https://indico.cern.ch/event/489973/contributions/2148893/attachments/12718
> 19/1884943/11_AdvancedSources2016.pdf
>
> Cheers
> Vasilis
>
>
> ____________________________________________________________________________
> From: owner-fluka-discuss_at_mi.infn.it
> [owner-fluka-discuss_at_mi.infn.it] on behalf of Ševčik Aleksandras
> [aleksandras.sevcik_at_ktu.edu]
> Sent: Friday, June 30, 2017 16:19
> To: fluka-discuss_at_fluka.org
> Subject: [fluka-discuss]: regarding x-ray beam simulation of
> dental machines
>
>
>  
>
> Dear experts,
>
>
> I am interested in the calculation of effective doses during the
> dental -xray  procedures. Is it possible to generate a beam only
> of  the specific piece of x-ray equipment when manufacturer's
> parameters are known, or should one model all the internal
> components of the machine i.e. the whole x-ray generation
> process in order to have valid results?
>
>
> Sorry if this question too obvious. In addition, I would
> appreciate any links/comments/references that would help me in
> pursue of this matter,
>
>
> Regards
>
> Aleksandras
>
> Dear Alex,
>
>
> It looks like you didn't compile a source.f routine and didn't run the
> modified FLUKA executive file with the modified source routine (see
> FM.pdf, Chapter 13, User routines). That is why your simulation uses
> only monoenergetic value from the BEAM card!
>
>
> There is no need to use OPEN card in an input file, when you open file
> for reading in source.f (OAUXFI or OPEN calls).
>
>
> Just use either
>
> 1) OAUXFI call in source.f and disable OPEN card
>
> or
>
> 2) OPEN card and do not open file for reading in source.f once again.
>
>
> I prefer first approach nowadays.
>
>
> As a result your scored spectrum should be OK.
>
>
>
> Best regards,
>
> Mikhail
>
> Dear Alex,
>
>
> Since you have been sampling an energy distribution in your custom source,
> you must sample a spatial distribution as well!
>
>
> For rectangular beam shape:
>
>
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+
> ....
> BEAM           7E-05                           3.0       5.0      
> 1.0PHOTON   
>
>
> source file:
>
>
>       ZFLK   (NPFLKA) = ZBEAM
>       XFLK   (NPFLKA) = XBEAM + (TWOTWO * FLRNDM(XDUMMY) - ONEONE)
>      &  * XSPOT
>       YFLK   (NPFLKA) = YBEAM + (TWOTWO * FLRNDM(YDUMMY) - ONEONE)
>      &  * YSPOT
>
>
>
> Best regards,
>
>
> Mikhail
>
>
>

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Jul 28 2017 - 14:57:06 CEST

This archive was generated by hypermail 2.3.0 : Fri Jul 28 2017 - 14:57:26 CEST