Message: Why an Ion of H3 does not generate any secondary? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Why an Ion of H3 does not generate any secondary? 

Keywords: Why an Ion of H3 does not generate any secondary?
Forum: Event and Track Management
Date: 22 Nov, 2006
From: Jorge Cabello <Jorge Cabello>

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.

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

1 None: Re: Why an Ion of H3 does not generate any secondary?   (Makoto Asai - 30 Nov, 2006)
 Add Message Add Message
to: "Why an Ion of H3 does not generate any secondary?"

 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 ]