include/Geant4/G4ComponentAntiNuclNuclearXS.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 //  Calculation of the total, elastic and inelastic cross-sections
0028 //  of anti-nucleon and anti-nucleus interactions with nuclei
0029 //  based on Glauber approach,  V. Grishine formulaes for
0030 //  interpolations (ref. V.M.Grichine, Eur.Phys.J., C62(2009) 399;
0031 //  NIM, B267 (2009) 2460) and our parametrization of hadron-nucleon
0032 //  cross-sections
0033 //
0034 //
0035 //   Created by A.Galoyan and V. Uzhinsky, 18.11.2010
0036
0037
0038 #ifndef G4ComponentAntiNuclNuclearXS_h
0039 #define G4ComponentAntiNuclNuclearXS_h
0040
0041 #include <CLHEP/Units/PhysicalConstants.h>  // pi, fermi,..
0042
0043 #include "globals.hh"
0044 #include "G4Proton.hh"
0045 #include "G4AntiProton.hh"
0046 #include "G4AntiNeutron.hh"
0047 #include "G4AntiDeuteron.hh"
0048 #include "G4AntiHe3.hh"
0049 #include "G4AntiTriton.hh"
0050 #include "G4AntiAlpha.hh"
0051 #include "G4Nucleus.hh"
0052 #include "G4Pow.hh"
0053 #include "G4VComponentCrossSection.hh"
0054
0055
0056 class G4ParticleDefinition;
0057
0058
0059 class G4ComponentAntiNuclNuclearXS : public G4VComponentCrossSection {
0060
0061   public:
0062     G4ComponentAntiNuclNuclearXS ();
0063     virtual ~G4ComponentAntiNuclNuclearXS ();
0064     virtual G4double GetTotalIsotopeCrossSection(const G4ParticleDefinition* aParticle,
0065                                  G4double kinEnergy, G4int Z, G4int A);
0066     virtual G4double GetTotalElementCrossSection(const G4ParticleDefinition* aParticle,
0067                                  G4double kinEnergy, G4int Z, G4double A);
0068     virtual G4double GetInelasticIsotopeCrossSection(const G4ParticleDefinition* aParticle,
0069                                  G4double kinEnergy, G4int Z, G4int A);
0070     virtual G4double GetInelasticElementCrossSection(const G4ParticleDefinition* aParticle,
0071                                  G4double kinEnergy, G4int Z, G4double A);
0072     virtual G4double GetElasticElementCrossSection(const G4ParticleDefinition* aParticle,
0073                                G4double kinEnergy, G4int Z, G4double A);
0074     virtual G4double GetElasticIsotopeCrossSection(const G4ParticleDefinition* aParticle,
0075                                G4double kinEnergy, G4int Z, G4int A);
0076     virtual void BuildPhysicsTable(const G4ParticleDefinition&) {}
0077     virtual void DumpPhysicsTable(const G4ParticleDefinition&)  {}
0078     virtual void CrossSectionDescription(std::ostream&) const;
0079     // Method for calculation of Anti-Hadron Nucleon Total Cross-section
0080     G4double GetAntiHadronNucleonTotCrSc(const G4ParticleDefinition* aParticle, G4double kinEnergy);
0081     // Method for calculation of Anti-Hadron Nucleon Elastic Cross-section
0082     G4double GetAntiHadronNucleonElCrSc(const G4ParticleDefinition* aParticle, G4double kinEnergy);
0083
0084   private:
0085     G4double fRadiusEff;  // Effective Radius for AntiNucleus
0086     G4double fTotalXsc, fElasticXsc, fInelasticXsc;
0087     G4double fAntiHadronNucleonTotXsc, fAntiHadronNucleonElXsc;
0088     G4double Elab, S, SqrtS ;
0089     G4double Mn, b0, b2,  SqrtS0, S0, R0;  // Parameters for AntiHadron-Nucleon Xsc
0090     G4ParticleDefinition* theAProton;
0091     G4ParticleDefinition* theANeutron;
0092     G4ParticleDefinition* theADeuteron;
0093     G4ParticleDefinition* theATriton;
0094     G4ParticleDefinition* theAAlpha;
0095     G4ParticleDefinition* theAHe3;
0096
0097     const G4double ReffTot[5][5] =  { {0.000, 3.800, 3.300, 3.300, 2.376},    // Pbar   + p, d, t, He3, He4
0098                                       {3.800, 3.238, 3.144, 3.144, 2.544},    // Dbar   + p, d, t, He3, He4
0099                                       {3.300, 3.144, 3.075, 3.075, 3.589},    // Tbar   + p, d, t, He3, He4
0100                                       {3.300, 3.144, 3.075, 3.075, 2.589},    // He3bar + p, d, t, He3, He4
0101                                       {2.376, 2.544, 3.589, 3.598, 2.241} };  // He4bar + p, d, t, He3, He4
0102     const G4double ReffInel[5][5] = { {0.000, 3.582, 3.105, 3.105, 2.209},    // Pbar   + p, d, t, He3, He4
0103                                       {3.582, 3.169, 3.066, 3.066, 2.498},    // Dbar   + p, d, t, He3, He4
0104                                       {3.105, 3.066, 2.973, 2.973, 2.508},    // Tbar   + p, d, t, He3, He4
0105                                       {3.105, 3.066, 2.973, 2.973, 2.508},    // He3bar + p, d, t, He3, He4
0106                                       {2.209, 2.498, 2.508, 2.508, 2.158} };  // He4bar + p, d, t, He3, He4
0107
0108     const G4Pow* theG4Pow = G4Pow::GetInstance();
0109 };
0110
0111 #endif