Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 26 #include "G4DNAPTBExcitationStructure.hh" 27 28 #include "G4Material.hh" 29 #include "G4SystemOfUnits.hh" 30 31 G4DNAPTBExcitationStructure::G4DNAPTBExcitationStructure() 32 { 33 fpN2 = G4Material::GetMaterial("N2", false); 34 35 // taken directly from PTra code by MPietrzak 36 if (fpN2 != nullptr) { 37 auto index = fpN2->GetIndex(); 38 energyConstant[index].push_back(1.85 * eV); 39 energyConstant[index].push_back(2.15 * eV); 40 energyConstant[index].push_back(8.00 * eV); 41 energyConstant[index].push_back(8.50 * eV); 42 energyConstant[index].push_back(8.60 * eV); 43 energyConstant[index].push_back(11.05 * eV); 44 energyConstant[index].push_back(11.79 * eV); 45 energyConstant[index].push_back(11.90 * eV); 46 energyConstant[index].push_back(12.25 * eV); 47 energyConstant[index].push_back(12.50 * eV); 48 energyConstant[index].push_back(13.01 * eV); 49 energyConstant[index].push_back(13.19 * eV); 50 energyConstant[index].push_back(13.30 * eV); 51 energyConstant[index].push_back(14.33 * eV); 52 energyConstant[index].push_back(14.84 * eV); 53 energyConstant[index].push_back(15.18 * eV); 54 energyConstant[index].push_back(15.70 * eV); 55 energyConstant[index].push_back(15.75 * eV); 56 energyConstant[index].push_back(15.86 * eV); 57 energyConstant[index].push_back(17.36 * eV); 58 energyConstant[index].push_back(17.95 * eV); 59 energyConstant[index].push_back(19.77 * eV); 60 energyConstant[index].push_back(20.79 * eV); 61 energyConstant[index].push_back(20.87 * eV); 62 energyConstant[index].push_back(22.27 * eV); 63 energyConstant[index].push_back(22.83 * eV); 64 energyConstant[index].push_back(37.19 * eV); 65 energyConstant[index].push_back(38.67 * eV); 66 energyConstant[index].push_back(39.23 * eV); 67 } 68 69 for (const auto& [index, levels] : energyConstant) { 70 nExcLevels[index] = (G4int)levels.size(); 71 } 72 } 73 74 G4double G4DNAPTBExcitationStructure::ExcitationEnergy( 75 const G4int& ExcLevel, const size_t& materialID) 76 { 77 size_t matNameModif = ReplaceMaterial(materialID); 78 79 // check if the material exist in the map 80 if (energyConstant.find(matNameModif) == energyConstant.end()) { 81 std::ostringstream oss; 82 oss << "Material name was not found in energyConstantMap. Problematic material is: " 83 << matNameModif; 84 G4Exception( 85 "G4DNAPTBExcitationStructure::ExcitationEnergy", "em0002", FatalException, oss.str().c_str()); 86 } 87 88 G4double excitation = 0.; 89 90 if (ExcLevel >= 0 && ExcLevel < nExcLevels[matNameModif]) 91 excitation = energyConstant[matNameModif][ExcLevel]; 92 93 return excitation; 94 } 95 96 G4int G4DNAPTBExcitationStructure::NumberOfExcLevels(const size_t& matID) 97 { 98 auto matNameModif = ReplaceMaterial(matID); 99 100 // check if the material exist in the map 101 if (nExcLevels.find(matNameModif) == nExcLevels.end()) { 102 std::ostringstream oss; 103 oss << "Material name was not found in energyConstantMap. Problematic material is: " 104 << matNameModif; 105 G4Exception("G4DNAPTBNDExcitationStructure::NumberOfExcLevels", "em0002", FatalException, 106 oss.str().c_str()); 107 } 108 109 return nExcLevels[matNameModif]; 110 } 111 112 size_t G4DNAPTBExcitationStructure::ReplaceMaterial(const size_t& materialID) 113 { 114 auto output = materialID; 115 auto G4_N2 = G4Material::GetMaterial("G4_N2", false)->GetIndex(); 116 if (materialID == G4_N2) { 117 output = fpN2->GetIndex(); 118 } 119 120 return output; 121 } 122