diff --git a/DECAY_BELLE2_rel-10-00.DEC b/DECAY_BELLE2_rel-10-00.DEC new file mode 100644 index 0000000..973167c --- /dev/null +++ b/DECAY_BELLE2_rel-10-00.DEC @@ -0,0 +1,13333 @@ +# BelleII collaboration generic DECAY.DEC +# Based on EvtGen DECAY.DEC and Belle DECAY.DEC +# +# Changes +# +# October 2016, Umberto Tamponi (tamponi@to.infn.it) +# Light meson decays (rho, eta, phi, f'_0...) updated to PDG 2015 +# New Dalitz modes for eta and eta' +# +# Updates listed below are tabulated on https://xwiki.desy.de/xwiki/rest/p/4d75f +# See individual pull requests for details. +# +# December 2017, B. Oberhof (benjamin.oberhof@lnf.infn.it) +# Updated to new CPV models. See channels with #[w/CPV - BO 09/2017] +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1240/overview +# +# March 2018, A. Ermakov (ermakov@physik.uni-bonn.de) +# Updated charmless semileptonic BFs to PDG 2017 values +# https://stash.desy.de/projects/B2/repos/software/pull-requests/2087/overview +# +# March 2018, M. Merola (mario.merola@na.infn.it) +# Updated B -> tau nu BR to PDG2018 values +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1711/overview +# +# March 2018, R. Godang (godang@southalabama.edu) +# Updated the BRs of the semitauonic decays: B0, anti-B0, B-, B+ decays including pi+ tau- anti-nu_tau +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1769/overview +# https://github.com/belle2/basf2/blob/629dc58140ef04e18106dc7a7df9fff27974c98c/generators/evtgen/decayfiles/DECAY_BELLE2.DEC +# +# March 2018, T. Lueck (lueck@pi.infn.it) +# Updated B -> X_c l nu BR to PDG2018 values +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1739/overview +# +# March 2018, W. Sutcliffe (william.sutcliffe@desy.de) +# Updated B -> s gamma and B -> s l l +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1712/overview# +# +# March 2018, S. Jia & Y. Li (jiasen@buaa.edu.cn, liyubo@buaa.edu.cn) +# Updated BR of charmonium modes in B_s0/anti-B_s0 decays and chi_cJ decays to PDG2017 +# Added new charmonium modes in B_s0/anti-B_s0 decays and chi_cJ decays +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1816/overview +# +# April 2018, B. Oberhof (benjamin.oberhof@lnf.infn.it) +# Updated CP-violating decays to charmonium + X and CP-violating decays to D mesons. See channels with #[w/CPV - BO 09/2017] +# https://stash.desy.de/projects/B2/repos/software/pull-requests/1932/overview +# +# April 2018, P. Goldenzweig (pablo.goldenzweig@kit.edu) +# Corrected unwanted asymmetries and ordered decays to be symmetric between Decay X and Decay anti-X +# https://stash.desy.de/projects/B2/repos/software/pull-requests/2024/overview +# +# May 2018, S. Duell (s.duell@physik.uni-bonn.de) +# Updated B->D(*) Form Factor parameters to HFLAV 2016 values (1612.07233) +# https://stash.desy.de/projects/B2/repos/software/pull-requests/2147/overview +# +# March 2019, L. Cao (lu.cao@desy.de) +# Updated B -> D l nu form factor and parameters to the latest Belle measurement (1510.03657) +# Updated B -> pi l nu form factor to BTOXELNU BCL with parameters based on HFLAV 2016 (1612.07233) from M. Prim (markus.prim@kit.edu) +# Updated B - > rho/omega l nu from factor to BTOXELNU BCL with parameters from M. Prim +# https://stash.desy.de/projects/B2/repos/software/pull-requests/3721/overview +# +# March 2020, E. R. Oxford (emma.oxford@desy.de) +# **All changes are documented in BELLE2-NOTE-PH-2020-008 on the Belle II Document Server** +# Substantial updates to charm decays, such as: +# Updated D0/anti-D0 decay tables +# Updated D+/D- decay tables +# Updated D_s+/D_s- decay tables +# Updated charmed baryon decay tables +# Added ModelAliases for SLPOLE and SLBKPOLE models now used in D semi-leptonic decays +# Added two new "fake" resonances, phiR and f_2R, and improved documentation for K*0R and K*BR "fake" resonances +# Added decay table for sigma_0 particle (f_0(500) in the PDG) +# +# Sept.-Dec. 2020, P. Goldenzweig (pablo.goldenzweig@kit.edu) +# Updated hundreds of B0 and B+ BFs and some FFs as described in https://gitlab.desy.de/belle2/software/basf2/-/issues/7124. +# Further updates and corrections applied in https://stash.desy.de/projects/B2/repos/software/pull-requests/7501/overview. +# Added the remaining rare decays from the Belle mixedrare.dec and chargedrare.dec files. +# Pythia hadron contributions rescaled by S. Duell (https://stash.desy.de/projects/B2/repos/software/pull-requests/6690/overview) +# All Decay B0 & Decay anti-B0 channels have symmetric values. Same for B+ / B-. I.e., no Acp included. +# Updates flagged with #[BF update - PG 09/2020]. +# +# June 2021, U. Tamponi (tamponi@to.infn.it) +# Updated all the narrow bottomonia (M < Y(4S)) +# +# June 2021 Automated: Converted PYTHIA ids +# Automated conversion of all PYTHIA decay ids to Pythia 8 using conversion defined in evtgen/src/EvtGenExternal/EvtPythiaEngine.cpp +# +# December 2021, V. Bertacchi (bertacchi@cppm.in2p3.fr) +# fixed a_1+ and a_10 BR decay table to restore isospin symmetry and reimplement the resonant contribution (KK*) +# rescaled B->DKK* BR to take into account of resonant contribution B->D(a1->KK*), preserving the BR(NR+Res)= Belle measurement (https://doi.org/10.1016/S0370-2693(02)02373-0) +# https://gitlab.desy.de/belle2/software/basf2/-/merge_requests/581 +# +# March 2022, H. Svidras (henrikas.svidras@desy.de) +# Corrected the BTOXSGAMMA model branching fraction to the correct inclusive B->Xs gamma decay branching fraction in the PDG + +# June 2023, Niharika Rout (niharika.rout@ts.infn.it) +# Updated the B -> D(**)npi BFs according to latest PDG/paper measurements. + +# June 2023, Meihong Liu (liumh19@fudan.edu.cn) +# Updated the B0 -> D+(**)npi BFs according to latest PDG/paper measurements. + +Define qoverp_incohMix_B_s0 1.0 +Define dm_incohMix_B_s0 17.8e12 +Define qoverp_incohMix_B0 1.0 +Define dm_incohMix_B0 0.507e12 +# Old definition of dm still in some decay models +Define dm 0.507e12 +#Define dgamma 0 +#Define qoverp 1 +#Define phaseqoverp 0 +#Define values for B0s mixing +#Define dms 20.e12 +# DeltaGammas corresponds to DG/G = 10% +#Define dgammas 6.852e10 +# Activate incoherent Mixing +### TODO: find a way to give mixing parameters through decay file to EvtGen +#####yesIncoherentB0Mixing dm dgamma +#####yesIncoherentBsMixing dms dgammas +# define the values of the CKM angles (alpha=70, beta=40) +Define alpha 1.365 +Define beta 0.39 +Define gamma 1.387 +Define twoBetaPlusGamma 2.167 +Define betaPlusHalfGamma 1.0835 +Define minusGamma -1.387 +Define minusTwoBeta -0.78 + +# New definitions for psiKstar modes (Lange, July 26, 2000) +Define PKHplus 0.159 +Define PKHzero 0.775 +Define PKHminus 0.612 +Define PKphHplus 1.563 +Define PKphHzero 0.0 +Define PKphHminus 2.712 + +Define Aplus 0.490 +Define Azero 1.10 +Define Aminus 0.4 + +Define phAplus 2.5 +Define phAzero 0.0 +Define phAminus -0.17 + + +# These model aliases are used to specify parameters +# for SLPOLE and SLBKPOLE models used in D semi-leptonic +# decays (see EvtGen documentation for details) +ModelAlias SLBKPOLE_DtoKlnu SLBKPOLE 1.0 0.303 1.0 2.112; +ModelAlias SLBKPOLE_Dtopilnu SLBKPOLE 1.0 0.281 1.0 2.010; +ModelAlias SLBKPOLE_Dtoetalnu SLBKPOLE 1.0 0.281 1.0 2.010; +ModelAlias SLBKPOLE_Dtoetaplnu SLBKPOLE 1.0 0.281 1.0 2.010; +ModelAlias SLPOLE_DtoKstarlnu SLPOLE 1.00 -0.558 0.0 1.0 0.85 -0.558 0.0 1.0 1.50 -0.790 0.0 1.0 0.00 -0.558 0.0 1.0; +ModelAlias SLPOLE_Dtorholnu SLPOLE 1.00 -0.605 0.0 1.0 0.60 -0.605 0.0 1.0 1.50 -0.871 0.0 1.0 0.00 -0.605 0.0 1.0; +ModelAlias SLPOLE_Dtoomegalnu SLPOLE 1.00 -0.605 0.0 1.0 0.60 -0.605 0.0 1.0 1.50 -0.871 0.0 1.0 0.00 -0.605 0.0 1.0; + + +# These particle aliases are used in for CP violating decays +# in which the decay distributions depends on how the K* decayed. +Alias K*L K*0 +Alias K*S K*0 +Alias K*BL anti-K*0 +Alias K*BS anti-K*0 +Alias K*0T K*0 +Alias anti-K*0T anti-K*0 +Alias anti-K_0*0N anti-K_0*0 +Alias K_0*0N K_0*0 +ChargeConj K*L K*BL +ChargeConj K*S K*BS +ChargeConj K*0T anti-K*0T +ChargeConj K_0*0N anti-K_0*0N + + +# These particle aliases are used to avoid double-counting +# resonances included in D_DALITZ models +Alias phiR phi +Alias f_2R f_2 +Alias K*BR anti-K*0 +Alias K*0R K*0 +ChargeConj K*0R K*BR +Alias K*-R K*- +Alias K*+R K*+ +ChargeConj K*+R K*-R + +# PR LHCb: Alias for signal productions in LHCb +Alias B0sig B0 +Alias anti-B0sig anti-B0 +ChargeConj B0sig anti-B0sig +Alias B+sig B+ +Alias B-sig B- +ChargeConj B+sig B-sig +Alias B_s0sig B_s0 +Alias anti-B_s0sig anti-B_s0 +ChargeConj B_s0sig anti-B_s0sig +Alias B_c+sig B_c+ +Alias B_c-sig B_c- +ChargeConj B_c+sig B_c-sig +Alias eta_bsig eta_b +ChargeConj eta_bsig eta_bsig +Alias h_bsig h_b +ChargeConj h_bsig h_bsig +Alias Sigma_b-sig Sigma_b- +Alias anti-Sigma_b+sig anti-Sigma_b+ +ChargeConj Sigma_b-sig anti-Sigma_b+sig +Alias Lambda_b0sig Lambda_b0 +Alias anti-Lambda_b0sig anti-Lambda_b0 +ChargeConj Lambda_b0sig anti-Lambda_b0sig +Alias Omega_b-sig Omega_b- +Alias anti-Omega_b+sig anti-Omega_b+ +ChargeConj Omega_b-sig anti-Omega_b+sig +Alias Xi_b-sig Xi_b- +Alias anti-Xi_b+sig anti-Xi_b+ +ChargeConj Xi_b-sig anti-Xi_b+sig +Alias Xi_b0sig Xi_b0 +Alias anti-Xi_b0sig anti-Xi_b0 +ChargeConj Xi_b0sig anti-Xi_b0sig +Alias J/psisig J/psi +ChargeConj J/psisig J/psisig +Alias chi_c0sig chi_c0 +ChargeConj chi_c0sig chi_c0sig +Alias chi_c1sig chi_c1 +ChargeConj chi_c1sig chi_c1sig +Alias chi_c2sig chi_c2 +ChargeConj chi_c2sig chi_c2sig +Alias psi(2S)sig psi(2S) +ChargeConj psi(2S)sig psi(2S)sig +Alias psi(3770)sig psi(3770) +ChargeConj psi(3770)sig psi(3770)sig +Alias D*+sig D*+ +Alias D*-sig D*- +ChargeConj D*+sig D*-sig +Alias D*0sig D*0 +Alias anti-D*0sig anti-D*0 +ChargeConj D*0sig anti-D*0sig +Alias D0sig D0 +Alias anti-D0sig anti-D0 +ChargeConj D0sig anti-D0sig +Alias D+sig D+ +Alias D-sig D- +ChargeConj D+sig D-sig +Alias D_s+sig D_s+ +Alias D_s-sig D_s- +ChargeConj D_s+sig D_s-sig +Alias Lambda_c+sig Lambda_c+ +Alias anti-Lambda_c-sig anti-Lambda_c- +ChargeConj Lambda_c+sig anti-Lambda_c-sig +Alias tau+sig tau+ +Alias tau-sig tau- +ChargeConj tau+sig tau-sig +Alias Upsilonsig Upsilon +ChargeConj Upsilonsig Upsilonsig +Alias Upsilon(2S)sig Upsilon(2S) +ChargeConj Upsilon(2S)sig Upsilon(2S)sig +Alias Upsilon(3S)sig Upsilon(3S) +ChargeConj Upsilon(3S)sig Upsilon(3S)sig +Alias Upsilon(4S)sig Upsilon(4S) +ChargeConj Upsilon(4S)sig Upsilon(4S)sig +Alias Upsilon(5S)sig Upsilon(5S) +ChargeConj Upsilon(5S)sig Upsilon(5S)sig +Alias X_1(3872)sig X_1(3872) +ChargeConj X_1(3872)sig X_1(3872)sig +Alias h_csig h_c +ChargeConj h_csig h_csig +Alias Sigma+sig Sigma+ +Alias anti-Sigma-sig anti-Sigma- +ChargeConj Sigma+sig anti-Sigma-sig +Alias Lambda0sig Lambda0 +Alias anti-Lambda0sig anti-Lambda0 +ChargeConj Lambda0sig anti-Lambda0sig +Alias B_10sig B_10 +Alias anti-B_10sig anti-B_10 +ChargeConj B_10sig anti-B_10sig +Alias B_2*0sig B_2*0 +Alias anti-B_2*0sig anti-B_2*0 +ChargeConj B_2*0sig anti-B_2*0sig +Alias B_s10sig B_s10 +Alias anti-B_s10sig anti-B_s10 +ChargeConj B_s10sig anti-B_s10sig +Alias B_s2*0sig B_s2*0 +Alias anti-B_s2*0sig anti-B_s2*0 +ChargeConj B_s2*0sig anti-B_s2*0sig +Alias B_1+sig B_1+ +Alias B_1-sig B_1- +ChargeConj B_1+sig B_1-sig +Alias B_2*+sig B_2*+ +Alias B_2*-sig B_2*- +ChargeConj B_2*+sig B_2*-sig + + +#JetSet parameter modifications +#(Very important that there are no blank spaces in the parameter string!) +#Turn of B0-B0B mixing in JetSet: +#JetSetPar MSTJ(26)=0 +#Control of L=1 mesons (in order: J1S0 J0S1 J1S1 J2S1) - commented out by NB/WP +#JetSetPar PARJ(14)=0.05 +#JetSetPar PARJ(15)=0.05 +#JetSetPar PARJ(16)=0.05 +#JetSetPar PARJ(17)=0.05 +#cut-off parameter used to stop fragmentation process (should not be changed) +#JetSetPar PARJ(33)=0.3 SET NOW IN GAUSS + + + +# Turn on FSR for all decays +yesFSR + + + +################################################################################## +########## ########## +########## Upsilon(4S) ########## +########## ########## +################################################################################## +# +# Updated to PDG 2011 +# +# Use VSS_BMIX mixing decay model (DK,28-Oct-1999) +Decay Upsilon(4S) +0.515122645 B+ B- VSS; +0.483122645 B0 anti-B0 VSS_BMIX dm; +0.000015583 e+ e- VLL; +0.000015766 mu+ mu- VLL; +0.000015766 tau+ tau- VLL; +0.000084099 Upsilon(2S) pi+ pi- VVPIPI; +0.000044342 Upsilon(2S) pi0 pi0 VVPIPI; +0.000080123 Upsilon pi+ pi- VVPIPI; +0.000044342 Upsilon pi0 pi0 VVPIPI; +0.000194392 Upsilon eta PARTWAVE 0.0 0.0 1.0 0.0 0.0 0.0; +# BF ~ (2J+1)E^3_gamma; see PRL 94, 032001 +# V-> gamma S Partial wave (L,S)=(0,0) +0.000092625 gamma chi_b0(3P) HELAMP 1. 0. 1. 0.; +# V-> gamma V Partial wave (L,S)=(0,1) +0.000138938 gamma chi_b1(3P) HELAMP 1. 0. 1. 0. -1. 0. -1. 0.; +# V-> gamma T Partial wave (L,S)=(0,1) +0.000129084 gamma chi_b2(3P) HELAMP 2.4494897 0. 1.7320508 0. 1. 0. 1. 0. 1.7320508 0. 2.4494897 0.; +# V-> gamma S Partial wave (L,S)=(0,0) +0.000002956 gamma chi_b0(2P) HELAMP 1. 0. 1. 0.; +# V-> gamma V Partial wave (L,S)=(0,1) +0.000007883 gamma chi_b1(2P) HELAMP 1. 0. 1. 0. -1. 0. -1. 0.; +# V-> gamma T Partial wave (L,S)=(0,1) +0.000011825 gamma chi_b2(2P) HELAMP 2.4494897 0. 1.7320508 0. 1. 0. 1. 0. 1.7320508 0. 2.4494897 0.; +0.000837571 g g g PYTHIA 92; +0.000039415 gamma g g PYTHIA 92; +Enddecay + + + + + +################################################################################## +########## ########## +########## B0 anti-B0 ########## +########## ########## +################################################################################## + + +############################################# +# # +# B0 # +# # +############################################# + +Decay B0 + +# b -> c semileptonic +# +# Updated to PDG 2017: using the Isospin average between B+ and B0. +# +0.051100000 D*- e+ nu_e BTOXELNU BGL 0.02596 -0.06049 0.01311 0.01713 0.00753 -0.09346; #FF values taken from [arxiv:2008.09341] +0.021400000 D- e+ nu_e BTOXELNU BGL 0.0126 -0.094 0.34 -0.1 0.0115 -0.057 0.12 0.4; #[Belle measurement: 1510.03657] +0.007040000 D_1- e+ nu_e BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; +0.003620000 D_0*- e+ nu_e BTOXELNU LLSW 0.68 -0.2 0.3; +0.004010000 D'_1- e+ nu_e BTOXELNU LLSW 0.68 -0.2 0.3; +0.003470000 D_2*- e+ nu_e BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; +# The measured modes are already saturated by B->X_clnu; X_c->D(*)pi. +# The values below give a rough estimate of max. allowed BF taking into uncertainties. +0.000920000 anti-D*0 pi- e+ nu_e GOITY_ROBERTS; +0.000460000 D*- pi0 e+ nu_e GOITY_ROBERTS; +0.000920000 anti-D0 pi- e+ nu_e GOITY_ROBERTS; +0.000460000 D- pi0 e+ nu_e GOITY_ROBERTS; +# The B->Dpipilnu is the measured rate not covered by B->D_1lnu; D_1->Dpipi; modes with pi0 are inferred assuming Isospin symmetry. +0.000210000 D- pi+ pi- e+ nu_e PHSP; +0.000140000 anti-D0 pi- pi0 e+ nu_e PHSP; +0.000140000 D- pi0 pi0 e+ nu_e PHSP; +# B->D*pipilnu modes with pi0 are inferred assuming Isospin symmetry. +0.001050000 D*- pi+ pi- e+ nu_e PHSP; +0.000700000 anti-D*0 pi- pi0 e+ nu_e PHSP; +0.000700000 D*- pi0 pi0 e+ nu_e PHSP; +# Fill the gap between inclusive B->X_clnu and sum of exclusive with hypothetical decay (unmeasured). +0.002170000 D- eta e+ nu_e PHSP; +0.002170000 D*- eta e+ nu_e PHSP; + +0.051100000 D*- mu+ nu_mu BTOXELNU BGL 0.02596 -0.06049 0.01311 0.01713 0.00753 -0.09346; #FF values taken from [arxiv:2008.09341] +0.021400000 D- mu+ nu_mu BTOXELNU BGL 0.0126 -0.094 0.34 -0.1 0.0115 -0.057 0.12 0.4; #[Belle measurement: 1510.03657] +0.007040000 D_1- mu+ nu_mu BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; +0.003620000 D_0*- mu+ nu_mu BTOXELNU LLSW 0.68 -0.2 0.3; +0.004010000 D'_1- mu+ nu_mu BTOXELNU LLSW 0.68 -0.2 0.3; +0.003470000 D_2*- mu+ nu_mu BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; +# The measured modes are already saturated by B->X_clnu; X_c->D(*)pi. +# The values below give a rough estimate of max. allowed BF taking into uncertainties. +0.000920000 anti-D*0 pi- mu+ nu_mu GOITY_ROBERTS; +0.000460000 D*- pi0 mu+ nu_mu GOITY_ROBERTS; +0.000920000 anti-D0 pi- mu+ nu_mu GOITY_ROBERTS; +0.000460000 D- pi0 mu+ nu_mu GOITY_ROBERTS; +# The B->Dpipilnu is the measured rate not covered by B->D_1lnu; D_1->Dpipi; modes with pi0 are inferred assuming Isospin symmetry. +0.000210000 D- pi+ pi- mu+ nu_mu PHSP; +0.000140000 anti-D0 pi- pi0 mu+ nu_mu PHSP; +0.000140000 D- pi0 pi0 mu+ nu_mu PHSP; +# B->D*pipilnu ; modes with pi0 are inferred assuming isospin symmetry +0.001050000 D*- pi+ pi- mu+ nu_mu PHSP; +0.000700000 anti-D*0 pi- pi0 mu+ nu_mu PHSP; +0.000700000 D*- pi0 pi0 mu+ nu_mu PHSP; +# Fill the gap between inclusive B->X_clnu and sum of exclusive with hypothetical decay (unmeasured). +0.002170000 D- eta mu+ nu_mu PHSP; +0.002170000 D*- eta mu+ nu_mu PHSP; + +# b -> c tau nu +# +# D(*)0 tau nu: FF values taken from HFLAV 2016 (1612.07233); BF based on the HFLAV theory average of R(D/D*) of 0.299 and 0.258 and our D/D* l nu BFs in the MC. +0.013183800 D*- tau+ nu_tau BTOXELNU HQET3 0.912 1.205 1.15 1.404 0.854; #[BF update - PG 09/2020] +0.006398600 D- tau+ nu_tau BTOXELNU HQET3 1.0541 1.128 1.0; #[BF update - PG 09/2020] + +0.001300000 D_1- tau+ nu_tau BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; +0.001300000 D_0*- tau+ nu_tau BTOXELNU LLSW 0.68 -0.2 0.3; +0.002000000 D'_1- tau+ nu_tau BTOXELNU LLSW 0.68 -0.2 0.3; +0.002000000 D_2*- tau+ nu_tau BTOXELNU LLSW 0.71 -1.6 -0.5 2.9; + +# B -> cc= s +# +#$ forbid_pythia() +0.000000000 J/psi K0 PHSP; +0.000436500 J/psi K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[w/CPV - BO 09/2017] +0.000436500 J/psi K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; + +#$ forbid_pythia() +0.000000000 J/psi K*0 PHSP; +0.000211667 J/psi K*S SVV_CP beta dm 1 Aplus phAplus Azero phAzero Aminus phAminus; #J/psi K*0 is Gamma_200 PDGLive = 1.27e-3. #[Comment - PG 09/2020] #[w/CPV - BO 09/2017] +0.000211667 J/psi K*L SVV_CP beta dm -1 Aplus phAplus Azero phAzero Aminus phAminus; +0.000846667 J/psi K*0T SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus; + +# Requires subtraction of K* component. +#$ BF('S042:Desig=4') - BF('S042:Desig=22') * BF('M018:Desig=12') * 2 / 3 +0.000303333 J/psi K+ pi- PHSP; #[BF update - PG 09/2020] +# This branching fraction is not measured. The value is assigned +# from isospin relation. +#$ BF('B0 -> J/psi K+ pi-') / 2 +0.000151665 J/psi K0 pi0 PHSP; #[BF update - PG 09/2020] + +0.000016570 J/psi pi0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #[BF update - PG 09/2020] +0.000017700 J/psi omega SVV_CP beta dm 1 Aplus phAplus Azero phAzero Aminus phAminus; #[BF update - PG 09/2020] + +0.000025450 J/psi rho0 SVV_CP beta dm 1 Aplus phAplus Azero phAzero Aminus phAminus; #[BF update - PG 09/2020] +0.000003260 J/psi f_2 PHSP; #[BF update - PG 09/2020] +0.000020300 J/psi pi+ pi- PHSP; #Gamma_210 PDGLive => 0.000039440 includes rho0, f0(500), f0(980), f2, rho(1450). Only rho0 and f_2 included here. Add others and subtract from inclusive BF. #[BF update - PG 09/2020] + +0.000434530 J/psi K0 pi- pi+ PHSP; #[BF update - PG 09/2020] +#0.000350000 J/psi K0 pi0 pi0 PHSP; #No PDGLive value. Keep commented out. #[Comment - PG 09/2020] +#0.000350000 J/psi K+ pi- pi0 PHSP; #No PDGLive value. Keep commented out. #[Comment - PG 09/2020] + +0.001300000 J/psi K_10 SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0; +0.000100000 J/psi K'_10 SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0; +0.000500000 J/psi K_2*0 PHSP; +0.000048680 J/psi K0 phi PHSP; #[BF update - PG 09/2020] + +0.000002500 J/psi K+ K- PHSP; #Gamma_219 PDGLive includes a0, but we do not include J/psi a0, so no need to subtract here. #[BF update - PG 09/2020] +0.000108000 J/psi eta K0 PHSP; +0.000230000 J/psi omega K0 PHSP; +0.000010780 J/psi eta PHSP; #[BF update - PG 09/2020] +0.000007610 J/psi eta' PHSP; #[BF update - PG 09/2020] +0.000024890 J/psi K0 K+ K- PHSP; #[BF update - PG 09/2020] +0.000540000 J/psi K0 rho0 PHSP; +0.000770000 J/psi K*+ pi- PHSP; #[BF update - PG 09/2020] +0.000014240 J/psi pi+ pi- pi+ pi- PHSP; #Gamma_229 PDGLive. #[BF update - PG 09/2020] +0.000008440 J/psi f_1 PHSP; #Gamma_230 PDGLive. #[BF update - PG 09/2020] +0.000660000 J/psi K*0 pi+ pi- PHSP; + +#$ forbid_pythia() +0.000000000 psi(2S) K0 PHSP; +0.000290000 psi(2S) K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; +0.000290000 psi(2S) K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; + +#$ forbid_pythia() +0.000000000 psi(2S) K*0 PHSP; +0.000098333 psi(2S) K*S SVV_CP beta dm 1 Aplus phAplus Azero phAzero Aminus phAminus; #psi(2S) K*0 is Gamma_257 PDGLive = 5.9e-4. #[Comment - PG 09/2020] #[w/CPV - BO 09/2017] +0.000098333 psi(2S) K*L SVV_CP beta dm -1 Aplus phAplus Azero phAzero Aminus phAminus; +0.000393333 psi(2S) K*0T SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus; + +# Requires subtraction of K* component. +#$ BF('S042:Desig=70') - BF('S042:Desig=69') * BF('M018:Desig=12') * 2 / 3 +0.000186667 psi(2S) K+ pi- PHSP; #[BF update - PG 09/2020] +# This branching fraction is not measured. The value is assigned +# from isospin relation. +#$ BF('B0 -> psi(2S) K+ pi-') / 2 +0.000093333 psi(2S) K0 pi0 PHSP; #[BF update - PG 09/2020] +0.000200000 psi(2S) K0 pi- pi+ PHSP; +0.000100000 psi(2S) K0 pi0 pi0 PHSP; +0.000100000 psi(2S) K+ pi- pi0 PHSP; +0.000400000 psi(2S) K_10 PHSP; +0.000022090 psi(2S) pi+ pi- PHSP; #[BF update - PG 09/2020] +0.000011700 psi(2S) pi0 PHSP; + +#$ forbid_pythia() +0.000000000 eta_c K0 PHSP; +0.000400000 eta_c K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[w/CPV - BO 09/2017] +0.000400000 eta_c K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; + +#$ forbid_pythia() +0.000000000 eta_c K*0 PHSP; +0.000105000 K*S eta_c SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[w/CPV - BO 09/2017] +0.000105000 K*L eta_c SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; +0.000420000 K*0T eta_c SVS; + +0.000200000 eta_c K+ pi- PHSP; +0.000100000 eta_c K0 pi0 PHSP; +0.000200000 eta_c K0 pi- pi+ PHSP; +0.000100000 eta_c K0 pi0 pi0 PHSP; +0.000100000 eta_c K+ pi- pi0 PHSP; + +#$ forbid_pythia() +0.000000000 eta_c(2S) K0 PHSP; +# Br(B0 -> eta_c K_S0) = Br(B+ -> eta_c K+) * Br(B0 -> eta_c (-> p anti-p pi+ pi-) K_S0) / Br (B+ -> eta_c (-> p anti-p pi+ pi-) K+) = 4.4e-4 * (4.2e-7 / 1.12e-6) = 0.00016500. #[Comment - KL 09/2020] +0.000165000 eta_c(2S) K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[BF update - PG 09/2020] #[w/CPV - BO 09/2017] +0.000165000 eta_c(2S) K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #[BF update - PG 09/2020] + +#$ forbid_pythia() +0.000000000 eta_c(2S) K*0 PHSP; +0.000110000 K*S eta_c(2S) SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[w/CPV - BO 09/2017] +0.000110000 K*L eta_c(2S) SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; +0.000440000 K*0T eta_c(2S) SVS; + +0.000080000 eta_c(2S) K+ pi- PHSP; +0.000050000 eta_c(2S) K0 pi0 PHSP; +0.000080000 eta_c(2S) K0 pi- pi+ PHSP; +0.000050000 eta_c(2S) K0 pi0 pi0 PHSP; +0.000050000 eta_c(2S) K+ pi- pi0 PHSP; + +#$ forbid_pythia() +0.000000000 chi_c0 K0 PHSP; +0.000073500 chi_c0 K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #[w/CPV - BO 09/2017] +0.000073500 chi_c0 K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; + +#$ forbid_pythia() +0.000000000 chi_c0 K*0 PHSP; +0.000028000 K*S chi_c0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #[w/CPV - BO 09/2017] +0.000028000 K*L chi_c0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; +0.000113000 K*0T chi_c0 SVS; + +0.000200000 chi_c0 K+ pi- PHSP; +0.000100000 chi_c0 K0 pi0 PHSP; +0.000200000 chi_c0 K0 pi- pi+ PHSP; +0.000100000 chi_c0 K0 pi0 pi0 PHSP; +0.000100000 chi_c0 K+ pi- pi0 PHSP; + +# B0 -> chi_c1 X + +# B0 -> chi_c1 K0 (PDG 2019 data) +#$ forbid_pythia() +0.000000000 chi_c1 K0 PHSP; +0.000196500 chi_c1 K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #[w/CPV - BO 09/2017] +0.000196500 chi_c1 K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; + +#$ forbid_pythia() +0.000000000 chi_c1 K*0 PHSP; +0.000039570 chi_c1 K*S SVV_CP beta dm 1 Aplus phAplus Azero phAzero Aminus phAminus; #[w/CPV - BO 09/2017] +0.000039570 chi_c1 K*L SVV_CP beta dm -1 Aplus phAplus Azero phAzero Aminus phAminus; +0.000158280 chi_c1 K*0T SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus; #[w/CPV - BO 09/2017] + +# Requires subtraction of K* component. +#$ BF('S042:Desig=441') - BF('S042:Desig=127') * BF('M018:Desig=12') * 2 / 3 +0.000338720 chi_c1 K+ pi- PHSP; #[Comment - PG 09/2020] +# This branching fraction is not measured. The value is assigned +# from isospin relation. +#$ BF('B0 -> chi_c1 K+ pi-') / 2 +0.000169360 chi_c1 K0 pi0 PHSP; + +0.000316000 chi_c1 K0 pi- pi+ PHSP; #[BF update - PG 09/2020] +0.000200000 chi_c1 K0 pi0 pi0 PHSP; +0.000352000 chi_c1 K+ pi- pi0 PHSP; #[BF update - PG 09/2020] +0.000011200 chi_c1 pi0 PHSP; + +#$ forbid_pythia() +0.000000000 chi_c2 K0 PHSP; +0.000007500 chi_c2 K_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0; #Taken from UL of <1.5e-5. #[w/CPV - BO 09/2017] +0.000007500 chi_c2 K_L0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; + +0.000048820 chi_c2 K*0 PHSP; #[BF update - PG 09/2020] +0.000072000 chi_c2 K+ pi- PHSP; #No PDGLive value. Keep original value in .DEC. #[Comment - PG 09/2020] +0.000036000 chi_c2 K0 pi0 PHSP; #Assume BF(chi_c2 K0 pi0) = BF(chi_c2 K+ pi-) /2. #[BF update - PG 09/2020] + +0.000010000 chi_c2 K0 pi- pi+ PHSP; #No PDGLive value. Set BF