Message: Re: Simulation of an AmBe-neutron source Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Simulation of an AmBe-neutron source 

Forum: Hadronic Processes
Re: Question Simulation of an AmBe-neutron source (Martin Kirchner)
Re: None Re: Simulation of an AmBe-neutron source (Artem Zontikov)
Re: More Re: Simulation of an AmBe-neutron source (Martin Kirchner)
Re: None Re: Simulation of an AmBe-neutron source (Artem Zontikov)
Date: 14 Jun, 2016
From: Martin Kirchner <Martin Kirchner>

Thanks again for your reply and sorry for the late response!

ad 3.: I tried this, and other configurations of the physics list, but as long as G4ParticleHPInelastic was in use for alphas the problem persisted, the Sim did not run really stable and the results didn't seem to improve. So I reverted to the initial physics list.

ad 1. + 2.: I added the JENDL-3.3 data files and changed to BeO, which improved the results (attached). But i'm not really satisfied.

I tried to implement AmO2 grains spread around the BeO, but it would have needed far too many such grains(~10^9) to be realistic.

Now I wanted to try it with a realistic amount AmO2 added to the BeO. But with a Material definition like this:

	Cresst::Material::MaterialProperty SrcMat;
        SrcMat.Name = "BeOAmO2";
	SrcMat.Material = new G4Material(SrcMat.Name, 3.71 * g / cm3, 3); 		
	SrcMat.Material->AddMaterial(BeO.Material, 92.5*perCent);
	SrcMat.Material->AddMaterial(AmO2.Material, 7.5*perCent);

... the Sim produces a segmentation violation:

=========================================================== There was a crash. This is the entire stack trace of all threads: ===========================================================

#0  0x00000034300ac65e in waitpid () from /lib64/libc.so.6
#1  0x000000343003e609 in do_system () from /lib64/libc.so.6
#2  0x00007fc7092085aa in TUnixSystem::StackTrace() () from /afs/hephy.at/project/cresst/software/root-6.04.00/lib/root/libCore.so.6.05
#3  0x00007fc70920a54c in TUnixSystem::DispatchSignals(ESignals) () from /afs/hephy.at/project/cresst/software/root-6.04.00/lib/root/libCore.so.6.05
#4  <signal handler called>
#5  0x00007fc70a7c82df in G4Physics2DVector::Value(double, double, unsigned long&, unsigned long&) const () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/global/management/src/G4Physics2DVector.cc:182
#6  0x00007fc70c413a45 in G4SeltzerBergerModel::ComputeDXSectionPerAtom(double) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/standard/src/G4SeltzerBergerModel.cc:215
#7  0x00007fc70c432765 in G4eBremsstrahlungRelModel::ComputeXSectionPerAtom(double) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/standard/src/G4eBremsstrahlungRelModel.cc:333
#8  0x00007fc70c432826 in G4eBremsstrahlungRelModel::ComputeCrossSectionPerAtom(G4ParticleDefinition const*, double, double, double, double, double) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/standard/src/G4eBremsstrahlungRelModel.cc:296
#9  0x00007fc70c44f530 in G4EmElementSelector::Initialise(G4ParticleDefinition const*, double) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/utils/include/G4VEmModel.hh:538
#10 0x00007fc70c475d3d in G4VEmModel::InitialiseElementSelectors(G4ParticleDefinition const*, G4DataVector const&) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/utils/src/G4VEmModel.cc:205
#11 0x00007fc70c4310c6 in G4eBremsstrahlungRelModel::Initialise(G4ParticleDefinition const*, G4DataVector const&) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/standard/src/G4eBremsstrahlungRelModel.cc:182
#12 0x00007fc70c45329c in G4EmModelManager::Initialise(G4ParticleDefinition const*, G4ParticleDefinition const*, double, int) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/utils/src/G4EmModelManager.cc:545
#13 0x00007fc70c48743c in G4VEnergyLossProcess::PreparePhysicsTable(G4ParticleDefinition const&) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/processes/electromagnetic/utils/src/G4VEnergyLossProcess.cc:643
#14 0x00007fc70e302f3d in G4VUserPhysicsList::PreparePhysicsTable(G4ParticleDefinition*) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4VUserPhysicsList.cc:755
#15 0x00007fc70e30322f in G4VUserPhysicsList::BuildPhysicsTable() () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4VUserPhysicsList.cc:530
#16 0x00007fc70e2ebb92 in G4RunManagerKernel::BuildPhysicsTables(bool) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4RunManagerKernel.cc:707
#17 0x00007fc70e2ec09e in G4RunManagerKernel::RunInitialization(bool) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4RunManagerKernel.cc:609
#18 0x00007fc70e2dbcb7 in G4RunManager::RunInitialization() () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4RunManager.cc:313
#19 0x00007fc70e2d618d in G4RunManager::BeamOn(int, char const*, int) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4RunManager.cc:272
#20 0x00007fc70e2f227d in G4RunMessenger::SetNewValue(G4UIcommand*, G4String) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/run/src/G4RunMessenger.cc:376
#21 0x00007fc70aa0fe12 in G4UIcommand::DoIt(G4String) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UIcommand.cc:230
#22 0x00007fc70aa285dd in G4UImanager::ApplyCommand(char const*) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UImanager.cc:523
#23 0x00007fc70a9ff997 in G4UIbatch::ExecCommand(G4String const&) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UIbatch.cc:170
#24 0x00007fc70aa01043 in G4UIbatch::SessionStart() () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UIbatch.cc:215
#25 0x00007fc70aa227b3 in G4UImanager::ExecuteMacroFile(char const*) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UImanager.cc:296
#26 0x00007fc70aa1945c in G4UIcontrolMessenger::SetNewValue(G4UIcommand*, G4String) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UIcontrolMessenger.cc:313
#27 0x00007fc70aa0fe12 in G4UIcommand::DoIt(G4String) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UIcommand.cc:230
#28 0x00007fc70aa285dd in G4UImanager::ApplyCommand(char const*) () at /afs/hephy.at/project/cresst/software/geant4.10.02.p01/source/intercoms/src/G4UImanager.cc:523
#29 0x0000000000402895 in main ()
===========================================================

For the small grains AmO2 worked as Material, but here it is making problems.

Any ideas how to fix this problem or to improve the simulation in any other way? Best regards, Martin

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/06/14/16.38-21674-BeO_run0527_outofcube.jpg

 Add Message Add Message
to: "Re: Simulation of an AmBe-neutron source"

 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 ]