Message: Re: Ion decay. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Ion decay. 

Forum: Particles
Re: None Ion decay. (J.Ostariz)
Re: None Re: Ion decay. (michel maire)
Date: 05 Apr, 2011
From: J.Ostariz <J.Ostariz>

Thank you for your quick response. I have checked the environment variables and they seem to be fine.

I think that the problem arises because the recoiling ion has energy different from zero. In fact, when I set the energy of the primary ion different from zero it does not decay and I get no track. I suspect that something is lacking in my "PhysicsList":

Here, I construct the ions...

#include "G4MesonConstructor.hh"
#include "G4BaryonConstructor.hh"
#include "G4IonConstructor.hh"

void PhysicsList::ConstructHadrons()
{
    G4MesonConstructor mConstructor;  mConstructor.ConstructParticle();
    G4BaryonConstructor bConstructor; bConstructor.ConstructParticle();
    G4IonConstructor iConstructor;    iConstructor.ConstructParticle();
}

Here, I add transportation

void PhysicsList::AddTransportation()
{
    G4VUserPhysicsList::AddTransportation();

    theParticleIterator->reset();
    while((*theParticleIterator)())
    {
	G4ParticleDefinition* particle = theParticleIterator->value();
	G4ProcessManager* pmanager = particle->GetProcessManager();
	G4String particleName = particle->GetParticleName();
	pmanager->AddProcess(new G4StepLimiter, -1, -1, 1);
  }		      
}

Electromagnetic process for ions...

if(particleName == "GenericIon"||
		(particleName=="nucleus" && (particle->GetPDGCharge() != 0.0) ))
	{
	      G4MultipleScattering* aMultipleScattering 
		= new G4MultipleScattering(); 
	      pmanager->AddProcess(aMultipleScattering,-1,1,1);
	      //pmanager->AddProcess(new G4ionLowEnergyIonisation(), -1,2,2); 
              pmanager->AddProcess(new G4hIonisation(),-1,2,2);
	}

Should I add some method to my PhysicsList in order to handle ions in motion? As the RadioactiveDecay is applicable for ions with null energy, I'm not sure if the decay ion products will ever decay.

Cheers

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

1 None: Re: Ion decay.   (Sarah Bender - 03 May, 2011)
(_ None: Re: Ion decay.   (J.Ostariz - 03 May, 2011)
(_ None: Re: Ion decay.   (Sarah Bender - 04 May, 2011)
(_ None: Re: Ion decay.   (J.Ostariz - 05 May, 2011)
(_ None: Re: Ion decay.   (Sarah Bender - 05 May, 2011)
(_ Idea: Re: Ion decay.   (Vladimir Ivanchenko - 12 May, 2011)
(_ None: Re: Ion decay.   (Sarah Bender - 13 May, 2011)
 Add Message Add Message
to: "Re: Ion decay."

 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 ]