Message: Yet another radioactive decay question. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Yet another radioactive decay question. 

Keywords: co 57 decay source radioactive
Forum: Physics List
Date: 16 Dec, 2007
From: Burcin Donmez <burcind@umich.edu>

I know this has been asked a lot and there are answers in the forum. But I tried the suggestions and I still can not do simple radioactive decay of Co 57.

I looked the exRDM example and add the following to my PhysicsList file

-------------
#include "G4RadioactiveDecay.hh"
#include "G4GenericIon.hh"
void PhysicsList::ConstructGeneral()
{
  // Add Decay Process
  G4Decay* fDecayProcess = new G4Decay();
  G4RadioactiveDecay*  theRadioactiveDecay = new G4RadioactiveDecay();
  G4GenericIon* ion = G4GenericIon::GenericIon();

  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();

    if (particle == ion) {
      pmanager->AddProcess(theRadioactiveDecay, 0, -1, 3);

    } else if (fDecayProcess->IsApplicable(*particle)) { 

      pmanager ->AddProcess(fDecayProcess);

      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(fDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(fDecayProcess, idxAtRest);

    }
  } 
}
--------------
and tried the following macro (only part of it is shown)

-------------------- ## GPS

/gps/particle ion /gps/ion 27 57 0 0

#/gps/ene/type Mono
#/gps/energy 122.0 keV
/gps/energy 0 keV /gps/source/intensity 1.
# second isotope Cs137 with strength of 1.
#/gps/ion 55 137 0 662
#/gps/direction 0 0 1

#/gps/pos/type Plane
#/gps/pos/shape Square
#/gps/pos/centre 10.0 10.0 -1. mm
#/gps/pos/halfx 0.01 cm
#/gps/pos/halfy 0.01 cm

/gps/source/list

/run/beamOn 10000

---------------

and i'm getting Bus error when i tried to run after successful compile.

What am I doing wrong?

Thanks a lot.

Burçin

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

1 None: Re: Yet another radioactive decay question.   (Vladimir IVANTCHENKO - 17 Dec, 2007)
1 None: Re: Yet another radioactive decay question.   (Burcin Donmez - 17 Dec, 2007)
2 None: Re: Yet another radioactive decay question.   (Kevin Xie - 04 Jan, 2008)
 Add Message Add Message
to: "Yet another radioactive decay question."

 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 ]