phel_scattering_matrix.h

#ifndef PHEL_SCATTERING_MATRIX_H
#define PHEL_SCATTERING_MATRIX_H
 
#include "electron_h0_wannier.h"
#include "phonon_h0.h"
#include "base_ph_scattering_matrix.h"
#include "vector_bte.h"
#include "interaction_elph.h"
 
// TODO realistically this should be just absorbed
// as a diagonal only case for ph scattering matrix!
class PhElScatteringMatrix : public BasePhScatteringMatrix {
public:
 
  PhElScatteringMatrix(Context &context_, StatisticsSweep &statisticsSweep_,
                       BaseBandStructure &phBandStructure_,
                       PhononH0 *phononH0_, 
                       ElectronH0Wannier *electronH0_);
 
  PhElScatteringMatrix(const PhElScatteringMatrix &that);
 
  PhElScatteringMatrix &operator=(const PhElScatteringMatrix &that);
 
protected:
 
  PhononH0* phononH0; 
  ElectronH0Wannier *electronH0;
 
  void builder(std::shared_ptr<VectorBTE> linewidth, std::vector<VectorBTE> &inPopulations,
               std::vector<VectorBTE> &outPopulations) override;
 
  // to prevent mistakes in which these two outer and inner BS could
  // be accidentally swapped
  BaseBandStructure& getPhBandStructure() { return outerBandStructure; };
  BaseBandStructure& getElBandStructure() { return innerBandStructure; };
 
  friend void addPhElScattering(BasePhScatteringMatrix &matrix, Context &context,
              std::vector<VectorBTE> &inPopulations,
              std::vector<VectorBTE> &outPopulations,
              BaseBandStructure &phBandStructure,
              BaseBandStructure &elBandStructure,
              StatisticsSweep& statisticsSweep, 
              InteractionElPhWan &couplingElPhWan,
              std::shared_ptr<VectorBTE> linewidth);
 
};
 
#endif