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 // 23-Jan-2009 V.Ivanchenko make the class to 28 // 17-Aug-2012 V.Ivanchenko added hadronic mod 29 30 #include "G4NuclearPolarizationStore.hh" 31 #include "G4SystemOfUnits.hh" 32 33 G4ThreadLocal G4NuclearPolarizationStore* 34 G4NuclearPolarizationStore::instance = nullptr 35 36 G4NuclearPolarizationStore* G4NuclearPolarizat 37 { 38 if(nullptr == instance) { 39 static G4ThreadLocalSingleton<G4NuclearPol 40 instance = inst.Instance(); 41 } 42 return instance; 43 } 44 45 G4NuclearPolarizationStore::G4NuclearPolarizat 46 { 47 for(G4int i=0; i<maxNumStates; ++i) { nuclis 48 oldIdx = 0; 49 } 50 51 G4NuclearPolarizationStore::~G4NuclearPolariza 52 { 53 for(G4int i=0; i<maxNumStates; ++i) { 54 delete nuclist[i]; 55 nuclist[i] = nullptr; 56 } 57 } 58 59 void G4NuclearPolarizationStore::Register(G4Nu 60 { 61 G4int idx = -1; 62 for(G4int i=0; i<maxNumStates; ++i) { 63 if(ptr == nuclist[i]) { return; } 64 if(nullptr == nuclist[i]) { idx = i; } 65 } 66 if(idx >= 0) { 67 nuclist[idx] = ptr; 68 return; 69 } 70 // delete oldest object 71 delete nuclist[oldIdx]; 72 nuclist[oldIdx] = ptr; 73 // redefine oldIdx 74 ++oldIdx; 75 if(oldIdx >= maxNumStates) { oldIdx = 0; } 76 } 77 78 G4NuclearPolarization* 79 G4NuclearPolarizationStore::FindOrBuild(G4int 80 { 81 static const G4double tolerance = 10.*CLHEP: 82 for(G4int i=0; i<maxNumStates; ++i) { 83 auto nucp = nuclist[i]; 84 if(nucp && Z == nucp->GetZ() && A == nucp- 85 std::abs(Eexc - nucp->GetExcitationEner 86 return nucp; 87 } 88 } 89 G4NuclearPolarization* ptr = new G4NuclearPo 90 Register(ptr); 91 return ptr; 92 } 93 94 void G4NuclearPolarizationStore::RemoveMe(G4Nu 95 { 96 for(G4int i=0; i<maxNumStates; ++i) { 97 if(ptr == nuclist[i]) { 98 delete ptr; 99 nuclist[i] = nullptr; 100 // do we need redefine oldIdx? 101 if(i == oldIdx) { 102 for(G4int j=0; j<maxNumStates; ++j) { 103 if(j != i && nullptr != nuclist[j]) { 104 oldIdx = j; 105 return; 106 } 107 } 108 oldIdx = i; 109 } 110 return; 111 } 112 } 113 } 114