Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // >> 26 // >> 27 // $Id: G4TauMinus.cc 95951 2016-03-03 10:44:16Z gcosmo $ >> 28 // >> 29 // 26 // ------------------------------------------- 30 // ---------------------------------------------------------------------- 27 // GEANT 4 class implementation file 31 // GEANT 4 class implementation file 28 // 32 // 29 // History: first implementation, based o 33 // History: first implementation, based on object model of 30 // 4th April 1996, G.Cosmo 34 // 4th April 1996, G.Cosmo 31 // 7 July 1996 H.Kurash << 35 // 7 July 1996 H.Kurashige 32 // ******************************************* 36 // ********************************************************************** 33 // New impelemenataion as an utility class M 37 // New impelemenataion as an utility class M.Asai, 26 July 2004 34 // ------------------------------------------- 38 // ---------------------------------------------------------------------- 35 39 36 #include "G4TauMinus.hh" 40 #include "G4TauMinus.hh" 37 << 38 #include "G4DecayTable.hh" << 39 #include "G4ParticleTable.hh" << 40 #include "G4PhaseSpaceDecayChannel.hh" << 41 #include "G4PhysicalConstants.hh" 41 #include "G4PhysicalConstants.hh" 42 #include "G4String.hh" << 43 #include "G4SystemOfUnits.hh" 42 #include "G4SystemOfUnits.hh" >> 43 #include "G4ParticleTable.hh" >> 44 >> 45 #include "G4PhaseSpaceDecayChannel.hh" 44 #include "G4TauLeptonicDecayChannel.hh" 46 #include "G4TauLeptonicDecayChannel.hh" 45 #include "G4Types.hh" << 47 #include "G4DecayTable.hh" 46 #include "G4VDecayChannel.hh" << 47 48 48 G4TauMinus* G4TauMinus::theInstance = nullptr; << 49 // ###################################################################### >> 50 // ### TAUMINUS ### >> 51 // ###################################################################### >> 52 G4TauMinus* G4TauMinus::theInstance = 0; 49 53 50 G4TauMinus* G4TauMinus::Definition() 54 G4TauMinus* G4TauMinus::Definition() 51 { 55 { 52 if (theInstance != nullptr) return theInstan << 56 if (theInstance !=0) return theInstance; 53 const G4String name = "tau-"; 57 const G4String name = "tau-"; 54 // search in particle table] 58 // search in particle table] 55 G4ParticleTable* pTable = G4ParticleTable::G 59 G4ParticleTable* pTable = G4ParticleTable::GetParticleTable(); 56 G4ParticleDefinition* anInstance = pTable->F 60 G4ParticleDefinition* anInstance = pTable->FindParticle(name); 57 if (anInstance == nullptr) { << 61 if (anInstance ==0) 58 // create particle << 62 { 59 // << 63 // create particle 60 // Arguments for constructor are as fol << 64 // 61 // name mass << 65 // Arguments for constructor are as follows 62 // 2*spin parity C- << 66 // name mass width charge 63 // 2*Isospin 2*Isospin3 << 67 // 2*spin parity C-conjugation 64 // type lepton number ba << 68 // 2*Isospin 2*Isospin3 G-parity 65 // stable lifetime << 69 // type lepton number baryon number PDG encoding 66 // shortlived subType << 70 // stable lifetime decay table 67 << 71 // shortlived subType anti_encoding 68 // clang-format off << 69 anInstance = new G4ParticleDefinition( 72 anInstance = new G4ParticleDefinition( 70 name, 1.77686*GeV, 2.267 73 name, 1.77686*GeV, 2.267e-9*MeV, -1.*eplus, 71 1, 0, 0, 74 1, 0, 0, 72 0, 0, 0, 75 0, 0, 0, 73 "lepton", 1, 76 "lepton", 1, 0, 15, 74 false, 290.3e-6*ns, nullptr, << 77 false, 290.3e-6*ns, NULL, 75 false, "tau" 78 false, "tau" 76 ); 79 ); 77 // clang-format on << 78 << 79 // Bohr Magnetron << 80 G4double muB = -0.5 * eplus * hbar_Planck << 81 80 82 anInstance->SetPDGMagneticMoment(muB * 1.0 << 81 // Bohr Magnetron >> 82 G4double muB = -0.5*eplus*hbar_Planck/(anInstance->GetPDGMass()/c_squared) ; >> 83 >> 84 anInstance->SetPDGMagneticMoment( muB * 1.00118); >> 85 >> 86 //create Decay Table >> 87 G4DecayTable* table = new G4DecayTable(); >> 88 >> 89 // create decay channels >> 90 G4VDecayChannel* mode; >> 91 // tau- -> mu- + anti_nu_mu + nu_tau >> 92 mode = new G4TauLeptonicDecayChannel("tau-",0.1736,"mu-"); >> 93 table->Insert(mode); >> 94 // tau- -> e- + anti_nu_e + nu_tau >> 95 mode = new G4TauLeptonicDecayChannel("tau-",0.1784,"e-"); >> 96 table->Insert(mode); >> 97 // tau- -> pi- + nu_tau >> 98 mode = new G4PhaseSpaceDecayChannel("tau-",0.1106,2,"pi-","nu_tau"); >> 99 table->Insert(mode); >> 100 // tau- -> pi0 + pi- + nu_tau >> 101 mode = new G4PhaseSpaceDecayChannel("tau-",0.2541,3,"pi0","pi-","nu_tau"); >> 102 table->Insert(mode); >> 103 // tau- -> pi0 + pi0 + pi- + nu_tau >> 104 mode = new G4PhaseSpaceDecayChannel(); >> 105 mode->SetParent("tau-"); >> 106 mode->SetBR(0.0917); >> 107 mode->SetNumberOfDaughters(4); >> 108 mode->SetDaughter(0,"pi0"); >> 109 mode->SetDaughter(1,"pi0"); >> 110 mode->SetDaughter(2,"pi-"); >> 111 mode->SetDaughter(3,"nu_tau"); >> 112 table->Insert(mode); >> 113 // tau- -> pi- + pi- + pi+ + nu_tau >> 114 mode = new G4PhaseSpaceDecayChannel(); >> 115 mode->SetParent("tau-"); >> 116 mode->SetBR(0.0946); >> 117 mode->SetNumberOfDaughters(4); >> 118 mode->SetDaughter(0,"pi-"); >> 119 mode->SetDaughter(1,"pi-"); >> 120 mode->SetDaughter(2,"pi+"); >> 121 mode->SetDaughter(3,"nu_tau"); >> 122 table->Insert(mode); 83 123 84 // create Decay Table << 124 anInstance->SetDecayTable(table); 85 auto table = new G4DecayTable(); << 86 << 87 // create decay channels << 88 G4VDecayChannel* mode; << 89 // tau- -> mu- + anti_nu_mu + nu_tau << 90 mode = new G4TauLeptonicDecayChannel("tau- << 91 table->Insert(mode); << 92 // tau- -> e- + anti_nu_e + nu_tau << 93 mode = new G4TauLeptonicDecayChannel("tau- << 94 table->Insert(mode); << 95 // tau- -> pi- + nu_tau << 96 mode = new G4PhaseSpaceDecayChannel("tau-" << 97 table->Insert(mode); << 98 // tau- -> pi0 + pi- + nu_tau << 99 mode = new G4PhaseSpaceDecayChannel("tau-" << 100 table->Insert(mode); << 101 // tau- -> pi0 + pi0 + pi- + nu_tau << 102 mode = new G4PhaseSpaceDecayChannel(); << 103 mode->SetParent("tau-"); << 104 mode->SetBR(0.0917); << 105 mode->SetNumberOfDaughters(4); << 106 mode->SetDaughter(0, "pi0"); << 107 mode->SetDaughter(1, "pi0"); << 108 mode->SetDaughter(2, "pi-"); << 109 mode->SetDaughter(3, "nu_tau"); << 110 table->Insert(mode); << 111 // tau- -> pi- + pi- + pi+ + nu_tau << 112 mode = new G4PhaseSpaceDecayChannel(); << 113 mode->SetParent("tau-"); << 114 mode->SetBR(0.0946); << 115 mode->SetNumberOfDaughters(4); << 116 mode->SetDaughter(0, "pi-"); << 117 mode->SetDaughter(1, "pi-"); << 118 mode->SetDaughter(2, "pi+"); << 119 mode->SetDaughter(3, "nu_tau"); << 120 table->Insert(mode); << 121 << 122 anInstance->SetDecayTable(table); << 123 } 125 } 124 theInstance = static_cast<G4TauMinus*>(anIns << 126 theInstance = reinterpret_cast<G4TauMinus*>(anInstance); 125 return theInstance; 127 return theInstance; 126 } 128 } 127 129 128 G4TauMinus* G4TauMinus::TauMinusDefinition() << 130 G4TauMinus* G4TauMinus::TauMinusDefinition() 129 { 131 { 130 return Definition(); 132 return Definition(); 131 } 133 } 132 134 133 G4TauMinus* G4TauMinus::TauMinus() << 135 G4TauMinus* G4TauMinus::TauMinus() 134 { 136 { 135 return Definition(); 137 return Definition(); 136 } 138 } >> 139 >> 140 137 141