import pylab as p
import numpy as np
import csv
from tCal import timeToChannel, channelToTime
Hs=[]
t=[channelToTime(i)*1e-9 for i in range(8190)]
for i in range(2):
with open('muonAndPhoton'+['','2'][i]) as f:
Hs.append([])
for line in f:
Hs[-1].append(int(line.split()[1]))
if len(Hs[-1])==8190:
break
H=np.sum(Hs,axis=0)
channelSize=[(t[i+1]-t[i-1])/2 for i in range(1,len(t)-1)]
channelSize=[channelSize[0]]+channelSize+[channelSize[-1]]
pdfPoints=[(0. if h!=h or h==float('inf') else h) for h in H/channelSize]
noiseT1=3.75e-8
noiseT2=3.99e-8
noise=np.mean([pdfPoints[i] for i in range(len(pdfPoints)) if noiseT1<t[i]<noiseT2])
pdfPoints=pdfPoints-noise
pdfPoints=[0]+pdfPoints/sum(pdfPoints*channelSize)+[0]
pdf=lambda it: np.interp(it, [t[0]-channelSize[0]]+t+[t[-1]+channelSize[-1]], [0]+list(pdfPoints)+[0])
if __name__=='__main__':
#n=200
#ts=np.linspace(0,50e-9,n+1)
import pylab as p
import numpy as np
import csv
import tCal
with open('dCal') as f:
H=[]
for line in f:
H.append(int(line.split()[1]))
if len(H)==8190:
break
n=200
ts=np.linspace(0,50,n+1)
N=[sum(H[j] for j in range(len(H)) if ts[i]<tCal.channelToTime(j)-6<ts[i+1]) for i in range(n)]
p.plot(reduce(lambda a,b:a+b,[[t]*2 for t in ts]),[0]+reduce(lambda a,b:a+b,[[iN]*2 for iN in N])+[0])
p.show()
import pylab as p
import numpy as np
import csv
from tCal import timeToChannel, channelToTime
delay=8
ct=[]
with open('muonSpeed2') as f:
H=[]
for line in f:
H.append(int(line.split()[1]))
ct.append(channelToTime(len(ct))-delay)
if len(H)==8190:
break
n=200
ts=np.linspace(0,50,n+1)
N=[sum(H[j] for j in range(len(H)) if ts[i]<ct[j]<ts[i+1]) for i in range(n)]
#p.plot(reduce(lambda a,b:a+b,[[t]*2 for t in ts]),[0]+reduce(lambda a,b:a+b,[[iN]*2 for iN in N])+[0])
p.plot(N)
p.show()
