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

---

Detailed Description

Definition at line 37 of file G4XAqmTotal.hh.

---

Constructor & Destructor Documentation

G4XAqmTotal::G4XAqmTotal

(

)

Definition at line 57 of file G4XAqmTotal.cc.

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

G4XAqmTotal::~G4XAqmTotal

(

)

[virtual]

Definition at line 64 of file G4XAqmTotal.cc.

{ }

---

Member Function Documentation

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

Implements G4VCrossSectionSource.

Definition at line 80 of file G4XAqmTotal.cc.

References G4ParticleDefinition::GetAntiQuarkContent(), G4KineticTrack::GetDefinition(), and G4ParticleDefinition::GetQuarkContent().

Referenced by G4XAqmElastic::CrossSection().

{
  G4double sigma = 0.;

  // Get strangeness content
  G4ParticleDefinition* def1 = trk1.GetDefinition();
  G4int sTrk1 = def1->GetQuarkContent(3) + def1->GetAntiQuarkContent(3);
  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. * std::pow((2.0/3.0),nMesons) * (1. - 0.4 * sRatio1) * (1. - 0.4 * sRatio2) * millibarn;

  return sigma;
}

virtual const G4CrossSectionVector* G4XAqmTotal::GetComponents

(

)

const [inline, virtual]

Implements G4VCrossSectionSource.

Definition at line 51 of file G4XAqmTotal.hh.

{ return 0; }

G4bool G4XAqmTotal::IsValid

(

G4double

e

)

const [virtual]

Reimplemented from G4VCrossSectionSource.

Definition at line 132 of file G4XAqmTotal.cc.

References G4VCrossSectionSource::InLimits().

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

  return answer;
}

G4String G4XAqmTotal::Name

(

)

const [virtual]

Implements G4VCrossSectionSource.

Definition at line 124 of file G4XAqmTotal.cc.

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

G4bool G4XAqmTotal::operator!=

(

const G4XAqmTotal &

right

)

const

Definition at line 74 of file G4XAqmTotal.cc.

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

G4bool G4XAqmTotal::operator==

(

const G4XAqmTotal &

right

)

const

Definition at line 68 of file G4XAqmTotal.cc.

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

---

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

• G4XAqmTotal.hh
• G4XAqmTotal.cc

---