# 
# Mark Williams' non-resonant B->Xs,gamma penguin decays where the
# B is a B0, anti-B0, B+ or B- and the Xs is Xsd, anti-Xsd, Xsu or
# anti-Xsu, respectively. It uses the new BTOXSGAMMA model.
# The mode for JETSET decays is 11 which gives phase-space production
# of hadrons from the available quarks.
# (Note that this could be modified for use with B_s->Xss,gamma decays
# as the BTOXSGAMMA model already allows for B_s radiative penguins.) 
#

# Reset some Jetset parameters here rather than recompiling

#   Xsd meson  (sbar-d system, introduced for b->s gamma decays)
# Set Xsd 3*Charge:
JetSetPar KCHG(455,1)=0
# Set Xsd color information:
JetSetPar KCHG(455,2)=0
# Set Xsd particle/antiparticle distinction:
JetSetPar KCHG(455,3)=1
# Set Xsd mass:
JetSetPar PMAS(455,1)=1.600
# Set Xsd width (sigma):
JetSetPar PMAS(455,2)=0.500
# Set Xsd truncation:
JetSetPar PMAS(455,3)=2.000
# Set Xsd name:
JetSetPar CHAF(455)=Xsd
# Set Xsd so it can decay:
JetSetPar MDCY(455,1)=1
# Set decay table entry pt for Xsd: 
JetSetPar MDCY(455,2)=1154
# Number of decay channels for Xsd:                
JetSetPar MDCY(455,3)=1
# Switch on Xsd decay
JetSetPar MDME(1154,1)=1
# Phase space decays into hadrons from available quarks
JetSetPar MDME(1154,2)=11
# Xsd decays into two quarks a d and an anti-s
JetSetPar KFDP(1154,1)=-3
JetSetPar KFDP(1154,2)=1

#   Xsu meson  (sbar-u system, introduced for b->s gamma decays)
# Set Xsu 3*Charge:
JetSetPar KCHG(456,1)=3
# Set Xsu color information:
JetSetPar KCHG(456,2)=0
# Set Xsu particle/antiparticle distinction:
JetSetPar KCHG(456,3)=1
# Set Xsu mass:
JetSetPar PMAS(456,1)=1.600
# Set Xsu width (sigma):
JetSetPar PMAS(456,2)=0.500
# Set Xsu truncation:
JetSetPar PMAS(456,3)=2.000
# Set Xsu name:
JetSetPar CHAF(456)=Xsu 
# Set Xsu so it can decay:
JetSetPar MDCY(456,1)=1
# Set decay table entry pt for Xsu: 
JetSetPar MDCY(456,2)=1155
# Number of decay channels for Xsu:                
JetSetPar MDCY(456,3)=1
# Switch on Xsu decay                                  
JetSetPar MDME(1155,1)=1 
# Phase space decays into hadrons from available quarks
JetSetPar MDME(1155,2)=11
# Xsu decays into two quarks a u and an anti-s
JetSetPar KFDP(1155,1)=-3
JetSetPar KFDP(1155,2)=2

Alias myB0 B0
Alias myantiB0 anti-B0
ChargeConj myB0 myantiB0

Alias myB+ B+
Alias myB- B-
ChargeConj myB+ myB-

Decay Upsilon(4S)
0.21000 myB0 anti-B0                      VSS_MIX  dm;
0.21000 B0 myantiB0                       VSS_MIX  dm;
0.04000 anti-B0 myantiB0                  VSS_MIX  dm;
0.04000 B0 myB0                           VSS_MIX dm;
#0.25 myB+ B- VSS;
#0.25 myB- B+ VSS;
Enddecay

#coeffs: 0=model, 1=fermi model, 2=mB, 3=mb, 4=mu, 5=lam1, 6=delta, 7=z, 8=nintervalS, 9=nIntervalmH  
Decay myB0
#1.00000      Xsd    gamma     BTOXSGAMMA 2 ;
#1.00000      Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.95 4.95 -0.14 0.9 0.0841 100 80;
#1.00000      Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.65 4.65 -0.52 0.9 0.0841 100 80;
1.00000      Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000      Xsd    gamma     BTOXSGAMMA 2 2 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000      Xsd    gamma     BTOXSGAMMA 2 3 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
Enddecay

Decay myantiB0
#1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 ;
#1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.95 4.95 -0.14 0.9 0.0841 100 80;
#1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.65 4.65 -0.52 0.9 0.0841 100 80;
1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 1 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 2 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     anti-Xsd    gamma     BTOXSGAMMA 2 3 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
Enddecay

#Decay myB+
#1.00000     Xsu    gamma     BTOXSGAMMA 2 1 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     Xsu    gamma     BTOXSGAMMA 2 2 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     Xsu    gamma     BTOXSGAMMA 2 3 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#Enddecay

#Decay myB-
#1.00000     anti-Xsu    gamma     BTOXSGAMMA 2 1 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     anti-Xsu    gamma     BTOXSGAMMA 2 2 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#1.00000     anti-Xsu    gamma     BTOXSGAMMA 2 3 5.27885 4.8 4.8 -0.3 0.9 0.0841 100 80;
#Enddecay

End