Message: 6Li and neutrons Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None 6Li and neutrons 

Forum: Hadronic Processes
Date: 12 Mar, 2010
From: Kareem Kazkaz <Kareem Kazkaz>

Hello everyone.

I am running a simulation using ZnS / LiF as a neutron absorber. The simulation works perfectly if I don't include 6Li in the material. Unfortunately, 6Li is the important isotope in this case, so I can't simply leave it out. With the 6Li in place, the simulation crashes with this message:

G4NucleiProperties::GetNuclearMass: Wrong values for A = 1 and Z = 2

I ran the simulation using gdb, and got this backtrace:

#0  0x0321431b in G4ReactionProduct::SetDefinition ()
#1  0x034288aa in G4NeutronHPInelasticCompFS::CompositeApply ()
#2  0x03451382 in G4NeutronHPTInelasticFS::ApplyYourself ()
#3  0x033f750c in G4NeutronHPChannel::ApplyYourself ()
#4  0x033f90a0 in G4NeutronHPChannelList::ApplyYourself ()
#5  0x0341d667 in G4NeutronHPInelastic::ApplyYourself ()
#6  0x031d10d7 in G4HadronicProcess::PostStepDoIt ()
#7  0x02ef49ab in G4SteppingManager::InvokePSDIP ()
#8  0x02ef4cb2 in G4SteppingManager::InvokePostStepDoItProcs ()
#9  0x02ef2775 in G4SteppingManager::Stepping ()
#10 0x02ef9ad2 in G4TrackingManager::ProcessOneTrack ()
#11 0x02e554c4 in G4EventManager::DoProcessing ()
#12 0x02dfed92 in G4RunManager::DoEventLoop ()
#13 0x02dffff1 in G4RunManager::BeamOn ()
#14 0x02e02ed5 in G4RunMessenger::SetNewValue ()
#15 0x056f87b1 in G4UIcommand::DoIt ()
#16 0x05700da6 in G4UImanager::ApplyCommand ()
#17 0x056ec094 in G4UIbatch::ExecCommand ()
#18 0x056ec2eb in G4UIbatch::SessionStart ()
#19 0x056ff8cd in G4UImanager::ExecuteMacroFile ()
#20 0x056fed02 in G4UIcontrolMessenger::SetNewValue ()
#21 0x056f87b1 in G4UIcommand::DoIt ()
#22 0x05700da6 in G4UImanager::ApplyCommand ()
#23 0x000038ea in main ()

I am running GEANT4.9.3, and I am using the G4PhysicsListFactory to register the hadronic models. I get this crash and error message with both the QGSP_BERT_HP and QGSP_BIC_HP lists. Here is the relevant bit of code from the physics list:

	G4PhysListFactory factory;
	G4VModularPhysicsList *phys = factory.GetReferencePhysList( "QGSP_BIC_HP" );
	for( G4int i=0; ; ++i ) {
		G4VPhysicsConstructor *elem =
				const_cast<G4VPhysicsConstructor*> (phys->GetPhysics(i));
		if( elem == NULL ) break;
		G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
		RegisterPhysics( elem );
	}

Any ideas? I saw an earlier post on the subject (thread #1041 in this forum), but there weren't any responses.

Thanks for any and all guidance.

		Kareem

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

1 None: Re: 6Li and neutrons   (Kareem Kazkaz - 18 Mar, 2010)
 Add Message Add Message
to: "6Li and neutrons"

 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 ]