Now the rooted triplet distribution works when branch length information needs to be kept (as for MPest). Had to change the birth death simulate function, because it did not keep species information at the tips, and added clone and "=" function to the Taxon class. Not sure the "==" function in the Taxon class works by the way (I have not changed it though).