Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm8/src/PhysicsList.cc

Timestamp:

Nov 5, 2010, 4:08:39 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/examples/extended/electromagnetic/TestEm8/src/PhysicsList.cc (modified) (8 diffs)

trunk/examples/extended/electromagnetic/TestEm8/src/PhysicsList.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-//
-// $Id: PhysicsList.cc,v 1.19 2010/06/04 17:08:15 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: PhysicsList.cc,v 1.21 2010/10/27 14:52:07 vnivanch Exp $
+// GEANT4 tag $Name: examples-V09-03-09 $
+//
+//---------------------------------------------------------------------------
+//
+// ClassName:   PhysicsList
+//
+// Description: EM physics with a possibility to add PAI model
+//
+// Author:      V.Ivanchenko 01.09.2010
+//
+//----------------------------------------------------------------------------
 //
 
@@ -39 +48 @@
 #include "G4EmStandardPhysics_option2.hh"
 #include "G4EmStandardPhysics_option3.hh"
+#include "G4EmLivermorePhysics.hh"
+#include "G4EmPenelopePhysics.hh"
 #include "G4DecayPhysics.hh"
 
@@ -47 +58 @@
 #include "G4Electron.hh"
 #include "G4Positron.hh"
+#include "G4Proton.hh"
 
 #include "G4UnitsTable.hh"
@@ -67 +79 @@
   cutForElectron  = defaultCutValue;
   cutForPositron  = defaultCutValue;
-
-  stepMaxProcess  = 0;
+  cutForProton    = defaultCutValue;
 
   pMessenger = new PhysicsListMessenger(this);
+
+  stepMaxProcess = new StepMax();
 
   // Decay Physics is always defined
@@ -140 +153 @@
     emPhysicsList = new G4EmStandardPhysics_option3();
 
+  } else if (name == "emlivermore") {
+
+    emName = name;
+    delete emPhysicsList;
+    emPhysicsList = new G4EmLivermorePhysics();
+
+  } else if (name == "empenelope") {
+
+    emName = name;
+    delete emPhysicsList;
+    emPhysicsList = new G4EmPenelopePhysics();
+
   } else if (name == "pai") {
 
@@ -163 +188 @@
 {
   // Step limitation seen as a process
-  stepMaxProcess = new StepMax();
 
   theParticleIterator->reset();
@@ -195 +219 @@
   SetCutValue(cutForElectron, "e-");
   SetCutValue(cutForPositron, "e+");
+  SetCutValue(cutForProton, "proton");
 
   if ( verboseLevel > 0 ) { DumpCutValuesTable(); }
@@ -221 +246 @@
   cutForPositron = cut;
   SetParticleCuts(cutForPositron, G4Positron::Positron());
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
+
+void PhysicsList::SetCutForProton(G4double cut)
+{
+  cutForPositron = cut;
+  SetParticleCuts(cutForProton, G4Proton::Proton());
 }