G4ParticleHPInterpolator Class Reference

#include <G4ParticleHPInterpolator.hh>

Collaboration diagram for G4ParticleHPInterpolator:

Public Member Functions
	G4ParticleHPInterpolator ()
	~G4ParticleHPInterpolator ()
G4double	Lin (G4double x, G4double x1, G4double x2, G4double y1, G4double y2)
G4double	Interpolate (G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	Interpolate2 (G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	GetBinIntegral (const G4InterpolationScheme &aScheme, const G4double x1, const G4double x2, const G4double y1, const G4double y2)
G4double	GetWeightedBinIntegral (const G4InterpolationScheme &aScheme, const G4double x1, const G4double x2, const G4double y1, const G4double y2)

Private Member Functions
G4double	Histogram (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	LinearLinear (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	LinearLogarithmic (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	LogarithmicLinear (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	LogarithmicLogarithmic (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
G4double	Random (G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const

Detailed Description

Definition at line 44 of file G4ParticleHPInterpolator.hh.

Constructor & Destructor Documentation

G4ParticleHPInterpolator()

G4ParticleHPInterpolator::G4ParticleHPInterpolator ( )

inline

Definition at line 48 of file G4ParticleHPInterpolator.hh.

{}

~G4ParticleHPInterpolator()

G4ParticleHPInterpolator::~G4ParticleHPInterpolator ( )

inline

Definition at line 49 of file G4ParticleHPInterpolator.hh.

   {
    //  G4cout <<"deleted the interpolator"<<G4endl;
   }

Member Function Documentation

GetBinIntegral()

G4double G4ParticleHPInterpolator::GetBinIntegral	(	const G4InterpolationScheme &	aScheme,
		const G4double	x1,
		const G4double	x2,
		const G4double	y1,
		const G4double	y2
	)

Definition at line 36 of file G4ParticleHPInterpolator.cc.

  { // inline again later on @@@@
    G4double result = 0;
    if(aScheme==HISTO||aScheme==CHISTO||aScheme==UHISTO)
    {
      result = y1*(x2-x1);
    }
    else if(aScheme==LINLIN||aScheme==CLINLIN||aScheme==ULINLIN)
    {
      result = 0.5*(y2+y1)*(x2-x1);
    }
    else if(aScheme==LINLOG||aScheme==CLINLOG||aScheme==ULINLOG)
    {
      if(x1==0) result = y1;
      else if(x2==0) result = y2;
      else
      {
        G4double b = (y2-y1)/(G4Log(x2)-G4Log(x1));
        G4double a = y1 - b*G4Log(x1);
        result = (a-b)*(x2-x1) + b*(x2*G4Log(x2)-x1*G4Log(x1));
      }
    }
    else if(aScheme==LOGLIN||aScheme==CLOGLIN||aScheme==ULOGLIN)
    {
      if ( y1==0||y2==0 ) { 
         result =0;
      } else {
        // G4double b = (std::log(y2)-std::log(y1))/(x2-x1);
        // G4double a = std::log(y1) - b*x1;
        //***************************************************************
        //EMendoza:
        //result = (std::exp(a)/b)*(std::exp(b*x2)-std::exp(b*x1));
        //***************************************************************
        if ( y1!=y2 ) {
           result = (x2-x1)*(y2-y1)/G4Log(y2/y1);
        } else { 
          result = y2*(x2-x1);
        }
        //***************************************************************
      }
    }
    else if(aScheme==LOGLOG||aScheme==CLOGLOG||aScheme==ULOGLOG)
    {
      if(x1==0) result = y1;
      else if(x2==0) result = y2;
      else if(y1==0||y2==0) result =0;
      else
      {      
        G4double b = (G4Log(y2)-G4Log(y1))/(G4Log(x2)-G4Log(x1));
        G4double a = G4Log(y1) - b*G4Log(x1);;
        result = (G4Exp(a)/(b+1))*(G4Pow::GetInstance()->powA(x2,b+1)-G4Pow::GetInstance()->powA(x1,b+1));
      }
    }
    else
    {
      throw G4HadronicException(__FILE__, __LINE__, "Unknown interpolation scheme in G4ParticleHPVector::Integrate");
    }
    return result;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetWeightedBinIntegral()

G4double G4ParticleHPInterpolator::GetWeightedBinIntegral	(	const G4InterpolationScheme &	aScheme,
		const G4double	x1,
		const G4double	x2,
		const G4double	y1,
		const G4double	y2
	)

Definition at line 97 of file G4ParticleHPInterpolator.cc.

  { // inline again later on @@@@
    G4double result = 0;
    if(aScheme==HISTO||aScheme==CHISTO||aScheme==UHISTO)
    {
      result = 0.5*y1*(x2*x2-x1*x1);
    }
    else if(aScheme==LINLIN||aScheme==CLINLIN||aScheme==ULINLIN)
    {
      //        G4double b = (y2-y1)/(x2-x1);
      //        G4double a = y1 - b*x1;
      //        result = 0.5*a*(x2*x2-x1*x1) + (b/3.)*(x2*x2*x2-x1*x1*x1);
      //  Factor out x2-x1 to avoid divide by zero
      
      result = (y1*x2 - y2*x1)*(x2 + x1)/2. + (y2-y1)*(x2*x2 + x2*x1 + x1*x1)/3.;
    }
    else if(aScheme==LINLOG||aScheme==CLINLOG||aScheme==ULINLOG)
    {
      if(x1==0) result = y1;
      else if(x2==0) result = y2;
      else
      {
        G4double b = (y2-y1)/(G4Log(x2)-G4Log(x1));
        G4double a = y1 - b*G4Log(x1);
        result = ( x2*x2/2. * (a-b/2.+b*G4Log(x2)) )
                -( x1*x1/2. * (a-b/2.+b*G4Log(x1)) );
      }
    }
    else if(aScheme==LOGLIN||aScheme==CLOGLIN||aScheme==ULOGLIN)
    {
      if(y1==0||y2==0) result = 0;
      else
      {
        G4double b = (G4Log(y2)-G4Log(y1))/(x2-x1);
        G4double a = G4Log(y1) - b*x1;
        result = G4Exp(a)/(b*b)*( G4Exp(b*x2)*(b*x2-1.) - G4Exp(b*x1)*(b*x1-1.) );
      }
    }
    else if(aScheme==LOGLOG||aScheme==CLOGLOG||aScheme==ULOGLOG)
    {
      if(x1==0) result = y1;
      else if(x2==0) result = y2;
      else if(y1==0||y2==0) result = 0;
      else
      {
        G4double b = (G4Log(y2)-G4Log(y1))/(G4Log(x2)-G4Log(x1));
        G4double a = G4Log(y1) - b*G4Log(x1);;
        result = G4Exp(a)/(b+2.)*( G4Pow::GetInstance()->powA(x2, b+2.) - G4Pow::GetInstance()->powA(x1, b+2) );
      }
    }
    else
    {
      throw G4HadronicException(__FILE__, __LINE__, "Unknown interpolation scheme in G4ParticleHPVector::Integrate");
    }
    return result;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

Histogram()

G4double G4ParticleHPInterpolator::Histogram ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 163 of file G4ParticleHPInterpolator.hh.

{
  G4double result;
  result = y1;
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Interpolate()

G4double G4ParticleHPInterpolator::Interpolate ( G4InterpolationScheme  aScheme, G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inline

Definition at line 91 of file G4ParticleHPInterpolator.hh.

{
  G4double result(0);
  G4int theScheme = aScheme;
  theScheme = theScheme%CSTART_;
  switch(theScheme)
  {
    case 1:
      //080809
      //result = Histogram(x, x1, x2, y1, y2);
      result = LinearLinear(x, x1, x2, y1, y2);
      break;
    case 2:
      result = LinearLinear(x, x1, x2, y1, y2);
      break;
    case 3:
      result = LinearLogarithmic(x, x1, x2, y1, y2);
      break;
    case 4:
      result = LogarithmicLinear(x, x1, x2, y1, y2);
      break;
    case 5:
      result = LogarithmicLogarithmic(x, x1, x2, y1, y2);
      break;
    case 6:
      result = Random(x, x1, x2, y1, y2);
      break;
    default:
      G4cout << "theScheme = "<<theScheme<<G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPInterpolator::Carthesian Invalid InterpolationScheme");
      break;
  }
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Interpolate2()

G4double G4ParticleHPInterpolator::Interpolate2 ( G4InterpolationScheme  aScheme, G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inline

Definition at line 128 of file G4ParticleHPInterpolator.hh.

{
  G4double result(0);
  G4int theScheme = aScheme;
  theScheme = theScheme%CSTART_;
  switch(theScheme)
  {
    case 1:
      result = Histogram(x, x1, x2, y1, y2);
      break;
    case 2:
      result = LinearLinear(x, x1, x2, y1, y2);
      break;
    case 3:
      result = LinearLogarithmic(x, x1, x2, y1, y2);
      break;
    case 4:
      result = LogarithmicLinear(x, x1, x2, y1, y2);
      break;
    case 5:
      result = LogarithmicLogarithmic(x, x1, x2, y1, y2);
      break;
    case 6:
      result = Random(x, x1, x2, y1, y2);
      break;
    default:
      G4cout << "theScheme = "<<theScheme<<G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPInterpolator::Carthesian Invalid InterpolationScheme");
      break;
  }
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Lin()

G4double G4ParticleHPInterpolator::Lin ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  )

inline

Definition at line 54 of file G4ParticleHPInterpolator.hh.

  {
    G4double slope=0, off=0;
    if(x2-x1==0) return (y2+y1)/2.;
    slope = (y2-y1)/(x2-x1);
    off = y2-x2*slope;
    G4double y = x*slope+off;
    return y;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

LinearLinear()

G4double G4ParticleHPInterpolator::LinearLinear ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 171 of file G4ParticleHPInterpolator.hh.

{
  G4double slope=0, off=0;
  if(x2-x1==0) return (y2+y1)/2.;
  slope = (y2-y1)/(x2-x1);
  off = y2-x2*slope;
  G4double y = x*slope+off;
  return y;
}

Here is the call graph for this function:

Here is the caller graph for this function:

LinearLogarithmic()

G4double G4ParticleHPInterpolator::LinearLogarithmic ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 182 of file G4ParticleHPInterpolator.hh.

{
  G4double result;
  if(x==0) result = y1+y2/2.;
  else if(x1==0) result = y1;
  else if(x2==0) result = y2;
  else result = LinearLinear(G4Log(x), G4Log(x1), G4Log(x2), y1, y2);
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

LogarithmicLinear()

G4double G4ParticleHPInterpolator::LogarithmicLinear ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 193 of file G4ParticleHPInterpolator.hh.

{
  G4double result;
  if(y1==0||y2==0) result = 0;
  else 
  {
    result = LinearLinear(x, x1, x2, G4Log(y1), G4Log(y2));
    result = G4Exp(result);
  }
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

LogarithmicLogarithmic()

G4double G4ParticleHPInterpolator::LogarithmicLogarithmic ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 206 of file G4ParticleHPInterpolator.hh.

{
  if(x==0) return y1+y2/2.;
  else if(x1==0) return y1;
  else if(x2==0) return y2;
  G4double result;
  if(y1==0||y2==0) result = 0;
  else 
  {
    result = LinearLinear(G4Log(x), G4Log(x1), G4Log(x2), G4Log(y1), G4Log(y2));
    result = G4Exp(result);
  }
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Random()

G4double G4ParticleHPInterpolator::Random ( G4double  x, G4double  x1, G4double  x2, G4double  y1, G4double  y2  ) const

inlineprivate

Definition at line 222 of file G4ParticleHPInterpolator.hh.

{
  G4double result;
  result = y1+G4UniformRand()*(y2-y1);
  return result;
}

Here is the caller graph for this function:

---

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

• G4ParticleHPInterpolator.hh
• G4ParticleHPInterpolator.cc