Message: problem simulating radioactive decay Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question problem simulating radioactive decay 

Forum: Hadronic Processes
Date: 20 Oct, 2012
From: Cejo George <Cejo George>

Dear all,
  I'm having a problem simulating radioactive decay process of Cs-137 source.
The following is the terminal output I obtain.

 ======================== run summary ======================
		 The run was 20 Cs137[662.0] of 0 eV
 ===========================================================
Nb of generated particles: 
     Cs137[0.0]:      20  Emean =      0 eV 	( 0 eV  --> 0 eV )
   Cs137[662.0]:      20  Emean =      0 eV 	( 0 eV  --> 0 eV )             
e-:       1  Emean =  189.1 keV	( 189.1 keV --> 189.1 keV)
         gamma:      20  Emean =    662 keV	( 662 keV --> 662 keV)
  opticalphoton:    7156  Emean =  3.006 eV 	( 2.07 eV  --> 4.13 eV )
 Energy and momentum balance : final state - initial state
 (excluding gamma de-excitation from momentum balance) : 
  Energy:   mean =  -28.1 MeV	( -28.1 MeV --> -28.1 MeV)
  Momentum: mean =      0 eV 	( 0 eV  --> 0 eV )
 Time of life : mean =      0 ps   half-life =      0 ps    ( 0 ps  --> 0 ps )
 Activity =    inf Bq/g   (inf Ci/g) 

1.Why is Activity above inf?
2.What does these two line mean?   
   Cs137[0.0]:      20  Emean =      0 eV 	( 0 eV  --> 0 eV )
   Cs137[662.0]:      20  Emean =      0 eV 	( 0 eV  --> 0 eV )             

Below is my primary generator action file.Please correct me if there are any mistakes.

#include "PrimaryGeneratorAction.hh"
#include "G4Event.hh"
#include "G4ParticleTable.hh"
#include "G4ParticleDefinition.hh"
#include "G4Geantino.hh"
#include "Randomize.hh"
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

PrimaryGeneratorAction::PrimaryGeneratorAction() { G4int n_particle = 1; particleGun = new G4ParticleGun(n_particle); particleGun->SetParticleEnergy(0*keV); particleGun->SetParticlePosition(G4ThreeVector(0.*cm,0.*cm,2.*cm)); particleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,-1.)); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... PrimaryGeneratorAction::~PrimaryGeneratorAction() { delete particleGun; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent) { if (particleGun->GetParticleDefinition() == G4Geantino::Geantino()) { // G4int A,Z; G4int Z = 55, A = 137; G4double ionCharge = 0.*eplus; G4double excitEnergy = 662.*keV; G4ParticleDefinition* ion = G4ParticleTable::GetParticleTable()->GetIon(Z,A,excitEnergy); particleGun->SetParticleDefinition(ion); particleGun->SetParticleCharge(ionCharge); } //create vertex particleGun->GeneratePrimaryVertex(anEvent); } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... Also i didn't understand what do these parameters describe particleGun->SetParticleEnergy(0*keV);G4double ionCharge = 0.*eplus;G4double excitEnergy =662.*keV;

Best Regards, Cejo George

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

1 None: Re: problem simulating radioactive decay   (michel maire - 20 Oct, 2012)
(_ Question: Re: problem simulating radioactive decay   (Cejo George - 22 Oct, 2012)
(_ None: Re: problem simulating radioactive decay   (michel maire - 24 Oct, 2012)
(_ Question: Re: problem simulating radioactive decay   (Cejo George - 25 Oct, 2012)
(_ None: Re: problem simulating radioactive decay   (Michael H. Kelsey - 25 Oct, 2012)
 Add Message Add Message
to: "problem simulating radioactive 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 ]