G4ComponentAntiNuclNuclearXS Class Reference

#include <G4ComponentAntiNuclNuclearXS.hh>

Inheritance diagram for G4ComponentAntiNuclNuclearXS:

Public Member Functions
	G4ComponentAntiNuclNuclearXS ()
virtual	~G4ComponentAntiNuclNuclearXS ()
virtual G4double	GetTotalIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
virtual G4double	GetTotalElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
virtual G4double	GetInelasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
virtual G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
virtual G4double	GetElasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
virtual G4double	GetElasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
G4double	GetAntiHadronNucleonTotCrSc (const G4ParticleDefinition *aParticle, G4double kinEnergy)
G4double	GetAntiHadronNucleonElCrSc (const G4ParticleDefinition *aParticle, G4double kinEnergy)
Public Member Functions inherited from G4VComponentCrossSection
	G4VComponentCrossSection (const G4String &nam="")
virtual	~G4VComponentCrossSection ()
G4double	GetTotalElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
virtual G4double	GetTotalElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)=0
virtual G4double	GetTotalIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)=0
G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
virtual G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)=0
virtual G4double	GetInelasticIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)=0
G4double	GetElasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
virtual G4double	GetElasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)=0
virtual G4double	GetElasticIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)=0
virtual G4double	ComputeQuasiElasticRatio (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	Description () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Detailed Description

Definition at line 57 of file G4ComponentAntiNuclNuclearXS.hh.

Constructor & Destructor Documentation

G4ComponentAntiNuclNuclearXS()

G4ComponentAntiNuclNuclearXS::G4ComponentAntiNuclNuclearXS

(

)

Definition at line 48 of file G4ComponentAntiNuclNuclearXS.cc.

: G4VComponentCrossSection("AntiAGlauber"), fUpperLimit( 10000 * GeV ),
  fLowerLimit( 10 * MeV )
{
  theAProton       = G4AntiProton::AntiProton();
  theANeutron      = G4AntiNeutron::AntiNeutron();
  theADeuteron     = G4AntiDeuteron::AntiDeuteron();
  theATriton       = G4AntiTriton::AntiTriton();
  theAAlpha        = G4AntiAlpha::AntiAlpha();
  theAHe3          = G4AntiHe3::AntiHe3();
 
 Mn       = 0.93827231;           // GeV
 b0       = 11.92;                // GeV^(-2)
 b2       = 0.3036;               // GeV^(-2)
 SqrtS0   = 20.74;                // GeV
 S0       = 33.0625;              // GeV^2
 
}

~G4ComponentAntiNuclNuclearXS()

G4ComponentAntiNuclNuclearXS::~G4ComponentAntiNuclNuclearXS ( )

virtual

Definition at line 71 of file G4ComponentAntiNuclNuclearXS.cc.

{
}

Member Function Documentation

BuildPhysicsTable()

virtual void G4ComponentAntiNuclNuclearXS::BuildPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 95 of file G4ComponentAntiNuclNuclearXS.hh.

 {}

CrossSectionDescription()

void G4ComponentAntiNuclNuclearXS::CrossSectionDescription ( std::ostream &  outFile ) const

virtual

Definition at line 78 of file G4ComponentAntiNuclNuclearXS.cc.

{
  outFile << "G4ComponentAntiNuclNuclearXS describes the total, elastic\n"
          << "and inelastic cross sections for the scattering of light\n"
          << "antinuclei from nuclei using the Glauber approach and Grichine\n"
          << "parameterization.  It is valid from 10 MeV to 10 TeV incident\n"
          << "energies and applies to antiprotons, antineutrons, antideuterons,\n"
          << "anti3He, antitritons and antialphas.  This is a cross section\n"
          << "component which is to be used as part of a cross section data\n"
          << "set.\n";
} 

DumpPhysicsTable()

virtual void G4ComponentAntiNuclNuclearXS::DumpPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 99 of file G4ComponentAntiNuclNuclearXS.hh.

 {}

GetAntiHadronNucleonElCrSc()

G4double G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonElCrSc	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy
	)

Definition at line 335 of file G4ComponentAntiNuclNuclearXS.cc.

{
 G4double xsection;
 
 G4double   SigAss;
 G4double   C, d1, d2, d3  ;
 
 GetAntiHadronNucleonTotCrSc(aParticle,kinEnergy);
 
 SigAss   = 4.5 + 0.101*std::log(S/S0)*std::log(S/S0);            //mb
  
 C        = 59.27;
 d1       = -6.95;
 d2       = 23.54;
 d3       = -25.34;
 
 xsection = SigAss* (1 + 1. / (std::sqrt(S-4.*Mn*Mn)) / (std::pow(R0, 3.))
  *C* ( 1+d1/SqrtS+d2/(std::pow(SqrtS,2.))+d3/(std::pow(SqrtS,3.)) ));  
 
//  xsection *= millibarn;
 
 fAntiHadronNucleonElXsc = xsection;
 return fAntiHadronNucleonElXsc;
}

Referenced by GetInelasticElementCrossSection(), and GetTotalElementCrossSection().

GetAntiHadronNucleonTotCrSc()

G4double G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy
	)

Definition at line 294 of file G4ComponentAntiNuclNuclearXS.cc.

{
  G4double xsection, Pmass, Energy, momentum;
  const G4ParticleDefinition* theParticle = aParticle;
  Pmass=theParticle->GetPDGMass();
  Energy=Pmass+kinEnergy;
  momentum=std::sqrt(Energy*Energy-Pmass*Pmass)/std::abs(theParticle->GetBaryonNumber());
  G4double Plab = momentum / GeV/std::abs(aParticle->GetBaryonNumber()); // Uzhi 24 Nov. 2011
 
  if(Plab < 0.1) { Plab = 0.1; }    // Uzhi 24 Nov. 2011
 
  G4double   B, SigAss;
  G4double   C, d1, d2, d3  ;
 
  Elab     = std::sqrt(Mn*Mn + Plab*Plab);   // GeV
  S        = 2.*Mn*Mn + 2. *Mn*Elab;         // GeV^2
  SqrtS    = std::sqrt(S);                   // GeV 
 
  B        = b0+b2*std::log(SqrtS/SqrtS0)*std::log(SqrtS/SqrtS0); //GeV^(-2)
  SigAss   = 36.04 +0.304*std::log(S/S0)*std::log(S/S0);          //mb 
  R0       = std::sqrt(0.40874044*SigAss - B);                   //GeV^(-2)
 
  C        = 13.55;
  d1       = -4.47;
  d2       = 12.38;
  d3       = -12.43;
  xsection = SigAss*(1 + 1./(std::sqrt(S-4.*Mn*Mn)) / (std::pow(R0, 3.))
  *C* (1+d1/SqrtS+d2/(std::pow(SqrtS,2.))+d3/(std::pow(SqrtS,3.)) ));  
 
//  xsection *= millibarn;
 
  fAntiHadronNucleonTotXsc = xsection;
  return fAntiHadronNucleonTotXsc;
}

Referenced by GetInelasticElementCrossSection(), GetTotalElementCrossSection(), and G4AntiNuclElastic::SampleInvariantT().

GetElasticElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetElasticElementCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4double  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 271 of file G4ComponentAntiNuclNuclearXS.cc.

{
 fElasticXsc = GetTotalElementCrossSection(aParticle, kinEnergy, Z, A)-
   GetInelasticElementCrossSection(aParticle, kinEnergy, Z, A);
 
 if (fElasticXsc < 0.) fElasticXsc = 0.;
 
 return fElasticXsc;
}

Referenced by GetElasticIsotopeCrossSection(), and G4AntiNuclElastic::SampleInvariantT().

GetElasticIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetElasticIsotopeCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4int  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 286 of file G4ComponentAntiNuclNuclearXS.cc.

{ return GetElasticElementCrossSection(aParticle, kinEnergy, Z, (G4double) A); }

GetInelasticElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4double  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 181 of file G4ComponentAntiNuclNuclearXS.cc.

{
  G4double  inelxsection,  sigmaTotal, sigmaElastic;
 
  const G4ParticleDefinition* theParticle = aParticle;
 
  sigmaTotal        = GetAntiHadronNucleonTotCrSc(theParticle,kinEnergy);
  sigmaElastic      = GetAntiHadronNucleonElCrSc(theParticle,kinEnergy);
  
// calculation of sqr of radius NN-collision
  fRadiusNN2=sigmaTotal*sigmaTotal*0.1/(8.*sigmaElastic*pi);   // fm^2   
 
 
// calculation of effective nuclear radius for Pbar and Nbar interaction (can be changed)
 
  if ( (theParticle == theAProton) || (theParticle == theANeutron) )
  {
   if (A==1)
      { fInelasticXsc = (sigmaTotal - sigmaElastic) * millibarn;
        return fInelasticXsc;  
      } 
   fRadiusEff = 1.31*std::pow(A, 0.22)+0.9/std::pow(A, 1./3.);  //fm
    
   if( (Z==1) && (A==2) ) fRadiusEff = 3.582;               //fm
   if( (Z==1) && (A==3) ) fRadiusEff = 3.105;               
   if( (Z==2) && (A==3) ) fRadiusEff = 3.105;
   if( (Z==2) && (A==4) ) fRadiusEff = 2.209;
  }
 
//calculation of effective nuclear radius for AntiDeuteron interaction (can be changed)
 
  if (theParticle ==theADeuteron) 
  { 
    fRadiusEff = 1.38*std::pow(A, 0.21)+1.55/std::pow(A, 1./3.);
  
    if( (Z==1) && (A==2) ) fRadiusEff = 3.169;            //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 3.066;
    if( (Z==2) && (A==3) ) fRadiusEff = 3.066;
    if( (Z==2) && (A==4) ) fRadiusEff = 2.498;
  }
 
//calculation of effective nuclear radius for AntiHe3 interaction (can be changed)
 
  if( (theParticle ==theAHe3) || (theParticle ==theATriton) )
  {
    fRadiusEff = 1.34 * std::pow(A, 0.21)+1.51/std::pow(A, 1./3.);
  
    if( (Z==1) && (A==2) ) fRadiusEff = 3.066;           //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 2.973;
    if( (Z==2) && (A==3) ) fRadiusEff = 2.973;
    if( (Z==2) && (A==4) ) fRadiusEff = 2.508;
  
  }
 
//calculation of effective nuclear radius for AntiAlpha interaction (can be changed)
 
  if (theParticle == theAAlpha) 
  {
    fRadiusEff = 1.3*std::pow(A, 0.21)+1.05/std::pow(A, 1./3.);
    
    if( (Z==1) && (A==2) ) fRadiusEff = 2.498;            //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 2.508;
    if( (Z==2) && (A==3) ) fRadiusEff = 2.508;
    if( (Z==2) && (A==4) ) fRadiusEff = 2.158;
  }
  G4double R2 = fRadiusEff*fRadiusEff;
  G4double REf2  = R2+fRadiusNN2;
  G4double  ApAt= std::abs(theParticle->GetBaryonNumber())  *  A;
 
  inelxsection  = pi*REf2 *10* std::log(1+(ApAt*sigmaTotal/(pi*REf2*10.))); //mb
  inelxsection  = inelxsection * millibarn;  
  fInelasticXsc =  inelxsection; 
  return fInelasticXsc;
}

Referenced by GetElasticElementCrossSection(), and GetInelasticIsotopeCrossSection().

GetInelasticIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetInelasticIsotopeCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4int  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 261 of file G4ComponentAntiNuclNuclearXS.cc.

{return GetInelasticElementCrossSection(aParticle, kinEnergy, Z, (G4double) A); }

GetTotalElementCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4double  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 95 of file G4ComponentAntiNuclNuclearXS.cc.

{
  G4double xsection,   sigmaTotal, sigmaElastic;
 
 const G4ParticleDefinition* theParticle = aParticle;
 
  sigmaTotal        = GetAntiHadronNucleonTotCrSc(theParticle,kinEnergy);
  sigmaElastic      = GetAntiHadronNucleonElCrSc(theParticle,kinEnergy);
 
// calculation of squared radius of  NN-collision
  fRadiusNN2=sigmaTotal*sigmaTotal*0.1/(8.*sigmaElastic*pi) ;  //fm^2   
 
// calculation of effective nuclear radius for Pbar and Nbar interactions (can be changed)
 
  if ( (theParticle == theAProton) || (theParticle == theANeutron) )
  {   
     if(A==1)
     { fTotalXsc = sigmaTotal * millibarn;
        return fTotalXsc;  }
 
     fRadiusEff = 1.34*std::pow(A,0.23)+1.35/std::pow(A,1./3.);   //fm
  
     if( (Z==1) && (A==2) ) fRadiusEff = 3.800;     //fm
     if( (Z==1) && (A==3) ) fRadiusEff = 3.300;  
     if( (Z==2) && (A==3) ) fRadiusEff = 3.300;  
     if( (Z==2) && (A==4) ) fRadiusEff = 2.376;     
  }
      
//calculation of effective nuclear radius for AntiDeuteron interaction (can be changed)
  if (theParticle == theADeuteron) 
  { fRadiusEff = 1.46 * std::pow(A,0.21) + 1.45 / std::pow(A,1./3.);
 
    if( (Z==1) && (A==2) ) fRadiusEff = 3.238;     //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 3.144;     
    if( (Z==2) && (A==3) ) fRadiusEff = 3.144;      
    if( (Z==2) && (A==4) ) fRadiusEff = 2.544;     
  }
// calculation of effective nuclear radius for AntiHe3 interaction (can be changed)
 
  if( (theParticle ==theAHe3) || (theParticle ==theATriton) )
  { fRadiusEff = 1.40* std::pow(A,0.21)+1.63/std::pow(A,1./3.);
 
    if( (Z==1) && (A==2) ) fRadiusEff = 3.144;     //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 3.075;  
    if( (Z==2) && (A==3) ) fRadiusEff = 3.075;  
    if( (Z==2) && (A==4) ) fRadiusEff = 2.589;  
  }
 
//calculation of effective nuclear radius for AntiAlpha interaction (can be changed)
 
  if (theParticle == theAAlpha) 
  {
    fRadiusEff = 1.35* std::pow(A,0.21)+1.1/std::pow(A,1./3.);
  
    if( (Z==1) && (A==2) ) fRadiusEff = 2.544;     //fm
    if( (Z==1) && (A==3) ) fRadiusEff = 2.589;   
    if( (Z==2) && (A==3) ) fRadiusEff = 2.589;   
    if( (Z==2) && (A==4) ) fRadiusEff = 2.241;    
  
  }
 
  G4double R2 = fRadiusEff*fRadiusEff;
  G4double REf2  = R2+fRadiusNN2;
  G4double ApAt = std::abs(theParticle->GetBaryonNumber())  *  A;
 
  xsection = 2*pi*REf2*10.*std::log(1+(ApAt*sigmaTotal/(2*pi*REf2*10.)));  //mb
  xsection =xsection *millibarn; 
  fTotalXsc   = xsection;
 
  return fTotalXsc; 
}

Referenced by GetElasticElementCrossSection(), GetTotalIsotopeCrossSection(), and G4AntiNuclElastic::SampleInvariantT().

GetTotalIsotopeCrossSection()

G4double G4ComponentAntiNuclNuclearXS::GetTotalIsotopeCrossSection ( const G4ParticleDefinition *  aParticle, G4double  kinEnergy, G4int  Z, G4int  A  )

virtual

Implements G4VComponentCrossSection.

Definition at line 173 of file G4ComponentAntiNuclNuclearXS.cc.

{ return GetTotalElementCrossSection(aParticle, kinEnergy, Z, (G4double) A);  }

---

The documentation for this class was generated from the following files:

• G4ComponentAntiNuclNuclearXS.hh
• G4ComponentAntiNuclNuclearXS.cc