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

None Re: Problem with low energy EM physics list 

Forum: Physics List
Re: None Problem with low energy EM physics list (Peter Kimstrand)
Date: 21 Apr, 2006
From: Vladimir IVANTCHENKO <vnivanch@mail.cern.ch>

On Fri, 21 Apr 2006, Peter Kimstand wrote:

> *** Discussion title: Physics List
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/phys-list/253"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> 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
> 

It is not garanteed any of your PhysicsList are correct.

Please, look into some official examples.

$G4INSTALL/examples/extended/electromagnetic/TestEm7/src/PhysListEmStandard.cc

shows optimal physics for standard. If you need low-energy option for 
hadrons/ion, then replace both (!) G4hIonisation and G4ionIonisation by 
G4hLowEnergyIonisation.

VI


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   (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)
1 Sad: Re: Problem with low energy EM physics list   (Mathieu Trocmé - 19 Sep, 2009)
2 Sad: Re: Problem with low energy EM physics list   (Mathieu Trocmé - 19 Sep, 2009)
... 6 Message(s)
 Add Message Add Message
to: "Re: 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 ]