Geant4  10.02.p03
Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
G4XNNElasticLowE Class Reference

#include <G4XNNElasticLowE.hh>

Inheritance diagram for G4XNNElasticLowE:
Collaboration 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
 

Private Member Functions

 G4XNNElasticLowE (const G4XNNElasticLowE &right)
 
const G4XNNElasticLowE & operator= (const G4XNNElasticLowE &right)
 

Private Attributes

std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > > xMap
 
G4double _eMin
 
G4double _eMax
 

Static Private Attributes

static const G4double _lowLimit = 0.
 
static const G4double _highLimit = 3.*GeV
 
static const G4double ppTable [101]
 
static const G4double npTable [101]
 
static const G4int tableSize = 101
 
static const G4double _eMinTable = 1.8964808
 
static const G4double _eStepLog = 0.01
 

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 39 of file G4XNNElasticLowE.hh.

Constructor & Destructor Documentation

◆ G4XNNElasticLowE() [1/2]

G4XNNElasticLowE::G4XNNElasticLowE ( )

Definition at line 89 of file G4XNNElasticLowE.cc.

90 {
91  // Cross-sections are available in the range (_eMin,_eMax)
92 
93  _eMin = _eMinTable * GeV;
94  _eMax = G4Exp(G4Log(_eMinTable) + tableSize * _eStepLog) * GeV;
95  if (_eMin < _lowLimit)
96  throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");
97  if (_highLimit > _eMax)
98  throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - High energy limit not valid");
99  G4PhysicsVector* pp = new G4PhysicsLnVector(_eMin,_eMax,tableSize);
100 
101  _eMin = G4Exp(G4Log(_eMinTable)-_eStepLog)*GeV;
102  if (_eMin < _lowLimit)
103  throw G4HadronicException(__FILE__, __LINE__, "G4XNNElasticLowE::G4XNNElasticLowE - Low energy limit not valid");
104  G4PhysicsVector* np = new G4PhysicsLnVector(_eMin,_eMax,tableSize);
105 
106  G4int i;
107  for (i=0; i<tableSize; i++)
108  {
109  G4double value = ppTable[i] * millibarn;
110  pp->PutValue(i,value);
111  value = npTable[i] * millibarn;
112  np->PutValue(i,value);
113  }
114  xMap[G4Proton::ProtonDefinition()] = pp;
115  xMap[G4Neutron::NeutronDefinition()] = np;
116 }
G4XNNElasticLowE::xMap
std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > > xMap
Definition: G4XNNElasticLowE.hh:80
G4XNNElasticLowE::_highLimit
static const G4double _highLimit
Definition: G4XNNElasticLowE.hh:73
G4XNNElasticLowE::_eMinTable
static const G4double _eMinTable
Definition: G4XNNElasticLowE.hh:77
G4XNNElasticLowE::tableSize
static const G4int tableSize
Definition: G4XNNElasticLowE.hh:76
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4XNNElasticLowE::npTable
static const G4double npTable[101]
Definition: G4XNNElasticLowE.hh:75
G4XNNElasticLowE::_lowLimit
static const G4double _lowLimit
Definition: G4XNNElasticLowE.hh:72
G4XNNElasticLowE::_eStepLog
static const G4double _eStepLog
Definition: G4XNNElasticLowE.hh:78
G4int
int G4int
Definition: G4Types.hh:78
G4XNNElasticLowE::ppTable
static const G4double ppTable[101]
Definition: G4XNNElasticLowE.hh:74
G4PhysicsVector::PutValue
void PutValue(size_t index, G4double theValue)
GeV
static const double GeV
Definition: G4SIunits.hh:214
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:230
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
millibarn
static const double millibarn
Definition: G4SIunits.hh:105
G4double
double G4double
Definition: G4Types.hh:76
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ~G4XNNElasticLowE()

G4XNNElasticLowE::~G4XNNElasticLowE ( )
virtual

Definition at line 119 of file G4XNNElasticLowE.cc.

120 {
121  delete xMap[G4Proton::ProtonDefinition()];
122  delete xMap[G4Neutron::NeutronDefinition()];
123 }
G4XNNElasticLowE::xMap
std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > > xMap
Definition: G4XNNElasticLowE.hh:80
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
Here is the call graph for this function:

◆ G4XNNElasticLowE() [2/2]

G4XNNElasticLowE::G4XNNElasticLowE ( const G4XNNElasticLowE &  right)
private

Member Function Documentation

◆ CrossSection()

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

Implements G4VCrossSectionSource.

Definition at line 139 of file G4XNNElasticLowE.cc.

140 {
141  G4double sigma = 0.;
142  G4double sqrtS = (trk1.Get4Momentum() + trk2.Get4Momentum()).mag();
143  G4bool dummy = false;
144 
145  const G4ParticleDefinition * key = FindKeyParticle(trk1,trk2);
146 
147  typedef std::map <const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > StringPhysMap;
148 
149  if (xMap.find(key)!= xMap.end())
150  {
151 
152  StringPhysMap::const_iterator iter;
153  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
154  {
155  const G4ParticleDefinition * str = (*iter).first;
156  if (str == key)
157  {
158  G4PhysicsVector* physVector = (*iter).second;
159  // G4PhysicsVector* physVector = xMap[key];
160  if (sqrtS >= _eMin && sqrtS <= _eMax)
161  {
162  sigma = physVector->GetValue(sqrtS,dummy);
163  } else if ( sqrtS < _eMin )
164  {
165  sigma = physVector->GetValue(_eMin,dummy);
166  }
167  //G4cout << " sqrtS / sigma " << sqrtS/GeV << " / " <<
168  // sigma/millibarn << G4endl;
169  }
170  }
171  }
172  return sigma;
173 }
G4XNNElasticLowE::xMap
std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > > xMap
Definition: G4XNNElasticLowE.hh:80
G4VCrossSectionSource::FindKeyParticle
const G4ParticleDefinition * FindKeyParticle(const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
Definition: G4VCrossSectionSource.cc:51
G4bool
bool G4bool
Definition: G4Types.hh:79
G4PhysicsVector::GetValue
G4double GetValue(G4double theEnergy, G4bool &isOutRange) const
G4double
double G4double
Definition: G4Types.hh:76
G4KineticTrack::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4KineticTrack.hh:240
Here is the call graph for this function:

◆ GetComponents()

virtual const G4CrossSectionVector* G4XNNElasticLowE::GetComponents ( ) const
inlinevirtual

Implements G4VCrossSectionSource.

Definition at line 53 of file G4XNNElasticLowE.hh.

53 { return 0; }
Here is the call graph for this function:

◆ HighLimit()

virtual G4double G4XNNElasticLowE::HighLimit ( ) const
inlinevirtual

Reimplemented from G4VCrossSectionSource.

Definition at line 61 of file G4XNNElasticLowE.hh.

61 { return _highLimit; }
G4XNNElasticLowE::_highLimit
static const G4double _highLimit
Definition: G4XNNElasticLowE.hh:73
Here is the call graph for this function:

◆ IsValid()

G4bool G4XNNElasticLowE::IsValid ( G4double  e) const
virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 247 of file G4XNNElasticLowE.cc.

248 {
249  G4bool answer = InLimits(e,_lowLimit,_highLimit);
250 
251  return answer;
252 }
G4XNNElasticLowE::_highLimit
static const G4double _highLimit
Definition: G4XNNElasticLowE.hh:73
G4XNNElasticLowE::_lowLimit
static const G4double _lowLimit
Definition: G4XNNElasticLowE.hh:72
G4bool
bool G4bool
Definition: G4Types.hh:79
G4VCrossSectionSource::InLimits
G4bool InLimits(G4double e, G4double eLow, G4double eHigh) const
Definition: G4VCrossSectionSource.cc:134
Here is the call graph for this function:
Here is the caller graph for this function:

◆ Name()

G4String G4XNNElasticLowE::Name ( ) const
virtual

Implements G4VCrossSectionSource.

Definition at line 239 of file G4XNNElasticLowE.cc.

240 {
241  G4String name("NNElasticLowE");
242  return name;
243 }
name
G4String name
Definition: TRTMaterials.hh:40
Here is the caller graph for this function:

◆ operator!=()

G4bool G4XNNElasticLowE::operator!= ( const G4XNNElasticLowE &  right) const

Definition at line 132 of file G4XNNElasticLowE.cc.

133 {
134 
135  return (this != (G4XNNElasticLowE *) &right);
136 }

◆ operator=()

const G4XNNElasticLowE& G4XNNElasticLowE::operator= ( const G4XNNElasticLowE &  right)
private
Here is the caller graph for this function:

◆ operator==()

G4bool G4XNNElasticLowE::operator== ( const G4XNNElasticLowE &  right) const

Definition at line 126 of file G4XNNElasticLowE.cc.

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

◆ Print()

void G4XNNElasticLowE::Print ( void  ) const
virtual

Reimplemented from G4VCrossSectionSource.

Definition at line 176 of file G4XNNElasticLowE.cc.

177 {
178  // Dump the pp cross-section table
179 
180  G4cout << Name() << ", pp cross-section: " << G4endl;
181 
182  G4bool dummy = false;
183  G4int i;
184  const G4ParticleDefinition * key = G4Proton::ProtonDefinition();
185  G4PhysicsVector* pp = 0;
186 
187  typedef std::map <const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > StringPhysMap;
188  StringPhysMap::const_iterator iter;
189 
190  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
191  {
192  const G4ParticleDefinition * str = (*iter).first;
193  if (str == key)
194  {
195  pp = (*iter).second;
196  }
197  }
198 
199  if (pp != 0)
200  {
201  for (i=0; i<tableSize; i++)
202  {
203  G4double e = pp->GetLowEdgeEnergy(i);
204  G4double sigma = pp->GetValue(e,dummy) / millibarn;
205  G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
206  }
207  }
208 
209  // Dump the np cross-section table
210 
211  G4cout << Name() << ", np cross-section: " << G4endl;
212 
213  key = G4Neutron::NeutronDefinition();
214  G4PhysicsVector* np = 0;
215  for (iter = xMap.begin(); iter != xMap.end(); ++iter)
216  {
217  const G4ParticleDefinition * str = (*iter).first;
218  if (str == key)
219  {
220  np = (*iter).second;
221  }
222  }
223 
224  // G4PhysicsVector* np = xMap[G4Neutron::NeutronDefinition()->GetParticleName()];
225 
226  if (np != 0)
227  {
228  for (i=0; i<tableSize; i++)
229  {
230  G4double e = np->GetLowEdgeEnergy(i);
231  G4double sigma = np->GetValue(e,dummy) / millibarn;
232  G4cout << i << ") e = " << e / GeV << " GeV ---- Cross section = " << sigma << " mb " << G4endl;
233  }
234  }
235  G4VCrossSectionSource::Print();
236 }
G4XNNElasticLowE::xMap
std::map< const G4ParticleDefinition *, G4PhysicsVector *, std::less< const G4ParticleDefinition * > > xMap
Definition: G4XNNElasticLowE.hh:80
G4VCrossSectionSource::Print
virtual void Print() const
Definition: G4VCrossSectionSource.cc:87
G4XNNElasticLowE::tableSize
static const G4int tableSize
Definition: G4XNNElasticLowE.hh:76
G4XNNElasticLowE::Name
virtual G4String Name() const
Definition: G4XNNElasticLowE.cc:239
G4Proton::ProtonDefinition
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
G4int
int G4int
Definition: G4Types.hh:78
G4PhysicsVector::GetLowEdgeEnergy
G4double GetLowEdgeEnergy(size_t binNumber) const
Definition: G4PhysicsVector.cc:153
G4cout
G4GLOB_DLL std::ostream G4cout
G4bool
bool G4bool
Definition: G4Types.hh:79
GeV
static const double GeV
Definition: G4SIunits.hh:214
G4PhysicsVector::GetValue
G4double GetValue(G4double theEnergy, G4bool &isOutRange) const
millibarn
static const double millibarn
Definition: G4SIunits.hh:105
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4Neutron::NeutronDefinition
static G4Neutron * NeutronDefinition()
Definition: G4Neutron.cc:99
Here is the call graph for this function:
Here is the caller graph for this function:

Member Data Documentation

◆ _eMax

G4double G4XNNElasticLowE::_eMax
private

Definition at line 83 of file G4XNNElasticLowE.hh.

◆ _eMin

G4double G4XNNElasticLowE::_eMin
private

Definition at line 82 of file G4XNNElasticLowE.hh.

◆ _eMinTable

const G4double G4XNNElasticLowE::_eMinTable = 1.8964808
staticprivate

Definition at line 77 of file G4XNNElasticLowE.hh.

◆ _eStepLog

const G4double G4XNNElasticLowE::_eStepLog = 0.01
staticprivate

Definition at line 78 of file G4XNNElasticLowE.hh.

◆ _highLimit

const G4double G4XNNElasticLowE::_highLimit = 3.*GeV
staticprivate

Definition at line 73 of file G4XNNElasticLowE.hh.

◆ _lowLimit

const G4double G4XNNElasticLowE::_lowLimit = 0.
staticprivate

Definition at line 72 of file G4XNNElasticLowE.hh.

◆ npTable

const G4double G4XNNElasticLowE::npTable
staticprivate
Initial value:
=
{
1500.00,
248.20, 93.38, 55.26, 44.50, 41.33, 38.48, 37.20, 35.98,
35.02, 34.47, 32.48, 30.76, 29.46, 28.53, 27.84, 27.20,
26.53, 25.95, 25.59, 25.46, 25.00, 24.49, 24.08, 23.86,
23.17, 22.70, 21.88, 21.48, 20.22, 19.75, 18.97, 18.39,
17.98, 17.63, 17.21, 16.72, 16.68, 16.58, 16.42, 16.22,
15.98, 15.71, 15.42, 15.14, 14.87, 14.65, 14.44, 14.26,
14.10, 13.95, 13.80, 13.64, 13.47, 13.29, 13.09, 12.89,
12.68, 12.47, 12.27, 12.06, 11.84, 11.76, 11.69, 11.60,
11.50, 11.41, 11.29, 11.17, 11.06, 10.93, 10.81, 10.68,
10.56, 10.44, 10.33, 10.21, 10.12, 10.03, 9.96, 9.89,
9.83, 9.80, 9.77, 9.75, 9.74, 9.74, 9.74, 9.76,
9.73, 9.70, 9.68, 9.65, 9.63, 9.60, 9.57, 9.55,
9.52, 9.49, 9.46, 9.43
}

Definition at line 75 of file G4XNNElasticLowE.hh.

◆ ppTable

const G4double G4XNNElasticLowE::ppTable
staticprivate
Initial value:
=
{
60.00,
33.48, 26.76, 25.26, 24.55, 23.94, 23.77, 23.72, 23.98,
25.48, 27.52, 27.72, 27.21, 25.80, 26.00, 24.32, 23.81,
24.37, 24.36, 23.13, 22.43, 21.71, 21.01, 20.83, 20.74,
20.25, 20.10, 20.59, 20.04, 20.83, 20.84, 21.07, 20.83,
20.79, 21.88, 21.15, 20.92, 19.00, 18.60, 17.30, 17.00,
16.70, 16.50, 16.20, 15.80, 15.57, 15.20, 15.00, 14.60,
14.20, 14.00, 13.80, 13.60, 13.40, 13.20, 13.00, 12.85,
12.70, 12.60, 12.50, 12.40, 12.30, 12.20, 12.10, 12.00,
11.90, 11.80, 11.75, 11.70, 11.64, 11.53, 11.41, 11.31,
11.22, 11.13, 11.05, 10.97, 10.89, 10.82, 10.75, 10.68,
10.61, 10.54, 10.48, 10.41, 10.35, 10.28, 10.22, 10.16,
10.13, 10.10, 10.08, 10.05, 10.02, 9.99, 9.96, 9.93,
9.90, 9.87, 9.84, 9.80
}

Definition at line 74 of file G4XNNElasticLowE.hh.

◆ tableSize

const G4int G4XNNElasticLowE::tableSize = 101
staticprivate

Definition at line 76 of file G4XNNElasticLowE.hh.

◆ xMap

std::map<const G4ParticleDefinition *, G4PhysicsVector*, std::less<const G4ParticleDefinition *> > G4XNNElasticLowE::xMap
private

Definition at line 80 of file G4XNNElasticLowE.hh.

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