G4UPiNuclearCrossSection Class Reference

#include <G4UPiNuclearCrossSection.hh>

Inheritance diagram for G4UPiNuclearCrossSection:

Public Member Functions
	G4UPiNuclearCrossSection ()
	~G4UPiNuclearCrossSection () override
G4bool	IsElementApplicable (const G4DynamicParticle *aParticle, G4int Z, const G4Material *) final
G4double	GetElasticCrossSection (const G4DynamicParticle *aParticle, G4int Z, G4int A) const
G4double	GetInelasticCrossSection (const G4DynamicParticle *aParticle, G4int Z, G4int A) const
void	BuildPhysicsTable (const G4ParticleDefinition &) final
void	DumpPhysicsTable (const G4ParticleDefinition &) final
void	CrossSectionDescription (std::ostream &) const final
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
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	GetElementCrossSection (const G4DynamicParticle *, G4int Z, 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 const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4int	GetVerboseLevel () const
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

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 45 of file G4UPiNuclearCrossSection.hh.

Constructor & Destructor Documentation

G4UPiNuclearCrossSection()

G4UPiNuclearCrossSection::G4UPiNuclearCrossSection ( )

explicit

Definition at line 61 of file G4UPiNuclearCrossSection.cc.

 : G4VCrossSectionDataSet("G4UPiNuclearCrossSection")
{
  isMaster = false;
  piPlus  = G4PionPlus::PionPlus();
  piMinus = G4PionMinus::PionMinus();
 
  aPower  = 0.75;
  elow    = 20.0*CLHEP::MeV;
}

~G4UPiNuclearCrossSection()

G4UPiNuclearCrossSection::~G4UPiNuclearCrossSection ( )

override

Definition at line 72 of file G4UPiNuclearCrossSection.cc.

{
  if(isMaster) {
    if(piPlusElastic) {
      piPlusElastic->clearAndDestroy();
      delete piPlusElastic;
      piPlusElastic = nullptr;
    }
    if(piPlusInelastic) {
      piPlusInelastic->clearAndDestroy();
      delete piPlusInelastic;
      piPlusInelastic = nullptr;
    }
    if(piMinusElastic) {
      piMinusElastic->clearAndDestroy();
      delete piMinusElastic;
      piMinusElastic = nullptr;
    }
    if(piMinusInelastic) {
      piMinusInelastic->clearAndDestroy();
      delete piMinusInelastic;
      piMinusInelastic = nullptr;
    }
  }
}

Member Function Documentation

BuildPhysicsTable()

void G4UPiNuclearCrossSection::BuildPhysicsTable ( const G4ParticleDefinition &  p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 169 of file G4UPiNuclearCrossSection.cc.

{
  //G4cout << "G4UPiNuclearCrossSection::BuildPhysicsTable for " 
  //     << p.GetParticleName() << "  " << piPlusElastic << G4endl; 
  if(piPlusElastic) { return; }
 
  if(&p != piPlus && &p != piMinus) { 
    G4ExceptionDescription ed;
    ed << "This cross section is applicable only to pions and not to " 
       << p.GetParticleName() << G4endl; 
    G4Exception("G4UPiNuclearCrossSection::BuildPhysicsTable", "had001", 
        FatalException, ed);
    return;
  }
 
  if(!piPlusElastic) { 
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&pionUXSMutex);
    if(!piPlusElastic) { 
#endif
      isMaster = true;
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&pionUXSMutex);
#endif
  }
  if(!isMaster) { return; }
  
  G4NistManager* nist = G4NistManager::Instance();
  G4Pow* g4pow = G4Pow::GetInstance();
  for(G4int i=0; i<NZ; ++i) {
    theA[i] = nist->GetAtomicMassAmu(theZ[i]);
  }
  for(G4int i=1; i<93; ++i) {
    APower[i] = g4pow->powA(nist->GetAtomicMassAmu(i),aPower);
  }
  G4int idx = 1;
  for(G4int i=3; i<93; ++i) {
    if(theZ[idx] == i) {
      idxZ[i] = -idx;
      ++idx;
    } else {
      idxZ[i] = idx;
    }
  }
 
  piPlusElastic    = new G4PhysicsTable();
  piPlusInelastic  = new G4PhysicsTable();
  piMinusElastic   = new G4PhysicsTable();
  piMinusInelastic = new G4PhysicsTable();
 
  LoadData();
  //G4cout << "### Initialisation done " << G4endl;
}

Referenced by G4BGGPionElasticXS::BuildPhysicsTable(), and G4BGGPionInelasticXS::BuildPhysicsTable().

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 590 of file G4UPiNuclearCrossSection.cc.

{
  outFile << "G4UPiNuclearCrossSection calculates the total, elastic and\n"
          << "inelastic cross sections for pion scattering from nuclei\n"
          << "heavier than hydrogen.  It is based on the Barashenkov\n"
          << "parameterization and is valid for all incident energies.\n";
}

DumpPhysicsTable()

void G4UPiNuclearCrossSection::DumpPhysicsTable ( const G4ParticleDefinition &  p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 154 of file G4UPiNuclearCrossSection.cc.

{
  if(&p == piPlus) {
    G4cout << "### G4UPiNuclearCrossSection Elastic data for pi+" << G4endl;
    G4cout << *piPlusElastic << G4endl;
    G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi+" << G4endl;
    G4cout << *piPlusInelastic << G4endl;
  } else if(&p == piMinus) {
    G4cout << "### G4UPiNuclearCrossSection Elastic data for pi-" << G4endl;
    G4cout << *piMinusElastic << G4endl;
    G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi-" << G4endl;
    G4cout << *piMinusInelastic << G4endl;
  }
}

GetElasticCrossSection()

G4double G4UPiNuclearCrossSection::GetElasticCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z, G4int  A  ) const

inline

Definition at line 106 of file G4UPiNuclearCrossSection.hh.

{
  const G4PhysicsTable* table = 
    (dp->GetDefinition() == piPlus) ? piPlusElastic : piMinusElastic; 
  return Interpolate(Z, A, dp->GetKineticEnergy(), table);
}

Referenced by G4BGGPionElasticXS::BuildPhysicsTable(), and G4BGGPionElasticXS::GetElementCrossSection().

GetInelasticCrossSection()

G4double G4UPiNuclearCrossSection::GetInelasticCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z, G4int  A  ) const

inline

Definition at line 115 of file G4UPiNuclearCrossSection.hh.

{
  const G4PhysicsTable* table = 
    (dp->GetDefinition() == piPlus) ? piPlusInelastic : piMinusInelastic; 
  return Interpolate(Z, A, dp->GetKineticEnergy(), table);
}

Referenced by G4BGGPionInelasticXS::BuildPhysicsTable(), and G4BGGPionInelasticXS::GetElementCrossSection().

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 99 of file G4UPiNuclearCrossSection.cc.

{
  return (1 < Z);
}

---

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

• G4UPiNuclearCrossSection.hh
• G4UPiNuclearCrossSection.cc