Add functions to read in file and create density plot