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: Hadronic Processes
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 physiclistsfEM 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 canft get any information about pionminus or pionplus Could I use hadronphysic or Bertini cascade models? but I don't kown how to write the physicslist. has someone give some instructions or can 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 ]