eparticle.cpp

Go to the documentation of this file.

#include "wcpplib/particle/eparticle.h"
#include "wcpplib/util/emul_new_stand.h"
/*
1998 - 2004,   I. Smirnov
*/
 
namespace Heed {
 
void field_map(const point& pt, vec& Efield, vec& Hfield, vfloat& mrange);
// defined anywhere outside
 
eparticle::eparticle(manip_absvol* primvol, const point& pt, const vec& vel,
                     vfloat time, particle_def* fpardef)
    : mparticle(), particle_type(fpardef) {
  gparticle gp(primvol, pt, vel, time);
  statcast(mparticle&, *this) = mparticle(gp, fpardef->mass);
}
 
eparticle::eparticle(manip_absvol* primvol, const point& pt, const vec& vel,
                     vfloat time, particle_def* fpardef, double gamma_1)
    : mparticle(primvol, pt, vel, time, fpardef->mass, gamma_1),
      particle_type(fpardef) {}
 
int eparticle::force(const point& pt, vec& f, vec& f_perp, vfloat& mrange) {
  vec efield;
  vec hfield;
  field_map(pt, efield, hfield, mrange);
  f = pardef->charge * efield;
  f_perp = pardef->charge * hfield;
  return 1;
}
 
void eparticle::print(std::ostream& file, int l) const {
  if (l >= 0) {
    Ifile << "eparticle: particle is ";
    print_notation(file);
    file << '\n';
    mparticle::print(file, l);
  }
}
 
}