From 3707b1b10ac909421a7b4527e59d85a98840639e Mon Sep 17 00:00:00 2001
From: adrien-matta <a.matta@surrey.ac.uk>
Date: Tue, 21 Oct 2014 10:49:39 +0100
Subject: [PATCH] * Fixing issue in the Physics list with proton energy loss *
 The proton were registered as ion and did not interect any mor

---
 NPSimulation/src/PhysicsList.cc | 41 +--------------------------------
 1 file changed, 1 insertion(+), 40 deletions(-)

diff --git a/NPSimulation/src/PhysicsList.cc b/NPSimulation/src/PhysicsList.cc
index cf4c2f89e..5407585c1 100644
--- a/NPSimulation/src/PhysicsList.cc
+++ b/NPSimulation/src/PhysicsList.cc
@@ -45,26 +45,6 @@
 #include "G4eBremsstrahlung.hh"
 #include "G4eplusAnnihilation.hh"
 
-// Penelope
-//#include "G4PenelopeIonisation.hh"
-//#include "G4PenelopeBremsstrahlung.hh"
-//#include "G4PenelopeAnnihilation.hh"
-
-//#include "G4PenelopeCompton.hh"
-//#include "G4PenelopeGammaConversion.hh"
-//#include "G4PenelopePhotoElectric.hh"
-//#include "G4PenelopeRayleigh.hh"
-
-// Low energy ~ Penelope
-//#include "G4LowEnergyIonisation.hh"
-//#include "G4LowEnergyBremsstrahlung.hh"
-
-//#include "G4LowEnergyCompton.hh"
-//#include "G4LowEnergyGammaConversion.hh"
-//#include "G4LowEnergyPhotoElectric.hh"
-//#include "G4LowEnergyRayleigh.hh"
-
-
 #include "G4MuIonisation.hh"
 #include "G4MuMultipleScattering.hh"
 #include "G4MuBremsstrahlung.hh"
@@ -73,7 +53,6 @@
 #include "G4hIonisation.hh"
 #include "G4ionIonisation.hh"
 #include "G4hMultipleScattering.hh"
-//#include "G4hLowEnergyIonisation.hh"
 
 #include "G4EmProcessOptions.hh"
 #include "G4ProcessManager.hh"
@@ -209,17 +188,6 @@ void PhysicsList::ConstructEM(){
       pmanager->AddDiscreteProcess(new G4PhotoElectricEffect) ;
       pmanager->AddDiscreteProcess(new G4ComptonScattering) ;
       pmanager->AddDiscreteProcess(new G4GammaConversion) ;
-      //Low energy
-      //pmanager->AddDiscreteProcess(new G4LowEnergyPhotoElectric) ;
-      //pmanager->AddDiscreteProcess(new G4LowEnergyCompton) ;
-      //pmanager->AddDiscreteProcess(new G4LowEnergyGammaConversion) ;
-      //pmanager->AddDiscreteProcess(new G4LowEnergyRayleigh) ;
-      // Penelope
-      //pmanager->AddDiscreteProcess(new G4PenelopePhotoElectric) ;
-      //pmanager->AddDiscreteProcess(new G4PenelopeCompton) ;
-      //pmanager->AddDiscreteProcess(new G4PenelopeGammaConversion) ;
-      //pmanager->AddDiscreteProcess(new G4PenelopeRayleigh) ;
-      
     }
     else if (particleName == "e-") {
       //electron
@@ -227,12 +195,6 @@ void PhysicsList::ConstructEM(){
       //standard geant4:
       pmanager->AddProcess(new G4eIonisation         , -1,  2, 2)     ;
       pmanager->AddProcess(new G4eBremsstrahlung     , -1, -1, 3)     ;
-      // Low energy:
-      //pmanager->AddProcess(new G4LowEnergyIonisation         , -1,  2, 2)     ;
-      //pmanager->AddProcess(new G4LowEnergyBremsstrahlung     , -1, -1, 3)     ;
-      // Penelope:
-      // pmanager->AddProcess(new G4PenelopeIonisation         , -1,  2, 2)     ;
-      // pmanager->AddProcess(new G4PenelopeBremsstrahlung     , -1, -1, 3)     ;
     }
     
     else if (particleName == "e+") {
@@ -257,8 +219,7 @@ void PhysicsList::ConstructEM(){
               || particleName == "He3"
               || particleName == "alpha"
               || particleName == "deuteron"
-              || particleName == "triton"
-              || particleName == "proton") {
+              || particleName == "triton") {
       pmanager->AddProcess(new G4hMultipleScattering(), -1, 1, 1) ;
       G4ionIonisation* iI = new G4ionIonisation ;
       iI->ActivateStoppingData(true) ;
-- 
GitLab