G4BGGPionElasticXS Class Reference

#include <G4BGGPionElasticXS.hh>

Inheritance diagram for G4BGGPionElasticXS:

Public Member Functions
	G4BGGPionElasticXS (const G4ParticleDefinition *)
	~G4BGGPionElasticXS () final
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *) final
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm, const G4Material *mat) final
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat) final
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr) final
void	BuildPhysicsTable (const G4ParticleDefinition &) final
void	CrossSectionDescription (std::ostream &) const final
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
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	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	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
G4String	name

Detailed Description

Definition at line 64 of file G4BGGPionElasticXS.hh.

Constructor & Destructor Documentation

G4BGGPionElasticXS()

G4BGGPionElasticXS::G4BGGPionElasticXS ( const G4ParticleDefinition * p )

explicit

Definition at line 68 of file G4BGGPionElasticXS.cc.

 : G4VCrossSectionDataSet("BarashenkovGlauberGribov") 
{
  verboseLevel = 0;
  fGlauberEnergy = 91.*GeV;
  fLowEnergy = 20.*MeV;
  fLowestEnergy = 1.*MeV;
  SetMinKinEnergy(0.0);
  SetMaxKinEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
 
  fPion = nullptr;
  fGlauber = nullptr;
  fHadron  = nullptr;
 
  fG4pow   = G4Pow::GetInstance();
 
  theProton= G4Proton::Proton();
  thePiPlus= G4PionPlus::PionPlus();
  isPiplus = (p == thePiPlus);
  isMaster = false;
  SetForAllAtomsAndEnergies(true);
}

~G4BGGPionElasticXS()

G4BGGPionElasticXS::~G4BGGPionElasticXS ( )

final

Definition at line 93 of file G4BGGPionElasticXS.cc.

{
  delete fHadron;
}

Member Function Documentation

BuildPhysicsTable()

void G4BGGPionElasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 175 of file G4BGGPionElasticXS.cc.

{
  if(fPion) { return; }
  if(verboseLevel > 1) {
    G4cout << "G4BGGPionElasticXS::BuildPhysicsTable for " 
           << p.GetParticleName() << G4endl;
  } 
  if(&p == G4PionPlus::PionPlus() || &p == G4PionMinus::PionMinus()) {
    isPiplus = (&p == G4PionPlus::PionPlus());
  } else {
    G4ExceptionDescription ed;
    ed << "This BGG cross section is applicable only to pions and not to " 
       << p.GetParticleName() << G4endl; 
    G4Exception("G4BGGPionElasticXS::BuildPhysicsTable", "had001", 
                FatalException, ed);
    return;
  }
 
  fPion    = new G4UPiNuclearCrossSection();
  fGlauber = new G4ComponentGGHadronNucleusXsc();
  fHadron  = new G4HadronNucleonXsc();
 
  fPion->BuildPhysicsTable(p);
 
  if(0 == theA[0]) { 
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&pionElasticXSMutex);
    if(0 == theA[0]) { 
#endif
      isMaster = true;
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&pionElasticXSMutex);
#endif
  } else {
    return;
  }
 
  if(isMaster && 0 == theA[0]) {
 
    theA[0] = theA[1] = 1;
    G4ThreeVector mom(0.0,0.0,1.0);
    G4DynamicParticle dp(thePiPlus, mom, fGlauberEnergy);
 
    G4NistManager* nist = G4NistManager::Instance();
 
    G4double csup, csdn;
    for(G4int iz=2; iz<93; ++iz) {
 
      G4int A = G4lrint(nist->GetAtomicMassAmu(iz));
      theA[iz] = A;
 
      csup = fGlauber->GetElasticGlauberGribov(&dp, iz, A);
      csdn = fPion->GetElasticCrossSection(&dp, iz, A);
      theGlauberFacPiPlus[iz] = csdn/csup;
    }
 
    dp.SetDefinition(G4PionMinus::PionMinus());
    for(G4int iz=2; iz<93; ++iz) {
      csup = fGlauber->GetElasticGlauberGribov(&dp, iz, theA[iz]);
      csdn = fPion->GetElasticCrossSection(&dp, iz, theA[iz]);
      theGlauberFacPiMinus[iz] = csdn/csup;
 
      if(verboseLevel > 0) {
        G4cout << "Z= " << iz <<  "  A= " << theA[iz] 
               << " factorPiPlus= " << theGlauberFacPiPlus[iz] 
               << " factorPiMinus= " << theGlauberFacPiMinus[iz] 
               << G4endl; 
      }
    }
    theCoulombFacPiPlus[1] = 1.0;
    theCoulombFacPiMinus[1]= 1.0;
    dp.SetKineticEnergy(fLowEnergy);
    dp.SetDefinition(thePiPlus);
    for(G4int iz=2; iz<93; ++iz) {
      theCoulombFacPiPlus[iz] = fPion->GetElasticCrossSection(&dp, iz, theA[iz])
        /CoulombFactorPiPlus(fLowEnergy, iz);
    }
    dp.SetDefinition(G4PionMinus::PionMinus());
    for(G4int iz=2; iz<93; ++iz) {
      theCoulombFacPiMinus[iz] = fPion->GetElasticCrossSection(&dp, iz, theA[iz])
        /FactorPiMinus(fLowEnergy);
 
      if(verboseLevel > 0) {
        G4cout << "Z= " << iz <<  "  A= " << theA[iz] 
               << " CoulombFactorPiPlus= " << theCoulombFacPiPlus[iz] 
               << " CoulombFactorPiMinus= " << theCoulombFacPiMinus[iz] 
               << G4endl; 
      }
    }
  }
}

Referenced by G4LightIonQMDReaction::G4LightIonQMDReaction(), and G4QMDReaction::G4QMDReaction().

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 286 of file G4BGGPionElasticXS.cc.

{
  outFile << "The Barashenkov-Glauber-Gribov cross section handles elastic\n"
          << "scattering of pions from nuclei at all energies. The\n"
          << "Barashenkov parameterization is used below 91 GeV and the\n"
          << "Glauber-Gribov parameterization is used above 91 GeV.\n";
}

GetElementCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 119 of file G4BGGPionElasticXS.cc.

{
  // this method should be called only for Z > 1
  G4double cross = 0.0;
  G4double ekin = std::max(dp->GetKineticEnergy(), fLowestEnergy);
  G4int Z = std::min(ZZ, 92);
  if(1 == Z) {
    cross = 1.0115*GetIsoCrossSection(dp,1,1);
  } else {
    if(ekin <= fLowEnergy) {
      cross = (isPiplus) ? theCoulombFacPiPlus[Z]*CoulombFactorPiPlus(ekin, Z)
        : theCoulombFacPiMinus[Z]*FactorPiMinus(ekin);
    } else if(ekin > fGlauberEnergy) {
      cross = (isPiplus) ? theGlauberFacPiPlus[Z] : theGlauberFacPiMinus[Z];
      cross *= fGlauber->GetElasticGlauberGribov(dp, Z, theA[Z]);
    } else {
      cross = fPion->GetElasticCrossSection(dp, Z, theA[Z]);
    }
  }
  if(verboseLevel > 1) {
    G4cout << "G4BGGPionElasticXS::GetElementCrossSection  for "
           << dp->GetDefinition()->GetParticleName()
           << "  Ekin(GeV)= " << dp->GetKineticEnergy()
           << " in nucleus Z= " << Z << "  A= " << theA[Z]
           << " XS(b)= " << cross/barn 
           << G4endl;
  }
  return cross;
}

Referenced by G4LightIonQMDReaction::ApplyYourself(), and G4QMDReaction::ApplyYourself().

GetIsoCrossSection()

G4double G4BGGPionElasticXS::GetIsoCrossSection ( const G4DynamicParticle * dp, G4int Z, G4int A, const G4Isotope * iso = nullptr, const G4Element * elm = nullptr, const G4Material * mat = nullptr )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 151 of file G4BGGPionElasticXS.cc.

{
  // this method should be called only for Z = 1
  fHadron->HadronNucleonXscNS(dp->GetDefinition(), theProton, 
                              dp->GetKineticEnergy());
  G4double cross = A*fHadron->GetElasticHadronNucleonXsc();
 
  if(verboseLevel > 1) {
    G4cout << "G4BGGPionElasticXS::GetIsoCrossSection  for "
           << dp->GetDefinition()->GetParticleName()
           << "  Ekin(GeV)= " << dp->GetKineticEnergy()
           << " in nucleus Z= " << Z << "  A= " << A
           << " XS(b)= " << cross/barn 
           << G4endl;
  }
  return cross;
}

Referenced by GetElementCrossSection().

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 101 of file G4BGGPionElasticXS.cc.

{
  return true;
}

IsIsoApplicable()

G4bool G4BGGPionElasticXS::IsIsoApplicable ( const G4DynamicParticle * , G4int Z, G4int A, const G4Element * elm, const G4Material * mat )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 109 of file G4BGGPionElasticXS.cc.

{
  return (1 == Z);
}

---

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

• G4BGGPionElasticXS.hh
• G4BGGPionElasticXS.cc