Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/cross_sections/include/G4GammaNuclearXS.hh

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 ]

  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 //
 27 // -------------------------------------------------------------------
 28 //
 29 // GEANT4 Class header file
 30 //
 31 //
 32 // File name:    G4GammaNuclearXS
 33 //
 34 // Authors  V.Ivantchenko, Geant4, 20 October 2020
 35 //          B.Kutsenko, BINP/NSU, 10 August 2021
 36 //
 37 //
 38 // Class Description:
 39 // This is a base class for gamma-nuclear cross section based on
 40 // data files from IAEA Evaluated Photonuclear Data Library (IAEA/PD-2019) 
 41 // https://www-nds.iaea.org/photonuclear/
 42 // Database review:
 43 // https://www.sciencedirect.com/science/article/pii/S0090375219300699?via%3Dihub
 44 //
 45 
 46 #ifndef G4GammaNuclearXS_h
 47 #define G4GammaNuclearXS_h 1
 48 
 49 #include "G4VCrossSectionDataSet.hh"
 50 #include "globals.hh"
 51 #include <vector>
 52 
 53 class G4DynamicParticle;
 54 class G4ParticleDefinition;
 55 class G4Element;
 56 class G4ElementData;
 57 class G4PhysicsVector;
 58 class G4PhotoNuclearCrossSection;
 59 
 60 class G4GammaNuclearXS final : public G4VCrossSectionDataSet
 61 {
 62 public: 
 63 
 64   G4GammaNuclearXS();
 65 
 66   ~G4GammaNuclearXS() override = default;
 67     
 68   static const char* Default_Name() { return "GammaNuclearXS"; }
 69 
 70   G4bool IsElementApplicable(const G4DynamicParticle*, 
 71            G4int Z, const G4Material*) final;
 72 
 73   G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A,
 74        const G4Element*, const G4Material* mat) final;
 75 
 76   G4double GetElementCrossSection(const G4DynamicParticle*, 
 77                 G4int Z, 
 78                                   const G4Material* mat = nullptr) final; 
 79 
 80   G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A,
 81                               const G4Isotope* iso = nullptr,
 82                               const G4Element* elm = nullptr,
 83                               const G4Material* mat = nullptr) final;
 84 
 85   const G4Isotope* SelectIsotope(const G4Element*, 
 86                                  G4double kinEnergy, G4double logE) final;
 87 
 88   void BuildPhysicsTable(const G4ParticleDefinition&) final;
 89 
 90   G4double IsoCrossSection(const G4double ekin, const G4int Z, const G4int A);
 91 
 92   G4double ElementCrossSection(const G4double ekin, const G4int Z);
 93 
 94   G4double LowEnergyCrossSection(G4double ekin, G4int Z);
 95 
 96   void CrossSectionDescription(std::ostream&) const final;
 97       
 98   G4GammaNuclearXS & operator=(const G4GammaNuclearXS &right) = delete;
 99   G4GammaNuclearXS(const G4GammaNuclearXS&) = delete;
100   
101 private: 
102 
103   void Initialise(G4int Z);
104 
105   const G4String& FindDirectoryPath();
106 
107   G4PhysicsVector* RetrieveVector(std::ostringstream& in, G4bool warn, G4int Z);
108   
109   G4PhotoNuclearCrossSection* ggXsection = nullptr;
110   const G4ParticleDefinition* gamma;
111 
112   // Cache
113   G4double fXS = 0.0;  
114   G4double fEkin = 0.0;  
115   G4int fZ = 0;  
116 
117   static const G4int MAXZGAMMAXS = 95;
118   static const G4int MAXNFREE = 11;
119   static G4ElementData* data;
120   // The list of elements with non-linear parametrisation for better precision 
121   static const G4int freeVectorException[MAXNFREE];
122   // CHIPS photonuclear model had a problem with high energy parametrisation 
123   // for isotopes of H and He, coefficient is needed to connect isotope cross 
124   // section with element cross-section on high energies.
125   static G4double coeff[3][3]; 
126   static G4double xs150[MAXZGAMMAXS];
127   static G4String gDataDirectory;
128 
129   std::vector<G4double> temp;
130 };
131 
132 #endif
133