G4CrossSectionPatch Class Reference

#include <G4CrossSectionPatch.hh>

Inheritance diagram for G4CrossSectionPatch:

Public Member Functions
	G4CrossSectionPatch ()
virtual	~G4CrossSectionPatch ()
G4bool	operator== (const G4CrossSectionPatch &right) const
G4bool	operator!= (const G4CrossSectionPatch &right) const
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual const G4CrossSectionVector *	GetComponents () const =0
virtual G4bool	IsValid (G4double e) const
Public Member Functions inherited from G4VCrossSectionSource
	G4VCrossSectionSource ()
virtual	~G4VCrossSectionSource ()
G4bool	operator== (const G4VCrossSectionSource &right) const
G4bool	operator!= (const G4VCrossSectionSource &right) const
virtual G4String	Name () const =0
virtual void	Print () const
virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual G4double	HighLimit () const
virtual G4double	LowLimit () const

Protected Member Functions
G4double	Transition (const G4KineticTrack &trk1, const G4KineticTrack &trk2, const G4VCrossSectionSource *comp1, const G4VCrossSectionSource *comp2) const
G4double	Transition (G4double ecm, G4double sigma1, G4double sigma2, G4double e1, G4double e2) const
Protected Member Functions inherited from G4VCrossSectionSource
G4bool	InLimits (G4double e, G4double eLow, G4double eHigh) const
const G4ParticleDefinition *	FindLightParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4double	FcrossX (G4double e, G4double e0, G4double sigma, G4double eParam, G4double power) const
const G4ParticleDefinition *	FindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Detailed Description

Definition at line 50 of file G4CrossSectionPatch.hh.

Constructor & Destructor Documentation

G4CrossSectionPatch()

G4CrossSectionPatch::G4CrossSectionPatch

(

)

Definition at line 48 of file G4CrossSectionPatch.cc.

{ }

Referenced by GetComponents(), operator!=(), and operator==().

~G4CrossSectionPatch()

G4CrossSectionPatch::~G4CrossSectionPatch ( )

virtual

Definition at line 52 of file G4CrossSectionPatch.cc.

{ }

Member Function Documentation

CrossSection()

G4double G4CrossSectionPatch::CrossSection ( const G4KineticTrack & trk1, const G4KineticTrack & trk2 ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 68 of file G4CrossSectionPatch.cc.

{
  // The cross section is provided by one of the components, according to their energy
  // validity range
 
  G4double crossSection = 0.;
  G4double ecm = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
 
  const G4CrossSectionVector* components = GetComponents();
  if (components != 0)
    {
      std::size_t nComponents = this->GetComponents()->size();
      
      for (std::size_t i=0; i<nComponents; ++i)
    {
          G4CrossSectionSourcePtr componentPtr = (*components)[i];
      G4VCrossSectionSource* component = componentPtr();
      if (component->IsValid(ecm))
        {
          crossSection = component->CrossSection(trk1,trk2);
        }
      else if (i < (nComponents - 1) )
        {
          G4CrossSectionSourcePtr nextPtr = (*components)[i+1];
          G4VCrossSectionSource* next = nextPtr();
          if (ecm > component->HighLimit() && ecm < next->LowLimit())
        {
          // Merge cross-sections in transition region between two validity ranges
          crossSection = Transition(trk1,trk2,component,next);
        }
        }
    }
    }
 
  return crossSection;
}

GetComponents()

virtual const G4CrossSectionVector * G4CrossSectionPatch::GetComponents ( ) const

pure virtual

Implements G4VCrossSectionSource.

Implemented in G4XNNElastic, G4XNNTotal, G4XnpElastic, and G4XnpTotal.

Referenced by CrossSection(), and IsValid().

IsValid()

G4bool G4CrossSectionPatch::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 107 of file G4CrossSectionPatch.cc.

{
  // The Patch is valid if any of its components are valid
  G4bool answer = false;
  const G4CrossSectionVector* components = GetComponents();
  if (components != 0)
    {
      std::size_t n = components->size();
      for (std::size_t i=0; i<n; ++i)
    {
      G4CrossSectionSourcePtr componentPtr = (*components)[i];
          G4VCrossSectionSource* component = componentPtr();
      if (component->IsValid(e)) 
        {
          answer = true;
          break;
        }
    }
    }
  return answer;
}

operator!=()

G4bool G4CrossSectionPatch::operator!=

(

const G4CrossSectionPatch &

right

)

const

Definition at line 62 of file G4CrossSectionPatch.cc.

{
  return (this != (G4CrossSectionPatch*) &right);
}

operator==()

G4bool G4CrossSectionPatch::operator==

(

const G4CrossSectionPatch &

right

)

const

Definition at line 56 of file G4CrossSectionPatch.cc.

{
  return (this == (G4CrossSectionPatch*) &right);
}

Transition() [1/2]

G4double G4CrossSectionPatch::Transition ( const G4KineticTrack & trk1, const G4KineticTrack & trk2, const G4VCrossSectionSource * comp1, const G4VCrossSectionSource * comp2 ) const

protected

Definition at line 130 of file G4CrossSectionPatch.cc.

{
  //Merge two cross sections in the transition region between their validity ranges
 
  G4double crossSection = 0.;
 
  G4double ecm = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4double sigma1 = comp1->CrossSection(trk1,trk2);
  G4double sigma2 = comp2->CrossSection(trk1,trk2);
  G4double denom = comp2->LowLimit() - comp1->HighLimit();
  G4double diff = ecm - comp1->HighLimit();
  if (denom > 0. && diff > 0.)
    {
      G4double ratio = diff / denom;
      crossSection = (1.- ratio) * sigma1 + ratio * sigma2;
    }
 
  return crossSection;
}

Referenced by CrossSection().

Transition() [2/2]

G4double G4CrossSectionPatch::Transition ( G4double ecm, G4double sigma1, G4double sigma2, G4double e1, G4double e2 ) const

protected

Definition at line 153 of file G4CrossSectionPatch.cc.

{
  //Merge two cross sections in the transition region between their validity ranges
 
  G4double crossSection = 0.;
  
  G4double denom = e2 - e1;
  G4double diff = ecm - e1;
  if (denom > 0. && diff > 0.)
    {
      G4double ratio = diff / denom;
      crossSection = (1.- ratio) * sigma1 + ratio * sigma2;
    }
 
  return crossSection;
}

---

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

• G4CrossSectionPatch.hh
• G4CrossSectionPatch.cc