ContactHandler.h

Go to the documentation of this file.

/*
 *  This is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  any later version.
 *  The software is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *  You should have received a copy of the GNU General Public License
 *  along with the software. If not, see <http://www.gnu.org/licenses/>.
 *
 *  Copyright 2017, Kuylen E, Willem L, Broeckhove J
 */

#pragma once

#include <cmath>
#include <functional>

namespace stride {

class ContactHandler
{
public:
        explicit ContactHandler(std::function<double()> gen) : m_uniform01_generator(std::move(gen)) {}

        double operator()() { return m_uniform01_generator(); }

        bool HasContactAndTransmission(double contact_rate, double transmission_probability)
        {
                return m_uniform01_generator() < transmission_probability * RateToProbability(contact_rate);
        }

        bool HasContact(double contact_rate) { return m_uniform01_generator() < RateToProbability(contact_rate); }

        bool HasTransmission(double transmission_probability)
        {
                return m_uniform01_generator() < transmission_probability;
        }

private:
        double RateToProbability(double rate) { return 1.0 - std::exp(-rate); }

private:
        std::function<double()> m_uniform01_generator; 
};

} // namespace stride