G4EmStandardFysik.cc

Go to the documentation of this file.

#include "G4EmStandardFysik.hh"

#include "G4EmStandardPhysics.hh"
 #include "G4ParticleDefinition.hh"
 #include "G4ProcessManager.hh"
 #include "G4LossTableManager.hh"
 #include "G4EmProcessOptions.hh"
 
 #include "G4ComptonScattering.hh"
 #include "G4GammaConversion.hh"
 #include "G4PhotoElectricEffect.hh"
 
 #include "G4eMultipleScattering.hh"
 #include "G4MuMultipleScattering.hh"
 #include "G4hMultipleScattering.hh"
 #include "G4CoulombScattering.hh"
 #include "G4WentzelVIModel.hh"
 
 #include "G4eIonisation.hh"
 #include "G4eBremsstrahlung.hh"
 #include "G4eplusAnnihilation.hh"
 
 #include "G4MuIonisation.hh"
 #include "G4MuBremsstrahlung.hh"
 #include "G4MuPairProduction.hh"
 #include "G4hBremsstrahlung.hh"
 #include "G4hPairProduction.hh"
 
 #include "G4hIonisation.hh"
 #include "G4ionIonisation.hh"
 #include "G4alphaIonisation.hh"
 
 #include "G4Gamma.hh"
 #include "G4Electron.hh"
 #include "G4Positron.hh"
 #include "G4MuonPlus.hh"
 #include "G4MuonMinus.hh"
 #include "G4PionPlus.hh"
 #include "G4PionMinus.hh"
 #include "G4KaonPlus.hh"
 #include "G4KaonMinus.hh"
 #include "G4Proton.hh"
 #include "G4AntiProton.hh"
 #include "G4Deuteron.hh"
 #include "G4Triton.hh"
 #include "G4He3.hh"
 #include "G4Alpha.hh"
#include "G4GenericIon.hh"

G4EmStandardFysik::G4EmStandardFysik(G4int ver)
{
        verbose = ver;
}

void G4EmStandardFysik::ConstructProcess()
{ 
     theParticleIterator->reset();
     while( (*theParticleIterator)() ){
     G4ParticleDefinition* particle = theParticleIterator->value();
     G4ProcessManager* pmanager = particle->GetProcessManager();
     G4String particleName = particle->GetParticleName();
     if(verbose > 1)
       G4cout << "### " << GetPhysicsName() << " instantiates for " 
              << particleName << G4endl;
 
     if (particleName == "gamma") {
 
       pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
      // pmanager->AddDiscreteProcess(new G4ComptonScattering);
       pmanager->AddDiscreteProcess(new G4GammaConversion);
 
     } else if (particleName == "e-") {
 
       pmanager->AddProcess(new G4eMultipleScattering(), -1, 1, 1);
       pmanager->AddProcess(new G4eIonisation,           -1, 2, 2);
       pmanager->AddProcess(new G4eBremsstrahlung(),     -1,-3, 3);
 
     } else if (particleName == "e+") {
 
       pmanager->AddProcess(new G4eMultipleScattering(), -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-"    ) {
 
       G4MuMultipleScattering* msc = new G4MuMultipleScattering();
       msc->AddEmModel(0, new G4WentzelVIModel());
       pmanager->AddProcess(msc,                     -1, 1, 1);
       pmanager->AddProcess(new G4MuIonisation,      -1, 2, 2);
       pmanager->AddProcess(new G4MuBremsstrahlung,  -1,-3, 3);
       pmanager->AddProcess(new G4MuPairProduction,  -1,-4, 4);
       pmanager->AddDiscreteProcess(new G4CoulombScattering());
 
     } else if (particleName == "alpha" ||
                particleName == "He3") {
 
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       pmanager->AddProcess(new G4ionIonisation,       -1, 2, 2);
 
     } else if (particleName == "GenericIon") {
 
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       pmanager->AddProcess(new G4ionIonisation,       -1, 2, 2);
 
     } else if (particleName == "pi+" ||
                particleName == "pi-" ||
                particleName == "kaon+" ||
                particleName == "kaon-" ||
                particleName == "proton" ) {
 
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       pmanager->AddProcess(new G4hIonisation,         -1, 2, 2);
       pmanager->AddProcess(new G4hBremsstrahlung,     -1,-3, 3);
       pmanager->AddProcess(new G4hPairProduction,     -1,-4, 4);
 
     } else if (particleName == "B+" ||
                particleName == "B-" ||
                particleName == "D+" ||
                particleName == "D-" ||
                particleName == "Ds+" ||
                particleName == "Ds-" ||
                particleName == "anti_He3" ||
                particleName == "anti_alpha" ||
                particleName == "anti_deuteron" ||
                particleName == "anti_lambda_c+" ||
                particleName == "anti_omega-" ||
                particleName == "anti_proton" ||
                particleName == "anti_sigma_c+" ||
                particleName == "anti_sigma_c++" ||
                particleName == "anti_sigma+" ||
                particleName == "anti_sigma-" ||
                particleName == "anti_triton" ||
                particleName == "anti_xi_c+" ||
                particleName == "anti_xi-" ||
                particleName == "deuteron" ||
                particleName == "lambda_c+" ||
                particleName == "omega-" ||
                particleName == "sigma_c+" ||
                particleName == "sigma_c++" ||
                particleName == "sigma+" ||
                particleName == "sigma-" ||
                particleName == "tau+" ||
                particleName == "tau-" ||
                particleName == "triton" ||
                particleName == "xi_c+" ||
                particleName == "xi-" ) {
 
       pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
       pmanager->AddProcess(new G4hIonisation,         -1, 2, 2);
     }
   }
   G4EmProcessOptions opt;
   opt.SetVerbose(verbose);
   opt.SetPolarAngleLimit(0.2);
}