#include <G4CrossSectionPairGG.hh>

Inheritance diagram for G4CrossSectionPairGG:

Collaboration diagram for G4CrossSectionPairGG:

Public Member Functions
	G4CrossSectionPairGG (G4VCrossSectionDataSet *low, G4double Etransit)

virtual	~G4CrossSectionPairGG ()

virtual void	CrossSectionDescription (std::ostream &) const

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat=0)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=0)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm=0, const G4Material *mat=0)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

virtual G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=0, const G4Element elm=0, const G4Material mat=0)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

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 51 of file G4CrossSectionPairGG.hh.

Constructor & Destructor Documentation

G4CrossSectionPairGG::G4CrossSectionPairGG	(	G4VCrossSectionDataSet *	low,
		G4double	Etransit
	)

Definition at line 50 of file G4CrossSectionPairGG.cc.

                          :
       G4VCrossSectionDataSet("G4CrossSectionPairGG"), theLowX(low), ETransition(
             Etransit) {
    NistMan = G4NistManager::Instance();
    theHighX = new G4ComponentGGHadronNucleusXsc();
    verboseLevel = 0;
 }

Here is the call graph for this function:

G4CrossSectionPairGG::~G4CrossSectionPairGG ( )

virtual

Definition at line 59 of file G4CrossSectionPairGG.cc.

                                             {
    delete theHighX;
    // The cross section registry will delete theLowX
 }

Member Function Documentation

void G4CrossSectionPairGG::BuildPhysicsTable ( const G4ParticleDefinition & pDef )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 128 of file G4CrossSectionPairGG.cc.

                                                                              {
    theLowX->BuildPhysicsTable(pDef);
    theHighX->BuildPhysicsTable(pDef);
 
    if (verboseLevel > 0) {
       G4cout << "G4CrossSectionPairGG::BuildPhysicsTable "
             << theLowX->GetName() << "  " << theHighX->GetName() << G4endl;
    }
 
    const G4ParticleDefinition * myDef = &pDef;
    std::vector<ParticleXScale>::iterator iter;
    iter = scale_factors.begin();
    while (iter != scale_factors.end() && (*iter).first != myDef)  /* Loop checking, 08.01.2016, W. Pokorski */
      {
        ++iter;
      }
 
    //  new particle, initialise
 
    G4Material* mat = 0;
 
    if (iter == scale_factors.end()) {
       XS_factors factors(93);
       G4ThreeVector mom(0.0, 0.0, 1.0);
       G4DynamicParticle DynPart(myDef, mom, ETransition); // last is kinetic Energy
 
       if (verboseLevel > 0) {
          G4cout << "G4CrossSectionPairGG::BuildPhysicsTable for particle "
                << pDef.GetParticleName() << G4endl;
       }
       for (G4int aZ = 1; aZ < 93; ++aZ) {
          factors[aZ] = 1.; // default, to give reasonable value if only high is applicable
          G4int AA = G4lrint(NistMan->GetAtomicMassAmu(aZ));
          G4bool isApplicable = theLowX->IsElementApplicable(&DynPart, aZ,
                mat) && (aZ > 1);
 
          if (isApplicable) {
             factors[aZ] = theLowX->GetElementCrossSection(&DynPart, aZ, mat)
                               / theHighX->GetInelasticGlauberGribov(&DynPart, aZ, AA);
 
          }
          if (verboseLevel > 0) {
             G4cout << "Z=" << aZ << ",  A=" << AA << ", scale="
                   << factors[aZ];
             if (verboseLevel == 1) {
                G4cout << G4endl;
             } else {
                if (isApplicable) {
                   G4cout << ",  low / high "
                         << theLowX->GetElementCrossSection(&DynPart, aZ,
                               mat) << "  "
                               << theHighX->GetInelasticGlauberGribov(&DynPart,
                                     aZ, AA) << G4endl;
                } else {
                   G4cout << ",   N/A" << G4endl;
                }
             }
          }
       }
       ParticleXScale forPart(myDef, factors);
       scale_factors.push_back(forPart);
    }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 64 of file G4CrossSectionPairGG.cc.

                                  {
    outFile << "G4CrossSectionPairGG is used to add the relativistic rise to\n"
          << "hadronic cross section data sets above a given energy.  In this\n"
          << "case, the Glauber-Gribov cross section is used above 91 GeV.\n"
          << "At this energy the low energy cross section is smoothly joined\n"
          << "to the high energy cross section.  Below 91 GeV the Barashenkov\n"
          << "cross section is used for pions (G4PiNuclearCrossSection), the\n"
          << "Axen-Wellisch cross section is used for protons\n"
          << "(G4ProtonInelasticCrossSection), and the Wellisch-Laidlaw cross\n"
          << "section is used for neutrons (G4NeutronInelasticCrossSection).\n";
 }

void G4CrossSectionPairGG::DumpPhysicsTable ( const G4ParticleDefinition & )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 202 of file G4CrossSectionPairGG.cc.

                                                                        {
    G4cout << std::setw(24) << " " << " G4CrossSectionPairGG: "
          << theLowX->GetName() << " cross sections " << G4endl;
    G4cout << std::setw(27) << " " << "below " << ETransition / GeV
          << " GeV, Glauber-Gribov above " << G4endl;
 }

Here is the call graph for this function:

G4double G4CrossSectionPairGG::GetElementCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		const G4Material *	mat = `0`
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4CrossSectionPairGG.cc.

 {
    G4double Xsec(0.);
 
    if (aParticle->GetKineticEnergy() < ETransition)
    {
       Xsec = theLowX->GetElementCrossSection(aParticle, ZZ, mat);
    } else {
 
       std::vector<ParticleXScale>::iterator iter = scale_factors.begin();
       const G4ParticleDefinition * pDef = aParticle->GetDefinition();
       while (iter != scale_factors.end() && (*iter).first != pDef)  /* Loop checking, 08.01.2016, W. Pokorski */
       {
          ++iter;
       }
       if (iter != scale_factors.end() )
       {
          G4int AA = G4lrint(NistMan->GetAtomicMassAmu(ZZ));
          Xsec = theHighX->GetInelasticGlauberGribov(aParticle, ZZ, AA)
                             * (*iter).second[ZZ];
          if (verboseLevel > 2)
          {
             G4cout << " scaling .." << ZZ << " " << AA << " "
                   << (*iter).second[ZZ] << " "
                   << theHighX->GetInelasticGlauberGribov(aParticle, ZZ, AA)
                   << "  " << Xsec << G4endl;
          }
       } else {
          // BuildPhysicsTable not done for pDef=aParticle->GetDefinition
          //  build table, and recurse
          BuildPhysicsTable(*pDef);
          Xsec=GetElementCrossSection(aParticle, ZZ, mat);
       }
    }
 
    return Xsec;
 }

Here is the call graph for this function:

G4bool G4CrossSectionPairGG::IsElementApplicable	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		const G4Material *	mat = `0`
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 77 of file G4CrossSectionPairGG.cc.

                                                                           {
    G4bool isApplicable(false);
    G4double Ekin = aParticle->GetKineticEnergy();
    if (Ekin <= ETransition) {
       isApplicable = theLowX->IsElementApplicable(aParticle, Z, mat);
    } else if (Z > 1) {
       isApplicable = true;
    }
    return isApplicable;
 }

Here is the call graph for this function:

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

source/geant4.10.03.p03/source/processes/hadronic/cross_sections/include/G4CrossSectionPairGG.hh
source/geant4.10.03.p03/source/processes/hadronic/cross_sections/src/G4CrossSectionPairGG.cc

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation