Message: Re: G4Exception : 001 issued by : G4HadronicProcess Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: G4Exception : 001 issued by : G4HadronicProcess 

Keywords: G4Exception : 001 issued by : G4HadronicProcess
Forum: Hadronic Processes
Re: None G4Exception : 001 issued by : G4HadronicProcess (Lam YiHua)
Re: Feedback Re: G4Exception : 001 issued by : G4HadronicProcess (Gunter Folger)
Date: 01 Apr, 2007
From: Lam YiHua <LamYiHua@gmail.com>

Dear Gunter,

Please inform me if there is any error in the neutron Physics List. Thanks...

These are the dumped outputs:

_______________________________________________ ...

..

/particle/select neutron
/particle/process/dump -1
G4ProcessManager:particle[neutron]
[0]=== process[Transportation :Transportation] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                 -1:      -1:       0:       0:       5:       0:
  parameter             -1:      -1:       0:       0:       0:       0:
[1]=== process[Decay :Decay] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                  0:       0:      -1:      -1:       4:       1:
  parameter           1000:    1000:      -1:      -1:    1000:    1000:
[2]=== process[LElastic :Hadronic] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                 -1:      -1:      -1:      -1:       3:       2:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:
[3]=== process[NeutronInelastic :Hadronic] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                 -1:      -1:      -1:      -1:       2:       3:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:
[4]=== process[LFission :Hadronic] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                 -1:      -1:      -1:      -1:       1:       4:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:
[5]=== process[LCapture :Hadronic] Active
  Ordering::             AtRest             AlongStep          PostStep
                    GetPIL/    DoIt    GetPIL/    DoIt    GetPIL/    DoIt
  Ordering::
  index                 -1:      -1:      -1:      -1:       0:       5:
  parameter             -1:      -1:      -1:      -1:    1000:    1000:

..

/run/beamOn 1

msc:  Model variant of multiple scattering for e-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1
G4AugerData for Element no. 6 are loaded
G4AugerData for Element no. 7 are loaded
G4AugerData for Element no. 8 are loaded
G4AugerData for Element no. 16 are loaded
G4AugerData for Element no. 17 are loaded
G4AugerData for Element no. 29 are loaded
AugerTransitionTable complete

LowEnergyIoni:  Total cross sections from EEDL database.
      Gamma energy sampled from a parametrised formula.
      Implementation of the continuous dE/dx part.
      At present it can be used for electrons in the energy range [250eV,100GeV].
      The process must work with G4LowEnergyBremsstrahlung.

LowEnBrem:  Total cross sections from EEDL database.
      Gamma energy sampled from a parameterised formula.
      Implementation of the continuous dE/dx part.
      At present it can be used for electrons in the energy range [250eV,100GeV].
      The process must work with G4LowEnergyIonisation.

eIoni:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Delta cross sections from Moller+Bhabha, good description from 1 KeV to 100 GeV.
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

eBrem:   tables are built for  e+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Total cross sections from a parametrisation based on the EEDL data library.
      Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV.

annihil:       Sampling according eplus2gg model
      tables are built for  e+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.

msc:  Model variant of multiple scattering for proton
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1

hLowEIoni:    Knock-on electron cross sections .
        Good description above the mean excitation energy.
        Delta ray energy sampled from  differential Xsection.
        PhysicsTables from 10 eV  to 0.1 TeV  in 360 bins.
        Electronic stopping power model is  ICRU_R49p
        from 1 keV  to 2 MeV .

        Parametrisation model for antiprotons is  ICRU_R49p
        from 25 keV  to 2 MeV .
        Parametrization of the Barkas effect is switched on.
        Nuclear stopping power model is ICRU_R49

msc:  Model variant of multiple scattering for GenericIon
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1

hLowEIoni:    Knock-on electron cross sections .
        Good description above the mean excitation energy.
        Delta ray energy sampled from  differential Xsection.
        PhysicsTables from 10 eV  to 0.1 TeV  in 360 bins.
        Electronic stopping power model is  ICRU_R49p
        from 1 keV  to 2 MeV .

        Parametrisation model for antiprotons is  ICRU_R49p
        from 25 keV  to 2 MeV .
        Parametrization of the Barkas effect is switched on.
        Nuclear stopping power model is ICRU_R49

msc:  Model variant of multiple scattering for mu+
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1

muIoni:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

muBrems:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model

muPairProd:   tables are built for  mu+
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model

muIoni:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below,
      radiative corrections for E > 1 GeV
      Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1

muBrems:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model

muPairProd:   tables are built for  mu-
      dE/dx and range tables from 100 eV  to 100 TeV in 120 bins.
      Lambda tables from threshold to 100 TeV in 120 bins.
      Parametrised model

msc:  Model variant of multiple scattering for pi-
      Lambda tables from 100 eV  to 100 TeV in 120 bins.
      Boundary/stepping algorithm is active with facrange= 0.02  Step limitation 1

Region <DefaultRegionForTheWorld> -- appears in <World_phys> world volume
 Materials : Air CellDNA
 Production cuts :  gamma 1 mm     e- 1 mm     e+ 1 mm

========= Table of registered couples ==============================

Index : 0     used in the geometry : Yes     recalculation needed : No
 Material : Air
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 990 eV     e- 990 eV     e+ 990 eV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

Index : 1     used in the geometry : Yes     recalculation needed : No
 Material : CellDNA
 Range cuts        :  gamma 1 mm     e- 1 mm     e+ 1 mm
 Energy thresholds :  gamma 3.59178 keV    e- 496.074 keV    e+ 478.087 keV
 Region(s) which use this couple :
    DefaultRegionForTheWorld

====================================================================

Start closing geometry. G4GeometryManager::ReportVoxelStats -- Voxel Statistics

    Total memory consumed for geometry optimisation:   0 kByte
    Total CPU time elapsed for geometry optimisation: 0 seconds
### NAA001RunAction::BeginOfRunAction ###
                          Run 0 start...
/vis/scene/notifyHandlers
G4VisManager::SetCurrentViewer: viewer now viewer-0 (OpenGLImmediateXm)
G4VisManager::SetCurrentSceneHandler: scene handler now "scene-handler-0"
Traversing scene data...
Viewer "viewer-0 (OpenGLImmediateXm)" of scene handler "scene-handler-0"
  refreshed at request of scene "FusionScrene".
G4VisManager::SetCurrentViewer: viewer now viewer-0 (OpenGLImmediateXm)
G4VisManager::SetCurrentSceneHandler: scene handler now "scene-handler-0"
Start Run processing.
 NumberOfParticlesToBeGenerated: 1
Raw position -0.074,-0.74,0
Rotated Position (-0.074,-0.74,0)
Rotated and Translated position (-0.074,-0.74,0)
Resultant Planar wave  momentum vector (0,0,-1)
MaxwellBoltzman Energy is 2.9e-08
Creating primaries and assigning to vertex
Particle name: neutron
       Energy: 2.9e-08
     Position: (-0.074,-0.74,0)
    Direction: (0,0,-1)
 Primary Vetex generated !

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0  -0.0738   -0.741        0  2.89e-08        0        0         0 InnerPhantom_phys initStep
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
MeanFreePath = 91[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
InteractionLength= 2.1e+02[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
MeanFreePath = 1.8e+307[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
InteractionLength= 1.8e+307[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
MeanFreePath = 21[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
InteractionLength= 62[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
MeanFreePath = 0.81[cm]
G4VDiscreteProcess::PostStepGetPhysicalInteractionLength [ ]
 Particle type - neutron
   mass:        0.94[GeV]
   charge:      0[e]
   Direction x: 0, y: 0, z: -1
   Total Momentum = 7.4e-06[GeV]
   Momentum: 0[GeV], y: 0[GeV], z: -7.4e-06[GeV]
   Total Energy   = 0.94[GeV]
   Kinetic Energy = 2.9e-11[GeV]
   ProperTime     = 0[ns]
 in Material  CellDNA
InteractionLength= 0.2[cm]

 >>AlongStepDoIt (after all invocations):
    ++List of invoked processes
      1) Transportation

    ++G4Step Information
      Address of G4Track    : 0xb8319e0
      Step Length (mm)      : 1.998139590382948
      Energy Deposit (MeV)  : 0
      -----------------------------------------------------------------------
        StepPoint Information               PreStep            PostStep
      -----------------------------------------------------------------------
         Position - x (mm)   : -0.07376025683046769-0.07376025683046769
         Position - y (mm)   :  -0.7406760636610703 -0.7406760636610703
         Position - z (mm)   :                    0  -1.998139590382948
         Global Time (ns)    :                    0   849.9640007192295
         Local Time (ns)     :                    0   849.9640007192295
         Proper Time (ns)    :                    0   849.9640006930971
         Momentum Direct - x :                    0                   0
         Momentum Direct - y :                    0                   0
         Momentum Direct - z :                   -1                  -1
         Momentum - x (MeV/c):                    0                   0
         Momentum - y (MeV/c):                    0                   0
         Momentum - z (MeV/c): -0.007367696227668771-0.007367696227668771
         Total Energy (MeV)  :    939.5656300288873   939.5656300288873
         Kinetic Energy (MeV): 2.888725703087403e-082.888725703087403e-08
         Velocity (mm/ns)    : 0.0023508520227823130.002350852022782313
         Volume Name         :    InnerPhantom_phys   InnerPhantom_phys
         Safety (mm)         :                    0   0.744339711938141
         Polarization - x    :                    0                   0
         Polarization - y    :                    0                   0
         Polarization - Z    :                    0                   0
         Weight              :                    1                   1
         Step Status         :            Undefined       PostStep Proc
         Process defined Step:            Undefined            LElastic
      -----------------------------------------------------------------------

    ++List of secondaries generated (x,y,z,kE,t,PID):  No. of secodaries = 0

*** Geant4 Hadronic Reaction Information ***
    Nucleus A, Z = 16 8
    Projectile was a neutron
    projectile momentum (px, py, pz) = (9.02253e-19, 0, 0.0073677)
    Projectile energy = 939.566
*** End of Geant4 Hadronic Reaction Information ***

*** G4Exception : 001
      issued by : G4HadronicProcess
segmentation fault
*** Fatal Exception *** core dump ***

*** G4Exception: Aborting execution *** Aborted

[xxx@localhost YYY001]#

_______________________________________________

The neutron physics list:

				pManager = G4Neutron::Neutron()->GetProcessManager();

  // ELASTIC SCATTERING-----------------------------------------------------
  		 theNeutronElasticProcess = new G4HadronElasticProcess();

  G4LElastic	  *theElasticNeutronModel = new G4LElastic();
  		   			theElasticNeutronModel->SetMinEnergy(20.*MeV);  

  G4NeutronHPElastic *theElasticHPNeutron = new G4NeutronHPElastic();  
  		     		   theElasticHPNeutron->SetMinEnergy(0.0000000*MeV);  
  		     		   theElasticHPNeutron->SetMaxEnergy(19.999999999*MeV);

  		 		  theNeutronElasticProcess->RegisterMe(theElasticNeutronModel);
  		 		  theNeutronElasticProcess->RegisterMe(theElasticHPNeutron);

  G4NeutronHPElasticData *theNeutronElasticData = new G4NeutronHPElasticData();

  		 		  theNeutronElasticProcess->AddDataSet(theNeutronElasticData);

pManager->AddDiscreteProcess(theNeutronElasticProcess); G4cout<<"I have configured and tabulated Neutron Elastic Cross Section Data"<<G4endl;

  // INELASTIC SCATTERING--------------------------------------------------
  			   	   theNeutronInelasticProcess = new G4NeutronInelasticProcess();

  G4LENeutronInelastic 	   *theLENeutronModel = new G4LENeutronInelastic();			 
			  				 theLENeutronModel->SetMinEnergy(20.0*MeV);

  G4NeutronHPInelastic *theInelasticHPNeutron = new G4NeutronHPInelastic;
  			 			 theInelasticHPNeutron->SetMinEnergy(0.0000*MeV);
  			 			 theInelasticHPNeutron->SetMaxEnergy(19.999999999*MeV);

  					theNeutronInelasticProcess->RegisterMe(theLENeutronModel);
  					theNeutronInelasticProcess->RegisterMe(theInelasticHPNeutron);

  // binary
  G4BinaryCascade *neutronBC = new G4BinaryCascade;
  					neutronBC->SetMinEnergy(19.*MeV);
  					neutronBC->SetMaxEnergy(10.5*GeV);
   theNeutronInelasticProcess->RegisterMe(neutronBC);  

  // higher energy
   theNeutronInelasticProcess->RegisterMe(theTheoModel);

  // now the Neutron Inelastic Cross-sections Data //Last in First Out Rule
  G4NeutronHPInelasticData 	 	 	 *theNeutronInelasticData = new G4NeutronHPInelasticData();  
  G4NeutronInelasticCrossSection *theNeutronInelasticXsecData = new G4NeutronInelasticCrossSection();

  									theNeutronInelasticProcess->AddDataSet(theNeutronInelasticXsecData);
  									theNeutronInelasticProcess->AddDataSet(theNeutronInelasticData);

  	pManager->AddDiscreteProcess(theNeutronInelasticProcess);
  // FISSION----------------------------------------------------------------------------------------------------
  			  	 theNeutronFissionProcess = new G4HadronFissionProcess();

  G4LFission 	*theNeutronFissionModel   = new G4LFission();
  			  	 	theNeutronFissionModel->SetMinEnergy(20.0*MeV);

  G4NeutronHPFission *theFissionHPNeutron = new G4NeutronHPFission();
  				 	   theFissionHPNeutron->SetMinEnergy(0.000000*MeV);
  				 	   theFissionHPNeutron->SetMaxEnergy(19.999999999*MeV);

  				  theNeutronFissionProcess->RegisterMe(theNeutronFissionModel);
  				  theNeutronFissionProcess->RegisterMe(theFissionHPNeutron); 

  G4NeutronHPFissionData *theNeutronFissionData = new G4NeutronHPFissionData();
				  theNeutronFissionProcess->AddDataSet(theNeutronFissionData);

  pManager->AddDiscreteProcess(theNeutronFissionProcess);

  // CAPTURE----------------------------------------------------------------------------------------------------
    			theNeutronCaptureProcess = new G4HadronCaptureProcess();

  G4LCapture 	 *theNeutronCaptureModel = new G4LCapture();
  				   theNeutronCaptureModel->SetMinEnergy(20.0*MeV);

  G4NeutronHPCapture *theCaptureHPNeutron = new G4NeutronHPCapture();
  					   theCaptureHPNeutron->SetMinEnergy(0.000000*MeV);
  					   theCaptureHPNeutron->SetMaxEnergy(19.999999999*MeV);

  			theNeutronCaptureProcess->RegisterMe(theNeutronCaptureModel);
  			theNeutronCaptureProcess->RegisterMe(theCaptureHPNeutron);

  G4NeutronHPCaptureData *theNeutronCaptureData = new G4NeutronHPCaptureData();
  						theNeutronCaptureProcess->AddDataSet(theNeutronCaptureData);

  pManager->AddDiscreteProcess(theNeutronCaptureProcess);

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

 Add Message Add Message
to: "Re: G4Exception : 001 issued by : G4HadronicProcess"

 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 ]