Changeset 1230 for trunk/examples/novice/N06/src/ExN06PhysicsList.cc

Timestamp:

Jan 8, 2010, 3:02:48 PM (14 years ago)

Author:

garnier

Message:

update to geant4.9.3

File:

• trunk/examples/novice/N06/src/ExN06PhysicsList.cc (modified) (8 diffs)

trunk/examples/novice/N06/src/ExN06PhysicsList.cc

@@ -25 +25 @@
 //
 //
-// $Id: ExN06PhysicsList.cc,v 1.14 2007/09/30 22:51:03 gum Exp $
-// GEANT4 tag $Name:  $
+// $Id: ExN06PhysicsList.cc,v 1.17 2009/11/10 05:16:23 gum Exp $
+// GEANT4 tag $Name: geant4-09-03-cand-01 $
 //
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
@@ -53 +53 @@
 #include "G4OpRayleigh.hh"
 #include "G4OpBoundaryProcess.hh"
+
+#include "G4LossTableManager.hh"
+#include "G4EmSaturation.hh"
  
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
@@ -177 +180 @@
 #include "G4PhotoElectricEffect.hh"
 
-#include "G4MultipleScattering.hh"
+#include "G4eMultipleScattering.hh"
+#include "G4MuMultipleScattering.hh"
+#include "G4hMultipleScattering.hh"
 
 #include "G4eIonisation.hh"
@@ -209 +214 @@
     //electron
       // Construct processes for electron
-      pmanager->AddProcess(new G4MultipleScattering(),-1, 1, 1);
+      pmanager->AddProcess(new G4eMultipleScattering(),-1, 1, 1);
       pmanager->AddProcess(new G4eIonisation(),       -1, 2, 2);
       pmanager->AddProcess(new G4eBremsstrahlung(),   -1, 3, 3);
@@ -216 +221 @@
     //positron
       // Construct processes for positron
-      pmanager->AddProcess(new G4MultipleScattering(),-1, 1, 1);
+      pmanager->AddProcess(new G4eMultipleScattering(),-1, 1, 1);
       pmanager->AddProcess(new G4eIonisation(),       -1, 2, 2);
       pmanager->AddProcess(new G4eBremsstrahlung(),   -1, 3, 3);
@@ -225 +230 @@
     //muon
      // Construct processes for muon
-     pmanager->AddProcess(new G4MultipleScattering(),-1, 1, 1);
+     pmanager->AddProcess(new G4MuMultipleScattering(),-1, 1, 1);
      pmanager->AddProcess(new G4MuIonisation(),      -1, 2, 2);
      pmanager->AddProcess(new G4MuBremsstrahlung(),  -1, 3, 3);
@@ -234 +239 @@
           (particle->GetParticleName() != "chargedgeantino")) {
      // all others charged particles except geantino
-       pmanager->AddProcess(new G4MultipleScattering(),-1,1,1);
+       pmanager->AddProcess(new G4hMultipleScattering(),-1,1,1);
        pmanager->AddProcess(new G4hIonisation(),       -1,2,2);
      }
@@ -258 +263 @@
   SetVerbose(1);
   
-  theCerenkovProcess->SetMaxNumPhotonsPerStep(300);
+  theCerenkovProcess->SetMaxNumPhotonsPerStep(20);
+  theCerenkovProcess->SetMaxBetaChangePerStep(10.0);
   theCerenkovProcess->SetTrackSecondariesFirst(true);
   
   theScintillationProcess->SetScintillationYieldFactor(1.);
   theScintillationProcess->SetTrackSecondariesFirst(true);
+
+  // Use Birks Correction in the Scintillation process
+
+  G4EmSaturation* emSaturation = G4LossTableManager::Instance()->EmSaturation();
+  theScintillationProcess->AddSaturation(emSaturation);
 
   G4OpticalSurfaceModel themodel = unified;