Skip to content
Snippets Groups Projects

Add spruce lines for B2HPi0 with Dalitz pi0 decays

Open Tom Hadavizadeh requested to merge thadaviz/add_spruce_B2HPi0_Dalitz into bnoc_run3
Files
6
@@ -533,7 +533,8 @@ def make_b2hh(
pt_min=1200 * MeV,
dira_min=0.99,
ipchi2_max=9,
fdchi2_min=100,
fdchi2_min=None,
maxipchi2tracks_min=None,
AllowDiffInputsForSameIDChildren=False,
):
combination_cut = require_all(
@@ -541,6 +542,11 @@ def make_b2hh(
in_range(am_min, F.MASS, am_max),
F.MAXSDOCACHI2CUT(docachi2),
)
if maxipchi2tracks_min is not None:
combination_cut = require_all(
combination_cut, F.MAX(F.OWNPVIPCHI2) > maxipchi2tracks_min
)
composite_cut = require_all(
F.PT > pt_min,
F.OWNPVDIRA > dira_min,
@@ -1071,15 +1077,29 @@ def make_bds2ppbarppbar(
sdocachi2_max=20.0,
b_min_lft=0.2 * picosecond,
b_min_fdchi2=100,
MinimalPidProtonCut=0.01,
):
MinimalPid4Proton = (
F.CHILD(1, F.PROBNN_P)
* F.CHILD(2, F.PROBNN_P)
* F.CHILD(3, F.PROBNN_P)
* F.CHILD(4, F.PROBNN_P)
) > MinimalPidProtonCut
combination12_code = F.require_all(
F.MAXSDOCACHI2CUT(sdocachi2_max), F.MASS < am_max_mass
F.MAXSDOCACHI2CUT(sdocachi2_max),
F.MASS < am_max_mass,
)
combination123_code = F.MAXSDOCACHI2CUT(sdocachi2_max)
combination123_code = F.require_all(
F.MAXSDOCACHI2CUT(sdocachi2_max),
F.MASS < am_max_mass,
)
combination_code = F.require_all( ## form the ppbarppbar
in_range(am_min_mass, F.MASS, am_max_mass), F.MAXSDOCACHI2CUT(sdocachi2_max)
in_range(am_min_mass, F.MASS, am_max_mass),
F.MAXSDOCACHI2CUT(sdocachi2_max),
MinimalPid4Proton,
)
vertex_code = F.require_all( ## fit to the B0(s)->ppbarppbar
@@ -1120,11 +1140,13 @@ def make_bds2ppbar(
vtx_max_ipchi2=16,
b0_dira=0.9997,
b_min_lft=0.2 * picosecond,
MinimalPidProtonCut=0.01,
):
combination_code = F.require_all(
in_range(am_min_mass, F.MASS, am_max_mass),
F.MAX(F.PT) > comb_max_pt,
F.MAX(F.OWNPVIPCHI2) > am_max_ipchi2,
F.CHILD(1, F.PROBNN_P) * F.CHILD(2, F.PROBNN_P) > MinimalPidProtonCut,
)
vertex_code = F.require_all(
Loading