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 // 27 // File: G4ElementSelector 28 // 29 // Author: V.Ivanchenko (Vladimir.Ivanc 30 // 31 // Creation date: 2 April 2000 32 // 33 // Modifications: 34 // 18/08/2000 V.Ivanchenko Update description 35 // 17/05/2006 V.Ivanchenko Cleanup 36 // 02/10/2007 V.Ivanchenko Fixed typo in comp 37 // proposed by Victor 38 // 39 //-------------------------------------------- 40 41 #include "G4ElementSelector.hh" 42 #include "Randomize.hh" 43 #include "G4Material.hh" 44 #include "G4Nucleus.hh" 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 48 G4ElementSelector::G4ElementSelector() 49 {} 50 51 //....oooOO0OOooo........oooOO0OOooo........oo 52 53 G4ElementSelector::~G4ElementSelector() 54 {} 55 56 //....oooOO0OOooo........oooOO0OOooo........oo 57 58 const G4Element* 59 G4ElementSelector::SelectZandA(const G4Track& 60 { 61 // Fermi-Teller Z-low of mu- capture and exc 62 // for halogens and oxigen. 63 // N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1. 64 65 std::size_t i = 0; 66 const G4Material* mat = track.GetMaterial(); 67 std::size_t numberOfElements = mat->GetNumbe 68 const G4ElementVector* theElementVector = ma 69 70 if(1 < numberOfElements) { 71 if(numberOfElements > prob.size()) { prob. 72 73 const G4double* theAtomNumDensity = mat->G 74 75 G4double sum = 0.0; 76 for (i=0; i < numberOfElements; ++i) { 77 78 G4int Z = (*theElementVector)[i]->GetZas 79 80 // Halogens 81 if( (9 == Z) || (17 == Z) || (35 == Z) | 82 sum += 0.66 * Z * theAtomNumDensity[i]; 83 84 // Oxigen 85 } else if( 8 == Z ) { 86 sum += 0.56 * Z * theAtomNumDensity[i]; 87 88 // Others 89 } else { 90 sum += Z * theAtomNumDensity[i]; 91 } 92 prob[i] = sum; 93 } 94 95 sum *= G4UniformRand(); 96 for (i=0; i < numberOfElements; ++i) { 97 if(sum <= prob[i]) { break; } 98 } 99 } 100 101 const G4Element* elm = (*theElementVector)[i 102 G4int Z = elm->GetZasInt(); 103 104 // select isotope 105 const G4IsotopeVector* isv = elm->GetIsotope 106 std::size_t ni = isv->size(); 107 i = 0; 108 109 if(1 < ni) { 110 111 const G4double* ab = elm->GetRelativeAbund 112 G4double y = G4UniformRand(); 113 for(i=0; i<ni; ++i) { 114 y -= ab[i]; 115 if(y <= 0.0) { break; } 116 } 117 } 118 119 G4int A = elm->GetIsotope((G4int)i)->GetN(); 120 target->SetParameters(A, Z); 121 122 return elm; 123 } 124