|Message: Re: How to get SOBP output with Hadrontherapy example||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
On Sun, 04 Apr 2010 07:25:16 GMT, Jongseok wrote:
> Hello! I'm trying to get SOBP output. > > I changed the angle of the modulator whenever I simulated > proton_therapy.mac > > Like this > > /event/printEventNumber 1000 > > /modulator/angle 1 deg > > /analysis/setAnalysisFile deg1.root > > /run/beamOn 1000 /modulator/angle 1 deg > > /analysis/setAnalysisFile deg2.root > > /run/beamOn 1000 /modulator/angle 1 deg > > /analysis/setAnalysisFile deg3.root > > /run/beamOn 1000 /modulator/angle 1 deg > > /analysis/setAnalysisFile deg4.root > > /run/beamOn 1000 > > However, I couldn't get SOBP, and it's break peak points are almost > same!! > > I heard that I should change the angle of the modulator as a function of > time. > > Is it correct? > > If that is correct, what should I do for getting SOBP output. > > Please give me some advices. > > Thanks. >
Hadrontherapy to collect the total dose released by a beam into the phantom, doesn't matter if the beam is a mono-energetic or a modulated one, sum all the contribute at the end of a run. So you must not create a different file name for each angle but you have only change the various angles of the modulator. At the end, in each voxel/slice, the total energy due to the different contributes (i.e. the modulated dose distribution) will be available. I can also suggest you, in a first trial, to use the ASCII file generation of the Bragg peak that you can check in a more immediate way.
Moreover, the modulator wheel is already designed to get a very particular modulation and the only thing you have to do is to change the angles. The way as it is constructed contains the *time dependency* you mentioned as the different sectors have a different angular aperture so that beam intercept them for different periods of time.
I hope my explanation can help you.
Regards, Pabo Cirrone
|Inline Depth:||Outline Depth:||Add message:|