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

None Re: Yet another radioactive decay question. 

Forum: Physics List
Re: None Yet another radioactive decay question. (Burcin Donmez)
Date: 17 Dec, 2007
From: Vladimir IVANTCHENKO <vnivanch@mail.cern.ch>

On Sun, 16 Dec 2007, Burcin Donmez wrote:

> *** Discussion title: Physics List
> Email replies to PublicHyperNews@slac.stanford.edu must include:
>   In-Reply-To: <"/phys-list/446"@geant4-hn.slac.stanford.edu>
>   Subject: ...change this to be about your reply.
> 
> 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
> 

Hello. Bus error may not be Geant4. Try out any other G4 example 
(novice/exampleN02).

VI

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.   (Burcin Donmez - 17 Dec, 2007)
 Add Message Add Message
to: "Re: 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 ]