EDepSim::ArbElecField Class Reference

#include <EDepSimArbElecField.hh>

Inheritance diagram for EDepSim::ArbElecField:

Public Member Functions
	ArbElecField ()
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 59 of file EDepSimArbElecField.hh.

Constructor & Destructor Documentation

EDepSim::ArbElecField::ArbElecField

(

)

Definition at line 38 of file EDepSimArbElecField.cc.

{
}

Member Function Documentation

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

virtual

Definition at line 116 of file EDepSimArbElecField.cc.

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

    EDepSim::Cubic interp;
    field[3] = 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[4] = 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[5] = 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::ArbElecField::PrintInfo

(

)

const

Definition at line 128 of file EDepSimArbElecField.cc.

{
    EDepSimLog("Printing values for electric 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::ArbElecField::ReadFile

(

const std::string &

fname

)

Definition at line 42 of file EDepSimArbElecField.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 * (volt/cm));
            m_field_y[xcount][ycount].push_back(fy * (volt/cm));
            m_field_z[xcount][ycount].push_back(fz * (volt/cm));

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

    return true;
}

Member Data Documentation

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

private

Definition at line 71 of file EDepSimArbElecField.hh.

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

private

Definition at line 72 of file EDepSimArbElecField.hh.

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

private

Definition at line 73 of file EDepSimArbElecField.hh.

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

private

Definition at line 74 of file EDepSimArbElecField.hh.

std::string EDepSim::ArbElecField::m_filename

private

Definition at line 69 of file EDepSimArbElecField.hh.

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

private

Definition at line 70 of file EDepSimArbElecField.hh.

---

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

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