Refactored simulation code to use np arrays from the start.  Added reindexing and new column to track the actual lineage of events.
TODO: complete the func for p, the true connectivity.
Confused about the norming - is this how Mohler did it?
Runner now working in a loop BUT producing strange results.