Message: physics list of single pion photoproduction??---geant4 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None physics list of single pion photoproduction??---geant4  

Forum: Physics List
Date: 25 Nov, 2007
From: yuchun zhang <zhangyc@ihep.ac.cn>

Hello! I am a novice at geant4.recently,I have simulated the behavior of 1-2 GeV gammas hiting a thin carbon-targets in order to get pion (that is single pion photoproduction +n/p(-/+)+p/n ) (using geant4.9.0) however I could't find any similar example ,I modified N02 and used physiclistsfEM as follows:

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

   if (particleName == "gamma") { 
     // gamma         
     pmanager->AddDiscreteProcess(new G4PhotoElectricEffect); 
     pmanager->AddDiscreteProcess(new G4ComptonScattering); 
     pmanager->AddDiscreteProcess(new G4GammaConversion); 

   } else if (particleName == "e-") { 
     //electron 
     pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); 
     pmanager->AddProcess(new G4eIonisation,       -1, 2,2); 
     pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);       

   } else if (particleName == "e+") { 
     //positron 
     pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); 
     pmanager->AddProcess(new G4eIonisation,       -1, 2,2); 
     pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3); 
     pmanager->AddProcess(new G4eplusAnnihilation,  0,-1,4); 

   } else if( particleName == "mu+" || 
              particleName == "mu-"    ) { 
     //muon   
     pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); 
     pmanager->AddProcess(new G4MuIonisation,      -1, 2,2); 
     pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 3,3); 
     pmanager->AddProcess(new G4MuPairProduction,  -1, 4,4);       

   } else if ((!particle->IsShortLived()) && 
           (particle->GetPDGCharge() != 0.0) && 
           (particle->GetParticleName() != "chargedgeantino")) { 
     //all others charged particles except geantino 
     pmanager->AddProcess(new G4MultipleScattering,-1, 1,1); 
     pmanager->AddProcess(new G4hIonisation,       -1, 2,2); 
     //step limit 
     pmanager->AddProcess(new G4StepLimiter,       -1,-1,3);         
     ///pmanager->AddProcess(new G4UserSpecialCuts,   -1,-1,4);   
   } 
 } 
} 

but I canft get any information about pionminus or pionplus Could I use hadronphysic or Bertini cascade models? but I don't kown how to write the that physicslist. has someone give some instructions or give me some examples in this aspects(single pion photoproduction)?

Thanks and best regards, yuchun zhang beijing ihep

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

 Add Message Add Message
to: "physics list of single pion photoproduction??---geant4 "

 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 ]