36#include "RandomCellKiller.hpp"
41 mProbabilityOfDeathInAnHour(probabilityOfDeathInAnHour)
45 EXCEPTION(
"Probability of death must be between zero and one");
52 return mProbabilityOfDeathInAnHour;
71 double death_prob_this_timestep = 1.0 - pow((1.0 - mProbabilityOfDeathInAnHour),
SimulationTime::Instance()->GetTimeStep());
73 if (!pCell->HasApoptosisBegun() &&
77 this->mpCellPopulation->StartApoptosisOnCell(pCell,
"RandomCellKiller");
85 cell_iter != this->mpCellPopulation->End();
88 CheckAndLabelSingleCellForApoptosis(*cell_iter);
95 *rParamsFile <<
"\t\t\t<ProbabilityOfDeathInAnHour>" << mProbabilityOfDeathInAnHour <<
"</ProbabilityOfDeathInAnHour>\n";
#define EXCEPTION(message)
#define EXPORT_TEMPLATE_CLASS_SAME_DIMS(CLASS)
virtual void OutputCellKillerParameters(out_stream &rParamsFile)=0
double GetDeathProbabilityInAnHour() const
void OutputCellKillerParameters(out_stream &rParamsFile)
RandomCellKiller(AbstractCellPopulation< DIM > *pCellPopulation, double probabilityOfDeathInAnHour)
double mProbabilityOfDeathInAnHour
void CheckAndLabelCellsForApoptosisOrDeath()
void CheckAndLabelSingleCellForApoptosis(CellPtr pCell)
static RandomNumberGenerator * Instance()
static SimulationTime * Instance()