wirecell.gen.test.plot_impactzipper Namespace Reference

Functions
def	bilogify (hist, lmin=None)
def	set_palette (which="custom")

Variables
	stops = array('d',[ 0.00, 0.45, 0.50, 0.55, 1.00 ])
	reds = array('d',[ 0.00, 0.00, 1.00, 1.00, 0.51 ])
	greens = array('d',[ 0.00, 0.81, 1.00, 0.20, 0.00 ])
	blues = array('d',[ 0.51, 1.00, 1.00, 0.00, 0.00 ])
	fp = ROOT.TFile.Open("build/gen/test_impactzipper-uvw.root")
dictionary	hists = {u:fp.Get("h%d"%n) for n,u in enumerate("uvw")}
list	limits = [1,1,2]
list	lmins = [-3, -3, -3]
string	pdffile = "plot_impactzipper.pdf"
	c = ROOT.TCanvas()
	lminout = bilogify(h, lmin)
	title = h.GetTitle()

Function Documentation

def wirecell.gen.test.plot_impactzipper.bilogify	(		hist,
			lmin = None
	)

Modify histogram so that it is rescaled in a "bilog" manner.
Negative values are -log10(-z).  An offset in log10 is returned
and represents the 0.

Definition at line 6 of file plot_impactzipper.py.

def bilogify(hist, lmin = None):
    """
    Modify histogram so that it is rescaled in a "bilog" manner.
    Negative values are -log10(-z).  An offset in log10 is returned
    and represents the 0.
    """
    nx,ny = hist.GetNbinsX(), hist.GetNbinsY()
    zmax = max(abs(hist.GetMinimum()), abs(hist.GetMaximum()))

    if lmin is None:
        zmin = sum([abs(hist.GetBinContent(99,iy)) for iy in range(ny)]) / ny
        lmin = int(log10(zmin))

    lmax = 1+int(log10(zmax))

    for ix in range(hist.GetNbinsX()):
        for iy in range(hist.GetNbinsY()):
            val = hist.GetBinContent(ix, iy)
            if val == 0.0:
                hist.SetBinContent(ix, iy, 0.0)
                continue
            
            sign = 1.0
            if val > 0:                   # opposite sign to match Xin 
                sign = -1.0
            lval = log10(abs(val))
            if lval < lmin:
                lval = 0.0
            else:
                lval -= lmin
                lval *= sign
            hist.SetBinContent(ix, iy, lval)
            continue
        continue

    lhmax = lmax - lmin

    hist.SetMaximum(lhmax)
    hist.SetMinimum(-lhmax)
    return lmin

def wirecell.gen.test.plot_impactzipper.set_palette

(

which = "custom"

)

Definition at line 58 of file plot_impactzipper.py.

def set_palette(which = "custom"):
    if not which or which == "custom":
        ROOT.TColor.CreateGradientColorTable(len(stops), stops, reds, greens, blues, 100)
        return
    ROOT.gStyle.SetPalette(which)

Variable Documentation

wirecell.gen.test.plot_impactzipper.blues = array('d',[ 0.51, 1.00, 1.00, 0.00, 0.00 ])

Definition at line 52 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.c = ROOT.TCanvas()

Definition at line 71 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.fp = ROOT.TFile.Open("build/gen/test_impactzipper-uvw.root")

Definition at line 64 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.greens = array('d',[ 0.00, 0.81, 1.00, 0.20, 0.00 ])

Definition at line 51 of file plot_impactzipper.py.

dictionary wirecell.gen.test.plot_impactzipper.hists = {u:fp.Get("h%d"%n) for n,u in enumerate("uvw")}

Definition at line 65 of file plot_impactzipper.py.

list wirecell.gen.test.plot_impactzipper.limits = [1,1,2]

Definition at line 66 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.lminout = bilogify(h, lmin)

Definition at line 89 of file plot_impactzipper.py.

list wirecell.gen.test.plot_impactzipper.lmins = [-3, -3, -3]

Definition at line 67 of file plot_impactzipper.py.

string wirecell.gen.test.plot_impactzipper.pdffile = "plot_impactzipper.pdf"

Definition at line 69 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.reds = array('d',[ 0.00, 0.00, 1.00, 1.00, 0.51 ])

Definition at line 50 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.stops = array('d',[ 0.00, 0.45, 0.50, 0.55, 1.00 ])

Definition at line 49 of file plot_impactzipper.py.

wirecell.gen.test.plot_impactzipper.title = h.GetTitle()

Definition at line 90 of file plot_impactzipper.py.