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 // neutron_hp -- source file 27 // J.P. Wellisch, Nov-1996 28 // A prototype of the low energy neutron trans 29 // 30 // P. Arce, June-2014 Conversion neutron_hp to 31 // 32 #include "G4ParticleHPData.hh" 33 34 #include "G4PhysicsFreeVector.hh" 35 36 G4ParticleHPData::G4ParticleHPData(G4ParticleD 37 : theProjectile(projectile) 38 { 39 // const char* theDataDirVariable; 40 if (projectile == G4Neutron::Neutron()) { 41 theDataDirVariable = "G4NEUTRONHPDATA"; 42 } 43 else if (projectile == G4Proton::Proton()) { 44 theDataDirVariable = "G4PROTONHPDATA"; 45 } 46 else if (projectile == G4Deuteron::Deuteron( 47 theDataDirVariable = "G4DEUTERONHPDATA"; 48 } 49 else if (projectile == G4Triton::Triton()) { 50 theDataDirVariable = "G4TRITONHPDATA"; 51 } 52 else if (projectile == G4He3::He3()) { 53 theDataDirVariable = "G4HE3HPDATA"; 54 } 55 else if (projectile == G4Alpha::Alpha()) { 56 theDataDirVariable = "G4ALPHAHPDATA"; 57 } 58 59 numEle = (G4int)G4Element::GetNumberOfElemen 60 for (G4int i = 0; i < numEle; ++i) { 61 theData.push_back(new G4ParticleHPElementD 62 } 63 for (G4int i = 0; i < numEle; ++i) { 64 (*theData[i]).Init((*(G4Element::GetElemen 65 } 66 } 67 68 G4ParticleHPData::~G4ParticleHPData() 69 { 70 for (auto it = theData.cbegin(); it != theDa 71 delete *it; 72 theData.clear(); 73 } 74 75 G4ParticleHPData* G4ParticleHPData::Instance(G 76 { 77 static G4ThreadLocal G4ParticleHPData* theCr 78 if (theCrossSectionData_G4MT_TLS_ == nullptr 79 theCrossSectionData_G4MT_TLS_ = new G4Part 80 G4ParticleHPData& theCrossSectionData = *the 81 return &theCrossSectionData; 82 } 83 84 G4PhysicsVector* G4ParticleHPData::DoPhysicsVe 85 { 86 G4int len = theVector->GetVectorLength(); 87 if (len == 0) return new G4PhysicsFreeVector 88 G4double emin = theVector->GetX(0); 89 G4double emax = theVector->GetX(len - 1); 90 91 auto theResult = new G4PhysicsFreeVector(len 92 for (G4int i = 0; i < len; ++i) { 93 theResult->PutValues(i, theVector->GetX(i) 94 } 95 return theResult; 96 } 97 98 void G4ParticleHPData::addPhysicsVector() 99 { 100 for (G4int i = numEle; i < (G4int)G4Element: 101 theData.push_back(new G4ParticleHPElementD 102 (*theData[i]).Init((*(G4Element::GetElemen 103 } 104 numEle = (G4int)G4Element::GetNumberOfElemen 105 } 106