ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4KaonPlus.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4KaonPlus.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 //
26 //
27 //
28 //
29 // ----------------------------------------------------------------------
30 // GEANT 4 class implementation file
31 //
32 // History: first implementation, based on object model of
33 // 4th April 1996, G.Cosmo
34 // H.Kurashige 7 Jul 96
35 // **********************************************************************
36 // New impelemenataion as an utility class M.Asai, 26 July 2004
37 // ----------------------------------------------------------------------
38 
39 #include "G4KaonPlus.hh"
40 #include "G4SystemOfUnits.hh"
41 #include "G4ParticleTable.hh"
42 
44 #include "G4KL3DecayChannel.hh"
45 #include "G4DecayTable.hh"
46 
47 // ######################################################################
48 // ### KAONPLUS ###
49 // ######################################################################
50 
52 
54 {
55  if (theInstance !=0) return theInstance;
56  const G4String name = "kaon+";
57  // search in particle table]
59  G4ParticleDefinition* anInstance = pTable->FindParticle(name);
60  if (anInstance ==0)
61  {
62  // create particle
63  //
64  // Arguments for constructor are as follows
65  // name mass width charge
66  // 2*spin parity C-conjugation
67  // 2*Isospin 2*Isospin3 G-parity
68  // type lepton number baryon number PDG encoding
69  // stable lifetime decay table
70  // shortlived subType anti_encoding
71 
72  anInstance = new G4ParticleDefinition(
73  name, 0.493677*GeV, 5.317e-14*MeV, +1.*eplus,
74  0, -1, 0,
75  1, +1, 0,
76  "meson", 0, 0, 321,
77  false, 12.380*ns, NULL,
78  false, "kaon");
79 
80  //create Decay Table
81  G4DecayTable* table = new G4DecayTable();
82 
83  // create decay channels
84  G4VDecayChannel** mode = new G4VDecayChannel*[6];
85  // kaon+ -> mu+ + nu_mu
86  mode[0] = new G4PhaseSpaceDecayChannel("kaon+",0.6355,2,"mu+","nu_mu");
87  // kaon+ -> pi+ + pi0
88  mode[1] = new G4PhaseSpaceDecayChannel("kaon+",0.2066,2,"pi+","pi0");
89  // kaon+ -> pi+ + pi+ + pi-
90  mode[2] = new G4PhaseSpaceDecayChannel("kaon+",0.0559,3,"pi+","pi+","pi-");
91  // kaon+ -> pi+ + pi0 + pi0
92  mode[3] = new G4PhaseSpaceDecayChannel("kaon+",0.01761,3,"pi+","pi0","pi0");
93  // kaon+ -> pi0 + e+ + nu_e (Ke3)
94  mode[4] = new G4KL3DecayChannel("kaon+",0.0507,"pi0","e+","nu_e");
95  // kaon+ -> pi0 + mu+ + nu_mu (Kmu3)
96  mode[5] = new G4KL3DecayChannel("kaon+",0.0335,"pi0","mu+","nu_mu");
97 
98  for (G4int index=0; index <6; index++ ) table->Insert(mode[index]);
99  delete [] mode;
100 
101  anInstance->SetDecayTable(table);
102  }
103  theInstance = reinterpret_cast<G4KaonPlus*>(anInstance);
104  return theInstance;
105 }
106 
108 {
109  return Definition();
110 }
111 
113 {
114  return Definition();
115 }
116