G4XTankPhysicsListCustomization.cc

Go to the documentation of this file.

#include "G4XTankPhysicsListCustomization.h"
#include <G4UImanager.hh>
#include "G4XTankSimulator.h"

#include "G4XTankOpBoundaryProcess.h"
#include "G4XTankFastCerenkov.h"

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

using namespace G4XTankSimulatorAG;

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

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


G4XTankPhysicsListCustomization::~G4XTankPhysicsListCustomization()
{
}


void G4XTankPhysicsListCustomization::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 G4XTankPhysicsListCustomization::ActivateCustomProcesses(){
  G4cerr << "G4XTankPhysicsListCustomization::ActivateCustomProcesses" << G4endl;

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

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

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

  return true;
}

bool G4XTankPhysicsListCustomization::InactivateCustomProcesses(){
  G4cerr << "G4XTankPhysicsListCustomization::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
G4XTankPhysicsListCustomization::ConstructProcess()
{
  const G4int maxNumPhotons = 3;

  G4XTankSimulatorAG::G4XTankFastCerenkov* const fastCerenkov = new G4XTankSimulatorAG::G4XTankFastCerenkov("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();

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