window.h

#ifndef WINDOW_H
#define WINDOW_H
 
#include "context.h"
#include "eigen.h"
#include "particle.h"
 
class Window {
 public:
  Window(Context &context, Particle &particle_, const double &temperatureMin =
  0., const double &temperatureMax = 0.,
         const double &chemicalPotentialMin = 0.,
         const double &chemicalPotentialMax = 0.);
 
  std::tuple<std::vector<double>, std::vector<int>> apply(
      Eigen::VectorXd &energies);
 
  // Constants that identify the kind of filter to be used
  static const int nothing = 0;
 
  static const int population = 1;
 
  static const int energy = 2;
 
  static const int muEnergy = 3;
 
  int getMethodUsed() const;
 private:
  Particle &particle;
 
  // parameters for window
  double temperatureMin, temperatureMax;
  double chemicalPotentialMin, chemicalPotentialMax;
  double populationThreshold = 0.;
  double minEnergy = 0., maxEnergy = 0.;
 
  int method;
 
  int numBands = 0;
 
  // internal method to apply the window on population
  std::tuple<std::vector<double>, std::vector<int>> internalPopWindow(
      const Eigen::VectorXd &energies, const Eigen::VectorXd &popMin,
      const Eigen::VectorXd &popMax) const;
  // internal method to apply the window on energy
  std::tuple<std::vector<double>, std::vector<int>> internalEnWindow(
      const Eigen::VectorXd &energies) const;
};
 
#endif