#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 void	Print () const

virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) 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

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

G4double	GetTransversePionMass () const

G4double	GetMinStringMass () const

Detailed Description

Definition at line 46 of file G4XAnnihilationChannel.hh.

Constructor & Destructor Documentation

G4XAnnihilationChannel::G4XAnnihilationChannel ( )

Definition at line 37 of file G4XAnnihilationChannel.cc.

References DBL_MAX.

                                               : 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::G4XAnnihilationChannel	(	const G4ParticleDefinition *	resDefinition,
		const G4ResonanceWidth &	resWidths,
		const G4ResonancePartialWidth &	resPartWidths,
		const G4String &	partWidthLabel
	)

Definition at line 47 of file G4XAnnihilationChannel.cc.

References DBL_MAX, G4ParticleDefinition::GetParticleName(), G4ResonancePartialWidth::MassDependentWidth(), G4ResonanceWidth::MassDependentWidth(), and G4ResonanceNames::ShortName().

   : 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 ( )

virtual

Definition at line 70 of file G4XAnnihilationChannel.cc.

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

Member Function Documentation

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

virtual

Implements G4VCrossSectionSource.

Definition at line 91 of file G4XAnnihilationChannel.cc.

References G4KineticTrack::Get4Momentum(), G4KineticTrack::GetDefinition(), G4ParticleDefinition::GetPDGiSpin(), G4ParticleDefinition::GetPDGMass(), python.hepunit::hbarc_squared, python.hepunit::pi, and width.

 {
   G4double sigma = 0.;
   G4double eCM = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
 
   G4ParticleDefinition* def1 = trk1.GetDefinition();
   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;
 }

virtual const G4CrossSectionVector* G4XAnnihilationChannel::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 64 of file G4XAnnihilationChannel.hh.

64 { return 0; }

G4bool G4XAnnihilationChannel::IsValid ( G4double e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 144 of file G4XAnnihilationChannel.cc.

References G4VCrossSectionSource::InLimits().

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

G4String G4XAnnihilationChannel::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 136 of file G4XAnnihilationChannel.cc.

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

G4bool G4XAnnihilationChannel::operator!= ( const G4XAnnihilationChannel & right ) const

Definition at line 85 of file G4XAnnihilationChannel.cc.

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

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:

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation