D**

5/25/2004 AN D** mass resolution in MC

MC may be biased to lower pt value becasue the trigger efficiencies are not taken into account. To correct for this we can apply
scaling (weights while filling histograms). The plot below shows the ratio of Pt (lepton) distributions for tight mu D0 candidates in MC and data
(i.e. MC distribution is divided by data distribution). The dependence can be described by exp+const.

D**(2420) mass before and after scaling

D**(2420) fit with relativistic Breit-Wigner from PDG. BW is convoluted with a gaussian. Amplitude, mass, BW width
and gaussian sigma were free parameters of the fit. MC uses PDG values for mass and width (2.422 GeV and 19 MeV).
It looks the width of BW correlates with sigma of gaussian.

D**(2460) fit with relativistic Breit-Wigner from PDG. BW is convoluted with a gaussian. Amplitude, mass, BW width
and gaussian sigma were free parameters of the fit.

11/17/2003 AN D** study in MC

used MC B-> muon + D0 (->K pion) + X (requests 8681 and 8949) + standard D0 particles with no additional cuts. total # of entries = 27460. each entry has a reconstructed D0 candidate. Total # of reconstructed B->l D0 matched to MC particles is 25633. in the below table no oscillation was required for B0.

channel	D meson ID	total candidates for mu + D	reco'd D0->K pi	reco'd D*+->(Kpi)pi+	comment	k-factor
B- (ID=-521) => mu D^0	421	3180	2890	0	main channel	D0 k-factor
B- (ID=-521) => mu D^*0	423	9172	8402	0	main channel	D*0 k-factor
B- (ID=-521) => mu D^*0_2	425	202	180	34	D** channel	D*0_2 k-factor
B- (ID=-521) => mu D^*-	413	143	115	84	D** channel
B- (ID=-521) => mu D^*0_0	10421	69	62	0	D** channel
B- (ID=-521) => mu D^0_1	10423	504	435	161	D** channel
B- (ID=-521) => mu D^*0_1	20423	305	274	120	D** channel

B0 (ID=-511) => mu D^*-	413	6082	5504	3994	main channel
B0 (ID=-511) => mu D^*-_2	415	328	285	26	D** channel
B0 (ID=-511) => mu D^0	421	730	624	0	D** channel
B0 (ID=-511) => mu D^*0	423	185	153	0	D** channel
B0 (ID=-511) => mu D^*-_0	10411	136	120	0	D** channel
B0 (ID=-511) => mu D^-_1	10413	575	515	98	D** channel
B0 (ID=-511) => mu D^*-_1	20413	370	319	55	D** channel

observations :
1) xxx% of B0 decay to D0 channels without forming D*-
2) xxx% of B- decay to D*- channel

k-factor is defined as Pt (reco'ed D(*)+muon) / Pt (MC B). Distributions of the k-factors for 4 most important cases is shown below

Here are the text files with the above four histograms. the format is the following
(bin #, left edge of the bin, right edge of the bin, content of the bin) :
B+ -> D0 or D*0
B0 -> D*-
B -> D** -> D0
B -> D** -> D*-

below are Branching Fractions for B0 and B- used in MC

Decay anti-B0
#                   b -> c semileptonic
#
0.0560   D*+    mu- anti-nu_mu         PHOTOS HQET 0.92 1.18 0.72;
0.0056   D_1+   mu- anti-nu_mu         PHOTOS ISGW2;
0.3333    D*+ pi0
0.6667    D*0 pi+
0.0020   D_0*+   mu- anti-nu_mu         PHOTOS ISGW2;
0.6667   D0 pi+
0.0037   D'_1+   mu- anti-nu_mu         PHOTOS ISGW2;
0.3333    D*+ pi0
0.6667    D*0 pi+
0.0037   D_2*+   mu- anti-nu_mu         PHOTOS ISGW2;
0.1030    D*+ pi0
0.2090    D*0 pi+
0.4590    D0 pi+
0.0010   D*+ pi0   mu- anti-nu_mu    PHOTOS   GOITY_ROBERTS;
0.0020   D*0 pi+   mu- anti-nu_mu    PHOTOS   GOITY_ROBERTS;
0.0060   D0   pi+   mu- anti-nu_mu    PHOTOS   GOITY_ROBERTS;

i.e. Br(B0 => D* mu X ) = 5.6+0.56*0.33+0.37*0.33+0.37*0.10+0.1 = 5.6 (direct D*) + 0.44 (D* from D**) = 6.04%
   Br(B0 => D0 pi+ mu X ) = 0.56*0.67+0.20*0.67+0.37*0.67+0.37*0.67+0.2+0.6 = 1.81% (D0 from D**)

Decay B-
#
0.0560   D*0   mu- anti-nu_mu         PHOTOS   HQET 0.92 1.18 0.72;
0.0210   D0    mu- anti-nu_mu         PHOTOS   ISGW2;
0.0056   D_10 mu- anti-nu_mu         PHOTOS   ISGW2;
0.3333    D*0 pi0
0.6667    D*+ pi-
0.0020   D_0*0 mu- anti-nu_mu        PHOTOS    ISGW2;
0.3333   D0 pi0
0.0037   D'_10 mu- anti-nu_mu        PHOTOS    ISGW2;
0.6667    D*+ pi-
0.3333    D*0 pi0
0.0037   D_2*0 mu- anti-nu_mu        PHOTOS    ISGW2;
0.2090    D*+ pi-
0.1030    D*0 pi0
0.2290    D0 pi0
0.0020   D*+   pi-   mu-   anti-nu_mu PHOTOS   GOITY_ROBERTS;
0.0010   D*0   pi0   mu-   anti-nu_mu PHOTOS   GOITY_ROBERTS;
0.0030   D0    pi0   mu-   anti-nu_mu PHOTOS   GOITY_ROBERTS;

i.e Br(B- => D0 mu X ) = 5.6+2.1+0.56*0.33+0.2*0.33+0.37*0.33+0.37*0.33+0.1+0.3 = 7.7 (direct D0 and from D*0) + 0.90 (D0 from D**) = 8.6%
Br(B- => D*+ pi- mu X ) = 0.56*0.67+0.37*0.67+0.37*0.21+0.2 = 0.90% (D* from D**)

thus at the generator level the composition of the samples is

	main B- channels	main B0 channels	B- through D**	B0 through D**	total B-	total B0
D0 sample	74.0%	0	8.6%	17.4%	82.6%	17.4%
*D sample**	0	80.7%	13.0%	6.3%	13.0%	87%

D** decays in EvtGen/QQ tables, see

/d0dist/dist/releases/<version>/EvtGen/DECAY.DEC

decays contributing to D0 or D*- cases are highlighted

Decay D_0*+
0.3333   D+ pi0                         PHSP;
0.6667   D0 pi+                         PHSP;

Decay D_0*0
0.3333   D0 pi0                         PHSP;
0.6667   D+ pi-                         PHSP;

Decay D_1+
0.3333    D*+ pi0                        VVS_PWAVE 0.0 0.0 0.0 0.0 1.0 0.0;
0.6667    D*0 pi+                        VVS_PWAVE 0.0 0.0 0.0 0.0 1.0 0.0;

Decay D_10
0.3333    D*0 pi0                        VVS_PWAVE 0.0 0.0 0.0 0.0 1.0 0.0;
0.6667    D*+ pi-                        VVS_PWAVE 0.0 0.0 0.0 0.0 1.0 0.0;

Decay D'_1+
0.3333    D*+ pi0                        VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
0.6667    D*0 pi+                        VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;

Decay D'_10
0.6667    D*+ pi-                        VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
0.3333    D*0 pi0                        VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;

Decay D_2*+
0.1030    D*+ pi0                        TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
0.2090    D*0 pi+                        TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
0.2290    D+ pi0                        TSS;
0.4590    D0 pi+                        TSS;

Decay D_2*0
0.2090    D*+ pi-                        TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
0.1030    D*0 pi0                        TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
0.2290    D0 pi0                        TSS;
0.4590    D+ pi-                        TSS;