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

(

)

Definition at line 39 of file G4HadronNucleonXsc.cc.

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

virtual

Definition at line 79 of file G4HadronNucleonXsc.cc.

{}

Member Function Documentation

CalcMandelstamS()

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

Definition at line 1125 of file G4HadronNucleonXsc.cc.

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

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

CalculateEcmValue()

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

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

CrossSectionDescription()

void G4HadronNucleonXsc::CrossSectionDescription

(

std::ostream &

outFile

)

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

DumpPhysicsTable()

virtual void G4HadronNucleonXsc::DumpPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Definition at line 65 of file G4HadronNucleonXsc.hh.

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

GetCoulombBarrier()

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

Definition at line 1140 of file G4HadronNucleonXsc.cc.

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

Referenced by GetHadronNucleonXscNS().

GetElasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetElasticHadronNucleonXsc ( )

inline

Definition at line 89 of file G4HadronNucleonXsc.hh.

{ return fElasticXsc;   }; 

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

GetHadronNucleonXscEL()

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

Definition at line 141 of file G4HadronNucleonXsc.cc.

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

GetHadronNucleonXscNS()

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

pi- ////////////////////////////////////////////

Definition at line 363 of file G4HadronNucleonXsc.cc.

{
  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 ) ////// proton //////////////////////////////////////////////
  {
    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 ) /////////////////// p_bar ///////////////////////////
  {
    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 ) /// pi- ////////////////////////////////////////////
  {
    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;
}

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

GetHadronNucleonXscPDG()

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

Definition at line 219 of file G4HadronNucleonXsc.cc.

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

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

GetHadronNucleonXscVU()

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

Definition at line 961 of file G4HadronNucleonXsc.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;
 
  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;    
}

GetInelasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetInelasticHadronNucleonXsc ( )

inline

Definition at line 90 of file G4HadronNucleonXsc.hh.

{ return fInelasticXsc; }; 

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

GetKmNeutronTotXscVector()

G4double G4HadronNucleonXsc::GetKmNeutronTotXscVector ( G4double  logEnergy )

inline

Definition at line 97 of file G4HadronNucleonXsc.hh.

{ return fKmNeutronTotXscVector.Value(logEnergy); };

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

GetKmProtonTotXscVector()

G4double G4HadronNucleonXsc::GetKmProtonTotXscVector ( G4double  logEnergy )

inline

Definition at line 96 of file G4HadronNucleonXsc.hh.

{ return fKmProtonTotXscVector.Value(logEnergy); };

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

GetKpNeutronTotXscVector()

G4double G4HadronNucleonXsc::GetKpNeutronTotXscVector ( G4double  logEnergy )

inline

Definition at line 95 of file G4HadronNucleonXsc.hh.

{ return fKpNeutronTotXscVector.Value(logEnergy); };

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

GetKpProtonTotXscVector()

G4double G4HadronNucleonXsc::GetKpProtonTotXscVector ( G4double  logEnergy )

inline

Definition at line 94 of file G4HadronNucleonXsc.hh.

{ return fKpProtonTotXscVector.Value(logEnergy); };

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

GetTotalHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetTotalHadronNucleonXsc ( )

inline

Definition at line 88 of file G4HadronNucleonXsc.hh.

{ return fTotalXsc;     }; 

InitialiseKaonNucleonTotXsc()

void G4HadronNucleonXsc::InitialiseKaonNucleonTotXsc

(

)

Definition at line 1192 of file G4HadronNucleonXsc.cc.

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

Referenced by G4HadronNucleonXsc().

IsApplicable()

G4bool G4HadronNucleonXsc::IsApplicable ( const G4DynamicParticle *  aDP, const G4Element *  anElement  )

virtual

Definition at line 93 of file G4HadronNucleonXsc.cc.

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

IsIsoApplicable()

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

virtual

Definition at line 105 of file G4HadronNucleonXsc.cc.

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

Referenced by IsApplicable().

---

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

• G4HadronNucleonXsc.hh
• G4HadronNucleonXsc.cc