mag::Interpolator Class Reference

#include <GArMagneticField.h>

Public Member Functions
	Interpolator ()
float	interpolate (const float *point, const std::vector< std::vector< std::vector< float >>> &g, const float *delta, const float *offset) const
float	interpolate (float x, float y, float z, const std::vector< std::vector< std::vector< float >>> &g, float hx, float hy, float hz, float xo, float yo, float zo) const

Private Member Functions
float	conv_kernel (float s) const

Detailed Description

Definition at line 67 of file GArMagneticField.h.

Constructor & Destructor Documentation

mag::Interpolator::Interpolator

(

)

Definition at line 464 of file GArMagneticField.cxx.

{}

Member Function Documentation

float mag::Interpolator::conv_kernel ( float  s ) const

private

Definition at line 518 of file GArMagneticField.cxx.

{
    float v = 0;
    float z = std::abs(s);

    if(0 <= z && z < 1)
        v = 1 + (0.5) * z * z * (3 * z - 5);

    else if(1 < z && z < 2)
        v = 2 - z * (4 + 0.5 * z * (z - 5));

    return v;
}

float mag::Interpolator::interpolate	(	const float *	point,
		const std::vector< std::vector< std::vector< float >>> &	g,
		const float *	delta,
		const float *	offset
	)		const

Definition at line 467 of file GArMagneticField.cxx.

{
    return interpolate(point[0], point[1], point[2], g, delta[0], delta[1], delta[2], offset[0],
                       offset[1], offset[2]);
}

float mag::Interpolator::interpolate	(	float	x,
		float	y,
		float	z,
		const std::vector< std::vector< std::vector< float >>> &	g,
		float	hx,
		float	hy,
		float	hz,
		float	xo,
		float	yo,
		float	zo
	)		const

Definition at line 476 of file GArMagneticField.cxx.

{
    float v = 0;

    const int xsize = g.size();
    const int ysize = g[0].size();
    const int zsize = g[0][0].size();

    const float xp = (x - xo) / hx;
    const float yp = (y - yo) / hy;
    const float zp = (z - zo) / hz;

    const int x_idx = std::floor(xp);
    const int y_idx = std::floor(yp);
    const int z_idx = std::floor(zp);

    for(auto i : {x_idx - 1, x_idx, x_idx + 1, x_idx + 2})
    {
        for(auto j : {y_idx - 1, y_idx, y_idx + 1, y_idx + 2})
        {
            for(auto k : {z_idx - 1, z_idx, z_idx + 1, z_idx + 2})
            {
                if(i < 0 || j < 0 || k < 0)
                    continue;

                if(i >= xsize || j >= ysize || k >= zsize)
                    continue;

                const float x_c = conv_kernel(xp - i);
                const float y_c = conv_kernel(yp - j);
                const float z_c = conv_kernel(zp - k);
                v += g[i][j][k] * x_c * y_c * z_c;
            }
        }
    }

    return v;
}

---

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

• garsoft/DetectorInfo/GArMagneticField.h
• garsoft/DetectorInfo/GArMagneticField.cxx