#include <G4HadronNucleonXsc.hh>

Public Member Functions
	G4HadronNucleonXsc ()

	~G4HadronNucleonXsc ()=default

G4double	HadronNucleonXsc (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	HadronNucleonXscPDG (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	HadronNucleonXscNS (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	KaonNucleonXscNS (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	KaonNucleonXscGG (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	KaonNucleonXscVG (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	HyperonNucleonXscNS (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	SCBMesonNucleonXscNS (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	HadronNucleonXscVU (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	HadronNucleonXscEL (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	CoulombBarrier (const G4ParticleDefinition theParticle, const G4ParticleDefinition nucleon, G4double ekin)

G4double	GetHadronNucleonXscEL (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetHadronNucleonXscPDG (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetHadronNucleonXscNS (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetKaonNucleonXscGG (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetHyperonNucleonXscNS (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetHadronNucleonXscVU (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetCoulombBarrier (const G4DynamicParticle dp, const G4ParticleDefinition p)

G4double	GetTotalHadronNucleonXsc () const

G4double	GetElasticHadronNucleonXsc () const

G4double	GetInelasticHadronNucleonXsc () const

void	CrossSectionDescription (std::ostream &) const

Detailed Description

Definition at line 51 of file G4HadronNucleonXsc.hh.

Constructor & Destructor Documentation

◆ G4HadronNucleonXsc()

G4HadronNucleonXsc::G4HadronNucleonXsc ( )

Definition at line 66 of file G4HadronNucleonXsc.cc.

{
  // basic hadrons
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  thePiPlus   = G4PionPlus::PionPlus();
 
  // basic strange mesons
  theKPlus    = G4KaonPlus::KaonPlus();
  theKMinus   = G4KaonMinus::KaonMinus();
  theK0S      = G4KaonZeroShort::KaonZeroShort();
  theK0L      = G4KaonZeroLong::KaonZeroLong();
 
  g4calc = G4Pow::GetInstance();
}

◆ ~G4HadronNucleonXsc()

G4HadronNucleonXsc::~G4HadronNucleonXsc ( )

default

Member Function Documentation

◆ CoulombBarrier()

G4double G4HadronNucleonXsc::CoulombBarrier	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 1328 of file G4HadronNucleonXsc.cc.

{
  G4double tR = 0.895*CLHEP::fermi;
  G4double pR = 0.5*CLHEP::fermi;
 
  if     ( theParticle == theProton ) pR = 0.895*CLHEP::fermi;
  else if( theParticle == thePiPlus ) pR = 0.663*CLHEP::fermi;
  else if( theParticle == theKPlus )  pR = 0.340*CLHEP::fermi;
 
  G4double pZ = theParticle->GetPDGCharge();
  G4double tZ = nucleon->GetPDGCharge();
 
  G4double pM = theParticle->GetPDGMass(); 
  G4double tM = nucleon->GetPDGMass();
 
  G4double pElab = ekin + 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/(2.*(pR + tR));
 
  //G4cout<<"##CB: ekin = "<<ekin<<"; pElab= " << pElab
  //    <<"; bC = "<<bC<<"; Tcm = "<< totTcm <<G4endl;
 
  G4double ratio = (totTcm > bC) ? 1. - bC/totTcm : 0.0;
 
  // G4cout <<"ratio = "<<ratio<<G4endl;
  return ratio;
}

Referenced by GetCoulombBarrier(), and HadronNucleonXscPDG().

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

◆ GetCoulombBarrier()

G4double G4HadronNucleonXsc::GetCoulombBarrier	(	const G4DynamicParticle *	dp,
		const G4ParticleDefinition *	p )

inline

Definition at line 117 of file G4HadronNucleonXsc.hh.

119 { return CoulombBarrier(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4DynamicParticle::GetDefinition

G4ParticleDefinition * GetDefinition() const

G4DynamicParticle::GetKineticEnergy

G4double GetKineticEnergy() const

G4HadronNucleonXsc::CoulombBarrier

G4double CoulombBarrier(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:1328

◆ GetElasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetElasticHadronNucleonXsc ( ) const

inline

Definition at line 123 of file G4HadronNucleonXsc.hh.

123{ return fElasticXsc; };

◆ GetHadronNucleonXscEL()

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

inline

Definition at line 93 of file G4HadronNucleonXsc.hh.

95 { return HadronNucleonXscEL(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::HadronNucleonXscEL

G4double HadronNucleonXscEL(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:1263

◆ GetHadronNucleonXscNS()

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

inline

Definition at line 101 of file G4HadronNucleonXsc.hh.

103 { return HadronNucleonXscNS(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::HadronNucleonXscNS

G4double HadronNucleonXscNS(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:336

◆ GetHadronNucleonXscPDG()

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

inline

Definition at line 97 of file G4HadronNucleonXsc.hh.

99 { return HadronNucleonXscPDG(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::HadronNucleonXscPDG

G4double HadronNucleonXscPDG(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:144

◆ GetHadronNucleonXscVU()

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

inline

Definition at line 113 of file G4HadronNucleonXsc.hh.

115 { return HadronNucleonXscVU(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::HadronNucleonXscVU

G4double HadronNucleonXscVU(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:1118

◆ GetHyperonNucleonXscNS()

G4double G4HadronNucleonXsc::GetHyperonNucleonXscNS	(	const G4DynamicParticle *	dp,
		const G4ParticleDefinition *	p )

inline

Definition at line 109 of file G4HadronNucleonXsc.hh.

111 { return HyperonNucleonXscNS(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::HyperonNucleonXscNS

G4double HyperonNucleonXscNS(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:993

◆ GetInelasticHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetInelasticHadronNucleonXsc ( ) const

inline

Definition at line 124 of file G4HadronNucleonXsc.hh.

124{ return fInelasticXsc; };

◆ GetKaonNucleonXscGG()

G4double G4HadronNucleonXsc::GetKaonNucleonXscGG	(	const G4DynamicParticle *	dp,
		const G4ParticleDefinition *	p )

inline

Definition at line 105 of file G4HadronNucleonXsc.hh.

107 { return KaonNucleonXscGG(dp->GetDefinition(), p, dp->GetKineticEnergy()); }

G4HadronNucleonXsc::KaonNucleonXscGG

G4double KaonNucleonXscGG(const G4ParticleDefinition *theParticle, const G4ParticleDefinition *nucleon, G4double ekin)

Definition G4HadronNucleonXsc.cc:778

◆ GetTotalHadronNucleonXsc()

G4double G4HadronNucleonXsc::GetTotalHadronNucleonXsc ( ) const

inline

Definition at line 122 of file G4HadronNucleonXsc.hh.

122{ return fTotalXsc; };

◆ HadronNucleonXsc()

G4double G4HadronNucleonXsc::HadronNucleonXsc	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 92 of file G4HadronNucleonXsc.cc.

{
  G4double xsc(0.);
  G4int pdg = std::abs(theParticle->GetPDGEncoding());
  
  // p, n, pi+-, pbar, nbar
  if ( pdg == 2212 || pdg == 2112 || pdg == 211 ) { 
    xsc = HadronNucleonXscNS(theParticle, nucleon, ekin);
  }
  else if ( pdg == 22 ) // gamma
  {
    xsc = HadronNucleonXscPDG(theParticle, nucleon, ekin);
  }
  else if ( pdg == 321 || pdg == 310 || pdg == 130 ) // K+-, K0L, K0S
  {
    xsc = KaonNucleonXscNS(theParticle, nucleon, ekin);
  }
  else if (pdg > 3000) 
  {
    if (pdg == 3122 || pdg == 3222 || pdg == 3112 || pdg == 3212 || pdg == 3322 || pdg == 3312 || pdg == 3324 ||
    pdg == 4122 || pdg == 4332 || pdg == 4122 || pdg == 4212 || pdg == 4222 || pdg == 4112 || pdg == 4232 || pdg == 4132 ||
        pdg == 5122 || pdg == 5332 || pdg == 5122 || pdg == 5112 || pdg == 5222 || pdg == 5212 || pdg == 5132 || pdg == 5232
       ) // heavy s-,c-,b-hyperons
    {
      xsc = HyperonNucleonXscNS(theParticle, nucleon, ekin);
    }
    else 
    {
      xsc = HadronNucleonXscPDG(theParticle, nucleon, ekin);
    }
  } else if (pdg > 220) {
    if (pdg == 511 || pdg == 421 || pdg == 531 || pdg == 541 || pdg == 431 || pdg == 411 || pdg == 521 ||
    pdg == 221 || pdg == 331 || pdg == 441 || pdg == 443 || pdg == 543) // s-,c-,b-mesons
    {
      xsc = SCBMesonNucleonXscNS(theParticle, nucleon, ekin);
    }
    else 
    {
      xsc = HadronNucleonXscPDG(theParticle, nucleon, ekin);
    }
  } else {
    xsc = HadronNucleonXscPDG(theParticle, nucleon, ekin);
  }
  return xsc;
}

◆ HadronNucleonXscEL()

G4double G4HadronNucleonXsc::HadronNucleonXscEL	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 1263 of file G4HadronNucleonXsc.cc.

{
  G4int pdg = theParticle->GetPDGEncoding();
  G4double xsection(0.);
  static const G4double targ_mass = 
    0.5*(CLHEP::proton_mass_c2 + CLHEP::neutron_mass_c2);
  G4double sMand = 
    CalcMandelstamS(ekin, theParticle->GetPDGMass(), targ_mass)*invGeV2;
 
  G4double x1 = G4Exp(G4Log(sMand)*0.0808);
  G4double x2 = G4Exp(G4Log(-sMand)*0.4525);
 
  if(pdg == 22) 
  {
    xsection = 0.0677*x1 + 0.129*x2;
  } 
  else if(theParticle == theNeutron)  
  {
    xsection = 21.70*x1 + 56.08*x2;
  } 
  else if(theParticle == theProton) 
  {
    xsection = 21.70*x1 + 56.08*x2;
  } 
  // pbar
  else if(pdg == -2212) 
  {
    xsection = 21.70*x1 + 98.39*x2;
  } 
  else if(theParticle == thePiPlus) 
  {
    xsection = 13.63*x1 + 27.56*x2;
  } 
  // pi-
  else if(pdg == -211) 
  {
    xsection = 13.63*x1 + 36.02*x2;
  } 
  else if(theParticle == theKPlus) 
  {
    xsection = 11.82*x1 + 8.15*x2;
  } 
  else if(theParticle == theKMinus) 
  {
    xsection = 11.82*x1 + 26.36*x2;
  }
  else if(theParticle == theK0S || theParticle == theK0L) 
  {
    xsection = 11.82*x1 + 17.25*x2;
  }
  else  
  {
    xsection = 21.70*x1 + 56.08*x2;
  }
  fTotalXsc     = xsection*CLHEP::millibarn;
  fInelasticXsc = 0.83*fTotalXsc;
  fElasticXsc   = fTotalXsc - fInelasticXsc;
  return fTotalXsc;
}

Referenced by GetHadronNucleonXscEL().

◆ HadronNucleonXscNS()

G4double G4HadronNucleonXsc::HadronNucleonXscNS	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 336 of file G4HadronNucleonXsc.cc.

{  
  const G4double ekin = std::max(ekin0, ekinmin);
  G4int pdg = theParticle->GetPDGEncoding();
  /*  
  G4cout<< "HadronNucleonXscNS: Ekin(GeV)= " << ekin/GeV << "  "
        << theParticle->GetParticleName() << " + " 
        << nucleon->GetParticleName() << G4endl;
  */
  if(pdg == -2212 || pdg == -2112) {
    return HadronNucleonXscPDG(theParticle, nucleon, ekin);
  }
 
  G4double pM   = theParticle->GetPDGMass();
  G4double tM   = nucleon->GetPDGMass();
  G4double pE   = ekin + pM; 
  G4double pLab = std::sqrt(ekin*(ekin + 2*pM));
 
  G4double sMand = CalcMandelstamS(ekin, pM, tM)*invGeV2;
 
  pLab *= invGeV;
  pE   *= invGeV;
 
  if(pLab >= 10.) {
    fTotalXsc = HadronNucleonXscPDG(theParticle, nucleon, ekin)/CLHEP::millibarn;
  } else { fTotalXsc = 0.0; }
  fElasticXsc = 0.0;
  //G4cout << "Stot(mb)= " << fTotalXsc << " pLab= " << pLab 
  //     << " Smand= " << sMand <<G4endl;
  G4double logP = G4Log(pLab);
 
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);
 
  if( theParticle == theNeutron) 
  {
    if( pLab >= 373.)
    {
      fElasticXsc = 6.5 + 0.308*G4Exp(G4Log(G4Log(sMand/400.)*1.65)) 
    + 9.19*G4Exp(-G4Log(sMand)*0.458);
    }
    else if( pLab >= 100.)
    {
      fElasticXsc = 5.53 + 0.308*G4Exp(G4Log(G4Log(sMand/28.9))*1.1) 
    + 9.19*G4Exp(-G4Log(sMand)*0.458);
    }
    else if( pLab >= 10.)
    {
      fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
    }
    else  // pLab < 10 GeV/c
    {
      if( neutron )      // nn to be pp
      {
        G4double x = G4Log(pLab/0.73);
        if( pLab < 0.4 )
        {
          fTotalXsc = 23 + 50*std::sqrt(g4calc->powN(-x, 7));
          fElasticXsc = fTotalXsc;
        }
        else if( pLab < 0.73 )
        {
          fTotalXsc = 23 + 50*std::sqrt(g4calc->powN(-x, 7));
          fElasticXsc = fTotalXsc; 
        }
        else if( pLab < 1.05  )
        {
          fTotalXsc = 23 + 40*x*x;
          fElasticXsc = 23 + 20*x*x;
        }
        else    // 1.05 - 10 GeV/c
        {
          fTotalXsc = 39.0+75*(pLab - 1.2)/(g4calc->powN(pLab,3) + 0.15);
          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
      }
      if( proton )   // pn to be np
      {
        if( pLab < 0.02 )
        {
          fTotalXsc = 4100+30*G4Exp(G4Log(G4Log(1.3/pLab))*3.6); 
      fElasticXsc = fTotalXsc;
        }      
        else if( pLab < 0.8 )
        {
          fTotalXsc = 33+30*g4calc->powN(G4Log(pLab/1.3),4);
      fElasticXsc = fTotalXsc;
        }      
        else if( pLab < 1.4 )
        {
          fTotalXsc = 33+30*g4calc->powN(G4Log(pLab/0.95),2);
          G4double x = G4Log(0.511/pLab);
          fElasticXsc =  6 + 52/( x*x + 1.6 );
        }
        else    // 1.4 < pLab < 10.  )
        {
          fTotalXsc = 33.3 + 20.8*(pLab*pLab - 1.35)
        /(std::sqrt(g4calc->powN(pLab,5)) + 0.95);
          fElasticXsc =  6 + 20/( (logP-0.182)*(logP-0.182) + 1.0 );
        }
      }
    }
  }
  ////// proton //////////////////////////////////////////////
  else if(theParticle == theProton) 
  {
    if( pLab >= 373.) // pdg due to TOTEM data
    {
      fElasticXsc = 6.5 + 0.308*G4Exp(G4Log(G4Log(sMand/400.))*1.65) 
    + 9.19*G4Exp(-G4Log(sMand)*0.458);     
    }
    else if( pLab >= 100.)
    {
      fElasticXsc = 5.53 + 0.308*G4Exp(G4Log(G4Log(sMand/28.9))*1.1) 
    + 9.19*G4Exp(-G4Log(sMand)*0.458);     
    }
    else if( pLab >= 10.)
    {
      fElasticXsc = 6. + 20./( (logP-0.182)*(logP-0.182) + 1.0 );
    }
    else
    {
      // pp
      if( proton )
      {
        if( pLab < 0.73 )
        {
          fTotalXsc = 23 + 50*std::sqrt(g4calc->powN(G4Log(0.73/pLab),7));
          fElasticXsc = fTotalXsc;
        }
        else if( pLab < 1.05  )
        {
      G4double x = G4Log(pLab/0.73);
          fTotalXsc = 23 + 40*x*x;
          fElasticXsc = 23 + 20*x*x;
        }
        else    // 1.05 - 10 GeV/c
        {
          fTotalXsc = 39.0+75*(pLab - 1.2)/(g4calc->powN(pLab,3) + 0.15);
          fElasticXsc = 6. + 20./( (logP-0.182)*(logP-0.182) + 1.0 );
        }
      }
      else if( neutron )     // pn to be np
      {
        if( pLab < 0.02 )
        {
          fTotalXsc = 4100+30*G4Exp(G4Log(G4Log(1.3/pLab))*3.6); 
      fElasticXsc = fTotalXsc;
        }      
        else if( pLab < 0.8 )
        {
          fTotalXsc = 33+30*g4calc->powN(G4Log(pLab/1.3),4);
      fElasticXsc = fTotalXsc;
        }      
        else if( pLab < 1.4 )
        {
          G4double x1 = G4Log(pLab/0.95);
          G4double x2 = G4Log(0.511/pLab);
          fTotalXsc = 33+30*x1*x1;
          fElasticXsc =  6 + 52/(x2*x2 + 1.6);
        }
        else    // 1.4 < pLab < 10.  )
        {
          fTotalXsc = 33.3 + 20.8*(pLab*pLab - 1.35)
        /(std::sqrt(g4calc->powN(pLab,5)) + 0.95);
          fElasticXsc = 6. + 20./( (logP-0.182)*(logP-0.182) + 1.0 );
        }
      }
    }    
  }
  // pi+ p; pi- n 
  else if((pdg == 211 && proton) || (pdg == -211 && neutron))  
  {
    if( pLab < 0.28 )
    {
      fTotalXsc = 10./((logP + 1.273)*(logP + 1.273) + 0.05);
      fElasticXsc = fTotalXsc;
    }
    else if( pLab < 0.68 )
    {
      fTotalXsc = 14./((logP + 1.273)*(logP + 1.273) + 0.07);
      fElasticXsc = fTotalXsc;
    }
    else if( pLab < 0.85 )
    {
      G4double x = G4Log(pLab/0.77);
      fTotalXsc = 88.*x*x + 14.9;
      fElasticXsc = fTotalXsc*G4Exp(-3.*(pLab - 0.68));  
    }
    else if( pLab < 1.15 )
    {
      G4double x = G4Log(pLab/0.77);
      fTotalXsc = 88.*x*x + 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);
      fTotalXsc       = 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);
      fTotalXsc     = 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);
      fTotalXsc       = Ex1 + Ex2 + 27.5;
      fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
    }
    else if( pLab < 10. ) // my
    {
      fTotalXsc = 10.6 + 2.*G4Log(pE) + 25*G4Exp(-G4Log(pE)*0.43 ); 
      fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
    }
    else //  pLab > 10 // my
    {
      fElasticXsc = 3.0 + 6.20/( (logP - 0.336)*(logP - 0.336) + 0.8);    
    }
  }
  // pi+ n; pi- p
  else if((pdg == 211 && neutron) || (pdg == -211 && proton))  
  {
    if( pLab < 0.28 ) 
    {
      fTotalXsc = 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
    {
      fTotalXsc = 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 )
    {
      G4double y = G4Log(pLab/0.48);
      fTotalXsc = 26 + 110*y*y;
      fElasticXsc = 0.37*fTotalXsc;
    }
    else if( pLab < 0.65 )
    {
      G4double x = G4Log(pLab/0.48);
      fTotalXsc = 26. + 110.*x*x;
      fElasticXsc = 0.95/((pLab - 0.72)*(pLab - 0.72) + 0.049);
    }
    else if( pLab < 0.72 )
    {
      fTotalXsc = 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 )
    {
      fTotalXsc = 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 )
    {
      fTotalXsc = 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 )
    {
      fTotalXsc = 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 )
    {
      fTotalXsc = 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 < 10.) // < 3.0) // ns original
    {
      fTotalXsc = 36.1 + 0.079-4.313*logP+
    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 
    {
      fElasticXsc = 3. + 13./pLab;
    }
  }
  else if( (theParticle == theKMinus) && proton )   // K-p
  {
    if( pLab < pMin)
    {
      G4double psp = pLab*std::sqrt(pLab);
      fElasticXsc  = 5.2/psp;
      fTotalXsc    = 14./psp;
    }
    else if( pLab > pMax )
    {
      G4double ld  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
      G4double ld  = logP - 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 - 1.01;
      G4double hd  = lh*lh + .011;
 
      G4double lm  = pLab - .39;
      G4double md  = lm*lm + .000356;
      G4double lh1  = pLab - 0.78;
      G4double hd1  = lh1*lh1 + .00166;
      G4double lh2  = pLab - 1.63;
      G4double hd2  = lh2*lh2 + .007;
 
      // small peaks were added but commented out now
      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; 
 
      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) && neutron )  // K-n
  {
    if( pLab > pMax )
    {
      G4double ld  = logP - 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 sqrLogPlab = logP*logP;
 
      fElasticXsc = 5.0 +  8.1*G4Exp(-logP*1.8) + 0.16*sqrLogPlab 
    - 1.3*logP + .15/hd;
      fTotalXsc = 25.2 +  0.38*sqrLogPlab - 2.9*logP + 0.30/hd;
    }
  }
  else if( (theParticle == theKPlus) && proton  )  // K+p
  {
    // VI: modified low-energy part
    if( pLab < 0.631 )
    {
      fElasticXsc = fTotalXsc = 12.03;
    }
    else if( pLab > pMax )
    {
      G4double ld  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = logP - 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;
      // VI: tuned elastic
      fElasticXsc  = LE + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) 
    + 2./((pLab - 0.8)*(pLab - 0.8) + 0.652);
      fTotalXsc    = LE + (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) 
    + 2.6/((pLab - 1.)*(pLab - 1.) + 0.392);
    }
  }
  else if(  (theParticle == theKPlus) && neutron) // K+n  
  {
    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  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = logP - 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;
    }
  }
  fTotalXsc   *= CLHEP::millibarn; 
  fElasticXsc *= CLHEP::millibarn;
  fElasticXsc  = std::min(fElasticXsc, fTotalXsc);
 
  if( proton && theParticle->GetPDGCharge() > 0. && ekin < ekinmaxQB )
  {
    G4double cB = G4NuclearRadii::CoulombFactor(theParticle, nucleon, ekin);
    fTotalXsc   *= cB;
    fElasticXsc *= cB; 
  }
  fInelasticXsc = std::max(fTotalXsc - fElasticXsc,0.0);
  /*
  G4cout<< "HNXsc: Ekin(GeV)= " << ekin/GeV << "; tot(mb)= " << fTotalXsc/millibarn
        <<"; el(mb)= " <<fElasticXsc/millibarn
        <<"; inel(mb)= " <<fInelasticXsc/millibarn<<"  "
        << theParticle->GetParticleName() << " + " 
        << nucleon->GetParticleName() << G4endl;
  */
  return fTotalXsc;
}

Referenced by GetHadronNucleonXscNS(), HadronNucleonXsc(), HyperonNucleonXscNS(), KaonNucleonXscNS(), and SCBMesonNucleonXscNS().

◆ HadronNucleonXscPDG()

G4double G4HadronNucleonXsc::HadronNucleonXscPDG	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 144 of file G4HadronNucleonXsc.cc.

{
  static const G4double M    = 2.1206; // in Gev
  static const G4double eta1 = 0.4473;
  static const G4double eta2 = 0.5486;
  static const G4double H    = 0.272;
 
  G4int pdg = theParticle->GetPDGEncoding();
 
  G4double mass1 = (pdg == 22) ? 770. : theParticle->GetPDGMass();
  G4double mass2 = nucleon->GetPDGMass();
 
  G4double sMand = CalcMandelstamS(ekin, mass1, mass2)*invGeV2;
  G4double x = (mass1 + mass2)*invGeV + M;
  G4double blog = G4Log(sMand/(x*x));
 
  G4double P(0.0), R1(0.0), R2(0.0), del(1.0);
 
  G4bool proton  = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);
  
  if(theParticle == theNeutron)
  {
    if ( proton )
    {
      P  = 34.71;
      R1 = 12.52;
      R2 = -6.66;
    }
    else
    {
      P  = 34.41;
      R1 = 13.07;
      R2 = -7.394;
    }
  } 
  else if(theParticle == theProton) 
  {
    if ( neutron )
    {
      P  = 34.71;
      R1 = 12.52;
      R2 = -6.66;
    }
    else
    {
      P  = 34.41;
      R1 = 13.07;
      R2 = -7.394;
    }
  } 
  // pbar
  else if(pdg == -2212) 
  {
    if ( neutron )
    {
      P  = 34.71;
      R1 = 12.52;
      R2 = 6.66;
    }
    else
    {
      P  = 34.41;
      R1 = 13.07;
      R2 = 7.394;
    }
  } 
  // nbar
  else if(pdg == -2112) 
  {
    if ( proton )
    {
      P  = 34.71;
      R1 = 12.52;
      R2 = 6.66;
    }
    else
    {
      P  = 34.41;
      R1 = 13.07;
      R2 = 7.394;
    }
  } 
  // pi+
  else if(pdg == 211) 
  {
    P  = 18.75;
    R1 = 9.56;
    R2 = -1.767;
  } 
  // pi-
  else if(pdg == -211) 
  {
    P  = 18.75;
    R1 = 9.56;
    R2 = 1.767;
  } 
  else if(theParticle == theKPlus) 
  {
    if ( proton )
    {
      P  = 16.36;
      R1 = 4.29;
      R2 = -3.408;
    }
    else
    {
      P  = 16.31;
      R1 = 3.7;
      R2 = -1.826;
    }
  } 
  else if(theParticle == theKMinus) 
  {
    if ( proton )
    {
      P  = 16.36;
      R1 = 4.29;
      R2 = 3.408;
    }
    else
    {
      P  = 16.31;
      R1 = 3.7;
      R2 = 1.826;
    }
  }
  else if(theParticle == theK0S || theParticle == theK0L) 
  {
    P  = 16.36;
    R1 = 2.5;
    R2 = 0.;
  }
  // sigma-
  else if(pdg == 3112) 
  {
    P  = 34.7;
    R1 = -46.;
    R2 = 48.;
  } 
  // gamma
  else if(pdg == 22) // modify later on
  {
    del= 0.003063;
    P  = 34.71*del;
    R1 = (neutron) ? 0.0231 : 0.0139;
    R2 = 0.;
  } 
  else  // as proton ??? 
  {
    if ( neutron )
    {
      P  = 34.71;
      R1 = 12.52;
      R2 = -6.66;
    }
    else
    {
      P  = 34.41;
      R1 = 13.07;
      R2 = -7.394;
    }
  } 
  fTotalXsc = CLHEP::millibarn* 
    (del*(H*blog*blog + P) + R1*G4Exp(-eta1*blog) + R2*G4Exp(-eta2*blog));
  fInelasticXsc = 0.75*fTotalXsc;
  fElasticXsc   = fTotalXsc - fInelasticXsc;
 
  if( proton && theParticle->GetPDGCharge() > 0. && ekin < ekinmaxQB )
  {
    G4double cB = CoulombBarrier(theParticle, nucleon, ekin);
    fTotalXsc   *= cB;
    fElasticXsc *= cB; 
    fInelasticXsc *= cB; 
  }
  /*
  G4cout << "## HadronNucleonXscPDG: ekin(MeV)= " << ekin 
     << " tot= " << fTotalXsc/CLHEP::millibarn 
     << " inel= " << fInelasticXsc/CLHEP::millibarn
     << " el= " << fElasticXsc/CLHEP::millibarn
         << G4endl;
  */
  return fTotalXsc;
}

Referenced by GetHadronNucleonXscPDG(), HadronNucleonXsc(), and HadronNucleonXscNS().

◆ HadronNucleonXscVU()

G4double G4HadronNucleonXsc::HadronNucleonXscVU	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 1118 of file G4HadronNucleonXsc.cc.

{
  G4int PDGcode = theParticle->GetPDGEncoding();
  G4int absPDGcode = std::abs(PDGcode);
  G4double mass = theParticle->GetPDGMass();
  G4double Plab = std::sqrt(ekin*(ekin + 2.*mass))*invGeV;
 
  G4double logPlab = G4Log( Plab );
  G4double sqrLogPlab = logPlab * logPlab;
 
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);
   
  if( absPDGcode > 1000)  //------Projectile is baryon -
  {
    if(proton)
    {
      fTotalXsc   = 48.0 + 0.522*sqrLogPlab - 4.51*logPlab;
      fElasticXsc = 11.9 + 26.9*G4Exp(-logPlab*1.21) + 0.169*sqrLogPlab - 1.85*logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 47.3  + 0.513*sqrLogPlab - 4.27*logPlab;
      fElasticXsc = 11.9 + 26.9*G4Exp(-logPlab*1.21) + 0.169*sqrLogPlab - 1.85*logPlab;
    }
  }
  else if( PDGcode ==  211)  //------Projectile is PionPlus ----
  {
    if(proton)
    {
      fTotalXsc  = 16.4 + 19.3 *G4Exp(-logPlab*0.42) + 0.19 *sqrLogPlab - 0.0 *logPlab;
      fElasticXsc =  0.0 + 11.4*G4Exp(-logPlab*0.40) + 0.079*sqrLogPlab - 0.0 *logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   =  33.0 + 14.0 *G4Exp(-logPlab*1.36) + 0.456*sqrLogPlab - 4.03*logPlab;
      fElasticXsc = 1.76 + 11.2*G4Exp(-logPlab*0.64) + 0.043*sqrLogPlab - 0.0 *logPlab;
    }
  }
  else if( PDGcode == -211)  //------Projectile is PionMinus ----
  {
    if(proton)
    {
      fTotalXsc   = 33.0 + 14.0 *G4Exp(-logPlab*1.36) + 0.456*sqrLogPlab - 4.03*logPlab;
      fElasticXsc = 1.76 + 11.2*G4Exp(-logPlab*0.64) + 0.043*sqrLogPlab - 0.0 *logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 16.4 + 19.3 *G4Exp(-logPlab*0.42) + 0.19 *sqrLogPlab - 0.0 *logPlab;
      fElasticXsc =  0.0 + 11.4*G4Exp(-logPlab*0.40) + 0.079*sqrLogPlab - 0.0 *logPlab;
    }
  }
  else if( PDGcode ==  111)  //------Projectile is PionZero  --
  {
    if(proton)
    {
      fTotalXsc = (16.4 + 19.3 *G4Exp(-logPlab*0.42) + 0.19 *sqrLogPlab - 0.0 *logPlab +   //Pi+
           33.0 + 14.0 *G4Exp(-logPlab*1.36) + 0.456*sqrLogPlab - 4.03*logPlab)/2; //Pi-
 
      fElasticXsc = (0.0 + 11.4*G4Exp(-logPlab*0.40) + 0.079*sqrLogPlab - 0.0 *logPlab +   //Pi+
             1.76 + 11.2*G4Exp(-logPlab*0.64) + 0.043*sqrLogPlab - 0.0 *logPlab)/2;//Pi-
 
    }
    if(neutron)
    {    
      fTotalXsc = (33.0 + 14.0 *G4Exp(-logPlab*1.36) + 0.456*sqrLogPlab - 4.03*logPlab +   //Pi+
           16.4 + 19.3 *G4Exp(-logPlab*0.42) + 0.19 *sqrLogPlab - 0.0 *logPlab)/2; //Pi-
      fElasticXsc = (1.76 + 11.2*G4Exp(-logPlab*0.64) + 0.043*sqrLogPlab - 0.0 *logPlab +  //Pi+
             0.0  + 11.4*G4Exp(-logPlab*0.40) + 0.079*sqrLogPlab - 0.0 *logPlab)/2;//Pi-
    }
  }
  else if( PDGcode == 321 )    //------Projectile is KaonPlus --
  {
    if(proton)
    {
      fTotalXsc   = 18.1 +  0.26 *sqrLogPlab - 1.0 *logPlab;
      fElasticXsc =  5.0 +  8.1*G4Exp(-logPlab*1.8 ) + 0.16 *sqrLogPlab - 1.3 *logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 18.7  + 0.21 *sqrLogPlab - 0.89*logPlab;
      fElasticXsc =  7.3  + 0.29 *sqrLogPlab - 2.4 *logPlab;
    }
  }
  else if( PDGcode ==-321 )  //------Projectile is KaonMinus ----
  {
    if(proton)
    {
      fTotalXsc   = 32.1 + 0.66*sqrLogPlab - 5.6*logPlab;
      fElasticXsc =  7.3 + 0.29*sqrLogPlab - 2.4*logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 25.2 + 0.38*sqrLogPlab - 2.9*logPlab;
      fElasticXsc =  5.0 +  8.1*G4Exp(-logPlab*1.8 ) + 0.16*sqrLogPlab - 1.3*logPlab;
    }
  }
  else if( PDGcode == 311 )  //------Projectile is KaonZero -----
  {
    if(proton)
    {
      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*G4Exp(-logPlab*1.8 ) + 0.16 *sqrLogPlab - 1.3 *logPlab + //K+
                         7.3 + 0.29 *sqrLogPlab - 2.4 *logPlab)/2; //K-
    }
    if(neutron)
    {    
      fTotalXsc   = ( 18.7 + 0.21 *sqrLogPlab - 0.89*logPlab +   //K+
                      25.2 + 0.38 *sqrLogPlab - 2.9 *logPlab)/2; //K-
      fElasticXsc = (  7.3 + 0.29 *sqrLogPlab - 2.4 *logPlab +   //K+
                       5.0 + 8.1*G4Exp(-logPlab*1.8 ) + 0.16 *sqrLogPlab - 1.3 *logPlab)/2; //K-
    }
  }
  else  //------Projectile is undefined, Nucleon assumed
  {
    if(proton)
    {
      fTotalXsc   = 48.0 + 0.522*sqrLogPlab - 4.51*logPlab;
      fElasticXsc = 11.9 + 26.9*G4Exp(-logPlab*1.21) + 0.169*sqrLogPlab - 1.85*logPlab;
    }
    if(neutron)
    {    
      fTotalXsc   = 47.3 + 0.513*sqrLogPlab - 4.27*logPlab;
      fElasticXsc = 11.9 + 26.9*G4Exp(-logPlab*1.21) + 0.169*sqrLogPlab - 1.85*logPlab;
    }
  }
 
  fTotalXsc   *= CLHEP::millibarn;
  fElasticXsc *= CLHEP::millibarn;
  fElasticXsc = std::min(fElasticXsc, fTotalXsc);
  fInelasticXsc = fTotalXsc - fElasticXsc;
 
  return fTotalXsc;    
}

Referenced by GetHadronNucleonXscVU().

◆ HyperonNucleonXscNS()

G4double G4HadronNucleonXsc::HyperonNucleonXscNS	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 993 of file G4HadronNucleonXsc.cc.

{
  G4double coeff = 1.0;
  G4int pdg = std::abs(theParticle->GetPDGEncoding());
    
  // lambda, sigma+-0 and anti-hyperons
  if( pdg == 3122 || pdg == 3112 || pdg == 3212 || pdg == 3222 )
  {
    coeff = 0.88;
  } 
  // Xi hyperons and anti-hyperons
  else if( pdg == 3312 || pdg == 3322 )
  {
    coeff = 0.76;
  }
  // omega, anti_omega
  else if( pdg == 3334 )
  {
    coeff = 0.64;
  }
  // lambdaC, sigmaC+-0 and anti-hyperonsC
  else if( pdg == 4122 || pdg == 4112 || pdg == 4212 || pdg == 4222 ) 
  {
    coeff = 0.784378;
  }
  // omegaC0, anti_omegaC0
  else if( pdg == 4332 )
  {
    coeff = 0.544378;
  }
  // XiC+0 and anti-hyperonC
  else if( pdg == 4132 || pdg == 4232 )
  {
    coeff = 0.664378;
  }
  // lambdaB, sigmaB+-0 and anti-hyperonsB
  else if( pdg == 5122 || pdg == 5112 || pdg == 5212 || pdg == 5222 )
  {
    coeff = 0.740659;
  }
  // omegaB0, anti_omegaB0
  else if( pdg == 5332 )
  {
    coeff = 0.500659;
  }
  // XiB+0 and anti-hyperonB
  else if( pdg == 5132 || pdg == 5232 )
  {
    coeff = 0.620659;
  } 
  fTotalXsc = coeff*HadronNucleonXscNS( theProton, nucleon, ekin);
  fInelasticXsc *= coeff;
  fElasticXsc *= coeff;
  
  return fTotalXsc;
}

Referenced by GetHyperonNucleonXscNS(), and HadronNucleonXsc().

◆ KaonNucleonXscGG()

G4double G4HadronNucleonXsc::KaonNucleonXscGG	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 778 of file G4HadronNucleonXsc.cc.

{
  fTotalXsc = fElasticXsc = fInelasticXsc = 0.0;
  if(theParticle == theKMinus || theParticle == theKPlus) {
    KaonNucleonXscVG(theParticle, nucleon, ekin);
 
  } else if(theParticle == theK0S || theParticle == theK0L) {
    G4double stot  = KaonNucleonXscVG(theKMinus, nucleon, ekin);
    G4double sel   = fElasticXsc;
    G4double sinel = fInelasticXsc;
    stot  += KaonNucleonXscVG(theKPlus, nucleon, ekin);
    sel   += fElasticXsc;
    sinel += fInelasticXsc;
    fTotalXsc = stot*0.5; 
    fElasticXsc = sel*0.5; 
    fInelasticXsc = sinel*0.5;
  }
  return fTotalXsc;
}

Referenced by GetKaonNucleonXscGG().

◆ KaonNucleonXscNS()

G4double G4HadronNucleonXsc::KaonNucleonXscNS	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 804 of file G4HadronNucleonXsc.cc.

{
  fTotalXsc = fElasticXsc = fInelasticXsc = 0.0;
  if(theParticle == theKMinus || theParticle == theKPlus) {
    HadronNucleonXscNS(theParticle, nucleon, ekin);
 
  } else if(theParticle == theK0S || theParticle == theK0L) {
    G4double fact = 0.5;
    G4double stot = 0.0;
    G4double sel  = 0.0;
    G4double sinel= 0.0;
    if(ekin > ekinmaxQB) {
      stot  = HadronNucleonXscNS(theKMinus, nucleon, ekin);
      sel   = fElasticXsc;
      sinel = fInelasticXsc;
      stot  += HadronNucleonXscNS(theKPlus, nucleon, ekin);
      sel   += fElasticXsc;
      sinel += fInelasticXsc;
    } else {
      fact *= std::sqrt(ekinmaxQB/std::max(ekin, ekinmin));
      stot  = HadronNucleonXscNS(theKMinus, nucleon, ekinmaxQB);
      sel   = fElasticXsc;
      sinel = fInelasticXsc;
      stot  += HadronNucleonXscNS(theKPlus, nucleon, ekinmaxQB);
      sel   += fElasticXsc;
      sinel += fInelasticXsc;
    }
    fTotalXsc = stot*fact; 
    fElasticXsc = sel*fact; 
    fInelasticXsc = sinel*fact;
  }
  return fTotalXsc;
}

Referenced by HadronNucleonXsc().

◆ KaonNucleonXscVG()

G4double G4HadronNucleonXsc::KaonNucleonXscVG	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 844 of file G4HadronNucleonXsc.cc.

{
  G4double pM   = theParticle->GetPDGMass();
  G4double pLab = std::sqrt(ekin*(ekin + 2*pM));
 
  pLab *= invGeV;
  G4double logP = G4Log(pLab);
 
  fTotalXsc = 0.0;
 
  G4bool proton = (nucleon == theProton);
  G4bool neutron = (nucleon == theNeutron);
 
  if( (theParticle == theKMinus) && proton )   // K-p
  {
    if( pLab < pMin)
    {
      G4double psp = pLab*std::sqrt(pLab);
      fElasticXsc  = 5.2/psp;
      fTotalXsc    = 14./psp;
    }
    else if( pLab > pMax )
    {
      G4double ld  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = 1.1*cofLogT*ld2 + 19.7;
    }
    else
    {
      G4double ld  = logP - 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 - 1.01;
      G4double hd  = lh*lh + .011;
      fElasticXsc  = 5.2/psp + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .075/p4) + .15/hd;
      fTotalXsc    = 14./psp + (1.1*cofLogT*ld2 + 19.5)/(1. - .21/sp + .52/p4) + .60/hd;
    }
  }
  else if( (theParticle == theKMinus) && neutron )  // K-n
  {
    if( pLab > pMax )
    {
      G4double ld  = logP - 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;    // vg version
      G4double sqrLogPlab = logP*logP;
 
      fElasticXsc = 5.0 +  8.1*G4Exp(-logP*1.8) + 0.16*sqrLogPlab 
    - 1.3*logP + .15/hd;
      fTotalXsc = 25.2 +  0.38*sqrLogPlab - 2.9*logP + 0.60/hd;  // vg version
    }
  }
  else if( (theParticle == theKPlus) && proton  )  // K+p
  {
    // VI: modified low-energy part
    if( pLab < 0.631 )
    {
      fElasticXsc = fTotalXsc = 12.03;
    }
    else if( pLab > pMax )
    {
      G4double ld  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = logP - 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;
      // VI: tuned elastic
      fElasticXsc  = LE + (cofLogE*ld2 + 2.23)/(1. - .7/sp + .1/p4) 
    + 2./((pLab - 0.8)*(pLab - 0.8) + 0.652);
      fTotalXsc    = LE + (cofLogT*ld2 + 19.5)/(1. + .46/sp + 1.6/p4) 
    + 2.6/((pLab - 1.)*(pLab - 1.) + 0.392);
    }
  }
  else if(  (theParticle == theKPlus) && neutron) // K+n  
  {
    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  = logP - minLogP;
      G4double ld2 = ld*ld;
      fElasticXsc  = cofLogE*ld2 + 2.23;
      fTotalXsc    = cofLogT*ld2 + 19.2;
    }
    else
    {
      G4double ld  = logP - 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;
    }
  }
 
  fTotalXsc   *= CLHEP::millibarn; 
  fElasticXsc *= CLHEP::millibarn;
 
  if( proton && theParticle->GetPDGCharge() > 0. )
  {
    G4double cB = G4NuclearRadii::CoulombFactor(theParticle, nucleon, ekin);
    fTotalXsc   *= cB;
    fElasticXsc *= cB; 
  }
  fElasticXsc = std::min(fElasticXsc, fTotalXsc);
  fInelasticXsc = std::max(fTotalXsc - fElasticXsc,0.0);
  /*
  G4cout << "HNXscVG: E= " << ekin << " " << theParticle->GetParticleName() 
     << " P: " << proton << " xtot(b)= " << fTotalXsc/barn
     << " xel(b)= " <<  fElasticXsc/barn << " xinel(b)= " << fInelasticXsc/barn
     << G4endl;
  */
  return fTotalXsc;
}

Referenced by KaonNucleonXscGG().

◆ SCBMesonNucleonXscNS()

G4double G4HadronNucleonXsc::SCBMesonNucleonXscNS	(	const G4ParticleDefinition *	theParticle,
		const G4ParticleDefinition *	nucleon,
		G4double	ekin )

Definition at line 1056 of file G4HadronNucleonXsc.cc.

{
  G4double coeff(1.0);
  G4int pdg = std::abs(theParticle->GetPDGEncoding());
 
  // B+-0 anti
  if( pdg == 511 || pdg == 521 )
  {
    coeff = 0.610989;
  }
  // D+-0 anti
  else if( pdg == 411 || pdg == 421 )
  {
    coeff = 0.676568;
  }
  // Bs, antiBs
  else if( pdg == 531 )
  {
    coeff = 0.430989;
  }
  // Bc+-
  else if( pdg == 541 )
  {
    coeff = 0.287557;
  }
  // Ds+-
  else if( pdg == 431 )
  {
    coeff = 0.496568;
  }
  // etaC, J/Psi
  else if( pdg == 441 || pdg == 443 )
  {
    coeff = 0.353135;
  }
  // Upsilon
  else if( pdg == 553 )
  {
    coeff = 0.221978;
  }
  // eta
  else if( pdg == 221 )
  {
    coeff = 0.76;
  }
  // eta'
  else if( pdg == 331 )
  {
    coeff = 0.88;
  }
  fTotalXsc = coeff*HadronNucleonXscNS(thePiPlus, nucleon, ekin);
  fElasticXsc *= coeff;
  fInelasticXsc *= coeff;
  return fTotalXsc;
}

Referenced by HadronNucleonXsc().

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4HadronNucleonXsc()

◆ ~G4HadronNucleonXsc()

Member Function Documentation

◆ CoulombBarrier()

◆ CrossSectionDescription()

◆ GetCoulombBarrier()

◆ GetElasticHadronNucleonXsc()

◆ GetHadronNucleonXscEL()

◆ GetHadronNucleonXscNS()

◆ GetHadronNucleonXscPDG()

◆ GetHadronNucleonXscVU()

◆ GetHyperonNucleonXscNS()

◆ GetInelasticHadronNucleonXsc()

◆ GetKaonNucleonXscGG()

◆ GetTotalHadronNucleonXsc()

◆ HadronNucleonXsc()

◆ HadronNucleonXscEL()

◆ HadronNucleonXscNS()

◆ HadronNucleonXscPDG()

◆ HadronNucleonXscVU()

◆ HyperonNucleonXscNS()

◆ KaonNucleonXscGG()

◆ KaonNucleonXscNS()

◆ KaonNucleonXscVG()

◆ SCBMesonNucleonXscNS()