source: trunk/examples/extended/exoticphysics/monopole/src/G4Monopole.cc@ 1036

Last change on this file since 1036 was 807, checked in by garnier, 17 years ago

update

File size: 4.8 KB
RevLine 
[807]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//
26//---------------------------------------------------------------------------
27//
28// ClassName: G4Monopole
29//
30// Description:
31//
32// Authors: 21.03.05 V.Ivanchenko
33//
34// Modified:
35//
36//----------------------------------------------------------------------------
37//
38
39//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
40//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
41
42#include "G4Monopole.hh"
43#include "G4ParticleTable.hh"
44
45// ######################################################################
46// ### Monopole ###
47// ######################################################################
48
49G4Monopole* G4Monopole::theMonopole = 0;
50G4double G4Monopole::magCharge = 0.0;
51
52G4Monopole::G4Monopole(
53 const G4String& aName, G4double mass,
54 G4double width, G4double charge,
55 G4int iSpin, G4int iParity,
56 G4int iConjugation, G4int iIsospin,
57 G4int iIsospin3, G4int gParity,
58 const G4String& pType, G4int lepton,
59 G4int baryon, G4int encoding,
60 G4bool stable, G4double lifetime,
61 G4DecayTable *decaytable)
62 : G4ParticleDefinition( aName, mass, width, charge, iSpin, iParity,
63 iConjugation, iIsospin, iIsospin3, gParity, pType,
64 lepton, baryon, encoding, stable, lifetime, decaytable )
65{}
66
67//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
68
69G4Monopole::~G4Monopole()
70{}
71
72//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
73//
74// Arguments for constructor are as follows
75// name mass width charge
76// 2*spin parity C-conjugation
77// 2*Isospin 2*Isospin3 G-parity
78// type lepton number baryon number PDG encoding
79// stable lifetime decay table
80//
81//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
82
83G4Monopole* G4Monopole::MonopoleDefinition(G4double mass, G4int mCharge, G4int eCharge)
84{
85 if(!theMonopole) {
86 magCharge = eplus * G4double(mCharge) / fine_structure_const * 0.5;
87 theMonopole = new G4Monopole(
88 "monopole", mass, 0.0*MeV, eplus*G4double(eCharge),
89 0, 0, 0,
90 0, 0, 0,
91 "boson", 0, 0, 0,
92 true, -1.0, 0);
93
94
95 G4cout << "Monopole is created: m(GeV)= " << theMonopole->GetPDGMass()/GeV
96 << " Qel= " << theMonopole->GetPDGCharge()/eplus
97 << " Qmag= " << magCharge/eplus
98 << G4endl;
99 }
100 return theMonopole;
101}
102
103
104//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
105G4Monopole* G4Monopole::Monopole()
106{
107 if(!theMonopole) theMonopole = MonopoleDefinition();
108 return theMonopole;
109}
110
111//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
112
Note: See TracBrowser for help on using the repository browser.