1 #include "SimG4Core/Physics/interface/G4Monopole.hh"
2 #include "G4ParticleTable.hh"
4 G4Monopole::G4Monopole(
const G4String aName, G4int pdgEncoding, G4double mass,
5 G4int mCharge, G4int eCharge) :
6 G4ParticleDefinition(aName, mass, 0.0*MeV, eplus*eCharge, 0, 0, 0,
7 0, 0, 0,
"boson", 0, 0, pdgEncoding,
true, -1.0, 0) {
9 magCharge = eplus * G4double(mCharge) / fine_structure_const * 0.5;
10 G4cout <<
"Monopole is created: m(GeV)= " << GetPDGMass()/GeV
11 <<
" Qel= " << GetPDGCharge()/eplus
12 <<
" Qmag= " << magCharge/eplus
13 <<
" PDG encoding = " << pdgEncoding << G4endl;
16 G4Monopole::~G4Monopole() {}