Skip to content
Snippets Groups Projects

remove GhostProb cuts from standard particles

Merged Aniol Lobo Salvia requested to merge alobo_GhostProb_standardparticles into master
@@ -533,12 +533,13 @@ def filter_leptons_loose(particles,
pid_mu=None,
ismuon=False,
minipchi2=-999,
trghostprob=999):
trghostprob=None):
"""Returns loosely preselected leptons """
pvs = pv_maker()
code = F.require_all(F.PT > pt_min, F.P > p_min,
F.MINIPCHI2(pvs) > minipchi2,
F.GHOSTPROB < trghostprob)
F.MINIPCHI2(pvs) > minipchi2)
if trghostprob is not None:
code = code & (F.GHOSTPROB < trghostprob)
if lepton == 'muon':
if probnn_mu is not None:
code = (F.PROBNN_MU > probnn_mu) & code
@@ -570,7 +571,7 @@ def make_detached_dielectron(name='make_detached_dielectron_{hash}',
probnn_e=None,
pt_e=0.25 * GeV,
minipchi2=9.,
trghostprob=0.25,
trghostprob=None,
adocachi2cut=30,
bpvvdchi2=30,
vfaspfchi2ndof=10):
@@ -637,7 +638,7 @@ def make_detached_dielectron_with_brem(
pt_e=0.25 * GeV,
p_e=0. * GeV,
minipchi2=9.,
trghostprob=0.25,
trghostprob=None,
dielectron_ID="J/psi(1S)",
pt_diE=0 * MeV,
m_diE_min=0 * MeV,
@@ -653,8 +654,9 @@ def make_detached_dielectron_with_brem(
pvs = pv_maker()
particles = electron_maker()
code_e = F.require_all(F.PT > pt_e, F.P > p_e,
F.MINIPCHI2(pvs) > minipchi2,
F.GHOSTPROB < trghostprob)
F.MINIPCHI2(pvs) > minipchi2)
if trghostprob is not None:
code_e = code_e & (F.GHOSTPROB < trghostprob)
if PIDe_min is not None:
code_e = code_e & (F.PID_E > PIDe_min)
if probnn_e is not None:
@@ -687,7 +689,7 @@ def make_detached_mue(name='make_detached_mue_{hash}',
probnn_e=None,
pt_e=0.25 * GeV,
minipchi2=9.,
trghostprob=0.25,
trghostprob=None,
adocachi2cut=30,
bpvvdchi2=30,
vfaspfchi2ndof=10):
@@ -740,7 +742,7 @@ def make_detached_mue_with_brem(name='make_detached_mue_with_brem_{hash}',
probnn_e=None,
min_pt_e=0.25 * GeV,
minipchi2_track=9.,
max_trghostprob=0.25,
max_trghostprob=None,
max_adocachi2=30,
min_bpvvdchi2=30,
max_vchi2ndof=10):
@@ -793,7 +795,7 @@ def make_detached_mumu(name='make_detached_mumu_{hash}',
probnn_mu=None,
pt_mu=0. * GeV,
minipchi2=9.,
trghostprob=0.25,
trghostprob=None,
adocachi2cut=30,
bpvvdchi2=30,
vfaspfchi2ndof=10):
Loading