<div dir="ltr"><div>Dear Chunhu,</div><div> </div><div>It happened because the amount of yttrium in the Gate materials in very high. </div><div>for example for 10% yttrium you have to make change as below:</div><div> </div> <div>LYSO: d=7.105 g/cm3; n=4 ; state=Solid +el: name=Lutetium ; f=0.714609801 +el: name=Yttrium ; f=0.040381125 +el: name=Silicon; f=0.063520871 +el: name=Oxygen; f=0.181488203</div> <div> </div><div> </div><div>Regards,</div><div>Tahereh</div></div><div class="gmail_extra"> <div class="gmail_quote">On Sun, Mar 9, 2014 at 3:33 AM, <<a href="mailto:gate-users-request@lists.opengatecollaboration.org" target="_blank">gate-users-request@lists.opengatecollaboration.org</a>> wrote: <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Gate-users mailing list submissions to <a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a> To subscribe or unsubscribe via the World Wide Web, visit <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a> or, via email, send a message with subject or body 'help' to <a href="mailto:gate-users-request@lists.opengatecollaboration.org">gate-users-request@lists.opengatecollaboration.org</a> You can reach the person managing the list at <a href="mailto:gate-users-owner@lists.opengatecollaboration.org">gate-users-owner@lists.opengatecollaboration.org</a> When replying, please edit your Subject line so it is more specific than "Re: Contents of Gate-users digest..." Today's Topics: 1. Re: Unable to use phase space file from IAEA with Gate 6.2 (Stathis Kamperis) 2. Re: Unable to use phase space file from IAEA with Gate 6.2 (Stathis Kamperis) 3. About the performance comparison between LSO ans LYSO (???) ---------------------------------------------------------------------- Message: 1 Date: Sat, 8 Mar 2014 22:13:46 +0200 From: Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> To: gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>> Subject: Re: [Gate-users] Unable to use phase space file from IAEA with Gate 6.2 Message-ID: <CACYCkUtbj8p9gE72CC2No7B0n03wXrF=3q9BaQ=<a href="mailto:edMbEfKae1Q@mail.gmail.com">edMbEfKae1Q@mail.gmail.com</a>> Content-Type: text/plain; charset=ISO-8859-1 I'm not sure why Gate can't parse the phase space file. I tried to extract the information with Mathematica and it seems that the phase space file is ok (not sure about the header though, especially where the record contents are declared). For example I get only 1,2,3 (gammas, e-, e+) for particles: file = BinaryReadList[ "~/Dropbox/foo", Join[ Flatten@ {"Byte", PadLeft[{}, 6, ToString@Real32], "Integer32", "Integer32"}], ByteOrdering -> -1]; h = HistogramList[file[[All, 1]]] {{1, 2, 3, 4}, {1543049, 8306, 161}} And the ratios seem to be in alignment with what is declared at the header: First@Last@h/Total@Last@h // N 0.994543 125806703/126491813 // N 0.994584 I've only checked 50 MB worth of data but I doubt that there's a wrong particle number in the rest (will check though). So... has anyone ever successfully used the IAEA phase space files with Gate 6.2 ? Best regards, Stathis On Fri, Mar 7, 2014 at 11:40 PM, Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> wrote: > Relevant: > > <a href="http://comments.gmane.org/gmane.comp.science.opengate.user/3115" target="_blank">http://comments.gmane.org/gmane.comp.science.opengate.user/3115</a> > > Best regards, > Stathis > > On Fri, Mar 7, 2014 at 11:26 PM, Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> wrote: >> Greetings everyone, >> >> I downloaded it from here: >> <a href="https://www-nds.iaea.org/phsp/photon1/" target="_blank">https://www-nds.iaea.org/phsp/photon1/</a> >> >> ELEKTA_PRECISE_6mv_part1.IAEAphsp and its associated header >> ELEKTA_PRECISE_6mv_part1.IAEAheader >> >> When I try to use it from within Gate with: >> >> /gate/source/addSource mybeam phaseSpace >> /gate/source/mybeam/addPhaseSpaceFile >> /home/stathis/Downloads/ELEKTA_PRECISE_6mv_part1.IAEAphsp >> >> I get the following error: >> /home/stathis/Downloads/gate_v6.2/source/physics/src/GateSourcePhaseSpace.cc >> (l.252): Source phase space: particle not available in IAEA phase >> space format. >> >> I checked the source code: >> <a href="http://www.opengatecollaboration.org/lxr/source/6.2.0/source/physics/src/GateSourcePhaseSpace.cc#L231" target="_blank">http://www.opengatecollaboration.org/lxr/source/6.2.0/source/physics/src/GateSourcePhaseSpace.cc#L231</a> >> >> and it appears as if Gate can't recognize some particle type that the >> header file is including. But, >> <a href="https://www-nds.iaea.org/phsp/photon1/ELEKTA_PRECISE_6mv_part1.IAEAheader" target="_blank">https://www-nds.iaea.org/phsp/photon1/ELEKTA_PRECISE_6mv_part1.IAEAheader</a> >> has only photons, electrons, positrons that Gate should recognize. >> >> Any thoughts? >> >> Best regards, >> Stathis ------------------------------ Message: 2 Date: Sat, 8 Mar 2014 22:34:36 +0200 From: Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> To: gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>> Subject: Re: [Gate-users] Unable to use phase space file from IAEA with Gate 6.2 Message-ID: <<a href="mailto:CACYCkUuFkakuzc-WGG9%2BfpZsjreKK0oaw7MGgfFbwDZkeiB2HA@mail.gmail.com">CACYCkUuFkakuzc-WGG9+fpZsjreKK0oaw7MGgfFbwDZkeiB2HA@mail.gmail.com</a>> Content-Type: text/plain; charset=ISO-8859-1 On Sat, Mar 8, 2014 at 10:13 PM, Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> wrote: > I'm not sure why Gate can't parse the phase space file. I tried to > extract the information with Mathematica and it seems that the phase > space file is ok (not sure about the header though, especially where > the record contents are declared). > > For example I get only 1,2,3 (gammas, e-, e+) for particles: > > file = > BinaryReadList[ > "~/Dropbox/foo", > Join[ > Flatten@ > {"Byte", PadLeft[{}, 6, ToString@Real32], "Integer32", "Integer32"}], > ByteOrdering -> -1]; > > h = HistogramList[file[[All, 1]]] > {{1, 2, 3, 4}, {1543049, 8306, 161}} > > And the ratios seem to be in alignment with what is declared at the header: > > First@Last@h/Total@Last@h // N > 0.994543 > > 125806703/126491813 // N > 0.994584 > > I've only checked 50 MB worth of data but I doubt that there's a wrong > particle number in the rest (will check though). > The whole phase space file is ok. I checked it with the following code snippet: [stathis@fedora]~/Downloads% cat read.c 22:32:29 #include <assert.h> #include <stdio.h> struct rec { char particle; char pad[32]; } __attribute__((packed)); int main(int argc, char *argv[]) { FILE *fp; struct rec r; int cnt, rv; fp = fopen(argv[1], "r"); assert(fp); cnt = 0; do { ++cnt; rv = fread(&r, sizeof(struct rec), 1, fp); assert(rv > 0 || (rv == 0 && feof(fp))); if (r.particle != 1 && r.particle != 2 && r.particle != 3) { printf("Wrong particle = %d at record = %d\n", r.particle, cnt); } } while(!feof(fp) && !ferror(fp)); fclose(fp); return 0; } [stathis@fedora]~/Downloads% gcc -Wall -ansi -pedantic read.c [stathis@fedora]~/Downloads% ./a.out ELEKTA_PRECISE_6mv_part1.IAEAphsp [stathis@fedora]~/Downloads% [stathis@fedora]~/Downloads% echo $? 0 [stathis@fedora]~/Downloads% > So... has anyone ever successfully used the IAEA phase space files > with Gate 6.2 ? > > Best regards, > Stathis > > > > > On Fri, Mar 7, 2014 at 11:40 PM, Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> wrote: >> Relevant: >> >> <a href="http://comments.gmane.org/gmane.comp.science.opengate.user/3115" target="_blank">http://comments.gmane.org/gmane.comp.science.opengate.user/3115</a> >> >> Best regards, >> Stathis >> >> On Fri, Mar 7, 2014 at 11:26 PM, Stathis Kamperis <<a href="mailto:ekamperi@gmail.com">ekamperi@gmail.com</a>> wrote: >>> Greetings everyone, >>> >>> I downloaded it from here: >>> <a href="https://www-nds.iaea.org/phsp/photon1/" target="_blank">https://www-nds.iaea.org/phsp/photon1/</a> >>> >>> ELEKTA_PRECISE_6mv_part1.IAEAphsp and its associated header >>> ELEKTA_PRECISE_6mv_part1.IAEAheader >>> >>> When I try to use it from within Gate with: >>> >>> /gate/source/addSource mybeam phaseSpace >>> /gate/source/mybeam/addPhaseSpaceFile >>> /home/stathis/Downloads/ELEKTA_PRECISE_6mv_part1.IAEAphsp >>> >>> I get the following error: >>> /home/stathis/Downloads/gate_v6.2/source/physics/src/GateSourcePhaseSpace.cc >>> (l.252): Source phase space: particle not available in IAEA phase >>> space format. >>> >>> I checked the source code: >>> <a href="http://www.opengatecollaboration.org/lxr/source/6.2.0/source/physics/src/GateSourcePhaseSpace.cc#L231" target="_blank">http://www.opengatecollaboration.org/lxr/source/6.2.0/source/physics/src/GateSourcePhaseSpace.cc#L231</a> >>> >>> and it appears as if Gate can't recognize some particle type that the >>> header file is including. But, >>> <a href="https://www-nds.iaea.org/phsp/photon1/ELEKTA_PRECISE_6mv_part1.IAEAheader" target="_blank">https://www-nds.iaea.org/phsp/photon1/ELEKTA_PRECISE_6mv_part1.IAEAheader</a> >>> has only photons, electrons, positrons that Gate should recognize. >>> >>> Any thoughts? >>> >>> Best regards, >>> Stathis ------------------------------ Message: 3 Date: Sun, 9 Mar 2014 10:07:39 +0800 From: ??? <<a href="mailto:arnoxiaolei86@gmail.com">arnoxiaolei86@gmail.com</a>> To: "<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>" <<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>> Subject: [Gate-users] About the performance comparison between LSO ans LYSO Message-ID: <CAGMGNVN=Rx8YpX=<a href="mailto:Qg5Rq%2ByDbO%2B1YjY0tM47fDDHMYW1bRgVaFA@mail.gmail.com">Qg5Rq+yDbO+1YjY0tM47fDDHMYW1bRgVaFA@mail.gmail.com</a>> Content-Type: text/plain; charset="iso-8859-1" Hello everyone: I did two experiments using GateV6.2, The conditions of these two experiments are same with each other except the material of the crystal (using LSO in first experiment and LYSO in secone one), but the Coincidences number recorded in the first experiment is treble as many as that recorded in the second experiment. What a bad performance of LYSO! Can anyone tell me why? The macros are in the attachment. Thanks a lot!--Chunhui Zhang -------------- next part -------------- A non-text attachment was scrubbed... Name: PET.zip Type: application/zip Size: 113372 bytes Desc: not available URL: <<a href="http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20140309/f3829bbd/attachment.zip" target="_blank">http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20140309/f3829bbd/attachment.zip</a>> ------------------------------ _______________________________________________ Gate-users mailing list <a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a> <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a> End of Gate-users Digest, Vol 94, Issue 15 ****************************************** </blockquote></div> </div>