G4ChargeExchangeXS Class Reference

#include <G4ChargeExchangeXS.hh>

Inheritance diagram for G4ChargeExchangeXS:

Public Member Functions
	G4ChargeExchangeXS ()
	~G4ChargeExchangeXS () override=default
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *) final
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *) final
void	CrossSectionDescription (std::ostream &) const final
const G4ParticleDefinition *	SampleSecondaryType (const G4ParticleDefinition *, const G4Material *, G4int Z, G4int A, G4double etot)
G4double	GetPartialPionXS (G4int idx) const
void	SetEnergyLimit (G4double val)
void	SetCrossSectionFactor (G4double val)
G4double	GetCrossSectionFactor () const
G4bool	isPion () const
G4double	SampleTforPion (const G4double etot, const G4double tmax) const
G4ChargeExchangeXS &	operator= (const G4ChargeExchangeXS &right)=delete
	G4ChargeExchangeXS (const G4ChargeExchangeXS &)=delete
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
bool	ForAllAtomsAndEnergies () const
void	SetForAllAtomsAndEnergies (G4bool val)
const G4String &	GetName () const
void	SetName (const G4String &nam)
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel {0}
G4String	name

Detailed Description

Definition at line 65 of file G4ChargeExchangeXS.hh.

Constructor & Destructor Documentation

G4ChargeExchangeXS() [1/2]

G4ChargeExchangeXS::G4ChargeExchangeXS

(

)

Definition at line 77 of file G4ChargeExchangeXS.cc.

{
  if (verboseLevel > 1) {
    G4cout  << "G4ChargeExchangeXS::G4ChargeExchangeXS" << G4endl;
  }
  fMassPi = G4PionPlus::PionPlus()->GetPDGMass();
  g4calc = G4Pow::GetInstance();
  auto table = G4ParticleTable::GetParticleTable();
  const G4String nam[5] = {"pi0", "eta", "eta_prime", "omega", "f2(1270)"};
  for (G4int i=0; i<5; ++i) {
    fPionSecPD[i] = table->FindParticle(nam[i]);
    if (nullptr == fPionSecPD[i]) {
      G4ExceptionDescription ed;
      ed << "### meson " << nam[i] << " is not found out in the particle table";
      G4Exception("G4ChargeExchangeXS::G4ChargeExchangeXS()","had044",
                  FatalException, ed,"");
    }
  }
}

Referenced by G4ChargeExchangeXS(), and operator=().

~G4ChargeExchangeXS()

G4ChargeExchangeXS::~G4ChargeExchangeXS ( )

overridedefault

G4ChargeExchangeXS() [2/2]

G4ChargeExchangeXS::G4ChargeExchangeXS ( const G4ChargeExchangeXS & )

delete

Member Function Documentation

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 98 of file G4ChargeExchangeXS.cc.

{
  outFile << "G4ChargeExchangeXS calculates charge exchange cross section for "
          << "pi+, pi-, K+, K-, KL\n";
}

GetCrossSectionFactor()

G4double G4ChargeExchangeXS::GetCrossSectionFactor ( ) const

inline

Definition at line 91 of file G4ChargeExchangeXS.hh.

{ return fFactor; };

GetElementCrossSection()

G4double G4ChargeExchangeXS::GetElementCrossSection ( const G4DynamicParticle * dp, G4int Z, const G4Material * mat )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 111 of file G4ChargeExchangeXS.cc.

{
  G4double pE = dp->GetTotalEnergy();
  return (pE > fEnergyLimit) ?
    GetCrossSection(dp->GetDefinition(), mat, Z, pE) : 0.0;
}

GetPartialPionXS()

G4double G4ChargeExchangeXS::GetPartialPionXS

(

G4int

idx

)

const

Definition at line 314 of file G4ChargeExchangeXS.cc.

{
  G4double res = 0.0;
  if (0 == idx) { res = fXSecPion[0]; }
  else if (0 < idx && 5 > idx) {
    res = fXSecPion[idx] - fXSecPion[idx - 1];
  }
  return res;
}

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 104 of file G4ChargeExchangeXS.cc.

{
  return true;
}

isPion()

G4bool G4ChargeExchangeXS::isPion ( ) const

inline

Definition at line 93 of file G4ChargeExchangeXS.hh.

{ return findex >= 0; };

operator=()

G4ChargeExchangeXS & G4ChargeExchangeXS::operator= ( const G4ChargeExchangeXS & right )

delete

SampleSecondaryType()

const G4ParticleDefinition * G4ChargeExchangeXS::SampleSecondaryType	(	const G4ParticleDefinition *	part,
		const G4Material *	mat,
		G4int	Z,
		G4int	A,
		G4double	etot )

Definition at line 220 of file G4ChargeExchangeXS.cc.

{
  // index of pion partial x-section
  findex = -1;
  // recompute x-section for the element in complex material
  GetCrossSection(part, mat, Z, etot);
  
  const G4ParticleDefinition* pd = nullptr;
  G4int pdg = std::abs(part->GetPDGEncoding());  
  G4cout << pdg << G4endl;
  // pi- + p /  pi+ + n  
  if (pdg == 211) {
    pd = fPionSecPD[0];
    G4double x = fXSecPion[4]*G4UniformRand();
    for (findex = 0; findex < 5; ++findex) {
      if (x <= fXSecPion[findex]) {
        pd = fPionSecPD[findex];
    break;
      }
    }
  }
 
  // K- + p /  K+ + n 
  // Equal opportunity of producing k-short and k-long
  else if (pdg == 321) {
    if (G4UniformRand() >= 0.5) {
      pd = G4KaonZeroLong::KaonZeroLong();
    }
    else {
      pd = G4KaonZeroShort::KaonZeroShort();
    }
  }
 
  // KL + nucleus
  else if (pdg == 130) {
    G4double prob = (G4double)Z/(G4double)A;
    if (G4UniformRand() >= prob) {
      pd = G4KaonMinus::KaonMinus();
    }
    else {
      pd = G4KaonPlus::KaonPlus();
    }
  }
  if (verboseLevel > 1) {
    G4cout << "G4ChargeExchangeXS::SampleSecondaryType for "
       << pd->GetParticleName() << "  findex=" << findex
       << G4endl;
  }
  return pd;
} 

SampleTforPion()

G4double G4ChargeExchangeXS::SampleTforPion	(	const G4double	etot,
		const G4double	tmax ) const

Definition at line 273 of file G4ChargeExchangeXS.cc.

{
  G4double tmax = ltmax/GeV2;
  G4double tM = CLHEP::proton_mass_c2;
  G4double logX = G4Log((tM*tM + 2*tM*etot +  fMassPi*fMassPi)*inv1e7);
 
  G4double gl = pG0[findex] + pG1[findex]*logX;
  G4double cl = pC0[findex] + pC1[findex]*logX;
  G4double gc = gl*cl;
 
  G4double t{0};
  G4double sigmaMax = (gc > 0.0) ? gl*G4Exp(-(gl - 1.0)/gl) : 1.0;
  for (G4int i = 0; i < 100000; ++i) {
    t = tmax*G4UniformRand();
    G4double sigma = (1.0 + gc*t)*G4Exp(-cl*t);
    if (G4UniformRand()*sigmaMax <= sigma) {
      return t*GeV2;
    }
  }
  return 0.0;
}

SetCrossSectionFactor()

void G4ChargeExchangeXS::SetCrossSectionFactor ( G4double val )

inline

Definition at line 89 of file G4ChargeExchangeXS.hh.

{ fFactor = val; };

SetEnergyLimit()

void G4ChargeExchangeXS::SetEnergyLimit ( G4double val )

inline

Definition at line 87 of file G4ChargeExchangeXS.hh.

{ fEnergyLimit = val; };

---

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

• G4ChargeExchangeXS.hh
• G4ChargeExchangeXS.cc