Message: Re: Generating a Beam with a given Energy Probability Distribution Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: Generating a Beam with a given Energy Probability Distribution 

Keywords: Generating a Beam with a given Energy Probability Distribution
Forum: Particles
Re: Question Generating a Beam with a given Energy Probability Distribution
Re: None Re: Generating a Beam with a given Energy Probability Distribution (Maria Grazia Pia)
Re: None Re: Generating a Beam with a given Energy Probability Distribution (Giovanni Santin)
Re: Question Re: Generating a Beam with a given Energy Probability Distribution
Re: None Re: Generating a Beam with a given Energy Probability Distribution (Giovanni Santin)
Date: 02 Aug, 2006
From: <maryam@jlab.org>

Dear Giovanni,

Thanks for your great explanation. It's clear to me now how a decay is performed in geant4. I added the following to my PhysicsListGeneral.cc to associate the decay process to Sr90 ion:

   const G4IonTable *theIonTable = G4ParticleTable::GetParticleTable()->GetIonTable();
   G4RadioactiveDecay*  theRadioactiveDecay = new G4RadioactiveDecay();
  for (G4int i=0; i<theIonTable->Entries(); i++)
    {
      G4String particleName = theIonTable->GetParticle(i)->GetParticleName();
            if (particleName == "Sr90")
         {
          G4ProcessManager* pmanager = theIonTable->GetParticle(i)->GetProcessManager();
          pmanager->SetVerboseLevel(0);
          pmanager ->AddProcess(theRadioactiveDecay);
          pmanager ->SetProcessOrdering(theRadioactiveDecay, idxPostStep);
          pmanager ->SetProcessOrdering(theRadioactiveDecay, idxAtRest);
	    }
    } 

However when I start the run I don't see any betas emerge from the source. There seems to be no decay. (I need to use the beta tracks as a beam with a certain energy distribution and see the signals in a fiber beam monitor.) Setting the tacking verbose to 2, I get the following:

Event number: 0

*********************************************************************************************************
* G4Track Information:   Particle = Sr90[0.0],   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0         0        0        0         0       HallB initStep
    1        0        0        0         0        0        0         0       HallB User Limit

It looks like it's trying to track the source (?) Do you think I'm missing something? Is there anywhere else the decay process should be defined?

Thanks very very much again for all your help,

Best regards,

Maryam

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

1 None: Re: Generating a Beam with a given Energy Probability Distribution   (Vladimir IVANTCHENKO - 02 Aug, 2006)
(_ Sad: Re: Generating a Beam with a given Energy Probability Distribution   (maryam@jlab.org - 02 Aug, 2006)
(_ Ok: Re: Generating a Beam with a given Energy Probability Distribution   (maryam@jlab.org - 04 Aug, 2006)
 Add Message Add Message
to: "Re: Generating a Beam with a given Energy Probability Distribution"

 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 ]