G4HadronNucleonXsc Class Reference

#include <G4HadronNucleonXsc.hh>

Collaboration diagram for G4HadronNucleonXsc:

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

Private Attributes
const G4double	fLowerLimit
G4double	fTotalXsc
G4double	fElasticXsc
G4double	fInelasticXsc
G4ParticleDefinition *	theGamma
G4ParticleDefinition *	theProton
G4ParticleDefinition *	theNeutron
G4ParticleDefinition *	theAProton
G4ParticleDefinition *	theANeutron
G4ParticleDefinition *	thePiPlus
G4ParticleDefinition *	thePiMinus
G4ParticleDefinition *	thePiZero
G4ParticleDefinition *	theKPlus
G4ParticleDefinition *	theKMinus
G4ParticleDefinition *	theK0S
G4ParticleDefinition *	theK0L
G4ParticleDefinition *	theL
G4ParticleDefinition *	theAntiL
G4ParticleDefinition *	theSPlus
G4ParticleDefinition *	theASPlus
G4ParticleDefinition *	theSMinus
G4ParticleDefinition *	theASMinus
G4ParticleDefinition *	theS0
G4ParticleDefinition *	theAS0
G4ParticleDefinition *	theXiMinus
G4ParticleDefinition *	theXi0
G4ParticleDefinition *	theAXiMinus
G4ParticleDefinition *	theAXi0
G4ParticleDefinition *	theOmega
G4ParticleDefinition *	theAOmega
G4ParticleDefinition *	theD
G4ParticleDefinition *	theT
G4ParticleDefinition *	theA
G4ParticleDefinition *	theHe3

Detailed Description

Definition at line 51 of file G4HadronNucleonXsc.hh.

Constructor & Destructor Documentation

G4HadronNucleonXsc()

G4HadronNucleonXsc::G4HadronNucleonXsc

(

)

Definition at line 41 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();
}

Here is the call graph for this function:

~G4HadronNucleonXsc()

G4HadronNucleonXsc::~G4HadronNucleonXsc ( )

virtual

Definition at line 83 of file G4HadronNucleonXsc.cc.

{}

Member Function Documentation

CalcMandelstamS()

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

Definition at line 1418 of file G4HadronNucleonXsc.cc.

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

  return sMand;
}

Here is the caller graph for this function:

CalculateEcmValue()

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

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

Here is the caller graph for this function:

CrossSectionDescription()

void G4HadronNucleonXsc::CrossSectionDescription

(

std::ostream &

outFile

)

const

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

Here is the caller graph for this function:

DumpPhysicsTable()

virtual void G4HadronNucleonXsc::DumpPhysicsTable ( const G4ParticleDefinition &  )

inlinevirtual

Definition at line 65 of file G4HadronNucleonXsc.hh.

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

Here is the call graph for this function:

GetCoulombBarrier()

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetElasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetElasticHadronNucleonXsc ( )

inline

Definition at line 91 of file G4HadronNucleonXsc.hh.

{ return fElasticXsc;   }; 

Here is the caller graph for this function:

GetHadronNucleonXscEL()

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

Definition at line 145 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*G4Log(sMand)*G4Log(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*G4Log(sMand)*G4Log(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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetHadronNucleonXscNS()

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

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

Definition at line 367 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 = G4Log(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;
  G4double minLogP = 3.5;       // min of (lnP-minLogP)^2 
  G4double cofLogE = .0557;     // elastic (lnP-minLogP)^2 
  G4double cofLogT = .3;        // total (lnP-minLogP)^2 
  G4double pMin = .1;        // fast LE calculation 
  G4double pMax = 1000.;     // fast HE calculation 


  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(G4Log(sMand/400.),1.65) + 9.19*std::pow(sMand,-0.458);
    
      fTotalXsc = xsection;
    
    }
    else if( pLab >= 100.)
    {
      B0 = 7.5;
      A0 = 100. - B0*G4Log(3.0e7);

      xsection = A0 + B0*G4Log(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(G4Log(sMand/28.9),1.1) + 9.19*std::pow(sMand,-0.458);
      
      fTotalXsc = xsection;
    }
    else if( pLab >= 10.)
    {
        B0 = 7.5;
        A0 = 100. - B0*G4Log(3.0e7);

        xsection = A0 + B0*G4Log(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( G4Log(0.73/pLab), 3.5 ) );
          fElasticXsc = hnXsc;
        }
        else if( pLab < 0.73 )
        {
          hnXsc = 23 + 50*( std::pow( G4Log(0.73/pLab), 3.5 ) );
          fElasticXsc = hnXsc; 
        }
        else if( pLab < 1.05  )
        {
          hnXsc = 23 + 40*(G4Log(pLab/0.73))*
                         (G4Log(pLab/0.73));
          fElasticXsc = 23 + 20*(G4Log(pLab/0.73))*
                         (G4Log(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(G4Log(1.3/pLab),3.6); // was as pLab < 0.8
      fElasticXsc = hpXsc;
        }      
        else if( pLab < 0.8 )
        {
          hpXsc = 33+30*std::pow(G4Log(pLab/1.3),4.0);
      fElasticXsc = hpXsc;
        }      
        else if( pLab < 1.05 )
        {
          hpXsc = 33+30*std::pow(G4Log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( G4Log(0.511/pLab)*G4Log(0.511/pLab) + 1.6 );
        }
        else if( pLab < 1.4 )
        {
          hpXsc = 33+30*std::pow(G4Log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( G4Log(0.511/pLab)*G4Log(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(G4Log(sMand/400.),1.65) + 9.19*std::pow(sMand,-0.458);
     
      fTotalXsc = xsection;
    }
    else if( pLab >= 100.)
    {
      B0 = 7.5;
      A0 = 100. - B0*G4Log(3.0e7);

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

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

      xsection = A0 + B0*G4Log(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( G4Log(0.73/pLab), 3.5 ) );
          fElasticXsc = hpXsc;
        }
        else if( pLab < 0.73 )
        {
          hpXsc = 23 + 50*( std::pow( G4Log(0.73/pLab), 3.5 ) );
          fElasticXsc = hpXsc; 
        }
        else if( pLab < 1.05  )
        {
          hpXsc = 23 + 40*(G4Log(pLab/0.73))*
                         (G4Log(pLab/0.73));
          fElasticXsc = 23 + 20*(G4Log(pLab/0.73))*
                         (G4Log(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(G4Log(1.3/pLab),3.6); // was as pLab < 0.8
      fElasticXsc = hnXsc;
        }      
        else if( pLab < 0.8 )
        {
          hnXsc = 33+30*std::pow(G4Log(pLab/1.3),4.0);
      fElasticXsc = hnXsc;
        }      
        else if( pLab < 1.05 )
        {
          hnXsc = 33+30*std::pow(G4Log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( G4Log(0.511/pLab)*G4Log(0.511/pLab) + 1.6 );
        }
        else if( pLab < 1.4 )
        {
          hnXsc = 33+30*std::pow(G4Log(pLab/0.95),2.0);
          fElasticXsc =  6 + 52/( G4Log(0.511/pLab)*G4Log(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(G4Log(sMand/s0),2.) 
                          + 42.53*std::pow(sMand,-eta1) + 33.34*std::pow(sMand,-eta2);
    }
    if( neutron ) // ???
    {
      xsection = 35.80 + B*std::pow(G4Log(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*(G4Log(pLab/0.77))*(G4Log(pLab/0.77));
        hpXsc        = Ex4 + 14.9;
        fElasticXsc = hpXsc*G4Exp(-3.*(pLab - 0.68));  
      }
      else if( pLab < 1.15 )
      {
        G4double Ex4 = 88*(G4Log(pLab/0.77))*(G4Log(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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.*G4Log(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*(G4Log(pLab/0.48))*(G4Log(pLab/0.48));
        fElasticXsc = 0.37*hnXsc;
      }
      else if( pLab < 0.65 )
      {
        hnXsc       = 26 + 110*(G4Log(pLab/0.48))*(G4Log(pLab/0.48));
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.72 )
      {
        hnXsc = 36.1 + 10*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Log(pLab)+
                3*G4Exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*G4Exp(-(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*G4Log(pLab)+
                3*G4Exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*G4Exp(-(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*(G4Log(pLab/0.77))*(G4Log(pLab/0.77));
        hnXsc        = Ex4 + 14.9;
        fElasticXsc = hnXsc*G4Exp(-3.*(pLab - 0.68));  
      }
      else if( pLab < 1.15 )
      {
        G4double Ex4 = 88*(G4Log(pLab/0.77))*(G4Log(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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*G4Exp(-(pLab-2.55)*(pLab-2.55)/0.55/0.55);
        G4double Ex2 = 12*G4Exp(-(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.*G4Log(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*(G4Log(pLab/0.48))*(G4Log(pLab/0.48));
        fElasticXsc = 0.37*hpXsc;
      }
      else if( pLab < 0.65 )
      {
        hpXsc       = 26 + 110*(G4Log(pLab/0.48))*(G4Log(pLab/0.48));
        fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
      }
      else if( pLab < 0.72 )
      {
        hpXsc = 36.1+
                10*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Exp(-(pLab-0.72)*(pLab-0.72)/0.06/0.06)+
                24*G4Exp(-(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*G4Log(pLab)+
                3*G4Exp(-(pLab-2.1)*(pLab-2.1)/0.4/0.4)+
                1.5*G4Exp(-(pLab-1.4)*(pLab-1.4)/0.12/0.12);
        fElasticXsc = 3. +13./pLab; // *G4Log(pLab*6.79);
      }
      else   // mb
      {
        hpXsc = GetHadronNucleonXscPDG(aParticle, nucleon )/millibarn;
        fElasticXsc = 3. + 13./pLab;
      }
      fTotalXsc = hpXsc;
    }
  } 
  else if( (theParticle == theKMinus || theParticle == theK0S) && proton )   // Kmp/K0p /////////////////////////////////
  {
    if( pLab < pMin)
    {
      G4double psp = pLab*std::sqrt(pLab);
      fElasticXsc  = 5.2/psp;
      fTotalXsc    = 14./psp;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double sp  = std::sqrt(pLab);
      G4double psp = pLab*sp;
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
      G4double lm  = pLab - .39;
      G4double md  = lm*lm + .000356;

      G4double lh1  = pLab - 0.78;
      G4double hd1  = lh1*lh1 + .00166;
 
      G4double lh  = pLab - 1.01;
      G4double hd  = lh*lh + .011;

      G4double lh2  = pLab - 1.63;
      G4double hd2  = lh2*lh2 + .007;

      fElasticXsc  = 5.2/psp + (1.1*cofLogE*ld2 + 2.23)/(1. - .7/sp + .075/p4) 
                 + .004/md + 0.005/hd1+ 0.01/hd2 +.15/hd; // small peaks were added

      fTotalXsc    = 14./psp + (1.1*cofLogT*ld2 + 19.5)/(1. - .21/sp + .52/p4) 
                 + .006/md  + 0.01/hd1+ 0.02/hd2 + .20/hd ;
    }
  }
  else if( (theParticle == theKMinus || theParticle == theK0S) && neutron )   // Kmn/K0n //////////////////////////////
  {
    if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
 
      G4double lh  = pLab - 0.98;
      G4double hd  = lh*lh + .021;

      G4double LogPlab = G4Log( pLab );
      G4double sqrLogPlab = LogPlab * LogPlab;

      fElasticXsc  = // 5.2/psp + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .075/p4) + .004/md 
                     5.0 +  8.1*std::pow(pLab,-1.8 ) + 0.16*sqrLogPlab - 1.3*LogPlab + .15/hd;
      fTotalXsc    = // 14./psp + 
                     //  (1.1*cofLogT*ld2 + 19.5)/(1. - .21/sp + .52/p4) 
    // WP                     25.2 +  0. *std::pow(pLab, 0.  ) + 0.38*sqrLogPlab - 2.9*LogPlab               
                     25.2 +  0.38*sqrLogPlab - 2.9*LogPlab               
                     //       + .006/md  + 0.01/hd1+ 0.02/hd2 
                        + 0.30/hd ;
    }
  }
  else if(  (theParticle == theKPlus || theParticle == theK0L) && proton  )  // Kpp/aKp ////////////////////////
  {
    if( pLab < pMin )
    {
      G4double lr = pLab - .38;
      G4double lm = pLab - 1.;
      G4double md = lm*lm + .392;   
      fElasticXsc = .7/(lr*lr + .076) + 2./md;
      fTotalXsc   = .7/(lr*lr + .076) + 2.6/md;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double lr  = pLab - .38;
      G4double LE  = .7/(lr*lr + .076);
      G4double sp  = std::sqrt(pLab);
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
      G4double lm  = pLab - 1.;
      G4double md  = lm*lm + .392;
      fElasticXsc  = LE + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) + 2./md;
      fTotalXsc    = LE + (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) + 2.6/md;
    }
  }
  else if(  (theParticle == theKPlus || theParticle == theK0L) && neutron  )  // Kpn/aKn ///////////////////////
  {
    if( pLab < pMin )
    {
      G4double lm = pLab - 0.94;
      G4double md = lm*lm + .392;   
      fElasticXsc = 2./md;
      fTotalXsc   = 4.6/md;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double sp  = std::sqrt(pLab);
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
      G4double lm  = pLab - 0.94;
      G4double md  = lm*lm + .392;
      fElasticXsc  = (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) + 2./md;
      fTotalXsc    = (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) + 4.6/md;
    }
  }
  else if( theParticle == theSMinus && pORn ) 
  {
    xsection  = 35.20 + B*std::pow(G4Log(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(G4Log(sMand/s0),2.) 
      + 0.032*std::pow(sMand,-eta1); // WP - 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(G4Log(sMand/s0),2.) 
                          + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2);
    }
    if( neutron )
    {
      xsection += 35.80 + B*std::pow(G4Log(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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetHadronNucleonXscPDG()

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

Definition at line 223 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(G4Log(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(G4Log(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(G4Log(sMand/s0),2.) 
                          + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 35.80 + B*std::pow(G4Log(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(G4Log(sMand/s0),2.) 
                          + 42.53*std::pow(sMand,-eta1) + 33.34*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 35.80 + B*std::pow(G4Log(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(G4Log(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(G4Log(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(G4Log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) - 13.45*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 17.87 + B*std::pow(G4Log(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(G4Log(sMand/s0),2.) 
                          + 7.14*std::pow(sMand,-eta1) + 13.45*std::pow(sMand,-eta2));
    }
    if ( neutron )
    {
      xsection = Nt*( 17.87 + B*std::pow(G4Log(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(G4Log(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(G4Log(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(G4Log(sMand/s0),2.) 
                       + 42.53*std::pow(sMand,-eta1) - 33.34*std::pow(sMand,-eta2) );
    }
    if ( neutron )
    {
      xsection = Nt*( 35.80 + B*std::pow(G4Log(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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetHadronNucleonXscVU()

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

Definition at line 1238 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 = G4Log( 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)
    {
      // WP      fTotalXsc   = 48.0 +  0. *std::pow(Plab, 0.  ) + 0.522*sqrLogPlab - 4.51*LogPlab;
      fTotalXsc   = 48.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)
    {    
      // WP      fTotalXsc   = 47.3 +  0. *std::pow(Plab, 0.  ) + 0.513*sqrLogPlab - 4.27*LogPlab;
      fTotalXsc   = 47.3  + 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)
    {
      // WP      fTotalXsc   = 18.1 +  0. *std::pow(Plab, 0.  ) + 0.26 *sqrLogPlab - 1.0 *LogPlab;
      fTotalXsc   = 18.1 +  0.26 *sqrLogPlab - 1.0 *LogPlab;
      fElasticXsc =  5.0 +  8.1*std::pow(Plab,-1.8 ) + 0.16 *sqrLogPlab - 1.3 *LogPlab;
    }
    if(neutron)
    {    
      // WP      fTotalXsc   = 18.7 +  0. *std::pow(Plab, 0.  ) + 0.21 *sqrLogPlab - 0.89*LogPlab;
      // WP      fElasticXsc =  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab;
      fTotalXsc   = 18.7  + 0.21 *sqrLogPlab - 0.89*LogPlab;
      fElasticXsc =  7.3  + 0.29 *sqrLogPlab - 2.4 *LogPlab;
    }
  }
  else if( PDGcode ==-321  && pORn )  //------Projectile is KaonMinus ----
  {
    if(proton)
    {
      // WP      fTotalXsc   = 32.1 +  0. *std::pow(Plab, 0.  ) + 0.66*sqrLogPlab - 5.6*LogPlab;
      // WP      fElasticXsc =  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29*sqrLogPlab - 2.4*LogPlab;
      fTotalXsc   = 32.1 + 0.66*sqrLogPlab - 5.6*LogPlab;
      fElasticXsc =  7.3 + 0.29*sqrLogPlab - 2.4*LogPlab;
    }
    if(neutron)
    {    
      // WP      fTotalXsc   = 25.2 +  0. *std::pow(Plab, 0.  ) + 0.38*sqrLogPlab - 2.9*LogPlab;
      fTotalXsc   = 25.2 + 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)
    {
      // WP      fTotalXsc   = ( 18.1 +  0. *std::pow(Plab, 0.  ) + 0.26 *sqrLogPlab - 1.0 *LogPlab +   //K+
      // WP                      32.1 +  0. *std::pow(Plab, 0.  ) + 0.66 *sqrLogPlab - 5.6 *LogPlab)/2; //K-
      fTotalXsc   = ( 18.1 + 0.26 *sqrLogPlab - 1.0 *LogPlab +   //K+
                      32.1 + 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.29 *sqrLogPlab - 2.4 *LogPlab)/2; //K-
      // WP                         7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab)/2; //K-
    }
    if(neutron)
    {    
      // WP      fTotalXsc   = ( 18.7 +  0. *std::pow(Plab, 0.  ) + 0.21 *sqrLogPlab - 0.89*LogPlab +   //K+
      // WP                      25.2 +  0. *std::pow(Plab, 0.  ) + 0.38 *sqrLogPlab - 2.9 *LogPlab)/2; //K-
      fTotalXsc   = ( 18.7 + 0.21 *sqrLogPlab - 0.89*LogPlab +   //K+
                      25.2 + 0.38 *sqrLogPlab - 2.9 *LogPlab)/2; //K-
      // WP      fElasticXsc = (  7.3 +  0. *std::pow(Plab,-0.  ) + 0.29 *sqrLogPlab - 2.4 *LogPlab +   //K+
      fElasticXsc = (  7.3 + 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)
    {
      // WP      fTotalXsc   = 48.0 +  0. *std::pow(Plab, 0.  ) + 0.522*sqrLogPlab - 4.51*LogPlab;
      fTotalXsc   = 48.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)
    {    
      // WP      fTotalXsc   = 47.3 +  0. *std::pow(Plab, 0.  ) + 0.513*sqrLogPlab - 4.27*LogPlab;
      fTotalXsc   = 47.3 + 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;    
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetInelasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetInelasticHadronNucleonXsc ( )

inline

Definition at line 92 of file G4HadronNucleonXsc.hh.

{ return fInelasticXsc; }; 

Here is the caller graph for this function:

GetKaonNucleonXscGG()

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

Definition at line 1072 of file G4HadronNucleonXsc.cc.

{
  G4double pLab = aParticle->GetMomentum().mag();

  pLab /= GeV;
  G4double LogPlab = G4Log( pLab );
  G4double sqrLogPlab = LogPlab * LogPlab;

  G4double minLogP = 3.5;       // min of (lnP-minLogP)^2 
  G4double cofLogE = .0557;     // elastic (lnP-minLogP)^2 
  G4double cofLogT = .3;        // total (lnP-minLogP)^2 
  G4double pMin = .1;        // fast LE calculation 
  G4double pMax = 1000.;     // fast HE calculation 

  const G4ParticleDefinition* theParticle = aParticle->GetDefinition();

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

  if(  (theParticle == theKMinus || theParticle == theK0S) && proton ) // (K-,K0)on p ////////////////////////////
  {

    if( pLab < pMin)
    {
      G4double psp = pLab*std::sqrt(pLab);
      fElasticXsc  = 5.2/psp;
      fTotalXsc    = 14./psp;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double sp  = std::sqrt(pLab);
      G4double psp = pLab*sp;
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
 
      G4double lh  = pLab - 0.98;
      G4double hd  = lh*lh + .045;


      fElasticXsc  = 5.2/psp + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .075/p4) // + .004/md 
               + .15/hd;
      fTotalXsc    = 14./psp + (1.1*cofLogT*ld2 + 19.5)/(1. - .21/sp + .52/p4) 
                  //  + .006/md  + 0.01/hd1 + 0.02/hd2 
                     + .60/hd;
    }
  }
  else if( (theParticle == theKMinus || theParticle == theK0S) && neutron )   // Kmn/K0n /////////////////////////////
  {
    if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
 
      G4double lh  = pLab - 0.98;
      G4double hd  = lh*lh + .045;

      fElasticXsc  = // 5.2/psp + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .075/p4) + .004/md 
                     5.0 +  8.1*std::pow(pLab,-1.8 ) + 0.16*sqrLogPlab - 1.3*LogPlab + .15/hd;
      fTotalXsc    = // 14./psp + 
                     //  (1.1*cofLogT*ld2 + 19.5)/(1. - .21/sp + .52/p4) 
    // WP                     25.2 +  0. *std::pow(pLab, 0.  ) + 0.38*sqrLogPlab - 2.9*LogPlab               
                     25.2 + 0.38*sqrLogPlab - 2.9*LogPlab                
                     //       + .006/md  + 0.01/hd1+ 0.02/hd2 
                        + 0.60/hd ;
    }
  }
  else if(  (theParticle == theKPlus || theParticle == theK0L) && proton )  // Kpp/aKp //////////////////////
  {
    if( pLab < pMin )
    {
      G4double lr = pLab - .38;
      G4double lm = pLab - 1.;
      G4double md = lm*lm + .392;   
      fElasticXsc = .7/(lr*lr + .076) + 2./md;
      fTotalXsc   = // .7/(lr*lr + .076) + 
                2.6/md;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double lr  = pLab - .38;
      G4double LE  = .7/(lr*lr + .076);
      G4double sp  = std::sqrt(pLab);
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
      G4double lm  = pLab - 0.8;
      G4double md  = lm*lm + .652;
      fElasticXsc  = LE + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) + 2./md;
      fTotalXsc    = (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) + 7.6/md; // + LE;
    }
  }
  else if( (theParticle == theKPlus || theParticle == theK0L) && neutron )  // Kpn/aKn //////////////////////////////////
  {
    if( pLab < pMin )
    {
      G4double lm = pLab - 0.94;
      G4double md = lm*lm + .392;   
      fElasticXsc = 2./md;
      fTotalXsc   = 4.6/md;
    }
    else if( pLab > pMax )
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc           = cofLogE*ld2 + 2.23;
      fTotalXsc           = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = G4Log(pLab) - minLogP;
      G4double ld2 = ld*ld;
      G4double sp  = std::sqrt(pLab);
      G4double p2  = pLab*pLab;
      G4double p4  = p2*p2;
      G4double lm  = pLab - 0.8;
      G4double md  = lm*lm + .652;
      fElasticXsc  = (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) + 2./md;
      fTotalXsc    = (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) + 7.6/md;
    }
  }
  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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetTotalHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetTotalHadronNucleonXsc ( )

inline

Definition at line 90 of file G4HadronNucleonXsc.hh.

{ return fTotalXsc;     }; 

IsApplicable()

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

virtual

Definition at line 97 of file G4HadronNucleonXsc.cc.

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

Here is the call graph for this function:

IsIsoApplicable()

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

virtual

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

fElasticXsc

G4double G4HadronNucleonXsc::fElasticXsc

private

Definition at line 107 of file G4HadronNucleonXsc.hh.

fInelasticXsc

G4double G4HadronNucleonXsc::fInelasticXsc

private

Definition at line 107 of file G4HadronNucleonXsc.hh.

fLowerLimit

const G4double G4HadronNucleonXsc::fLowerLimit

private

Definition at line 92 of file G4HadronNucleonXsc.hh.

fTotalXsc

G4double G4HadronNucleonXsc::fTotalXsc

private

Definition at line 107 of file G4HadronNucleonXsc.hh.

theA

G4ParticleDefinition* G4HadronNucleonXsc::theA

private

Definition at line 157 of file G4HadronNucleonXsc.hh.

theANeutron

G4ParticleDefinition* G4HadronNucleonXsc::theANeutron

private

Definition at line 133 of file G4HadronNucleonXsc.hh.

theAntiL

G4ParticleDefinition* G4HadronNucleonXsc::theAntiL

private

Definition at line 142 of file G4HadronNucleonXsc.hh.

theAOmega

G4ParticleDefinition* G4HadronNucleonXsc::theAOmega

private

Definition at line 154 of file G4HadronNucleonXsc.hh.

theAProton

G4ParticleDefinition* G4HadronNucleonXsc::theAProton

private

Definition at line 132 of file G4HadronNucleonXsc.hh.

theAS0

G4ParticleDefinition* G4HadronNucleonXsc::theAS0

private

Definition at line 148 of file G4HadronNucleonXsc.hh.

theASMinus

G4ParticleDefinition* G4HadronNucleonXsc::theASMinus

private

Definition at line 146 of file G4HadronNucleonXsc.hh.

theASPlus

G4ParticleDefinition* G4HadronNucleonXsc::theASPlus

private

Definition at line 144 of file G4HadronNucleonXsc.hh.

theAXi0

G4ParticleDefinition* G4HadronNucleonXsc::theAXi0

private

Definition at line 152 of file G4HadronNucleonXsc.hh.

theAXiMinus

G4ParticleDefinition* G4HadronNucleonXsc::theAXiMinus

private

Definition at line 151 of file G4HadronNucleonXsc.hh.

theD

G4ParticleDefinition* G4HadronNucleonXsc::theD

private

Definition at line 155 of file G4HadronNucleonXsc.hh.

theGamma

G4ParticleDefinition* G4HadronNucleonXsc::theGamma

private

Definition at line 129 of file G4HadronNucleonXsc.hh.

theHe3

G4ParticleDefinition* G4HadronNucleonXsc::theHe3

private

Definition at line 158 of file G4HadronNucleonXsc.hh.

theK0L

G4ParticleDefinition* G4HadronNucleonXsc::theK0L

private

Definition at line 140 of file G4HadronNucleonXsc.hh.

theK0S

G4ParticleDefinition* G4HadronNucleonXsc::theK0S

private

Definition at line 139 of file G4HadronNucleonXsc.hh.

theKMinus

G4ParticleDefinition* G4HadronNucleonXsc::theKMinus

private

Definition at line 138 of file G4HadronNucleonXsc.hh.

theKPlus

G4ParticleDefinition* G4HadronNucleonXsc::theKPlus

private

Definition at line 137 of file G4HadronNucleonXsc.hh.

theL

G4ParticleDefinition* G4HadronNucleonXsc::theL

private

Definition at line 141 of file G4HadronNucleonXsc.hh.

theNeutron

G4ParticleDefinition* G4HadronNucleonXsc::theNeutron

private

Definition at line 131 of file G4HadronNucleonXsc.hh.

theOmega

G4ParticleDefinition* G4HadronNucleonXsc::theOmega

private

Definition at line 153 of file G4HadronNucleonXsc.hh.

thePiMinus

G4ParticleDefinition* G4HadronNucleonXsc::thePiMinus

private

Definition at line 135 of file G4HadronNucleonXsc.hh.

thePiPlus

G4ParticleDefinition* G4HadronNucleonXsc::thePiPlus

private

Definition at line 134 of file G4HadronNucleonXsc.hh.

thePiZero

G4ParticleDefinition* G4HadronNucleonXsc::thePiZero

private

Definition at line 136 of file G4HadronNucleonXsc.hh.

theProton

G4ParticleDefinition* G4HadronNucleonXsc::theProton

private

Definition at line 130 of file G4HadronNucleonXsc.hh.

theS0

G4ParticleDefinition* G4HadronNucleonXsc::theS0

private

Definition at line 147 of file G4HadronNucleonXsc.hh.

theSMinus

G4ParticleDefinition* G4HadronNucleonXsc::theSMinus

private

Definition at line 145 of file G4HadronNucleonXsc.hh.

theSPlus

G4ParticleDefinition* G4HadronNucleonXsc::theSPlus

private

Definition at line 143 of file G4HadronNucleonXsc.hh.

theT

G4ParticleDefinition* G4HadronNucleonXsc::theT

private

Definition at line 156 of file G4HadronNucleonXsc.hh.

theXi0

G4ParticleDefinition* G4HadronNucleonXsc::theXi0

private

Definition at line 150 of file G4HadronNucleonXsc.hh.

theXiMinus

G4ParticleDefinition* G4HadronNucleonXsc::theXiMinus

private

Definition at line 149 of file G4HadronNucleonXsc.hh.

---

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

• G4HadronNucleonXsc.hh
• G4HadronNucleonXsc.cc