Skip to content

Instantly share code, notes, and snippets.

@mattearnshaw
Last active August 29, 2015 14:14
Show Gist options
  • Select an option

  • Save mattearnshaw/6f1525d580e7fc966fac to your computer and use it in GitHub Desktop.

Select an option

Save mattearnshaw/6f1525d580e7fc966fac to your computer and use it in GitHub Desktop.
Parkinson Cuts
/* Contributions from B0s→φµ+µ−decays
Candidates are rejected if
- 5321<mK(π→K)µµ<5411MeV/c2 and either:
- 1010<mK(π→K)<1030MeV/c2 and the pion satisfies DLLKπ>−10;
- or 1030<mK(π→K)<1075MeV/c2 and the pion satisfies DLLKπ>10. */
RooFormulaVar B0stoPhimm("B0stoPhimm", "!((5321<@0<5411) && "
"((1010<@1<1030 && @2>-10) ||"
"(1030<@1<1075 && PI_PIDK>10)))",
RooArgList(B0_MM, B0_kpiswapMass, Pi_PIDK));
/* Contributions from B0→K∗0µ+µ−decays
Candidates are rejected if the kaon and pion satisfy
kaon DLL_Kπ + 10 < pion DLL_Kπ and 792 < m(K↔π)< 992MeV/c2.*/
RooFormulaVar B0toKstmm("B0toKstmm", "!(((@0 + 10) < @1) && (792<@2<992))",
RooArgList(K_PIDK, Pi_PIDK, B0_kpiswapMass));
/* Contributions from B0→J/ψK∗0 decays pp. 73
Candidates are rejected if
3036<m(π→µ)µ<3156MeV/c2 and the pion satisfies either
the IsMuon criteria or DLLµπ > 5.0;
or if 3036<m_(K→µ)_µ<3156MeV/c2 and the kaon satisfies either
the IsMuon criteria or DLL_µπ> 5.0. */
RooFormulaVar B0toJpsiK_1("B0toJpsiK_1",
"!((3036<@0<3156) && (@1 || @2 > 5.0))",
RooArgList(B0_pimuMass, Pi_isMuon, Pi_PIDmu));
RooFormulaVar B0toJpsiK_2("B0toJpsiK_2",
"!((3036<@0<3156) && (@1 || @2 > 5.0))",
RooArgList(B0_kmuMass, K_isMuon, Pi_PIDmu));
/* Contributions from Λ0b→Λ∗0(1520)µ + µ−decays pp. 73
Candidates are rejected if 5575 < mK(π→p)µµ< 5665MeV/c 2,
1490 < mK(π→p)<1550MeV/c2and the pion satisfies DLLpπ> 15.
Candidates are rejected if
5575<m(K→p)(π→K)µµ<5665MeV/c2, 1500<m(K→p)(π→K)<1550MeV/c2
and the pion satisfies DLLKπ> 10. */
RooFormulaVar L0btoLs0mm_1("L0btoLs0mm_1",
"!((5575<@0<5665) && (1490<@1<1550) && (@2>15))",
RooArgList(B0_MM, B0_Lambda_M_pK, Pi_PIDK));
RooFormulaVar L0btoLs0mm_2("L0btoLs0mm_2",
"!((5575<@0<5665) && (1500<@1<1550) && (@2>10))",
RooArgList(B0_MM, B0_Lambda_M_pK2, Pi_PIDK));
/* Contributions from B+→K+µ+µ−decays
candidates are removed by requiring mKπµµ> 5380MeV/c2 and
5220<mKµµ<5340MeV/c2, where mKµµis the Kµµ invariant mass. */
RooFormulaVar BptoKpmm("BptoKpmm", "!(@0>5380 && (5220<@1<5340))",
RooArgList(B0_MM, B0_kmuMass));
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment