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
G4ParticleDefinition *	FindKeyParticle (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4double	GetTransversePionMass () const
G4double	GetMinStringMass () const

Detailed Description

Definition at line 38 of file G4VCrossSectionSource.hh.

Constructor & Destructor Documentation

G4VCrossSectionSource::G4VCrossSectionSource

(

)

Definition at line 41 of file G4VCrossSectionSource.cc.

{ }

G4VCrossSectionSource::~G4VCrossSectionSource ( )

virtual

Definition at line 45 of file G4VCrossSectionSource.cc.

{ }

Member Function Documentation

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

pure virtual

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

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

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

protected

Definition at line 174 of file G4VCrossSectionSource.cc.

References G4INCL::Math::max().

{
  G4double result = 0.;

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

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

protected

Definition at line 50 of file G4VCrossSectionSource.cc.

References G4KineticTrack::GetDefinition(), G4Neutron::Neutron(), and G4Proton::Proton().

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

{
  G4ParticleDefinition * result;
  
  G4ParticleDefinition * p1 = trk1.GetDefinition();
  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;
}

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

protected

Definition at line 158 of file G4VCrossSectionSource.cc.

References G4KineticTrack::GetDefinition(), and G4ParticleDefinition::GetPDGMass().

Referenced by G4XMesonBaryonElastic::CrossSection().

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

virtual const G4CrossSectionVector* G4VCrossSectionSource::GetComponents ( ) const

pure virtual

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

Referenced by Print(), and PrintAll().

G4double G4VCrossSectionSource::GetMinStringMass ( ) const

protected

Definition at line 196 of file G4VCrossSectionSource.cc.

References python.hepunit::GeV.

{
  // Parameter from UrQMD
  const G4double minStringMass = 0.52 * GeV;
  return minStringMass;
}

G4double G4VCrossSectionSource::GetTransversePionMass ( ) const

protected

Definition at line 188 of file G4VCrossSectionSource.cc.

References python.hepunit::GeV.

{
  // Parameter from UrQMD
  const G4double transversePionMass = 0.3 * GeV;
  return transversePionMass;
}

G4double G4VCrossSectionSource::HighLimit ( ) const

virtual

Reimplemented in G4XnpElasticLowE, G4XnpTotalLowE, G4XNNElasticLowE, and G4XNNTotalLowE.

Definition at line 146 of file G4VCrossSectionSource.cc.

References DBL_MAX.

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

{
  return DBL_MAX; 
}

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

protected

Definition at line 133 of file G4VCrossSectionSource.cc.

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

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

G4bool G4VCrossSectionSource::IsValid ( G4double  e ) const

virtual

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

Definition at line 151 of file G4VCrossSectionSource.cc.

References HighLimit(), and LowLimit().

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

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

G4double G4VCrossSectionSource::LowLimit ( ) const

virtual

Reimplemented in G4XPDGElastic, and G4XPDGTotal.

Definition at line 140 of file G4VCrossSectionSource.cc.

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

{
  return 0.; 
}

virtual G4String G4VCrossSectionSource::Name ( ) const

pure virtual

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

Referenced by Print(), and PrintAll().

G4bool G4VCrossSectionSource::operator!=

(

const G4VCrossSectionSource &

right

)

const

Definition at line 80 of file G4VCrossSectionSource.cc.

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

G4bool G4VCrossSectionSource::operator==

(

const G4VCrossSectionSource &

right

)

const

Definition at line 74 of file G4VCrossSectionSource.cc.

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

void G4VCrossSectionSource::Print ( void  ) const

virtual

Reimplemented in G4XnpElasticLowE, G4XnpTotalLowE, and G4XNNElasticLowE.

Definition at line 86 of file G4VCrossSectionSource.cc.

References G4cout, G4endl, GetComponents(), Name(), nComponents, and Print().

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

{
  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();
    }
}

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

virtual

Definition at line 107 of file G4VCrossSectionSource.cc.

References CrossSection(), G4cout, G4endl, G4KineticTrack::Get4Momentum(), GetComponents(), python.hepunit::GeV, python.hepunit::millibarn, Name(), nComponents, and PrintAll().

Referenced by PrintAll().

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

---

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

• G4VCrossSectionSource.hh
• G4VCrossSectionSource.cc