G4VCrossSectionSource Class Reference

#include <G4VCrossSectionSource.hh>

Inheritance diagram for G4VCrossSectionSource:

Public Member Functions
	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

Protected Member Functions
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 38 of file G4VCrossSectionSource.hh.

Constructor & Destructor Documentation

G4VCrossSectionSource()

G4VCrossSectionSource::G4VCrossSectionSource

(

)

Definition at line 41 of file G4VCrossSectionSource.cc.

{ }

~G4VCrossSectionSource()

G4VCrossSectionSource::~G4VCrossSectionSource ( )

virtual

Definition at line 45 of file G4VCrossSectionSource.cc.

{ }

Member Function Documentation

CrossSection()

virtual G4double G4VCrossSectionSource::CrossSection ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

pure virtual

Implemented in G4CrossSectionComposite, G4CrossSectionPatch, G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, G4XnpTotalLowE, G4XPDGElastic, G4XPDGTotal, G4XpimNTotal, G4XpipNTotal, and G4XResonance.

Referenced by G4CollisionComposite::CrossSection(), G4CollisionNN::CrossSection(), G4CrossSectionComposite::CrossSection(), G4CrossSectionPatch::CrossSection(), G4VCollision::CrossSection(), PrintAll(), and G4CrossSectionPatch::Transition().

FcrossX()

G4double G4VCrossSectionSource::FcrossX ( G4double  e, G4double  e0, G4double  sigma, G4double  eParam, G4double  power  ) const

protected

Definition at line 174 of file G4VCrossSectionSource.cc.

{
  G4double result = 0.;
 
  G4double denom = eParam*eParam + (e-e0)*(e-e0);
  if (denom > 0.) 
  {
    G4double value = (2.* eParam * sigma * (e-e0) / denom) * G4Pow::GetInstance()->powA(((e0 + eParam) / e), power);
    result = std::max(0., value);
  }
  return result;
}     

FindKeyParticle()

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

protected

Definition at line 49 of file G4VCrossSectionSource.cc.

{
  const G4ParticleDefinition * result;
  
  const G4ParticleDefinition * p1 = trk1.GetDefinition();
  const G4ParticleDefinition * p2 = trk2.GetDefinition();
  
  if( (p1==G4Proton::Proton() && p2==G4Proton::Proton() ) ||
      (p1==G4Neutron::Neutron() && p2==G4Neutron::Neutron()) )
  {
    result = G4Proton::Proton();
  }
  else if( (p1==G4Neutron::Neutron() && p2==G4Proton::Proton()) ||
           (p2==G4Neutron::Neutron() && p1==G4Proton::Proton()) )
  {
    result = G4Neutron::Neutron();
  }
  else
  {
    throw G4HadronicException(__FILE__, __LINE__, "G4VCrossSectionSource: unklnown particles in FindKeyParticle");
  }
  return result;
}

Referenced by G4XNNElasticLowE::CrossSection(), and G4XNNTotalLowE::CrossSection().

FindLightParticle()

const G4ParticleDefinition * G4VCrossSectionSource::FindLightParticle ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

protected

Definition at line 158 of file G4VCrossSectionSource.cc.

{
  G4double mass1 = trk1.GetDefinition()->GetPDGMass();
  G4double mass2 = trk2.GetDefinition()->GetPDGMass();
  if (mass1 < mass2)
    {
      return trk1.GetDefinition();
    }
  else
    {
      return trk2.GetDefinition();
    }
}

Referenced by G4XMesonBaryonElastic::CrossSection().

GetComponents()

virtual const G4CrossSectionVector * G4VCrossSectionSource::GetComponents ( ) const

pure virtual

Implemented in G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElastic, G4XNNElasticLowE, G4XNNTotal, G4XNNTotalLowE, G4XnpElastic, G4XnpElasticLowE, G4XnpTotal, G4XnpTotalLowE, G4XPDGElastic, G4XPDGTotal, G4XpimNTotal, G4XpipNTotal, G4XResonance, G4CrossSectionComposite, and G4CrossSectionPatch.

Referenced by Print(), and PrintAll().

HighLimit()

G4double G4VCrossSectionSource::HighLimit ( ) const

virtual

Reimplemented in G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 146 of file G4VCrossSectionSource.cc.

{
  return DBL_MAX; 
}

Referenced by G4CrossSectionPatch::CrossSection(), IsValid(), and G4CrossSectionPatch::Transition().

InLimits()

G4bool G4VCrossSectionSource::InLimits ( G4double  e, G4double  eLow, G4double  eHigh  ) const

protected

Definition at line 133 of file G4VCrossSectionSource.cc.

{
  G4bool answer = false;
  if (e >= eLow && e <= eHigh) answer = true;
  return answer;
}

Referenced by G4XAnnihilationChannel::IsValid(), G4XAqmElastic::IsValid(), G4XAqmTotal::IsValid(), G4XMesonBaryonElastic::IsValid(), G4XNNElasticLowE::IsValid(), G4XnpElasticLowE::IsValid(), G4XnpTotalLowE::IsValid(), G4XPDGElastic::IsValid(), and G4XPDGTotal::IsValid().

IsValid()

G4bool G4VCrossSectionSource::IsValid ( G4double  e ) const

virtual

Reimplemented in G4CrossSectionComposite, G4CrossSectionPatch, G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, G4XnpTotalLowE, G4XPDGElastic, and G4XPDGTotal.

Definition at line 151 of file G4VCrossSectionSource.cc.

{
  G4bool answer = false;
  if (e >= LowLimit() && e <= HighLimit()) answer = true;
  return answer;
}

Referenced by G4CrossSectionComposite::CrossSection(), G4CrossSectionPatch::CrossSection(), G4CrossSectionComposite::IsValid(), and G4CrossSectionPatch::IsValid().

LowLimit()

G4double G4VCrossSectionSource::LowLimit ( ) const

virtual

Reimplemented in G4XPDGElastic, and G4XPDGTotal.

Definition at line 140 of file G4VCrossSectionSource.cc.

{
  return 0.; 
}

Referenced by G4CrossSectionPatch::CrossSection(), IsValid(), and G4CrossSectionPatch::Transition().

Name()

virtual G4String G4VCrossSectionSource::Name ( ) const

pure virtual

Implemented in G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElastic, G4XNNElasticLowE, G4XNNTotal, G4XNNTotalLowE, G4XnpElastic, G4XnpElasticLowE, G4XnpTotal, G4XnpTotalLowE, G4XPDGElastic, G4XPDGTotal, G4XpimNTotal, G4XpipNTotal, and G4XResonance.

Referenced by Print(), and PrintAll().

operator!=()

G4bool G4VCrossSectionSource::operator!=

(

const G4VCrossSectionSource &

right

)

const

Definition at line 80 of file G4VCrossSectionSource.cc.

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

operator==()

G4bool G4VCrossSectionSource::operator==

(

const G4VCrossSectionSource &

right

)

const

Definition at line 74 of file G4VCrossSectionSource.cc.

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

Print()

void G4VCrossSectionSource::Print ( ) const

virtual

Reimplemented in G4XNNElasticLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 86 of file G4VCrossSectionSource.cc.

{
  G4int nComponents = 0;
  const G4CrossSectionVector* components = GetComponents();
  if (components)
    {
      nComponents = components->size();
    }
  G4cout << "---- " << this->Name() << " ---- has " << nComponents << " components" <<G4endl;
  G4int i;
  for (i=0; i<nComponents; i++)
    {
      G4cout << "-" <<  this->Name() << " - Component " << i << ": " <<G4endl;
 
      G4CrossSectionSourcePtr componentPtr = (*components)[i];
      G4VCrossSectionSource* component = componentPtr();
      component->Print();
    }
}

Referenced by G4VCollision::Print(), Print(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), and G4XnpTotalLowE::Print().

PrintAll()

void G4VCrossSectionSource::PrintAll ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

virtual

Definition at line 107 of file G4VCrossSectionSource.cc.

{
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4double sigma = CrossSection(trk1,trk2) / millibarn;
  G4cout << "---- " << Name() << ": "
     << "Ecm = " << sqrtS / GeV << " GeV -  " 
     << " Cross section = " << sigma << " mb "
     << G4endl;
 
  G4int nComponents = 0;
  const G4CrossSectionVector* components = GetComponents();
  if (components != 0)
    {
      nComponents = components->size();
    }
  G4int i;
  for (i=0; i<nComponents; i++)
    {
      G4cout << "* Component " << i << ": ";
      G4CrossSectionSourcePtr componentPtr = (*components)[i];
      G4VCrossSectionSource* component = componentPtr();
      component->PrintAll(trk1,trk2);
    }
}

Referenced by PrintAll().

---

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

• G4VCrossSectionSource.hh
• G4VCrossSectionSource.cc