Skip to content
Snippets Groups Projects

Draft: Update diphoton.py

Open Alejandro Novo Cal requested to merge diphotonHlt2_try into qee_upgrade
1 unresolved thread
1 file
+ 10
10
Compare changes
  • Side-by-side
  • Inline
@@ -33,6 +33,7 @@ from Hlt2Conf.standard_particles import make_photons
all_lines = {}
@acasaisv
@configurable
def filter_qee_photons(cl_min,
@@ -41,19 +42,18 @@ def filter_qee_photons(cl_min,
name='filter_qee_photons_{hash}',
make_particles=make_photons):
"""Photon filter, currently missing CL functor"""
code = F.require_all(F.PT > et_min, F.P > e_min)
code = F.require_all(F.PT > et_min, F.P > e_min, F.IS_NOT_H > cl_min)
return ParticleFilter(make_particles(), F.FILTER(code), name=name)
@configurable
def make_diphoton(
min_cl_gamma=0.3, # CL is not implemented currently
min_pt_gamma=5 * GeV,
min_p_gamma=16. * GeV,
min_sumpt=12. * GeV,
min_B_pt=7. * GeV,
min_B_mass=2. * GeV,
min_cl_gamma=0.3, # CL is now implemented
min_pt_gamma=3 * GeV,
min_p_gamma=5. * GeV,
min_sumpt=6. * GeV,
min_B_pt=3. * GeV,
min_B_mass=4.3 * GeV,
min_pt_asym=0.1):
"""
Return Diphotons from B/ALPs around a few GeV.
@@ -62,7 +62,7 @@ def make_diphoton(
descriptor = 'B_s0 -> gamma gamma'
photons = filter_qee_photons(
cl_min=min_cl_gamma, et_min=min_pt_gamma,
e_min=min_p_gamma) #CL is not implemented at the moment
e_min=min_p_gamma) #CL is now implemented
sum_pt = F.CHILD(1, F.PT) + F.CHILD(2, F.PT)
pt_asym = (F.CHILD(1, F.PT) - F.CHILD(2, F.PT)) / (
Loading