[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