BesMuonPhysics.cc

Go to the documentation of this file.

 
#include "PhySim/BesMuonPhysics.hh"
 
#include "G4ios.hh"
#include "globals.hh"
#include <iomanip>
 
BesMuonPhysics::BesMuonPhysics( const G4String& name ) : G4VPhysicsConstructor( name ) {}
 
BesMuonPhysics::~BesMuonPhysics() {}
 
#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"
 
#include "G4AntiNeutrinoMu.hh"
#include "G4AntiNeutrinoTau.hh"
#include "G4MuonMinus.hh"
#include "G4MuonPlus.hh"
#include "G4NeutrinoMu.hh"
#include "G4NeutrinoTau.hh"
#include "G4TauMinus.hh"
#include "G4TauPlus.hh"
 
void BesMuonPhysics::ConstructParticle() {
  // Mu
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();
  G4NeutrinoMu::NeutrinoMuDefinition();
  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
 
  // Tau
  G4TauMinus::TauMinusDefinition();
  G4TauPlus::TauPlusDefinition();
  G4NeutrinoTau::NeutrinoTauDefinition();
  G4AntiNeutrinoTau::AntiNeutrinoTauDefinition();
}
 
#include "G4ProcessManager.hh"
 
void BesMuonPhysics::
    ConstructProcess() { /*
                          G4ProcessManager * pManager = 0;
 
                          // Muon Plus Physics
                          pManager = G4MuonPlus::MuonPlus()->GetProcessManager();
 
                          pManager->AddProcess(&fMuPlusMultipleScattering,-1,  1, 1);
                          pManager->AddProcess(&fMuPlusIonisation,        -1,  2, 2);
                          pManager->AddProcess(&fMuPlusBremsstrahlung,    -1,  3, 3);
                          pManager->AddProcess(&fMuPlusPairProduction,    -1,  4, 4);
 
                          // Muon Minus Physics
                          pManager = G4MuonMinus::MuonMinus()->GetProcessManager();
 
                           pManager->AddProcess(&fMuMinusMultipleScattering,-1,  1, 1);
                           pManager->AddProcess(&fMuMinusIonisation,        -1,  2, 2);
                           pManager->AddProcess(&fMuMinusBremsstrahlung,    -1,  3, 3);
                           pManager->AddProcess(&fMuMinusPairProduction,    -1,  4, 4);
 
                             pManager->AddRestProcess(&fMuMinusCaptureAtRest);
 
                          // Tau Plus Physics
                          pManager = G4TauPlus::TauPlus()->GetProcessManager();
 
                          pManager->AddProcess(&fTauPlusMultipleScattering, -1, 1, 1);
                          pManager->AddProcess(&fTauPlusIonisation,         -1, 2, 2);
 
                          // Tau Minus Physics
                          pManager = G4TauMinus::TauMinus()->GetProcessManager();
 
                          pManager->AddProcess(&fTauMinusMultipleScattering, -1, 1, 1);
                          pManager->AddProcess(&fTauMinusIonisation,         -1, 2, 2);
                        */
}