Message: Problem with low energy EM physics list Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Problem with low energy EM physics list 

Keywords: low energy physics
Forum: Physics List
Date: 21 Apr, 2006
From: Peter Kimstrand <Peter Kimstrand>

Hello,
I'm having trouble using the low energy electromagnetic physics list for hadrons, G4hLowEnergyIonisation. 
The physics list for Em processes of my code looks like this (with standard package for EM interaction):
.
.
  else if (particleName == "proton"     ||
	     particleName == "alpha"      ||
	     particleName == "deuteron"   ||
	     particleName == "triton"     ||
	     particleName == "He3"        ||
	     particleName == "GenericIon" ) {

       G4MultipleScattering* msc = new G4MultipleScattering();
      G4hIonisation* ionisation = new G4hIonisation();

      pmanager -> AddProcess(msc, -1,1,1);   
      pmanager -> AddProcess(ionisation,-1,2,2);
      pmanager -> AddProcess(new G4StepLimiter(),-1,-1, 3);
 }

This works fine, no trouble at all. But if I change the code to:

.
.
.
  else if (particleName == "proton"     ||
	     particleName == "alpha"      ||
	     particleName == "deuteron"   ||
	     particleName == "triton"     ||
	     particleName == "He3"        ||
	     particleName == "GenericIon" ) {

       G4MultipleScattering* msc = new G4MultipleScattering();
      G4hLowEnergyIonisation* ionisation = new G4hLowEnergyIonisation();

      pmanager -> AddProcess(msc, -1,1,1);   
      pmanager -> AddProcess(ionisation,-1,2,2);
      pmanager -> AddProcess(new G4StepLimiter(),-1,-1, 3);
}

Compiling this code is OK, but when I try run the code I get a segmentation fault, the output looks like this:

Idle> 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. 18 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 10 0 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:       Heilter model of formula of annihilation into 2 photons
      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
Segmentation fault

Can anyone explain this? I've tried looking at the examples that use the Low Energy physics list (I have no trouble in compiling and runnig the examples), but I can't find the error. Help would be much apprieciated. / Peter

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

1 None: Re: Problem with low energy EM physics list   (Vladimir IVANTCHENKO - 21 Apr, 2006)
1 None: Re: Problem with low energy EM physics list   (Peter Kimstrand - 21 Apr, 2006)
(_ None: Re: Problem with low energy EM physics list   (Vladimir IVANTCHENKO - 21 Apr, 2006)
(_ None: Re: Problem with low energy EM physics list   (Peter Kimstrand - 24 Apr, 2006)
(_ None: Re: Problem with low energy EM physics list   (Gabriel Sawakuchi - 24 Apr, 2006)
1 None: Re: Problem with low energy EM physics list   (Peter Kimstrand - 24 Apr, 2006)
...
3 None: Re: Problem with low energy EM physics list   (Vladimir IVANTCHENKO - 24 Apr, 2006)
... 8 Message(s)
2 None: Re: Problem with low energy EM physics list   (Matt Enright - 06 Jul, 2006)
4 None: SegFault problem with physics list.   (Ryan Romero - 08 Aug, 2006)
1 None: Re: SegFault problem with physics list.   (Vladimir IVANTCHENKO - 09 Aug, 2006)
 Add Message Add Message
to: "Problem with low energy EM 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 ]