ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4Monopole.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4Monopole.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
28 //
29 //
30 //---------------------------------------------------------------------------
31 //
32 // ClassName: G4Monopole
33 //
34 // Description:
35 //
36 // Authors: 21.03.05 V.Ivanchenko
37 //
38 // Modified:
39 //
40 // 12.07.10 S.Burdin (changed the magnetic and electric charge variables from integer to double)
41 //----------------------------------------------------------------------------
42 //
43 
44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46 
47 #include "G4Monopole.hh"
48 #include "G4ParticleTable.hh"
49 #include "G4PhysicalConstants.hh"
50 #include "G4SystemOfUnits.hh"
51 
52 // ######################################################################
53 // ### Monopole ###
54 // ######################################################################
55 
56 G4Monopole* G4Monopole::theMonopole = 0;
57 G4double G4Monopole::magCharge = 0.0;
58 
59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
60 
61 G4Monopole::G4Monopole(
62  const G4String& aName, G4double mass,
63  G4double width, G4double charge,
64  G4int iSpin, G4int iParity,
65  G4int iConjugation, G4int iIsospin,
66  G4int iIsospin3, G4int gParity,
67  const G4String& pType, G4int lepton,
68  G4int baryon, G4int encoding,
69  G4bool stable, G4double lifetime,
70  G4DecayTable *decaytable)
71  : G4ParticleDefinition( aName, mass, width, charge, iSpin, iParity,
72  iConjugation, iIsospin, iIsospin3, gParity, pType,
73  lepton, baryon, encoding, stable, lifetime, decaytable )
74 {}
75 
76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
77 
78 G4Monopole::~G4Monopole()
79 {}
80 
81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
82 //
83 // Arguments for constructor are as follows
84 // name mass width charge
85 // 2*spin parity C-conjugation
86 // 2*Isospin 2*Isospin3 G-parity
87 // type lepton number baryon number PDG encoding
88 // stable lifetime decay table
89 //
90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
91 
92 G4Monopole* G4Monopole::MonopoleDefinition(G4double mass, G4double mCharge,
93  G4double eCharge)
94 {
95  if(!theMonopole) {
96  magCharge = eplus * mCharge / fine_structure_const * 0.5;
97  theMonopole = new G4Monopole(
98  "monopole", mass, 0.0*MeV, eplus*eCharge,
99  0, 0, 0,
100  0, 0, 0,
101  "boson", 0, 0, 0,
102  true, -1.0, 0);
103 
104 
105  G4cout << "Monopole is created: m(GeV)= " << theMonopole->GetPDGMass()/GeV
106  << " Qel= " << theMonopole->GetPDGCharge()/eplus
107  << " Qmag= " << magCharge/eplus
108  << G4endl;
109  }
110  return theMonopole;
111 }
112 
113 
114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
115 
116 G4Monopole* G4Monopole::Monopole()
117 {
118  if(!theMonopole) { theMonopole = MonopoleDefinition(); }
119  return theMonopole;
120 }
121 
122 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
123 
124 G4double G4Monopole::MagneticCharge() const
125 {
126  return magCharge;
127 }
128 
129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
130