<div dir="ltr">Hi,<br>I've attached my main macro, let me know if you have any problems to see it.<br><br clear="all"><div><div><div class="gmail_signature">Débora Salvado<br><br><br>------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 16 Jan 2015 12:18:30 +0100<br>
From: narbor <<a href="mailto:nicolas.arbor@iphc.cnrs.fr">nicolas.arbor@iphc.cnrs.fr</a>><br>
To: D?bora Salvado <<a href="mailto:debbiesal@gmail.com">debbiesal@gmail.com</a>><br>
Cc: Gate Users <<a href="mailto:gate-users@lists.opengatecollaboration.org">gate-users@lists.opengatecollaboration.org</a>><br>
Subject: Re: [Gate-users] Ion source: empty output<br>
Message-ID: <<a href="mailto:54B8F386.1060407@iphc.cnrs.fr">54B8F386.1060407@iphc.cnrs.fr</a><div id=":1vj" class="">><br>
Content-Type: text/plain; charset="windows-1252"; Format="flowed"<br>
<br>
Hi D?bora,<br>
Could you send us your main.mac file just to have a look ?<br>
Thanks.<br>
Nicolas Arbor<br>
<br>
***********************************************************************<br>
Nicolas ARBOR<br>
Ma?tre de Conf?rences<br>
Universit? de Strasbourg - UFR Physique et Ing?nierie<br>
Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien<br>
***********************************************************************<br>
<br>
Le 16/01/2015 12:09, D?bora Salvado a ?crit :<br>
> I increased the verbosity of my simulation to:<br>
><br>
> /event/verbose 1<br>
><br>
> and now when I run my simulation, it keeps giving me this message:<br>
><br>
> =====================================<br>
>   G4EventManager::ProcessOneEvent()<br>
> =====================================<br>
> G4PrimaryTransformer::PrimaryVertex (0(mm),0(mm),0(mm),1.51525e+08(nsec))<br>
> 1 primaries are passed from G4EventTransformer.<br>
> !!!!!!! Now start processing an event !!!!!!!<br>
> Track (trackID 1, parentID 0) is processed with stopping code 2<br>
> Track (trackID 3, parentID 1) is processed with stopping code 2<br>
> Track (trackID 2, parentID 1) is processed with stopping code 2<br>
> NULL returned from G4StackManager.<br>
> Terminate current event processing.<br>
><br>
> I've read that some users solved this problem by adding the<br>
> RadioactiveDecay process to their physics list. But I have it in my<br>
> macro... so can someone help me?<br>
><br>
> Thanks,<br>
> D?bora Salvado<br>
><br>
><br>
><br>
> 2015-01-15 19:22 GMT+00:00 Marc Chamberland<br>
> <<a href="mailto:MarcChamberland@cmail.carleton.ca">MarcChamberland@cmail.carleton.ca</a><br>
> <mailto:<a href="mailto:MarcChamberland@cmail.carleton.ca">MarcChamberland@cmail.carleton.ca</a>>>:<br>
><br>
>     Hmm? I don?t have experience with the sinogram output, so I?d try<br>
>     looking at the Root singles output. Is it also empty?<br>
><br>
>     Perhaps also increasing the verbosity might help you diagnose<br>
>     what?s going on. I don?t know which verbosity to increase, though!<br>
>     Perhaps:<br>
><br>
>     /event/verbose 1 (or 2 or 3)<br>
>     See if this produces any messages that might point you in the<br>
>     right direction.<br>
><br>
>     Marc<br>
><br>
><br>
><br>
>     __________________________<br>
><br>
>     Marc Chamberland, MSc<br>
>     PhD candidate<br>
>     Department of Physics<br>
>     Carleton University<br>
>     Ottawa (ON)<br>
><br>
>>     On Jan 15, 2015, at 14:04, D?bora Salvado <<a href="mailto:debbiesal@gmail.com">debbiesal@gmail.com</a><br>
>>     <mailto:<a href="mailto:debbiesal@gmail.com">debbiesal@gmail.com</a>>> wrote:<br>
>><br>
>>     Thanks Marc,<br>
>><br>
>>     I add this line to the physics list:<br>
>>     /gate/physics/addProcess RadioactiveDecay<br>
>><br>
>>     and tried to run: again no errors but an empty .sin file at the end.<br>
>><br>
>>     Any more suggestions?<br>
>><br>
>>     D?bora Salvado<br>
>><br>
>>     2015-01-15 18:36 GMT+00:00 Marc Chamberland<br>
>>     <<a href="mailto:MarcChamberland@cmail.carleton.ca">MarcChamberland@cmail.carleton.ca</a><br>
>>     <mailto:<a href="mailto:MarcChamberland@cmail.carleton.ca">MarcChamberland@cmail.carleton.ca</a>>>:<br>
>><br>
>>         Hi D?bora,<br>
>><br>
>>         The first thing I would check is that your physics list has<br>
>>         the radioactive decay process enabled.<br>
>><br>
>>         Marc<br>
>><br>
>><br>
>><br>
>>         __________________________<br>
>><br>
>>         Marc Chamberland, MSc<br>
>>         PhD candidate<br>
>>         Department of Physics<br>
>>         Carleton University<br>
>>         Ottawa (ON)<br>
>><br>
>>>         On Jan 15, 2015, at 13:30, D?bora Salvado<br>
>>>         <<a href="mailto:debbiesal@gmail.com">debbiesal@gmail.com</a> <mailto:<a href="mailto:debbiesal@gmail.com">debbiesal@gmail.com</a>>> wrote:<br>
>>><br>
>>>         Dear Gate users,<br>
>>><br>
>>>         I have my main .mac file parameterized for different sources.<br>
>>><br>
>>>         It works for the following definition of source:<br>
>>><br>
>>>         /gate/source/addSource source1<br>
>>>         /gate/source/source1/gps/type Point<br>
>>>         /gate/source/source1/gps/centre 0. 0. 0. mm<br>
>>>         /gate/source/source1/gps/particle gamma<br>
>>>         /gate/source/source1/gps/energy 159. keV<br>
>>>         /gate/source/source1/setActivity 3000000. Bq<br>
>>>         /gate/source/source1/gps/angtype iso<br>
>>><br>
>>>         However, when I try to run it for an ion source, with this<br>
>>>         definition:<br>
>>><br>
>>>         /gate/source/addSource source1<br>
>>>         /gate/source/source1/setActivity 3000000. Bq<br>
>>>         /gate/source/source1/gps/type Point<br>
>>>         /gate/source/source1/gps/centre 0. 0. 0. mm<br>
>>>         /gate/source/source1/gps/particle ion<br>
>>>         /gate/source/source1/gps/ion 53 123 0 0<br>
>>>         /gate/source/source1/setForcedUnstableFlag true<br>
>>>         /gate/source/source1/setForcedHalfLife 47520. s<br>
>>>         /gate/source/source1/gps/energytype Mono<br>
>>>         /gate/source/source1/gps/monoenergy 0. keV<br>
>>>         /gate/source/source1/gps/number 1<br>
>>>         /gate/source/source1/gps/angtype iso<br>
>>>         /gate/source/source1/gps/mintheta 0. deg<br>
>>>         /gate/source/source1/gps/maxtheta 180. deg<br>
>>>         /gate/source/source1/gps/minphi 0. deg<br>
>>>         /gate/source/source1/gps/maxphi 360. deg<br>
>>><br>
>>>         It runs but the output projection file .sin is empty...<br>
>>>         When I run with the first source instead, I get a .sin file<br>
>>>         with some data.<br>
>>>         I don't really have an error, it's just that my simulation<br>
>>>         doesn't generate anything. Can someone help me figure out<br>
>>>         what's happening?<br>
>>><br>
>>>         Regards,<br>
>>>         D?bora Salvado<br>
>>>         _______________________________________________<br>
>>>         Gate-users mailing list<br>
>>>         <a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br>
>>>         <mailto:<a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a>><br>
>>>         <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>
>><br>
>><br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> Gate-users mailing list<br>
> <a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br>
> <a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></div><br></div></div>
</div></div>