Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/lowenergy/src/G4ANSTOecpssrKxsModel.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/lowenergy/src/G4ANSTOecpssrKxsModel.cc (Version 11.3.0) and /processes/electromagnetic/lowenergy/src/G4ANSTOecpssrKxsModel.cc (Version 8.2.p1)


  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 // History:                                       
 27 // -----------                                    
 28 //  10 Nov 2021   S. Guatelli & S. Bakr, Wollo    
 29 //                                                
 30 // Class description                              
 31 // ----------------                               
 32 //  Computation of K, L & M shell ECPSSR ionis    
 33 //  Based on the work of                          
 34 //  - S. Bakr et al. (2021) NIM B, 507:11-19.     
 35 //  - S. Bakr et al (2018), NIMB B, 436: 285-2    
 36 // -------------------------------------------    
 37                                                   
 38 #include <fstream>                                
 39 #include <iomanip>                                
 40                                                   
 41 #include "G4ANSTOecpssrKxsModel.hh"               
 42                                                   
 43 #include "globals.hh"                             
 44 #include "G4SystemOfUnits.hh"                     
 45 #include "G4ios.hh"                               
 46                                                   
 47 #include "G4EMDataSet.hh"                         
 48 #include "G4LogLogInterpolation.hh"               
 49 #include "G4Proton.hh"                            
 50 #include "G4Alpha.hh"                             
 51                                                   
 52 //....oooOO0OOooo........oooOO0OOooo........oo    
 53                                                   
 54 G4ANSTOecpssrKxsModel::G4ANSTOecpssrKxsModel()    
 55 {                                                 
 56                                                   
 57   G4cout << "Using ANSTO K Cross Sections! "<<    
 58                                                   
 59   interpolation = new G4LogLogInterpolation();    
 60   for (G4int i=6; i<93; i++)                      
 61   {                                               
 62       protonDataSetMap[i] = new G4EMDataSet(i,    
 63       protonDataSetMap[i]->LoadData("pixe_ANST    
 64   }                                               
 65                                                   
 66   for (G4int i=6; i<93; i++)                      
 67   {                                               
 68       alphaDataSetMap[i] = new G4EMDataSet(i,i    
 69       alphaDataSetMap[i]->LoadData("pixe_ANSTO    
 70   }                                               
 71                                                   
 72 }                                                 
 73                                                   
 74 //....oooOO0OOooo........oooOO0OOooo........oo    
 75                                                   
 76 G4ANSTOecpssrKxsModel::~G4ANSTOecpssrKxsModel(    
 77 {                                                 
 78   protonDataSetMap.clear();                       
 79   alphaDataSetMap.clear();                        
 80   delete interpolation;                           
 81 }                                                 
 82                                                   
 83 //....oooOO0OOooo........oooOO0OOooo........oo    
 84                                                   
 85 G4double G4ANSTOecpssrKxsModel::CalculateCross    
 86 {                                                 
 87   G4Proton* aProton = G4Proton::Proton();         
 88   G4Alpha* aAlpha = G4Alpha::Alpha();             
 89   G4double sigma = 0;                             
 90                                                   
 91   if (massIncident == aProton->GetPDGMass())      
 92     {                                             
 93       if (energyIncident > 0.2*MeV && energyIn    
 94                                                   
 95           sigma = protonDataSetMap[zTarget]->F    
 96           if (sigma !=0 && energyIncident > pr    
 97         }                                         
 98     }                                             
 99                                                   
100   else if (massIncident == aAlpha->GetPDGMass(    
101         {                                         
102           if (energyIncident > 0.2*MeV && ener    
103                                                   
104           sigma = alphaDataSetMap[zTarget]->Fi    
105           if (sigma !=0 && energyIncident > al    
106           }                                       
107          }                                        
108                                                   
109    else                                           
110           {                                       
111             sigma = 0.;                           
112           }                                       
113                                                   
114   // sigma is in internal units: it has been c    
115   // the input file in barns bt the EmDataset     
116   return sigma;                                   
117 }                                                 
118