G4GGNuclNuclCrossSection Class Reference

#include <G4GGNuclNuclCrossSection.hh>

Inheritance diagram for G4GGNuclNuclCrossSection:

Public Member Functions
	G4GGNuclNuclCrossSection ()
virtual	~G4GGNuclNuclCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
G4double	GetZandACrossSection (const G4DynamicParticle *, G4int Z, G4int A)
G4double	GetCoulombBarier (const G4DynamicParticle *, G4double Z, G4double A, G4double pR, G4double tR)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
G4double	GetRatioSD (const G4DynamicParticle *, G4double At, G4double Zt)
G4double	GetRatioQE (const G4DynamicParticle *, G4double At, G4double Zt)
G4double	GetHadronNucleonXsc (const G4DynamicParticle *, const G4Element *)
G4double	GetHadronNucleonXsc (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHadronNucleonXscPDG (G4ParticleDefinition *, G4double sMand, G4ParticleDefinition *)
G4double	GetHadronNucleonXscNS (G4ParticleDefinition *, G4double pTkin, G4ParticleDefinition *)
G4double	GetHNinelasticXscVU (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	CalculateEcmValue (const G4double, const G4double, const G4double)
G4double	CalcMandelstamS (const G4double, const G4double, const G4double)
G4double	GetElasticGlauberGribov (const G4DynamicParticle *, G4int Z, G4int A)
G4double	GetInelasticGlauberGribov (const G4DynamicParticle *, G4int Z, G4int A)
G4double	GetTotalGlauberGribovXsc ()
G4double	GetElasticGlauberGribovXsc ()
G4double	GetInelasticGlauberGribovXsc ()
G4double	GetProductionGlauberGribovXsc ()
G4double	GetDiffractionGlauberGribovXsc ()
G4double	GetRadiusConst ()
G4double	GetNucleusRadius (const G4DynamicParticle *, const G4Element *)
G4double	GetNucleusRadius (G4double Zt, G4double At)
G4double	GetNucleusRadiusGG (G4double At)
G4double	GetNucleusRadiusDE (G4double Z, G4double A)
G4double	GetNucleusRadiusRMS (G4double Z, G4double A)
void	SetEnergyLowerLimit (G4double E)
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 50 of file G4GGNuclNuclCrossSection.hh.

Constructor & Destructor Documentation

G4GGNuclNuclCrossSection()

G4GGNuclNuclCrossSection::G4GGNuclNuclCrossSection

(

)

Definition at line 46 of file G4GGNuclNuclCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()),
   fUpperLimit(100000*GeV), fLowerLimit(0.1*MeV),
   fRadiusConst(1.08*fermi),  // 1.1, 1.3 ?
   fTotalXsc(0.0), fElasticXsc(0.0), fInelasticXsc(0.0), fProductionXsc(0.0),
   fDiffractionXsc(0.0), fHadronNucleonXsc(0.0)
{
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  hnXsc = new G4HadronNucleonXsc();
}

~G4GGNuclNuclCrossSection()

G4GGNuclNuclCrossSection::~G4GGNuclNuclCrossSection ( )

virtual

Definition at line 59 of file G4GGNuclNuclCrossSection.cc.

{}

Member Function Documentation

BuildPhysicsTable()

virtual void G4GGNuclNuclCrossSection::BuildPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 75 of file G4GGNuclNuclCrossSection.hh.

  {}

CalcMandelstamS()

G4double G4GGNuclNuclCrossSection::CalcMandelstamS	(	const G4double	mp,
		const G4double	mt,
		const G4double	Plab
	)

Definition at line 766 of file G4GGNuclNuclCrossSection.cc.

{
  G4double Elab = std::sqrt ( mp * mp + Plab * Plab );
  G4double sMand  = mp*mp + mt*mt + 2*Elab*mt ;
 
  return sMand;
}

Referenced by GetHadronNucleonXsc(), and GetHadronNucleonXscNS().

CalculateEcmValue()

G4double G4GGNuclNuclCrossSection::CalculateEcmValue	(	const G4double	mp,
		const G4double	mt,
		const G4double	Plab
	)

Definition at line 749 of file G4GGNuclNuclCrossSection.cc.

{
  G4double Elab = std::sqrt ( mp * mp + Plab * Plab );
  G4double Ecm  = std::sqrt ( mp * mp + mt * mt + 2 * Elab * mt );
  // G4double Pcm  = Plab * mt / Ecm;
  // G4double KEcm = std::sqrt ( Pcm * Pcm + mp * mp ) - mp;
 
  return Ecm ; // KEcm;
}

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 63 of file G4GGNuclNuclCrossSection.cc.

{
  outFile << "G4GGNuclNuclCrossSection calculates total, inelastic and\n"
          << "elastic cross sections for nucleus-nucleus collisions using\n"
          << "the Glauber model with Gribov corrections.  It is valid for\n"
          << "all incident energies above 100 keV./n";
}

Default_Name()

static const char * G4GGNuclNuclCrossSection::Default_Name ( )

inlinestatic

Definition at line 57 of file G4GGNuclNuclCrossSection.hh.

{return "Glauber-Gribov nucleus nucleus";}

Referenced by G4IonPhysics::ConstructProcess().

DumpPhysicsTable()

virtual void G4GGNuclNuclCrossSection::DumpPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 79 of file G4GGNuclNuclCrossSection.hh.

  {G4cout << "G4NuclNuclCrossSection: uses Glauber-Gribov formula"<<G4endl;}

GetCoulombBarier()

G4double G4GGNuclNuclCrossSection::GetCoulombBarier	(	const G4DynamicParticle *	aParticle,
		G4double	Z,
		G4double	A,
		G4double	pR,
		G4double	tR
	)

Definition at line 202 of file G4GGNuclNuclCrossSection.cc.

{
  G4double ratio;
  G4double pZ = aParticle->GetDefinition()->GetPDGCharge();
 
  G4double pTkin = aParticle->GetKineticEnergy();
  // G4double pPlab = aParticle->GetTotalMomentum();
  G4double pM    = aParticle->GetDefinition()->GetPDGMass();
  // G4double tM    = tZ*proton_mass_c2 + (tA-tZ)*neutron_mass_c2; // ~ 1% accuracy
  G4double tM    = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIonMass( G4int(tZ), G4int(tA) );
  G4double pElab = pTkin + pM;
  G4double totEcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM);
  // G4double pPcm  = pPlab*tM/totEcm;
  // G4double pTcm  = std::sqrt(pM*pM + pPcm*pPcm) - pM;
  G4double totTcm  = totEcm - pM -tM;
 
  G4double bC    = fine_structure_const*hbarc*pZ*tZ;
           bC   /= pR + tR;
           bC   /= 2.;  // 4., 2. parametrisation cof ??? vmg
 
       // G4cout<<"pTkin = "<<pTkin/GeV<<"; pPlab = "
       // <<pPlab/GeV<<"; bC = "<<bC/GeV<<"; pTcm = "<<pTcm/GeV<<G4endl;
 
  if( totTcm <= bC ) ratio = 0.;
  else             ratio = 1. - bC/totTcm;
 
  // if(ratio < DBL_MIN) ratio = DBL_MIN;
  if( ratio < 0.) ratio = 0.;
 
  // G4cout <<"ratio = "<<ratio<<G4endl;
  return ratio;
}

Referenced by GetZandACrossSection().

GetDiffractionGlauberGribovXsc()

G4double G4GGNuclNuclCrossSection::GetDiffractionGlauberGribovXsc ( )

inline

Definition at line 105 of file G4GGNuclNuclCrossSection.hh.

{ return fDiffractionXsc; }; 

GetElasticGlauberGribov()

G4double G4GGNuclNuclCrossSection::GetElasticGlauberGribov ( const G4DynamicParticle *  dp, G4int  Z, G4int  A  )

inline

Definition at line 136 of file G4GGNuclNuclCrossSection.hh.

{
  GetZandACrossSection(dp, Z, A);
  return fElasticXsc;
}

GetElasticGlauberGribovXsc()

G4double G4GGNuclNuclCrossSection::GetElasticGlauberGribovXsc ( )

inline

Definition at line 102 of file G4GGNuclNuclCrossSection.hh.

{ return fElasticXsc;   }; 

GetElementCrossSection()

G4double G4GGNuclNuclCrossSection::GetElementCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 90 of file G4GGNuclNuclCrossSection.cc.

{
  G4int A = G4lrint(G4NistManager::Instance()->GetAtomicMassAmu(Z));
  return GetZandACrossSection(aParticle, Z, A);
}

GetHadronNucleonXsc() [1/2]

G4double G4GGNuclNuclCrossSection::GetHadronNucleonXsc	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement
	)

Definition at line 328 of file G4GGNuclNuclCrossSection.cc.

{
  G4int At = G4lrint(anElement->GetN());  // number of nucleons 
  G4int Zt = G4lrint(anElement->GetZ());  // number of protons
  return GetHadronNucleonXsc(aParticle, At, Zt);
}

Referenced by GetHadronNucleonXsc().

GetHadronNucleonXsc() [2/2]

G4double G4GGNuclNuclCrossSection::GetHadronNucleonXsc	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt
	)

Definition at line 347 of file G4GGNuclNuclCrossSection.cc.

{
  G4double xsection = 0.;
 
  G4double targ_mass = G4ParticleTable::GetParticleTable()->
  GetIonTable()->GetIonMass(Zt, At);
  targ_mass = 0.939*GeV;  // ~mean neutron and proton ???
 
  G4double proj_mass = aParticle->GetMass();
  G4double proj_momentum = aParticle->GetMomentum().mag();
  G4double sMand = CalcMandelstamS ( proj_mass , targ_mass , proj_momentum );
 
  sMand /= GeV*GeV;  // in GeV for parametrisation
  proj_momentum /= GeV;
  const G4ParticleDefinition* pParticle = aParticle->GetDefinition();
 
  if(pParticle == theNeutron) // as proton ??? 
  {
    xsection = G4double(At)*(21.70*std::pow(sMand,0.0808) + 56.08*std::pow(sMand,-0.4525));
  } 
  else if(pParticle == theProton) 
  {
    xsection = G4double(At)*(21.70*std::pow(sMand,0.0808) + 56.08*std::pow(sMand,-0.4525));
  } 
 
  xsection *= millibarn;
  return xsection;
}

GetHadronNucleonXscNS()

G4double G4GGNuclNuclCrossSection::GetHadronNucleonXscNS	(	G4ParticleDefinition *	pParticle,
		G4double	pTkin,
		G4ParticleDefinition *	tParticle
	)

Definition at line 443 of file G4GGNuclNuclCrossSection.cc.

{
  G4double xsection(0);
  // G4double Delta;   DHW 19 May 2011: variable set but not used
  G4double A0, B0;
  G4double hpXscv(0);
  G4double hnXscv(0);
 
  G4double targ_mass = tParticle->GetPDGMass();
  G4double proj_mass = pParticle->GetPDGMass(); 
 
  G4double proj_energy   = proj_mass + pTkin; 
  G4double proj_momentum = std::sqrt(pTkin*(pTkin+2*proj_mass));
 
  G4double sMand = CalcMandelstamS ( proj_mass , targ_mass , proj_momentum );
 
  sMand         /= GeV*GeV;  // in GeV for parametrisation
  proj_momentum /= GeV;
  proj_energy   /= GeV;
  proj_mass     /= GeV;
 
  // General PDG fit constants
 
  //  G4double s0   = 5.38*5.38; // in Gev^2
  //  G4double eta1 = 0.458;
  //  G4double eta2 = 0.458;
  //  G4double B    = 0.308;
  
  if( proj_momentum >= 373.)
  {
      return GetHadronNucleonXscPDG(pParticle,sMand,tParticle);
  }
  else if( proj_momentum >= 10. ) // high energy: pp = nn = np
    // if( proj_momentum >= 2.)
  {
    //  Delta = 1.;    // DHW 19 May 2011: variable set but not used
    //  if (proj_energy < 40.) Delta = 0.916+0.0021*proj_energy;
 
    if (proj_momentum >= 10.) {
      B0 = 7.5;
      A0 = 100. - B0*std::log(3.0e7);
 
      xsection = A0 + B0*std::log(proj_energy) - 11
                + 103*std::pow(2*0.93827*proj_energy + proj_mass*proj_mass+
                  0.93827*0.93827,-0.165);        //  mb
    }
  }
  else // low energy pp = nn != np
  {
      if(pParticle == tParticle) // pp or nn      // nn to be pp
      {
        if( proj_momentum < 0.73 )
        {
          hnXscv = 23 + 50*( std::pow( std::log(0.73/proj_momentum), 3.5 ) );
        }
        else if( proj_momentum < 1.05  )
        {
          hnXscv = 23 + 40*(std::log(proj_momentum/0.73))*
                         (std::log(proj_momentum/0.73));
        }
        else  // if( proj_momentum < 10.  )
        {
          hnXscv = 39.0 + 
              75*(proj_momentum - 1.2)/(std::pow(proj_momentum,3.0) + 0.15);
        }
        xsection = hnXscv;
      }
      else  // pn to be np
      {
        if( proj_momentum < 0.8 )
        {
          hpXscv = 33+30*std::pow(std::log(proj_momentum/1.3),4.0);
        }      
        else if( proj_momentum < 1.4 )
        {
          hpXscv = 33+30*std::pow(std::log(proj_momentum/0.95),2.0);
        }
        else    // if( proj_momentum < 10.  )
        {
          hpXscv = 33.3+
              20.8*(std::pow(proj_momentum,2.0)-1.35)/
                 (std::pow(proj_momentum,2.50)+0.95);
        }
        xsection = hpXscv;
      }
  }
  xsection *= millibarn; // parametrised in mb
  return xsection;
}

Referenced by GetRatioQE(), and GetRatioSD().

GetHadronNucleonXscPDG()

G4double G4GGNuclNuclCrossSection::GetHadronNucleonXscPDG	(	G4ParticleDefinition *	pParticle,
		G4double	sMand,
		G4ParticleDefinition *	tParticle
	)

Definition at line 386 of file G4GGNuclNuclCrossSection.cc.

{
  G4double xsection = 0.;
  // G4bool pORn = (tParticle == theProton || nucleon == theNeutron  );  
  G4bool proton = (tParticle == theProton);
  G4bool neutron = (tParticle == theNeutron);
 
  // General PDG fit constants
 
  G4double s0   = 5.38*5.38; // in Gev^2
  G4double eta1 = 0.458;
  G4double eta2 = 0.458;
  G4double B    = 0.308;
 
  // const G4ParticleDefinition* pParticle = aParticle->GetDefinition();
  
  if(pParticle == theNeutron) // proton-neutron fit 
  {
    if ( proton )
    {
      xsection = ( 35.80 + B*std::pow(std::log(sMand/s0),2.) 
                  + 40.15*std::pow(sMand,-eta1) - 30.*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = (35.45 + B*std::pow(std::log(sMand/s0),2.) 
         + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2)); // pp for nn
    }
  } 
  else if(pParticle == theProton) 
  {
    if ( proton )
    {
      xsection = (35.45 + B*std::pow(std::log(sMand/s0),2.) 
                + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2));
 
    }
    if ( neutron )
    {
      xsection = (35.80 + B*std::pow(std::log(sMand/s0),2.) 
                  + 40.15*std::pow(sMand,-eta1) - 30.*std::pow(sMand,-eta2));
    }
  } 
  xsection *= millibarn; // parametrised in mb
  return xsection;
}

Referenced by GetHadronNucleonXscNS().

GetHNinelasticXscVU()

G4double G4GGNuclNuclCrossSection::GetHNinelasticXscVU	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt
	)

Definition at line 540 of file G4GGNuclNuclCrossSection.cc.

{
  G4int PDGcode = aParticle->GetDefinition()->GetPDGEncoding();
  G4int absPDGcode = std::abs(PDGcode);
  G4double Elab = aParticle->GetTotalEnergy();              
                          // (s - 2*0.88*GeV*GeV)/(2*0.939*GeV)/GeV;
  G4double Plab = aParticle->GetMomentum().mag();            
                          // std::sqrt(Elab * Elab - 0.88);
 
  Elab /= GeV;
  Plab /= GeV;
 
  G4double LogPlab    = std::log( Plab );
  G4double sqrLogPlab = LogPlab * LogPlab;
 
  //G4cout<<"Plab = "<<Plab<<G4endl;
 
  G4double NumberOfTargetProtons  = Zt; 
  G4double NumberOfTargetNucleons = At;
  G4double NumberOfTargetNeutrons = NumberOfTargetNucleons - NumberOfTargetProtons;
 
  if(NumberOfTargetNeutrons < 0.) NumberOfTargetNeutrons = 0.;
 
  G4double Xtotal = 0., Xelastic = 0., Xinelastic =0.;
 
  if( absPDGcode > 1000 )  //------Projectile is baryon --------
  {
       G4double XtotPP = 48.0 +  0. *std::pow(Plab, 0.  ) +
                         0.522*sqrLogPlab - 4.51*LogPlab;
 
       G4double XtotPN = 47.3 +  0. *std::pow(Plab, 0.  ) +
                         0.513*sqrLogPlab - 4.27*LogPlab;
 
       G4double XelPP  = 11.9 + 26.9*std::pow(Plab,-1.21) +
                         0.169*sqrLogPlab - 1.85*LogPlab;
 
       G4double XelPN  = 11.9 + 26.9*std::pow(Plab,-1.21) +
                         0.169*sqrLogPlab - 1.85*LogPlab;
 
       Xtotal          = ( NumberOfTargetProtons  * XtotPP +
                           NumberOfTargetNeutrons * XtotPN  );
 
       Xelastic        = ( NumberOfTargetProtons  * XelPP  +
                           NumberOfTargetNeutrons * XelPN   );
  }
 
  Xinelastic = Xtotal - Xelastic;
  if(Xinelastic < 0.) Xinelastic = 0.;
 
  return Xinelastic*= millibarn;
}

GetInelasticGlauberGribov()

G4double G4GGNuclNuclCrossSection::GetInelasticGlauberGribov ( const G4DynamicParticle *  dp, G4int  Z, G4int  A  )

inline

Definition at line 146 of file G4GGNuclNuclCrossSection.hh.

{
  GetZandACrossSection(dp, Z, A);
  return fInelasticXsc;
}

GetInelasticGlauberGribovXsc()

G4double G4GGNuclNuclCrossSection::GetInelasticGlauberGribovXsc ( )

inline

Definition at line 103 of file G4GGNuclNuclCrossSection.hh.

{ return fInelasticXsc; }; 

GetNucleusRadius() [1/2]

G4double G4GGNuclNuclCrossSection::GetNucleusRadius	(	const G4DynamicParticle *	,
		const G4Element *	anElement
	)

Definition at line 598 of file G4GGNuclNuclCrossSection.cc.

{
  G4double At = anElement->GetN();
  G4double oneThird = 1.0/3.0;
  G4double cubicrAt = std::pow (At, oneThird); 
 
  G4double R;  // = fRadiusConst*cubicrAt;
  R = fRadiusConst*cubicrAt;
 
  G4double meanA  = 21.;
  G4double tauA1  = 40.; 
  G4double tauA2  = 10.; 
  G4double tauA3  = 5.; 
 
  G4double a1 = 0.85;
  G4double b1 = 1. - a1;
 
  G4double b2 = 0.3;
  G4double b3 = 4.;
 
  if (At > 20.)   // 20.
  {
    R *= ( a1 + b1*std::exp( -(At - meanA)/tauA1) ); 
  }
  else if (At > 3.5)
  {
    R *= ( 1.0 + b2*( 1. - std::exp( (At - meanA)/tauA2) ) ); 
  }
  else 
  {
    R *= ( 1.0 + b3*( 1. - std::exp( (At - meanA)/tauA3) ) ); 
  }
 
  return R;
}

Referenced by GetRatioQE(), GetRatioSD(), and GetZandACrossSection().

GetNucleusRadius() [2/2]

G4double G4GGNuclNuclCrossSection::GetNucleusRadius	(	G4double	Zt,
		G4double	At
	)

Definition at line 640 of file G4GGNuclNuclCrossSection.cc.

{
  G4double R;
  R = GetNucleusRadiusDE(Zt,At);
  // R = GetNucleusRadiusRMS(Zt,At);
 
  return R;
}

GetNucleusRadiusDE()

G4double G4GGNuclNuclCrossSection::GetNucleusRadiusDE	(	G4double	Z,
		G4double	A
	)

Definition at line 680 of file G4GGNuclNuclCrossSection.cc.

{
  // algorithm from diffuse-elastic
 
  G4double R, r0, a11, a12, a13, a2, a3;
 
  a11 = 1.26;  // 1.08, 1.16
  a12 = 1.;  // 1.08, 1.16
  a13 = 1.12;  // 1.08, 1.16
  a2 = 1.1;
  a3 = 1.;
 
  // Special rms radii for light nucleii
 
  if (A < 50.)
  {
    if     (std::abs(A-1.) < 0.5)                         return 0.89*fermi; // p
    else if(std::abs(A-2.) < 0.5)                         return 2.13*fermi; // d
    else if(std::abs(Z-1.) < 0.5 && std::abs(A-3.) < 0.5) return 1.80*fermi; // t
 
    else if(std::abs(Z-2.) < 0.5 && std::abs(A-3.) < 0.5) return 1.96*fermi; // He3
    else if(std::abs(Z-2.) < 0.5 && std::abs(A-4.) < 0.5) return 1.68*fermi; // He4
 
    else if(std::abs(Z-3.) < 0.5)                         return 2.40*fermi; // Li7
    else if(std::abs(Z-4.) < 0.5)                         return 2.51*fermi; // Be9
 
    else if( 10.  < A && A <= 16. ) r0  = a11*( 1 - std::pow(A, -2./3.) )*fermi;   // 1.08*fermi;
    else if( 15.  < A && A <= 20. ) r0  = a12*( 1 - std::pow(A, -2./3.) )*fermi;
    else if( 20.  < A && A <= 30. ) r0  = a13*( 1 - std::pow(A, -2./3.) )*fermi;
    else                            r0  = a2*fermi;
 
    R = r0*std::pow( A, 1./3. );
  }
  else
  {
    r0 = a3*fermi;
 
    R  = r0*std::pow(A, 0.27);
  }
  return R;
}

Referenced by GetNucleusRadius().

GetNucleusRadiusGG()

G4double G4GGNuclNuclCrossSection::GetNucleusRadiusGG

(

G4double

At

)

Definition at line 652 of file G4GGNuclNuclCrossSection.cc.

{
  G4double oneThird = 1.0/3.0;
  G4double cubicrAt = std::pow (At, oneThird); 
 
  G4double R;  // = fRadiusConst*cubicrAt;  
  R = fRadiusConst*cubicrAt;
 
  G4double meanA = 20.;
  G4double tauA  = 20.;
 
  if ( At > 20.)   // 20.
  {
    R *= ( 0.8 + 0.2*std::exp( -(At - meanA)/tauA) ); 
  }
  else
  {
    R *= ( 1.0 + 0.1*( 1. - std::exp( (At - meanA)/tauA) ) ); 
  }
 
  return R;
}

GetNucleusRadiusRMS()

G4double G4GGNuclNuclCrossSection::GetNucleusRadiusRMS	(	G4double	Z,
		G4double	A
	)

Definition at line 728 of file G4GGNuclNuclCrossSection.cc.

{
 
  if     (std::abs(A-1.) < 0.5)                         return 0.89*fermi; // p
  else if(std::abs(A-2.) < 0.5)                         return 2.13*fermi; // d
  else if(std::abs(Z-1.) < 0.5 && std::abs(A-3.) < 0.5) return 1.80*fermi; // t
 
  else if(std::abs(Z-2.) < 0.5 && std::abs(A-3.) < 0.5) return 1.96*fermi; // He3
  else if(std::abs(Z-2.) < 0.5 && std::abs(A-4.) < 0.5) return 1.68*fermi; // He4
 
  else if(std::abs(Z-3.) < 0.5)                         return 2.40*fermi; // Li7
  else if(std::abs(Z-4.) < 0.5)                         return 2.51*fermi; // Be9
 
  else                               return 1.24*std::pow(A, 0.28 )*fermi; // A > 9
}

GetProductionGlauberGribovXsc()

G4double G4GGNuclNuclCrossSection::GetProductionGlauberGribovXsc ( )

inline

Definition at line 104 of file G4GGNuclNuclCrossSection.hh.

{ return fProductionXsc; }; 

GetRadiusConst()

G4double G4GGNuclNuclCrossSection::GetRadiusConst ( )

inline

Definition at line 106 of file G4GGNuclNuclCrossSection.hh.

{ return fRadiusConst;  }; 

GetRatioQE()

G4double G4GGNuclNuclCrossSection::GetRatioQE	(	const G4DynamicParticle *	aParticle,
		G4double	At,
		G4double	Zt
	)

Definition at line 282 of file G4GGNuclNuclCrossSection.cc.

{
  G4double sigma, cofInelastic = 2.4, cofTotal = 2.0, nucleusSquare, ratio;
 
  G4double pZ = aParticle->GetDefinition()->GetPDGCharge();
  G4double pA = aParticle->GetDefinition()->GetBaryonNumber();
 
  G4double pTkin = aParticle->GetKineticEnergy();  
  pTkin /= pA;
 
  G4double pN = pA - pZ;
  if( pN < 0. ) pN = 0.;
 
  G4double tN = tA - tZ;
  if( tN < 0. ) tN = 0.;
 
  G4double tR = GetNucleusRadius(tZ,tA);  
  G4double pR = GetNucleusRadius(pZ,pA); 
 
  sigma = (pZ*tZ+pN*tN)*GetHadronNucleonXscNS(theProton, pTkin, theProton) +
          (pZ*tN+pN*tZ)*GetHadronNucleonXscNS(theProton, pTkin, theNeutron);
 
  nucleusSquare = cofTotal*pi*( pR*pR + tR*tR );   // basically 2piRR
  ratio = sigma/nucleusSquare;
  fInelasticXsc = nucleusSquare*std::log(1. + cofInelastic*ratio)/cofInelastic;
 
  //  sigma = GetHNinelasticXsc(aParticle, tA, tZ);
  ratio = sigma/nucleusSquare;
  fProductionXsc = nucleusSquare*std::log(1. + cofInelastic*ratio)/cofInelastic;
 
  if (fInelasticXsc > fProductionXsc) ratio = (fInelasticXsc-fProductionXsc)/fInelasticXsc;
  else                                ratio = 0.;
  if ( ratio < 0. )                   ratio = 0.;
 
  return ratio; 
}

GetRatioSD()

G4double G4GGNuclNuclCrossSection::GetRatioSD	(	const G4DynamicParticle *	aParticle,
		G4double	At,
		G4double	Zt
	)

Definition at line 242 of file G4GGNuclNuclCrossSection.cc.

{
  G4double sigma, cofInelastic = 2.4, cofTotal = 2.0, nucleusSquare, ratio;
 
  G4double pZ = aParticle->GetDefinition()->GetPDGCharge();
  G4double pA = aParticle->GetDefinition()->GetBaryonNumber();
 
  G4double pTkin = aParticle->GetKineticEnergy();  
  pTkin /= pA;
 
  G4double pN = pA - pZ;
  if( pN < 0. ) pN = 0.;
 
  G4double tN = tA - tZ;
  if( tN < 0. ) tN = 0.;
 
  G4double tR = GetNucleusRadius(tZ,tA);  
  G4double pR = GetNucleusRadius(pZ,pA); 
 
  sigma = (pZ*tZ+pN*tN)*GetHadronNucleonXscNS(theProton, pTkin, theProton) +
          (pZ*tN+pN*tZ)*GetHadronNucleonXscNS(theProton, pTkin, theNeutron);
 
  nucleusSquare = cofTotal*pi*( pR*pR + tR*tR );   // basically 2piRR
  ratio = sigma/nucleusSquare;
  fInelasticXsc = nucleusSquare*std::log(1. + cofInelastic*ratio)/cofInelastic;
  G4double difratio = ratio/(1.+ratio);
 
  fDiffractionXsc = 0.5*nucleusSquare*( difratio - std::log( 1. + difratio ) );
 
  if (fInelasticXsc > 0.) ratio = fDiffractionXsc/fInelasticXsc;
  else                    ratio = 0.;
 
  return ratio; 
}

GetTotalGlauberGribovXsc()

G4double G4GGNuclNuclCrossSection::GetTotalGlauberGribovXsc ( )

inline

Definition at line 101 of file G4GGNuclNuclCrossSection.hh.

{ return fTotalXsc;     }; 

GetZandACrossSection()

G4double G4GGNuclNuclCrossSection::GetZandACrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		G4int	A
	)

Definition at line 107 of file G4GGNuclNuclCrossSection.cc.

{
  G4double xsection;
  G4double sigma;
  G4double cofInelastic = 2.4;
  G4double cofTotal = 2.0;
  G4double nucleusSquare;
  G4double cB;
  G4double ratio;
 
  G4double pZ = aParticle->GetDefinition()->GetPDGCharge();
  G4double pA = aParticle->GetDefinition()->GetBaryonNumber();
 
  G4double pTkin = aParticle->GetKineticEnergy();  
  pTkin /= pA;
 
  G4double pN = pA - pZ;
  if( pN < 0. ) pN = 0.;
 
  G4double tN = tA - tZ;
  if( tN < 0. ) tN = 0.;
 
  G4double tR = GetNucleusRadius( G4double(tZ),G4double(tA) );  
  G4double pR = GetNucleusRadius(pZ,pA); 
 
  cB = GetCoulombBarier(aParticle, G4double(tZ), G4double(tA), pR, tR);
 
  if ( cB > 0. ) 
  {
    G4DynamicParticle* dProton = new G4DynamicParticle(theProton,
                                              G4ParticleMomentum(1.,0.,0.), 
                                              pTkin);
 
    G4DynamicParticle* dNeutron = new G4DynamicParticle(theNeutron,
                                              G4ParticleMomentum(1.,0.,0.), 
                                              pTkin);
 
    sigma = (pZ*tZ+pN*tN)*hnXsc->GetHadronNucleonXscNS(dProton, theProton);
 
    G4double ppInXsc = hnXsc->GetInelasticHadronNucleonXsc();
 
    sigma += (pZ*tN+pN*tZ)*hnXsc->GetHadronNucleonXscNS(dNeutron, theProton);
 
    G4double npInXsc = hnXsc->GetInelasticHadronNucleonXsc();
 
    // G4cout<<"ppInXsc = "<<ppInXsc/millibarn<<"; npInXsc = "<<npInXsc/millibarn<<G4endl;
    // G4cout<<"npTotXsc = "<<hnXsc->GetTotalHadronNucleonXsc()/millibarn<<"; npElXsc = "
    //                      <<hnXsc->GetElasticHadronNucleonXsc()/millibarn<<G4endl;
 
    nucleusSquare = cofTotal*pi*( pR*pR + tR*tR );   // basically 2piRR
 
    ratio      = sigma/nucleusSquare;
    xsection   = nucleusSquare*std::log( 1. + ratio );
    fTotalXsc  = xsection;
    fTotalXsc *= cB;
 
    fInelasticXsc = nucleusSquare*std::log( 1. + cofInelastic*ratio )/cofInelastic;
 
    fInelasticXsc *= cB;
    fElasticXsc   = fTotalXsc - fInelasticXsc;
 
    // if (fElasticXsc < DBL_MIN) fElasticXsc = DBL_MIN;
    /*    
    G4double difratio = ratio/(1.+ratio);
 
    fDiffractionXsc = 0.5*nucleusSquare*( difratio - std::log( 1. + difratio ) );
    */
    // production to be checked !!! edit MK xsc
 
    //sigma = (pZ*tZ+pN*tN)*GetHadronNucleonXscMK(theProton, pTkin, theProton) +
    //      (pZ*tN+pN*tZ)*GetHadronNucleonXscMK(theProton, pTkin, theNeutron);
 
    sigma = (pZ*tZ+pN*tN)*ppInXsc + (pZ*tN+pN*tZ)*npInXsc;
 
    ratio = sigma/nucleusSquare;
    fProductionXsc = nucleusSquare*std::log( 1. + cofInelastic*ratio )/cofInelastic;
 
    if (fElasticXsc < 0.) fElasticXsc = 0.;
  }
  else
  {
    fInelasticXsc  = 0.;
    fTotalXsc      = 0.;
    fElasticXsc    = 0.;
    fProductionXsc = 0.;
  }
  return fInelasticXsc;   // xsection; 
}

Referenced by GetElasticGlauberGribov(), GetElementCrossSection(), and GetInelasticGlauberGribov().

IsElementApplicable()

G4bool G4GGNuclNuclCrossSection::IsElementApplicable ( const G4DynamicParticle *  , G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 72 of file G4GGNuclNuclCrossSection.cc.

{
  G4bool applicable = true;
//  G4double kineticEnergy = aDP->GetKineticEnergy();
 
//  if (kineticEnergy >= fLowerLimit) applicable = true;
  return applicable;
}

SetEnergyLowerLimit()

void G4GGNuclNuclCrossSection::SetEnergyLowerLimit ( G4double  E )

inline

Definition at line 115 of file G4GGNuclNuclCrossSection.hh.

{fLowerLimit=E;};

---

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

• G4GGNuclNuclCrossSection.hh
• G4GGNuclNuclCrossSection.cc