#################
# Configuration #
#################
import ROOT
import numpy as np
import plotfactory as pf
from glob import glob
pf.setpfstyle()
tt = pf.makechain(False)
#file = ROOT.TFile('tree.root')
#tt = file.Get('tree')
ntries = tt.GetEntries()
print('Number of entries: ' + str(ntries))
#################
# Define x-axes #
#################
pTbins = np.arange(3., 65, 1)
dxybins = np.arange(0., 600, 10)
###################
# Create canvases #
###################
print('Preparing canvas')
t = ROOT.TCanvas('t', 'd')
t1 = ROOT.TCanvas('t1', 'd1')
t2 = ROOT.TCanvas('t2', 'd2')
t3 = ROOT.TCanvas('t3', 'd3')
import ROOT as rt import numpy as np import plotfactory as pf import ntup_dir as nt from glob import glob import sys from pdb import set_trace from os.path import normpath, basename #################################################################################################### outdir = '/eos/user/v/vstampf/ntuples/plots/id_efficiency/' indir = nt.getntupdir() #################################################################################################### ntdr = basename(normpath(indir)) #################################################################################################### t = pf.makechain(True) #################################################################################################### #################################################################################################### pf.setpfstyle() #################################################################################################### #################################################################################################### inc = 'is_in_acc == 1 & l0_matched_electron_pt > 30 & l1_matched_muon_pt > 4 & l2_matched_muon_pt > 4' loose = ' & l1_matched_muon_id_l & l2_matched_muon_id_l' globl = ' & l1_matched_muon_is_gl & l2_matched_muon_is_gl' emumu = ' & abs(l0_pdgId) == 11 & abs(l1_pdgId)==13 & abs(l2_pdgId) == 13 ' #################################################################################################### idgl = inc + globl + emumu idL = inc + loose + emumu denom = inc + emumu #################################################################################################### b_2disp = np.arange(0.,200, 5) ####################################################################################################
#################
# Configuration #
#################
import ROOT
import numpy as np
import plotfactory as pf
from glob import glob
pf.setpfstyle()
bewl = input("Makechain: True or False\n")
tt = pf.makechain(bewl)
#file = ROOT.TFile('tree.root')
#tt = file.Get('tree')
ntries = tt.GetEntries()
print('Number of entries: ' + str(ntries))
#################
# Define x-axes #
#################
pTbins = np.arange(5., 73, 5)
###################
# Create canvases #
###################
print('Preparing canvas')
t = ROOT.TCanvas('t', 'cap_d1')
#####################
# Create histograms #
fout = rt.TFile('histoseffpurmassd3.root', 'recreate')
#########################################
# Make Chain from selection of samples
#########################################
# Get the option from the command line, using 'True' as a fallback.
if len(sys.argv) > 1 and sys.argv[1] == 'test':
setting = False
print('Using a selection of samples')
else:
setting = True
print('Using all samples')
tt = pf.makechain(setting)
nentries = tt.GetEntries()
print('number of total entries in chain:\t\t\t%d' % (nentries))
#########################################
# Produce KPIs
#########################################
n_2OrMoreMuons = tt.GetEntries()
print('number of all events with 2 or more muons:\t\t%d' % (n_2OrMoreMuons))
n_reconstructable = tt.GetEntries(
'flag_hnl_reconstructable == 1 & abs(l1_pdgId) == 13 & abs(l1_eta) < 2.4 & abs(l2_pdgId) == 13 & abs(l2_eta) < 2.4'
)
print('number of events with reconstructable HNLs:\t\t%d' %
