Message: Re: How to gereate unstable Phosphor Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Idea Re: How to gereate unstable Phosphor 

Keywords: Radioactive decay
Forum: Particles
Re: Question How to gereate unstable Phosphor (Mario Schubert)
Date: 15 Feb, 2002
From: Alexander Dierlamm <Alexander Dierlamm>

To generate an unstable ion I did the following:

Define the particles and decay tables in the PhysicsList (here for Sr90): 

void ExN03PhysicsList::ConstructSource()
{
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4ParticleDefinition* particle;   


 // Sr90
   particle = particleTable->GetIon(38,90,0.0);
   particle->SetPDGStable(0);
   particle->SetPDGLifeTime(9.082e8);
 
  G4VDecayChannel* channel = new G4BetaMinusDecayChannel(2,particle,1.,.546,0,1,false,NULL);
  G4DecayTable* table=new G4DecayTable();
  table->Insert(channel);
  particle->SetDecayTable(table);
 
  //Y90
   particle = particleTable->GetIon(39,90,0.0);
   particle->SetPDGStable(0);
   particle->SetPDGLifeTime(2.308e5);
  
   table=new G4DecayTable();

   channel = new G4BetaMinusDecayChannel(2,particle,1,2.2801,0,1,false,NULL);
   table->Insert(channel);

   particle->SetDecayTable(table);
}

And for the generation in PrimaryGeneratorAction:


  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();

  G4ParticleDefinition* particle;
  particle= particleTable->GetIon(38,90,0.0);
  //particle->DumpTable(); 
 
  particleSource->SetParticleDefinition(particle);

This should generate Sr90 which decays into Y90 and this further in Zr90. 

There might be an easier solution and I would be glad to learn an easier way
by accessing the predefined RadioactiveDecay data or something else...

Cheers, Alexander

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

 Add Message Add Message
to: "Re: How to gereate unstable Phosphor"

 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 ]