[Gate-users] Read a .IAEAphsp file : bug report ?

[REUZE SYLVAIN]

Dear Luc,

I had the same problem and finally, we decided to create a C++ code to convert IAEA files to .root files.
Z value is read by Gate (not in the header, but in the phsp file) and your particles will be at Z = 24 if you attach your PSF to this volume. You have to take this distance into account in your simulations.

Best regards,


Sylvain Reuzé
Stagiaire physicien
Département de Radiothérapie
Service de physique médicale

Tél. : 01 42 11 54 06
www.gustaveroussy.fr<http://www.gustaveroussy.fr>
[cid:image001.gif at 01D088A4.DFDE3440]
________________________________
De : gate-users-bounces at lists.opengatecollaboration.org [mailto:gate-users-bounces at lists.opengatecollaboration.org] De la part de Luc SIMON
Envoyé : mercredi 6 mai 2015 23:24
À : gate-users at lists.opengatecollaboration.org
Objet : [Gate-users] Read a .IAEAphsp file : bug report ?

Dear Gaters
A bug report and a small question

I use Gate 7.1 for Radiotherapy
 I try to read a .IAEAphsp to use it as a source (TrueBeam) and it doesn't work. I found an old email on this topic on GATE users list from Stathis but it didn't help me.

I finally fixed this problem and would like to share my solution (that is maybe not a robust one but it works for me)
For a reason I don't understand, in GateIAEARecord.cc when Gate starts to read particles :

if( fread(&ctmp, sizeof(char),   1, p_file) != 1)

the p_file pointer was not at the begining of the IAEAphsp file (but at byte 37).
I added the following:
// set the pointer at the beginning of the file
fseek(p_file,0,SEEK_SET);

just before the fread line and now it works!
I don't know if it is a bug or a newbee error... but maybe it could help someone

Now a small question:
My phase space doesn't contains Z coordinates. It is a planar phase space and the z value is constant and stored in the .IAEAheader file:
$RECORD_CONSTANT:
     24.0       // Constant Z
     1.0000     // Constant Weight
My question is: Is this z value in the header read and used by gate (I suppose no)

And if "yes" where should I attach this phase space ?

I mean if I create a box named PhS_Plane at Z = 0

and use:
/gate/source/MyBeam/attachTo PhS_Plane
what is the Z value : 0 or 24 ?
Thank you for any help.
Luc SIMON



--
___________________

SIMON Luc, Physicien Médical, PhD
Département d'Ingénierie et de Physique Médicale - IUCT-Oncopole - Toulouse, FRANCE
INSERM UMR 1037 (CRCT) - Equipe 15 - Toulouse, FRANCE
-------------- section suivante --------------
Une pi?ce jointe HTML a ?t? nettoy?e...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150507/c666391c/attachment-0001.html>
-------------- section suivante --------------
Une pi?ce jointe autre que texte a ?t? nettoy?e...
Nom: image001.gif
Type: image/gif
Taille: 6696 octets
Desc: image001.gif
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20150507/c666391c/attachment-0001.gif>