eventaidsmortality.h

#ifndef EVENTAIDSMORTALITY_H

#define EVENTAIDSMORTALITY_H

#include "eventmortalitybase.h"
#include "eventvariablefiretime.h"

// AIDS mortality
class EventAIDSMortality : public EventMortalityBase
{
public:
        EventAIDSMortality(Person *pPerson);
        ~EventAIDSMortality();

        std::string getDescription(double tNow) const;
        void writeLogs(double tNow) const;
#ifndef NDEBUG
        void fire(State *pState, double t);
#endif // !NDEBUG

        static void processConfig(ConfigSettings &config);
        static void obtainConfig(ConfigWriter &config);

        static double getExpectedSurvivalTime(const Person *pPerson);
        static double getExpectedTimeOfDeath(const Person *pPerson);
private:
        double getNewInternalTimeDifference(GslRandomNumberGenerator *pRndGen, const State *pState);
        double calculateInternalTimeInterval(const State *pState, double t0, double dt);
        double solveForRealTimeInterval(const State *pState, double Tdiff, double t0);
        void checkFireTime(double t0);

        EventVariableFireTime_Helper m_eventHelper;

        // TODO: use access functions
        static double m_C;
        static double m_k;
};

inline double EventAIDSMortality::getExpectedSurvivalTime(const Person *pPerson)
{
        assert(pPerson);
        double Vsp = pPerson->getSetPointViralLoad();
        assert(Vsp > 0);

        double tSurvival = m_C/std::pow(Vsp, -m_k);
        assert(tSurvival > 0);

        return tSurvival;
}

inline double EventAIDSMortality::getExpectedTimeOfDeath(const Person *pPerson)
{
        assert(pPerson);
        assert(pPerson->isInfected());

        double expectedTimeOfDeath = getExpectedSurvivalTime(pPerson) + pPerson->getInfectionTime();
        
        return expectedTimeOfDeath;
}

#endif // EVENTMORTALITY_H