G4HadronNucleonXsc Class Reference

#include <G4HadronNucleonXsc.hh>

Public Member Functions
	G4HadronNucleonXsc ()
virtual	~G4HadronNucleonXsc ()
virtual G4bool	IsApplicable (const G4DynamicParticle *aDP, const G4Element *)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *aDP, G4int Z, G4int A)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
void	CrossSectionDescription (std::ostream &) const
G4double	GetHadronNucleonXscEL (const G4DynamicParticle *, const G4ParticleDefinition *)
G4double	GetHadronNucleonXscPDG (const G4DynamicParticle *, const G4ParticleDefinition *)
G4double	GetHadronNucleonXscNS (const G4DynamicParticle *, const G4ParticleDefinition *)
G4double	GetHadronNucleonXscVU (const G4DynamicParticle *, const G4ParticleDefinition *)
G4double	CalculateEcmValue (const G4double, const G4double, const G4double)
G4double	CalcMandelstamS (const G4double, const G4double, const G4double)
G4double	GetCoulombBarrier (const G4DynamicParticle *aParticle, const G4ParticleDefinition *nucleon)
G4double	GetTotalHadronNucleonXsc ()
G4double	GetElasticHadronNucleonXsc ()
G4double	GetInelasticHadronNucleonXsc ()
void	InitialiseKaonNucleonTotXsc ()
G4double	GetKpProtonTotXscVector (G4double logEnergy)
G4double	GetKpNeutronTotXscVector (G4double logEnergy)
G4double	GetKmProtonTotXscVector (G4double logEnergy)
G4double	GetKmNeutronTotXscVector (G4double logEnergy)

---

Detailed Description

Definition at line 51 of file G4HadronNucleonXsc.hh.

---

Constructor & Destructor Documentation

G4HadronNucleonXsc::G4HadronNucleonXsc

(

)

Definition at line 39 of file G4HadronNucleonXsc.cc.

References G4Alpha::Alpha(), G4AntiLambda::AntiLambda(), G4AntiNeutron::AntiNeutron(), G4AntiOmegaMinus::AntiOmegaMinus(), G4AntiProton::AntiProton(), G4AntiSigmaMinus::AntiSigmaMinus(), G4AntiSigmaPlus::AntiSigmaPlus(), G4AntiSigmaZero::AntiSigmaZero(), G4AntiXiMinus::AntiXiMinus(), G4AntiXiZero::AntiXiZero(), G4Deuteron::Deuteron(), G4Gamma::Gamma(), G4He3::He3(), InitialiseKaonNucleonTotXsc(), G4KaonMinus::KaonMinus(), G4KaonPlus::KaonPlus(), G4KaonZeroLong::KaonZeroLong(), G4KaonZeroShort::KaonZeroShort(), G4Lambda::Lambda(), G4Neutron::Neutron(), G4OmegaMinus::OmegaMinus(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), G4PionZero::PionZero(), G4Proton::Proton(), G4SigmaMinus::SigmaMinus(), G4SigmaPlus::SigmaPlus(), G4SigmaZero::SigmaZero(), G4Triton::Triton(), G4XiMinus::XiMinus(), and G4XiZero::XiZero().

: fUpperLimit( 10000 * GeV ),
  fLowerLimit( 0.03 * MeV ),
  fTotalXsc(0.0), fElasticXsc(0.0), fInelasticXsc(0.0), fHadronNucleonXsc(0.0)
{
  theGamma    = G4Gamma::Gamma();
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  theAProton  = G4AntiProton::AntiProton();
  theANeutron = G4AntiNeutron::AntiNeutron();
  thePiPlus   = G4PionPlus::PionPlus();
  thePiMinus  = G4PionMinus::PionMinus();
  thePiZero   = G4PionZero::PionZero();
  theKPlus    = G4KaonPlus::KaonPlus();
  theKMinus   = G4KaonMinus::KaonMinus();
  theK0S      = G4KaonZeroShort::KaonZeroShort();
  theK0L      = G4KaonZeroLong::KaonZeroLong();
  theL        = G4Lambda::Lambda();
  theAntiL    = G4AntiLambda::AntiLambda();
  theSPlus    = G4SigmaPlus::SigmaPlus();
  theASPlus   = G4AntiSigmaPlus::AntiSigmaPlus();
  theSMinus   = G4SigmaMinus::SigmaMinus();
  theASMinus  = G4AntiSigmaMinus::AntiSigmaMinus();
  theS0       = G4SigmaZero::SigmaZero();
  theAS0      = G4AntiSigmaZero::AntiSigmaZero();
  theXiMinus  = G4XiMinus::XiMinus();
  theXi0      = G4XiZero::XiZero();
  theAXiMinus = G4AntiXiMinus::AntiXiMinus();
  theAXi0     = G4AntiXiZero::AntiXiZero();
  theOmega    = G4OmegaMinus::OmegaMinus();
  theAOmega   = G4AntiOmegaMinus::AntiOmegaMinus();
  theD        = G4Deuteron::Deuteron();
  theT        = G4Triton::Triton();
  theA        = G4Alpha::Alpha();
  theHe3      = G4He3::He3();

  InitialiseKaonNucleonTotXsc();
}

G4HadronNucleonXsc::~G4HadronNucleonXsc

(

)

[virtual]

Definition at line 79 of file G4HadronNucleonXsc.cc.

{}

---

Member Function Documentation

G4double G4HadronNucleonXsc::CalcMandelstamS	(	const	G4double,
		const	G4double,
		const	G4double
	)

Definition at line 1125 of file G4HadronNucleonXsc.cc.

Referenced by GetHadronNucleonXscEL(), GetHadronNucleonXscNS(), and GetHadronNucleonXscPDG().

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

  return sMand;
}

G4double G4HadronNucleonXsc::CalculateEcmValue	(	const	G4double,
		const	G4double,
		const	G4double
	)

Definition at line 1108 of file G4HadronNucleonXsc.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 G4HadronNucleonXsc::CrossSectionDescription

(

std::ostream &

)

const

Definition at line 82 of file G4HadronNucleonXsc.cc.

{
  outFile << "G4HadronNucleonXsc calculates the total, inelastic and elastic\n"
          << "hadron-nucleon cross sections using several different\n"
          << "parameterizations within the Glauber-Gribov approach. It is\n"
          << "valid for all incident gammas and long-lived hadrons at\n"
          << "energies above 30 keV.  This is a cross section component which\n"
          << "is to be used to build a cross section data set.\n"; 
}

virtual void G4HadronNucleonXsc::DumpPhysicsTable

(

const G4ParticleDefinition &

)

[inline, virtual]

Definition at line 65 of file G4HadronNucleonXsc.hh.

References G4cout, and G4endl.

  {G4cout << "G4HadronNucleonXsc: uses parametrisation"<<G4endl;}

G4double G4HadronNucleonXsc::GetCoulombBarrier	(	const G4DynamicParticle *	aParticle,
		const G4ParticleDefinition *	nucleon
	)

Definition at line 1140 of file G4HadronNucleonXsc.cc.

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4ParticleDefinition::GetPDGCharge(), and G4ParticleDefinition::GetPDGMass().

Referenced by GetHadronNucleonXscNS().

{
  G4double ratio;

  G4double tR = 0.895*fermi, pR;

  if     ( aParticle->GetDefinition() == theProton ) pR = 0.895*fermi;
  else if( aParticle->GetDefinition() == thePiPlus ) pR = 0.663*fermi;
  else if( aParticle->GetDefinition() == theKPlus )  pR = 0.340*fermi;
  else                                               pR = 0.500*fermi;

  G4double pZ = aParticle->GetDefinition()->GetPDGCharge();
  G4double tZ = nucleon->GetPDGCharge();

  G4double pTkin = aParticle->GetKineticEnergy();
  
  G4double pM    = aParticle->GetDefinition()->GetPDGMass(); 
  G4double tM    = nucleon->GetPDGMass();

  G4double pElab = pTkin + pM;

  G4double totEcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM);

  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 G4HadronNucleonXsc::GetElasticHadronNucleonXsc

(

)

[inline]

Definition at line 89 of file G4HadronNucleonXsc.hh.

Referenced by G4BGGPionElasticXS::BuildPhysicsTable(), G4BGGNucleonElasticXS::BuildPhysicsTable(), G4NeutronElasticXS::GetElementCrossSection(), G4BGGPionElasticXS::GetIsoCrossSection(), and G4BGGNucleonElasticXS::GetIsoCrossSection().

{ return fElasticXsc;   }; 

G4double G4HadronNucleonXsc::GetHadronNucleonXscEL	(	const G4DynamicParticle *	,
		const G4ParticleDefinition *
	)

Definition at line 141 of file G4HadronNucleonXsc.cc.

References CalcMandelstamS(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetMass(), and G4DynamicParticle::GetMomentum().

{
  G4double xsection;


  G4double 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* theParticle = aParticle->GetDefinition();

  G4bool pORn = (nucleon == theProton || nucleon == theNeutron  );  
  

  if(theParticle == theGamma && pORn ) 
  {
    xsection = (0.0677*std::pow(sMand,0.0808) + 0.129*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == theNeutron && pORn ) // as proton ??? 
  {
    xsection = (21.70*std::pow(sMand,0.0808) + 56.08*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == theProton && pORn ) 
  {
    xsection = (21.70*std::pow(sMand,0.0808) + 56.08*std::pow(sMand,-0.4525));

    // xsection = At*( 49.51*std::pow(sMand,-0.097) + 0.314*std::log(sMand)*std::log(sMand) );
    // xsection = At*( 38.4 + 0.85*std::abs(std::pow(log(sMand),1.47)) );
  } 
  else if(theParticle == theAProton && pORn ) 
  {
    xsection = ( 21.70*std::pow(sMand,0.0808) + 98.39*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == thePiPlus && pORn ) 
  {
    xsection = (13.63*std::pow(sMand,0.0808) + 27.56*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == thePiMinus && pORn ) 
  {
    // xsection = At*( 55.2*std::pow(sMand,-0.255) + 0.346*std::log(sMand)*std::log(sMand) );
    xsection = (13.63*std::pow(sMand,0.0808) + 36.02*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == theKPlus && pORn ) 
  {
    xsection = (11.82*std::pow(sMand,0.0808) + 8.15*std::pow(sMand,-0.4525));
  } 
  else if(theParticle == theKMinus && pORn ) 
  {
    xsection = (11.82*std::pow(sMand,0.0808) + 26.36*std::pow(sMand,-0.4525));
  }
  else  // as proton ??? 
  {
    xsection = (21.70*std::pow(sMand,0.0808) + 56.08*std::pow(sMand,-0.4525));
  }
  xsection *= millibarn;

  fTotalXsc     = xsection;
  fInelasticXsc = 0.83*xsection;
  fElasticXsc   = fTotalXsc - fInelasticXsc;
  if (fElasticXsc < 0.)fElasticXsc = 0.;
 
  return xsection;
}

G4double G4HadronNucleonXsc::GetHadronNucleonXscNS	(	const G4DynamicParticle *	,
		const G4ParticleDefinition *
	)

Definition at line 363 of file G4HadronNucleonXsc.cc.

References CalcMandelstamS(), GetCoulombBarrier(), G4DynamicParticle::GetDefinition(), GetHadronNucleonXscPDG(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentum(), G4ParticleDefinition::GetPDGCharge(), G4DynamicParticle::GetTotalEnergy(), neutron, G4InuclParticleNames::proton, and G4InuclParticleNames::s0.

Referenced by G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4GlauberGribovCrossSection::GetIsoCrossSection(), G4ComponentGGHadronNucleusXsc::GetIsoCrossSection(), G4BGGNucleonInelasticXS::GetIsoCrossSection(), G4GGNuclNuclCrossSection::GetZandACrossSection(), and G4ComponentGGNuclNuclXsc::GetZandACrossSection().

{
  G4double xsection(0); 
  
  G4double A0, B0;
  G4double hpXsc(0);
  G4double hnXsc(0);


  G4double tM = 0.939*GeV;  // ~mean neutron and proton ???

  G4double pM   = aParticle->GetMass();
  G4double pE   = aParticle->GetTotalEnergy(); 
  G4double pLab = aParticle->GetMomentum().mag();

  G4double sMand = CalcMandelstamS ( pM , tM , pLab );

  sMand         /= GeV*GeV;  // in GeV for parametrisation
  pLab /= GeV;
  pE   /= GeV;
  pM     /= GeV;

  G4double logP = std::log(pLab);


  // 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* theParticle = aParticle->GetDefinition();

  G4bool pORn = (nucleon == theProton || nucleon == theNeutron  );  
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);

  if( theParticle == theNeutron && pORn ) 
  {
    if( pLab >= 373.)
    {
      xsection =  GetHadronNucleonXscPDG(aParticle, nucleon)/millibarn;

      fElasticXsc = 6.5 + 0.308*std::pow(std::log(sMand/400.),1.65) + 9.19*std::pow(sMand,-0.458);
    
      fTotalXsc = xsection;
    
    }
    else if( pLab >= 100.)
    {
      B0 = 7.5;
      A0 = 100. - B0*std::log(3.0e7);

      xsection = A0 + B0*std::log(pE) - 11
          // + 103*std::pow(2*0.93827*pE + pM*pM+0.93827*0.93827,-0.165);        //  mb
                  + 103*std::pow(sMand,-0.165);        //  mb

      fElasticXsc = 5.53 + 0.308*std::pow(std::log(sMand/28.9),1.1) + 9.19*std::pow(sMand,-0.458);
      
      fTotalXsc = xsection;
    }
    else if( pLab >= 10.)
    {
        B0 = 7.5;
        A0 = 100. - B0*std::log(3.0e7);

        xsection = A0 + B0*std::log(pE) - 11
                  + 103*std::pow(2*0.93827*pE + pM*pM+
                     0.93827*0.93827,-0.165);        //  mb      
      fTotalXsc = xsection;
      fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
    }
    else  // pLab < 10 GeV/c
    {
      if( neutron )      // nn to be pp
      {
        if( pLab < 0.4 )
        {
          hnXsc = 23 + 50*( std::pow( std::log(0.73/pLab), 3.5 ) );
          fElasticXsc = hnXsc;
        }
        else if( pLab < 0.73 )
        {
          hnXsc = 23 + 50*( std::pow( std::log(0.73/pLab), 3.5 ) );
          fElasticXsc = hnXsc; 
        }
        else if( pLab < 1.05  )
        {
          hnXsc = 23 + 40*(std::log(pLab/0.73))*
                         (std::log(pLab/0.73));
          fElasticXsc = 23 + 20*(std::log(pLab/0.73))*
                         (std::log(pLab/0.73));
        }
        else    // 1.05 - 10 GeV/c
        {
          hnXsc = 39.0+75*(pLab - 1.2)/(std::pow(pLab,3.0) + 0.15);

          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
        fTotalXsc = hnXsc;
      }
      if( proton )   // pn to be np
      {
        if( pLab < 0.02 )
        {
          hpXsc = 4100+30*std::pow(std::log(1.3/pLab),3.6); // was as pLab < 0.8
          fElasticXsc = hpXsc;
        }      
        else if( pLab < 0.8 )
        {
          hpXsc = 33+30*std::pow(std::log(pLab/1.3),4.0);
          fElasticXsc = hpXsc;
        }      
        else if( pLab < 1.05 )
        {
          hpXsc = 33+30*std::pow(std::log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( std::log(0.511/pLab)*std::log(0.511/pLab) + 1.6 );
        }
        else if( pLab < 1.4 )
        {
          hpXsc = 33+30*std::pow(std::log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( std::log(0.511/pLab)*std::log(0.511/pLab) + 1.6 );
        }
        else    // 1.4 < pLab < 10.  )
        {
          hpXsc = 33.3 + 20.8*(std::pow(pLab,2.0) - 1.35)/(std::pow(pLab,2.50) + 0.95);
          
          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
        fTotalXsc = hpXsc;
      }
    }
  } 
  else if( theParticle == theProton && pORn ) 
  {
    if( pLab >= 373.) // pdg due to TOTEM data
    {
      xsection =  GetHadronNucleonXscPDG(aParticle, nucleon)/millibarn;

      fElasticXsc = 6.5 + 0.308*std::pow(std::log(sMand/400.),1.65) + 9.19*std::pow(sMand,-0.458);
     
      fTotalXsc = xsection;
    }
    else if( pLab >= 100.)
    {
      B0 = 7.5;
      A0 = 100. - B0*std::log(3.0e7);

      xsection = A0 + B0*std::log(pE) - 11 + 103*std::pow(sMand,-0.165);        //  mb

      fElasticXsc = 5.53 + 0.308*std::pow(std::log(sMand/28.9),1.1) + 9.19*std::pow(sMand,-0.458);
      
      fTotalXsc = xsection;
    }
    else if( pLab >= 10.)
    {
      B0 = 7.5;
      A0 = 100. - B0*std::log(3.0e7);

      xsection = A0 + B0*std::log(pE) - 11 + 103*std::pow(sMand,-0.165);        //  mb

      fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
      
      fTotalXsc = xsection;
    }
    else
    {
      // pp

      if( proton )
      {
        if( pLab < 0.4 )
        {
          hpXsc = 23 + 50*( std::pow( std::log(0.73/pLab), 3.5 ) );
          fElasticXsc = hpXsc;
        }
        else if( pLab < 0.73 )
        {
          hpXsc = 23 + 50*( std::pow( std::log(0.73/pLab), 3.5 ) );
          fElasticXsc = hpXsc; 
        }
        else if( pLab < 1.05  )
        {
          hpXsc = 23 + 40*(std::log(pLab/0.73))*
                         (std::log(pLab/0.73));
          fElasticXsc = 23 + 20*(std::log(pLab/0.73))*
                         (std::log(pLab/0.73));
        }
        else    // 1.05 - 10 GeV/c
        {
          hpXsc = 39.0+75*(pLab - 1.2)/(std::pow(pLab,3.0) + 0.15);

          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
        fTotalXsc = hpXsc;
      }
      if( neutron )     // pn to be np
      {
        if( pLab < 0.02 )
        {
          hnXsc = 4100+30*std::pow(std::log(1.3/pLab),3.6); // was as pLab < 0.8
          fElasticXsc = hnXsc;
        }      
        else if( pLab < 0.8 )
        {
          hnXsc = 33+30*std::pow(std::log(pLab/1.3),4.0);
          fElasticXsc = hnXsc;
        }      
        else if( pLab < 1.05 )
        {
          hnXsc = 33+30*std::pow(std::log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( std::log(0.511/pLab)*std::log(0.511/pLab) + 1.6 );
        }
        else if( pLab < 1.4 )
        {
          hnXsc = 33+30*std::pow(std::log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( std::log(0.511/pLab)*std::log(0.511/pLab) + 1.6 );
        }
        else    // 1.4 < pLab < 10.  )
        {
          hnXsc = 33.3 + 20.8*(std::pow(pLab,2.0) - 1.35)/(std::pow(pLab,2.50) + 0.95);
          
          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
        fTotalXsc = hnXsc;
      }
    }    
  } 
  else if( theParticle == theAProton && pORn ) 
  {
    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);
    }
    fTotalXsc = xsection;
  } 
  else if( theParticle == thePiPlus && pORn ) // pi+ /////////////////////////////////////////////
  {
    if( proton ) // pi+ p
    {
      if( pLab < 0.28 )
      {
        hpXsc       = 10./((logP + 1.273)*(logP + 1.273) + 0.05);
        fElasticXsc = hpXsc;
      }
      else if( pLab < 0.4 )
      {
        hpXsc       = 14./( (logP + 1.273)*(logP + 1.273) + 0.07);
        fElasticXsc = hpXsc;
      }
      else if( pLab < 0.68 )
      {
        hpXsc       = 14./( (logP + 1.273)*(logP + 1.273) + 0.07);
        fElasticXsc = hpXsc;
      }
      else if( pLab < 0.85 )
      {
        G4double Ex4 = 88*(std::log(pLab/0.77))*(std::log(pLab/0.77));
        hpXsc        = Ex4 + 14.9;
        fElasticXsc = hpXsc*std::exp(-3.*(pLab - 0.68));  
      }
      else if( pLab < 1.15 )
      {
        G4double Ex4 = 88*(std::log(pLab/0.77))*(std::log(pLab/0.77));
        hpXsc        = Ex4 + 14.9;

        fElasticXsc = 6.0 + 1.4/(( pLab - 1.4)*( pLab - 1.4) + 0.1);
      }
      else if( pLab < 1.4) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hpXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 6.0 + 1.4/(( pLab - 1.4)*( pLab - 1.4) + 0.1);
      }
      else if( pLab < 2.0 ) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hpXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 3.0 + 1.36/( (logP - 0.336)*(logP - 0.336) + 0.08);    
      }
      else if( pLab < 3.5 ) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hpXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      else if( pLab < 200. ) // my
      {
        hpXsc = 10.6 + 2.*std::log(pE) + 25*std::pow(pE, -0.43 ); // ns original
        // hpXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      else //  pLab > 100 // my
      {
        hpXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      fTotalXsc = hpXsc;
    }    
    if( neutron )  // pi+ n = pi- p??
    {
      if( pLab < 0.28 ) 
      {
        hnXsc       = 0.288/((pLab - 0.28)*(pLab - 0.28) + 0.004);
        fElasticXsc = 1.8/((logP + 1.273)*(logP + 1.273) + 0.07);
      }
      else if( pLab < 0.395676 ) // first peak
      {
        hnXsc       = 0.648/((pLab - 0.28)*(pLab - 0.28) + 0.009);
        fElasticXsc = 0.257/((pLab - 0.28)*(pLab - 0.28) + 0.01);
       }
      else if( pLab < 0.5 )
      {
        hnXsc       = 26 + 110*(std::log(pLab/0.48))*(std::log(pLab/0.48));
        fElasticXsc = 0.37*hnXsc;
      }
      else if( pLab < 0.65 )
      {
        hnXsc       = 26 + 110*(std::log(pLab/0.48))*(std::log(pLab/0.48));
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.72 )
      {
        hnXsc = 36.1 + 10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.88 )
      {
        hnXsc = 36.1 + 10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 1.03 )
      {
        hnXsc = 36.1 + 10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 2.0 + 0.4/((pLab - 1.03)*(pLab - 1.03) + 0.016);
      }
      else if( pLab < 1.15 )
      {
        hnXsc = 36.1 + 10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 2.0 + 0.4/((pLab - 1.03)*(pLab - 1.03) + 0.016);
      }
      else if( pLab < 1.3 )
      {
        hnXsc = 36.1 + 10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 3. + 13./pLab;
      }
      else if( pLab < 2.6 ) // < 3.0) // ns original
      {
        hnXsc = 36.1 + 0.079-4.313*std::log(pLab)+
                3*std::exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*std::exp(-(pLab-1.4)*(pLab-1.4)/0.12/0.12);
        fElasticXsc = 3. + 13./pLab; 
      }
      else if( pLab < 20. ) // < 3.0) // ns original
      {
        hnXsc = 36.1 + 0.079 - 4.313*std::log(pLab)+
                3*std::exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*std::exp(-(pLab-1.4)*(pLab-1.4)/0.12/0.12);
        fElasticXsc = 3. + 13./pLab; 
      }
      else   // mb 
      {
        hnXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3. + 13./pLab;
      }
      fTotalXsc = hnXsc;
    }
  } 
  else if( theParticle == thePiMinus && pORn ) 
  {
    if( neutron )     // pi- n = pi+ p??
    {
      if( pLab < 0.28 )
      {
        hnXsc       = 10./((logP + 1.273)*(logP + 1.273) + 0.05);
        fElasticXsc = hnXsc;
      }
      else if( pLab < 0.4 )
      {
        hnXsc       = 14./( (logP + 1.273)*(logP + 1.273) + 0.07);
        fElasticXsc = hnXsc;
      }
      else if( pLab < 0.68 )
      {
        hnXsc       = 14./( (logP + 1.273)*(logP + 1.273) + 0.07);
        fElasticXsc = hnXsc;
      }
      else if( pLab < 0.85 )
      {
        G4double Ex4 = 88*(std::log(pLab/0.77))*(std::log(pLab/0.77));
        hnXsc        = Ex4 + 14.9;
        fElasticXsc = hnXsc*std::exp(-3.*(pLab - 0.68));  
      }
      else if( pLab < 1.15 )
      {
        G4double Ex4 = 88*(std::log(pLab/0.77))*(std::log(pLab/0.77));
        hnXsc        = Ex4 + 14.9;

        fElasticXsc = 6.0 + 1.4/(( pLab - 1.4)*( pLab - 1.4) + 0.1);
      }
      else if( pLab < 1.4) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hnXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 6.0 + 1.4/(( pLab - 1.4)*( pLab - 1.4) + 0.1);
      }
      else if( pLab < 2.0 ) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hnXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 3.0 + 1.36/( (logP - 0.336)*(logP - 0.336) + 0.08);    
      }
      else if( pLab < 3.5 ) // ns original
      {
        G4double Ex1 = 3.2*std::exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*std::exp(-(pLab-1.47)*(pLab-1.47)/0.225/0.225);
        hnXsc        = Ex1 + Ex2 + 27.5;
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      else if( pLab < 200. ) // my
      {
        hnXsc = 10.6 + 2.*std::log(pE) + 25*std::pow(pE, -0.43 ); // ns original
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      else //  pLab > 100 // my
      {
        hnXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
      }
      fTotalXsc = hnXsc;
    }
    if( proton )    // pi- p
    {
      if( pLab < 0.28 ) 
      {
        hpXsc       = 0.288/((pLab - 0.28)*(pLab - 0.28) + 0.004);
        fElasticXsc = 1.8/((logP + 1.273)*(logP + 1.273) + 0.07);
      }
      else if( pLab < 0.395676 ) // first peak
      {
        hpXsc       = 0.648/((pLab - 0.28)*(pLab - 0.28) + 0.009);
        fElasticXsc = 0.257/((pLab - 0.28)*(pLab - 0.28) + 0.01);
       }
      else if( pLab < 0.5 )
      {
        hpXsc       = 26 + 110*(std::log(pLab/0.48))*(std::log(pLab/0.48));
        fElasticXsc = 0.37*hpXsc;
      }
      else if( pLab < 0.65 )
      {
        hpXsc       = 26 + 110*(std::log(pLab/0.48))*(std::log(pLab/0.48));
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.72 )
      {
        hpXsc = 36.1+
                10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.88 )
      {
        hpXsc = 36.1+
                10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 1.03 )
      {
        hpXsc = 36.1+
                10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 2.0 + 0.4/((pLab - 1.03)*(pLab - 1.03) + 0.016);
      }
      else if( pLab < 1.15 )
      {
        hpXsc = 36.1+
                10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 2.0 + 0.4/((pLab - 1.03)*(pLab - 1.03) + 0.016);
      }
      else if( pLab < 1.3 )
      {
        hpXsc = 36.1+
                10*std::exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*std::exp(-(pLab-1.015)*(pLab-1.015)/0.075/0.075);
        fElasticXsc = 3. + 13./pLab;
      }
      else if( pLab < 2.6 ) // < 3.0) // ns original
      {
        hpXsc = 36.1+0.079-4.313*std::log(pLab)+
                3*std::exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*std::exp(-(pLab-1.4)*(pLab-1.4)/0.12/0.12);
        fElasticXsc = 3. +13./pLab; // *std::log(pLab*6.79);
      }
      else   // mb
      {
        hpXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3. + 13./pLab;
      }
      fTotalXsc = hpXsc;
    }
  } 
  else if( theParticle == theKPlus && pORn ) 
  {
    if( proton )
    {
      xsection  = 17.91 + B*std::pow(std::log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) - 13.45*std::pow(sMand,-eta2);
    }
    if( neutron )
    {
      xsection = 17.87 + B*std::pow(std::log(sMand/s0),2.) 
                          + 5.17*std::pow(sMand,-eta1) - 7.23*std::pow(sMand,-eta2);
    }
    fTotalXsc = xsection;
  } 
  else if( theParticle == theKMinus && pORn ) 
  {
    if( proton )
    {
      xsection  = 17.91 + B*std::pow(std::log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) + 13.45*std::pow(sMand,-eta2);
    }
    if( neutron )
    {
      xsection = 17.87 + B*std::pow(std::log(sMand/s0),2.) 
                          + 5.17*std::pow(sMand,-eta1) + 7.23*std::pow(sMand,-eta2);
    }
    fTotalXsc = xsection;
  }
  else if( theParticle == theSMinus && pORn ) 
  {
    xsection  = 35.20 + B*std::pow(std::log(sMand/s0),2.) 
                          - 199.*std::pow(sMand,-eta1) + 264.*std::pow(sMand,-eta2);
  } 
  else if( theParticle == theGamma && pORn ) // modify later on
  {
    xsection  = 0.0 + B*std::pow(std::log(sMand/s0),2.) 
                          + 0.032*std::pow(sMand,-eta1) - 0.0*std::pow(sMand,-eta2);
    fTotalXsc = xsection;   
  } 
  else  // other then p,n,pi+,pi-,K+,K- as proton ??? 
  {
    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);
    }
    fTotalXsc = xsection;
  } 
  fTotalXsc   *= millibarn; // parametrised in mb
  fElasticXsc *= millibarn; // parametrised in mb

  if( proton && aParticle->GetDefinition()->GetPDGCharge() > 0. )
  {
    G4double cB = GetCoulombBarrier(aParticle, nucleon);
    fTotalXsc   *= cB;
    fElasticXsc *= cB; 
  }
  fInelasticXsc = fTotalXsc - fElasticXsc;
  if( fInelasticXsc < 0. ) fInelasticXsc = 0.;

  // G4cout<<fTotalXsc/millibarn<<"; "<<fElasticXsc/millibarn<<"; "<<fInelasticXsc/millibarn<<G4endl;

  return fTotalXsc;
}

G4double G4HadronNucleonXsc::GetHadronNucleonXscPDG	(	const G4DynamicParticle *	,
		const G4ParticleDefinition *
	)

Definition at line 219 of file G4HadronNucleonXsc.cc.

References CalcMandelstamS(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentum(), neutron, G4InuclParticleNames::proton, and G4InuclParticleNames::s0.

Referenced by G4BGGPionInelasticXS::BuildPhysicsTable(), G4BGGPionElasticXS::BuildPhysicsTable(), G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4BGGNucleonElasticXS::BuildPhysicsTable(), G4NeutronInelasticXS::GetElementCrossSection(), G4NeutronElasticXS::GetElementCrossSection(), GetHadronNucleonXscNS(), G4BGGPionInelasticXS::GetIsoCrossSection(), G4BGGPionElasticXS::GetIsoCrossSection(), G4BGGNucleonInelasticXS::GetIsoCrossSection(), and G4BGGNucleonElasticXS::GetIsoCrossSection().

{
  G4double xsection(0);
  G4int Zt=1, Nt=1, At=1;

   G4double 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

  // 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* theParticle = aParticle->GetDefinition();

  G4bool pORn = (nucleon == theProton || nucleon == theNeutron  );  
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);
  
  if(theParticle == theNeutron) // proton-neutron fit 
  {
    if ( proton )
    {
      xsection = Zt*( 35.80 + B*std::pow(std::log(sMand/s0),2.) 
                 + 40.15*std::pow(sMand,-eta1) - 30.*std::pow(sMand,-eta2));// on p
    }
    if ( neutron )
    {
      xsection  = Nt*( 35.45 + B*std::pow(std::log(sMand/s0),2.) 
                      + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2)); // on n pp for nn
    }
  } 
  else if(theParticle == theProton) 
  {
    if ( proton )
    {      
      xsection  = Zt*( 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 = Nt*( 35.80 + B*std::pow(std::log(sMand/s0),2.) 
                          + 40.15*std::pow(sMand,-eta1) - 30.*std::pow(sMand,-eta2));
    }
  } 
  else if(theParticle == theAProton) 
  {
    if ( proton )
    {      
      xsection  = Zt*( 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 = Nt*( 35.80 + B*std::pow(std::log(sMand/s0),2.) 
                          + 40.15*std::pow(sMand,-eta1) + 30.*std::pow(sMand,-eta2));
    }
  } 
  else if(theParticle == thePiPlus && pORn ) 
  {
    xsection  = At*( 20.86 + B*std::pow(std::log(sMand/s0),2.) 
                          + 19.24*std::pow(sMand,-eta1) - 6.03*std::pow(sMand,-eta2));
  } 
  else if(theParticle == thePiMinus && pORn ) 
  {
    xsection  = At*( 20.86 + B*std::pow(std::log(sMand/s0),2.) 
                          + 19.24*std::pow(sMand,-eta1) + 6.03*std::pow(sMand,-eta2));
  } 
  else if(theParticle == theKPlus) 
  {
    if ( proton )
    {      
      xsection  = Zt*( 17.91 + B*std::pow(std::log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) - 13.45*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 17.87 + B*std::pow(std::log(sMand/s0),2.) 
                          + 5.17*std::pow(sMand,-eta1) - 7.23*std::pow(sMand,-eta2));
    }
  } 
  else if(theParticle == theKMinus) 
  {
    if ( proton )
    {      
      xsection  = Zt*( 17.91 + B*std::pow(std::log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) + 13.45*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 17.87 + B*std::pow(std::log(sMand/s0),2.) 
                          + 5.17*std::pow(sMand,-eta1) + 7.23*std::pow(sMand,-eta2) );
    }
  }
  else if(theParticle == theSMinus && pORn ) 
  {
    xsection  = At*( 35.20 + B*std::pow(std::log(sMand/s0),2.) 
                          - 199.*std::pow(sMand,-eta1) + 264.*std::pow(sMand,-eta2) );
  } 
  else if(theParticle == theGamma && pORn ) // modify later on
  {
    xsection  = At*( 0.0 + B*std::pow(std::log(sMand/s0),2.) 
                          + 0.032*std::pow(sMand,-eta1) - 0.0*std::pow(sMand,-eta2) );
   
  } 
  else  // as proton ??? 
  {
    if ( proton )
    {      
      xsection  = Zt*( 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 = Nt*( 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

  fTotalXsc     = xsection;
  fInelasticXsc = 0.75*xsection;
  fElasticXsc   = fTotalXsc - fInelasticXsc;
  if (fElasticXsc < 0.) fElasticXsc = 0.;

  return xsection;
}

G4double G4HadronNucleonXsc::GetHadronNucleonXscVU	(	const G4DynamicParticle *	,
		const G4ParticleDefinition *
	)

Definition at line 961 of file G4HadronNucleonXsc.cc.

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetMomentum(), G4ParticleDefinition::GetPDGEncoding(), G4DynamicParticle::GetTotalEnergy(), neutron, and G4InuclParticleNames::proton.

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

  G4bool pORn = (nucleon == theProton || nucleon == theNeutron  );  
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);

   
  if( absPDGcode > 1000 && pORn )  //------Projectile is baryon -
  {
    if(proton)
    {
      fTotalXsc   = 48.0 +  0. *std::pow(Plab, 0.  ) + 0.522*sqrLogPlab - 4.51*LogPlab;
      fElasticXsc = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 47.3 +  0. *std::pow(Plab, 0.  ) + 0.513*sqrLogPlab - 4.27*LogPlab;
      fElasticXsc = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
    }
  }
  else if( PDGcode ==  211  && pORn )  //------Projectile is PionPlus ----
  {
    if(proton)
    {
      fTotalXsc  = 16.4 + 19.3 *std::pow(Plab,-0.42) + 0.19 *sqrLogPlab - 0.0 *LogPlab;
      fElasticXsc =  0.0 + 11.4*std::pow(Plab,-0.40) + 0.079*sqrLogPlab - 0.0 *LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   =  33.0 + 14.0 *std::pow(Plab,-1.36) + 0.456*sqrLogPlab - 4.03*LogPlab;
      fElasticXsc = 1.76 + 11.2*std::pow(Plab,-0.64) + 0.043*sqrLogPlab - 0.0 *LogPlab;
    }
  }
  else if( PDGcode == -211  && pORn )  //------Projectile is PionMinus ----
  {
    if(proton)
    {
      fTotalXsc   = 33.0 + 14.0 *std::pow(Plab,-1.36) + 0.456*sqrLogPlab - 4.03*LogPlab;
      fElasticXsc = 1.76 + 11.2*std::pow(Plab,-0.64) + 0.043*sqrLogPlab - 0.0 *LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 16.4 + 19.3 *std::pow(Plab,-0.42) + 0.19 *sqrLogPlab - 0.0 *LogPlab;
      fElasticXsc =  0.0 + 11.4*std::pow(Plab,-0.40) + 0.079*sqrLogPlab - 0.0 *LogPlab;
    }
  }
  else if( PDGcode ==  111  && pORn )  //------Projectile is PionZero  --
  {
    if(proton)
    {
      fTotalXsc   = (16.4 + 19.3 *std::pow(Plab,-0.42) + 0.19 *sqrLogPlab - 0.0 *LogPlab +   //Pi+
                        33.0 + 14.0 *std::pow(Plab,-1.36) + 0.456*sqrLogPlab - 4.03*LogPlab)/2; //Pi-

      fElasticXsc = ( 0.0 + 11.4*std::pow(Plab,-0.40) + 0.079*sqrLogPlab - 0.0 *LogPlab +    //Pi+
                         1.76 + 11.2*std::pow(Plab,-0.64) + 0.043*sqrLogPlab - 0.0 *LogPlab)/2; //Pi-

    }
    if(neutron)
    {    
      fTotalXsc   = (33.0 + 14.0 *std::pow(Plab,-1.36) + 0.456*sqrLogPlab - 4.03*LogPlab +   //Pi+
                        16.4 + 19.3 *std::pow(Plab,-0.42) + 0.19 *sqrLogPlab - 0.0 *LogPlab)/2; //Pi-
      fElasticXsc = ( 1.76 + 11.2*std::pow(Plab,-0.64) + 0.043*sqrLogPlab - 0.0 *LogPlab +   //Pi+
                         0.0  + 11.4*std::pow(Plab,-0.40) + 0.079*sqrLogPlab - 0.0 *LogPlab)/2; //Pi-
    }
  }
  else if( PDGcode == 321  && pORn )    //------Projectile is KaonPlus --
  {
    if(proton)
    {
      fTotalXsc   = 18.1 +  0. *std::pow(Plab, 0.  ) + 0.26 *sqrLogPlab - 1.0 *LogPlab;
      fElasticXsc =  5.0 +  8.1*std::pow(Plab,-1.8 ) + 0.16 *sqrLogPlab - 1.3 *LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 18.7 +  0. *std::pow(Plab, 0.  ) + 0.21 *sqrLogPlab - 0.89*LogPlab;
      fElasticXsc =  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab;
    }
  }
  else if( PDGcode ==-321  && pORn )  //------Projectile is KaonMinus ----
  {
    if(proton)
    {
      fTotalXsc   = 32.1 +  0. *std::pow(Plab, 0.  ) + 0.66*sqrLogPlab - 5.6*LogPlab;
      fElasticXsc =  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29*sqrLogPlab - 2.4*LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 25.2 +  0. *std::pow(Plab, 0.  ) + 0.38*sqrLogPlab - 2.9*LogPlab;
      fElasticXsc =  5.0 +  8.1*std::pow(Plab,-1.8 ) + 0.16*sqrLogPlab - 1.3*LogPlab;
    }
  }
  else if( PDGcode == 311  && pORn )  //------Projectile is KaonZero -----
  {
    if(proton)
    {
      fTotalXsc   = ( 18.1 +  0. *std::pow(Plab, 0.  ) + 0.26 *sqrLogPlab - 1.0 *LogPlab +   //K+
                        32.1 +  0. *std::pow(Plab, 0.  ) + 0.66 *sqrLogPlab - 5.6 *LogPlab)/2; //K-
      fElasticXsc = (  5.0 +  8.1*std::pow(Plab,-1.8 ) + 0.16 *sqrLogPlab - 1.3 *LogPlab +   //K+
                         7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab)/2; //K-
    }
    if(neutron)
    {    
      fTotalXsc   = ( 18.7 +  0. *std::pow(Plab, 0.  ) + 0.21 *sqrLogPlab - 0.89*LogPlab +   //K+
                         25.2 +  0. *std::pow(Plab, 0.  ) + 0.38 *sqrLogPlab - 2.9 *LogPlab)/2; //K-
      fElasticXsc = (  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab +   //K+
                         5.0 +  8.1*std::pow(Plab,-1.8 ) + 0.16 *sqrLogPlab - 1.3 *LogPlab)/2; //K-
    }
  }
  else  //------Projectile is undefined, Nucleon assumed
  {
    if(proton)
    {
      fTotalXsc   = 48.0 +  0. *std::pow(Plab, 0.  ) + 0.522*sqrLogPlab - 4.51*LogPlab;
      fElasticXsc = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 47.3 +  0. *std::pow(Plab, 0.  ) + 0.513*sqrLogPlab - 4.27*LogPlab;
      fElasticXsc = 11.9 + 26.9*std::pow(Plab,-1.21) + 0.169*sqrLogPlab - 1.85*LogPlab;
    }
  }
  fTotalXsc   *= millibarn;
  fElasticXsc *= millibarn;
  fInelasticXsc   = fTotalXsc - fElasticXsc;
  if (fInelasticXsc < 0.) fInelasticXsc = 0.;

  return fTotalXsc;    
}

G4double G4HadronNucleonXsc::GetInelasticHadronNucleonXsc

(

)

[inline]

Definition at line 90 of file G4HadronNucleonXsc.hh.

Referenced by G4BGGPionInelasticXS::BuildPhysicsTable(), G4BGGNucleonInelasticXS::BuildPhysicsTable(), G4NeutronInelasticXS::GetElementCrossSection(), G4GlauberGribovCrossSection::GetIsoCrossSection(), G4ComponentGGHadronNucleusXsc::GetIsoCrossSection(), G4BGGPionInelasticXS::GetIsoCrossSection(), G4BGGNucleonInelasticXS::GetIsoCrossSection(), G4GGNuclNuclCrossSection::GetZandACrossSection(), and G4ComponentGGNuclNuclXsc::GetZandACrossSection().

{ return fInelasticXsc; }; 

G4double G4HadronNucleonXsc::GetKmNeutronTotXscVector

(

G4double

logEnergy

)

[inline]

Definition at line 97 of file G4HadronNucleonXsc.hh.

References G4PhysicsVector::Value().

Referenced by G4GlauberGribovCrossSection::GetKaonNucleonXscVector(), and G4ComponentGGHadronNucleusXsc::GetKaonNucleonXscVector().

{ return fKmNeutronTotXscVector.Value(logEnergy); };

G4double G4HadronNucleonXsc::GetKmProtonTotXscVector

(

G4double

logEnergy

)

[inline]

Definition at line 96 of file G4HadronNucleonXsc.hh.

References G4PhysicsVector::Value().

Referenced by G4GlauberGribovCrossSection::GetKaonNucleonXscVector(), and G4ComponentGGHadronNucleusXsc::GetKaonNucleonXscVector().

{ return fKmProtonTotXscVector.Value(logEnergy); };

G4double G4HadronNucleonXsc::GetKpNeutronTotXscVector

(

G4double

logEnergy

)

[inline]

Definition at line 95 of file G4HadronNucleonXsc.hh.

References G4PhysicsVector::Value().

Referenced by G4GlauberGribovCrossSection::GetKaonNucleonXscVector(), and G4ComponentGGHadronNucleusXsc::GetKaonNucleonXscVector().

{ return fKpNeutronTotXscVector.Value(logEnergy); };

G4double G4HadronNucleonXsc::GetKpProtonTotXscVector

(

G4double

logEnergy

)

[inline]

Definition at line 94 of file G4HadronNucleonXsc.hh.

References G4PhysicsVector::Value().

Referenced by G4GlauberGribovCrossSection::GetKaonNucleonXscVector(), and G4ComponentGGHadronNucleusXsc::GetKaonNucleonXscVector().

{ return fKpProtonTotXscVector.Value(logEnergy); };

G4double G4HadronNucleonXsc::GetTotalHadronNucleonXsc

(

)

[inline]

Definition at line 88 of file G4HadronNucleonXsc.hh.

{ return fTotalXsc;     }; 

void G4HadronNucleonXsc::InitialiseKaonNucleonTotXsc

(

)

Definition at line 1192 of file G4HadronNucleonXsc.cc.

References G4LPhysicsFreeVector::PutValues(), and G4PhysicsVector::SetSpline().

Referenced by G4HadronNucleonXsc().

{
  G4int i = 0, iMax;
  G4double tmpxsc[106];

  // Kp-proton tot xsc

  iMax = 66;
  fKpProtonTotXscVector = G4LPhysicsFreeVector(iMax, fKpProtonTotTkin[0], fKpProtonTotTkin[iMax-1]);
  fKpProtonTotXscVector.SetSpline(true);

  for( i = 0; i < iMax; i++)
  {
    tmpxsc[i] = 0.;

    if( i == 0 || i == iMax-1 ) tmpxsc[i] = fKpProtonTotXsc[i];
    else                        tmpxsc[i] = 0.5*(fKpProtonTotXsc[i-1]+fKpProtonTotXsc[i+1]);

    fKpProtonTotXscVector.PutValues(size_t(i), fKpProtonTotTkin[i], tmpxsc[i]*millibarn);
  }

  // Kp-neutron tot xsc

  iMax = 75;
  fKpNeutronTotXscVector = G4LPhysicsFreeVector(iMax, fKpNeutronTotTkin[0], fKpNeutronTotTkin[iMax-1]);
  fKpNeutronTotXscVector.SetSpline(true);

  for( i = 0; i < iMax; i++)
  {
    tmpxsc[i] = 0.;
    if( i == 0 || i == iMax-1 ) tmpxsc[i] = fKpNeutronTotXsc[i];
    else                        tmpxsc[i] = 0.5*(fKpNeutronTotXsc[i-1]+fKpNeutronTotXsc[i+1]);

    fKpNeutronTotXscVector.PutValues(size_t(i), fKpNeutronTotTkin[i], tmpxsc[i]*millibarn);
  }

  // Km-proton tot xsc

  iMax = 106;
  fKmProtonTotXscVector = G4LPhysicsFreeVector(iMax, fKmProtonTotTkin[0], fKmProtonTotTkin[iMax-1]);
  fKmProtonTotXscVector.SetSpline(true);

  for( i = 0; i < iMax; i++)
  {
    tmpxsc[i] = 0.;

    if( i == 0 || i == iMax-1 ) tmpxsc[i] = fKmProtonTotXsc[i];
    else                        tmpxsc[i] = 0.5*(fKmProtonTotXsc[i-1]+fKmProtonTotXsc[i+1]);

    fKmProtonTotXscVector.PutValues(size_t(i), fKmProtonTotTkin[i], tmpxsc[i]*millibarn);
  }

  // Km-neutron tot xsc

  iMax = 68;
  fKmNeutronTotXscVector = G4LPhysicsFreeVector(iMax, fKmNeutronTotTkin[0], fKmNeutronTotTkin[iMax-1]);
  fKmNeutronTotXscVector.SetSpline(true);

  for( i = 0; i < iMax; i++)
  {
    tmpxsc[i] = 0.;
    if( i == 0 || i == iMax-1 ) tmpxsc[i] = fKmNeutronTotXsc[i];
    else                        tmpxsc[i] = 0.5*(fKmNeutronTotXsc[i-1]+fKmNeutronTotXsc[i+1]);

    fKmNeutronTotXscVector.PutValues(size_t(i), fKmNeutronTotTkin[i], tmpxsc[i]*millibarn);
  }
}

G4bool G4HadronNucleonXsc::IsApplicable	(	const G4DynamicParticle *	aDP,
		const G4Element *
	)			[virtual]

Definition at line 93 of file G4HadronNucleonXsc.cc.

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

{
  G4int Z = G4lrint(anElement->GetZ());
  G4int A = G4lrint(anElement->GetN());
  return IsIsoApplicable(aDP, Z, A);
} 

G4bool G4HadronNucleonXsc::IsIsoApplicable	(	const G4DynamicParticle *	aDP,
		G4int	Z,
		G4int	A
	)			[virtual]

Definition at line 105 of file G4HadronNucleonXsc.cc.

References G4DynamicParticle::GetDefinition(), and G4DynamicParticle::GetKineticEnergy().

Referenced by IsApplicable().

{
  G4bool applicable = false;
  // G4int baryonNumber     = aDP->GetDefinition()->GetBaryonNumber();
  G4double kineticEnergy = aDP->GetKineticEnergy();

  const G4ParticleDefinition* theParticle = aDP->GetDefinition();
 
  if ( ( kineticEnergy  >= fLowerLimit &&
         Z > 1 &&      // >=  He
       ( theParticle == theAProton   ||
         theParticle == theGamma     ||
         theParticle == theKPlus     ||
         theParticle == theKMinus    || 
         theParticle == theSMinus)      )    ||  

       ( kineticEnergy  >= 0.1*fLowerLimit &&
         Z > 1 &&      // >=  He
       ( theParticle == theProton    ||
         theParticle == theNeutron   ||   
         theParticle == thePiPlus    ||
         theParticle == thePiMinus       ) )    ) applicable = true;

  return applicable;
}

---

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

• G4HadronNucleonXsc.hh
• G4HadronNucleonXsc.cc

---