G4HadronNucleonXsc.hh

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 // Calculation of the total, elastic and inelastic cross-sections
00027 // based on parametrisations of (proton, pion, kaon, photon) nucleon
00028 // cross-sections and the hadron-nucleous cross-section model in 
00029 // the framework of Glauber-Gribov approach
00030 //
00031 //
00032 //
00033 // 14.03.07 V. Grichine - first implementation
00034 // 04.11.11 V. Grichine - update for kaon-(p,n) xsc, vector spline
00035 // 21.02.12 V. Grichine - update for pion-(p,n) xsc, NS fit++, vector spline
00036 //
00037 //
00038 
00039 #ifndef G4HadronNucleonXsc_h
00040 #define G4HadronNucleonXsc_h
00041 
00042 #include "globals.hh"
00043 #include "G4Proton.hh"
00044 #include "G4Nucleus.hh"
00045 #include "G4LPhysicsFreeVector.hh"
00046 
00047 
00048 class G4ParticleDefinition;
00049 class G4DynamicParticle;
00050 
00051 class G4HadronNucleonXsc 
00052 {
00053 public:
00054 
00055   G4HadronNucleonXsc ();
00056   virtual ~G4HadronNucleonXsc ();
00057    
00058   virtual
00059   G4bool IsApplicable(const G4DynamicParticle* aDP, const G4Element*);
00060 
00061   virtual
00062   G4bool IsIsoApplicable(const G4DynamicParticle* aDP, G4int Z, G4int A);
00063 
00064   virtual
00065   void DumpPhysicsTable(const G4ParticleDefinition&) 
00066   {G4cout << "G4HadronNucleonXsc: uses parametrisation"<<G4endl;}
00067 
00068   void CrossSectionDescription(std::ostream&) const;
00069 
00070   // Xsc parametrisations
00071 
00072   G4double GetHadronNucleonXscEL(const G4DynamicParticle*, const G4ParticleDefinition*);
00073 
00074   G4double GetHadronNucleonXscPDG(const G4DynamicParticle*, const G4ParticleDefinition*);
00075 
00076   G4double GetHadronNucleonXscNS(const G4DynamicParticle*, const G4ParticleDefinition*);
00077 
00078   G4double GetHadronNucleonXscVU(const G4DynamicParticle*, const G4ParticleDefinition*);
00079 
00080   // kinematics and set/get
00081 
00082   G4double CalculateEcmValue ( const G4double , const G4double , const G4double ); 
00083 
00084   G4double CalcMandelstamS( const G4double , const G4double , const G4double );
00085 
00086   G4double GetCoulombBarrier(const G4DynamicParticle* aParticle, const G4ParticleDefinition* nucleon );
00087 
00088   G4double GetTotalHadronNucleonXsc()    { return fTotalXsc;     }; 
00089   G4double GetElasticHadronNucleonXsc()  { return fElasticXsc;   }; 
00090   G4double GetInelasticHadronNucleonXsc(){ return fInelasticXsc; }; 
00091 
00092   void InitialiseKaonNucleonTotXsc();
00093 
00094   G4double GetKpProtonTotXscVector(G4double logEnergy){ return fKpProtonTotXscVector.Value(logEnergy); };
00095   G4double GetKpNeutronTotXscVector(G4double logEnergy){ return fKpNeutronTotXscVector.Value(logEnergy); };
00096   G4double GetKmProtonTotXscVector(G4double logEnergy){ return fKmProtonTotXscVector.Value(logEnergy); };
00097   G4double GetKmNeutronTotXscVector(G4double logEnergy){ return fKmNeutronTotXscVector.Value(logEnergy); };
00098 
00099 private:
00100 
00101   const G4double fUpperLimit;
00102   const G4double fLowerLimit; 
00103 
00104   G4double fTotalXsc, fElasticXsc, fInelasticXsc;
00105   G4double fHadronNucleonXsc;
00106 
00107   // K-nucleon tot xsc (mb) fit data, std::log(Tkin(MeV))
00108 
00109   static const G4double fKpProtonTotXsc[66];
00110   static const G4double fKpProtonTotTkin[66];
00111 
00112   static const G4double fKpNeutronTotXsc[75];
00113   static const G4double fKpNeutronTotTkin[75];
00114 
00115   static const G4double fKmProtonTotXsc[106];
00116   static const G4double fKmProtonTotTkin[106];
00117 
00118   static const G4double fKmNeutronTotXsc[68];
00119   static const G4double fKmNeutronTotTkin[68];
00120 
00121   G4LPhysicsFreeVector fKpProtonTotXscVector;
00122   G4LPhysicsFreeVector fKpNeutronTotXscVector;
00123   G4LPhysicsFreeVector fKmProtonTotXscVector;
00124   G4LPhysicsFreeVector fKmNeutronTotXscVector;
00125  
00126   G4ParticleDefinition* theGamma;
00127   G4ParticleDefinition* theProton;
00128   G4ParticleDefinition* theNeutron;
00129   G4ParticleDefinition* theAProton;
00130   G4ParticleDefinition* theANeutron;
00131   G4ParticleDefinition* thePiPlus;
00132   G4ParticleDefinition* thePiMinus;
00133   G4ParticleDefinition* thePiZero;
00134   G4ParticleDefinition* theKPlus;
00135   G4ParticleDefinition* theKMinus;
00136   G4ParticleDefinition* theK0S;
00137   G4ParticleDefinition* theK0L;
00138   G4ParticleDefinition* theL;
00139   G4ParticleDefinition* theAntiL;
00140   G4ParticleDefinition* theSPlus;
00141   G4ParticleDefinition* theASPlus;
00142   G4ParticleDefinition* theSMinus;
00143   G4ParticleDefinition* theASMinus;
00144   G4ParticleDefinition* theS0;
00145   G4ParticleDefinition* theAS0;
00146   G4ParticleDefinition* theXiMinus;
00147   G4ParticleDefinition* theXi0;
00148   G4ParticleDefinition* theAXiMinus;
00149   G4ParticleDefinition* theAXi0;
00150   G4ParticleDefinition* theOmega;
00151   G4ParticleDefinition* theAOmega;
00152   G4ParticleDefinition* theD;
00153   G4ParticleDefinition* theT;
00154   G4ParticleDefinition* theA;
00155   G4ParticleDefinition* theHe3;
00156 
00157 };
00158 
00159 
00160 #endif

Generated on Mon May 27 17:48:26 2013 for Geant4 by  doxygen 1.4.7