[Gate-users] ION source

li zhi lizhi_1125 at hotmail.com
Mon Nov 20 11:07:23 CET 2006 

Dear everyone:

I have also this problem when I use ION source before. Then we changed  
GateDecayPhysics, after line #72.  It working and we do not have the long 
error now. 
But I can not see the hits visually and there is no dataset be written in 
file Gatehits.dat.  While for using Flour-18 isotope as a positron 
particle, everything works fine.
 
Could someone have some solutions or give me hints of this? 

Thank you very much!


Kind Regards
Zhi LI

 
The source I test is the same as Srinivas's:


>>----------------------------------------------
>>/gate/source/addSource F18
>>/gate/source/F18/setActivity 3700000. becquerel
>>/gate/source/F18/gps/particle e+
>>/gate/source/F18/setForcedUnstableFlag true
>>/gate/source/F18/setForcedHalfLife 6586.2 s
>>/gate/source/F18/gps/energytype Fluor18
>>/gate/source/F18/gps/type Volume
>>/gate/source/F18/gps/shape Cylinder
>>/gate/source/F18/gps/radius 40 mm
>>/gate/source/F18/gps/halfz 40 mm
>>/gate/source/F18/gps/angtype iso
>>/gate/source/F18/gps/centre 10. 0. 0. mm
>>/gate/source/F18/gps/confine NULL
>>/gate/source/F18/gps/mintheta 0. deg
>>/gate/source/F18/gps/maxtheta 180. deg
>>/gate/source/F18/gps/minphi 0. deg
>>/gate/source/F18/gps/maxphi 360. deg
>>------------------------------------------------
>>
>>Flour-18 isotope as an ion particle (SIMULATION CRASHES) :
>>
>>------------------------------------------------
>>/gate/source/addSource F18
>>/gate/source/F18/setActivity 3700000. becquerel
>>/gate/source/F18/gps/particle ion
>>/gate/source/F18/gps/ion 9 18 0 0
>>/gate/source/F18/gps/monoenergy 0. keV
>>/gate/source/F18/setForcedUnstableFlag true
>>/gate/source/F18/setForcedHalfLife 6586.2 s
>>/gate/source/F18/gps/type Volume
>>/gate/source/F18/gps/shape Cylinder
>>/gate/source/F18/gps/radius 40 mm
>>/gate/source/F18/gps/halfz 40 mm
>>/gate/source/F18/gps/angtype iso
>>/gate/source/F18/gps/centre 10. 0. 0. mm
>>/gate/source/F18/gps/confine NULL
>>/gate/source/F18/gps/mintheta 0. deg
>>/gate/source/F18/gps/maxtheta 180. deg
>>/gate/source/F18/gps/minphi 0. deg
>>/gate/source/F18/gps/maxphi 360. deg
>>---------------------------------------------------
>>


























**************************************************************************************************************




>
>----------------------------------------------------------------------
>
>Message: 1
>Date: Fri, 17 Nov 2006 11:48:42 +0000
>From: "Rebecca Gregory" <Rebecca.Gregory at icr.ac.uk>
>Subject: Re: [Gate-users] Ion particle simulation Issue
>To: "Sebastien JAN" <sebastien.jan at cea.fr>
>Cc: Gate-users at lists.healthgrid.org
>Message-ID: <s55da1b2.011 at groupwise.icr.ac.uk>
>Content-Type: text/plain; charset=US-ASCII
>
>Great thanks,
>
>It does seem to work now, the simulation runs without the Bus error in
>this case (Mac OS X 10.4 Tiger
>on intel core Duo). I still need to test this on the SUN, which gets a
>segmentation violation.
>
>I'm not that familiar with C++. I made the change in
>$GATEHOME/src/GateDecayPhysics.cc then did a make clean then make for
>the whole of GATE and had trouble with out of date libraries etc. so it
>took some time. Is there a simpler way?
>
>Becky
>
>
> >>> Sebastien JAN <sebastien.jan at cea.fr> 16 November 2006 12:32:19 >>>
>Hi,
>
>Could you test this patch
>
>In GateDecayPhysics, after line #72, insert the following line:
>
>theIonTable->GetParticle(i)->SetAtomicMass(1);     // GEANT4 Bug Fix
>#846
>
>and let me known if it works fine
>
>Cheers
>Seb
>
>
>
>Rebecca Gregory wrote:
> > Hi,
> >
> > I have also had this problem with version 3.3.
> >
> > Have you made any progress fixing this? I couldn't see a reply on
>the
> > mailing list.
> >
> > Thanks,
> > Becky
> >
> >
> > Rebecca Gregory
> > Nuclear Medicine Physicist
> > Physics Department,
> > Royal Marsden NHS Foundation Trust,
> > Downs Road, Sutton,
> > Surrey SM2 5PT, UK
> > Tel:   +44(0)20 8661 3673
> >
> >
> >>>>Sebastien JAN <jan at shfj.cea.fr> 22 August 2006 10:35:59 >>>
> >
> > Hi,
> >
> > I will check this problem asap
> > Thanks for your report
> >
> > Sebastien
> >
> > sd48 at buffalo.edu wrote:
> >
> >>Hi All,
> >>      I'am facing the following problem with regards to the
> >
> > simulation
> >
> >>of an "ion particle" isotope.
> >>
> >>My simulation crashes with the following "crash messages" when my
> >
> > source
> >
> >>isotope is an "ion particle".
> >>
> >>
> >
> >
>-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

> >
> >
> >>*** Break *** segmentation violation
> >>Generating stack trace...
> >>/usr/bin/addr2line: 'Gate': No such file
> >>/usr/bin/addr2line: 'Gate': No such file
> >>0x000000374ca90663 in std::basic_string<char,
>std::char_traits<char>,
> >
> >
> >>std::allocator<char> >::assign(std::basic_string<char,
> >>std::char_traits<char>, std::allocator<char> > const&) + 0x23 from
> >>/usr/lib64/libstdc++.so.6
> >>0x0000002a957d7c8e in GateToRoot::RecordStep(G4Step const*) + 0x18e
> >
> > from
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/gate_v3.0.0/lib/libGate.so
> >
> >>0x0000002a9574a03e in GateOutputMgr::RecordStep(G4Step const*) +
>0x7e
> >
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/gate_v3.0.0/lib/libGate.so
> >
> >>0x0000002a957b0161 in GateSteppingAction::UserSteppingAction(G4Step
> >>const*) + 0x351 from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/gate_v3.0.0/lib/libGate.so
> >
> >>0x0000002a9b21d9f3 in G4SteppingManager::Stepping() + 0x443 from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4tracking.so

> >
> >
> >>0x0000002a9b22442c in G4TrackingManager::ProcessOneTrack(G4Track*) +
> >
> >
> >>0x15c from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4tracking.so

> >
> >
> >>0x0000002a9ae7cdb9 in G4EventManager::DoProcessing(G4Event*) + 0x339
> >
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4event.so

> >
> >
> >>0x0000002a9ad5096a in G4RunManager::DoEventLoop(int, char const*,
> >
> > int) +
> >
> >>0xfa from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4run.so

> >
> >
> >>0x0000002a9ad4f061 in G4RunManager::BeamOn(int, char const*, int) +
> >
> > 0x91
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4run.so

> >
> >
> >>0x0000002a9564366f in GateApplicationMgr::StartDAQ() + 0x13f from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/gate_v3.0.0/lib/libGate.so
> >
> >>0x0000002a9564289a in
> >>GateApplicationMgrMessenger::SetNewValue(G4UIcommand*, G4String) +
> >
> > 0x14a
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/gate_v3.0.0/lib/libGate.so
> >
> >>0x0000002a9cd41e7d in G4UIcommand::DoIt(G4String) + 0x45d from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4ebeb in G4UImanager::ApplyCommand(char const*) + 0x4fb
> >
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd35354 in G4UIbatch::SessionStart() + 0x494 from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4b765 in G4UImanager::ExecuteMacroFile(char const*) +
> >
> > 0x45
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4981f in
>G4UIcontrolMessenger::SetNewValue(G4UIcommand*,
> >
> >
> >>G4String) + 0x11f from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd41e7d in G4UIcommand::DoIt(G4String) + 0x45d from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4ebeb in G4UImanager::ApplyCommand(char const*) + 0x4fb
> >
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd35354 in G4UIbatch::SessionStart() + 0x494 from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4b765 in G4UImanager::ExecuteMacroFile(char const*) +
> >
> > 0x45
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4981f in
>G4UIcontrolMessenger::SetNewValue(G4UIcommand*,
> >
> >
> >>G4String) + 0x11f from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd41e7d in G4UIcommand::DoIt(G4String) + 0x45d from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a9cd4ebeb in G4UImanager::ApplyCommand(char const*) + 0x4fb
> >
> >
> >>from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4intercoms.so

> >
> >
> >>0x0000002a98c4978d in G4UIterminal::ExecuteCommand(G4String) + 0x8d
> >
> > from
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4UIbasic.so

> >
> >
> >>0x0000002a98c4abfe in G4UIterminal::SessionStart() + 0x7e from
> >>
> >
> >
>/san1/projects/nucmed/rutaoyao/gate/gate_3.0/geant4.8.0.p01/lib/Linux-g++/libG4UIbasic.so

> >
> >
> >>0x0000000000407947 in main + 0x3b7 from Gate
> >>0x000000374b51c3fb in __libc_start_main + 0xdb from
> >
> > /lib64/tls/libc.so.6
> >
> >>0x000000000040693a in __gxx_personality_v0 + 0x7a from Gate
> >>Abort
> >>
> >
> >
>----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

> >
> >
> >>
> >>        However, if i retain all the other environmental settings
>the
> >
> >
> >>same and change the source isotope(in the source.mac) into a
> >
> > different
> >
> >>type of particle (eg: a positron), the simulation goes thru'
> >
> > properly.
> >
> >>Below are the macro files that i have used for simulation :
> >>
> >>Flour-18 isotope as a positron particle (SIMULATION RUNS WITHOUT
> >>CRASHING) :
> >>
> >>----------------------------------------------
> >>/gate/source/addSource F18
> >>/gate/source/F18/setActivity 3700000. becquerel
> >>/gate/source/F18/gps/particle e+
> >>/gate/source/F18/setForcedUnstableFlag true
> >>/gate/source/F18/setForcedHalfLife 6586.2 s
> >>/gate/source/F18/gps/energytype Fluor18
> >>/gate/source/F18/gps/type Volume
> >>/gate/source/F18/gps/shape Cylinder
> >>/gate/source/F18/gps/radius 40 mm
> >>/gate/source/F18/gps/halfz 40 mm
> >>/gate/source/F18/gps/angtype iso
> >>/gate/source/F18/gps/centre 10. 0. 0. mm
> >>/gate/source/F18/gps/confine NULL
> >>/gate/source/F18/gps/mintheta 0. deg
> >>/gate/source/F18/gps/maxtheta 180. deg
> >>/gate/source/F18/gps/minphi 0. deg
> >>/gate/source/F18/gps/maxphi 360. deg
> >>------------------------------------------------
> >>
> >>Flour-18 isotope as an ion particle (SIMULATION CRASHES) :
> >>
> >>------------------------------------------------
> >>/gate/source/addSource F18
> >>/gate/source/F18/setActivity 3700000. becquerel
> >>/gate/source/F18/gps/particle ion
> >>/gate/source/F18/gps/ion 9 18 0 0
> >>/gate/source/F18/gps/monoenergy 0. keV
> >>/gate/source/F18/setForcedUnstableFlag true
> >>/gate/source/F18/setForcedHalfLife 6586.2 s
> >>/gate/source/F18/gps/type Volume
> >>/gate/source/F18/gps/shape Cylinder
> >>/gate/source/F18/gps/radius 40 mm
> >>/gate/source/F18/gps/halfz 40 mm
> >>/gate/source/F18/gps/angtype iso
> >>/gate/source/F18/gps/centre 10. 0. 0. mm
> >>/gate/source/F18/gps/confine NULL
> >>/gate/source/F18/gps/mintheta 0. deg
> >>/gate/source/F18/gps/maxtheta 180. deg
> >>/gate/source/F18/gps/minphi 0. deg
> >>/gate/source/F18/gps/maxphi 360. deg
> >>---------------------------------------------------
> >>
> >>              I'am using the gate_v3.0.0 version for the simulation.
> >
> >
> >>Could anyone kindly provide some hints for this issue. Thank you.
> >>
> >>
> >>regards,
> >>Srinivas
> >>
> >>
> >>---------------------
> >>Srinivas Prasad D.G.
> >>Research Assistant,
> >>Department of Nuclear Medicine.
> >>University at Buffalo, State University of New York,
> >>Buffalo, NY.
> >>
> >>
> >>_______________________________________________
> >>Gate-users mailing list
> >>Gate-users at lists.healthgrid.org
> >>http://lists.healthgrid.org/mailman/listinfo/gate-users
> >>
> >>
> >
> > _______________________________________________
> > Gate-users mailing list
> > Gate-users at lists.healthgrid.org
> > http://lists.healthgrid.org/mailman/listinfo/gate-users
> > _______________________________________________
> > Gate-users mailing list
> > Gate-users at lists.healthgrid.org
> > http://lists.healthgrid.org/mailman/listinfo/gate-users
> >
> >
>
>
>------------------------------
>
>_______________________________________________
>Gate-users mailing list
>Gate-users at lists.healthgrid.org
>http://lists.healthgrid.org/mailman/listinfo/gate-users
>
>
>End of Gate-users Digest, Vol 6, Issue 14
>*****************************************

_________________________________________________________________
享用世界上最大的电子邮件系统― MSN Hotmail。 http://www.hotmail.com

More information about the Gate-users mailing list