Message: simulate gamma + proton -> pion + proton Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question simulate gamma + proton -> pion + proton 

Forum: Hadronic Processes
Date: 14 Jun, 2010
From: farzad.totonchi <farzad.totonchi>

Hi I am trying to simulate the pi0 photoproduction from gamma near proton with threshold energy 144.7 MeV (gamma + proton -> pion + proton ). I use following code for gamma process :

############################

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

if (particleName == "gamma") {

     //EM
     pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
     pmanager->AddDiscreteProcess(new G4ComptonScattering);
     pmanager->AddDiscreteProcess(new G4GammaConversion);

    //Hadronic
   G4PhotoNuclearProcess* photoProc = new G4PhotoNuclearProcess("PhotonInelastic");
   G4GammaNuclearReaction* gammaNuclear = new G4GammaNuclearReaction();
   photoProc->RegisterMe(gammaNuclear);
   pmanager->AddDiscreteProcess(photoProc);
}

#############################

when run program do'nt produce the pi0 i think this my method is wrong if class can produce this process please introduce to me.

 Add Message Add Message
to: "simulate gamma + proton -> pion + proton"

 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 ]