###########################################################################
# This is a DECAY.DEC file for B0(B0bar) -> omega KS  :                   #
#   Upsilon(4S) -> B0   B0bar                                +CC          #
#                  |    |                                                 #
#                  |     ->  omega KS                                     #
#                  |          |    |                                      #
#                  |          |     -> pi+ pi-                            #
#                   -> X      |                                           #
#                              -> pi+ pi- pi0                             #
#                                                                         #
#   Modify so that TD gives S=0.51 and C=-0.55 (values from Run1-4)        #
#                                                                         #
#   Jim Smith, Keith Ulmer, Enrico Robutti     12-FEB-2006                #
#                                                                         #
###########################################################################
#
# Create our own versions of B0, B0bar which will decay to omega KS
#
Alias myB0 B0
Alias myanti-B0 anti-B0
#
# Define CP parameters
#
#
#
Decay Upsilon(4S)
1.000 B0 anti-B0 myB0 myanti-B0 VSS_BMIX dm;
Enddecay
#
# Create our own omega to force it to decay to three pions
#
Alias myOmega  omega
#
Decay myOmega
1.000    pi-  pi+  pi0         OMEGA_DALITZ;
Enddecay
#
# Create our own Kshort to force it to decay to pi+ pi-
#
Alias myK_S0  K_S0
#
Decay myK_S0
1.0000  pi+ pi-                   PHSP;
Enddecay
#
# Force myB0 and myB0bar to always decay to myOmega myK_S0
#
# The last four args of SSD_CP are |A_f|, arg(f), |Abar_f|, arg(Abar_f) (see
# BAD 522).  |Abar_f/A_f|=|lambda|=sqrt((1-C)/(1+C)) and then normalize so
# A_f^2+Abar_f^2=1.  The sign of Abar_f is negative for omegaKs which has
# negative CP.  Then set arg(A_f)=0 and
# arg(Abar_f)=asin(S*(1+lambda^2)/(2*lambda))+2*beta # where beta=0.39.
# lambda=1.856 for C=-0.55.
Decay myB0
1.0000 myOmega myK_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
#1.0000 myOmega myK_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 0.4743 0.0000 -0.8803 -0.1231;
#1.0000 myOmega myK_S0      SSD_CP dm 0.0 1.0 -0.78 0.4743 0.0000 0.8803 -0.1231;
#1.0000 myOmega myK_S0      SSD_CP dm 0.0 1.0 0.78 0.4743 0.0000 0.8803 1.4369;
Enddecay
#
Decay myanti-B0
1.0000 myOmega myK_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
#1.0000 myOmega myK_S0 SSD_CP dm 0.0 1.0 minusTwoBeta 0.4743 0.0000 -0.8803 -0.1231;
#1.0000 myOmega myK_S0      SSD_CP dm 0.0 1.0 -0.78 0.4743 0.0000 0.8803 -0.1231;
#1.0000 myOmega myK_S0      SSD_CP dm 0.0 1.0 0.78 0.4743 0.0000 0.8803 1.4369;
Enddecay
#
End