G4XAnnihilationChannel Class Reference

#include <G4XAnnihilationChannel.hh>

Inheritance diagram for G4XAnnihilationChannel:

Collaboration 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

Private Member Functions
	G4XAnnihilationChannel (const G4XAnnihilationChannel &right)
const G4XAnnihilationChannel &	operator= (const G4XAnnihilationChannel &right)
G4double	Branch (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4double	VariableWidth (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4double	VariablePartialWidth (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4double	NormalizedClebsch (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Private Attributes
G4double	lowLimit
G4double	highLimit
G4Clebsch	clebsch
G4ResonanceNames	theNames
G4PhysicsVector *	widthTable
G4PhysicsVector *	partWidthTable
const G4ParticleDefinition *	resonance

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/3]

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

Here is the caller graph for this function:

G4XAnnihilationChannel() [2/3]

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

Here is the call graph for this function:

~G4XAnnihilationChannel()

G4XAnnihilationChannel::~G4XAnnihilationChannel ( )

virtual

Definition at line 70 of file G4XAnnihilationChannel.cc.

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

G4XAnnihilationChannel() [3/3]

G4XAnnihilationChannel::G4XAnnihilationChannel ( const G4XAnnihilationChannel &  right )

private

Member Function Documentation

Branch()

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

private

Definition at line 152 of file G4XAnnihilationChannel.cc.

{
  G4double w=VariableWidth(trk1,trk2);
  if(w==0) return 0;
  return VariablePartialWidth(trk1,trk2) / VariableWidth(trk1,trk2);
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the call graph for this function:

GetComponents()

virtual const G4CrossSectionVector* G4XAnnihilationChannel::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 64 of file G4XAnnihilationChannel.hh.

{ return 0; }

Here is the call graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

Name()

G4String G4XAnnihilationChannel::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 136 of file G4XAnnihilationChannel.cc.

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

Here is the caller graph for this function:

NormalizedClebsch()

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

private

Definition at line 199 of file G4XAnnihilationChannel.cc.

{
  G4double cleb = 0.;
  const G4ParticleDefinition* def1 = trk1.GetDefinition();
  const G4ParticleDefinition* def2 = trk2.GetDefinition();

  G4int iso31 = def1->GetPDGiIsospin3();
  G4int iso32 = def2->GetPDGiIsospin3();
  G4int iso3 = iso31 + iso32;
  G4int iso1 = def1->GetPDGiIsospin();
  G4int iso2 = def2->GetPDGiIsospin();

  G4int isoRes = resonance->GetPDGiIsospin();
  
  if (isoRes < iso3) return 0.;
  if ((iso1*iso2) == 0) return 1.;

  cleb = clebsch.NormalizedClebschGordan(isoRes,iso3,iso1,iso2,iso31,iso32);

  // Special case: particle-antiparticle, charge-conjugated states have the same weight
  G4String type1 = def1->GetParticleType();
  G4String type2 = def2->GetParticleType();
  G4int anti = def1->GetPDGEncoding() * def2->GetPDGEncoding();
  G4int strangeness = resonance->GetQuarkContent(3) + resonance->GetAntiQuarkContent(3);
  if ( ((type1 == "baryon" && type2 == "baryon") ||(type1 == "meson" && type2 == "meson")) &&
       anti < 0 && strangeness == 0) 
    {
      if (def1->GetPDGEncoding() != -(def2->GetPDGEncoding())) cleb = 0.5 * cleb;
    }
       
  return cleb;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

G4bool G4XAnnihilationChannel::operator!=

(

const G4XAnnihilationChannel &

right

)

const

Definition at line 85 of file G4XAnnihilationChannel.cc.

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

operator=()

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

private

Here is the caller graph for this function:

operator==()

G4bool G4XAnnihilationChannel::operator==

(

const G4XAnnihilationChannel &

right

)

const

Definition at line 79 of file G4XAnnihilationChannel.cc.

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

VariablePartialWidth()

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

private

Definition at line 176 of file G4XAnnihilationChannel.cc.

{
  // Calculate mass dependent partial width of resonance, 
  // based on UrQMD tabulations

  G4double width(0);

  if (partWidthTable != 0)
  {
    G4double sqrtS = 0;
    G4bool dummy = false;
    sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
    width = partWidthTable->GetValue(sqrtS,dummy);
  }
  else
  {
    width = resonance->GetPDGWidth();
  }
  return width;
}

Here is the call graph for this function:

Here is the caller graph for this function:

VariableWidth()

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

private

Definition at line 160 of file G4XAnnihilationChannel.cc.

{
  // actual production width of resonance, depending on available energy.

  G4double width = resonance->GetPDGWidth();
  G4bool dummy = false;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  if (widthTable != 0) 
    {
      width = widthTable->GetValue(sqrtS,dummy);
    }
  return width;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

clebsch

G4Clebsch G4XAnnihilationChannel::clebsch

private

Definition at line 93 of file G4XAnnihilationChannel.hh.

highLimit

G4double G4XAnnihilationChannel::highLimit

private

Definition at line 91 of file G4XAnnihilationChannel.hh.

lowLimit

G4double G4XAnnihilationChannel::lowLimit

private

Definition at line 90 of file G4XAnnihilationChannel.hh.

partWidthTable

G4PhysicsVector* G4XAnnihilationChannel::partWidthTable

private

Definition at line 98 of file G4XAnnihilationChannel.hh.

resonance

const G4ParticleDefinition* G4XAnnihilationChannel::resonance

private

Definition at line 101 of file G4XAnnihilationChannel.hh.

theNames

G4ResonanceNames G4XAnnihilationChannel::theNames

private

Definition at line 94 of file G4XAnnihilationChannel.hh.

widthTable

G4PhysicsVector* G4XAnnihilationChannel::widthTable

private

Definition at line 97 of file G4XAnnihilationChannel.hh.

---

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

• G4XAnnihilationChannel.hh
• G4XAnnihilationChannel.cc