#include <G4LogLogInterpolator.hh>

Inheritance diagram for G4LogLogInterpolator:

Public Member Functions
G4double	Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data) const

virtual G4IInterpolator *	Clone () const

	G4LogLogInterpolator ()

	~G4LogLogInterpolator ()

Detailed Description

Definition at line 49 of file G4LogLogInterpolator.hh.

Constructor & Destructor Documentation

◆ G4LogLogInterpolator()

G4LogLogInterpolator::G4LogLogInterpolator ( )

Definition at line 41 of file G4LogLogInterpolator.cc.

42{ }

Referenced by Clone().

◆ ~G4LogLogInterpolator()

G4LogLogInterpolator::~G4LogLogInterpolator ( )

Definition at line 46 of file G4LogLogInterpolator.cc.

47{ }

Member Function Documentation

◆ Calculate()

G4double G4LogLogInterpolator::Calculate	(	G4double	point,
		G4int	bin,
		const G4DataVector &	energies,
		const G4DataVector &	data
	)		const

virtual

Implements G4IInterpolator.

Definition at line 53 of file G4LogLogInterpolator.cc.

{
  G4int nBins = data.size() - 1;
  // ---- MGP ---- To be verified; value = 0. is redundant here
  G4double value = 0.;
  if (x < points[0] || x == 0.)
    {
      value = 0;
    }
  else if (bin < nBins)
    {
      G4double e1 = points[bin];
      G4double e2 = points[bin+1];
      G4double d1 = data[bin];     
      G4double d2 = data[bin+1];
      
      if (d1 != 0. && d2 != 0.)
    {
      value = (std::log10(d1)*std::log10(e2/x) + std::log10(d2)*std::log10(x/e1)) / std::log10(e2/e1);
      value = std::pow(10.,value);
    }
      else
    {
      value = 0.;
    }
    }
  else
    {
      value = data[nBins];
    }
 
  return value;
}