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()
'../../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()
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()
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()