[Gate-users] /random/resetEngineFrom Error
Daniel Woodsworth
danieljw at phas.ubc.ca
Fri Jun 29 20:48:33 CEST 2007
Hi,
I was trying to use the the command /random/resetEngineFrom seed.rndm. In
the directory I have a a file seed.rndm which has a number 17 in it. It
appears that GATE manages to retrieve this seed value (see output below).
However, it then has an error which I don't understand at all. I have ran
many simulations with GATE and have not had any problems: I am pretty sure
that my installation of GATE is running properly. Here is the error output
from my simulation (I included a bit of the preceding output for context):
----------------------------------------
G4AugerData for Element no. 6 are loaded
G4AugerData for Element no. 7 are loaded
G4AugerData for Element no. 8 are loaded
G4AugerData for Element no. 11 are loaded
G4AugerData for Element no. 14 are loaded
G4AugerData for Element no. 18 are loaded
G4AugerData for Element no. 20 are loaded
G4AugerData for Element no. 53 are loaded
G4AugerData for Element no. 82 are loaded
AugerTransitionTable complete
*** Break *** segmentation violation
Generating stack trace...
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
0x41c00567 in CLHEP::HepJamesRandom::put() const at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/JamesRandom.cc:359
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c00fab in CLHEP::HepJamesRandom::put(std::basic_ostream<char,
std::char_traits<char> >&) const at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/JamesRandom.cc:330
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c1cde3 in CLHEP::operator<<(std::basic_ostream<char,
std::char_traits<char> >&, CLHEP::HepRandomEngine const&) at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/RandomEngine.cc:104
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c1d9a5 in CLHEP::HepRandom::saveFullState(std::basic_ostream<char,
std::char_traits<char> >&) at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/Random.cc:164 from
/home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x085e4c74 in G4RunManager::GenerateEvent(int) + 0x54 from Gate
0x085e5126 in G4RunManager::DoEventLoop(int, char const*, int) + 0xa6
from Gate
0x085e32c4 in G4RunManager::BeamOn(int, char const*, int) + 0x44 from Gate
0x081a5faf in GateApplicationMgr::StartDAQ() + 0xcf from Gate
0x081f9360 in GateApplicationMgrMessenger::SetNewValue(G4UIcommand*,
G4String)+ 0x120 from Gate
0x08859624 in G4UIcommand::DoIt(G4String) + 0x74 from Gate
0x088613eb in G4UImanager::ApplyCommand(char const*) + 0x3ab from Gate
0x0886c803 in G4UIbatch::SessionStart() + 0x2d3 from Gate
0x08861793 in G4UImanager::ExecuteMacroFile(char const*) + 0x33 from Gate
0x081e5a47 in GateUIcontrolMessenger::LaunchMacroFile(G4String) + 0x97
from Gate
0x081e5b24 in GateUIcontrolMessenger::SetNewValue(G4UIcommand*, G4String)
+ 0x84 from Gate
0x08859624 in G4UIcommand::DoIt(G4String) + 0x74 from Gate
0x088613eb in G4UImanager::ApplyCommand(char const*) + 0x3ab from Gate
0x081a4013 in main + 0x463 from Gate
0x41de2b10 in __libc_start_main + 0xe0 from /lib/tls/libc.so.6
0x081a2ac1 in TApplicationImp::ShowMembers(TMemberInspector&, char*) +
0x95 from Gate
Abort
Here is the part of my macro in which I deal with the random generator:
# RANDOM
/gate/output/root/setSaveRndmFlag 0
/random/setEngineFrom seed.rndm
As a test I commented out the first line and ran the simulation- with the
same results. I then commented out the second line and changed 0 to 1 in
the first so that I had
# RANDOM
/gate/output/root/setSaveRndmFlag 1
This gave normal GATE output.
I have pasted the complete output from my simulation below. Is this a bug
in GATE? Has anyone else encountered this?
Thanks for any help.
Daniel Woodsworth
enchmarkSPECT/st> Gate SPECTbench.mac
*************************************************************
Geant4 version Name: geant4-08-01-patch-02 (10-November-2006)
Copyright : Geant4 Collaboration
Reference : NIM A 506 (2003), 250-303
WWW : http://cern.ch/geant4
*************************************************************
I/O warning : failed to load external entity "Materials.xml"
GateMaterialDatabase: did not find the Materials.xml file: no
properties read for material 'Air'
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
Visualization Manager instantiating...
Visualization Manager initialising...
Registering graphics systems...
You have successfully registered the following graphics systems.
Current available graphics systems are:
ASCIITree (ATree)
DAWNFILE (DAWNFILE)
GAGTree (GAGTree)
G4HepRep (HepRepXML)
G4HepRepFile (HepRepFile)
RayTracer (RayTracer)
VRML1FILE (VRML1FILE)
VRML2FILE (VRML2FILE)
OpenGLImmediateX (OGLIX)
OpenGLStoredX (OGLSX)
RayTracerX (RayTracerX)
Registering model factories...
You have successfully registered the following model factories.
Registered model factories:
generic
drawByCharge
drawByParticleID
drawByOriginVolume
Registered models:
None
Registered filter factories:
chargeFilter
particleFilter
originVolumeFilter
Registered filters:
None
[GateSystemListManager::CheckScannerAutoCreation:
Creating new system based on volume inserter 'SPECThead'
I/O warning : failed to load external entity "Materials.xml"
I/O warning : failed to load external entity "Materials.xml"
I/O warning : failed to load external entity "Materials.xml"
I/O warning : failed to load external entity "Materials.xml"
GateMaterialDatabase: did not find the Materials.xml file: no
properties read for material 'Lead'
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
GateMaterialDatabase: did not find the Materials.xml file: no
properties read for material 'NaI'
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
GateMaterialDatabase: did not find the Materials.xml file: no
properties read for material 'Glass'
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
GateMaterialDatabase: did not find the Materials.xml file: no
properties read for material 'Water'
This is only a problem when OPTICAL PHOTONS are transported in this
material.'
GATE object: 'systems/SPECThead'
Components:
GATE object: 'systems/SPECThead/base'
Attached to volume: SPECThead
Nb of children: 1
GATE object: 'systems/SPECThead/crystal'
Attached to volume: crystal
Nb of children: 1
GATE object: 'systems/SPECThead/pixel'
Attached to volume: ---
Nb of children: 0
Selected processes for gamma physics:
Rayleigh effect: lowenergy
Photoelectric effect: lowenergy
Compton scattering: lowenergy
Gamma conversion: inactive
WARNING !
SetCutValueForOthers became obsolete.
It is harmless to remove this invokation without any side effects.
This dummy method implementation will be removed soon.
Confirming confinement to volume 'movsource_P'...
Confinement to volume 'movsource_P' confirmed
----- HepJamesRandom engine status -----
Initial seed = 17
u[] = 0.174685 0.883674 0.958512 0.690013 0.0250537 0.142294 0.783649
0.3560810.214774 0.69287 0.700302 0.465109 0.22215 0.768752 0.461089
0.397602 0.522215 0.121546 0.960227 0.510892 0.0116568 0.0286447 0.238777
0.585219 0.780612 0.709182 0.082106 0.828601 0.61446 0.0798209 0.688002
0.00978529 0.370831 0.376435 0.710184 0.864017 0.167338 0.21114 0.349235
0.522479 0.943713 0.0702708 0.778742 0.639458 0.506555 0.151751 0.990123
0.564688 0.430497 0.889832 0.995398 0.725368 0.710512 0.0550886 0.0540335
0.0159176 0.728379 0.947141 0.34922 0.20551 0.3510290.672539 0.118976
0.655784 0.195375 0.0380308 0.293556 0.237497 0.79723 0.158375 0.029761
0.602215 0.616657 0.995638 0.734827 0.943425 0.130007 0.277064 0.306344
0.973337 0.183926 0.547536 0.914344 0.666572 0.210074 0.300351 0.83602
0.267525 0.248549 0.533332 0.1754 0.299529 0.538308 0.555472 0.283894
0.485254 0.222273
c = 0.0216029, cd = 0.456233, cm = 1
i97 = -95, u[i97] = 1.4854e-313
j97 = -1073750680, u[j97] = 1.12088e-300
----------------------------------------
G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e-
The cut in range [1e+06 (mm)] is too big for material idx=1
The cut in energy is set1.79769e+305GeV
G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e-
The cut in range [1e+06 (mm)] is too big for material idx=2
The cut in energy is set1.79769e+305GeV
G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e-
The cut in range [1e+06 (mm)] is too big for material idx=3
The cut in energy is set1.79769e+305GeV
G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e-
The cut in range [1e+06 (mm)] is too big for material idx=4
The cut in energy is set1.79769e+305GeV
eBrem: tables are built for e-
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Total cross sections from a parametrisation based on the EEDL data
library.
Good description from 1 KeV to 100 GeV, log scale extrapolation
above 100GeV.
eIoni: tables are built for e-
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Delta cross sections from Moller+Bhabha, good description from 1 KeV
to 100 GeV.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc: Model variant of multiple scattering for e-
Lambda tables from 100 eV to 100 TeV in 120 bins.
Boundary/stepping algorithm is active with facrange= 0.02 Step
limitation 1
eBrem: tables are built for e+
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Total cross sections from a parametrisation based on the EEDL data
library.
Good description from 1 KeV to 100 GeV, log scale extrapolation
above 100GeV.
annihil: Heilter model of formula of annihilation into 2 photons
tables are built for e+
Lambda tables from 100 eV to 100 TeV in 120 bins.
annihil: Heilter model of formula of annihilation into 2 photons
tables are built for e+
Lambda tables from 100 eV to 100 TeV in 120 bins.
eIoni: tables are built for e+
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Delta cross sections from Moller+Bhabha, good description from 1 KeV
to 100 GeV.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
hIoni: tables are built for proton
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Scaling relation is used to proton dE/dx and range
Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for
protons below.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc: Model variant of multiple scattering for proton
Lambda tables from 100 eV to 100 TeV in 120 bins.
Boundary/stepping algorithm is active with facrange= 0.02 Step
limitation 1
msc: Model variant of multiple scattering for GenericIon
Boundary/stepping algorithm is active with facrange= 0.02 Step
limitation 1
hIoni: tables are built for anti_proton
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Scaling relation is used to proton dE/dx and range
Bether-Bloch model for Escaled > 2 MeV, ICRU49 parametrisation for
protons below.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
hIoni: tables are built for pi+
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Scaling relation is used to proton dE/dx and range
Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49
parametrisation forprotons below.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
hIoni: tables are built for pi-
dE/dx and range tables from 100 eV to 100 TeV in 120 bins.
Lambda tables from threshold to 100 TeV in 120 bins.
Scaling relation is used to proton dE/dx and range
Bether-Bloch model for Escaled > 0.297504 MeV, ICRU49
parametrisation forprotons below.
Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1
msc: Model variant of multiple scattering for pi-
Lambda tables from 100 eV to 100 TeV in 120 bins.
Boundary/stepping algorithm is active with facrange= 0.02 Step
limitation 1
========= Table of registered couples ==============================
Index : 0 used in the geometry : Yes recalculation needed : No
Material : Air
Range cuts : gamma 1 mm e- 1 km e+ 1 mm
Energy thresholds : gamma 990 eV e- 377.776 MeV e+ 990 eV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes recalculation needed : No
Material : Lead
Range cuts : gamma 1 mm e- 1 km e+ 1 mm
Energy thresholds : gamma 100.909 keV e- 2.88022e+295 J e+
1.28002 MeV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 2 used in the geometry : Yes recalculation needed : No
Material : NaI
Range cuts : gamma 1 mm e- 1 km e+ 1 mm
Energy thresholds : gamma 29.5738 keV e- 2.88022e+295 J e+
589.38 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 3 used in the geometry : Yes recalculation needed : No
Material : Glass
Range cuts : gamma 1 mm e- 1 km e+ 1 mm
Energy thresholds : gamma 5.708 keV e- 2.88022e+295 J e+ 568.011
keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 4 used in the geometry : Yes recalculation needed : No
Material : Water
Range cuts : gamma 1 mm e- 1 km e+ 1 mm
Energy thresholds : gamma 2.90186 keV e- 2.88022e+295 J e+
338.695 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
====================================================================
----- HepJamesRandom engine status -----
Initial seed = 17
u[] = 0.174685 0.883674 0.958512 0.690013 0.0250537 0.142294 0.783649
0.3560810.214774 0.69287 0.700302 0.465109 0.22215 0.768752 0.461089
0.397602 0.522215 0.121546 0.960227 0.510892 0.0116568 0.0286447 0.238777
0.585219 0.780612 0.709182 0.082106 0.828601 0.61446 0.0798209 0.688002
0.00978529 0.370831 0.376435 0.710184 0.864017 0.167338 0.21114 0.349235
0.522479 0.943713 0.0702708 0.778742 0.639458 0.506555 0.151751 0.990123
0.564688 0.430497 0.889832 0.995398 0.725368 0.710512 0.0550886 0.0540335
0.0159176 0.728379 0.947141 0.34922 0.20551 0.3510290.672539 0.118976
0.655784 0.195375 0.0380308 0.293556 0.237497 0.79723 0.158375 0.029761
0.602215 0.616657 0.995638 0.734827 0.943425 0.130007 0.277064 0.306344
0.973337 0.183926 0.547536 0.914344 0.666572 0.210074 0.300351 0.83602
0.267525 0.248549 0.533332 0.1754 0.299529 0.538308 0.555472 0.283894
0.485254 0.222273
c = 0.0216029, cd = 0.456233, cm = 1
i97 = -95, u[i97] = 1.4854e-313
j97 = -1073750680, u[j97] = 1.12088e-300
----------------------------------------
G4AugerData for Element no. 6 are loaded
G4AugerData for Element no. 7 are loaded
G4AugerData for Element no. 8 are loaded
G4AugerData for Element no. 11 are loaded
G4AugerData for Element no. 14 are loaded
G4AugerData for Element no. 18 are loaded
G4AugerData for Element no. 20 are loaded
G4AugerData for Element no. 53 are loaded
G4AugerData for Element no. 82 are loaded
AugerTransitionTable complete
*** Break *** segmentation violation
Generating stack trace...
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
/usr/bin/addr2line: 'Gate': No such file
0x41c00567 in CLHEP::HepJamesRandom::put() const at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/JamesRandom.cc:359
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c00fab in CLHEP::HepJamesRandom::put(std::basic_ostream<char,
std::char_traits<char> >&) const at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/JamesRandom.cc:330
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c1cde3 in CLHEP::operator<<(std::basic_ostream<char,
std::char_traits<char> >&, CLHEP::HepRandomEngine const&) at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/RandomEngine.cc:104
from /home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x41c1d9a5 in CLHEP::HepRandom::saveFullState(std::basic_ostream<char,
std::char_traits<char> >&) at
/home.mirg/danieljw/CLHEP/2.0.2.3/CLHEP/Random/src/Random.cc:164 from
/home.mirg/danieljw/CLHEP/lib/libCLHEP-2.0.2.3.so
0x085e4c74 in G4RunManager::GenerateEvent(int) + 0x54 from Gate
0x085e5126 in G4RunManager::DoEventLoop(int, char const*, int) + 0xa6
from Gate
0x085e32c4 in G4RunManager::BeamOn(int, char const*, int) + 0x44 from Gate
0x081a5faf in GateApplicationMgr::StartDAQ() + 0xcf from Gate
0x081f9360 in GateApplicationMgrMessenger::SetNewValue(G4UIcommand*,
G4String)+ 0x120 from Gate
0x08859624 in G4UIcommand::DoIt(G4String) + 0x74 from Gate
0x088613eb in G4UImanager::ApplyCommand(char const*) + 0x3ab from Gate
0x0886c803 in G4UIbatch::SessionStart() + 0x2d3 from Gate
0x08861793 in G4UImanager::ExecuteMacroFile(char const*) + 0x33 from Gate
0x081e5a47 in GateUIcontrolMessenger::LaunchMacroFile(G4String) + 0x97
from Gate
0x081e5b24 in GateUIcontrolMessenger::SetNewValue(G4UIcommand*, G4String)
+ 0x84 from Gate
0x08859624 in G4UIcommand::DoIt(G4String) + 0x74 from Gate
0x088613eb in G4UImanager::ApplyCommand(char const*) + 0x3ab from Gate
0x081a4013 in main + 0x463 from Gate
0x41de2b10 in __libc_start_main + 0xe0 from /lib/tls/libc.so.6
0x081a2ac1 in TApplicationImp::ShowMembers(TMemberInspector&, char*) +
0x95 from Gate
Abort
More information about the Gate-users
mailing list