Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/cross_sections/include/G4NeutronInelasticXS.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 ]

Diff markup

Differences between /processes/hadronic/cross_sections/include/G4NeutronInelasticXS.hh (Version 11.3.0) and /processes/hadronic/cross_sections/include/G4NeutronInelasticXS.hh (Version 10.7)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // -------------------------------------------     26 // -------------------------------------------------------------------
 27 //                                                 27 //
 28 // GEANT4 Class header file                        28 // GEANT4 Class header file
 29 //                                                 29 //
 30 //                                                 30 //
 31 // File name:    G4NeutronInelasticXS              31 // File name:    G4NeutronInelasticXS
 32 //                                                 32 //
 33 // Author  Ivantchenko, Geant4, 3-AUG-09           33 // Author  Ivantchenko, Geant4, 3-AUG-09
 34 //                                                 34 //
 35                                                    35 
 36 // Class Description:                              36 // Class Description:
 37 // This is a base class for neutron inelastic      37 // This is a base class for neutron inelastic hadronic cross section based on
 38 // data files from G4PARTICLEXSDATA data set       38 // data files from G4PARTICLEXSDATA data set 
 39 // Class Description - End                         39 // Class Description - End
 40                                                    40  
 41 #ifndef G4NeutronInelasticXS_h                     41 #ifndef G4NeutronInelasticXS_h
 42 #define G4NeutronInelasticXS_h 1                   42 #define G4NeutronInelasticXS_h 1
 43                                                    43 
 44 #include "G4VCrossSectionDataSet.hh"               44 #include "G4VCrossSectionDataSet.hh"
 45 #include "globals.hh"                              45 #include "globals.hh"
 46 #include "G4ElementData.hh"                        46 #include "G4ElementData.hh"
 47 #include "G4PhysicsVector.hh"                      47 #include "G4PhysicsVector.hh"
                                                   >>  48 #include "G4Threading.hh"
 48 #include <vector>                                  49 #include <vector>
 49                                                    50 
                                                   >>  51 const G4int MAXZINEL = 93;
                                                   >>  52 
 50 class G4DynamicParticle;                           53 class G4DynamicParticle;
 51 class G4ParticleDefinition;                        54 class G4ParticleDefinition;
 52 class G4Element;                                   55 class G4Element;
 53 class G4VComponentCrossSection;                    56 class G4VComponentCrossSection;
 54                                                    57 
 55 class G4NeutronInelasticXS final : public G4VC     58 class G4NeutronInelasticXS final : public G4VCrossSectionDataSet
 56 {                                                  59 {
 57 public:                                            60 public: 
 58                                                    61 
 59   G4NeutronInelasticXS();                      <<  62   explicit G4NeutronInelasticXS();
 60                                                    63 
 61   ~G4NeutronInelasticXS() override = default;  <<  64   ~G4NeutronInelasticXS() final;
 62                                                    65 
 63   static const char* Default_Name() { return " <<  66   static const char* Default_Name() {return "G4NeutronInelasticXS";}
 64                                                    67     
 65   G4bool IsElementApplicable(const G4DynamicPa     68   G4bool IsElementApplicable(const G4DynamicParticle*, G4int Z,
 66            const G4Material*) final;               69            const G4Material*) final;
 67                                                    70 
 68   G4bool IsIsoApplicable(const G4DynamicPartic     71   G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A,
 69        const G4Element*, const G4Material*) fi     72        const G4Element*, const G4Material*) final;
 70                                                    73 
 71   G4double GetElementCrossSection(const G4Dyna     74   G4double GetElementCrossSection(const G4DynamicParticle*, 
 72           G4int Z, const G4Material*) final;       75           G4int Z, const G4Material*) final;
 73                                                    76 
 74   G4double ComputeCrossSectionPerElement(G4dou << 
 75                                          const << 
 76                                          const << 
 77                                          const << 
 78                                                << 
 79   G4double ComputeIsoCrossSection(G4double kin << 
 80                                   const G4Part << 
 81                                   G4int Z, G4i << 
 82                                   const G4Isot << 
 83                                   const G4Elem << 
 84                                   const G4Mate << 
 85                                                << 
 86   G4double GetIsoCrossSection(const G4DynamicP     77   G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A,
 87                               const G4Isotope*     78                               const G4Isotope* iso,
 88                               const G4Element*     79                               const G4Element* elm,
 89                               const G4Material     80                               const G4Material* mat) final;
 90                                                    81 
 91   const G4Isotope* SelectIsotope(const G4Eleme     82   const G4Isotope* SelectIsotope(const G4Element*, 
 92                                  G4double kinE     83                                  G4double kinEnergy, G4double logE) final;
 93                                                    84 
 94   void BuildPhysicsTable(const G4ParticleDefin     85   void BuildPhysicsTable(const G4ParticleDefinition&) final;
 95                                                    86 
 96   void CrossSectionDescription(std::ostream&)      87   void CrossSectionDescription(std::ostream&) const final;
 97                                                    88 
 98   G4double ElementCrossSection(G4double kinEne << 
 99                                                << 
100   G4double IsoCrossSection(G4double ekin, G4do     89   G4double IsoCrossSection(G4double ekin, G4double logekin, G4int Z, G4int A);
101                                                    90 
102   G4NeutronInelasticXS & operator=(const G4Neu     91   G4NeutronInelasticXS & operator=(const G4NeutronInelasticXS &right) = delete;
103   G4NeutronInelasticXS(const G4NeutronInelasti     92   G4NeutronInelasticXS(const G4NeutronInelasticXS&) = delete;
104                                                    93 
105 private:                                           94 private: 
106                                                    95 
107   void Initialise(G4int Z);                        96   void Initialise(G4int Z);
108                                                    97 
109   void InitialiseOnFly(G4int Z);                   98   void InitialiseOnFly(G4int Z);
110                                                    99 
111   const G4String& FindDirectoryPath();            100   const G4String& FindDirectoryPath();
112                                                   101 
113   inline const G4PhysicsVector* GetPhysicsVect    102   inline const G4PhysicsVector* GetPhysicsVector(G4int Z);
114                                                   103 
115   G4PhysicsVector* RetrieveVector(std::ostring    104   G4PhysicsVector* RetrieveVector(std::ostringstream& in, G4bool warn);
116                                                   105   
117   G4VComponentCrossSection* ggXsection = nullp << 106   G4VComponentCrossSection* ggXsection;
118                                                   107 
119   const G4ParticleDefinition* neutron;            108   const G4ParticleDefinition* neutron;
120                                                   109 
121   std::vector<G4double> temp;                     110   std::vector<G4double> temp;
122                                                   111 
123   G4double elimit;                             << 112   G4bool  isMaster;
124   G4double lowElimit;                          << 
125   G4double loglowElimit;                       << 
126                                                << 
127   G4bool isInitializer{false};                 << 
128                                                   113 
129   static const G4int MAXZINEL = 93;            << 
130   static G4ElementData* data;                     114   static G4ElementData* data;
131   static G4double coeff[MAXZINEL];                115   static G4double coeff[MAXZINEL];
132   static G4double lowcoeff[MAXZINEL];          << 
133   static G4String gDataDirectory;                 116   static G4String gDataDirectory;
                                                   >> 117 
                                                   >> 118 #ifdef G4MULTITHREADED
                                                   >> 119   static G4Mutex neutronInelasticXSMutex;
                                                   >> 120 #endif
134 };                                                121 };
135                                                   122 
136 inline                                            123 inline
137 const G4PhysicsVector* G4NeutronInelasticXS::G    124 const G4PhysicsVector* G4NeutronInelasticXS::GetPhysicsVector(G4int Z)
138 {                                                 125 {
139   const G4PhysicsVector* pv = data->GetElement    126   const G4PhysicsVector* pv = data->GetElementData(Z);
140   if (pv == nullptr) {                         << 127   if(pv == nullptr) { 
141     InitialiseOnFly(Z);                           128     InitialiseOnFly(Z);
142     pv = data->GetElementData(Z);                 129     pv = data->GetElementData(Z);
143   }                                               130   }
144   return pv;                                      131   return pv;
145 }                                                 132 }
146                                                   133 
147 #endif                                            134 #endif
148                                                   135