particle.h

#ifndef PARTICLE_H
#define PARTICLE_H
 
#include "eigen.h"
 
 
class Particle {
public:
  static const int electron = -1;
  static const int phonon = -2;
 
  explicit Particle(int kind_);
 
  double getPopulation(const double &energy, const double &temperature,
                       const double &chemicalPotential = 0.) const;
 
  double getDndt(const double &energy, const double &temperature,
                 const double &chemicalPotential = 0.,
                 const bool &symmetrize = false) const;
 
  double getDnde(const double &energy, const double &temperature,
                 const double &chemicalPotential = 0.,
                 const bool &symmetrize = false) const;
 
  double getPopPopPm1(const double &energy, const double &temperature,
                      const double &chemicalPotential = 0.) const;
 
  bool isFermi() const;
 
  bool isBose() const;
 
  bool isElectron() const;
 
  bool isPhonon() const;
 
  int getParticleKind() const;
 
private:
  static const int fermi = 0;
  static const int bose = 1;
 
  // identifies whether the instance is a fermion of a boson
  int statistics;
  // identifies whether we have an electron or a phonon
  // Note: we keep particle != statistics because we might add particles
  int kind;
};
 
#endif