Skip to content
This repository was archived by the owner on Nov 6, 2019. It is now read-only.
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
11 changes: 6 additions & 5 deletions fbu/monitoring.py
Original file line number Diff line number Diff line change
@@ -1,14 +1,15 @@
import matplotlib.mlab as mlab
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

remove unused import

import matplotlib.pyplot as plt
from numpy import mean,std,arange,array
from scipy.stats import norm

def plothistandtrace(name,xx,lower,upper):
ax = plt.subplot(211)
mu = mean(xx) if 'truth' in name else 0.
sigma = std(xx) if 'truth' in name else 1.
n, bins, patches = plt.hist(xx, bins=50, normed=1, facecolor='green',
n, bins, patches = plt.hist(xx, bins=50, normed=1, facecolor='green',
alpha=0.5, histtype='stepfilled')
yy = mlab.normpdf(bins,mu,sigma)
yy = norm.pdf(bins,mu,sigma)
plt.plot(bins,yy,'r-')
plt.ylabel('Probability')
plt.xlabel('Bin content')
Expand Down Expand Up @@ -43,12 +44,12 @@ def plot(dirname,data,bkgd,resmat,trace,nuisancetrace,lower=[],upper=[]):
plt.close()

for name,nuisance in nuisancetrace.items():
plothistandtrace(dirname+name,nuisance,-5.,5.)
plothistandtrace(dirname+name,nuisance,-5.,5.)

nbins = len(trace)
for bin in range(nbins):
for bin in range(nbins):
plothistandtrace(dirname+'bin%d'%bin,trace[bin],lower[bin],upper[bin])

for name,nuisance in nuisancetrace.items():
plt.plot(trace[bin],nuisance,',')
plt.savefig(dirname+'%s_bin%d.png'%(name,bin))
Expand Down