Newer
Older
STAging / scripts / checkTrackSelection.py
@Elena Graverini Elena Graverini on 3 May 2017 3 KB [script] Draw the track quality cuts
#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Author: Elena Graverini
# @Date:   2015-10-27 18:26:49
# @Last Modified by:   Elena Graverini
# @Last Modified time: 2017-05-03 10:57:13

import ROOT as r
import os

location = {'TT': "/disk/data1/hep/elena/data/ST/Aging_Tuples/TT/TTaU/",
            'IT': "/disk/data1/hep/elena/data/ST/Aging_Tuples/IT/T3X2/"}
out_location = os.getcwd() + '/trackSel'
os.system('mkdir -p %s' % out_location)

fills = [2797, 3108, 3478, 3960, 4518]
dets = ['TT', 'IT']
val = {'TT': 5, 'IT': 7}
stat = 10000000

r.gROOT.ProcessLine('.X %s/../include/lhcbstyle.C' % os.getcwd())
#r.gStyle.SetPadRightMargin(0.16)
#r.gStyle.SetPadLeftMargin(r.gStyle.GetPadLeftMargin() * 0.7)
r.gStyle.SetPadLeftMargin(0.20)
r.gStyle.SetPadRightMargin(0.25)
r.gStyle.SetPadBottomMargin(0.15)
r.gStyle.SetPadTopMargin(0.05)


def printSelection(det, fill, stat):
    f = r.TFile(location[det] + '%s.root' % fill, 'read')
    t_sig = f.Get('STADCTrackMonitor/HitInfo/TTaU')
    # t_bg = f.Get('STADCTrackMonitor/NoiseInfo/TTaU')
    # t_sig.Draw("GhostP : TrChi2/TrNDoF>>hnew(100,0., 5,100, 0.,1.)", "", "colz", 10000, 0)

    t_sig.Draw("GhostP : TrChi2/TrNDoF>>hsig(100,0., 5,100, 0.,1.)", "val%s>12" % val[det], "colz", stat, 0)
    h_sig = r.gPad.GetPrimitive('hsig')
    t_sig.Draw("GhostP : TrChi2/TrNDoF>>hnoise(100,0., 5,100, 0.,1.)", "val%s<12" % val[det], "colz", stat, 0)
    h_noise = r.gPad.GetPrimitive('hnoise')

    c = r.TCanvas()
    h_sig.Divide(h_noise)
    h_sig.Draw("colz")
    h_sig.GetZaxis().SetRangeUser(0, 10)
    h_sig.SetTitle('TTaU S/B ratio for fill %s' % fill)
    h_sig.GetXaxis().SetTitle(r'#chi^{2}_{track}/ndf')
    h_sig.GetYaxis().SetTitle(r'Ghost Prob.')
    h_sig.GetZaxis().SetTitle(r'S/B ratio')
    h_sig.Draw("colz")

    c, line = putLine(c, det)

    label = r.TPaveText( 0.74 - r.gStyle.GetPadRightMargin(), 0.87 - r.gStyle.GetPadTopMargin(),
                         0.95 - r.gStyle.GetPadRightMargin(), 0.95 - r.gStyle.GetPadTopMargin(), "BRNDC")
    label.SetFillColor(0)
    label.SetTextAlign(12)
    label.SetBorderSize(0)
    label.SetTextFont(132)
    # label.SetTextFont(22)
    label.SetTextSize(0.06)
    label.SetTextAlign(22)
    # label.SetTextColor(r.kWhite)
    label.AddText('Fill %s' % fill)
    label.Draw('same')
    c.Modified()
    c.Update()

    c.SaveAs('%s/%s_fill%s.pdf' % (out_location, det, fill))
    f.Close()


def putLine(c, det):
    line = r.TLine()
    line.SetLineWidth(4)
    line.SetLineColor(kWhite)
    c.cd()
    if 'TT' in det:
        line.DrawLine(0.0, 0.01, 1.6, 0.16)
        line.DrawLine(3.0, 0.07, 1.6, 0.16)
        line.DrawLine(3.0, 0.07, 3.0, 0.0)
    else:
        line.DrawLine(0.0, 0.0, 3.5, 0.35)
        line.DrawLine(3.5, 0.35, 5.0, 0.20)
    return c, line


def printByStep(fill, stat):
    f = r.TFile(location_tt + '%s.root' % fill, 'read')
    t_sig = f.Get('STADCTrackMonitor/HitInfo/TTaU')

    # Look at Landau distribution
    for step in range(67):
        c = r.TCanvas()
        t_sig.Draw('val5>>hnew(70,-30,90)', 'odinStep==%s' % step)
        c.SaveAs()


    draw_sig = "GhostP : TrChi2/TrNDoF>>hsig(100, 0., 10., 100, 0., 1.)"
    draw_bg = "GhostP : TrChi2/TrNDoF>>hnoise(100, 0., 10., 100, 0., 1.)"
    f.Close()


if __name__ == '__main__':
    for det in dets:
        for fill in fills:
            printSelection(det, fill, stat)