APrimePhysics.cc

Go to the documentation of this file.

#include "SimG4Core/CustomPhysics/interface/APrimePhysics.h"
#include "SimG4Core/CustomPhysics/interface/G4APrime.h"
#include "SimG4Core/CustomPhysics/interface/G4muDarkBremsstrahlung.h"
// Geant 4
#include "G4Electron.hh"
#include "G4MuonMinus.hh"
#include "G4MuonPlus.hh"
#include "G4ProcessManager.hh"
#include "G4SystemOfUnits.hh"

APrimePhysics::APrimePhysics(double APMass, const G4String& scalefile, const G4double cxBias, const G4String& name)
    : G4VPhysicsConstructor(name), aprimeDef_(nullptr) {
  apmass = APMass;
  mgfile = scalefile;
  biasFactor = cxBias;
}

APrimePhysics::~APrimePhysics() {}

void APrimePhysics::ConstructParticle() {
  aprimeDef_ = G4APrime::APrime(apmass);
  //aprimeDef->SetProcessManager(new G4ProcessManager(aprimeDef));
}

void APrimePhysics::ConstructProcess() {
  G4ParticleDefinition* muonminus = G4MuonMinus::MuonMinusDefinition();
  G4ParticleDefinition* muonplus = G4MuonPlus::MuonPlusDefinition();
  G4ProcessManager* pmplus = muonplus->GetProcessManager();
  G4ProcessManager* pmminus = muonminus->GetProcessManager();
  pmplus->AddDiscreteProcess(new G4muDarkBremsstrahlung(mgfile, biasFactor), 6);
  pmminus->AddDiscreteProcess(new G4muDarkBremsstrahlung(mgfile, biasFactor), 6);
}