[Gate-users] problem with .IAEAphsp format

Tue Mar 14 12:44:13 CET 2017

Hi David, thanks for the tips! 

On one hand, regarding the IAEA issue that I explained below, finally I figured it out thanks to an old Gate-forum mail that I found: 

http://comments.gmane.org/gmane.comp.science.opengate.user/4333 

that had the same problem and explains how to solve it changing the GatePhaseSpaceActor.cc file. 

On the other hand, regarding the .root issue: I gather well the phase space in .root format, but when I use that file as a source qstaitself, it crashes and says "need to define the particle type", which I don't understand because all that info is already gathered in the .root file as I double check. Anyway, I'm gonna check the cc files that you say. 

Thanks! 
Consuelo. 

De: "David Sarrut" <David.Sarrut at creatis.insa-lyon.fr> 
À: "Consuelo Guardiola" <guardiola at imnc.in2p3.fr>, "gate-users" <gate-users at lists.opengatecollaboration.org> 
Cc: "Prezado Yolanda" <prezado at imnc.in2p3.fr> 
Envoyé: Lundi 13 Mars 2017 10:44:52 
Objet: Re: [Gate-users] problem with .IAEAphsp format 

Hello, 
- for IAEA phase-space -> sorry, I dont know 

- for root phase-space, please explain precisely what the issue is. When compiling ? When writing or when reading a phsp ? What is you goal ? 

Short information: 

To write a phsp, it is managed by the [ https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GatePhaseSpaceActor.cc | PhaseSpaceActor ] . There are a lot of options (see the macros in the [ https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GatePhaseSpaceActorMessenger.cc | messenger ] ). Do you setEnablePartName to true to store the particle name in the phsp ? 

To read a phsp as a source of particles, it is managed by the [ https://github.com/OpenGATE/Gate/blob/develop/source/physics/src/GateSourcePhaseSpace.cc | SourcePhaseSpace ] . Also, have a look at the code and the associated messengers to get the full list of options. 

Phase Space is a large piece of code and could be quite powerful, but it is complex. Please note that some bug fixes and modifications have been done since the last Gate version. You may have a better luck to give a try with the current Gate version (develop branch, which will be the next Gate 8.0 soon). 

David 
PS: please keep gate-user mailing list in copy, it could help other people ... and other people could help ;) 

On Mon, Mar 13, 2017 at 10:28 AM, Consuelo Guardiola < [ mailto:guardiola at imnc.in2p3.fr | guardiola at imnc.in2p3.fr ] > wrote: 

Hi David (I'm adding Yolanda Prezado in CC, my advisor here, since she told me she has recently discussed some v7.0 vs 7.1 discrepancies with you) : 

I lately tried with IAEA phsp format since in the Gate forum someone replied me this: 

-------------------------------------------------------------------------------------------------- 

Hi, 

I had the same problem few months ago and i think it was related to the root version. Version 5.34/36 didn't work but 5.34/26 worked. Otherwise, i'm finally using .IAEAphsp format to gather PhaseSpace files and its also working. 

Hope this would help, 

Tony 

when I asked this: 

BQ_BEGIN

Dear Gate users, 

I gathered a Phase Space file (in .root format) to be used as a source afterwards. When I try to run it as it is explained in the manual: 

[ http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.2:Phase_space_concept | http://wiki.opengatecollaboration.org/index.php/Users_Guide_V7.2:Phase_space_concept ] 

it fails and says that I need to define the particle type. Nevertheless, the .root file contains all the information (particle types, etc). I have not found more info neither in forums or around. Please could someone tell me if I am losing some info or what I may do to figure it out? 

Thanks in advance, 
Consuelo. 
----------------------------------------------------------------------------------------------------------------------- 

I have tried both things: first reinstalling everything with 5.34/26 root version, but did not work, later with IAEA phsp format... the point is that neither with Gate 7.0 nor 7.1 I can launch the phase space file WITHOUT specify the particle type (which would be completely wrong...). Do you have any idea how to figure this out? I would appreciate whatever advice. 

Thanks in advance 
Consuelo 

De: "David Sarrut" < [ mailto:David.Sarrut at creatis.insa-lyon.fr | David.Sarrut at creatis.insa-lyon.fr ] > 
À: "Consuelo Guardiola" < [ mailto:guardiola at imnc.in2p3.fr | guardiola at imnc.in2p3.fr ] > 
Envoyé: Lundi 13 Mars 2017 10:18:36 
Objet: Re: [Gate-users] problem with .IAEAphsp format 

Hello, 
to my knowledge, the IAEA phsp form in Gate in buggy: it only works for some type of particles. 

You may look at the source code: 
- [ https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAHeader.cc | https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAHeader.cc ] 
- [ https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAPhsp.cc | https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAPhsp.cc ] 
- [ https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEARecord.cc | https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEARecord.cc ] 
- [ https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAUtilities.cc | https://github.com/OpenGATE/Gate/blob/develop/source/general/src/GateIAEAUtilities.cc ] 

it could help to find what's wrong ... 

David 
PS: sorry I cannot help, I do not use such format 

On Wed, Mar 8, 2017 at 3:34 PM, Consuelo Guardiola < [ mailto:guardiola at imnc.in2p3.fr | guardiola at imnc.in2p3.fr ] > wrote: 

BQ_BEGIN

Dear Gate users, 

If I try to gather the phase space in .IAEAphsp format in both Gate v7.0 and v7.1, it does not work. Apparently, it does not get any error-file back, but in the running-file I find this: 

GatePhaseSpaceActor.cc (l.296): Actor phase space: particle not available in IAEA format. 

(no matter the version of Gate that I use or even running different jobs to discard a possible particle/material issue) 

But, unexpectedly, I get an empty (0 size) output.IAEAphsp file. 

Does anyone know what happens? 

Thanks in advance, 
Consuelo. 

_______________________________________________ 
Gate-users mailing list 
[ mailto:Gate-users at lists.opengatecollaboration.org | Gate-users at lists.opengatecollaboration.org ] 
[ http://lists.opengatecollaboration.org/mailman/listinfo/gate-users | http://lists.opengatecollaboration.org/mailman/listinfo/gate-users ] 

BQ_END

-- 
David Sarrut, Phd 
Directeur de recherche CNRS 
CREATIS, UMR CNRS 5220, Inserm U1206 
Centre de lutte contre le cancer Léon Bérard 
28 rue Laënnec, 69373 Lyon cedex 08 
Tel : [ tel:04%2078%2078%2051%2051 | 04 78 78 51 51 ] / [ tel:06%2074%2072%2005%2042 | 06 74 72 05 42 ] 
[ http://www.creatis.insa-lyon.fr/~dsarrut | http://www.creatis.insa-lyon.fr/~dsarrut ] 
_________________________________ 
"2 + 2 = 5, for extremely large values of 2" 
_________________________________ 

BQ_END

-- 
David Sarrut, Phd 
Directeur de recherche CNRS 
CREATIS, UMR CNRS 5220, Inserm U1206 
Centre de lutte contre le cancer Léon Bérard 
28 rue Laënnec, 69373 Lyon cedex 08 
Tel : 04 78 78 51 51 / 06 74 72 05 42 
[ http://www.creatis.insa-lyon.fr/~dsarrut | http://www.creatis.insa-lyon.fr/~dsarrut ] 
_________________________________ 
"2 + 2 = 5, for extremely large values of 2" 
_________________________________ 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20170314/2e05e552/attachment.html>