ComponentUser.cc

Go to the documentation of this file.

#include <iostream>
 
#include "ComponentUser.hh"
 
namespace Garfield {
 
ComponentUser::ComponentUser()
    : ComponentBase(),
      hasField(false),
      field(0),
      hasPotential(false),
      potential(0),
      hasWeightingField(false),
      wfield(0),
      hasWeightingPotential(false),
      wpot(0) {
 
  m_className = "ComponentUser";
}
 
void ComponentUser::ElectricField(const double x, const double y,
                                  const double z, double& ex, double& ey,
                                  double& ez, Medium*& m, int& status) {
 
  if (!hasField) {
    ex = ey = ez = 0.;
    m = 0;
    status = -10;
    return;
  }
 
  field(x, y, z, ex, ey, ez);
  m = GetMedium(x, y, z);
  if (m == NULL) {
    if (debug) {
      std::cerr << m_className << "::ElectricField:\n";
      std::cerr << "    (" << x << ", " << y << ", " << z << ")"
                << " is not inside a medium.\n";
    }
    status = -6;
    return;
  }
 
  if (m->IsDriftable()) {
    status = 0;
  } else {
    status = -5;
  }
}
 
void ComponentUser::ElectricField(const double x, const double y,
                                  const double z, double& ex, double& ey,
                                  double& ez, double& v, Medium*& m,
                                  int& status) {
 
  if (!hasField) {
    ex = ey = ez = v = 0.;
    m = 0;
    status = -10;
    return;
  }
  field(x, y, z, ex, ey, ez);
 
  if (hasPotential) {
    potential(x, y, z, v);
  } else {
    v = 0.;
  }
 
  m = GetMedium(x, y, z);
  if (m == NULL) {
    if (debug) {
      std::cerr << m_className << "::ElectricField:\n";
      std::cerr << "    (" << x << ", " << y << ", " << z << ")"
                << " is not inside a medium.\n";
    }
    status = -6;
    return;
  }
 
  if (m->IsDriftable()) {
    status = 0;
  } else {
    status = -5;
  }
}
 
bool ComponentUser::GetVoltageRange(double& vmin, double& vmax) {
 
  vmin = vmax = 0.;
  return false;
}
 
void ComponentUser::WeightingField(const double x, const double y,
                                   const double z, double& wx, double& wy,
                                   double& wz, const std::string label) {
 
  wx = wy = wz = 0.;
  if (!hasWeightingField) return;
  wfield(x, y, z, wx, wy, wz, label);
}
 
double ComponentUser::WeightingPotential(const double x, const double y,
                                         const double z,
                                         const std::string label) {
 
  double v = 0.;
  if (hasWeightingPotential) {
    wpot(x, y, z, v, label);
  }
  return v;
}
 
void ComponentUser::SetElectricField(void (*f)(const double, const double,
                                               const double, double&, double&,
                                               double&)) {
 
  if (f == 0) {
    std::cerr << m_className << "::SetElectricField:\n";
    std::cerr << "    Function pointer is null.\n";
    return;
  }
  field = f;
  hasField = true;
  ready = true;
}
 
void ComponentUser::SetPotential(void (*f)(const double, const double,
                                           const double, double&)) {
 
  if (f == 0) {
    std::cerr << m_className << "::SetPotential:\n";
    std::cerr << "    Function pointer is null.\n";
    return;
  }
  potential = f;
  hasPotential = true;
}
 
void ComponentUser::SetWeightingField(void (*f)(const double, const double,
                                                const double, double&, double&,
                                                double&, const std::string)) {
 
  if (f == 0) {
    std::cerr << m_className << "::SetWeightingField:\n";
    std::cerr << "    Function pointer is null.\n";
    return;
  }
  wfield = f;
  hasWeightingField = true;
}
 
void ComponentUser::SetWeightingPotential(void (*f)(const double, const double,
                                                    const double, double&,
                                                    const std::string)) {
 
  if (f == 0) {
    std::cerr << m_className << "::SetWeightingPotential:\n";
    std::cerr << "    Function pointer is null.\n";
    return;
  }
  wpot = f;
  hasWeightingPotential = true;
}
 
void ComponentUser::Reset() {
 
  field = 0;
  potential = 0;
  wfield = 0;
  wpot = 0;
  hasField = false;
  hasPotential = false;
  hasWeightingField = false;
  hasWeightingPotential = false;
  ready = false;
}
 
void ComponentUser::UpdatePeriodicity() {
 
  if (debug) {
    std::cerr << m_className << "::UpdatePeriodicity:\n";
    std::cerr << "    Periodicities are not supported.\n";
  }
}
}