Message: Re: How to add thermal neutron scattering to a production physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: How to add thermal neutron scattering to a production physics list 

Forum: Hadronic Processes
Re: Idea How to add thermal neutron scattering to a production physics list (Michael H. Kelsey)
Re: None Re: How to add thermal neutron scattering to a production physics list (Sanchit Sharma)
Re: None Re: How to add thermal neutron scattering to a production physics list (michel maire)
Date: Apr 12, 17:41
From: Sanchit Sharma <Sanchit Sharma>

Yes I did and after the implementation, I got the following error message: I am a novice Geant4 user kindly excuse me for my basic mistakes, I appreciate anyhelp in this regard:

G4WT26 > 
-------- 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 ***

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

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

I also get an output error message in my output file: G4WT29 > ASO 0.000191662 G4WT28 > ASO 1.37624e-06 G4WT10 > ASO 0.00030847 G4WT9 > ASO 3.60814e-05 G4WT1 > ASO 0.000686083 G4WT22 > ASO 0.000356742 G4WT26 > ASO 0.000739097 G4WT19 > ASO 0.000787554

G4WT2 >  ASO 6.75211e-05
G4WT2 >  ASO 4.56935e-05
G4WT26 > G4EnergyRangeManager:GetHadronicInteraction: counter=2, Ek=24.5718, Material = Water_TS, Element = Oxygen
G4WT26 > *0* low=4e-06, high=20
G4WT26 > *1* low=0, high=4e-06
G4WT26 > In /homes/sanchitsharma/Geant4/source/source/processes/hadronic/management/src/G4EnergyRangeManager.cc, line 128: 
G4WT26 > ===> GetHadronicInteraction: No Model found

My materials in detector construction is as follows:

WATER
  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 );
  H2O_TS -> AddElement(TSH,2);
  H2O_TS -> AddElement(elO,1);

HDPE
  G4Element* elTSH = new G4Element( "TS_H_of_Polyethylene" , "h_polyethylene" , 1.0 , 1.0079*g/mole );
  G4Material* HDPE = new G4Material( "Polyethylene_TS" , density = 0.970*g/cm3 , ncomponents = 2);
  HDPE -> AddElement(elTSH,2);
  HDPE -> AddElement(elC,1);

Aluminium
  G4Element* Al = new G4Element( "TS_Aluminium_Metal" , "al_metal" , 13. , 26.981539*g/mole );
  G4Material* Aluminum = new G4Material("Aluminum_TS", 2.7 * g/cm3 , 1);
  Aluminum->AddElement(Al, 100 * perCent);

 Add Message Add Message
to: "Re: How to add thermal neutron scattering to a production 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 ]