Changeset 1313 for trunk/examples/advanced/hadrontherapy/src/HadrontherapyPhysicsList.cc

Timestamp:

Jun 14, 2010, 3:54:58 PM (14 years ago)

Author:

garnier

Message:

geant4.9.4 beta rc0

File:

• trunk/examples/advanced/hadrontherapy/src/HadrontherapyPhysicsList.cc (modified) (7 diffs)

trunk/examples/advanced/hadrontherapy/src/HadrontherapyPhysicsList.cc

@@ -76 +76 @@
 //                        and local physic for ion-ion enelastic processes)
 
+#include "G4RunManager.hh"
+#include "G4Region.hh"
+#include "G4RegionStore.hh"
 #include "HadrontherapyPhysicsList.hh"
 #include "HadrontherapyPhysicsListMessenger.hh"
@@ -86 +89 @@
 #include "LocalINCLIonIonInelasticPhysic.hh"         // Physic dedicated to the ion-ion inelastic processes using INCL/ABLA
 
-// Physic lists (contained inside the Geant4 distribution)
+// Physic lists (contained inside the Geant4 source code, in the 'physicslists folder')
 #include "G4EmStandardPhysics_option3.hh"
 #include "G4EmLivermorePhysics.hh"
@@ -173 +176 @@
 {
   // transportation
-  //
   AddTransportation();
 
   // electromagnetic physics list
-  //
   emPhysicsList->ConstructProcess();
   em_config.AddModels();
 
   // decay physics list
-  //
   decPhysicsList->ConstructProcess();
 
@@ -207 +207 @@
   //   ELECTROMAGNETIC MODELS
   /////////////////////////////////////////////////////////////////////////////
-
-  if (name == "standard_opt3") {
+    if (name == "standard_opt3") {
     emName = name;
     delete emPhysicsList;
     emPhysicsList = new G4EmStandardPhysics_option3();
+    G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
     G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmStandardPhysics_option3" << G4endl;
 
@@ -218 +218 @@
     delete emPhysicsList;
     emPhysicsList = new G4EmLivermorePhysics();
+    G4RunManager::GetRunManager()-> PhysicsHasBeenModified();
     G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmLivermorePhysics" << G4endl;
 
@@ -308 +309 @@
   SetCutValue(cutForPositron, "e+");
 
+  // Set cuts for detector
+  SetDetectorCut(defaultCutValue); 
   if (verboseLevel>0) DumpCutValuesTable();
 }
@@ -331 +334 @@
   SetParticleCuts(cutForPositron, G4Positron::Positron());
 }
+
+void HadrontherapyPhysicsList::SetDetectorCut(G4double cut)
+{
+
+  G4String regionName = "DetectorLog";
+  G4Region* region = G4RegionStore::GetInstance()->GetRegion(regionName);
+
+  G4ProductionCuts* cuts = new G4ProductionCuts ;
+  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("gamma"));
+  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e-"));
+  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e+"));
+  region -> SetProductionCuts(cuts);
+}
+