Python Plotter Beispiele

Beispiel #1

Datei: plotter_tau_leg_of_di_tau_mc.py Projekt: rmanzoni/triggertools

for fname in filenames:
    if fname.endswith('.url'):
        with open(fname) as ff:
            fname = ff.readlines()[0].rstrip()
    t1.Add(fname)

nbins   = 40
bins    = [0., 10., 20., 25., 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200., 500., 1000.]
# bins    = [0., 10., 20., 22.5, 25., 27.5, 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200., 500., 1000.]

variables = [
#     Efficiency1D(tree=t1, name='bx'     , variable='bx'        , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='# bunch crossing'  , ylabel='L1 + HLT #tau efficiency', bins=125  , bini=0     , bine=3500   ),
#     Efficiency1D(tree=t1, name='run'    , variable='run'       , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='run'               , ylabel='L1 + HLT #tau efficiency', bins=1400 , bini=273100, bine=274500 ),
#     Efficiency1D(tree=t1, name='tau_eta', variable='l2_eta'    , histo_name='tau_eta', cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.   , bine=   3.  ),
#     Efficiency1D(tree=t1, name='tau_phi', variable='l2_phi'    , histo_name='tau_phi', cut_num='1', cut_den='1', xlabel='offline #tau #phi' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.15 , bine=   3.15),
#     Efficiency1D(tree=t1, name='npv'    , variable='n_vertices', histo_name='npv'    , cut_num='1', cut_den='1', xlabel='# PV'              , ylabel='L1 + HLT #tau efficiency', bins=10   , bini= 0    , bine=  60   ),
    Efficiency1D(tree=t1, name='tau_pt'      , variable='l2_pt'            , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline #tau p_{T}'       , ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
    Efficiency1D(tree=t1, name='tau_pt_up'   , variable='1.03 * l2_pt'     , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline 1.03 * #tau p_{T}', ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
    Efficiency1D(tree=t1, name='tau_pt_down' , variable='0.97 * l2_pt'     , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline 0.97 * #tau p_{T}', ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
#     Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis'   , cut_num='1', cut_den='1', xlabel='m_{#mu#tau}^{vis}' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0.   , bine= 100.  ),
]

HLTPlotter = Plotter(variables     = variables                         , 
                     #out_filename  = 'tauleg_of_di_tau_real_taus_skim_mc_v2.root'     , 
                     out_filename  = 'fuffa.root'     , 
                     sel_baseline  = baseline                          , 
                     sel_extra_den = [ weight, eta_bins, mt, zmass, sign, iso, decaymode ], 
                     sel_num       = [ TriggerSelection ]              )

HLTPlotter.run()

Beispiel #2

        tree=t1,
        name='tau_pt',
        variable='l2_pt',
        histo_name='tau_pt',
        cut_num='1',
        cut_den='1',
        xlabel='offline #tau p_{T}',
        ylabel='L1 + HLT #tau efficiency',
        bins=bins,
    ),
    Efficiency1D(tree=t1,
                 name='mvis',
                 variable='mvis',
                 histo_name='mvis',
                 cut_num='1',
                 cut_den='1',
                 xlabel='m_{#mu#tau}^{vis}',
                 ylabel='L1 + HLT #tau efficiency',
                 bins=nbins,
                 bini=0.,
                 bine=100.),
]

HLTPlotter = Plotter(variables=variables,
                     out_filename='mu_tau_singleL1_plots_ss.root',
                     sel_baseline=baseline,
                     sel_extra_den=[eta_bins, mt, zmass, sign, iso],
                     sel_num=[TriggerSelection])

HLTPlotter.run()

Beispiel #3

Datei: plotter_tau_leg_of_lepton_plus_tau_mc.py Projekt: rmanzoni/triggertools

for fname in filenames:
    if fname.endswith('.url'):
        with open(fname) as ff:
            fname = ff.readlines()[0].rstrip()
    t1.Add(fname)

nbins   = 40
# bins    = [0., 10., 20., 25., 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200., 500., 1000.]
bins    = [0., 10., 20., 22.5, 25., 27.5, 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200., 500., 1000.]

variables = [
#     Efficiency1D(tree=t1, name='bx'     , variable='bx'        , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='# bunch crossing'  , ylabel='L1 + HLT #tau efficiency', bins=125  , bini=0     , bine=3500   ),
#     Efficiency1D(tree=t1, name='run'    , variable='run'       , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='run'               , ylabel='L1 + HLT #tau efficiency', bins=1400 , bini=273100, bine=274500 ),
#     Efficiency1D(tree=t1, name='tau_eta', variable='l2_eta'    , histo_name='tau_eta', cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.   , bine=   3.  ),
#     Efficiency1D(tree=t1, name='tau_phi', variable='l2_phi'    , histo_name='tau_phi', cut_num='1', cut_den='1', xlabel='offline #tau #phi' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.15 , bine=   3.15),
#     Efficiency1D(tree=t1, name='npv'    , variable='n_vertices', histo_name='npv'    , cut_num='1', cut_den='1', xlabel='# PV'              , ylabel='L1 + HLT #tau efficiency', bins=10   , bini= 0    , bine=  60   ),
    Efficiency1D(tree=t1, name='tau_pt'      , variable='l2_pt'            , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline #tau p_{T}'       , ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
    Efficiency1D(tree=t1, name='tau_pt_up'   , variable='1.03 * l2_pt'     , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline 1.03 * #tau p_{T}', ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
    Efficiency1D(tree=t1, name='tau_pt_down' , variable='0.97 * l2_pt'     , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline 0.97 * #tau p_{T}', ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
#     Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis'   , cut_num='1', cut_den='1', xlabel='m_{#mu#tau}^{vis}' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0.   , bine= 100.  ),
]

HLTPlotter = Plotter(variables     = variables                         , 
                     out_filename  = 'tauleg_of_lepton_plus_tau_real_taus_skim_mc_v2.root'     , 
                     sel_baseline  = baseline                          , 
                     sel_extra_den = [ weight, eta_bins, mt, zmass, sign, iso, decaymode ], 
                     sel_num       = [ TriggerSelection ]              )

HLTPlotter.run()

Beispiel #4

Datei: plotter_tau_met_v4.py Projekt: tstreble/TauTriggerValidation

                 variable='n_vertices',
                 histo_name='npv',
                 cut_num='1',
                 cut_den='1',
                 xlabel='# PV',
                 ylabel='L1 + HLT #tau efficiency',
                 bins=10,
                 bini=0,
                 bine=30),
    Efficiency1D(
        tree=t1,
        name='tau_pt',
        variable='l2_pt',
        histo_name='tau_pt',
        cut_num='1',
        cut_den='1',
        xlabel='offline #tau p_{T}',
        ylabel='L1 + HLT #tau efficiency',
        bins=bins,
    ),
    #     Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis'   , cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0.   , bine= 100.  ),
]

HLTPlotter = Plotter(variables=variables,
                     out_filename='tau_met_plots.root',
                     sel_baseline=baseline,
                     sel_extra_den=[eta_bins, mt, zmass, sign],
                     sel_num=[TriggerSelection])

HLTPlotter.run()

Beispiel #5

Datei: L1_plotter.py Projekt: cms-tau-pog/TauTriggerValidation

    '../../SingleMuon_Run2016B_PromptReco_v2/H2TauTauTreeProducerTauMu/tree.root',    
]

t1 = ROOT.TChain('tree')

for fname in filenames:
    t1.Add(fname)

nbins   = 40
bins    = [0., 10., 20., 25., 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200.]

selections = ' (%s) * (%s) * (%s) * (%s) ' %(baseline, subtraction, eta_bins[''], decaymode[''])

variables = [
    Efficiency1D(tree=t1, name='tau_eta', variable='l2_eta'    , histo_name='tau_eta_iso', cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='iso L1 #tau efficiency', bins=nbins, bini=-3.  , bine=3.  ),
    Efficiency1D(tree=t1, name='tau_phi', variable='l2_phi'    , histo_name='tau_phi_iso', cut_num='1', cut_den='1', xlabel='offline #tau #phi' , ylabel='iso L1 #tau efficiency', bins=nbins, bini=-3.15, bine=3.15),
    Efficiency1D(tree=t1, name='npv'    , variable='n_vertices', histo_name='npv_iso'    , cut_num='1', cut_den='1', xlabel='# PV'              , ylabel='iso L1 #tau efficiency', bins=10   , bini= 0   , bine=30  ),
#     Efficiency1D(tree=t1, name='L1_pt'  , variable='l2_L1_pt'  , histo_name='L1_pt_iso'  , cut_num='1', cut_den='1', xlabel='L1 #tau E_{T}'     , ylabel='iso L1 #tau efficiency', bins=bins ,                      ),
    Efficiency1D(tree=t1, name='tau_pt' , variable='l2_pt'     , histo_name='tau_pt_iso' , cut_num='1', cut_den='1', xlabel='offline #tau p_{T}', ylabel='iso L1 #tau efficiency', bins=bins ,                      ),
    Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis_iso'   , cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='iso L1 #tau efficiency', bins=nbins, bini= 0.  , bine=100.),
]

L1Plotter = Plotter(variables     = variables                   , 
                    out_filename  = 'L1_plots.root'             , 
                    sel_baseline  = baseline                    , 
                    sel_extra_den = [ eta_bins, mt, zmass, sign], 
                    sel_num       = [ L1selection, L1pt ]       )

L1Plotter.run()

Beispiel #6

Datei: plotter_tau_leg_of_di_tau.py Projekt: rmanzoni/triggertools

variables = [
    #     Efficiency1D(tree=t1, name='bx'     , variable='bx'        , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='# bunch crossing'  , ylabel='L1 + HLT #tau efficiency', bins=125  , bini=0     , bine=3500   ),
    #     Efficiency1D(tree=t1, name='run'    , variable='run'       , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='run'               , ylabel='L1 + HLT #tau efficiency', bins=1400 , bini=273100, bine=274500 ),
    #     Efficiency1D(tree=t1, name='tau_eta', variable='l2_eta'    , histo_name='tau_eta', cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.   , bine=   3.  ),
    #     Efficiency1D(tree=t1, name='tau_phi', variable='l2_phi'    , histo_name='tau_phi', cut_num='1', cut_den='1', xlabel='offline #tau #phi' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.15 , bine=   3.15),
    #     Efficiency1D(tree=t1, name='npv'    , variable='n_vertices', histo_name='npv'    , cut_num='1', cut_den='1', xlabel='# PV'              , ylabel='L1 + HLT #tau efficiency', bins=10   , bini= 0    , bine=  60   ),
    Efficiency1D(
        tree=t1,
        name='tau_pt',
        variable='l2_pt',
        histo_name='tau_pt',
        cut_num='1',
        cut_den='1',
        xlabel='offline #tau p_{T}',
        ylabel='L1 + HLT #tau efficiency',
        bins=bins,
    ),
    #     Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis'   , cut_num='1', cut_den='1', xlabel='m_{#mu#tau}^{vis}' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0.   , bine= 100.  ),
]

HLTPlotter = Plotter(
    variables=variables,
    out_filename='tauleg_of_di_tau_real_taus_skim_v2.root',
    sel_baseline=baseline,
    sel_extra_den=[period, eta_bins, mt, zmass, sign, iso, decaymode],
    #                      sel_extra_den = [ period, previous_level, eta_bins, mt, zmass, sign, iso, pttau, decaymode ],
    sel_num=[TriggerSelection])

HLTPlotter.run()

Beispiel #7

Datei: L1_plotter.py Projekt: tstreble/TauTriggerValidation

        tree=t1,
        name='tau_pt',
        variable='l2_pt',
        histo_name='tau_pt_iso',
        cut_num='1',
        cut_den='1',
        xlabel='offline #tau p_{T}',
        ylabel='iso L1 #tau efficiency',
        bins=bins,
    ),
    Efficiency1D(tree=t1,
                 name='mvis',
                 variable='mvis',
                 histo_name='mvis_iso',
                 cut_num='1',
                 cut_den='1',
                 xlabel='offline #tau #eta',
                 ylabel='iso L1 #tau efficiency',
                 bins=nbins,
                 bini=0.,
                 bine=100.),
]

L1Plotter = Plotter(variables=variables,
                    out_filename='L1_plots.root',
                    sel_baseline=baseline,
                    sel_extra_den=[eta_bins, mt, zmass, sign],
                    sel_num=[L1selection, L1pt])

L1Plotter.run()

Beispiel #8

Datei: plotter.py Projekt: rmanzoni/triggertools

for fname in filenames:
    if fname.endswith('.url'):
        with open(fname) as ff:
            fname = ff.readlines()[0].rstrip()
    t1.Add(fname)

import pdb ; pdb.set_trace()

nbins   = 40
bins    = [0., 10., 20., 25., 30., 32.5, 35., 37.5, 40., 42.5, 45., 50., 55., 60., 70., 90., 120., 200., 500., 1000.]

variables = [
#     Efficiency1D(tree=t1, name='bx'     , variable='bx'        , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='# bunch crossing'  , ylabel='L1 + HLT #tau efficiency', bins=125  , bini=0     , bine=3500   ),
#     Efficiency1D(tree=t1, name='run'    , variable='run'       , histo_name='bx'     , cut_num='1', cut_den='1', xlabel='run'               , ylabel='L1 + HLT #tau efficiency', bins=1400 , bini=273100, bine=274500 ),
#     Efficiency1D(tree=t1, name='tau_eta', variable='l2_eta'    , histo_name='tau_eta', cut_num='1', cut_den='1', xlabel='offline #tau #eta' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.   , bine=   3.  ),
#     Efficiency1D(tree=t1, name='tau_phi', variable='l2_phi'    , histo_name='tau_phi', cut_num='1', cut_den='1', xlabel='offline #tau #phi' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini=-3.15 , bine=   3.15),
#     Efficiency1D(tree=t1, name='npv'    , variable='n_vertices', histo_name='npv'    , cut_num='1', cut_den='1', xlabel='# PV'              , ylabel='L1 + HLT #tau efficiency', bins=10   , bini= 0    , bine=  60   ),
    Efficiency1D(tree=t1, name='tau_pt' , variable='l2_pt'     , histo_name='tau_pt' , cut_num='1', cut_den='1', xlabel='offline #tau p_{T}', ylabel='L1 + HLT #tau efficiency', bins=bins ,                         ),
#     Efficiency1D(tree=t1, name='mvis'   , variable='mvis'      , histo_name='mvis'   , cut_num='1', cut_den='1', xlabel='m_{#mu#tau}^{vis}' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0.   , bine= 100.  ),
]

HLTPlotter = Plotter(variables     = variables                         , 
#                      out_filename  = 'mediumiso_plots_v3.root'     , 
                     out_filename  = 'mediumiso_plots_decaymode_and_eta.root'     , 
                     sel_baseline  = baseline                          , 
                     sel_extra_den = [ previous_level, eta_bins, mt, zmass, sign, iso, pttau, decaymode ], 
                     sel_num       = [ TriggerSelection ]              )

HLTPlotter.run()