#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.

42{ }

◆ ~G4VCrossSectionSource()

G4VCrossSectionSource::~G4VCrossSectionSource ( )

virtual

Definition at line 45 of file G4VCrossSectionSource.cc.

46{ }

Member Function Documentation

◆ CrossSection()

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

pure virtual

◆ FcrossX()

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

protected

Definition at line 172 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 156 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 G4CrossSectionComposite, G4CrossSectionPatch, G4XAnnihilationChannel, G4XAqmElastic, G4XAqmTotal, G4XMesonBaryonElastic, G4XNNElastic, G4XNNElasticLowE, G4XNNTotal, G4XNNTotalLowE, G4XnpElastic, G4XnpElasticLowE, G4XnpTotal, G4XnpTotalLowE, G4XPDGElastic, G4XPDGTotal, G4XpimNTotal, G4XpipNTotal, and G4XResonance.

Referenced by Print(), and PrintAll().

◆ HighLimit()

G4double G4VCrossSectionSource::HighLimit ( ) const

virtual

Reimplemented in G4XNNElasticLowE, G4XNNTotalLowE, G4XnpElasticLowE, and G4XnpTotalLowE.

Definition at line 144 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 131 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 149 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 138 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.

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

Referenced by G4VCollision::Print(), 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 106 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;
 
  std::size_t nComponents = 0;
  const G4CrossSectionVector* components = GetComponents();
  if (components != 0)
    {
      nComponents = components->size();
    }
  for (std::size_t 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:

Public Member Functions

Protected Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ G4VCrossSectionSource()

◆ ~G4VCrossSectionSource()

Member Function Documentation

◆ CrossSection()

◆ FcrossX()

◆ FindKeyParticle()

◆ FindLightParticle()

◆ GetComponents()

◆ HighLimit()

◆ InLimits()

◆ IsValid()

◆ LowLimit()

◆ Name()

◆ operator!=()

◆ operator==()

◆ Print()

◆ PrintAll()