Message: G4AntiProtonAnnihilationAtRest Issues Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question G4AntiProtonAnnihilationAtRest Issues 

Keywords: AntiProton
Forum: Hadronic Processes
Date: 06 Nov, 2007
From: <slaug634@erau.edu>

GEANT4 Users/Developers, I am a new user to GEANT4 and am encountering some issues. While running a simulation of annihilating antiprotons, I am having an error with the produced secondary particles. One example is my proton antiproton annihilation is producing 2 gamma-rays at 938 MeV! Obviously something is wrong with my code:

The only place I call upon the G4AntiProtonAnnihilationAtRest process is listed in my code below. Any help at all would be greatly appreciated.

void MyPhysicsList::ConstructEM()
{
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();

     if(particle->GetParticleName() == "anti_proton") {
       G4AntiProtonAnnihilationAtRest* apAnnihil = 
	  new G4AntiProtonAnnihilationAtRest();
      pmanager->AddRestProcess(apAnnihil);
	  G4AntiProtonInelasticProcess* theInelasticProcess = new G4AntiProtonInelasticProcess("inelastic");
      G4LEAntiProtonInelastic* theLEInelasticModel = new G4LEAntiProtonInelastic;
      theInelasticProcess->RegisterMe(theLEInelasticModel);
      G4HEAntiProtonInelastic* theHEInelasticModel = new G4HEAntiProtonInelastic;
      theInelasticProcess->RegisterMe(theHEInelasticModel);
      pmanager->AddDiscreteProcess(theInelasticProcess);

}

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

1 None: Re: G4AntiProtonAnnihilationAtRest Issues   (Vladimir IVANTCHENKO - 07 Nov, 2007)
 Add Message Add Message
to: "G4AntiProtonAnnihilationAtRest Issues"

 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 ]