coherent_elastic/include/G4AntiNuclElastic.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 //
0028 // Geant4 Header : G4AntiNuclElastic
0029 //
0030 // Author : A. Galoyan
0031 //
0032 // Class Description
0033 // Model for AntiNuclear Nuclear Elastic Scattering;
0034 // Class Description - End
0035
0036 #ifndef G4AntiNuclElastic_h
0037 #define G4AntiNuclElastic_h 1
0038
0039 #include "G4HadronElastic.hh"
0040 #include "globals.hh"
0041 #include "G4Nucleus.hh"
0042 #include "G4ComponentAntiNuclNuclearXS.hh"
0043
0044 class G4ParticleDefinition;
0045
0046 class G4AntiNuclElastic : public G4HadronElastic
0047 {
0048 public:
0049
0050   explicit G4AntiNuclElastic();
0051
0052   ~G4AntiNuclElastic() override;
0053
0054   G4double SampleInvariantT(const G4ParticleDefinition* p,
0055                             G4double plab, G4int Z, G4int A) override;
0056
0057  G4double SampleThetaCMS(const G4ParticleDefinition* p, G4double plab,
0058                          G4int Z, G4int A);
0059
0060  G4double SampleThetaLab(const G4ParticleDefinition* p,
0061                                 G4double plab, G4int Z, G4int A);
0062
0063
0064  G4double CalculateParticleBeta( const G4ParticleDefinition* particle,
0065                                     G4double momentum    );
0066
0067  G4double CalculateZommerfeld( G4double beta, G4double Z1, G4double Z2 );
0068
0069  G4double CalculateAm( G4double momentum, G4double n, G4double Z);
0070
0071  G4double DampFactor(G4double z);
0072
0073  G4double BesselJzero(G4double z);
0074
0075  G4double BesselJone(G4double z);
0076
0077  G4double BesselOneByArg(G4double z);
0078
0079  G4double GetcosTeta1( G4double plab, G4int A);
0080
0081  inline G4ComponentAntiNuclNuclearXS* GetComponentCrossSection();
0082
0083 private:
0084
0085   // Assignment operator and copy constructor
0086   G4AntiNuclElastic & operator=(const G4AntiNuclElastic &right);
0087   G4AntiNuclElastic(const G4AntiNuclElastic&);
0088
0089   G4ComponentAntiNuclNuclearXS* cs;    //cross section of antiA-A interaction
0090   const G4ParticleDefinition*  fParticle;
0091
0092   G4double fTetaCMS;        //  sampled Theta in CMS
0093   G4double fThetaLab;        //sampled Theta in Lab system
0094   G4double fWaveVector;
0095   G4double fBeta;         // velosity of projectile
0096   G4double fZommerfeld;   // parameter of Zommerfeld for calculation of Coulomb cross-section
0097   G4double fAm;           // parameter for calculation of Coulomb cross-section
0098   G4double fRa;           // Radius of target
0099   G4double fRef;          // Effective radiuse for Calculation of hadron cross-section
0100   G4double fceff;         //  Effective diffuse parameter
0101
0102   G4ThreeVector fbst;          // boost vector
0103   G4double fptot;         // momentum of projectile in CMS system
0104   G4double fTmax;
0105
0106   G4ParticleDefinition* theAProton;
0107   G4ParticleDefinition* theANeutron;
0108   G4ParticleDefinition* theADeuteron;
0109   G4ParticleDefinition* theATriton;
0110   G4ParticleDefinition* theAAlpha;
0111   G4ParticleDefinition* theAHe3;
0112
0113   G4ParticleDefinition* theProton;
0114   G4ParticleDefinition* theNeutron;
0115   G4ParticleDefinition* theDeuteron;
0116   G4ParticleDefinition* theAlpha;
0117
0118 };
0119
0120 inline G4ComponentAntiNuclNuclearXS*
0121 G4AntiNuclElastic::GetComponentCrossSection()
0122 {
0123   return cs;
0124 }
0125
0126 #endif
0127
0128