G4TankSimulatorOG/G4TankPhysicsListCustomization.cc

Go to the documentation of this file.

#include "G4TankPhysicsListCustomization.h"
#include <G4UImanager.hh>
#include "G4TankSimulator.h"

#include "G4TankOpBoundaryProcess.h"
#include "G4TankFastCerenkov.h"

#include <G4MuonMinusCaptureAtRest.hh>
#include <G4ParticleDefinition.hh>
#include <G4ParticleTypes.hh>
#include <G4ProcessManager.hh>

using namespace G4TankSimulatorOG;

G4TankPhysicsListCustomization::G4TankPhysicsListCustomization(bool fastCerenkov):
  tls::G4VPhysicsListCustomization()
{
  defaultCutValue = 0.1*mm;

  fFastCerenkov = fastCerenkov;
  fUImanager = G4UImanager::GetUIpointer();
}


G4TankPhysicsListCustomization::~G4TankPhysicsListCustomization()
{
}


void G4TankPhysicsListCustomization::SetCustomCuts(){
  const double cutForGamma    = 30*mm;
  const double cutForElectron = 0.5*mm;
  const double cutForPositron = 1e-5*mm;

  SetCutValue(cutForGamma, "gamma");
  SetCutValue(cutForElectron, "e-");
  SetCutValue(cutForPositron, "e+");
  SetCutValue(0, "proton"); //why??
}

bool G4TankPhysicsListCustomization::ActivateCustomProcesses(){
  G4cerr << "G4TankPhysicsListCustomization::ActivateCustomProcesses" << G4endl;

  if(fFastCerenkov) {
    G4cerr << "activating TankFastCerenkov" << G4endl;
    fUImanager->ApplyCommand("/process/activate TankFastCerenkov");
    G4cerr << "de-activating StandardCerenkov" << G4endl;
    fUImanager->ApplyCommand("/process/inactivate StandardCerenkov");
  }

  if(G4TankSimulator::MuCaptureOn()){
    G4cerr << "activating muMinusCaptureAtRest" << G4endl;
    fUImanager->ApplyCommand("/process/activate muMinusCaptureAtRest");
  }

  G4cerr << "activating TankOpBoundaryProcess" << G4endl;
  fUImanager->ApplyCommand("/process/activate TankOpBoundaryProcess");

  return true;
}

bool G4TankPhysicsListCustomization::InactivateCustomProcesses(){
  G4cerr << "G4TankPhysicsListCustomization::InactivateCustomProcesses" << G4endl;

  G4cerr << "de-activating TankFastCerenkov" << G4endl;
  fUImanager->ApplyCommand("/process/inactivate TankFastCerenkov");

  if(fFastCerenkov){
    G4cerr << "activating StandardCerenkov" << G4endl;
    fUImanager->ApplyCommand("/process/activate StandardCerenkov");
  }

  G4cerr << "de-activating muMinusCaptureAtRest" << G4endl;
  fUImanager->ApplyCommand("/process/inactivate muMinusCaptureAtRest");

  G4cerr << "de-activating TankOpBoundaryProcess" << G4endl;
  fUImanager->ApplyCommand("/process/inactivate TankOpBoundaryProcess");

  return true;
}


void
G4TankPhysicsListCustomization::ConstructProcess()
{
  const G4int maxNumPhotons = 3;

  G4TankSimulatorOG::G4TankFastCerenkov* const fastCerenkov = new G4TankSimulatorOG::G4TankFastCerenkov("TankFastCerenkov");
  fastCerenkov->SetTrackSecondariesFirst(true);
  fastCerenkov->SetMaxNumPhotonsPerStep(maxNumPhotons);
  G4Electron::Electron()->GetProcessManager()->AddContinuousProcess(fastCerenkov);
  G4Positron::Positron()->GetProcessManager()->AddContinuousProcess(fastCerenkov);
  G4MuonMinus::MuonMinusDefinition()->GetProcessManager()->AddContinuousProcess(fastCerenkov);
  G4MuonPlus::MuonPlusDefinition()->GetProcessManager()->AddContinuousProcess(fastCerenkov);

  G4MuonMinus::MuonMinusDefinition()->GetProcessManager()->AddProcess(new G4MuonMinusCaptureAtRest);

  G4ProcessManager& pManager =
    *G4OpticalPhoton::OpticalPhotonDefinition()->GetProcessManager();

  G4TankSimulatorOG::G4TankOpBoundaryProcess* const boundary = new G4TankSimulatorOG::G4TankOpBoundaryProcess("TankOpBoundaryProcess");
  const G4OpticalSurfaceModel model = unified;
  boundary->SetModel(model);
  pManager.AddDiscreteProcess(boundary);
}