Re: Re: RE: [Gate-users] [ no subject ]
Jasmine Schirmer
jasmine.schirmer at tum.de
Wed Dec 19 08:18:39 CET 2007
Dear Simon,
For clarification, we are using the following versions:
GATE 3.0.0
geant 4.8.0.p01
clhep 2.0.2.2
and
GATE 3.1.2
geant 4.9.0.p01
clhep 1.9.3.1
We look forward to hearing from you about your simulations using the different GATE versions.
Jasmine and Irene
Original Message:
Hi Irene and Jasmine,
<br />
<br />Could you tell me please how long is your simulation 'single_layer.mac' ?
<br />Because if it is not too much, I will run it will all different versions of
<br />Gate between 3.0.0 and 3.1.2 to see what are the differences and locate
<br />between which ones the problem occured.
<br />
<br />I would like to know too what version of Geant4 and clhep do you use with
<br />these different Gate versions. Because it has been demonstrated that with a
<br />wrong version of Geant4 and clhep, you could get wrong results. Do you
<br />really use the correct ones ?
<br />
<br />Hope we fix this problem as soon as possible,
<br />Simon
<br />
<br />2007/12/18, Jasmine Schirmer :
<br />>
<br />> Dear Min-Jae and GATE developers,
<br />>
<br />> Thank-you for looking into this problem. We are attaching an example of a
<br />> simple GATE simulation (single_layer.mac and the supplementary macros:
<br />> no_visu.mac, physics.mac, single_layer_allmod.mac, CylFOV_Ph.mac and
<br />> G_CapillariesS_5m_2D.mac ). You can run this simulation by using the command
<br />> ./control/execute/single_layer.mac. We have run this on two versions of
<br />> GATE. You will see that if you run this simulation using the newest version
<br />> 3.1.2, you will obtain for the singles and coincidences, energy spectra
<br />> which do not fall to zero. The output is as follows for a 400keV low energy
<br />> threshold:
<br />>
<br />> There are: 868838 Singles
<br />> => 35514 coincidences
<br />> => 4395 intercrystal coincidences
<br />> => 9888 randoms
<br />> => 2264 intercrystal randoms
<br />>
<br />> where intercrystal coincidences are true coincidences detected in the same
<br />> rsector. On the other hand, running the SAME macro under an older version of
<br />> GATE (we used 3.0.0), the energy spectra falls to zero as expected with
<br />> the following output:
<br />>
<br />> There are: 867717 Singles
<br />> => 50812 coincidences
<br />> => 0 intercrystal coincidences
<br />> => 9438 randoms
<br />> => 2199 intercrystal randoms
<br />>
<br />> This bug in the newest version of GATE means that we cannot use the output
<br />> data reliably and our reconstructed images are incorrect.
<br />>
<br />> We have noticed that this problem has also been detected by Georgi
<br />> Gerganov. As well, we have seen the same problem using an ECAT simulation
<br />> (0_ScriptSystem.mac and 0_ScriptSystem.mac). The resulting plots from this
<br />> simulation have been posted with the November 29th email to the GATE mailing
<br />> list.
<br />>
<br />> Please let us know if you have any further questions.
<br />>
<br />> Thanks in advance,
<br />> Irene and Jasmine
<br />>
<br />> Original Message:
<br />>
<br />> Dear Irene,
<br />> <br />
<br />> <br />Let me check what you are saying...
<br />> <br />
<br />> <br />There is no problem with benchmarkPET...
<br />> <br />[In both the old and new versions of GATE, the benchmark PET was
<br />> <br />successfully completed and no difference in the energy spectra was
<br />> <br />observed.]
<br />> <br />
<br />> <br />In your simulation code, uphold value is 950 keV...
<br />> <br />[However, the upper energy threshold of the benchmark is set to
<br />> <br />650keV while our simulations are performed at higher upper energy
<br />> <br />thresholds.]
<br />> <br />[/gate/digitizer/Singles/upholder/setUphold 950. keV]
<br />> <br />
<br />> <br />In fact, I can't get that the benchmark is same but your code is
<br />> <br />different...
<br />> <br />
<br />> <br />Can I take a look on your code and simulate your code with my gate?
<br />> <br />
<br />> <br />MJ
<br />> <br />
<br />> <br />-----Original Message-----
<br />> <br />From: gate-users-bounces at lists.healthgrid.org
<br />> <br />[mailto:gate-users-bounces at lists.healthgrid.org] On Behalf Of
<br />> Jasmine
<br />> <br />Schirmer
<br />> <br />Sent: Monday, December 17, 2007 11:24 PM
<br />> <br />To: gate-users at lists.healthgrid.org
<br />> <br />Subject: [Gate-users] [ no subject ]
<br />> <br />
<br />> <br />Dear GATE users and developers,
<br />> <br />
<br />> <br />A few weeks ago, we sent an email to the users list concerning a
<br />> problem we
<br />> <br />are having with the newest version of GATE. The email was dated the
<br />> 29th of
<br />> <br />November and is titled 'differences in energy spectra from different
<br />> GATE
<br />> <br />releases'.
<br />> <br />
<br />> <br />At the moment, we can no longer use reliably the newest version of
<br />> GATE
<br />> <br />until this issue is resolved. To recall your attention, we showed a
<br />> <br />difference in the energy spectra produced by two identical
<br />> simulations using
<br />> <br />two versions of GATE (version 3.0.0 and the newest 3.1.2). The
<br />> newest
<br />> <br />version gives a different and unexpected result for the energy
<br />> spectra as
<br />> <br />you can see from the figures attached in the previous email. (The
<br />> energy
<br />> <br />spectra does not fall to zero!)
<br />> <br />
<br />> <br />We kindly ask the developers and users to address this issue. In
<br />> addition,
<br />> <br />if anyone would like to perform a similiar test on their version of
<br />> GATE, we
<br />> <br />would be happy to provide the macro we used.
<br />> <br />
<br />> <br />Thanks in advance,
<br />> <br />Irene Torres and Jasmine Schirmer
<br />> <br />
<br />> <br />_______________________________________________
<br />> <br />Gate-users mailing list
<br />> <br />Gate-users at lists.healthgrid.org
<br />> <br />http://lists.healthgrid.org/mailman/listinfo/gate-users
<br />> <br />
<br />>
<br />> --
<br />>
<br />> _______________________________________________
<br />> Gate-users mailing list
<br />> Gate-users at lists.healthgrid.org
<br />> http://lists.healthgrid.org/mailman/listinfo/gate-users
<br />>
<br />>
<br />>
<br />
--
More information about the Gate-users
mailing list