#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	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	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 G4int	GetVerboseLevel () const

virtual 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 ( )

Definition at line 46 of file G4GGNuclNuclCrossSection.cc.

References G4Neutron::Neutron(), and G4Proton::Proton().

  : 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 ( )

virtual

Definition at line 61 of file G4GGNuclNuclCrossSection.cc.

 {
   delete hnXsc;
 }

Member Function Documentation

virtual void G4GGNuclNuclCrossSection::BuildPhysicsTable ( const G4ParticleDefinition & )

inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 75 of file G4GGNuclNuclCrossSection.hh.

76 {}

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

Definition at line 774 of file G4GGNuclNuclCrossSection.cc.

Referenced by GetHadronNucleonXsc(), and GetHadronNucleonXscNS().

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

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

Definition at line 757 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;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 67 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";
 }

static const char* G4GGNuclNuclCrossSection::Default_Name ( )

inlinestatic

Definition at line 57 of file G4GGNuclNuclCrossSection.hh.

Referenced by DMXPhysicsList::ConstructHad(), and GammaRayTelIonPhysics::ConstructProcess().

57 {return "Glauber-Gribov nucleus nucleus";}

virtual void G4GGNuclNuclCrossSection::DumpPhysicsTable ( const G4ParticleDefinition & )

inlinevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 79 of file G4GGNuclNuclCrossSection.hh.

References G4cout, and G4endl.

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

G4cout

G4GLOB_DLL std::ostream G4cout

G4endl

#define G4endl

Definition: G4ios.hh:61

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

Definition at line 211 of file G4GGNuclNuclCrossSection.cc.

References python.hepunit::fine_structure_const, G4DynamicParticle::GetDefinition(), G4IonTable::GetIonMass(), G4ParticleTable::GetIonTable(), G4DynamicParticle::GetKineticEnergy(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), and python.hepunit::hbarc.

Referenced by GetZandACrossSection().

 {
   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;
 }

G4double G4GGNuclNuclCrossSection::GetDiffractionGlauberGribovXsc ( )

inline

Definition at line 105 of file G4GGNuclNuclCrossSection.hh.

105 { return fDiffractionXsc; };

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

inline

Definition at line 136 of file G4GGNuclNuclCrossSection.hh.

References GetZandACrossSection().

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

G4double G4GGNuclNuclCrossSection::GetElasticGlauberGribovXsc ( )

inline

Definition at line 102 of file G4GGNuclNuclCrossSection.hh.

102 { return fElasticXsc; };

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 95 of file G4GGNuclNuclCrossSection.cc.

References G4lrint(), GetZandACrossSection(), and G4NistManager::Instance().

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

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

Definition at line 336 of file G4GGNuclNuclCrossSection.cc.

References G4lrint(), G4Element::GetN(), and G4Element::GetZ().

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

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

Definition at line 355 of file G4GGNuclNuclCrossSection.cc.

References CalcMandelstamS(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentum(), G4ParticleTable::GetParticleTable(), python.hepunit::GeV, CLHEP::Hep3Vector::mag(), and python.hepunit::millibarn.

 {
   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;
 }

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

Definition at line 451 of file G4GGNuclNuclCrossSection.cc.

References CalcMandelstamS(), GetHadronNucleonXscPDG(), G4ParticleDefinition::GetPDGMass(), python.hepunit::GeV, and python.hepunit::millibarn.

Referenced by GetRatioQE(), and GetRatioSD().

 {
   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;
 }

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

Definition at line 394 of file G4GGNuclNuclCrossSection.cc.

References python.hepunit::millibarn, neutron, G4InuclParticleNames::proton, and G4InuclParticleNames::s0.

Referenced by GetHadronNucleonXscNS().

 {
   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;
 }

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

Definition at line 548 of file G4GGNuclNuclCrossSection.cc.

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetMomentum(), G4ParticleDefinition::GetPDGEncoding(), G4DynamicParticle::GetTotalEnergy(), python.hepunit::GeV, CLHEP::Hep3Vector::mag(), and python.hepunit::millibarn.

 {
   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;
 }

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

inline

Definition at line 146 of file G4GGNuclNuclCrossSection.hh.

References GetZandACrossSection().

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

G4double G4GGNuclNuclCrossSection::GetInelasticGlauberGribovXsc ( )

inline

Definition at line 103 of file G4GGNuclNuclCrossSection.hh.

103 { return fInelasticXsc; };

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

Definition at line 606 of file G4GGNuclNuclCrossSection.cc.

References G4Element::GetN(), and G4INCL::Math::oneThird.

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

 {
   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;
 }

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

Definition at line 648 of file G4GGNuclNuclCrossSection.cc.

References GetNucleusRadiusDE().

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

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

Definition at line 688 of file G4GGNuclNuclCrossSection.cc.

References python.hepunit::fermi.

Referenced by GetNucleusRadius().

 {
   // 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;
 }

G4double G4GGNuclNuclCrossSection::GetNucleusRadiusGG ( G4double At )

Definition at line 660 of file G4GGNuclNuclCrossSection.cc.

References G4INCL::Math::oneThird.

 {
   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;
 }

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

Definition at line 736 of file G4GGNuclNuclCrossSection.cc.

References python.hepunit::fermi.

 {
 
   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
 }

G4double G4GGNuclNuclCrossSection::GetProductionGlauberGribovXsc ( )

inline

Definition at line 104 of file G4GGNuclNuclCrossSection.hh.

104 { return fProductionXsc; };

G4double G4GGNuclNuclCrossSection::GetRadiusConst ( )

inline

Definition at line 106 of file G4GGNuclNuclCrossSection.hh.

106 { return fRadiusConst; };

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

Definition at line 291 of file G4GGNuclNuclCrossSection.cc.

References G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), GetHadronNucleonXscNS(), G4DynamicParticle::GetKineticEnergy(), GetNucleusRadius(), G4ParticleDefinition::GetPDGCharge(), and python.hepunit::pi.

 {
   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; 
 }

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

Definition at line 251 of file G4GGNuclNuclCrossSection.cc.

References G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), GetHadronNucleonXscNS(), G4DynamicParticle::GetKineticEnergy(), GetNucleusRadius(), G4ParticleDefinition::GetPDGCharge(), and python.hepunit::pi.

 {
   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; 
 }

G4double G4GGNuclNuclCrossSection::GetTotalGlauberGribovXsc ( )

inline

Definition at line 101 of file G4GGNuclNuclCrossSection.hh.

101 { return fTotalXsc; };

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

Definition at line 112 of file G4GGNuclNuclCrossSection.cc.

References G4ParticleDefinition::GetBaryonNumber(), GetCoulombBarier(), G4DynamicParticle::GetDefinition(), G4HadronNucleonXsc::GetHadronNucleonXscNS(), G4HadronNucleonXsc::GetInelasticHadronNucleonXsc(), G4DynamicParticle::GetKineticEnergy(), GetNucleusRadius(), G4ParticleDefinition::GetPDGCharge(), and python.hepunit::pi.

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

 {
   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();
 
       delete dProton;
       delete dNeutron;
       
     // 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; 
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 76 of file G4GGNuclNuclCrossSection.cc.

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

void G4GGNuclNuclCrossSection::SetEnergyLowerLimit ( G4double E )

inline

Definition at line 115 of file G4GGNuclNuclCrossSection.hh.

115 {fLowerLimit=E;};

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

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation