#include <G4LindhardSorensenData.hh>

Public Member Functions
	G4LindhardSorensenData ()

	~G4LindhardSorensenData ()

G4double	GetDeltaL (G4int Z, G4double gamma) const

G4LindhardSorensenData &	operator= (const G4LindhardSorensenData &right)=delete

	G4LindhardSorensenData (const G4LindhardSorensenData &)=delete

Detailed Description

Definition at line 52 of file G4LindhardSorensenData.hh.

Constructor & Destructor Documentation

◆ G4LindhardSorensenData() [1/2]

G4LindhardSorensenData::G4LindhardSorensenData ( )

explicit

Definition at line 56 of file G4LindhardSorensenData.cc.

{
  g4calc = G4Pow::GetInstance();
  InitialiseData();
}

◆ ~G4LindhardSorensenData()

G4LindhardSorensenData::~G4LindhardSorensenData ( )

Definition at line 64 of file G4LindhardSorensenData.cc.

{
  for(G4int i=0; i<=LVECT; ++i) { delete data[i]; }
}

◆ G4LindhardSorensenData() [2/2]

G4LindhardSorensenData::G4LindhardSorensenData ( const G4LindhardSorensenData & )

delete

Member Function Documentation

◆ GetDeltaL()

G4double G4LindhardSorensenData::GetDeltaL	(	G4int	Z,
		G4double	gamma
	)		const

Definition at line 71 of file G4LindhardSorensenData.cc.

{
  G4int idx = 0;
  for(; idx<LVECT; ++idx) {
    if(Z < zlist[idx+1]) { break; }
  }
  idx = std::min(idx, LVECT);
 
  G4double x = G4Log(gamma - 1.0);
  G4double y = ComputeDeltaL(idx, x);
 
  // interpolation over Z if needed
  if(idx < LVECT && Z > zlist[idx]) {
    G4double y1 = ComputeDeltaL(idx+1, x);
    //G4cout << "idx= " << idx << " x= " << x << " y= " << y << " y1= " << y1 << G4endl;
    y += (y1 - y)*(Z - zlist[idx])/(G4double)(zlist[idx+1] - zlist[idx]);
  }
  y *= g4calc->Z23(Z);
  return y;
}

Referenced by G4LindhardSorensenIonModel::CorrectionsAlongStep().

◆ operator=()

G4LindhardSorensenData & G4LindhardSorensenData::operator= ( const G4LindhardSorensenData & right )