|Message: Re: Bragg peak of 330Mev/u of carbon using 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 Wed, 21 Apr 2010 08:38:56 GMT, Morad Hamad wrote:|
> Dear users/developers,
> when I ran default proton_therapy.mac I get a good Bragg peak which is
> agreeing with the experimental data includes inside the code. However,
> when I try to use ion_therapy.mac to compare it with 330 Mev/u of
> carbon12 from GSI, which has an experimentally a peak at about 20 cm in
> water I get a peak at about 6 cm, also I get the same depth when using
> energy of 135 Mev/u.
> I only changed the energy inside ion_therapy.mac to be 330 Mev/u once
> and 135 Mev/u other.
> I still don't know what is the physical meaning of sigmaenergy variable.
> Is the error because I don't use GSI synchrotron geometry? Or that I
> have to change other things inside either the code or the macro files?
> Please see the attached Bragg peak for 330 Mev/u of carbon-12.
> I will be appreciated for your assistance to solve the problem.
it simply seems to me that you are not changing correctly the energy of the carbon ions.
Firstly check the following points:
1.- The energy of your ions must be 330*12 = 3960 MeV as Geant4 requires the total kinetic energy and not the energy per nucleon
2.- Try to change the energy of the beam directly in the HadrontherapyPrimaryGeneratorAction.cc class file (where now the energy for a proton beam is defined as default), so to be sure that you are setting the correct energy value. I have in fact, the suspect that some default is changing the energy you set in the macro file.
Can you please report what version (both of Geant4 and Hadrontherapy) are you using?
Looking into the plot you show I must also recommend you to change the slice thickness of the phantom as, in the case of a carbon beam, the Bragg peak is very narrow and you need very thin slices to achieve a correct spatial resolution. I can suggest you to try to use a 50 micron slice with a 1 micron step max and 5 micron of production cut inside the detector.
Hoping this can help you,
|Inline Depth:||Outline Depth:||Add message:|