[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