Message: Wierd paricle getting generated in Geant4 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Wierd paricle getting generated in Geant4  

Forum: Hadronic Processes
Date: Apr 16, 13:40
From: Sanchit Sharma <Sanchit Sharma>

I am using Thermal scattering cross-sections. And my runs are getting stopped with a weird error, saying that a neitron of 24 MeV is interacting with oxygen, and we don't have crossections. My source is 14.1 MeV neutrons surrounded by water as a bulk media.

My gps is below: /run/initialize /gps/particle neutron /gps/ene/mono 14.1 MeV /gps/ang/type iso /gps/pos/type Point /gps/source/intensity 1E8 /gps/position 0.4600 0.5334 -0.1670 m

The error is as follows:
G4WT17 > 
-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had005
      issued by : G4HadronElasticProcess::PostStepDoIt
In /homes/sanchitsharma/Geant4/source/source/processes/hadronic/management/src/G4EnergyRangeManager.cc, line 128: 
===> GetHadronicInteraction: No Model found
Target element Oxygen  Z= 8  A= 16
Unrecoverable error in the method ChooseHadronicInteraction of hadElastic
TrackID= 2  ParentID= 1  neutron
Ekin(GeV)= 0.0245718;  direction= (0.231985,-0.92939,-0.287085)
Position(mm)= (575.535,378.204,-222.997);  material Water_TS
PhysicalVolume  <TESTING-MATERIAL>
 No HadronicInteraction found out

*** Fatal Exception *** core dump ***

G4WT17 > G4Track (0x7fdf24914250) - track ID = 2, parent ID = 1
G4WT17 >  Particle type : neutron - creator process : neutronInelastic, creator model : Undefined
G4WT17 >  Kinetic energy : 24.5718 MeV - Momentum direction : (0.231985,-0.92939,-0.287085)
G4WT17 >  Step length : 5.97031 cm  - total energy deposit : 0 eV 
G4WT17 >  Pre-step point : (561.685,433.692,-205.857) - Physical volume : TESTING-MATERIAL (Water_TS)
G4WT17 >  - defined by : not available
G4WT17 >  Post-step point : (575.535,378.204,-222.997) - Physical volume : TESTING-MATERIAL (Water_TS)
G4WT17 >  - defined by : hadElastic - step status : 4
G4WT17 >  *** Note: Step information might not be properly updated.
G4WT17 > 
-------- EEEE -------- G4Exception-END --------- EEEE -------

G4WT17 > 
G4WT17 > *** G4Exception: Aborting execution ***
/var/spool/slurm/d/job3400918/slurm_script: line 59: 20182 Aborted /homes/sanchitsharma/Geant4/CNEC/Thermal/build/task /homes/sanchitsharma/Geant4/CNEC/Thermal/build/macros/vis.mac

my material definition is as follows:
My material definition is as follows :
  a = 15.999 * g/mole;
  G4Element* elO = new G4Element(name="Oxygen", symbol="O",z= 8., a);
  G4Element* TSH = new G4Element( "TS_H_of_Water" , "h_water" , 1.0 , 1.0079*g/mole);
  G4Material* H2O_TS = new G4Material( "Water_TS" , density = 0.998207 * g/cm3 , ncomponents = 2, kStateLiquid, 293.15*kelvin, 1*atmosphere  );
  H2O_TS -> AddElement(TSH,2);
  H2O_TS -> AddElement(elO,1);

My physicslist.cc is attached.

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2019/04/16/13.38-47877-PhysicsList.cc

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

1 Idea: Re: Wierd paricle getting generated in Geant4   (Vladimir Ivanchenko - May 07, 09:38) newer
(_ None: Re: Wierd paricle getting generated in Geant4   (Sanchit Sharma - May 07, 09:52) newer
(_ Idea: Re: Wierd paricle getting generated in Geant4   (Vladimir Ivanchenko - May 07, 09:57) newer
 Add Message Add Message
to: "Wierd paricle getting generated in Geant4 "

 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 ]