Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/dna/utils/src/G4DNAPTBExcitationStructure.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/dna/utils/src/G4DNAPTBExcitationStructure.cc (Version 11.3.0) and /processes/electromagnetic/dna/utils/src/G4DNAPTBExcitationStructure.cc (Version 9.2.p4)


  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 #include "G4DNAPTBExcitationStructure.hh"         
 27                                                   
 28 #include "G4Material.hh"                          
 29 #include "G4SystemOfUnits.hh"                     
 30                                                   
 31 G4DNAPTBExcitationStructure::G4DNAPTBExcitatio    
 32 {                                                 
 33   fpN2 = G4Material::GetMaterial("N2", false);    
 34                                                   
 35   //  taken directly from PTra code by MPietrz    
 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] : energyCon    
 70     nExcLevels[index] = (G4int)levels.size();     
 71   }                                               
 72 }                                                 
 73                                                   
 74 G4double G4DNAPTBExcitationStructure::Excitati    
 75   const G4int& ExcLevel, const size_t& materia    
 76 {                                                 
 77   size_t matNameModif = ReplaceMaterial(materi    
 78                                                   
 79   // check if the material exist in the map       
 80   if (energyConstant.find(matNameModif) == ene    
 81     std::ostringstream oss;                       
 82     oss << "Material name was not found in ene    
 83         << matNameModif;                          
 84     G4Exception(                                  
 85       "G4DNAPTBExcitationStructure::Excitation    
 86   }                                               
 87                                                   
 88   G4double excitation = 0.;                       
 89                                                   
 90   if (ExcLevel >= 0 && ExcLevel < nExcLevels[m    
 91     excitation = energyConstant[matNameModif][    
 92                                                   
 93   return excitation;                              
 94 }                                                 
 95                                                   
 96 G4int G4DNAPTBExcitationStructure::NumberOfExc    
 97 {                                                 
 98   auto matNameModif = ReplaceMaterial(matID);     
 99                                                   
100   // check if the material exist in the map       
101   if (nExcLevels.find(matNameModif) == nExcLev    
102     std::ostringstream oss;                       
103     oss << "Material name was not found in ene    
104         << matNameModif;                          
105     G4Exception("G4DNAPTBNDExcitationStructure    
106       oss.str().c_str());                         
107   }                                               
108                                                   
109   return nExcLevels[matNameModif];                
110 }                                                 
111                                                   
112 size_t G4DNAPTBExcitationStructure::ReplaceMat    
113 {                                                 
114   auto output = materialID;                       
115   auto G4_N2 = G4Material::GetMaterial("G4_N2"    
116   if (materialID == G4_N2) {                      
117     output = fpN2->GetIndex();                    
118   }                                               
119                                                   
120   return output;                                  
121 }                                                 
122