G4ParticleHPInterpolator Class Reference

#include <G4ParticleHPInterpolator.hh>

Public Member Functions
	G4ParticleHPInterpolator ()
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)
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	Lin (G4double x, G4double x1, G4double x2, G4double y1, G4double y2)
	~G4ParticleHPInterpolator ()

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 35 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;
  }

References CHISTO, CLINLIN, CLINLOG, CLOGLIN, CLOGLOG, G4Exp(), G4Log(), G4Pow::GetInstance(), HISTO, LINLIN, LINLOG, LOGLIN, LOGLOG, G4Pow::powA(), UHISTO, ULINLIN, ULINLOG, ULOGLIN, and ULOGLOG.

Referenced by G4ParticleHPVector::GetMeanX(), G4ParticleHPVector::IntegrateAndNormalise(), and G4ParticleHPContAngularPar::Sample().

GetWeightedBinIntegral()

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

Definition at line 96 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;
  }

References CHISTO, CLINLIN, CLINLOG, CLOGLIN, CLOGLOG, G4Exp(), G4Log(), G4Pow::GetInstance(), HISTO, LINLIN, LINLOG, LOGLIN, LOGLOG, G4Pow::powA(), UHISTO, ULINLIN, ULINLOG, ULOGLIN, and ULOGLOG.

Referenced by G4ParticleHPVector::GetMeanX(), and G4ParticleHPContAngularPar::Sample().

Histogram()

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

inlineprivate

Definition at line 162 of file G4ParticleHPInterpolator.hh.

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

Referenced by Interpolate2().

Interpolate()

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

inline

Definition at line 90 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;
}

References CSTART_, G4cout, G4endl, LinearLinear(), LinearLogarithmic(), LogarithmicLinear(), LogarithmicLogarithmic(), and Random().

Referenced by G4ParticleHPContAngularPar::BuildByInterpolation(), G4ParticleHPAngularP::GetCosTh(), G4ParticleHPVector::GetXsec(), G4ParticleHPPartial::GetY(), G4ParticleHPVector::Merge(), G4ParticleHPDiscreteTwoBody::Sample(), G4ParticleHPLabAngularEnergy::Sample(), G4ParticleHPContAngularPar::Sample(), G4ParticleHPLegendreStore::Sample(), G4ParticleHPPartial::Sample(), G4ParticleHPLegendreStore::SampleDiscreteTwoBody(), G4ParticleHPLegendreStore::SampleElastic(), and G4ParticleHPLegendreStore::SampleMax().

Interpolate2()

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

inline

Definition at line 127 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;
}

References CSTART_, G4cout, G4endl, Histogram(), LinearLinear(), LinearLogarithmic(), LogarithmicLinear(), LogarithmicLogarithmic(), and Random().

Referenced by G4ParticleHPContAngularPar::BuildByInterpolation().

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;
  }

Referenced by G4ParticleHPVector::Get50percentBorder(), G4ParticleHPVector::SampleLin(), and G4ParticleHPVector::ThinOut().

LinearLinear()

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

inlineprivate

Definition at line 170 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;
}

Referenced by Interpolate(), Interpolate2(), LinearLogarithmic(), LogarithmicLinear(), and LogarithmicLogarithmic().

LinearLogarithmic()

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

inlineprivate

Definition at line 181 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;
}

References G4Log(), and LinearLinear().

Referenced by Interpolate(), and Interpolate2().

LogarithmicLinear()

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

inlineprivate

Definition at line 192 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;
}

References G4Exp(), G4Log(), and LinearLinear().

Referenced by Interpolate(), and Interpolate2().

LogarithmicLogarithmic()

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

inlineprivate

Definition at line 205 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;
}

References G4Exp(), G4Log(), and LinearLinear().

Referenced by Interpolate(), and Interpolate2().

Random()

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

inlineprivate

Definition at line 221 of file G4ParticleHPInterpolator.hh.

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

References G4UniformRand.

Referenced by Interpolate(), and Interpolate2().

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPInterpolator.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPInterpolator.cc