[Gate-users] Ion source: empty output

narbor nicolas.arbor at iphc.cnrs.fr
Fri Jan 16 12:18:30 CET 2015


Hi Débora,
Could you send us your main.mac file just to have a look ?
Thanks.
Nicolas Arbor

***********************************************************************
Nicolas ARBOR
Maître de Conférences
Université de Strasbourg - UFR Physique et Ingénierie
Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien
***********************************************************************

Le 16/01/2015 12:09, Débora Salvado a écrit :
> I increased the verbosity of my simulation to:
>
> /event/verbose 1
>
> and now when I run my simulation, it keeps giving me this message:
>
> =====================================
>   G4EventManager::ProcessOneEvent()
> =====================================
> G4PrimaryTransformer::PrimaryVertex (0(mm),0(mm),0(mm),1.51525e+08(nsec))
> 1 primaries are passed from G4EventTransformer.
> !!!!!!! Now start processing an event !!!!!!!
> Track (trackID 1, parentID 0) is processed with stopping code 2
> Track (trackID 3, parentID 1) is processed with stopping code 2
> Track (trackID 2, parentID 1) is processed with stopping code 2
> NULL returned from G4StackManager.
> Terminate current event processing.
>
> I've read that some users solved this problem by adding the 
> RadioactiveDecay process to their physics list. But I have it in my 
> macro... so can someone help me?
>
> Thanks,
> Débora Salvado
>
>
>
> 2015-01-15 19:22 GMT+00:00 Marc Chamberland 
> <MarcChamberland at cmail.carleton.ca 
> <mailto:MarcChamberland at cmail.carleton.ca>>:
>
>     Hmm… I don’t have experience with the sinogram output, so I’d try
>     looking at the Root singles output. Is it also empty?
>
>     Perhaps also increasing the verbosity might help you diagnose
>     what’s going on. I don’t know which verbosity to increase, though!
>     Perhaps:
>
>     /event/verbose 1 (or 2 or 3)
>     See if this produces any messages that might point you in the
>     right direction.
>
>     Marc
>
>
>
>     __________________________
>
>     Marc Chamberland, MSc
>     PhD candidate
>     Department of Physics
>     Carleton University
>     Ottawa (ON)
>
>>     On Jan 15, 2015, at 14:04, Débora Salvado <debbiesal at gmail.com
>>     <mailto:debbiesal at gmail.com>> wrote:
>>
>>     Thanks Marc,
>>
>>     I add this line to the physics list:
>>     /gate/physics/addProcess RadioactiveDecay
>>
>>     and tried to run: again no errors but an empty .sin file at the end.
>>
>>     Any more suggestions?
>>
>>     Débora Salvado
>>
>>     2015-01-15 18:36 GMT+00:00 Marc Chamberland
>>     <MarcChamberland at cmail.carleton.ca
>>     <mailto:MarcChamberland at cmail.carleton.ca>>:
>>
>>         Hi Débora,
>>
>>         The first thing I would check is that your physics list has
>>         the radioactive decay process enabled.
>>
>>         Marc
>>
>>
>>
>>         __________________________
>>
>>         Marc Chamberland, MSc
>>         PhD candidate
>>         Department of Physics
>>         Carleton University
>>         Ottawa (ON)
>>
>>>         On Jan 15, 2015, at 13:30, Débora Salvado
>>>         <debbiesal at gmail.com <mailto:debbiesal at gmail.com>> wrote:
>>>
>>>         Dear Gate users,
>>>
>>>         I have my main .mac file parameterized for different sources.
>>>
>>>         It works for the following definition of source:
>>>
>>>         /gate/source/addSource source1
>>>         /gate/source/source1/gps/type Point
>>>         /gate/source/source1/gps/centre 0. 0. 0. mm
>>>         /gate/source/source1/gps/particle gamma
>>>         /gate/source/source1/gps/energy 159. keV
>>>         /gate/source/source1/setActivity 3000000. Bq
>>>         /gate/source/source1/gps/angtype iso
>>>
>>>         However, when I try to run it for an ion source, with this
>>>         definition:
>>>
>>>         /gate/source/addSource source1
>>>         /gate/source/source1/setActivity 3000000. Bq
>>>         /gate/source/source1/gps/type Point
>>>         /gate/source/source1/gps/centre 0. 0. 0. mm
>>>         /gate/source/source1/gps/particle ion
>>>         /gate/source/source1/gps/ion 53 123 0 0
>>>         /gate/source/source1/setForcedUnstableFlag true
>>>         /gate/source/source1/setForcedHalfLife 47520. s
>>>         /gate/source/source1/gps/energytype Mono
>>>         /gate/source/source1/gps/monoenergy 0. keV
>>>         /gate/source/source1/gps/number 1
>>>         /gate/source/source1/gps/angtype iso
>>>         /gate/source/source1/gps/mintheta 0. deg
>>>         /gate/source/source1/gps/maxtheta 180. deg
>>>         /gate/source/source1/gps/minphi 0. deg
>>>         /gate/source/source1/gps/maxphi 360. deg
>>>
>>>         It runs but the output projection file .sin is empty...
>>>         When I run with the first source instead, I get a .sin file
>>>         with some data.
>>>         I don't really have an error, it's just that my simulation
>>>         doesn't generate anything. Can someone help me figure out
>>>         what's happening?
>>>
>>>         Regards,
>>>         Débora Salvado
>>>         _______________________________________________
>>>         Gate-users mailing list
>>>         Gate-users at lists.opengatecollaboration.org
>>>         <mailto:Gate-users at lists.opengatecollaboration.org>
>>>         http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>>
>>
>
>
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users

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


More information about the Gate-users mailing list