G4XAqmTotal Class Reference

#include <G4XAqmTotal.hh>

Inheritance diagram for G4XAqmTotal:

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

Additional Inherited Members
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 37 of file G4XAqmTotal.hh.

Constructor & Destructor Documentation

G4XAqmTotal()

G4XAqmTotal::G4XAqmTotal

(

)

Definition at line 58 of file G4XAqmTotal.cc.

{ 
  // As a first approximation the model is assumed to be valid over 
  // the entire energy range
}

~G4XAqmTotal()

G4XAqmTotal::~G4XAqmTotal ( )

virtual

Definition at line 65 of file G4XAqmTotal.cc.

{ }

Member Function Documentation

CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 81 of file G4XAqmTotal.cc.

{
  G4double sigma = 0.;
 
  // Get strangeness content
  const G4ParticleDefinition* def1 = trk1.GetDefinition();
  G4int sTrk1 = def1->GetQuarkContent(3) + def1->GetAntiQuarkContent(3);
  const G4ParticleDefinition* def2 = trk2.GetDefinition();
  G4int sTrk2 = def2->GetQuarkContent(3) + def2->GetAntiQuarkContent(3);
  
  // Get non-strange quark content
  G4int qTrk1 = def1->GetQuarkContent(1) + def1->GetAntiQuarkContent(1) +
                def1->GetQuarkContent(2) + def1->GetAntiQuarkContent(2) +
                def1->GetQuarkContent(4) + def1->GetAntiQuarkContent(4) +
                def1->GetQuarkContent(5) + def1->GetAntiQuarkContent(5) +
                def1->GetQuarkContent(6) + def1->GetAntiQuarkContent(6);
 
  G4int qTrk2 = def2->GetQuarkContent(1) + def2->GetAntiQuarkContent(1) +
                def2->GetQuarkContent(2) + def2->GetAntiQuarkContent(2) +
                def2->GetQuarkContent(4) + def2->GetAntiQuarkContent(4) +
                def2->GetQuarkContent(5) + def2->GetAntiQuarkContent(5) +
                def2->GetQuarkContent(6) + def2->GetAntiQuarkContent(6);
 
  G4double sRatio1 = 0.;
  if (qTrk1 != 0) sRatio1 = sTrk1 / qTrk1;
  
  G4double sRatio2 = 0.;
  if (qTrk2 != 0) sRatio2 = sTrk2 / qTrk2;
  
  // Calculate the number of colliding mesons
  G4int nMesons = 0;
  G4int nQ1 = sTrk1 + qTrk1;
 
  if (nQ1 == 2) nMesons++;
  G4int nQ2 = sTrk2 + qTrk2;
  if (nQ2 == 2) nMesons++;
 
  // Cross-section (units to be checked!)
  sigma = 40. * G4Pow::GetInstance()->powN((2.0/3.0),nMesons) * (1. - 0.4 * sRatio1) * (1. - 0.4 * sRatio2) * millibarn;
 
  return sigma;
}

Referenced by G4XAqmElastic::CrossSection().

GetComponents()

virtual const G4CrossSectionVector * G4XAqmTotal::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 51 of file G4XAqmTotal.hh.

{ return 0; }

IsValid()

G4bool G4XAqmTotal::IsValid ( G4double  e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 133 of file G4XAqmTotal.cc.

{
  G4bool answer = InLimits(e,_lowLimit,_highLimit);
 
  return answer;
}

Name()

G4String G4XAqmTotal::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 125 of file G4XAqmTotal.cc.

{
  G4String name("AqmTotalCrossSection");
  return name;
}

operator!=()

G4bool G4XAqmTotal::operator!=

(

const G4XAqmTotal &

right

)

const

Definition at line 75 of file G4XAqmTotal.cc.

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

operator==()

G4bool G4XAqmTotal::operator==

(

const G4XAqmTotal &

right

)

const

Definition at line 69 of file G4XAqmTotal.cc.

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

---

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

• G4XAqmTotal.hh
• G4XAqmTotal.cc