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
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

---

Detailed Description

Definition at line 50 of file G4CrossSectionPatch.hh.

---

Constructor & Destructor Documentation

G4CrossSectionPatch::G4CrossSectionPatch

(

)

Definition at line 48 of file G4CrossSectionPatch.cc.

{ }

G4CrossSectionPatch::~G4CrossSectionPatch

(

)

[virtual]

Definition at line 52 of file G4CrossSectionPatch.cc.

{ }

---

Member Function Documentation

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

Implements G4VCrossSectionSource.

Definition at line 68 of file G4CrossSectionPatch.cc.

References G4VCrossSectionSource::CrossSection(), G4KineticTrack::Get4Momentum(), GetComponents(), G4VCrossSectionSource::HighLimit(), G4VCrossSectionSource::IsValid(), G4VCrossSectionSource::LowLimit(), and Transition().

{
  // 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)
    {
      G4int nComponents = this->GetComponents()->size();
//      G4int nValid = 0;
      
      G4int i;
      for (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;
}

virtual const G4CrossSectionVector* G4CrossSectionPatch::GetComponents

(

)

const [pure virtual]

Implements G4VCrossSectionSource.

Implemented in G4XNNElastic, G4XNNTotal, G4XnpElastic, and G4XnpTotal.

Referenced by CrossSection(), and IsValid().

G4bool G4CrossSectionPatch::IsValid

(

G4double

e

)

const [virtual]

Reimplemented from G4VCrossSectionSource.

Definition at line 109 of file G4CrossSectionPatch.cc.

References GetComponents(), G4VCrossSectionSource::IsValid(), and CLHEP::detail::n.

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

G4bool G4CrossSectionPatch::operator!=

(

const G4CrossSectionPatch &

right

)

const

Definition at line 62 of file G4CrossSectionPatch.cc.

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

G4bool G4CrossSectionPatch::operator==

(

const G4CrossSectionPatch &

right

)

const

Definition at line 56 of file G4CrossSectionPatch.cc.

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

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

Definition at line 156 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;
}

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

Definition at line 133 of file G4CrossSectionPatch.cc.

References G4VCrossSectionSource::CrossSection(), G4KineticTrack::Get4Momentum(), G4VCrossSectionSource::HighLimit(), and G4VCrossSectionSource::LowLimit().

Referenced by CrossSection().

{
  //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;
}

---

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

• G4CrossSectionPatch.hh
• G4CrossSectionPatch.cc

---