[gate-users] gate-users Digest, Vol 8, Issue 10 SPECT benchmark problems of Ywyang
karine.assie at univ-rouen.fr
karine.assie at univ-rouen.fr
Thu Nov 18 10:11:50 CET 2004
Hi Ywyang and all,
I have checked your GateSingles.dat files. Everything seems to be right.
I have downloaded the last version of GATE on the GATE web site. I have
recompiled the programm benchmark_spectra.c. My results are the following:
Number of counts for the total energy spectrum = 35901
Number of counts for the Unscattered spectrum =12778
Number of counts for the total scatter spectrum =23123
Number of counts for the crystal scatter spectrum =2345
Number of counts for the shie scatter spectrum = 2
Number of counts for the coll scatter spectrum = 114
Number of counts for the comp scatter spectrum = 439
Number of counts for the tabl scatter spectrum = 1008
Number of counts for the movs scatter spectrum = 4944
Number of counts for the Phan scatter spectrum = 14271
Number of counts for the order1 scatter spectrum = 11202
Number of counts for the order2 scatter spectrum = 6046
Number of counts for the order3 scatter spectrum = 2902
Number of counts for the order4 scatter spectrum = 1543
Number of counts for scatter orders greater than 4 spectrum = 1430
##########################################################################
Number of detected counts from 0.020 to 0.190 MeV : 35901
Unscattered events (percentage of total counts)= 35.59
Scatter in the crystal (percentage of total counts) = 6.53
Scatter in the shie volume (percentage of total scatter) = 0.01
Scatter in the coll volume (percentage of total scatter) = 0.32
Scatter in the comp volume (percentage of total scatter) = 1.22
Scatter in the tabl volume (percentage of total scatter) = 2.81
Scatter in the movs volume (percentage of total scatter) = 13.77
Scatter in the Phan volume (percentage of total scatter) = 39.75
Scatter order1 (percentage of total scatter) = 48.45
Scatter order2 (percentage of total scatter) = 26.15
Scatter order3 (percentage of total scatter) = 12.55
Scatter order4 (percentage of total scatter) = 6.67
Scatter orders greater than 4 (percentage of total scatter) = 6.18
My results don't agree with yours. Moreover, my results agree with the mean
results presented in the results_benchmark.doc files.
I think that your output file is OK and that your simulation is right. Perhaps
you could have done a mistake during the compiling of the *.c files.
Hope this help,
Assié Karine
> Send gate-users mailing list submissions to
> gate-users at lphe1pet1.epfl.ch
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lphe1pet1.epfl.ch/mailman/listinfo/gate-users
>
> Please use the GATE's site webmaster e-mail to reach the person
> managing the list.
>
>
> Today's Topics:
>
> 1. Cygwin (Nehmeh, Sadek/Medical Physics)
> 2. PETbench works fine but SPECTbench spectra goes wrong (ywyang)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 12 Nov 2004 15:21:20 -0500
> From: "Nehmeh, Sadek/Medical Physics" <nehmehs at mskcc.org>
> Subject: [gate-users] Cygwin
> To: "GATE feedback and helpline for Users"
> <gate-users at lphe1pet1.epfl.ch>
> Message-ID:
> <2867578BB7767E45B3C9E3CBA9C5A65F7762C3 at smskpexmbx3.mskcc.root.mskcc.org>
>
> Content-Type: text/plain; charset=us-ascii
>
>
> Dear cygwin users,
>
> I would like to mention that running Gate on cygwin may not be reliable.
> Eventhough running the GATE benchmarks produced the same results as
> those reported by the GATE collaboration, "sometimes" results do not
> match those from running on RedHat. Another major problem, I could not
> let GATE decay a radioactive particle in cygwin. Running the same mac
> files on Linux was just fine. Hope this helps.
>
> Regards,
>
> -----------------------------------------
> Sadek A. Nehmeh, Ph.D.
> Assistant Attending Physicist
> Medical Physics Department
> Nuclear Medicine Service
> Memorial Sloan-Kettering Cancer Center
> Tel: (212)-639-2175
> Email: nehmehs at mskcc.org
> ------------------------------------------
>
>
> =====================================================================
>
> Please note that this e-mail and any files transmitted with it may be
> privileged, confidential, and protected from disclosure under
> applicable law. If the reader of this message is not the intended
> recipient, or an employee or agent responsible for delivering this
> message to the intended recipient, you are hereby notified that any
> reading, dissemination, distribution, copying, or other use of this
> communication or any of its attachments is strictly prohibited. If
> you have received this communication in error, please notify the
> sender immediately by replying to this message and deleting this
> message, any attachments, and all copies and backups from your
> computer.
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 13 Nov 2004 09:46:08 +0800
> From: "ywyang" <ywyang98 at ms25.hinet.net>
> Subject: [gate-users] PETbench works fine but SPECTbench spectra goes
> wrong
> To: "ywyang" <ywyang98 at ms25.hinet.net>, <gate-users at lphe1pet1.epfl.ch>
> Message-ID: <000c01c4c922$8befb4f0$0c0419ac at pet17585>
> Content-Type: text/plain; format=flowed; charset="iso-8859-1";
> reply-type=response
>
> Dear All :
>
> I am struggling with benchmarks. I followed the instructions of
> manuals and got help from this group. I used the mac files downloaded.
>
> It seems everything is OK when the results of PET bench are comparable
> with manuals. But benchSPECT did not work well. The scatters in medium
> and scatter order was different from the results showed in the manuals.
> Is it allright or I have some troubles in installation ?
>
> Followed are the messages.
>
> [gate at localhost benchmarkSPECT]$ ./benchmark_spectra gateSingles.dat
> Enter the lower energy threshold (MeV)
>
> 0.02
>
>
> Lower energy threshold is 0.020000 MeV
>
>
> Enter the upper energy threshold (MeV)
>
> 0.19
>
> Upper energy threshold is 0.190000 MeV
>
> Enter the number of energy windows
>
> 100
>
> Energy range will be divided into 100 energy windows
>
> Enter the number of scattering volume (max is 10)
>
> 6
>
> 6 scattering volumes will be considered
>
> Enter the first four letters of the scattering volume number 1
>
> shie
>
> Enter the first four letters of the scattering volume number 2
>
> coll
>
> Enter the first four letters of the scattering volume number 3
>
> comp
>
> Enter the first four letters of the scattering volume number 4
>
> tabl
>
> Enter the first four letters of the scattering volume number 5
>
> movs
>
> Enter the first four letters of the scattering volume number 6
>
> Phan
>
>
>
> The scattering volumes that will be considered are the following:
>
> shie
>
> coll
>
> comp
>
> tabl
>
> movs
>
> Phan
>
>
>
> Beginning of processing
>
> End of processing
>
>
>
> File fich_p correponds to the Unscattered spectrum
>
> File fich_c correponds to the photons which have scattered in the crystal
>
> File fich_t correponds to the total ernergy spectrum
>
> File fich_1 correspond to the photon whose last scatter interaction was in
> the shie volume
>
> File fich_2 correspond to the photon whose last scatter interaction was in
> the coll volume
>
> File fich_3 correspond to the photon whose last scatter interaction was in
> the comp volume
>
> File fich_4 correspond to the photon whose last scatter interaction was in
> the tabl volume
>
> File fich_5 correspond to the photon whose last scatter interaction was in
> the movs volume
>
> File fich_6 correspond to the photon whose last scatter interaction was in
> the Phan volume
>
> File name fich_o1 corresponds to the first order scatter
>
> File name fich_o2 corresponds to the second order scatter
>
> File name fich_o3 corresponds to the third order scatter
>
> File name fich_o4 corresponds to the fourth order scatter
>
> File name fich_o5 corresponds to the scatter order greater than 4
>
> Number of counts for the total energy spectrum = 35795
>
> Number of counts for the Unscattered spectrum =30077
>
> Number of counts for the total scatter spectrum =5718
>
> Number of counts for the crystal scatter spectrum =3786
>
> Number of counts for the shie scatter spectrum = 1
>
> Number of counts for the coll scatter spectrum = 705
>
> Number of counts for the comp scatter spectrum = 46
>
> Number of counts for the tabl scatter spectrum = 132
>
> Number of counts for the movs scatter spectrum = 275
>
> Number of counts for the Phan scatter spectrum = 773
>
> Number of counts for the order1 scatter spectrum = 11195
>
> Number of counts for the order2 scatter spectrum = 5940
>
> Number of counts for the order3 scatter spectrum = 3013
>
> Number of counts for the order4 scatter spectrum = 1460
>
> Number of counts for scatter orders greater than 4 spectrum = 1311
>
>
>
> ##########################################################################
>
>
>
> Number of detected counts from 0.020 to 0.190 MeV : 35795
>
> Unscattered events (percentage of total counts)= 84.03
>
> Scatter in the crystal (percentage of total counts) = 10.58
>
> Scatter in the shie volume (percentage of total scatter) = 0.00
>
> Scatter in the coll volume (percentage of total scatter) = 1.97
>
> Scatter in the comp volume (percentage of total scatter) = 0.13
>
> Scatter in the tabl volume (percentage of total scatter) = 0.37
>
> Scatter in the movs volume (percentage of total scatter) = 0.77
>
> Scatter in the Phan volume (percentage of total scatter) = 2.16
>
> Scatter order1 (percentage of total scatter) = 195.79
>
> Scatter order2 (percentage of total scatter) = 103.88
>
> Scatter order3 (percentage of total scatter) = 52.69
>
> Scatter order4 (percentage of total scatter) = 25.53
>
> Scatter orders greater than 4 (percentage of total scatter) = 22.93
>
>
>
>
>
>
>
> Thanks for your kind help.
>
>
>
> best regards,
>
> Yu-Wen Yang
>
>
>
> ------------------------------
>
> _______________________________________________
> gate-users mailing list
> gate-users at lphe1pet1.epfl.ch
> http://lphe1pet1.epfl.ch/mailman/listinfo/gate-users
>
> End of gate-users Digest, Vol 8, Issue 10
> *****************************************
>
--
Assié Karine
CLCC Henri Becquerel
1 Rue d'Amiens
76 000 Rouen
France
Tel : (33) (0)2 35 14 86 17
Fax : (33) (0)2 32 08 25 50
e-mail :
karine.assie at univ-rouen.fr
More information about the Gate-users
mailing list