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()
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., 500., 1000.]
bins = [0., 10., 25., 32.5, 40., 50., 60., 90., 200.]
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,
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., 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,
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= 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='m_{#mu#tau}^{vis}' , ylabel='L1 + HLT #tau efficiency', bins=nbins, bini= 0. , bine= 100. ),
]
HLTPlotter = Plotter(variables=variables,
out_filename='diTau_plots_good.root',
sel_baseline=baseline,
sel_extra_den=[eta_bins, mt, zmass, sign, iso],
sel_num=[TriggerSelection])
HLTPlotter.run()