Tue 20 Nov 17:26:02 GMT 2018

1 files changed, 0 insertions, 62 deletions

diff --git a/src/G4BeamTestEMPhysics.cxx b/src/G4BeamTestEMPhysics.cxx
index 1aa9a7c..55ecf99 100644
--- a/src/G4BeamTestEMPhysics.cxx
+++ b/src/G4BeamTestEMPhysics.cxx
@@ -55,66 +55,4 @@ void G4BeamTestEMPhysics::ConstructProcess()
   pManager->AddProcess(&positronIonisation,         -1, 2, 2);
   pManager->AddProcess(&positronBremsStrahlung,     -1, 3, 3);  
   pManager->AddProcess(&annihilation,                0,-1, 4);  
-
-  // Cerenkov Physics
-  ConstructOp();
-}
-
-void G4BeamTestEMPhysics::ConstructOp()
-{
-  G4Cerenkov* cerenkovProcess = &cerenkov;
-  G4int fMaxNumPhotonStep = 20;
-
-  cerenkovProcess->SetMaxNumPhotonsPerStep(fMaxNumPhotonStep);
-  G4cout << "MaxNumPhoton" << fMaxNumPhotonStep << G4endl; 
-
-  cerenkovProcess->SetMaxBetaChangePerStep(10.0);
-  cerenkovProcess->SetTrackSecondariesFirst(true);
- G4Scintillation* scintillationProcess = new G4Scintillation("Scintillation");
- scintillationProcess->SetScintillationYieldFactor(1.);
- scintillationProcess->SetTrackSecondariesFirst(true);
-  G4OpAbsorption* absorptionProcess = new G4OpAbsorption();
-  G4OpRayleigh* rayleighScatteringProcess = new G4OpRayleigh();
-  G4OpMieHG* mieHGScatteringProcess = new G4OpMieHG();
-  G4OpBoundaryProcess* boundaryProcess = new G4OpBoundaryProcess();
-
-/*   cerenkovProcess->SetVerboseLevel(fVerboseLebel); */
-/*  // scintillationProcess->SetVerboseLevel(fVerboseLebel); */
-/*   absorptionProcess->SetVerboseLevel(fVerboseLebel); */
-/*   rayleighScatteringProcess->SetVerboseLevel(fVerboseLebel); */
-/*   mieHGScatteringProcess->SetVerboseLevel(fVerboseLebel); */
-/*   boundaryProcess->SetVerboseLevel(fVerboseLebel); */
-  
-  // Use Birks Correction in the Scintillation process
-
- if(!G4Threading::IsWorkerThread())
- {
-   G4EmSaturation* emSaturation =
-             G4LossTableManager::Instance()->EmSaturation();
-     scintillationProcess->AddSaturation(emSaturation);
- }
-
-  auto theParticleIterator = GetParticleIterator();
-  theParticleIterator->reset();
-  while( (*theParticleIterator)() ){
-    G4ParticleDefinition* particle = theParticleIterator->value();
-    G4ProcessManager* pmanager = particle->GetProcessManager();
-    G4String particleName = particle->GetParticleName();
-    if (cerenkovProcess->IsApplicable(*particle)) {
-      pmanager->AddProcess(cerenkovProcess);
-      pmanager->SetProcessOrdering(cerenkovProcess,idxPostStep);
-    }
-   if (scintillationProcess->IsApplicable(*particle)) {
-     pmanager->AddProcess(scintillationProcess);
-     pmanager->SetProcessOrderingToLast(scintillationProcess, idxAtRest);
-     pmanager->SetProcessOrderingToLast(scintillationProcess, idxPostStep);
-   }
-    if (particleName == "opticalphoton") {
-      G4cout << " AddDiscreteProcess to OpticalPhoton " << G4endl;
-      pmanager->AddDiscreteProcess(absorptionProcess);
-      pmanager->AddDiscreteProcess(rayleighScatteringProcess);
-      pmanager->AddDiscreteProcess(mieHGScatteringProcess);
-      pmanager->AddDiscreteProcess(boundaryProcess);
-    }
-  }
 }