<br>Hi,<br>Thanks for your help.<br><br>I can run benchmark_spectra and benchmark_projections now.<br>But the results obtained was different from the results showed in the manual.<br>Is it allright, or something wrong somewhere?
<br>Followed are the messages:<br>
<p class="MsoNormal"><span style="font-size: 10pt;">stolen:~/gate/gate_v3.0.0/benchmarks/benchmarkSPECT>
./benchmark_spectra gateSingles.dat</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the lower energy
threshold (MeV)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">0.02</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Lower energy threshold is
0.020000 MeV</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the upper energy
threshold (MeV)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">0.19</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Upper energy threshold is
0.190000 MeV</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the number of energy
windows</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">100</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Energy range will be divided
into 100 energy windows </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the number of
scattering volume (max is 10)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">6</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">6 scattering volumes will be
considered</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 1 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">shie</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 2 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">coll</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 3 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">comp</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 4 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">tabl</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 5 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">movs</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the first four letters
of the scattering volume number 6 </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">phan </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">The scattering volumes that
will be considered are the following:</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">shie</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">coll</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">comp</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">tabl</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">movs</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">phan</span></p>
<br>
<p class="MsoNormal"><span style="font-size: 10pt;">Beginning of processing</span></p><span style="font-size: 10pt;"></span><span style="font-size: 10pt;">End of processing</span>
<p class="MsoNormal"><span style="font-size: 10pt;"> <br></span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_p correponds to the
Unscattered spectrum</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_c correponds to the
photons which have scattered in the crystal</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_t correponds to the
total ernergy spectrum</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_1 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the shie volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_2 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the coll volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_3 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the comp volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_4 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the tabl volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_5 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the movs volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File fich_6 correspond to the
photon whose last scatter interaction was<span style="">
</span>in the phan volume</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File name fich_o1 corresponds
to the first order scatter</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File name fich_o2 corresponds
to the second order scatter</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File name fich_o3 corresponds
to the third order scatter</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File name fich_o4 corresponds
to the fourth order scatter</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File name fich_o5 corresponds
to the scatter order greater than 4</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
total energy spectrum = 36041</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
Unscattered spectrum =12809</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
total scatter spectrum =23232</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
crystal scatter spectrum =2392</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the shie
scatter spectrum = 1</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the coll
scatter spectrum = 107</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the comp
scatter spectrum = 456</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the tabl
scatter spectrum = 1081</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the movs
scatter spectrum = 5123</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the phan
scatter spectrum = 0</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
order1 scatter spectrum = 11176</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
order2 scatter spectrum = 6120</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
order3 scatter spectrum = 3004</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for the
order4 scatter spectrum = 1500</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of counts for scatter
orders greater than 4 spectrum = 1432</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;"> </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">##########################################################################</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;"> </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of detected counts
from 0.020 to 0.190 MeV : 36041</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Unscattered events
(percentage of total counts)= 35.54</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the crystal
(percentage of total counts) = 6.64</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the shie volume
(percentage of total scatter) = 0.00</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the coll volume
(percentage of total scatter) = 0.30</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the comp volume
(percentage of total scatter) = 1.27</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the tabl volume
(percentage of total scatter) = 3.00</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the movs volume
(percentage of total scatter) = 14.21</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter in the phan volume
(percentage of total scatter) = 0.00</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter order1 (percentage of
total scatter) = 48.11</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter order2 (percentage of
total scatter) = 26.34</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter order3 (percentage of
total scatter) = 12.93</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter order4 (percentage of
total scatter) = 6.46</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Scatter orders greater than 4
(percentage of total scatter) = 6.16</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;"> </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">##########################################################################</span></p>
<br><br><p class="MsoNormal"><span style="font-size: 11pt;">For benchmark_projections;<br></span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">stolen:~/gate/gate_v3.0.0/benchmarks/benchmarkSPECT>
./benchmark_projections gateSingles.dat</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the number of rows and
columns of the images to be created</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">128</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Image size is 128 x 128</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the lower energy
threshold (MeV)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">0.02</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Lower energy threshold is
0.020000 MeV</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the upper energy
threshold (MeV)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">0.19</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Upper energy threshold is
0.190000 MeV</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Enter the size of the pixels
(mm)</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">0.904</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Pixels are 0.904000 mm x
0.904000 mm x</span></p><span style="font-size: 10pt;"><br></span><span style="font-size: 10pt;"></span><span style="font-size: 10pt;">Beginning of processing</span><span style="font-size: 10pt;"> </span>
<p class="MsoNormal"><span style="font-size: 10pt;">End of processing</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;"> <br>Number of detected events in
the image of primary and scattered photons = 17619</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of primary events in
the image of primary photons= 6201</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">Number of scattered
events<span style=""> </span>in the image of scattered photons
= 11418</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;"> </span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File Total correponds to the
image of primary and scattered photons</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File Primary correponds to
the image of primary photons</span></p>
<p class="MsoNormal"><span style="font-size: 10pt;">File Scatter correponds to
the image of scattered photons</span></p>
<br><br>Appreciate your reply and help.<br>Many thanks,<br>S<br><br><div><span class="gmail_quote">On 5/19/07, <b class="gmail_sendername">Daniel Woodsworth</b> <<a href="mailto:danieljw@phas.ubc.ca">danieljw@phas.ubc.ca
</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Suhaili,<br><br>See below for ideas.<br><br>Cheers,<br><br>Daniel<br>
<br>> Hi all,<br>> Can someone help me with these problems?<br>><br>> When I tried to run benchmark macros (i.e. SPECTbench.mac and<br>> benchmarkPET.mac), these problems occur:<br>><br>> 1. There is no file for
benchSPECT.root.<br>><br>Are you sure you loaded root using root benchSPECT.root ? That tells root<br>to 'mount' benchSPECT.root.<br>> When I typed root [0] .x benchSPECT.C<br>><br>> I have this message: Error in <TFile::TFile> file does not exist
<br>><br>> Error: illegal pointer to class object<br>> Singles 0´0 133 benchSPECT.C:21<br>><br>><br>><br>> 1. When tried to complie benchmark_spectra.c using:<br>>
<br>> gcc -o benchmark_spectra benchmark_spectra.c –lm<br>><br>> I got this message:<br>><br>> benchmark_spectra.c: In function â:<br>><br>> benchmark_spectra.c:768: error: missing terminating " character
<br>><br>> benchmark_spectra.c:769: error: â undeclared (first use in this function)<br>><br>> benchmark_spectra.c:769: error: (Each undeclared identifier is reported<br>> only<br>> once<br>><br>> benchmark_spectra.c:769: error: for each function it appears in.)
<br>><br>> benchmark_spectra.c:769: error: missing terminating " character<br>><br>> benchmark_spectra.c:771: error: syntax error before â token<br>><br>I had this exact same problem. If you open up benchmark_spectra.c in a
<br>text editor and go to probably line 768 (but check to be sure) you should<br>see that the line ends and then 769 is a new line. But it should be one<br>continuous piece of code. But the compiler is "seeing" a new line of code,
<br>which of course is nonsense. If you just delete out any spacing between<br>the two lines (but make sure to preserve the syntax) and recompile it<br>should work.<br><br>><br>><br>> 1. Same goes with benchmark_spectra.c
<br>><br>> benchmark_projections.c: In function â:<br>><br>> benchmark_projections.c:662: error: missing terminating " character<br>><br>> benchmark_projections.c:663: error: â undeclared (first use in this
<br>> function)<br>><br>> benchmark_projections.c:663: error: (Each undeclared identifier is<br>> reported<br>> only once<br>><br>> benchmark_projections.c:663: error: for each function it appears in.)
<br>><br>> benchmark_projections.c:663: error: syntax error before â<br>><br>> benchmark_projections.c:663: error: stray â in program<br>><br>> benchmark_projections.c:663: error: missing terminating " character
<br>><br>exactly the same problem and solution as for benchmark_spectra.c<br><br>><br>><br>> 1. Does x. file.C file always appear after simulation runs? And how to<br>> convert .root file to .C file?<br>
><br>><br>> 1. For benchmarkPET.mac, are there no files for benchmark_spectra.c<br>> and benchmark_projections.c? Why?<br>><br>><br><br>when I did this I just ran<br><br>danieljw$ root benchSPECT.root
<br><br>root [0] .x benchSPECT.C<br><br>it should work (just make sure you have appropriate files in your working<br>directory)<br><br><br>> I really appreciate your help and explanation. Looking forward to your<br>> replies.
<br>> Many thanks,<br>> S<br>><br>> * *<br>> _______________________________________________<br>> Gate-users mailing list<br>> <a href="mailto:Gate-users@lists.healthgrid.org">Gate-users@lists.healthgrid.org
</a><br>> <a href="http://lists.healthgrid.org/mailman/listinfo/gate-users">http://lists.healthgrid.org/mailman/listinfo/gate-users</a><br>><br><br><br></blockquote></div><br>