newmen/versatile-diamond

#include "base_events_container.h"

namespace vd
{

Reaction *BaseEventsContainer::selectEvent(double r)
{
    assert(_events.size() > 0);

    uint index = (uint)(r / _events.front()->rate());
    assert(index < _events.size());

#if defined(PRINT) || defined(MC_PRINT)
    debugPrint([&](IndentStream &os) {
        os << "BaseEventsContainer::selectEvent: " << index;
    });
#endif // PRINT || MC_PRINT

    return _events[index];
}

double BaseEventsContainer::oneRate() const
{
    assert(_events.front());
    return _events.front()->rate();
}

double BaseEventsContainer::commonRate() const
{
    if (_events.empty())
    {
        return 0;
    }
    else
    {
        return _events.size() * oneRate();
    }
}

}