Allow enabling DNA physics list for specific regions

mjordanst · May 25, 2023 · bafa8cd · bafa8cd
1 parent c64675e
commit bafa8cd
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 3 deletions.
diff --git a/source/physics/include/GatePhysicsList.hh b/source/physics/include/GatePhysicsList.hh
@@ -26,6 +26,8 @@
 #include "GateVProcess.hh"
 #include "GateUserLimits.hh"
 
+class G4EmDNAPhysicsActivator;
+
 //class GateVProcess;
 class GatePhysicsListMessenger;
 class GatePhysicsList: public G4VModularPhysicsList
@@ -140,6 +142,8 @@ protected:
   G4double mLowEnergyRangeLimit;
 
   G4EmParameters *emPar;
+
+	G4EmDNAPhysicsActivator* emDNAActivator;
 };
 
 

diff --git a/source/physics/src/GatePhysicsList.cc b/source/physics/src/GatePhysicsList.cc
@@ -55,6 +55,7 @@
 #include "G4EmDNAPhysics_option6.hh"
 #include "G4EmDNAPhysics_option7.hh"
 #include "G4EmDNAPhysics_option8.hh"
+#include "G4EmDNAPhysicsActivator.hh"
 #include "G4LossTableManager.hh"
 #include "G4UAtomicDeexcitation.hh"
 #include "G4RadioactiveDecayPhysics.hh"
@@ -111,6 +112,10 @@ GatePhysicsList::GatePhysicsList(): G4VModularPhysicsList()
 #if G4VERSION_MAJOR >= 10 && G4VERSION_MINOR >= 5
   mUseICRU90Data = false;
 #endif
+
+	// used to conditionally activate DNA physics list in specific regions
+	// while using other em physics list elsewhere
+	emDNAActivator = new G4EmDNAPhysicsActivator;
 }
 //-----------------------------------------------------------------------------------------
 
@@ -185,6 +190,8 @@ GatePhysicsList::~GatePhysicsList()
     if((*vect)[i]) G4cout<<"Process= "<<(*vect)[i]->GetProcessName()<< Gateendl;
     }
     }*/
+
+	delete emDNAActivator;
 }
 //-----------------------------------------------------------------------------------------
 
@@ -285,6 +292,9 @@ void GatePhysicsList::ConstructProcess()
      }
   */
 
+	// Construct G4 DNA activator process
+	emDNAActivator->ConstructProcess();
+
   mLoadState++;
 }
 //-----------------------------------------------------------------------------------------
@@ -467,11 +477,8 @@ void GatePhysicsList::ConstructParticle()
   //  Construct hybridino
   G4Hybridino::HybridinoDefinition();
 
-
-
   //Construct G4DNA particles
 
-
   G4DNAGenericIonsManager* dnagenericIonsManager;
   dnagenericIonsManager=G4DNAGenericIonsManager::Instance();
   dnagenericIonsManager->GetIon("hydrogen");
@@ -483,6 +490,9 @@ void GatePhysicsList::ConstructParticle()
   dnagenericIonsManager->GetIon("iron");
   dnagenericIonsManager->GetIon("oxygen");
 
+	// Construct G4 DNA activator particles
+	emDNAActivator->ConstructParticle();
+
  //Construct positroniums
  GateParaPositronium::ParaPositroniumDefinition();
  GateOrthoPositronium::OrthoPositroniumDefinition();