Hi there,
I have been simulating ions of C14(average energy 49keV) and S35
(average energy 48-53 keV) with the GPS with no problem so far. These
ions decay into 3 secondaries, 1 beta electron (the one that I am really
interested), 1 anti_nu_e and 1 ion of another element.
The point is that now I have to do the same simulations with H3 or triton
(average energy 5.7 keV) but the code does not generate absolutely any
secondary.
I thought that it was because of the Cuts. I changed them in the Physics
List:
void Exp2PhysicsList::SetCuts()
{
defaultCutValue = 0.001*um;
SetCutsWithDefault();
}
... but no success. Then I read that these cuts are only applied to
EM processes and in this case, to use the GPS and GenericIons, I had
to register the G4Decay process to the GenericIon. So these cuts
are not applied to this particle (correct me if I'm wrong).
Then I thought the it could be because I had to register the low energy
processes for the electrons instead of the normal ones:
if (particleName == "e-") {
//electron
/*pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
pmanager->AddProcess(new G4eIonisation, -1, 2,2);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 3,3);*/
loweIon = new G4LowEnergyIonisation("LowEnergyIoni");
loweBrem = new G4LowEnergyBremsstrahlung("LowEnBrem");
loweBrem->SetAngularGenerator("tsai");
pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
pmanager->AddProcess(loweIon, -1, 2,2);
pmanager->AddProcess(loweBrem, -1,3,3);
}
... but no success.
I think that the secondaries are being killed somewhere in the code
but I don't know where. If I only could find somewhere in the code
where the secondaries are still alive and try to enable something
to let them live, like the flag SetGoodForTrackingFlag or something
like that.
Please, help me to keep these secondaries safe from the code. They
have the right to live as well.
Thanks a lot.
|