dos.hpp

Go to the documentation of this file.

// Copyright 2015-2018 The ALMA Project Developers
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
// implied. See the License for the specific language governing
// permissions and limitations under the License.

#pragma once


#include <qpoint_grid.hpp>
#include <boost/math/distributions/normal.hpp>

namespace alma {
class Gaussian_for_DOS {
public:
    const double mu;
    const double sigma;
    const bool truncated;
    const double lbound;
    const double ubound;
    Gaussian_for_DOS(const Gamma_grid& grid,
                     std::size_t iq,
                     std::size_t im,
                     double scalebroad)
        : mu(grid.get_spectrum_at_q(iq).omega(im)),
          sigma(scalebroad *
                grid.base_sigma(grid.get_spectrum_at_q(iq).vg.col(im))),
          truncated(mu - constants::nsigma * sigma < 0.),
          lbound(truncated ? 0. : mu - constants::nsigma * sigma),
          ubound(mu + constants::nsigma * sigma) {
        if ((this->mu == 0.) || (this->sigma == 0.))
            this->dist = nullptr;
        else
            this->dist = make_unique<boost::math::normal>(mu, sigma);
    }


    double get_contribution(double omega) const {
        if (this->dist)
            return boost::math::pdf(*(this->dist), omega);
        else
            return 0.;
    }


private:
    std::unique_ptr<boost::math::normal> dist;
};
} // namespace alma