#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

G4ParticleDefinition *	FindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

G4double	GetTransversePionMass () const

G4double	GetMinStringMass () const

Detailed Description

Definition at line 50 of file G4CrossSectionPatch.hh.

Constructor & Destructor Documentation

G4CrossSectionPatch::G4CrossSectionPatch ( )

Definition at line 48 of file G4CrossSectionPatch.cc.

49 { }

G4CrossSectionPatch::~G4CrossSectionPatch ( )

virtual

Definition at line 52 of file G4CrossSectionPatch.cc.

53 { }

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(), nComponents, 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, G4XnpElastic, G4XnpTotal, and G4XNNTotal.

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 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	(	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;
 }

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

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

Public Member Functions

Protected Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation