[Gate-users] Ion source: empty output

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


Could it comes from vGATE itself ?
I never used it so I don't know how physics lists work on it...
Nicolas

***********************************************************************
Nicolas ARBOR
Maître de Conférences
Université de Strasbourg - UFR Physique et Ingénierie
Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien
23, rue du Loess BP28 67037 Strasbourg Cedex 2
mail : nicolas.arbor at iphc.cnrs.fr / tel : 03.88.10.64.27
***********************************************************************

Le 16/01/2015 15:24, Marc Chamberland a écrit :
> Hi Débora,
>
> I tried your macro on Gate v7.0 with Geant4.9.6 and it worked (see 
> Root output attached). I don’t know what could be going wrong on your 
> system.
> Marc
>
>
>
>> On Jan 16, 2015, at 08:32, Débora Salvado <debbiesal at gmail.com 
>> <mailto:debbiesal at gmail.com>> wrote:
>>
>> Dear GATE users,
>>
>> I've tried Uwe's suggestion and I keep getting lots of messages saying:
>>  =====================================
>>  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 also tried Nicolas' suggestion about the physics process but it 
>> says command not found (is it not included on vGATE? that's strange).
>>
>> Any more suggestions?
>>
>> Regards,
>> Débora Salvado
>>
>>
>> 2015-01-16 13:08 GMT+00:00 Débora Salvado <debbiesal at gmail.com 
>> <mailto:debbiesal at gmail.com>>:
>>
>>
>>     2015-01-16 12:37 GMT+00:00 Nicolas ARBOR
>>     <nicolas.arbor at iphc.cnrs.fr <mailto:nicolas.arbor at iphc.cnrs.fr>>:
>>
>>         This is strange...may be build physics list don't worked with
>>         vGate (?).
>>         Have you try modifications from Uwe ? Does it work ?
>>         Nicolas
>>
>>         ***********************************************************************
>>         Nicolas ARBOR
>>         Maître de Conférences
>>         Université de Strasbourg - UFR Physique et Ingénierie
>>         Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien
>>         23, rue du Loess BP28 67037 Strasbourg Cedex 2
>>         mail :nicolas.arbor at iphc.cnrs.fr  <mailto:nicolas.arbor at iphc.cnrs.fr>  / tel : 03.88.10.64.27
>>         ***********************************************************************
>>
>>         Le 16/01/15 13:27, Débora Salvado a écrit :
>>>         It gives me an error:
>>>         ***** COMMAND NOT FOUND </gate/physics/addPhysicsList
>>>         emstandard_opt3> *****
>>>
>>>         Débora Salvado
>>>
>>>         2015-01-16 11:50 GMT+00:00 narbor
>>>         <nicolas.arbor at iphc.cnrs.fr
>>>         <mailto:nicolas.arbor at iphc.cnrs.fr>>:
>>>
>>>             Just to check...does it changes something if you use for
>>>             the physics (instead of what you have) ?
>>>             /gate/physics/addPhysicsList emstandard_opt3
>>>             /gate/physics/addProcess RadioactiveDecay
>>>             Nicolas
>>>
>>>             ***********************************************************************
>>>             Nicolas ARBOR
>>>             Maître de Conférences
>>>             Université de Strasbourg - UFR Physique et Ingénierie
>>>             Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien
>>>             23, rue du Loess BP28 67037 Strasbourg Cedex 2
>>>             mail :nicolas.arbor at iphc.cnrs.fr  <mailto:nicolas.arbor at iphc.cnrs.fr>  / tel : 03.88.10.64.27
>>>             ***********************************************************************
>>>
>>>             Le 16/01/2015 12:40, Débora Salvado a écrit :
>>>>             Hi,
>>>>             I've attached my main macro, let me know if you have
>>>>             any problems to see it.
>>>>
>>>>             Débora Salvado
>>>>
>>>>
>>>>             ------------------------------
>>>>
>>>>             Message: 2
>>>>             Date: Fri, 16 Jan 2015 12:18:30 +0100
>>>>             From: narbor <nicolas.arbor at iphc.cnrs.fr
>>>>             <mailto:nicolas.arbor at iphc.cnrs.fr>>
>>>>             To: D?bora Salvado <debbiesal at gmail.com
>>>>             <mailto:debbiesal at gmail.com>>
>>>>             Cc: Gate Users
>>>>             <gate-users at lists.opengatecollaboration.org
>>>>             <mailto:gate-users at lists.opengatecollaboration.org>>
>>>>             Subject: Re: [Gate-users] Ion source: empty output
>>>>             Message-ID: <54B8F386.1060407 at iphc.cnrs.fr
>>>>             <mailto:54B8F386.1060407 at iphc.cnrs.fr>
>>>>             >
>>>>             Content-Type: text/plain; charset="windows-1252";
>>>>             Format="flowed"
>>>>
>>>>             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>
>>>>             > <mailto: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>
>>>>             >>  <mailto: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>
>>>>             >>  <mailto: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>
>>>>             <mailto: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>
>>>>             >>>  <mailto: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
>>>>             <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  <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 
>> <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/fd7bcb63/attachment-0001.html>


More information about the Gate-users mailing list