G4XNNElasticLowE Class Reference

#include <G4XNNElasticLowE.hh>

Inheritance diagram for G4XNNElasticLowE:

Public Member Functions
	G4XNNElasticLowE ()
virtual	~G4XNNElasticLowE ()
G4bool	operator== (const G4XNNElasticLowE &right) const
G4bool	operator!= (const G4XNNElasticLowE &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
virtual void	Print () const
virtual G4double	HighLimit () const
Public Member Functions inherited from G4VCrossSectionSource
	G4VCrossSectionSource ()
virtual	~G4VCrossSectionSource ()
G4bool	operator== (const G4VCrossSectionSource &right) const
G4bool	operator!= (const G4VCrossSectionSource &right) const
virtual void	PrintAll (const G4KineticTrack &trk1, const G4KineticTrack &trk2) 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 39 of file G4XNNElasticLowE.hh.

Constructor & Destructor Documentation

G4XNNElasticLowE::G4XNNElasticLowE

(

)

Definition at line 88 of file G4XNNElasticLowE.cc.

References python.hepunit::GeV, python.hepunit::millibarn, G4Neutron::NeutronDefinition(), G4InuclParticleNames::pp, G4Proton::ProtonDefinition(), and G4PhysicsVector::PutValue().

{ 
  // Cross-sections are available in the range (_eMin,_eMax)

  _eMin = _eMinTable * GeV;
  _eMax = std::exp(std::log(_eMinTable) + tableSize * _eStepLog) * GeV;
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");    
  if (_highLimit > _eMax)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - High energy limit not valid");    
  G4PhysicsVector* pp = new G4PhysicsLnVector(_eMin,_eMax,tableSize);

  _eMin = std::exp(std::log(_eMinTable)-_eStepLog)*GeV;
  if (_eMin < _lowLimit)
    throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");
  G4PhysicsVector* np = new G4PhysicsLnVector(_eMin,_eMax,tableSize);

  G4int i;
  for (i=0; i<tableSize; i++)
    {
      G4double value = ppTable[i] * millibarn;
      pp->PutValue(i,value);
      value = npTable[i] * millibarn;
      np->PutValue(i,value);
    }
  xMap[G4Proton::ProtonDefinition()] = pp;
  xMap[G4Neutron::NeutronDefinition()] = np;
}

G4XNNElasticLowE::~G4XNNElasticLowE ( )

virtual

Definition at line 118 of file G4XNNElasticLowE.cc.

References G4Neutron::NeutronDefinition(), and G4Proton::ProtonDefinition().

{
  delete xMap[G4Proton::ProtonDefinition()];
  delete xMap[G4Neutron::NeutronDefinition()];
}

Member Function Documentation

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

virtual

Implements G4VCrossSectionSource.

Definition at line 138 of file G4XNNElasticLowE.cc.

References G4VCrossSectionSource::FindKeyParticle(), G4KineticTrack::Get4Momentum(), and G4PhysicsVector::GetValue().

{
  G4double sigma = 0.;
  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
  G4bool dummy = false;

  G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);

  typedef std::map <G4ParticleDefinition *, G4PhysicsVector*, std::less<G4ParticleDefinition *> > StringPhysMap;

  if (xMap.find(key)!= xMap.end())
    {

      StringPhysMap::const_iterator iter;
      for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      G4ParticleDefinition * str = (*iter).first;
          if (str == key)
        {
          G4PhysicsVector* physVector = (*iter).second; 
          //     G4PhysicsVector* physVector = xMap[key];
          if (sqrtS >= _eMin && sqrtS <= _eMax)
        {
          sigma = physVector->GetValue(sqrtS,dummy);
        } else if ( sqrtS < _eMin )
                {
                  sigma = physVector->GetValue(_eMin,dummy);
                }
        //G4cout << " sqrtS / sigma " << sqrtS/GeV << " / " <<
        //          sigma/millibarn << G4endl;
        }
    }
    }
  return sigma;
}

virtual const G4CrossSectionVector* G4XNNElasticLowE::GetComponents ( ) const

inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 53 of file G4XNNElasticLowE.hh.

{ return 0; }

virtual G4double G4XNNElasticLowE::HighLimit ( ) const

inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 61 of file G4XNNElasticLowE.hh.

{ return _highLimit; }

G4bool G4XNNElasticLowE::IsValid ( G4double  e ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 246 of file G4XNNElasticLowE.cc.

References G4VCrossSectionSource::InLimits().

{
  G4bool answer = InLimits(e,_lowLimit,_highLimit);

  return answer;
}

G4String G4XNNElasticLowE::Name ( ) const

virtual

Implements G4VCrossSectionSource.

Definition at line 238 of file G4XNNElasticLowE.cc.

Referenced by Print().

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

G4bool G4XNNElasticLowE::operator!=

(

const G4XNNElasticLowE &

right

)

const

Definition at line 131 of file G4XNNElasticLowE.cc.

{

  return (this != (G4XNNElasticLowE *) &right);
}

G4bool G4XNNElasticLowE::operator==

(

const G4XNNElasticLowE &

right

)

const

Definition at line 125 of file G4XNNElasticLowE.cc.

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

void G4XNNElasticLowE::Print ( void  ) const

virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 175 of file G4XNNElasticLowE.cc.

References G4cout, G4endl, G4PhysicsVector::GetLowEdgeEnergy(), G4PhysicsVector::GetValue(), python.hepunit::GeV, python.hepunit::millibarn, Name(), G4Neutron::NeutronDefinition(), G4InuclParticleNames::pp, G4VCrossSectionSource::Print(), and G4Proton::ProtonDefinition().

{
  // Dump the pp cross-section table

  G4cout << Name() << ", pp cross-section: " << G4endl;

  G4bool dummy = false;
  G4int i;
  G4ParticleDefinition * key = G4Proton::ProtonDefinition();
  G4PhysicsVector* pp = 0;

  typedef std::map <G4ParticleDefinition *, G4PhysicsVector*, std::less<G4ParticleDefinition *> > StringPhysMap;
  StringPhysMap::const_iterator iter;

  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      G4ParticleDefinition * str = (*iter).first;
      if (str == key)
    {
      pp = (*iter).second; 
    }
    }
  
  if (pp != 0)
    {   
      for (i=0; i<tableSize; i++)
    {
      G4double e = pp->GetLowEdgeEnergy(i);
      G4double sigma = pp->GetValue(e,dummy) / millibarn;
      G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
    }
    }
  
  // Dump the np cross-section table

  G4cout << Name() << ", np cross-section: " << G4endl;

  key = G4Neutron::NeutronDefinition();
  G4PhysicsVector* np = 0;
  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
    {
      G4ParticleDefinition * str = (*iter).first;
      if (str == key)
    {
      np = (*iter).second; 
    }
    }
  
  //  G4PhysicsVector* np = xMap[G4Neutron::NeutronDefinition()->GetParticleName()];
  
  if (np != 0)
    {   
      for (i=0; i<tableSize; i++)
    {
      G4double e = np->GetLowEdgeEnergy(i);
      G4double sigma = np->GetValue(e,dummy) / millibarn;
      G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
    }
    }
  G4VCrossSectionSource::Print();
}

---

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

• G4XNNElasticLowE.hh
• G4XNNElasticLowE.cc