Add in code from massPlots.C. ALso update mu1PT and mu2PT variables so as not to cause confusion (they stay pT ordered)