# This python script plots histograms on the same plot 
# with 60 input files, of x y column format. Use like so:
#
# python plot60Hists.py distance_1.txt . . . distance_60.txt name_of_output_file.png


import matplotlib
matplotlib.use('Agg')
import numpy as np
import matplotlib.pyplot as plt
from pylab import *
import matplotlib.ticker as ticker
import sys

c1 = sys.argv[1]
c2 = sys.argv[2]
c3 = sys.argv[3]
c4 = sys.argv[4]
c5 = sys.argv[5]
c6 = sys.argv[6]
c7 = sys.argv[7]
c8 = sys.argv[8]
c9 = sys.argv[9]
c10 = sys.argv[10]
c11 = sys.argv[11]
c12 = sys.argv[12]
c13 = sys.argv[13]
c14 = sys.argv[14]
c15 = sys.argv[15]
c16 = sys.argv[16]
c17 = sys.argv[17]
c18 = sys.argv[18]
c19 = sys.argv[19]
c20 = sys.argv[20]
c21 = sys.argv[21]
c22 = sys.argv[22]
c23 = sys.argv[23]
c24 = sys.argv[24]
c25 = sys.argv[25]
c26 = sys.argv[26]
c27 = sys.argv[27]
c28 = sys.argv[28]
c29 = sys.argv[29]
c30 = sys.argv[30]
c31 = sys.argv[31]
c32 = sys.argv[32]
c33 = sys.argv[33]
c34 = sys.argv[34]
c35 = sys.argv[35]
c36 = sys.argv[36]
c37 = sys.argv[37]
c38 = sys.argv[38]
c39 = sys.argv[39]
c40 = sys.argv[40]
c41 = sys.argv[41]
c42 = sys.argv[42]
c43 = sys.argv[43]
c44 = sys.argv[44]
c45 = sys.argv[45]
c46 = sys.argv[46]
c47 = sys.argv[47]
c48 = sys.argv[48]
c49 = sys.argv[49]
c50 = sys.argv[50]
c51 = sys.argv[51]
c52 = sys.argv[52]
c53 = sys.argv[53]
c54 = sys.argv[54]
c55 = sys.argv[55]
c56 = sys.argv[56]
c57 = sys.argv[57]
c58 = sys.argv[58]
c59 = sys.argv[59]
c60 = sys.argv[60]

name = sys.argv[61]

x1,y1 = np.loadtxt(c1,skiprows=1).transpose()
x2,y2 = np.loadtxt(c2,skiprows=1).transpose()
x3,y3 = np.loadtxt(c3,skiprows=1).transpose()
x4,y4 = np.loadtxt(c4,skiprows=1).transpose()
x5,y5 = np.loadtxt(c5,skiprows=1).transpose()
x6,y6 = np.loadtxt(c6,skiprows=1).transpose()
x7,y7 = np.loadtxt(c7,skiprows=1).transpose()
x8,y8 = np.loadtxt(c8,skiprows=1).transpose()
x9,y9 = np.loadtxt(c9,skiprows=1).transpose()
x10,y10 = np.loadtxt(c10,skiprows=1).transpose()
x11,y11 = np.loadtxt(c11,skiprows=1).transpose()
x12,y12 = np.loadtxt(c12,skiprows=1).transpose()
x13,y13 = np.loadtxt(c13,skiprows=1).transpose()
x14,y14 = np.loadtxt(c14,skiprows=1).transpose()
x15,y15 = np.loadtxt(c15,skiprows=1).transpose()
x16,y16 = np.loadtxt(c16,skiprows=1).transpose()
x17,y17 = np.loadtxt(c17,skiprows=1).transpose()
x18,y18 = np.loadtxt(c18,skiprows=1).transpose()
x19,y19 = np.loadtxt(c19,skiprows=1).transpose()
x20,y20 = np.loadtxt(c20,skiprows=1).transpose()
x21,y21 = np.loadtxt(c21,skiprows=1).transpose()
x22,y22 = np.loadtxt(c22,skiprows=1).transpose()
x23,y23 = np.loadtxt(c23,skiprows=1).transpose()
x24,y24 = np.loadtxt(c24,skiprows=1).transpose()
x25,y25 = np.loadtxt(c25,skiprows=1).transpose()
x26,y26 = np.loadtxt(c26,skiprows=1).transpose()
x27,y27 = np.loadtxt(c27,skiprows=1).transpose()
x28,y28 = np.loadtxt(c28,skiprows=1).transpose()
x29,y29 = np.loadtxt(c29,skiprows=1).transpose()
x30,y30 = np.loadtxt(c30,skiprows=1).transpose()
x31,y31 = np.loadtxt(c31,skiprows=1).transpose()
x32,y32 = np.loadtxt(c32,skiprows=1).transpose()
x33,y33 = np.loadtxt(c33,skiprows=1).transpose()
x34,y34 = np.loadtxt(c34,skiprows=1).transpose()
x35,y35 = np.loadtxt(c35,skiprows=1).transpose()
x36,y36 = np.loadtxt(c36,skiprows=1).transpose()
x37,y37 = np.loadtxt(c37,skiprows=1).transpose()
x38,y38 = np.loadtxt(c38,skiprows=1).transpose()
x39,y39 = np.loadtxt(c39,skiprows=1).transpose()
x40,y40 = np.loadtxt(c40,skiprows=1).transpose()
x41,y41 = np.loadtxt(c41,skiprows=1).transpose()
x42,y42 = np.loadtxt(c42,skiprows=1).transpose()
x43,y43 = np.loadtxt(c43,skiprows=1).transpose()
x44,y44 = np.loadtxt(c44,skiprows=1).transpose()
x45,y45 = np.loadtxt(c45,skiprows=1).transpose()
x46,y46 = np.loadtxt(c46,skiprows=1).transpose()
x47,y47 = np.loadtxt(c47,skiprows=1).transpose()
x48,y48 = np.loadtxt(c48,skiprows=1).transpose()
x49,y49 = np.loadtxt(c49,skiprows=1).transpose()
x50,y50 = np.loadtxt(c50,skiprows=1).transpose()
x51,y51 = np.loadtxt(c51,skiprows=1).transpose()
x52,y52 = np.loadtxt(c52,skiprows=1).transpose()
x53,y53 = np.loadtxt(c53,skiprows=1).transpose()
x54,y54 = np.loadtxt(c54,skiprows=1).transpose()
x55,y55 = np.loadtxt(c55,skiprows=1).transpose()
x56,y56 = np.loadtxt(c56,skiprows=1).transpose()
x57,y57 = np.loadtxt(c57,skiprows=1).transpose()
x58,y58 = np.loadtxt(c58,skiprows=1).transpose()
x59,y59 = np.loadtxt(c59,skiprows=1).transpose()
x60,y60 = np.loadtxt(c60,skiprows=1).transpose()

ax = plt.axes()
plt.plot(x1,y1, linewidth=1)#, marker='.')
plt.plot(x2,y2, linewidth=1)#, marker='.')
plt.plot(x3,y3, linewidth=1)#, marker='.')
plt.plot(x4,y4, linewidth=1)#, marker='.')
plt.plot(x5,y5, linewidth=1)#, marker='.')
plt.plot(x6,y6, linewidth=1)#, marker='.')
plt.plot(x7,y7, linewidth=1)#, marker='.')
plt.plot(x8,y8, linewidth=1)#, marker='.')
plt.plot(x9,y9, linewidth=1)#, marker='.')
plt.plot(x10,y10, linewidth=1)#, marker='.')
plt.plot(x11,y11, linewidth=1)#, marker='.')
plt.plot(x12,y12, linewidth=1)#, marker='.')
plt.plot(x13,y13, linewidth=1)#, marker='.')
plt.plot(x14,y14, linewidth=1)#, marker='.')
plt.plot(x15,y15, linewidth=1)#, marker='.')
plt.plot(x16,y16, linewidth=1)#, marker='.')
plt.plot(x17,y17, linewidth=1)#, marker='.')
plt.plot(x18,y18, linewidth=1)#, marker='.')
plt.plot(x19,y19, linewidth=1)#, marker='.')
plt.plot(x20,y20, linewidth=1)#, marker='.')
plt.plot(x21,y21, linewidth=1)#, marker='.')
plt.plot(x22,y22, linewidth=1)#, marker='.')
plt.plot(x23,y23, linewidth=1)#, marker='.')
plt.plot(x24,y24, linewidth=1)#, marker='.')
plt.plot(x25,y25, linewidth=1)#, marker='.')
plt.plot(x26,y26, linewidth=1)#, marker='.')
plt.plot(x27,y27, linewidth=1)#, marker='.')
plt.plot(x28,y28, linewidth=1)#, marker='.')
plt.plot(x29,y29, linewidth=1)#, marker='.')
plt.plot(x30,y30, linewidth=1)#, marker='.')
plt.plot(x31,y31, linewidth=1)#, marker='.')
plt.plot(x32,y32, linewidth=1)#, marker='.')
plt.plot(x33,y33, linewidth=1)#, marker='.')
plt.plot(x34,y34, linewidth=1)#, marker='.')
plt.plot(x35,y35, linewidth=1)#, marker='.')
plt.plot(x36,y36, linewidth=1)#, marker='.')
plt.plot(x37,y37, linewidth=1)#, marker='.')
plt.plot(x38,y38, linewidth=1)#, marker='.')
plt.plot(x39,y39, linewidth=1)#, marker='.')
plt.plot(x40,y40, linewidth=1)#, marker='.')
plt.plot(x41,y41, linewidth=1)#, marker='.')
plt.plot(x42,y42, linewidth=1)#, marker='.')
plt.plot(x43,y43, linewidth=1)#, marker='.')
plt.plot(x44,y44, linewidth=1)#, marker='.')
plt.plot(x45,y45, linewidth=1)#, marker='.')
plt.plot(x46,y46, linewidth=1)#, marker='.')
plt.plot(x47,y47, linewidth=1)#, marker='.')
plt.plot(x48,y48, linewidth=1)#, marker='.')
plt.plot(x49,y49, linewidth=1)#, marker='.')
plt.plot(x50,y50, linewidth=1)#, marker='.')
plt.plot(x51,y51, linewidth=1)#, marker='.')
plt.plot(x52,y52, linewidth=1)#, marker='.')
plt.plot(x53,y53, linewidth=1)#, marker='.')
plt.plot(x54,y54, linewidth=1)#, marker='.')
plt.plot(x55,y55, linewidth=1)#, marker='.')
plt.plot(x56,y56, linewidth=1)#, marker='.')
plt.plot(x57,y57, linewidth=1)#, marker='.')
plt.plot(x58,y58, linewidth=1)#, marker='.')
plt.plot(x59,y59, linewidth=1)#, marker='.')
plt.plot(x60,y60, linewidth=1)#, marker='.')



plt.title("Histograms Title")
#ax.xaxis.set_major_locator(ticker.MultipleLocator(5))
#ax.xaxis.set_minor_locator(ticker.MultipleLocator(1))

plt.xlabel(r"Distance (A)")
plt.ylabel("Population")

figsave=name
plt.savefig(figsave,dpi=300)