G4AntiNuclElastic.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 //
00027 // $Id: G4AntiNuclElastic.hh, 2011-04-12 agaloyan
00028 // GEANT4 tag $Name: not supported by cvs2svn $
00029 //
00030 // Geant4 Header : G4AntiNuclElastic
00031 //
00032 // Author : A. Galoyan  
00033 //  
00034 // Class Description
00035 // Model for AntiNuclear Nuclear Elastic Scattering;
00036 // Class Description - End
00037 
00038 #ifndef G4AntiNuclElastic_h
00039 #define G4AntiNuclElastic_h 1
00040  
00041 #include "G4HadronElastic.hh"
00042 #include "globals.hh"
00043 #include "G4Nucleus.hh"
00044 #include "G4ComponentAntiNuclNuclearXS.hh"
00045 
00046 using namespace std;
00047         
00048 class G4ParticleDefinition;
00049 
00050 class G4AntiNuclElastic : public G4HadronElastic
00051 {
00052 public:
00053 
00054   G4AntiNuclElastic();
00055 
00056   virtual ~G4AntiNuclElastic();
00057  
00058   virtual G4double SampleInvariantT(const G4ParticleDefinition* p, 
00059                                     G4double plab,
00060                                     G4int Z, G4int A);
00061 
00062  G4double SampleThetaCMS(const G4ParticleDefinition* p, G4double plab, 
00063  G4int Z, G4int A);
00064 
00065  G4double SampleThetaLab(const G4ParticleDefinition* p, 
00066                                 G4double plab, G4int Z, G4int A);
00067 
00068  
00069  G4double CalculateParticleBeta( const G4ParticleDefinition* particle, 
00070                                         G4double momentum    );
00071  
00072  G4double CalculateZommerfeld( G4double beta, G4double Z1, G4double Z2 );
00073 
00074  G4double CalculateAm( G4double momentum, G4double n, G4double Z);
00075  
00076  G4double DampFactor(G4double z);
00077 
00078  G4double BesselJzero(G4double z);
00079 
00080  G4double BesselJone(G4double z);
00081 
00082  G4double BesselOneByArg(G4double z); 
00083 
00084  G4double GetcosTeta1( G4double plab, G4int A);
00085 
00086  inline G4ComponentAntiNuclNuclearXS* GetComponentCrossSection();
00087  
00088 private:
00089 
00090   // Assignment operator and copy constructor
00091   G4AntiNuclElastic & operator=(const G4AntiNuclElastic &right);
00092   G4AntiNuclElastic(const G4AntiNuclElastic&);
00093 
00094   G4ComponentAntiNuclNuclearXS* cs;    //cross section of antiA-A interaction
00095 
00096   G4double fTetaCMS;        //  sampled Theta in CMS 
00097   G4double fThetaLab;        //sampled Theta in Lab system
00098   const G4ParticleDefinition*  fParticle;
00099   G4double fWaveVector;
00100   G4double fBeta;         // velosity of projectile 
00101   G4double fZommerfeld;   // parameter of Zommerfeld for calculation of Coulomb cross-section
00102   G4double fAm;           // parameter for calculation of Coulomb cross-section
00103   G4double fRa;           // Radius of target
00104   G4double fRef;          // Effective radiuse for Calculation of hadron cross-section
00105   G4double fceff;         //  Effective diffuse parameter
00106 
00107   G4ThreeVector fbst;          // boost vector
00108   G4double fptot;         // momentum of projectile in CMS system
00109   G4double fTmax;          
00110 
00111   G4ParticleDefinition* theAProton;
00112   G4ParticleDefinition* theANeutron;
00113   G4ParticleDefinition* theADeuteron;
00114   G4ParticleDefinition* theATriton;
00115   G4ParticleDefinition* theAAlpha;
00116   G4ParticleDefinition* theAHe3;
00117  
00118   G4ParticleDefinition* theProton;
00119   G4ParticleDefinition* theNeutron;
00120   G4ParticleDefinition* theDeuteron;
00121   G4ParticleDefinition* theAlpha;
00122     
00123 };
00124 
00125 inline G4ComponentAntiNuclNuclearXS* 
00126 G4AntiNuclElastic::GetComponentCrossSection()
00127 {
00128   return cs;
00129 }
00130 
00131 #endif
00132 
00133 

Generated on Mon May 27 17:47:40 2013 for Geant4 by  doxygen 1.4.7