#include <G4CrossSectionPairGG.hh>

Inheritance 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 48 of file G4CrossSectionPairGG.cc.

References G4NistManager::Instance(), and G4VCrossSectionDataSet::verboseLevel.

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

G4CrossSectionPairGG::~G4CrossSectionPairGG ( )

virtual

Definition at line 57 of file G4CrossSectionPairGG.cc.

                                             {
    //   The cross section registry wil delete these
    //    delete theLowX;
    //    delete theHighX;
 }

Member Function Documentation

void G4CrossSectionPairGG::BuildPhysicsTable ( const G4ParticleDefinition & pDef )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 127 of file G4CrossSectionPairGG.cc.

References G4VCrossSectionDataSet::BuildPhysicsTable(), G4cout, G4endl, G4lrint(), G4NistManager::GetAtomicMassAmu(), G4VCrossSectionDataSet::GetElementCrossSection(), G4GlauberGribovCrossSection::GetInelasticGlauberGribov(), G4VCrossSectionDataSet::GetName(), G4ParticleDefinition::GetParticleName(), G4VCrossSectionDataSet::IsElementApplicable(), and G4VCrossSectionDataSet::verboseLevel.

Referenced by GetElementCrossSection().

                                                                              {
    theLowX->BuildPhysicsTable(pDef);
    theHighX->BuildPhysicsTable(pDef);
 
    if (verboseLevel > 0) {
       G4cout << "G4CrossSectionPairGG::BuildPhysicsTable "
             << theLowX->GetName() << "  " << theHighX->GetName() << G4endl;
    }
 
    G4ParticleDefinition * myDef = const_cast<G4ParticleDefinition*>(&pDef);
    std::vector<ParticleXScale>::iterator iter;
    iter = scale_factors.begin();
    while (iter != scale_factors.end() && (*iter).first != myDef) {
       ++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);
    }
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 63 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 200 of file G4CrossSectionPairGG.cc.

References G4cout, G4endl, G4VCrossSectionDataSet::GetName(), and python.hepunit::GeV.

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

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 88 of file G4CrossSectionPairGG.cc.

References BuildPhysicsTable(), G4cout, G4endl, G4lrint(), G4NistManager::GetAtomicMassAmu(), G4DynamicParticle::GetDefinition(), G4VCrossSectionDataSet::GetElementCrossSection(), G4GlauberGribovCrossSection::GetInelasticGlauberGribov(), G4DynamicParticle::GetKineticEnergy(), and G4VCrossSectionDataSet::verboseLevel.

 {
    G4double Xsec(0.);
 
    if (aParticle->GetKineticEnergy() < ETransition)
    {
       Xsec = theLowX->GetElementCrossSection(aParticle, ZZ, mat);
    } else {
 
       std::vector<ParticleXScale>::iterator iter = scale_factors.begin();
       G4ParticleDefinition * pDef = aParticle->GetDefinition();
       while (iter != scale_factors.end() && (*iter).first != pDef)
       {
          ++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;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 76 of file G4CrossSectionPairGG.cc.

References G4DynamicParticle::GetKineticEnergy(), and G4VCrossSectionDataSet::IsElementApplicable().

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

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation