[Gate-users] Unable to use phase space file from IAEA with Gate 6.2

[Stathis Kamperis]

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 at 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 at Last@h/Total at 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 <ekamperi at gmail.com> wrote:
> Relevant:
>
> http://comments.gmane.org/gmane.comp.science.opengate.user/3115
>
> Best regards,
> Stathis
>
> On Fri, Mar 7, 2014 at 11:26 PM, Stathis Kamperis <ekamperi at gmail.com> wrote:
>> Greetings everyone,
>>
>> I downloaded it from here:
>> https://www-nds.iaea.org/phsp/photon1/
>>
>> 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:
>> http://www.opengatecollaboration.org/lxr/source/6.2.0/source/physics/src/GateSourcePhaseSpace.cc#L231
>>
>> and it appears as if Gate can't recognize some particle type that the
>> header file is including. But,
>> https://www-nds.iaea.org/phsp/photon1/ELEKTA_PRECISE_6mv_part1.IAEAheader
>> has only photons, electrons, positrons that Gate should recognize.
>>
>> Any thoughts?
>>
>> Best regards,
>> Stathis