EDepSim::ArbMagField Class Reference

#include <EDepSimArbMagField.hh>

Inheritance diagram for EDepSim::ArbMagField:

Public Member Functions
	ArbMagField ()
bool	ReadFile (const std::string &fname)
void	PrintInfo () const
virtual void	GetFieldValue (const G4double pos[4], G4double *field) const

Private Attributes
std::string	m_filename
std::array< double, 3 >	m_offset
std::array< double, 3 >	m_delta
std::vector< std::vector< std::vector< double > > >	m_field_x
std::vector< std::vector< std::vector< double > > >	m_field_y
std::vector< std::vector< std::vector< double > > >	m_field_z

Detailed Description

Definition at line 60 of file EDepSimArbMagField.hh.

Constructor & Destructor Documentation

EDepSim::ArbMagField::ArbMagField

(

)

Definition at line 38 of file EDepSimArbMagField.cc.

{
}

Member Function Documentation

void EDepSim::ArbMagField::GetFieldValue ( const G4double  pos[4], G4double *  field  ) const

virtual

Definition at line 116 of file EDepSimArbMagField.cc.

{
    const double x = pos[0];
    const double y = pos[1];
    const double z = pos[2];

    EDepSim::Cubic interp;
    field[0] = interp.interpolate(x, y, z, m_field_x, m_delta[0], m_delta[1], m_delta[2], m_offset[0], m_offset[1], m_offset[2]);
    field[1] = interp.interpolate(x, y, z, m_field_y, m_delta[0], m_delta[1], m_delta[2], m_offset[0], m_offset[1], m_offset[2]);
    field[2] = interp.interpolate(x, y, z, m_field_z, m_delta[0], m_delta[1], m_delta[2], m_offset[0], m_offset[1], m_offset[2]);
}

void EDepSim::ArbMagField::PrintInfo

(

)

const

Definition at line 128 of file EDepSimArbMagField.cc.

{
    EDepSimLog("Printing values for magnetic field.");
    EDepSimLog("m_filename : " << m_filename
              << "\nm_offset   : " << m_offset[0] << ", " << m_offset[1] << ", " << m_offset[2]
              << "\nm_delta    : " << m_delta[0] << ", " << m_delta[1] << ", " << m_delta[2]);
}

bool EDepSim::ArbMagField::ReadFile

(

const std::string &

fname

)

Definition at line 42 of file EDepSimArbMagField.cc.

{
    m_filename = fname;
    std::fstream fin(fname, std::fstream::in);

    if(!fin.is_open())
    {
        EDepSimError("Can't read " << fname << std::endl);
        return false;
    }
    else
    {
        EDepSimLog("Reading " << fname << " ...");
        int xcount{-1}, ycount{-1}, zcount{-1};
        double xcurr{0}, ycurr{0}, zcurr{0};
        std::string line;

        while(std::getline(fin >> std::ws, line))
        {
            std::stringstream ss(line);

            if(ss.str().front() == '#')
                continue;

            ss >> m_offset[0] >> m_offset[1] >> m_offset[2]
                >> m_delta[0] >> m_delta[1] >> m_delta[2];
            break;
        }

        while(std::getline(fin >> std::ws, line))
        {
            double x{0}, y{0}, z{0}, fx{0}, fy{0}, fz{0}, f{0};
            std::stringstream ss(line);

            if(ss.str().front() == '#')
                continue;

            ss >> x >> y >> z >> fx >> fy >> fz >> f;

            if(std::abs(x - xcurr) > 0.0 || xcount < 0)
            {
                xcurr = x;
                xcount += 1;
                ycount = -1;
                m_field_x.emplace_back(std::vector<std::vector<double>>{});
                m_field_y.emplace_back(std::vector<std::vector<double>>{});
                m_field_z.emplace_back(std::vector<std::vector<double>>{});
            }

            if(std::abs(y - ycurr) > 0.0 || ycount < 0)
            {
                ycurr = y;
                ycount += 1;
                zcount = -1;
                m_field_x[xcount].emplace_back(std::vector<double>{});
                m_field_y[xcount].emplace_back(std::vector<double>{});
                m_field_z[xcount].emplace_back(std::vector<double>{});
            }

            m_field_x[xcount][ycount].push_back(fx * tesla);
            m_field_y[xcount][ycount].push_back(fy * tesla);
            m_field_z[xcount][ycount].push_back(fz * tesla);

            if(std::abs(z - zcurr) > 0.0 || zcount < 0)
            {
                zcurr = z;
                zcount += 1;
            }
        }
    }

    return true;
}

Member Data Documentation

std::array<double, 3> EDepSim::ArbMagField::m_delta

private

Definition at line 72 of file EDepSimArbMagField.hh.

std::vector<std::vector<std::vector<double> > > EDepSim::ArbMagField::m_field_x

private

Definition at line 73 of file EDepSimArbMagField.hh.

std::vector<std::vector<std::vector<double> > > EDepSim::ArbMagField::m_field_y

private

Definition at line 74 of file EDepSimArbMagField.hh.

std::vector<std::vector<std::vector<double> > > EDepSim::ArbMagField::m_field_z

private

Definition at line 75 of file EDepSimArbMagField.hh.

std::string EDepSim::ArbMagField::m_filename

private

Definition at line 70 of file EDepSimArbMagField.hh.

std::array<double, 3> EDepSim::ArbMagField::m_offset

private

Definition at line 71 of file EDepSimArbMagField.hh.

---

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

• edep-sim/src/EDepSimArbMagField.hh
• edep-sim/src/EDepSimArbMagField.cc