G4XAnnihilationChannel Class Reference

#include <G4XAnnihilationChannel.hh>

Inheritance diagram for G4XAnnihilationChannel:

Public Member Functions
	G4XAnnihilationChannel ()
	G4XAnnihilationChannel (const G4ParticleDefinition *resDefinition, const G4ResonanceWidth &resWidths, const G4ResonancePartialWidth &resPartWidths, const G4String &partWidthLabel)
virtual	~G4XAnnihilationChannel ()
G4bool	operator== (const G4XAnnihilationChannel &right) const
G4bool	operator!= (const G4XAnnihilationChannel &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 46 of file G4XAnnihilationChannel.hh.

Constructor & Destructor Documentation

G4XAnnihilationChannel() [1/2]

G4XAnnihilationChannel::G4XAnnihilationChannel

(

)

Definition at line 37 of file G4XAnnihilationChannel.cc.

                                              : resonance(0)
{
    // As a first approximation the model is assumed to be valid over
      // the entire energy range
      lowLimit = 0.;
      highLimit = DBL_MAX;
      widthTable = 0;
      partWidthTable = 0;
}

G4XAnnihilationChannel() [2/2]

G4XAnnihilationChannel::G4XAnnihilationChannel	(	const G4ParticleDefinition *	resDefinition,
		const G4ResonanceWidth &	resWidths,
		const G4ResonancePartialWidth &	resPartWidths,
		const G4String &	partWidthLabel
	)

Definition at line 47 of file G4XAnnihilationChannel.cc.

  : resonance(resDefinition)
{ 
  // Get the tabulated mass-dependent widths for the resonance
  G4String resName = resonance->GetParticleName();
  // cout << "HPW "<<resName<<endl;
  G4String shortName = theNames.ShortName(resName);
  // cout << "HPW "<<shortName<<endl;
  // cout << "HPW "<<partWidthLabel<<endl;
 
  widthTable = resWidths.MassDependentWidth(shortName);
  partWidthTable = resPartWidths.MassDependentWidth(partWidthLabel);
 
  // As a first approximation the model is assumed to be valid over 
  // the entire energy range
  lowLimit = 0.;
  highLimit = DBL_MAX;
}

~G4XAnnihilationChannel()

G4XAnnihilationChannel::~G4XAnnihilationChannel ( )

virtual

Definition at line 70 of file G4XAnnihilationChannel.cc.

{
  if (widthTable) delete widthTable;
  widthTable = 0;
  if (partWidthTable) delete partWidthTable;
  partWidthTable = 0;
 }

Member Function Documentation

CrossSection()

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

virtual

Implements G4VCrossSectionSource.

Definition at line 91 of file G4XAnnihilationChannel.cc.

{
  G4double sigma = 0.;
  G4double eCM = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
 
  const G4ParticleDefinition* def1 = trk1.GetDefinition();
  const G4ParticleDefinition* def2 = trk2.GetDefinition();
 
  G4int J1 = def1->GetPDGiSpin();
  G4int J2 = def2->GetPDGiSpin();
  G4double m_1 = def1->GetPDGMass();
  G4double m_2 = def2->GetPDGMass();
 
  G4int JRes = resonance->GetPDGiSpin();
  G4double mRes = resonance->GetPDGMass();
 
  G4double branch = Branch(trk1,trk2);
  G4double width = VariableWidth(trk1,trk2);
  G4double cleb = NormalizedClebsch(trk1,trk2);
 
  G4double S = eCM * eCM;
  if (S == 0.) throw G4HadronicException(__FILE__, __LINE__, "G4XAnnihilationChannel::CrossSection - eCM = 0");
 
  G4double pCM = std::sqrt((S-(m_1+m_2)*(m_1+m_2))*(S-(m_1-m_2)*(m_1-m_2))/(4.*S));
 
  sigma = ( (JRes + 1.) / ( (J1 + 1) * (J2 + 1) ) 
        * pi / (pCM * pCM) * branch * width * width / 
        ( (eCM - mRes) * (eCM - mRes) + width * width / 4.0) * cleb * hbarc_squared);
 
//   G4cout << "SS " << branch<<" "<<sigma<<" "
//          << J1 <<" "
//   <<J2<<" "
//   <<m1<<" "
//   <<m2<<" "
//   <<JRes<<" "
//   <<mRes<<" "
//   <<wRes<<" "
//   <<width<<" "
//   <<cleb<<" "
//   <<G4endl;
  return sigma;
}

GetComponents()

virtual const G4CrossSectionVector * G4XAnnihilationChannel::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 64 of file G4XAnnihilationChannel.hh.

{ return 0; }

IsValid()

G4bool G4XAnnihilationChannel::IsValid ( G4double  e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 144 of file G4XAnnihilationChannel.cc.

{
  G4bool answer = InLimits(e,lowLimit,highLimit);
 
  return answer;
}

Name()

G4String G4XAnnihilationChannel::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 136 of file G4XAnnihilationChannel.cc.

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

operator!=()

G4bool G4XAnnihilationChannel::operator!=

(

const G4XAnnihilationChannel &

right

)

const

Definition at line 85 of file G4XAnnihilationChannel.cc.

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

operator==()

G4bool G4XAnnihilationChannel::operator==

(

const G4XAnnihilationChannel &

right

)

const

Definition at line 79 of file G4XAnnihilationChannel.cc.

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

---

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

• G4XAnnihilationChannel.hh
• G4XAnnihilationChannel.cc