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

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

Forum: Hadronic Processes
Date: 01 Dec, 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?   (Vladimir IVANTCHENKO - 01 Dec, 2006)
(_ None: Re: Why an Ion of H3 does not generate any secondary?   (Jorge - 01 Dec, 2006)
(_ None: Re: Why an Ion of H3 does not generate any secondary?   (Vladimir IVANTCHENKO - 01 Dec, 2006)
(_ None: Re: Why an Ion of H3 does not generate any secondary?   (Vladimir IVANTCHENKO - 01 Dec, 2006)
1 None: Re: Why an Ion of H3 does not generate any secondary?   (Jorge - 01 Dec, 2006)
3 None: Re: Why an Ion of H3 does not generate any secondary?   (Jorge - 05 Dec, 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 ]