G4GGNuclNuclCrossSection.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 nucleus-nucleus total, inelastic, production, 
00027 // elastic and quasi-elastic  cross-sections
00028 // based on parametrisations of nucleon-nucleon
00029 // cross-sections  in 
00030 // the framework of simplified Glauber-Gribov approach
00031 //
00032 //
00033 // 24.11.08 V. Grichine - first implementation based on G4GlauberGribovCrossSection
00034 //
00035 //
00036 
00037 #ifndef G4GGNuclNuclCrossSection_h
00038 #define G4GGNuclNuclCrossSection_h
00039 
00040 #include "globals.hh"
00041 #include "G4Proton.hh"
00042 #include "G4Nucleus.hh"
00043 #include "G4NistManager.hh"
00044 
00045 #include "G4VCrossSectionDataSet.hh"
00046 
00047 class G4ParticleDefinition;
00048 class G4HadronNucleonXsc;
00049 
00050 class G4GGNuclNuclCrossSection : public G4VCrossSectionDataSet
00051 {
00052 public:
00053 
00054   G4GGNuclNuclCrossSection ();
00055   virtual ~G4GGNuclNuclCrossSection ();
00056     
00057   static const char* Default_Name() {return "Glauber-Gribov nucleus nucleus";}
00058 
00059    
00060   virtual
00061   G4bool IsElementApplicable(const G4DynamicParticle*, 
00062                              G4int Z, const G4Material*);
00063 
00064   virtual
00065   G4double GetElementCrossSection(const G4DynamicParticle*, 
00066                                   G4int Z, const G4Material*);
00067 
00068   G4double GetZandACrossSection(const G4DynamicParticle*, 
00069                                 G4int Z, G4int A);
00070 
00071   G4double GetCoulombBarier(const G4DynamicParticle*, 
00072                             G4double Z, G4double A, G4double pR, G4double tR);
00073 
00074   virtual
00075   void BuildPhysicsTable(const G4ParticleDefinition&)
00076   {}
00077 
00078   virtual
00079   void DumpPhysicsTable(const G4ParticleDefinition&) 
00080   {G4cout << "G4NuclNuclCrossSection: uses Glauber-Gribov formula"<<G4endl;}
00081 
00082   virtual void CrossSectionDescription(std::ostream&) const;
00083 
00084   G4double GetRatioSD(const G4DynamicParticle*, G4double At, G4double Zt);
00085   G4double GetRatioQE(const G4DynamicParticle*, G4double At, G4double Zt);
00086 
00087   G4double GetHadronNucleonXsc(const G4DynamicParticle*, const G4Element*);
00088   G4double GetHadronNucleonXsc(const G4DynamicParticle*, G4int At, G4int Zt);
00089 
00090  
00091   G4double GetHadronNucleonXscPDG(G4ParticleDefinition*,G4double sMand, G4ParticleDefinition*);
00092   G4double GetHadronNucleonXscNS(G4ParticleDefinition*,G4double pTkin, G4ParticleDefinition*);
00093 
00094   G4double GetHNinelasticXscVU(const G4DynamicParticle*, G4int At, G4int Zt);
00095   G4double CalculateEcmValue(const G4double, const G4double, const G4double); 
00096   G4double CalcMandelstamS( const G4double , const G4double , const G4double );
00097 
00098   G4double GetElasticGlauberGribov(const G4DynamicParticle*,G4int Z, G4int A);
00099   G4double GetInelasticGlauberGribov(const G4DynamicParticle*,G4int Z, G4int A);
00100 
00101   G4double GetTotalGlauberGribovXsc()    { return fTotalXsc;     }; 
00102   G4double GetElasticGlauberGribovXsc()  { return fElasticXsc;   }; 
00103   G4double GetInelasticGlauberGribovXsc(){ return fInelasticXsc; }; 
00104   G4double GetProductionGlauberGribovXsc(){ return fProductionXsc; }; 
00105   G4double GetDiffractionGlauberGribovXsc(){ return fDiffractionXsc; }; 
00106   G4double GetRadiusConst()              { return fRadiusConst;  }; 
00107 
00108   G4double GetNucleusRadius(const G4DynamicParticle*, const G4Element*);
00109 
00110   G4double GetNucleusRadius(G4double Zt, G4double At);
00111   G4double GetNucleusRadiusGG(G4double At);
00112   G4double GetNucleusRadiusDE(G4double Z, G4double A);
00113   G4double GetNucleusRadiusRMS(G4double Z, G4double A);
00114 
00115   inline void SetEnergyLowerLimit(G4double E ){fLowerLimit=E;};
00116 
00117 private:
00118 
00119   const G4double fUpperLimit;
00120   G4double fLowerLimit; 
00121   const G4double fRadiusConst;
00122  
00123   G4double fTotalXsc, fElasticXsc, fInelasticXsc, fProductionXsc, fDiffractionXsc;
00124   G4double fHadronNucleonXsc;
00125  
00126   G4ParticleDefinition* theProton;
00127   G4ParticleDefinition* theNeutron;
00128   G4HadronNucleonXsc* hnXsc;
00129 };
00130 
00132 //
00133 // Inlines
00134 
00135 inline G4double
00136 G4GGNuclNuclCrossSection::GetElasticGlauberGribov(const G4DynamicParticle* dp,
00137                                                   G4int Z, G4int A)
00138 {
00139   GetZandACrossSection(dp, Z, A);
00140   return fElasticXsc;
00141 }
00142 
00144 
00145 inline G4double
00146 G4GGNuclNuclCrossSection::GetInelasticGlauberGribov(const G4DynamicParticle* dp,
00147                                                     G4int Z, G4int A)
00148 {
00149   GetZandACrossSection(dp, Z, A);
00150   return fInelasticXsc;
00151 }
00152 
00153 #endif

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