np.seterr(divide='ignore', invalid='ignore', over='ignore')
parser = argparse.ArgumentParser(description="leptonjet ABCD")
parser.add_argument(
"--sync",
action='store_true',
help="issue rsync command to sync plots folder to lxplus web server")
args = parser.parse_args()
dml = DatasetMapLoader()
# dataDS, dataMAP = dml.fetch('data')
bkgDS, bkgMAP, bkgSCALE = dml.fetch('bkg')
sdml = SigDatasetMapLoader()
sigDS_2mu2e, sigSCALE_2mu2e = sdml.fetch('2mu2e')
sigDS_4mu, sigSCALE_4mu = sdml.fetch('4mu')
class LjABCDProcessor(processor.ProcessorABC):
def __init__(self, data_type='bkg'):
self.data_type = data_type
dataset_axis = hist.Cat('dataset', 'dataset')
iso_axis = hist.Bin('iso', 'min pfIso', 50, 0, 0.5)
bin_axis = hist.Bin('val', 'binary value', 3, 0, 3)
self._accumulator = processor.dict_accumulator({
'chan-4mu':
hist.Hist('Counts', dataset_axis, iso_axis, bin_axis),
'chan-2mu2e':
hist.Hist('Counts', dataset_axis, iso_axis, bin_axis),
parser = argparse.ArgumentParser(
description="leptonjet isolation profile of pt, eta")
parser.add_argument(
"--sync",
action='store_true',
help="issue rsync command to sync plots folder to lxplus web server")
parser.add_argument("--preserve",
action='store_true',
help="preserve plots in ROOT file")
args = parser.parse_args()
import ROOT
ROOT.gROOT.SetBatch()
sdml = SigDatasetMapLoader()
sigDS, sigSCALE = sdml.fetch()
dml = DatasetMapLoader()
bkgDS, bkgMAP, bkgSCALE = dml.fetch('bkg')
dataDS, dataMAP = dml.fetch('data')
@contextmanager
def _setIgnoreLevel(level):
originalLevel = ROOT.gErrorIgnoreLevel
ROOT.gErrorIgnoreLevel = level
yield
ROOT.gErrorIgnoreLevel = originalLevel
"""Leptonjet isolation, pt, """
description="find isolation cut for EGM-type leptonjet")
parser.add_argument(
"--sync",
action='store_true',
help="issue rsync command to sync plots folder to lxplus web server")
parser.add_argument("--preserve",
action='store_true',
help="preserve plots in ROOT file")
args = parser.parse_args()
import ROOT
ROOT.gROOT.SetBatch()
sdml = SigDatasetMapLoader()
sigDS, sigSCALE = sdml.fetch('2mu2e')
dml = DatasetMapLoader()
bkgDS, bkgMAP, bkgSCALE = dml.fetch('bkg')
dataDS, dataMAP = dml.fetch('data')
@contextmanager
def _setIgnoreLevel(level):
originalLevel = ROOT.gErrorIgnoreLevel
ROOT.gErrorIgnoreLevel = level
yield
ROOT.gErrorIgnoreLevel = originalLevel
def root_filling(output, datatype):
from FireHydrant.Tools.metfilter import MetFilters
from FireHydrant.Tools.trigger import Triggers
np.seterr(divide='ignore', invalid='ignore', over='ignore')
plt.rcParams["savefig.dpi"] = 120
plt.rcParams["savefig.bbox"] = "tight"
parser = argparse.ArgumentParser(description="[AN] leptonjet pair delta phi")
parser.add_argument(
"--sync",
action='store_true',
help="issue rsync command to sync plots folder to lxplus web server")
args = parser.parse_args()
sdml = SigDatasetMapLoader()
sigDS, sigSCALE = sdml.fetch('all')
class MuLJVtxProcessor(processor.ProcessorABC):
def __init__(self):
dataset_axis = hist.Cat('dataset', 'dataset')
reco_axis = hist.Cat('reco', 'reco type')
vxy_axis = hist.Bin('vxy', 'vxy [cm]', [0, 5, 70, 700])
lxy_axis = hist.Bin('lxy', 'lxy [cm]', 100, 0, 700)
bool_axis = hist.Bin('boolean', 'true/false', 2, 0, 2)
reso_axis = hist.Bin('reso', '(vxy(reco)-vxy(gen)) [cm]', 100, -25, 25)
error_axis = hist.Bin('error', '$\sigma_{lxy}$', 100, 0, 100)
cos_axis = hist.Bin('cos', r'$|cos(\theta)|$', 50, -1, 1)
sig_axis = hist.Bin('sig', 'lxy/$\sigma_{lxy}$', 50, 0, 50)
self._accumulator = processor.dict_accumulator({
'vertexgood':
plt.rcParams["savefig.bbox"] = "tight"
parser = argparse.ArgumentParser(
description="produce 2D leptonjet isolation map")
parser.add_argument(
"--sync",
action='store_true',
help="issue rsync command to sync plots folder to lxplus web server")
args = parser.parse_args()
dml = DatasetMapLoader()
bkgDS, bkgMAP, bkgSCALE = dml.fetch('bkg')
dataDS, dataMAP = dml.fetch('data')
sdml = SigDatasetMapLoader()
sigDS, sigSCALE = sdml.fetch('simple')
"""Leptonjet Isolation."""
class LeptonJetIsoProcessor(processor.ProcessorABC):
def __init__(self, dphi_control=False, data_type='bkg'):
self.dphi_control = dphi_control
self.data_type = data_type
dataset_axis = hist.Cat('dataset', 'dataset')
lj0iso_axis = hist.Bin('lj0iso', 'iso value', 20, 0, 1)
lj1iso_axis = hist.Bin('lj1iso', 'iso value', 20, 0, 1)
channel_axis = hist.Cat('channel', 'channel')
type_axis = hist.Cat('isotype', 'isotype')
njet_axis = hist.Bin('njet', 'njet', 10, 0, 10)
self._accumulator = processor.dict_accumulator({