|
| string | plotBeamSpotVxVal.__author__ = 'Peter Loscutoff' |
| |
| string | plotBeamSpotVxVal.__version__ = '$Id $' |
| |
| string | plotBeamSpotVxVal.__usage__ = 'plotBeamSpotVxVal.py [options] bs.nt.root split.nt.root' |
| |
| | plotBeamSpotVxVal.parser = OptionParser(usage=__usage__, version=__version__) |
| |
| | plotBeamSpotVxVal.dest |
| |
| | plotBeamSpotVxVal.action |
| |
| | plotBeamSpotVxVal.default |
| |
| | plotBeamSpotVxVal.False |
| |
| | plotBeamSpotVxVal.help |
| |
| | plotBeamSpotVxVal.options |
| |
| | plotBeamSpotVxVal.args |
| |
| | plotBeamSpotVxVal.canvas = TCanvas("PlotBsVxVal","PlotBsVxVal",1000,500) |
| |
| | plotBeamSpotVxVal.bsfile = TFile(args[0]) |
| |
| | plotBeamSpotVxVal.splitfile = TFile(args[1]) |
| |
| dictionary | plotBeamSpotVxVal.subplotSVs = {} |
| |
| int | plotBeamSpotVxVal.fitbins = 3001 |
| |
| int | plotBeamSpotVxVal.plotrange = int(options.plotrange) + 1 |
| |
| int | plotBeamSpotVxVal.outlierCut = 20 |
| |
| | plotBeamSpotVxVal.plotSV = TH1D("p(s)SV","Pull vs # Tracks",plotrange*2,-.25,plotrange-.25) |
| |
| string | plotBeamSpotVxVal.nameSV = "plotSV%s" % i |
| |
| | plotBeamSpotVxVal.tree = gDirectory.Get("splitVertex") |
| |
| | plotBeamSpotVxVal.entries = tree.GetEntriesFast() |
| |
| | plotBeamSpotVxVal.noprint = tree.GetEntry(entry) |
| |
| tuple | plotBeamSpotVxVal.binSV = (tree.tracks_odd+tree.tracks_even) |
| |
| dictionary | plotBeamSpotVxVal.fitsigma = subplotSVs[nameSV].GetRMS() |
| |
| tuple | plotBeamSpotVxVal.fserror = (subplotSVs[nameSV].GetRMS())/sqrt(subplotSVs[nameSV].Integral()) |
| |
| int | plotBeamSpotVxVal.bin = i+1 |
| |
| | plotBeamSpotVxVal.legend = TLegend(.60,.75,.92,.92) |
| |
| dictionary | plotBeamSpotVxVal.subplots = {} |
| |
| dictionary | plotBeamSpotVxVal.bstable = {} |
| |
| list | plotBeamSpotVxVal.splitWeight = [] |
| |
| list | plotBeamSpotVxVal.fullWeight = [] |
| |
| | plotBeamSpotVxVal.lumiedge = array('d') |
| |
| int | plotBeamSpotVxVal.lumibins = 0 |
| |
| int | plotBeamSpotVxVal.lumiend = 0 |
| |
| | plotBeamSpotVxVal.bstree = gDirectory.Get("BeamSpotNt") |
| |
| | plotBeamSpotVxVal.lumibegin = bstree.lbStart |
| |
| string | plotBeamSpotVxVal.name = str(lumibegin) + "to" + str(lumiend) |
| |
| | plotBeamSpotVxVal.plotX = TH1D("p(s)X","k vs Lumi Block",lumibins,lumiedge) |
| |
| | plotBeamSpotVxVal.plotBS = TH1D("p(s)bs","k vs Lumi Block",lumibins,lumiedge) |
| |
| | plotBeamSpotVxVal.sptree = gDirectory.Get("splitVertex") |
| |
| | plotBeamSpotVxVal.splitBin = int(.5*(sptree.tracks_odd + sptree.tracks_even) ) |
| |
| | plotBeamSpotVxVal.vxtree = gDirectory.Get('Vertices') |
| |
| | plotBeamSpotVxVal.fullBin = vxtree.nTracks |
| |
| | plotBeamSpotVxVal.lumi = sptree.lumi |
| |
| dictionary | plotBeamSpotVxVal.bs_x0 = bstable[lumi]["x0"] |
| |
| dictionary | plotBeamSpotVxVal.bs_y0 = bstable[lumi]["y0"] |
| |
| dictionary | plotBeamSpotVxVal.bs_a = bstable[lumi]["ax"] |
| |
| dictionary | plotBeamSpotVxVal.bs_b = bstable[lumi]["ay"] |
| |
| dictionary | plotBeamSpotVxVal.bs_z = bstable[lumi]["z"] |
| |
| dictionary | plotBeamSpotVxVal.bs_sx = bstable[lumi]["sx"] |
| |
| dictionary | plotBeamSpotVxVal.bs_sy = bstable[lumi]["sy"] |
| |
| dictionary | plotBeamSpotVxVal.bs_sz = bstable[lumi]["sz"] |
| |
| dictionary | plotBeamSpotVxVal.bs_rho = bstable[lumi]["rhoxy"] |
| |
| | plotBeamSpotVxVal.range = k.split("to") |
| |
| | plotBeamSpotVxVal.cov = array('d') |
| |
| | plotBeamSpotVxVal.covmat = TMatrixDSym(2) |
| |
| int | plotBeamSpotVxVal.min = .8 |
| |
| float | plotBeamSpotVxVal.max = 1.2 |
| |
| int | plotBeamSpotVxVal.sumw = 0 |
| |
| int | plotBeamSpotVxVal.sumwx = 0 |
| |
| int | plotBeamSpotVxVal.midpoint = int(name.split("to")[0])+.5 |
| |
| | plotBeamSpotVxVal.globalbin = plotBS.FindBin(midpoint) |
| |
| int | plotBeamSpotVxVal.vsavg = sumwx/sumw |
| |
| float | plotBeamSpotVxVal.vserr = 1.0/sumw |
| |
| int | plotBeamSpotVxVal.bsavg = sumwx/sumw |
| |
| float | plotBeamSpotVxVal.bserr = 1.0/sumw |
| |
1.8.18