#
# ccbar continuum production of the Selex state:
#
#    DsJ(2636)+ -->  D0 K+
#                    |
#                    +--> K- pi+ 
#
#
# $Id: ccbar_DsJ2636+_D0K+.dec,v 1.1.1.1 2005/03/29 17:12:43 steinke Exp $
#

#
# The following parameters attempt to harden the p* distribution
# to better match data
#
JetSetPar PARJ(41)=.10
JetSetPar PARJ(42)=.73


#
# Redefinition of particles to correctly define DsJ(2317)
# and DsJ(2458)
#
Particle D_s0*+ 2.3173   0.0001
Particle D_s0*- 2.3173   0.0001
Particle D_s1+  2.632    0.0001
Particle D_s1-  2.632    0.0001
Particle D'_s1+ 2.53535  0.001
Particle D'_s1- 2.53535  0.001

#
# Aliases we need to force desired D0 decay
#
Alias myD0 D0

#
# Charm only
#
Decay vpho 
1.00    JSCONT 4;
Enddecay

#
# Correct decay of D'_s states, now that we have "moved" them
#
Decay D'_s1+
0.5000   D*+ K0          VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
0.5000   D*0 K+          VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
0.0000   gamma D_s*+     PHSP;
0.0000   gamma D_s+      PHSP;
Enddecay

Decay D'_s1-
0.5000   D*- anti-K0     VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
0.5000   anti-D*0 K-     VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
0.0000   gamma D_s*-     PHSP;
0.0000   gamma D_s-      PHSP;
Enddecay

#
# Correct decays of DsJ(2317)
#
Decay D_s0*+
1.000   D_s+ pi0   PHSP;
Enddecay

Decay D_s0*-
1.000   D_s- pi0   PHSP;
Enddecay

#
# Forced decay of hypothetical DsJ(2636)+
#
Decay D_s1+
0.8   myD0 K+  PHSP;
Enddecay

Decay D_s1-
0.8   D0 K- PHSP;
Enddecay


#
# A D0 that decays only the way we want it to
#
Decay myD0
#
0.0385   K-  pi+                         PHSP;
Enddecay

End