import ROOT,argparse




pfilename="tpc/FOPI/par/padplane_FOPI.dat"

pfile=open(pfilename,'r')

xpos=[]
ypos=[]
hist=ROOT.TH2D("padp","padp",3000,-15,15,3000,-15,15)
for line in pfile:
    if line[0]=='P':
        continue
    words=line.split()
    xpos.append(float(words[3]))
    ypos.append(float(words[4]))
    hist.Fill(xpos[-1],ypos[-1])


hist.Draw("colz")

u=raw_input("bla")