Message: Heavy ion in G4HadFinalState Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Heavy ion in G4HadFinalState 

Forum: Particles
Date: 05 Jun, 2017
From: Simon Albright <Simon Albright>

I'm trying to have a G4HadFinalState return both a light ion and a heavy ion, but can only get one to work at a time.

I have the function

ApplyYourself(..., ...)
    {
    G4DynamicParticle* recoilNucleus = new G4DynamicParticle();
    recoilNucleus -> SetDefinition(ionTable -> GetIon(theTarget.GetZ_asInt(), theTarget.GetA_asInt()));
    G4DynamicParticle* returnDeut = new G4DynamicParticle();
    returnDeut -> SetDefinition(ionTable -> GetIon(1, 2));
    returnDeut -> SetMomentum(theTrack.Get4Momentum());
    theParticleChange.AddSecondary(returnDeut);
    theParticleChange.AddSecondary(recoilNucleus);
    return &theParticleChange;
    }

with

G4IonTable* ionTable = new G4IonTable();

in the header file

But this causes the ion mass to be correct and the deuteron mass to be effectively 0 (order 10^-310)

I'm obviously doing something fundamentally wrong, but have no idea what. Can anyone help?

Many thanks,

Simon Albright

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

1 Feedback: Re: Heavy ion in G4HadFinalState   (Gumplinger Peter - 05 Jun, 2017)
(_ Question: Re: Heavy ion in G4HadFinalState   (Simon Albright - 17 Jun, 2017)
 Add Message Add Message
to: "Heavy ion in G4HadFinalState"

 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 ]