- Timestamp:
- Jun 18, 2010, 11:42:07 AM (14 years ago)
- Location:
- trunk/source/processes/hadronic/models/de_excitation/fermi_breakup
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/de_excitation/fermi_breakup/include/G4FermiPhaseSpaceDecay.hh
r819 r1315 33 33 #include "G4LorentzVector.hh" 34 34 #include "G4ParticleMomentum.hh" 35 #include "Randomize.hh" 36 #include "G4Pow.hh" 35 37 #include <CLHEP/Random/RandGamma.h> 36 38 … … 41 43 { 42 44 public: 43 inline G4FermiPhaseSpaceDecay(); 44 inline ~G4FermiPhaseSpaceDecay(); 45 46 G4FermiPhaseSpaceDecay(); 47 ~G4FermiPhaseSpaceDecay(); 45 48 46 49 inline std::vector<G4LorentzVector*> * … … 48 51 49 52 private: 50 inline G4FermiPhaseSpaceDecay(const G4FermiPhaseSpaceDecay&); 51 inline const G4FermiPhaseSpaceDecay & operator=(const G4FermiPhaseSpaceDecay &); 52 inline G4bool operator==(const G4FermiPhaseSpaceDecay&); 53 inline G4bool operator!=(const G4FermiPhaseSpaceDecay&); 53 54 G4FermiPhaseSpaceDecay(const G4FermiPhaseSpaceDecay&); 55 const G4FermiPhaseSpaceDecay & operator=(const G4FermiPhaseSpaceDecay &); 56 G4bool operator==(const G4FermiPhaseSpaceDecay&); 57 G4bool operator!=(const G4FermiPhaseSpaceDecay&); 54 58 55 59 inline G4double PtwoBody(G4double E, G4double P1, G4double P2) const; … … 69 73 }; 70 74 71 #include "G4FermiPhaseSpaceDecay.icc" 75 inline G4double 76 G4FermiPhaseSpaceDecay::PtwoBody(G4double E, G4double P1, G4double P2) const 77 { 78 G4double res = -1.0; 79 G4double P = (E+P1+P2)*(E+P1-P2)*(E-P1+P2)*(E-P1-P2)/(4.0*E*E); 80 if (P>0.0) { res = std::sqrt(P); } 81 return res; 82 } 83 84 inline std::vector<G4LorentzVector*> * G4FermiPhaseSpaceDecay:: 85 Decay(const G4double parent_mass, const std::vector<G4double>& fragment_masses) const 86 { 87 return KopylovNBodyDecay(parent_mass,fragment_masses); 88 } 89 90 inline G4double 91 G4FermiPhaseSpaceDecay::BetaKopylov(const G4int K) const 92 { 93 //JMQ 250410 old algorithm has been commented 94 // Notice that alpha > beta always 95 // const G4double beta = 1.5; 96 // G4double alpha = 1.5*(K-1); 97 // G4double Y1 = CLHEP::RandGamma::shoot(alpha,1); 98 // G4double Y2 = CLHEP::RandGamma::shoot(beta,1); 99 100 // return Y1/(Y1+Y2); 101 102 G4Pow* g4pow = G4Pow::GetInstance(); 103 G4int N = 3*K - 5; 104 G4double xN = G4double(N); 105 G4double F; 106 //G4double Fmax = std::pow((3.*K-5.)/(3.*K-4.),(3.*K-5.)/2.)*std::sqrt(1-((3.*K-5.)/(3.*K-4.))); 107 // VI variant 108 G4double Fmax = std::sqrt(g4pow->powZ(N, xN/(xN + 1))/(xN + 1)); 109 G4double chi; 110 do 111 { 112 chi = G4UniformRand(); 113 F = std::sqrt(g4pow->powZ(N, chi)*(1-chi)); 114 } while ( Fmax*G4UniformRand() > F); 115 return chi; 116 117 } 72 118 73 119 #endif -
trunk/source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiConfiguration.cc
r1228 r1315 25 25 // 26 26 // 27 // $Id: G4FermiConfiguration.cc,v 1.1 2 2009/12/16 17:51:09 gunterExp $28 // GEANT4 tag $Name: geant4-09-0 3$27 // $Id: G4FermiConfiguration.cc,v 1.13 2010/04/26 11:14:28 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $ 29 29 // 30 30 // Hadronic Process: Nuclear De-excitations 31 31 // by V. Lara (Nov 1998) 32 32 // 33 // J. M.Quesada (12 October 2009) new implementation of Gamma function in configuration weight34 33 // J. M. Quesada (12 October 2009) new implementation of Gamma function in configuration weight 34 // J. M. Quesada (09 March 2010) Kappa is set to 6. 35 35 36 36 #include "G4FermiConfiguration.hh" … … 40 40 // Kappa = V/V_0 it is used in calculation of Coulomb energy 41 41 // Kappa is adimensional 42 const G4double G4FermiConfiguration::Kappa = 1.0; 42 // JMQ 090310 according to model developer (A. Botvina) no theoretical constraint for kappa below 10 43 // kappa values larger than 1 seem to provide better results. 6 is a good choice 44 // const G4double G4FermiConfiguration::Kappa = 1.0; 45 const G4double G4FermiConfiguration::Kappa = 6.0; 43 46 44 47 // r0 is the nuclear radius -
trunk/source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiFragmentsPool.cc
r1196 r1315 28 28 // by V. Lara 29 29 // 30 // J.M.Quesada (JMQ)July 2009, bug fixed in excitation energies:30 // J.M.Quesada, July 2009, bug fixed in excitation energies: 31 31 // ALL of them are in MeV instead of keV (as they were expressed previously) 32 32 // source: http://www.nndc.bnl.gov/chart 33 33 // Unknown excitation energies in He5 and Li5 have been suppressed 34 34 // Long lived levels (half-lives of the order ps-fs have been included) 35 35 // 36 // J. M. Quesada, April 2010: excitation energies according to tabulated values 37 // in PhotonEvaporatoion2.0. Fake photons eliminated. 36 38 37 39 #include "G4FermiFragmentsPool.hh" … … 51 53 static const G4StableFermiFragment Fragment04( 3, 2, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment04); 52 54 static const G4StableFermiFragment Fragment05( 4, 2, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment05); 53 //JMQ 30/06/09 unknown levels have been supressed55 //JMQ 30/06/09 unknown levels have been supressed 54 56 static const G4He5FermiFragment Fragment06( 5, 2, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment06);// He5 55 57 static const G4Li5FermiFragment Fragment07( 5, 3, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment07);// Li5 56 58 static const G4StableFermiFragment Fragment08( 6, 2, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment08); 57 59 static const G4StableFermiFragment Fragment09( 6, 3, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment09); 58 static const G4StableFermiFragment Fragment10( 6, 3, 1, 3.56*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment10); 60 // static const G4StableFermiFragment Fragment10( 6, 3, 1, 3.56*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment10); 61 static const G4StableFermiFragment Fragment10( 6, 3, 1, 3.562880*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment10); 59 62 static const G4StableFermiFragment Fragment11( 7, 3, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment11); 60 static const G4StableFermiFragment Fragment12( 7, 3, 2, 0.48*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment12); 63 // static const G4StableFermiFragment Fragment12( 7, 3, 2, 0.48*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment12); 64 static const G4StableFermiFragment Fragment12( 7, 3, 2, 0.4776120*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment12); 61 65 static const G4StableFermiFragment Fragment13( 7, 4, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment13); 62 static const G4StableFermiFragment Fragment14( 7, 4, 2, 0.43*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment14); 66 // static const G4StableFermiFragment Fragment14( 7, 4, 2, 0.43*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment14); 67 static const G4StableFermiFragment Fragment14( 7, 4, 2, 0.4290800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment14); 63 68 static const G4StableFermiFragment Fragment15( 8, 3, 5, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment15); 64 static const G4StableFermiFragment Fragment16( 8, 3, 3, 0.98*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment16); 69 // static const G4StableFermiFragment Fragment16( 8, 3, 3, 0.98*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment16); 70 static const G4StableFermiFragment Fragment16( 8, 3, 3, 0.9808000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment16); 65 71 static const G4Be8FermiFragment Fragment17( 8, 4, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment17); // Be8 66 72 static const G4StableFermiFragment Fragment18( 9, 4, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment18); 67 73 static const G4B9FermiFragment Fragment19( 9, 5, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment19); // B9 68 74 static const G4StableFermiFragment Fragment20( 10, 4, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment20); 69 static const G4StableFermiFragment Fragment21( 10, 4, 5, 3.37*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment21); 70 static const G4StableFermiFragment Fragment22( 10, 4, 8, 5.96*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment22); 71 static const G4StableFermiFragment Fragment23( 10, 4, 1, 6.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment23); 72 static const G4StableFermiFragment Fragment24( 10, 4, 5, 6.26*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment24); 75 // static const G4StableFermiFragment Fragment21( 10, 4, 5, 3.37*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment21); 76 static const G4StableFermiFragment Fragment21( 10, 4, 5, 3.368030*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment21); 77 // static const G4StableFermiFragment Fragment22( 10, 4, 8, 5.96*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment22); 78 static const G4StableFermiFragment Fragment22( 10, 4, 8, 5.958390*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment22); 79 // static const G4StableFermiFragment Fragment23( 10, 4, 1, 6.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment23); 80 static const G4StableFermiFragment Fragment23( 10, 4, 1, 6.179300*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment23); 81 // static const G4StableFermiFragment Fragment24( 10, 4, 5, 6.26*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment24); 82 static const G4StableFermiFragment Fragment24( 10, 4, 5, 6.263300*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment24); 73 83 static const G4StableFermiFragment Fragment25( 10, 5, 7, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment25); 74 static const G4StableFermiFragment Fragment26( 10, 5, 3, 0.72*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment26); 75 static const G4StableFermiFragment Fragment27( 10, 5, 1, 1.74*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment27); 76 static const G4StableFermiFragment Fragment28( 10, 5, 3, 2.15*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment28); 77 static const G4StableFermiFragment Fragment29( 10, 5, 5, 3.59*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment29); 78 84 // static const G4StableFermiFragment Fragment26( 10, 5, 3, 0.72*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment26); 85 static const G4StableFermiFragment Fragment26( 10, 5, 3, 0.7183500*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment26); 86 // static const G4StableFermiFragment Fragment27( 10, 5, 1, 1.74*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment27); 87 static const G4StableFermiFragment Fragment27( 10, 5, 1, 1.740150*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment27); 88 // static const G4StableFermiFragment Fragment28( 10, 5, 3, 2.15*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment28); 89 static const G4StableFermiFragment Fragment28( 10, 5, 3, 2.154300*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment28); 90 // static const G4StableFermiFragment Fragment29( 10, 5, 5, 3.59*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment29); 91 static const G4StableFermiFragment Fragment29( 10, 5, 5, 3.587100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment29); 79 92 static const G4StableFermiFragment Fragment30( 10, 6, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment30); 80 static const G4StableFermiFragment Fragment31( 10, 6, 5, 3.35*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment31); 93 // static const G4StableFermiFragment Fragment31( 10, 6, 5, 3.35*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment31); 94 static const G4StableFermiFragment Fragment31( 10, 6, 5, 3.353600*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment31); 81 95 static const G4StableFermiFragment Fragment32( 11, 5, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment32); 82 static const G4StableFermiFragment Fragment33( 11, 5, 2, 2.13*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment33); 83 static const G4StableFermiFragment Fragment34( 11, 5, 6, 4.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment34); 84 static const G4StableFermiFragment Fragment35( 11, 5, 4, 5.02*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment35); 85 static const G4StableFermiFragment Fragment36( 11, 5, 10, 6.76*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment36); 86 static const G4StableFermiFragment Fragment37( 11, 5, 6, 7.29*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment37); 87 static const G4StableFermiFragment Fragment38( 11, 5, 4, 7.98*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment38); 88 static const G4StableFermiFragment Fragment39( 11, 5, 6, 8.56*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment39); 89 90 static const G4StableFermiFragment Fragment40( 11, 6, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment40); 91 static const G4StableFermiFragment Fragment41( 11, 6, 2, 2.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment41); 92 static const G4StableFermiFragment Fragment42( 11, 6, 6, 4.32*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment42); 93 static const G4StableFermiFragment Fragment43( 11, 6, 4, 4.80*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment43); 94 static const G4StableFermiFragment Fragment44( 11, 6, 2, 6.34*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment44); 95 static const G4StableFermiFragment Fragment45( 11, 6, 8, 6.48*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment45); 96 static const G4StableFermiFragment Fragment46( 11, 6, 6, 6.90*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment46); 97 static const G4StableFermiFragment Fragment47( 11, 6, 4, 7.50*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment47); 98 static const G4StableFermiFragment Fragment48( 11, 6, 4, 8.10*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment48); 99 static const G4StableFermiFragment Fragment49( 11, 6, 6, 8.42*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment49); 100 101 static const G4StableFermiFragment Fragment50( 12, 5, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment50); 102 static const G4StableFermiFragment Fragment51( 12, 5, 5, 0.95*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment51); 103 static const G4StableFermiFragment Fragment52( 12, 5, 5, 1.67*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment52); 104 static const G4StableFermiFragment Fragment53( 12, 5, 4, 2.65*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment53); 105 static const G4StableFermiFragment Fragment54( 12, 6, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment54); 106 static const G4StableFermiFragment Fragment55( 12, 6, 5, 4.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment55); 107 static const G4StableFermiFragment Fragment56( 13, 6, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment56); 108 static const G4StableFermiFragment Fragment57( 13, 6, 2, 3.09*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment57); 109 static const G4StableFermiFragment Fragment58( 13, 6, 4, 3.68*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment58); 110 111 static const G4StableFermiFragment Fragment59( 13, 6, 6, 3.85*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment59); 112 static const G4StableFermiFragment Fragment60( 13, 7, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment60); 113 static const G4StableFermiFragment Fragment61( 14, 6, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment61); 114 static const G4StableFermiFragment Fragment62( 14, 6, 3, 6.09*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment62); 115 // JMQ 010709 corrected excitation energies for 64-66, according to http://www.nndc.bnl.gov/chart 116 static const G4StableFermiFragment Fragment63( 14, 6, 1, 6.59*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment63); 117 static const G4StableFermiFragment Fragment64( 14, 6, 7, 6.73*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment64); 118 static const G4StableFermiFragment Fragment65( 14, 6, 1, 6.90*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment65); 119 static const G4StableFermiFragment Fragment66( 14, 6, 5, 7.01*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment66); 120 static const G4StableFermiFragment Fragment67( 14, 6, 5, 7.34*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment67); 121 122 static const G4StableFermiFragment Fragment68( 14, 7, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment68); 123 static const G4StableFermiFragment Fragment69( 14, 7, 1, 2.31*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment69); 124 static const G4StableFermiFragment Fragment70( 14, 7, 3, 3.95*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment70); 125 126 static const G4StableFermiFragment Fragment71( 14, 7, 1, 4.92*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment71); 127 static const G4StableFermiFragment Fragment72( 14, 7, 5, 5.11*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment72); 128 static const G4StableFermiFragment Fragment73( 14, 7, 3, 5.69*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment73); 129 static const G4StableFermiFragment Fragment74( 14, 7, 7, 5.83*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment74); 130 static const G4StableFermiFragment Fragment75( 14, 7, 3, 6.20*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment75); 131 static const G4StableFermiFragment Fragment76( 14, 7, 7, 6.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment76); 132 static const G4StableFermiFragment Fragment77( 14, 7, 5, 7.03*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment77); 133 static const G4StableFermiFragment Fragment78( 15, 7, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment78); 134 // JMQ 010709 two very close levels instead of only one, with their own spins 135 static const G4StableFermiFragment Fragment79( 15, 7, 6, 5.27*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment79); 136 static const G4StableFermiFragment Fragment80( 15, 7, 2, 5.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment80); 137 static const G4StableFermiFragment Fragment81( 15, 7, 4, 6.32*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment81); 138 //JMQ 010709 new level and corrected energy and spins 139 static const G4StableFermiFragment Fragment82( 15, 7, 6, 7.15*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment82); 140 static const G4StableFermiFragment Fragment83( 15, 7, 4, 7.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment83); 141 static const G4StableFermiFragment Fragment84( 15, 7, 8, 7.57*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment84); 142 static const G4StableFermiFragment Fragment85( 15, 7, 2, 8.31*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment85); 143 static const G4StableFermiFragment Fragment86( 15, 7, 4, 8.57*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment86); 144 static const G4StableFermiFragment Fragment87( 15, 7, 2, 9.05*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment87); 145 //JMQ 010709 new levels for N15 146 static const G4StableFermiFragment Fragment88( 15, 7, 4, 9.151*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment88); 147 static const G4StableFermiFragment Fragment89( 15, 7, 6, 9.154*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment89); 148 static const G4StableFermiFragment Fragment90( 15, 7, 2, 9.22*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment90); 149 static const G4StableFermiFragment Fragment91( 15, 7, 6, 9.76*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment91); 150 static const G4StableFermiFragment Fragment92( 15, 7, 8, 9.83*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment92); 151 static const G4StableFermiFragment Fragment93( 15, 7, 4, 9.93*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment93); 152 static const G4StableFermiFragment Fragment94( 15, 7, 4, 10.07*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment94); 153 154 155 static const G4StableFermiFragment Fragment95( 15, 8, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment95); 156 //JMQ 010709 new level and spins 157 static const G4StableFermiFragment Fragment96( 15, 8, 2, 5.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment96); 158 static const G4StableFermiFragment Fragment97( 15, 8, 6, 5.24*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment97); 159 static const G4StableFermiFragment Fragment98( 15, 8, 4, 6.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment98); 160 static const G4StableFermiFragment Fragment99( 15, 8, 4, 6.79*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment99); 161 static const G4StableFermiFragment Fragment100( 15, 8, 6, 6.86*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment100); 162 static const G4StableFermiFragment Fragment101( 15, 8, 8, 7.28*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment101); 163 164 165 static const G4StableFermiFragment Fragment102( 16, 7, 5, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment102); 166 static const G4StableFermiFragment Fragment103( 16, 7, 1, 0.12*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment103); 167 static const G4StableFermiFragment Fragment104( 16, 7, 7, 0.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment104); 168 static const G4StableFermiFragment Fragment105( 16, 7, 3, 0.40*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment105); 169 170 //JMQ 010709 some energies and spins have been changed 171 static const G4StableFermiFragment Fragment106( 16, 8, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment106); 172 static const G4StableFermiFragment Fragment107( 16, 8, 1, 6.05*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment107); 173 static const G4StableFermiFragment Fragment108( 16, 8, 7, 6.13*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment108); 174 static const G4StableFermiFragment Fragment109( 16, 8, 5, 6.92*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment109); 175 static const G4StableFermiFragment Fragment110( 16, 8, 3, 7.12*MeV ); if(MapIsEmpty) fragment_pool.push_back 176 (&Fragment110); 96 // static const G4StableFermiFragment Fragment33( 11, 5, 2, 2.13*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment33); 97 static const G4StableFermiFragment Fragment33( 11, 5, 2, 2.124693*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment33); 98 // static const G4StableFermiFragment Fragment34( 11, 5, 6, 4.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment34); 99 static const G4StableFermiFragment Fragment34( 11, 5, 6, 4.444890*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment34); 100 // static const G4StableFermiFragment Fragment35( 11, 5, 4, 5.02*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment35); 101 static const G4StableFermiFragment Fragment35( 11, 5, 4, 5.020310*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment35); 102 // static const G4StableFermiFragment Fragment36( 11, 5, 10, 6.76*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment36); 103 static const G4StableFermiFragment Fragment36( 11, 5, 8, 6.742900*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment36); 104 //JMQ 190410 new level, fragment numbering shifted accordingly from here onwards 105 static const G4StableFermiFragment Fragment37( 11, 5, 2, 6.791800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment37); 106 // static const G4StableFermiFragment Fragment37( 11, 5, 6, 7.29*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment37); 107 static const G4StableFermiFragment Fragment38( 11, 5, 6, 7.285510*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment38); 108 // static const G4StableFermiFragment Fragment38( 11, 5, 4, 7.98*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment38); 109 static const G4StableFermiFragment Fragment39( 11, 5, 4, 7.977840*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment39); 110 // static const G4StableFermiFragment Fragment39( 11, 5, 6, 8.56*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment39); 111 static const G4StableFermiFragment Fragment40( 11, 5, 6, 8.560300*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment40); 112 static const G4StableFermiFragment Fragment41( 11, 6, 4, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment41); 113 static const G4StableFermiFragment Fragment42( 11, 6, 2, 2.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment42); 114 // static const G4StableFermiFragment Fragment42( 11, 6, 6, 4.32*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment42); 115 static const G4StableFermiFragment Fragment43( 11, 6, 6, 4.318800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment43); 116 // static const G4StableFermiFragment Fragment43( 11, 6, 4, 4.80*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment43); 117 static const G4StableFermiFragment Fragment44( 11, 6, 4, 4.804200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment44); 118 // static const G4StableFermiFragment Fragment44( 11, 6, 2, 6.34*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment44); 119 static const G4StableFermiFragment Fragment45( 11, 6, 2, 6.339200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment45); 120 // static const G4StableFermiFragment Fragment45( 11, 6, 8, 6.48*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment45); 121 static const G4StableFermiFragment Fragment46( 11, 6, 8, 6.478200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment46); 122 // static const G4StableFermiFragment Fragment46( 11, 6, 6, 6.90*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment46); 123 static const G4StableFermiFragment Fragment47( 11, 6, 6, 6.904800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment47); 124 // static const G4StableFermiFragment Fragment47( 11, 6, 4, 7.50*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment47); 125 static const G4StableFermiFragment Fragment48( 11, 6, 4, 7.499700*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment48); 126 // static const G4StableFermiFragment Fragment48( 11, 6, 4, 8.10*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment48); 127 static const G4StableFermiFragment Fragment49( 11, 6, 4, 8.104500*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment49); 128 // static const G4StableFermiFragment Fragment49( 11, 6, 6, 8.42*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment49); 129 static const G4StableFermiFragment Fragment50( 11, 6, 6, 8.420000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment50); 130 static const G4StableFermiFragment Fragment51( 12, 5, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment51); 131 // static const G4StableFermiFragment Fragment51( 12, 5, 5, 0.95*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment51); 132 static const G4StableFermiFragment Fragment52( 12, 5, 5, 0.9531400*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment52); 133 // static const G4StableFermiFragment Fragment52( 12, 5, 5, 1.67*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment52); 134 static const G4StableFermiFragment Fragment53( 12, 5, 5, 1.673650*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment53); 135 // static const G4StableFermiFragment Fragment53( 12, 5, 4, 2.65*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment53); 136 static const G4StableFermiFragment Fragment54( 12, 5, 3, 2.620800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment54); 137 static const G4StableFermiFragment Fragment55( 12, 6, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment55); 138 // static const G4StableFermiFragment Fragment55( 12, 6, 5, 4.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment55); 139 static const G4StableFermiFragment Fragment56( 12, 6, 5, 4.438910*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment56); 140 static const G4StableFermiFragment Fragment57( 13, 6, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment57); 141 // static const G4StableFermiFragment Fragment57( 13, 6, 2, 3.09*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment57); 142 static const G4StableFermiFragment Fragment58( 13, 6, 2, 3.089443*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment58); 143 // static const G4StableFermiFragment Fragment58( 13, 6, 4, 3.68*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment58); 144 static const G4StableFermiFragment Fragment59( 13, 6, 4, 3.684507*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment59); 145 // static const G4StableFermiFragment Fragment59( 13, 6, 6, 3.85*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment59); 146 static const G4StableFermiFragment Fragment60( 13, 6, 6, 3.853807*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment60); 147 static const G4StableFermiFragment Fragment61( 13, 7, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment61); 148 static const G4StableFermiFragment Fragment62( 14, 6, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment62); 149 // static const G4StableFermiFragment Fragment62( 14, 6, 3, 6.09*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment62); 150 static const G4StableFermiFragment Fragment63( 14, 6, 3, 6.093800*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment63); 151 // JMQ 010709 corrected excitation energies for 64-66, according to http://www.nndc.bnl.gov/chart 152 // static const G4StableFermiFragment Fragment63( 14, 6, 1, 6.59*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment63); 153 static const G4StableFermiFragment Fragment64( 14, 6, 1, 6.589400*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment64); 154 // static const G4StableFermiFragment Fragment64( 14, 6, 7, 6.73*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment64); 155 static const G4StableFermiFragment Fragment65( 14, 6, 7, 6.728200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment65); 156 // static const G4StableFermiFragment Fragment65( 14, 6, 1, 6.90*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment65); 157 static const G4StableFermiFragment Fragment66( 14, 6, 1, 6.902600*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment66); 158 // static const G4StableFermiFragment Fragment66( 14, 6, 5, 7.01*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment66); 159 static const G4StableFermiFragment Fragment67( 14, 6, 5, 7.012000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment67); 160 // static const G4StableFermiFragment Fragment67( 14, 6, 5, 7.34*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment67); 161 static const G4StableFermiFragment Fragment68( 14, 6, 5, 7.341000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment68); 162 163 static const G4StableFermiFragment Fragment69( 14, 7, 3, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment69); 164 // static const G4StableFermiFragment Fragment69( 14, 7, 1, 2.31*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment69); 165 static const G4StableFermiFragment Fragment70( 14, 7, 1, 2.312798*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment70); 166 // static const G4StableFermiFragment Fragment70( 14, 7, 3, 3.95*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment70); 167 static const G4StableFermiFragment Fragment71( 14, 7, 3, 3.948100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment71); 168 // static const G4StableFermiFragment Fragment71( 14, 7, 1, 4.92*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment71); 169 static const G4StableFermiFragment Fragment72( 14, 7, 1, 4.915100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment72); 170 // static const G4StableFermiFragment Fragment72( 14, 7, 5, 5.11*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment72); 171 static const G4StableFermiFragment Fragment73( 14, 7, 5, 5.105890*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment73); 172 // static const G4StableFermiFragment Fragment73( 14, 7, 3, 5.69*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment73); 173 static const G4StableFermiFragment Fragment74( 14, 7, 3, 5.691440*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment74); 174 // static const G4StableFermiFragment Fragment74( 14, 7, 7, 5.83*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment74); 175 static const G4StableFermiFragment Fragment75( 14, 7, 7, 5.834250*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment75); 176 // static const G4StableFermiFragment Fragment75( 14, 7, 3, 6.20*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment75); 177 static const G4StableFermiFragment Fragment76( 14, 7, 3, 6.203500*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment76); 178 // static const G4StableFermiFragment Fragment76( 14, 7, 7, 6.44*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment76); 179 static const G4StableFermiFragment Fragment77( 14, 7, 7, 6.446170*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment77); 180 // static const G4StableFermiFragment Fragment77( 14, 7, 5, 7.03*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment77); 181 static const G4StableFermiFragment Fragment78( 14, 7, 5, 7.029120*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment78); 182 static const G4StableFermiFragment Fragment79( 15, 7, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment79); 183 // JMQ 010709 two very close levels instead of only one, with their own spins 184 // static const G4StableFermiFragment Fragment79( 15, 7, 6, 5.27*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment79); 185 static const G4StableFermiFragment Fragment80( 15, 7, 6, 5.270155*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment80); 186 // static const G4StableFermiFragment Fragment80( 15, 7, 2, 5.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment80); 187 static const G4StableFermiFragment Fragment81( 15, 7, 2, 5.298822*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment81); 188 // static const G4StableFermiFragment Fragment81( 15, 7, 4, 6.32*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment81); 189 static const G4StableFermiFragment Fragment82( 15, 7, 4, 6.323780*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment82); 190 //JMQ 010709 new level and corrected energy and spins 191 // static const G4StableFermiFragment Fragment82( 15, 7, 6, 7.15*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment82); 192 static const G4StableFermiFragment Fragment83( 15, 7, 6, 7.155050*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment83); 193 // static const G4StableFermiFragment Fragment83( 15, 7, 4, 7.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment83); 194 static const G4StableFermiFragment Fragment84( 15, 7, 4, 7.300830*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment84); 195 // static const G4StableFermiFragment Fragment84( 15, 7, 8, 7.57*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment84); 196 static const G4StableFermiFragment Fragment85( 15, 7, 8, 7.567100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment85); 197 // static const G4StableFermiFragment Fragment85( 15, 7, 2, 8.31*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment85); 198 static const G4StableFermiFragment Fragment86( 15, 7, 2, 8.312620*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment86); 199 // static const G4StableFermiFragment Fragment86( 15, 7, 4, 8.57*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment86); 200 static const G4StableFermiFragment Fragment87( 15, 7, 4, 8.571400*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment87); 201 // static const G4StableFermiFragment Fragment87( 15, 7, 2, 9.05*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment87); 202 static const G4StableFermiFragment Fragment88( 15, 7, 2, 9.049710*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment88); 203 //JMQ 010709 new levels for N15 204 // static const G4StableFermiFragment Fragment88( 15, 7, 4, 9.151*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment88); 205 static const G4StableFermiFragment Fragment89( 15, 7, 4, 9.151900*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment89); 206 // static const G4StableFermiFragment Fragment89( 15, 7, 6, 9.154*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment89); 207 static const G4StableFermiFragment Fragment90( 15, 7, 6, 9.154900*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment90); 208 // static const G4StableFermiFragment Fragment90( 15, 7, 2, 9.22*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment90); 209 static const G4StableFermiFragment Fragment91( 15, 7, 2, 9.222100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment91); 210 // static const G4StableFermiFragment Fragment91( 15, 7, 6, 9.76*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment91); 211 static const G4StableFermiFragment Fragment92( 15, 7, 6, 9.760000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment92); 212 // static const G4StableFermiFragment Fragment92( 15, 7, 8, 9.83*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment92); 213 static const G4StableFermiFragment Fragment93( 15, 7, 8, 9.829000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment93); 214 // static const G4StableFermiFragment Fragment93( 15, 7, 4, 9.93*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment93); 215 static const G4StableFermiFragment Fragment94( 15, 7, 4, 9.925000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment94); 216 // static const G4StableFermiFragment Fragment94( 15, 7, 4, 10.07*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment94); 217 static const G4StableFermiFragment Fragment95( 15, 7, 4, 10.06600*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment95); 218 219 static const G4StableFermiFragment Fragment96( 15, 8, 2, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment96); 220 //JMQ 010709 new level and spins 221 // static const G4StableFermiFragment Fragment96( 15, 8, 2, 5.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment96); 222 static const G4StableFermiFragment Fragment97( 15, 8, 2, 5.183000*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment97); 223 // static const G4StableFermiFragment Fragment97( 15, 8, 6, 5.24*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment97); 224 static const G4StableFermiFragment Fragment98( 15, 8, 6, 5.240900*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment98); 225 // static const G4StableFermiFragment Fragment98( 15, 8, 4, 6.18*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment98); 226 static const G4StableFermiFragment Fragment99( 15, 8, 4, 6.176300*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment99); 227 // static const G4StableFermiFragment Fragment99( 15, 8, 4, 6.79*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment99); 228 static const G4StableFermiFragment Fragment100( 15, 8, 4, 6.793100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment100); 229 // static const G4StableFermiFragment Fragment100( 15, 8, 6, 6.86*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment100); 230 static const G4StableFermiFragment Fragment101( 15, 8, 6, 6.859400*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment101); 231 // static const G4StableFermiFragment Fragment101( 15, 8, 8, 7.28*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment101); 232 static const G4StableFermiFragment Fragment102( 15, 8, 8, 7.275900*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment102); 233 234 235 static const G4StableFermiFragment Fragment103( 16, 7, 5, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment103); 236 // static const G4StableFermiFragment Fragment103( 16, 7, 1, 0.12*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment103); 237 static const G4StableFermiFragment Fragment104( 16, 7, 1, 0.1204200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment104); 238 // static const G4StableFermiFragment Fragment104( 16, 7, 7, 0.30*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment104); 239 static const G4StableFermiFragment Fragment105( 16, 7, 7, 0.2982200*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment105); 240 // static const G4StableFermiFragment Fragment105( 16, 7, 3, 0.40*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment105); 241 static const G4StableFermiFragment Fragment106( 16, 7, 3, 0.3972700*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment106); 242 243 //JMQ 010709 some energies and spins have been changed 244 static const G4StableFermiFragment Fragment107( 16, 8, 1, 0.00*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment107); 245 // static const G4StableFermiFragment Fragment107( 16, 8, 1, 6.05*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment107); 246 static const G4StableFermiFragment Fragment108( 16, 8, 1, 6.049400*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment108); 247 // static const G4StableFermiFragment Fragment108( 16, 8, 7, 6.13*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment108); 248 static const G4StableFermiFragment Fragment109( 16, 8, 7, 6.129890*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment109); 249 // static const G4StableFermiFragment Fragment109( 16, 8, 5, 6.92*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment109); 250 static const G4StableFermiFragment Fragment110( 16, 8, 5, 6.917100*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment110); 251 // static const G4StableFermiFragment Fragment110( 16, 8, 3, 7.12*MeV ); if(MapIsEmpty) fragment_pool.push_back 252 //JMQ 180510 fixed fragment 111 253 static const G4StableFermiFragment Fragment111( 16, 8, 3, 7.116850*MeV ); if(MapIsEmpty) fragment_pool.push_back(&Fragment111); 177 254 178 255 static std::multimap<const std::pair<G4int,G4int>, const G4VFermiFragment* , std::less<const std::pair<G4int,G4int> > > 179 256 theMapOfFragments; 180 257 181 258 if (MapIsEmpty) 182 259 { 183 260 for(size_t i=0; i<fragment_pool.size(); i++) 184 185 186 const G4VFermiFragment* >(std::pair<G4int,G4int>(fragment_pool[i]->GetA(),261 { 262 theMapOfFragments.insert(std::pair<const std::pair<G4int,G4int>, 263 const G4VFermiFragment* >(std::pair<G4int,G4int>(fragment_pool[i]->GetA(), 187 264 fragment_pool[i]->GetZ()),fragment_pool[i])); 188 265 } 189 266 MapIsEmpty = false; 190 267 } -
trunk/source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiPhaseSpaceDecay.cc
r819 r1315 37 37 #include <functional> 38 38 39 G4FermiPhaseSpaceDecay::G4FermiPhaseSpaceDecay() 40 { 41 } 42 43 G4FermiPhaseSpaceDecay::~G4FermiPhaseSpaceDecay() 44 { 45 } 39 46 40 47 std::vector<G4LorentzVector*> *
Note: See TracChangeset
for help on using the changeset viewer.