Message: prob in pi0 decay.............. Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question prob in pi0 decay.............. 

Forum: Physics List
Date: 16 Feb, 2010
From: Sanjib Muhuri <Sanjib Muhuri>

I have a question about taking care of pi0 decay in Geant4 for getting response in an electromagnetic calorimeter.

I have modeled an electromagnetic calorimeter where I can get response of electrons or photons at different momenta. Next I would like to get the response of two photons decaying from pi0. I have turned on the decay processes - still I don't see the decay of pi0. I am attaching part of the config file below and will appreciate receiving any suggestion in this regard. I have started my job taking exampleno3 as the basic and for the decay process I have added following to the physicslist.

////////////////////////////////////////////////////////////////////////

#include "G4Decay.hh" void ExN03PhysicsList::ConstructGeneral() {

  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4String particleName = particle->GetParticleName(); //TG
    G4ProcessManager* pmanager = particle->GetProcessManager();
    if (theDecayProcess->IsApplicable(*particle)) {

if (particleName == "pi0") {
      pmanager ->AddProcess(theDecayProcess);
      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest); }
    }
  }
}

/////////////////////////////////////////////////////////////////////////

If anybody can suggest me wheather I am doing some mistake or have to add some more steps......will be a great help for me.

 Add Message Add Message
to: "prob in pi0 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 ]