Message: UpdateGeometry with a HP neutron physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question UpdateGeometry with a HP neutron physics list 

Forum: Event and Track Management
Date: 10 Dec, 2010
From: Cristian Bungau <Cristian Bungau>

Dear all,

During the Idle state I want to modify the geometry with the command line. The problem is that I am using a HP neutron physics list and the neutron data has to be reloaded somehow. The lines in my UpdateGeometry are:

void DetectorConstruction::UpdateGeometry() {

  // Force geometry open
   G4GeometryManager::GetInstance()->OpenGeometry();

  // clean-up previous geometry
   G4SolidStore::GetInstance()->Clean();
   G4LogicalVolumeStore::GetInstance()->Clean();
   G4PhysicalVolumeStore::GetInstance()->Clean();

  // rebuild and notify run manager
   G4RunManager::GetRunManager()->DefineWorldVolume(Construct());
   G4RunManager::GetRunManager()->GeometryHasBeenModified();
   G4RunManager::GetRunManager()->PhysicsHasBeenModified();

  //  reset the navigator
        G4TransportationManager::GetTransportationManager()
    ->GetNavigatorForTracking()->SetWorldVolume(physiWorld);
  G4TransportationManager::GetTransportationManager()
    ->GetNavigatorForTracking()->LocateGlobalPointAndSetup
    (G4ThreeVector(),0,false);

}

But I get a segmentation fault after the geometry update:

 *** Break *** segmentation violation

===========================================================
There was a crash (#7 0xb72dfe5d in SigHandler(ESignals) ()).
This is the entire stack trace of all threads:
===========================================================
#0  0x00ca0422 in __kernel_vsyscall ()
#1  0x027197d3 in waitpid () from /lib/tls/i686/cmov/libc.so.6
#2  0x026bade3 in ?? () from /lib/tls/i686/cmov/libc.so.6
#3  0x0222f27d in system () from /lib/tls/i686/cmov/libpthread.so.0
#4  0xb72d83ed in TUnixSystem::Exec(char const*) ()
   from /home/cristian/root/lib/libCore.so
#5  0xb72dea9d in TUnixSystem::StackTrace() ()
   from /home/cristian/root/lib/libCore.so
#6  0xb72dfd55 in TUnixSystem::DispatchSignals(ESignals) ()
   from /home/cristian/root/lib/libCore.so
#7  0xb72dfe5d in SigHandler(ESignals) ()
   from /home/cristian/root/lib/libCore.so
#8  0xb72d5772 in sighandler(int) () from /home/cristian/root/lib/libCore.so
#9  <signal handler called>
#10 0x06310c97 in G4NeutronHPData::DoPhysicsVector(G4NeutronHPVector*) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4had_neu_hp.so
#11 0x063171ea in G4NeutronHPElasticData::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4had_neu_hp.so
#12 0x00c1f76f in G4CrossSectionDataStore::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_xsect.so
#13 0x003fbd9f in G4HadronicProcess::BuildPhysicsTable(G4ParticleDefinition const&) () from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_mgt.so
#14 0x0060296c in G4UHadronElasticProcess::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_coherent_elastic.so
#15 0x00329e95 in G4VUserPhysicsList::BuildPhysicsTable(G4ParticleDefinition*)
    () from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#16 0x0032a6d5 in G4VUserPhysicsList::BuildPhysicsTable() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#17 0x0031a1c3 in G4RunManagerKernel::BuildPhysicsTables() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#18 0x0031a320 in G4RunManagerKernel::RunInitialization() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#19 0x00314a55 in G4RunManager::RunInitialization() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#20 0x00312e03 in G4RunManager::BeamOn(int, char const*, int) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#21 0x0031b94a in G4RunMessenger::SetNewValue(G4UIcommand*, G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#22 0x0286b72f in G4UIcommand::DoIt(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#23 0x0287933c in G4UImanager::ApplyCommand(char const*) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#24 0x02879536 in G4UImanager::ApplyCommand(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#25 0x004888d0 in G4UIterminal::ExecuteCommand(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4UIbasic.so
#26 0x00488753 in G4UIterminal::SessionStart() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4UIbasic.so
#27 0x080502e7 in main ()
===========================================================

The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#10 0x06310c97 in G4NeutronHPData::DoPhysicsVector(G4NeutronHPVector*) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4had_neu_hp.so
#11 0x063171ea in G4NeutronHPElasticData::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4had_neu_hp.so
#12 0x00c1f76f in G4CrossSectionDataStore::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_xsect.so
#13 0x003fbd9f in G4HadronicProcess::BuildPhysicsTable(G4ParticleDefinition const&) () from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_mgt.so
#14 0x0060296c in G4UHadronElasticProcess::BuildPhysicsTable(G4ParticleDefinition const&) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4hadronic_coherent_elastic.so
#15 0x00329e95 in G4VUserPhysicsList::BuildPhysicsTable(G4ParticleDefinition*)
    () from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#16 0x0032a6d5 in G4VUserPhysicsList::BuildPhysicsTable() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#17 0x0031a1c3 in G4RunManagerKernel::BuildPhysicsTables() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#18 0x0031a320 in G4RunManagerKernel::RunInitialization() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#19 0x00314a55 in G4RunManager::RunInitialization() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#20 0x00312e03 in G4RunManager::BeamOn(int, char const*, int) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#21 0x0031b94a in G4RunMessenger::SetNewValue(G4UIcommand*, G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4run.so
#22 0x0286b72f in G4UIcommand::DoIt(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#23 0x0287933c in G4UImanager::ApplyCommand(char const*) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#24 0x02879536 in G4UImanager::ApplyCommand(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4intercoms.so
#25 0x004888d0 in G4UIterminal::ExecuteCommand(G4String) ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4UIbasic.so
#26 0x00488753 in G4UIterminal::SessionStart() ()
   from /home/cristian/geant4/geant4.9.3/lib/Linux-g++/libG4UIbasic.so
#27 0x080502e7 in main ()
===========================================================

I would be very grateful if you could please let me know what I am missing in my UpdateGeometry() to make this work.

Thank you,

Best regards,

Cristian

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Feedback: Re: UpdateGeometry with a HP neutron physics list   (Gumplinger Peter - 13 Dec, 2010)
(_ Sad: Re: UpdateGeometry with a HP neutron physics list   (Cristian Bungau - 14 Dec, 2010)
(_ Feedback: Re: UpdateGeometry with a HP neutron physics list   (Gumplinger Peter - 14 Dec, 2010)
(_ None: Re: UpdateGeometry with a HP neutron physics list   (Cristian Bungau - 15 Dec, 2010)
 Add Message Add Message
to: "UpdateGeometry with a HP neutron physics list"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]