[Gate-users] Ion particle simulation Issue
Rebecca Gregory
Rebecca.Gregory at icr.ac.uk
Fri Nov 17 12:48:42 CET 2006
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
>
>
More information about the Gate-users
mailing list