Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // Author: 2021 Alberto Ribon 27 // 28 //-------------------------------------------- 29 30 #include "G4HyperHe5.hh" 31 32 #include "G4DecayTable.hh" 33 #include "G4ParticleTable.hh" 34 #include "G4PhaseSpaceDecayChannel.hh" 35 #include "G4PhysicalConstants.hh" 36 #include "G4String.hh" 37 #include "G4SystemOfUnits.hh" 38 #include "G4Types.hh" 39 #include "G4VDecayChannel.hh" 40 41 G4HyperHe5* G4HyperHe5::theInstance = nullptr; 42 43 G4HyperHe5* G4HyperHe5::Definition() 44 { 45 if (theInstance != nullptr) return theInstan 46 const G4String name = "hyperHe5"; 47 // search in particle table 48 G4ParticleTable* pTable = G4ParticleTable::G 49 auto anInstance = static_cast<G4Ions*>(pTabl 50 if (anInstance == nullptr) { 51 // create particle 52 // 53 // Arguments for constructor are as fol 54 // name 55 // 2*spin 56 // 2*Isospin 2 57 // type lept 58 // stable 59 // shortlived 60 // excitation 61 62 // clang-format off 63 anInstance = new G4Ions( name, 48 64 3, 65 0, 66 "nucleus", 67 false, 68 false, "static", 69 0.0, 0 70 // clang-format on 71 72 // Magnetic Moment 73 G4double mN = eplus * hbar_Planck / 2.0 / 74 anInstance->SetPDGMagneticMoment(2.9789624 75 76 // create Decay Table 77 auto table = new G4DecayTable; 78 // create decay channels 79 /* 80 // The decays "mode[1]" and "mode[3]" prod 81 // respectively, whose corresponding anti- 82 // we therefore skip them for the time bei 83 // of anti_hyperHe5. 84 const G4double half_br_lambda_to_p_pim = 0 85 const G4double half_br_lambda_to_n_piz = 0 86 G4VDecayChannel** mode = new G4VDecayChann 87 // lambda -> proton + pi- , with 50% proba 88 mode[0] = new G4PhaseSpaceDecayChannel( "h 89 "alpha", "proton", "pi-" ); 90 mode[1] = new G4PhaseSpaceDecayChannel( "h 91 "Li5", "pi-" ); 92 // lambda -> neutron + pi0 , with 50% prob 93 mode[2] = new G4PhaseSpaceDecayChannel( "h 94 "alpha", "neutron", "pi0" ); 95 mode[3] = new G4PhaseSpaceDecayChannel( "h 96 "He5", "pi0" ); 97 for ( G4int index = 0; index < 4; ++index 98 */ 99 // Replacement decay for the time being 100 const G4double br_lambda_to_p_pim = 0.639; 101 const G4double br_lambda_to_n_piz = 0.358; 102 auto mode = new G4VDecayChannel*[2]; 103 // lambda -> proton + pi- , with 0% probab 104 mode[0] = 105 new G4PhaseSpaceDecayChannel("hyperHe5", 106 // lambda -> neutron + pi0 , with 0% proba 107 mode[1] = 108 new G4PhaseSpaceDecayChannel("hyperHe5", 109 for (G4int index = 0; index < 2; ++index) 110 table->Insert(mode[index]); 111 //--- 112 delete[] mode; 113 anInstance->SetDecayTable(table); 114 } 115 theInstance = static_cast<G4HyperHe5*>(anIns 116 return theInstance; 117 } 118 119 G4HyperHe5* G4HyperHe5::HyperHe5Definition() 120 { 121 return Definition(); 122 } 123 124 G4HyperHe5* G4HyperHe5::HyperHe5() 125 { 126 return Definition(); 127 } 128