ParticleNumberOperator.hpp

Go to the documentation of this file.

#ifndef OBSERVABLE_PARTICLENUMBEROPERATOR_HPP
#define OBSERVABLE_PARTICLENUMBEROPERATOR_HPP

#include "Observable.hpp"

namespace solid
{

template <template <typename> class T1, typename T2>
class ParticleNumberOperator : public Observable<T1, T2>, public IOperatorParameters<T1, T2>
{
    using Observable<T1, T2>::_operator;
    using Observable<T1, T2>::parameters;

public:
    ParticleNumberOperator<T1, T2>(int _L)
    {
        L = _L;
        SelectTerms();
        FillParameters();
    }

private:
    int L;

    void SelectTerms()
    {
        _operator.termsEnabled.terms = {TermsEnum::FermionOneBodyInteraction};
    }

    void FillParameters() override
    {
        arma::SpMat<T2> mu = arma::eye<arma::SpMat<T2>>(L, L);
        parameters["mu"] = mu;
    }
};

} // namespace solid

#endif