Message: ionIoni at low energies Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question ionIoni at low energies 

Forum: Physics List
Date: 20 Jun, 2013
From: <matt_b>

Hi,

I am experiencing a Segmentation fault with ionisation at lower energies on generic ions, 8Li in this case.

The following is the output from running the simulation.

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

Step# Time X Y Z KineE dEStep Volume Process

......
   70   84.6 mus -1.93 cm   -7.8 cm  -1.28 cm 0.00941 eV 0.0131 eV   scintillator  hadElastic
   71   86.8 mus -2.14 cm  -7.97 cm  -1.38 cm 0.00923 eV 0.0323 eV   scintillator  hadElastic
   72   90.2 mus -2.24 cm  -8.25 cm  -1.71 cm      0 eV      0 eV   scintillator    nCapture

*********************************************************************************************************
* G4Track Information:   Particle = Li8[0.0],   Track ID = 14,   Parent ID = 1
*********************************************************************************************************

Step#   Time         X         Y         Z        KineE    dEStep     Volume     Process
    0   90.2 mus -2.24 cm  -8.25 cm  -1.71 cm   21.3 eV      0 eV   scintillator    initStep
    1   90.2 mus -2.24 cm  -8.25 cm  -1.71 cm      0 eV   21.3 eV   scintillator     ionIoni
Segmentation fault (core dumped)

I think the relevant piece of code is in my physics list:

    else if (particleName == "GenericIon")
      {
        // OBJECT may be dynamically created as either a GenericIon or nucleus
        // G4Nucleus exists and therefore has particle type nucleus
        // genericIon:

        //multiple scattering
        pmanager->AddProcess(new G4hMultipleScattering,-1,1,1);

        //ionisation
        G4ionIonisation* ionIoni = new G4ionIonisation();
        ionIoni->SetEmModel(new G4IonParametrisedLossModel());
        ionIoni->SetStepFunction(0.1, 20*um);
        pmanager->AddProcess(ionIoni,                   -1, 2, 2);
      } 

I based my list on the DMXPhysicsList example.

Can anyone suggest what I may be missing for handling G4ionIonisation at lower energies? Any other code that may be more applicable to look at to see where the fault could be?

I'm running 4.9.6 with no patches (took me a long time to get a working install of geant and I don't want to change anything if at all possible).

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

1 Idea: Re: ionIoni at low energies   (Vladimir Ivanchenko - 20 Jun, 2013)
1 None: Re: ionIoni at low energies   (matt_b - 21 Jun, 2013)
2 None: Re: ionIoni at low energies   (matt_b - 21 Jun, 2013)
1 Idea: Re: ionIoni at low energies   (Michael H. Kelsey - 21 Jun, 2013)
... 8 Message(s)
 Add Message Add Message
to: "ionIoni at low energies"

 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 ]