HeedParticle_BGM.h

Go to the documentation of this file.

#ifndef HEEDPARTICLE_BGM_H
#define HEEDPARTICLE_BGM_H
 
#include <vector>
#include "HeedCluster.h"
#include "wcpplib/particle/eparticle.h"
 
namespace Heed {
extern long last_particle_number;
 
/// Definition of the particle which can be traced through the geometry.
/// 2003, I. Smirnov
 
class HeedParticle_BGM : public eparticle {
 public:
  /// Default constructor.
  HeedParticle_BGM() : eparticle() {}
  /// Constructor.
  /// if fs_loss_only == 1 - only transfer energy and
  /// no other physics: no deposition of clusters,
  /// no generation of virtual photons.
  /// Thus it is just a PAI without even clusters
  HeedParticle_BGM(manip_absvol* primvol, const point& pt, const vec& vel,
                   vfloat time, particle_def* fpardef, HeedFieldMap* fieldmap,
                   bool fs_loss_only = false, bool fs_print_listing = false);
  /// Destructor
  virtual ~HeedParticle_BGM() {}
 
  virtual void physics(std::vector<gparticle*>& secondaries);
  virtual void print(std::ostream& file, int l) const;
  virtual HeedParticle_BGM* copy() const { return new HeedParticle_BGM(*this); }
 
 private:
  bool s_print_listing;
  long particle_number;
 
  bool s_loss_only;
  std::vector<double> etransf;
  std::vector<long> natom;
  std::vector<long> nshell;
 
  std::vector<HeedCluster> m_clusterBank;
};
}
 
#endif