Modifying EventGeneratorisotropic.cc to simulate other nucleus as a

source

Inputs/EventGenerator/proton.source

@@ -4,10 +4,10 @@
 %			   Energy are given in MeV , Position in mm				  %	
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Isotropic
-	EnergyLow=  80	
-	EnergyHigh= 80
+	EnergyLow=  20	
+	EnergyHigh= 20
 	HalfOpenAngleMin= 0
-	HalfOpenAngleMax= 90
+	HalfOpenAngleMax= 5
 	x0= 0	
 	y0= 0	
 	z0= 0	

NPSimulation/Core/EventGeneratorIsotropic.cc

@@ -165,6 +165,7 @@ void EventGeneratorIsotropic::ReadConfiguration(string Path,int){
         else if(m_particleName=="alpha") { m_particleName="4He" ; check_ExcitationEnergy = true ;}
         else if(m_particleName=="gamma") { check_ExcitationEnergy = true ;}
         else if(m_particleName=="neutron") { check_ExcitationEnergy = true ;}
+          else { check_ExcitationEnergy = true ;}
 
       }
       

NPSimulation/Core/MaterialManager.cc

@@ -123,6 +123,7 @@ G4Material* MaterialManager::GetMaterialFromLibrary(string Name){
       m_Material[Name]=material;
       return material; 
     }
+    
 
     else  if(Name == "Kapton"){
       G4Material* material = new G4Material(Name, 1.39*g/cm3,3);
@@ -317,7 +318,7 @@ G4Material* MaterialManager::GetMaterialFromLibrary(string Name){
       return material; 
     }
       
-    else  if(Name == "P10"){
+    else  if(Name == "P10_1atm"){
         G4Material* material = new G4Material(Name, 1.74*mg/cm3,3); //@ 0K, 1 atm
         material->AddElement(GetElementFromLibrary("Ar"),0.9222);
         material->AddElement(GetElementFromLibrary("C"),0.0623);
@@ -325,7 +326,29 @@ G4Material* MaterialManager::GetMaterialFromLibrary(string Name){
         m_Material[Name]=material;
         return material;
     }
+      
+    else  if(Name == "P10"){
+        G4Material* material = new G4Material(Name, 0.57*mg/cm3,3); //@ 0K, 1/3 atm
+        material->AddElement(GetElementFromLibrary("Ar"),0.9222);
+        material->AddElement(GetElementFromLibrary("C"),0.0623);
+        material->AddElement(GetElementFromLibrary("H"),0.0155);
+        m_Material[Name]=material;
+        return material;
+    }
 
+       else  if(Name == "Air_1atm"){ // 1 atm
+           G4Material* material = new G4Material("Air", 1.290*mg/cm3, 2);
+           material->AddElement(GetElementFromLibrary("N"), 0.7);
+           material->AddElement(GetElementFromLibrary("O"), 0.3);
+       }
+      
+       else  if(Name == "Air"){ // 1/3 atm
+           G4Material* material = new G4Material("Air", 1.290/3*mg/cm3, 2);
+           material->AddElement(GetElementFromLibrary("N"), 0.7);
+           material->AddElement(GetElementFromLibrary("O"), 0.3);
+       }
+      
+      
     else{
       G4cout << "ERROR: Material requested \""<< Name <<"\" is not available in the Material Librairy" << G4endl;
       exit(1);