forked from cms-externals/evtgen
-
Notifications
You must be signed in to change notification settings - Fork 0
/
DECAY_2010.DEC
10267 lines (10155 loc) · 524 KB
/
DECAY_2010.DEC
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
########################################################################
# Copyright 1998-2020 CERN for the benefit of the EvtGen authors #
# #
# This file is part of EvtGen. #
# #
# EvtGen is free software: you can redistribute it and/or modify #
# it under the terms of the GNU General Public License as published by #
# the Free Software Foundation, either version 3 of the License, or #
# (at your option) any later version. #
# #
# EvtGen is distributed in the hope that it will be useful, #
# but WITHOUT ANY WARRANTY; without even the implied warranty of #
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the #
# GNU General Public License for more details. #
# #
# You should have received a copy of the GNU General Public License #
# along with EvtGen. If not, see <https://www.gnu.org/licenses/>. #
########################################################################
# Updated to PDG 2010 by Tomas Pilar - T.Pilar@warwick.ac.uk
# Updated Mixing Parameters Mark Whitehead May 2009
#
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 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 K*BR anti-K*0
Alias K*0R 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
ChargeConj K*0R K*BR
#
# 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
# PR LHCb 10/01/2006
# Turn on PHOTOS for all decays
yesPhotos
#
#
#Decay vpho
#1.000 PYCONT;
#Enddecay
#
# use new VSS_BMIX mixing decay model (DK,28-Oct-1999)
Decay Upsilon(4S)
0.515122645 B+ B- VSS; #[Reconstructed PDG2011]
0.483122645 B0 anti-B0 VSS_BMIX dm; #[Reconstructed PDG2011]
0.000015583 e+ e- PHOTOS VLL; #[Reconstructed PDG2011]
0.000015766 mu+ mu- PHOTOS VLL; #[Reconstructed PDG2011]
0.000015766 tau+ tau- PHOTOS VLL; #[Reconstructed PDG2011]
0.000084099 Upsilon(2S) pi+ pi- VVPIPI; #[Reconstructed PDG2011]
0.000044342 Upsilon(2S) pi0 pi0 VVPIPI; #[Reconstructed PDG2011]
0.000080123 Upsilon pi+ pi- VVPIPI; #[Reconstructed PDG2011]
0.000044342 Upsilon pi0 pi0 VVPIPI; #[Reconstructed PDG2011]
0.000194392 Upsilon eta PARTWAVE 0.0 0.0 1.0 0.0 0.0 0.0; #[Reconstructed PDG2011]
# 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.; #[Reconstructed PDG2011]
# V-> gamma V Partial wave (L,S)=(0,1)
0.000138938 gamma chi_b1(3P) HELAMP 1. 0. 1. 0. -1. 0. -1. 0.; #[Reconstructed PDG2011]
# 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.; #[Reconstructed PDG2011]
# V-> gamma S Partial wave (L,S)=(0,0)
0.000002956 gamma chi_b0(2P) HELAMP 1. 0. 1. 0.; #[Reconstructed PDG2011]
# V-> gamma V Partial wave (L,S)=(0,1)
0.000007883 gamma chi_b1(2P) HELAMP 1. 0. 1. 0. -1. 0. -1. 0.; #[Reconstructed PDG2011]
# 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.; #[Reconstructed PDG2011]
0.000837571 g g g PYTHIA 4; #[Reconstructed PDG2011]
0.000039415 gamma g g PYTHIA 4; #[Reconstructed PDG2011]
Enddecay
#
#
#
Decay anti-B0
# Updated to PDG 2008
# b -> c semileptonic
#
0.050100000 D*+ e- anti-nu_e PHOTOS HQET 0.77 1.33 0.92; #[Reconstructed PDG2011]
0.021700000 D+ e- anti-nu_e PHOTOS ISGW2; #[Reconstructed PDG2011]
0.0054 D_1+ e- anti-nu_e PHOTOS ISGW2;
0.0020 D_0*+ e- anti-nu_e PHOTOS ISGW2;
0.0050 D'_1+ e- anti-nu_e PHOTOS ISGW2;
0.0022 D_2*+ e- anti-nu_e PHOTOS ISGW2;
0.0003 D*+ pi0 e- anti-nu_e PHOTOS GOITY_ROBERTS;
0.004900000 D*0 pi+ e- anti-nu_e PHOTOS GOITY_ROBERTS; #[Reconstructed PDG2011]
0.0010 D+ pi0 e- anti-nu_e PHOTOS GOITY_ROBERTS;
0.0000 D0 pi+ e- anti-nu_e PHOTOS GOITY_ROBERTS;
0.050100000 D*+ mu- anti-nu_mu PHOTOS HQET 0.77 1.33 0.92; #[Reconstructed PDG2011]
0.021700000 D+ mu- anti-nu_mu PHOTOS ISGW2; #[Reconstructed PDG2011]
0.0054 D_1+ mu- anti-nu_mu PHOTOS ISGW2;
0.0020 D_0*+ mu- anti-nu_mu PHOTOS ISGW2;
0.0050 D'_1+ mu- anti-nu_mu PHOTOS ISGW2;
0.0022 D_2*+ mu- anti-nu_mu PHOTOS ISGW2;
0.0003 D*+ pi0 mu- anti-nu_mu PHOTOS GOITY_ROBERTS;
0.004900000 D*0 pi+ mu- anti-nu_mu PHOTOS GOITY_ROBERTS; #[Reconstructed PDG2011]
0.0010 D+ pi0 mu- anti-nu_mu PHOTOS GOITY_ROBERTS;
0.0000 D0 pi+ mu- anti-nu_mu PHOTOS GOITY_ROBERTS;
#
# b -> c tau nu
#
0.015000000 D*+ tau- anti-nu_tau ISGW2; #[Reconstructed PDG2011]
0.011000000 D+ tau- anti-nu_tau ISGW2; #[Reconstructed PDG2011]
0.0013 D_1+ tau- anti-nu_tau ISGW2;
0.0013 D_0*+ tau- anti-nu_tau ISGW2;
0.0020 D'_1+ tau- anti-nu_tau ISGW2;
0.0020 D_2*+ tau- anti-nu_tau ISGW2;
#
# b -> u l nu
#
# 03/28/01 S.Menke included inclusive B -> X_u l nu decays for m(X_u) > 1.26 GeV.
# with some thresholds to gradually switch from excl. to incl.
# modeling without big jumps keeping the better known excl.
# Brs at their actual values and constrain the incl. Br and
# weights in a way that ~40% of the total BR belongs to
# the mass region > 1.5 GeV
# October 20, 2004 - Lange - update using sl awg input
# March 30, 2005 J. Dingfelder: Updated B -> X_u l nu hybrid model
# (supersedes older version by S. Menke). Inclusive
# decays are reweighted in three kinematic
# variables mX, q2, and El (=> 512 weights total).
# The list of weights is included as ModelAlias
# once per B charge state.
# NOTE: - The specified parameters for model VUBHYBRID are:
# mb, a, alpha_s, nbins_mX, nbins_q2, nbins_El,
# list of threshold values for mX, q2, El in ascending order.
#
# *** DO NOT CHANGE ANY OF THESE PARAMETERS! ***
# Otherwise a new set of weights would be needed, since they were
# generated for a certain set of values for mb, a, excl. BFs, non-res BF.
#
# - If no binning (nbins or thresholds) are specified after the first
# three parameters (mb,a,alpha_s) no hybrid weighting is performed.
# --- Hybrid weights for neutral B -> Xu l nu decays
# NOTE: Do NOT CHANGE any BFs without using the corresponding set of
# hybrid weights (and vice versa).
0.000134000 pi+ e- anti-nu_e PHOTOS ISGW2; #[Reconstructed PDG2011]
0.000247000 rho+ e- anti-nu_e PHOTOS ISGW2; #[Reconstructed PDG2011]
0.000000 D(2S)+ e- anti-nu_e PHOTOS ISGW2;
0.000000 D*(2S)+ e- anti-nu_e PHOTOS ISGW2;
0.001892 Xu+ e- anti-nu_e PHOTOS VUB 4.8 1.29 0.22 20 0.30 0.55 1.20 0.61 1.26 0.85 1.34 1.08 1.41 1.21 1.48 1.30 1.55 1.30 1.61 1.33 1.67 1.36 1.73 1.39 1.79 1.33 1.84 1.42 1.90 1.39 1.95 1.39 2.00 1.37 2.50 1.30 3.00 0.74 3.50 0.99 4.00 1.09 4.50 1.00;
#
0.000134000 pi+ mu- anti-nu_mu PHOTOS ISGW2; #[Reconstructed PDG2011]
0.000247000 rho+ mu- anti-nu_mu PHOTOS ISGW2; #[Reconstructed PDG2011]
0.000000 D(2S)+ mu- anti-nu_mu PHOTOS ISGW2;
0.000000 D*(2S)+ mu- anti-nu_mu PHOTOS ISGW2;
0.001892 Xu+ mu- anti-nu_mu PHOTOS VUB 4.8 1.29 0.22 20 0.30 0.55 1.20 0.61 1.26 0.85 1.34 1.08 1.41 1.21 1.48 1.30 1.55 1.30 1.61 1.33 1.67 1.36 1.73 1.39 1.79 1.33 1.84 1.42 1.90 1.39 1.95 1.39 2.00 1.37 2.50 1.30 3.00 0.74 3.50 0.99 4.00 1.09 4.50 1.00;
#
0.000060 pi+ tau- anti-nu_tau ISGW2;
0.000083 rho+ tau- anti-nu_tau ISGW2;
0.000090 a_1+ tau- anti-nu_tau ISGW2;
0.000008 b_1+ tau- anti-nu_tau ISGW2;
0.000008 a_0+ tau- anti-nu_tau ISGW2;
0.000008 a_2+ tau- anti-nu_tau ISGW2;
0.000000 rho(2S)+ tau- anti-nu_tau ISGW2;
0.000000 D(2S)+ tau- anti-nu_tau ISGW2;
0.000000 D*(2S)+ tau- anti-nu_tau ISGW2;
#
#
# b->u hadronic
# Ref. [B1]:
# Lange Nov 14, 2004 (flip D_s K -- is opposite of D_s pi)
0.000024000 D_s- pi+ PHSP; #[Reconstructed PDG2011]
0.000030000 D_s+ K- PHSP; #[Reconstructed PDG2011]
0.000021000 D_s*- pi+ SVS; #[Reconstructed PDG2011]
0.000021900 D_s*+ K- SVS; #[Reconstructed PDG2011]
#
0.000016 rho+ D_s- SVS;
0.000035000 K*- D_s+ SVS; #[Reconstructed PDG2011]
0.000041000 D_s*- rho+ SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0; #[Reconstructed PDG2011]
0.000032000 D_s*+ K*- SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0; #[Reconstructed PDG2011]
#
# b -> s gamma
#
0.000043300 anti-K*0 gamma HELAMP 1.0 0.0 1.0 0.0; #[Reconstructed PDG2011]
# PR LHCb add omega gamma
0.000000440 omega gamma HELAMP 1.0 0.0 1.0 0.0; #[Reconstructed PDG2011]
#0.0000135 anti-K_10 gamma HELAMP 1.0 0.0 1.0 0.0;
#0.0000065 anti-K'_10 gamma HELAMP 1.0 0.0 1.0 0.0;
#0.0000128 anti-K'*0 gamma HELAMP 1.0 0.0 1.0 0.0;
#0.0000166 anti-K_2*0 gamma HELAMP 1.0 0.0 1.0 0.0;
#0.0000017 anti-K''*0 gamma HELAMP 1.0 0.0 1.0 0.0;
0.0003118 anti-Xsd gamma BTOXSGAMMA 2 ;
#
0.000000160 anti-K0 e+ e- PHOTOS BTOSLLBALL; #[Reconstructed PDG2011]
0.000001030 anti-K*0 e+ e- PHOTOS BTOSLLBALL; #[Reconstructed PDG2011]
0.0000050 anti-Xsd e+ e- PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
0.000000450 anti-K0 mu+ mu- PHOTOS BTOSLLBALL; #[Reconstructed PDG2011]
0.000001050 anti-K*0 mu+ mu- PHOTOS BTOSLLBALL; #[Reconstructed PDG2011]
0.0000025 anti-Xsd mu+ mu- PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
0.0000001 anti-K0 tau+ tau- PHOTOS BTOSLLBALL;
0.0000002 anti-K*0 tau+ tau- PHOTOS BTOSLLBALL;
0.0000002 anti-Xsd tau+ tau- PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
#
#------------------------------------------------------------------------------
# 2- and 3-body modes revised Feb.2005 Markus Cristinziani, SLAC
#vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
# B0 -> CP eigenstates (some exclusive b -> u) sum=0.00036
#
0.000005130 pi+ pi- PHSP; #[Reconstructed PDG2011]
0.000001620 pi0 pi0 PHSP; #[Reconstructed PDG2011]
0.0000010 pi0 eta PHSP;
0.000001200 pi0 eta' PHSP; #[Reconstructed PDG2011]
0.000001 pi0 a_00 PHSP;
0.000001 pi0 f_0 PHSP;
# pi0 rho0 is with the 3-body modes
0.000001 omega pi0 SVS;
0.000001 a_10 pi0 SVS;
0.000001 b_10 pi0 SVS;
0.000001 eta eta PHSP;
0.000001 eta eta' PHSP;
0.000001 eta a_00 PHSP;
0.000001 eta f_0 PHSP;
0.000001 rho0 eta SVS;
0.000000940 omega eta SVS; #[Reconstructed PDG2011]
0.000001 a_10 eta SVS;
0.000001 b_10 eta SVS;
0.000001 eta' eta' PHSP;
0.000001 eta' a_00 PHSP;
0.000001 eta' f_0 PHSP;
0.000001 rho0 eta' SVS;
0.000001000 omega eta' SVS; #[Reconstructed PDG2011]
0.000001 a_10 eta' SVS;
0.000001 b_10 eta' SVS;
0.000001 a_00 a_00 PHSP;
0.000001 a_00 f_0 PHSP;
0.000001 rho0 a_00 SVS;
0.000001 omega a_00 SVS;
0.000001 a_10 a_00 SVS;
0.000001 b_10 a_00 SVS;
0.000001 f_0 f_0 PHSP;
0.000001 rho0 f_0 SVS;
0.000001 omega f_0 SVS;
0.000001 a_10 f_0 SVS;
0.000001 b_10 f_0 SVS;
# Penguin dominated modes are sin2beta, not sin2alpha
0.000004300 phi K_S0 SVS; #[Reconstructed PDG2011]
0.000004300 phi K_L0 SVS; #[Reconstructed PDG2011]
0.000000550 eta K_S0 PHSP; #[Reconstructed PDG2011]
0.000000550 eta K_L0 PHSP; #[Reconstructed PDG2011]
0.000033000 eta' K_S0 PHSP; #[Reconstructed PDG2011]
0.000033000 eta' K_L0 PHSP; #[Reconstructed PDG2011]
0.000002500 omega K_S0 SVS; #[Reconstructed PDG2011]
0.000002500 omega K_L0 SVS; #[Reconstructed PDG2011]
#don't allow jetset to simulate the same decay but to a K0 or anti-K0
#Lange 1/30/2003
0.0 phi K0 PHSP;
0.0 phi anti-K0 PHSP;
0.0 eta K0 PHSP;
0.0 eta anti-K0 PHSP;
0.0 eta' K0 PHSP;
0.0 eta' anti-K0 PHSP;
0.0 omega K0 PHSP;
0.0 omega anti-K0 PHSP;
# phipi0 has no simple Pengiun or tree so who knows - ditto phieta/eta'?
0.000001 phi pi0 SVS;
0.000001 phi eta SVS;
0.000001 phi eta' SVS;
#
0.000015900 anti-K*0 eta SVS; #[Reconstructed PDG2011]
0.000003800 anti-K*0 eta' SVS; #[Reconstructed PDG2011]
0.000002000 anti-K*0 omega SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0; #[Reconstructed PDG2011]
0.000001 phi omega SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.000001 omega omega SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.000001 omega rho0 SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
#
#
#2-body modes from Jim Olsen
0.000019400 K- pi+ PHSP; #[Reconstructed PDG2011]
0.0000004 K+ K- PHSP;
# K_S0 K_L0 is 0
# see http://babar-hn.slac.stanford.edu:5090/HyperNews/get/pubboard14/156.html
0.000000240 K_S0 K_S0 PHSP; #[Reconstructed PDG2011]
0.0000000 K_S0 K_L0 PHSP;
0.000000240 K_L0 K_L0 PHSP; #[Reconstructed PDG2011]
# declare K0 pi0 not to have it through Pythia
0.000000 anti-K0 pi0 PHSP;
# Split into K_S0 K_L0 (LHCb - P. Robbe)
0.000004750 K_S0 pi0 PHSP; #[Reconstructed PDG2011]
0.000004750 K_L0 pi0 PHSP; #[Reconstructed PDG2011]
# 0 to keep PYTHIA from generating it
0.000000 K0 anti-K0 PHSP;
#
# 3-body decays
# this model generates interferences and higher mass rho
#
# 3-body by John Back (jback@slac.stanford.edu) - Oct 15, 2002
# JGS intersperses modes with pi0->eta,eta'
# anti-B0 modes
# Model generates rho(770), rho(1450) and rho(1700) resonances
# with interference (uses CKM angle alpha)
0.000024 pi- pi+ pi0 BTO3PI_CP dm alpha;
# Put in with 0 BR to keep PYTHIA from generating these (see BTO3PI_CP above)
0.000000 rho0 pi0 SVS;
0.000000 rho- pi+ SVS;
0.000000 rho+ pi- SVS;
0.000000 rho(2S)- pi+ SVS;
0.000000 rho(2S)+ pi- SVS;
0.000000 rho(3S)- pi+ SVS;
0.000000 rho(3S)+ pi- SVS;
#
# rho0 3-body modes
0.000002350 rho0 K_S0 SVS; #[Reconstructed PDG2011]
0.000002350 rho0 K_L0 SVS; #[Reconstructed PDG2011]
#
# rho+ 3-body modes
0.000008400 rho+ K- SVS; #[Reconstructed PDG2011]
#
# rho(1450) 3-body modes
0.0000005 rho(2S)0 K_S0 SVS;
0.0000005 rho(2S)0 K_L0 SVS;
#
#
# f0,a0(980) 3-body modes
0.0000055 f_0 anti-K0 PHSP;
0.000001 a_00 anti-K0 PHSP;
0.000001 a_0+ K- PHSP;
0.000003 a_0+ pi- PHSP;
0.000001 a_0- pi+ PHSP;
#
# K*0(892) 3-body modes
0.000003600 anti-K*0 pi0 SVS; #[Reconstructed PDG2011]
0.000003 anti-K*0 K0 SVS;
#
# K*0(1430) 3-body modes
0.000001 anti-K_0*0 pi0 PHSP;
0.000003 anti-K_0*0 K0 PHSP;
#
# K*0(1680) 3-body modes
0.000001 anti-K''*0 pi0 SVS;
0.000001 anti-K''*0 K0 SVS;
#
# K*-(892) 3-body modes
0.000009400 K*- pi+ SVS; #[Reconstructed PDG2011]
0.000002 K*- K+ SVS;
#
# K*-(1430) 3-body modes
0.000033000 K_0*- pi+ PHSP; #[Reconstructed PDG2011]
0.000002 K_0*- K+ PHSP;
#
# K*0(1680) 3-body modes
0.000001 K''*- pi+ SVS;
0.000001 K''*- K+ SVS;
#
#
# Non-resonant 3-body left-overs
# Most modes are set at 1e-6
#
# pi-pi+pi0: high mass + f0(400-1200) = (1+1)e-6
####0402270.000002 pi- pi+ pi0 PHSP;
0.000017 pi- pi+ eta PHSP;
0.000001 pi- pi+ eta' PHSP;
#
# pi- pi+ anti-K0: total-rhoK0-f0K0-piK*(892)-piK*(1430)=45-5-6-13./3-11./3
0.000000000 pi- pi+ anti-K0 PHSP; #[Reconstructed PDG2011]
#
# K- pi+ pi0
0.000027500 K- pi+ pi0 PHSP; #[Reconstructed PDG2011]
0.000020 K- pi+ eta PHSP;
0.000001 K- pi+ eta' PHSP;
#
# pi- K+ K0.
# Both B0 and B0bar can decay to these => 0.5e-6 each
0.0000005 pi- K+ anti-K0 PHSP;
0.0000005 pi+ K- K0 PHSP;
#
# K- K+ anti-K0
# PR LHCb 09 Apr 2004 split into KL/KS
0.000000 K- K+ anti-K0 PHSP;
0.000012350 K- K+ K_S0 PHSP; #[Reconstructed PDG2011]
0.000012350 K- K+ K_L0 PHSP; #[Reconstructed PDG2011]
#
# K- K+ pi0
0.000001 K- K+ pi0 PHSP;
0.000001 K- K+ eta PHSP;
0.000001 K- K+ eta' PHSP;
#
# pi0 anti-K0 K0
0.000002 anti-K0 K0 pi0 PHSP;
0.000001 anti-K0 K0 eta PHSP;
0.000001 anti-K0 K0 eta' PHSP;
#
# pi0 pi0 anti-K0
0.000002 anti-K0 pi0 pi0 PHSP;
0.000001 anti-K0 pi0 eta PHSP;
0.000001 anti-K0 pi0 eta' PHSP;
0.000001 anti-K0 eta eta PHSP;
0.000001 anti-K0 eta eta' PHSP;
0.000001 anti-K0 eta' eta' PHSP;
#
# anti-K0 K0 anti-K0: 8*BF(KsKsKs)~=48 (8 is a guess)
0.000048 anti-K0 K0 anti-K0 PHSP;
#
# pi0 pi0 pi0
0.000001 pi0 pi0 pi0 PHSP;
0.000001 pi0 pi0 eta PHSP;
0.000001 pi0 pi0 eta' PHSP;
0.000001 pi0 eta eta PHSP;
0.000001 pi0 eta eta' PHSP;
#
#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
# 2- and 3-body modes revised Feb.2005 Markus Cristinziani, SLAC
#---------------------------------------------------------------------------
#
#4-body modes from Andrei Gritsan
#--- 4-body rho-rho, rho-pi-pi, pi-pi-pi-pi --------------------------
# PR LHCb 21 Apr 2004 : set polarisation to longitudinal
0.000000730 rho0 rho0 SVV_HELAMP 0.0 0.0 1.0 0.0 0.0 0.0; #[Reconstructed PDG2011]
0.000024200 rho+ rho- SVV_HELAMP 0.0 0.0 1.0 0.0 0.0 0.0; #[Reconstructed PDG2011]
0.000003 rho0 pi+ pi- PHSP;
0.000003 omega pi+ pi- PHSP;
0.000001 rho0 pi0 pi0 PHSP;
0.000001 omega pi0 pi0 PHSP;
0.0000005 rho0 pi0 eta PHSP;
0.0000002 rho0 eta eta PHSP;
0.0000002 rho0 pi0 eta' PHSP;
0.0000001 rho0 eta eta' PHSP;
0.0000005 omega pi0 eta PHSP;
0.0000002 omega eta eta PHSP;
0.0000002 omega pi0 eta' PHSP;
0.0000001 omega eta eta' PHSP;
0.000010 rho+ pi- pi0 PHSP;
0.000005 rho+ pi- eta PHSP;
0.000010 rho- pi+ pi0 PHSP;
0.000005 rho- pi+ eta PHSP;
0.000010 pi+ pi- pi+ pi- PHSP;
0.000010 pi+ pi- pi0 pi0 PHSP;
0.000005 pi+ pi- eta pi0 PHSP;
0.000002 pi+ pi- eta eta PHSP;
0.000002 pi+ pi- eta' pi0 PHSP;
0.000001 pi+ pi- eta' eta PHSP;
0.000016500 a_1+ pi- SVS; #[Reconstructed PDG2011]
0.000016500 a_1- pi+ SVS; #[Reconstructed PDG2011]
0.0000055 b_1+ pi- SVS;
0.0000055 b_1- pi+ SVS;
0.000002 rho- a_0+ SVS;
0.000010 rho+ a_0- SVS;
0.000001 f_0 pi+ pi- PHSP;
0.000001 f_0 pi0 pi0 PHSP;
0.000001 a_0- pi+ pi0 PHSP;
0.000001 a_00 pi0 pi0 PHSP;
0.000001 a_0+ pi- pi0 PHSP;
#--- 4-body rho-K*, rho-pi-K, K*-pi-pi, pi-pi-pi-K -------------------
0.000003400 rho0 anti-K*0 SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0; #[Reconstructed PDG2011]
0.000010 rho+ K*- SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.000050600 pi- pi+ anti-K*0 PHSP; #[Reconstructed PDG2011]
0.000010 pi+ pi0 K*- PHSP;
0.000005 pi+ eta K*- PHSP;
0.000002 pi+ eta' K*- PHSP;
0.000010 rho+ K- pi0 PHSP;
0.000005 rho+ K- eta PHSP;
0.000002 rho+ K- eta' PHSP;
0.000010 rho+ anti-K0 pi- PHSP;
0.000010 rho- anti-K0 pi+ PHSP;
0.000005 rho0 anti-K0 pi0 PHSP;
0.000002 rho0 anti-K0 eta PHSP;
0.000001 rho0 anti-K0 eta' PHSP;
0.000010 rho0 K- pi+ PHSP;
0.000005 omega anti-K0 pi0 PHSP;
0.000002 omega anti-K0 eta PHSP;
0.000001 omega anti-K0 eta' PHSP;
0.000005100 omega K- pi+ PHSP; #[Reconstructed PDG2011]
0.000010 pi+ pi0 pi- anti-K0 PHSP;
0.000005 pi+ eta pi- anti-K0 PHSP;
0.000010 pi0 pi0 pi0 anti-K0 PHSP;
0.000005 pi0 pi0 eta anti-K0 PHSP;
0.000002 pi0 eta eta anti-K0 PHSP;
0.000002 pi0 pi0 eta' anti-K0 PHSP;
0.000001 pi0 eta eta' anti-K0 PHSP;
0.000010 pi+ pi- pi+ K- PHSP;
0.000010 pi0 pi0 pi+ K- PHSP;
0.000005 pi0 eta pi+ K- PHSP;
0.000002 eta eta pi+ K- PHSP;
0.000002 pi0 eta' pi+ K- PHSP;
0.000001 eta eta' pi+ K- PHSP;
0.000010 rho0 anti-K_0*0 PHSP;
0.000020 rho+ K_0*- PHSP;
0.000010 pi- pi+ anti-K_0*0 PHSP;
0.000010 pi0 pi0 anti-K_0*0 PHSP;
0.000010 pi+ pi0 K_0*- PHSP;
0.000005 anti-K*0 f_0 SVS;
0.000010 a_10 anti-K0 SVS;
0.000016000 a_1+ K- SVS; #[Reconstructed PDG2011]
0.000010 b_10 anti-K0 SVS;
0.0000074 b_1+ K- SVS;
0.000005 anti-K*0 a_00 SVS;
0.000005 K*- a_0+ SVS;
#--- 4-body K*-K*, rho-K-K, K*-K-pi, phi-pi-pi, pi-pi-K-K ------------
0.000001280 K*0 anti-K*0 SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0; #[Reconstructed PDG2011]
0.000001 K*+ K*- SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.0000001 phi rho0 SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.0000001 phi pi- pi+ PHSP;
0.000001 K+ K- rho0 PHSP;
0.000001 K0 anti-K0 rho0 PHSP;
0.000001 K+ K- omega PHSP;
0.000001 K0 anti-K0 omega PHSP;
0.000001 K+ anti-K0 rho- PHSP;
0.000001 K- K0 rho+ PHSP;
0.000001 K+ K- pi- pi+ PHSP;
0.000001 K+ K- pi0 pi0 PHSP;
0.0000005 K+ K- pi0 eta PHSP;
0.0000002 K+ K- eta eta PHSP;
0.0000002 K+ K- pi0 eta' PHSP;
0.0000001 K+ K- eta eta' PHSP;
0.000001 K0 anti-K0 pi- pi+ PHSP;
0.000001 K0 anti-K0 pi0 pi0 PHSP;
0.0000005 K0 anti-K0 pi0 eta PHSP;
0.0000002 K0 anti-K0 eta eta PHSP;
0.0000002 K0 anti-K0 pi0 eta' PHSP;
0.0000001 K0 anti-K0 eta eta' PHSP;
0.000001 K*+ K- pi0 PHSP;
0.000001 K*- K+ pi0 PHSP;
0.000001 K*0 anti-K0 pi0 PHSP;
0.000001 K0 anti-K*0 pi0 PHSP;
0.000003320 K+ anti-K*0 pi- PHSP; #[Reconstructed PDG2011]
0.000001 K*+ anti-K0 pi- PHSP;
0.0000033 K- K*0 pi+ PHSP;
0.000001 K*- K0 pi+ PHSP;
#--- 4-body phi-K*, phi-K-pi, K-K-K*, pi-K-K-K -----------------------
0.000009800 phi anti-K*0 SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0; #[Reconstructed PDG2011]
0.000010 phi K- pi+ PHSP;
0.000010 phi anti-K0 pi0 PHSP;
0.000007 phi anti-K0 eta PHSP;
0.000017700 K- K+ anti-K*0 PHSP; #[Reconstructed PDG2011]
0.000010 K- K*+ anti-K0 PHSP;
0.000010 K*- K+ anti-K0 PHSP;
0.000010 K+ K- K- pi+ PHSP;
0.000010 K+ K- anti-K0 pi0 PHSP;
0.000005 K+ K- anti-K0 eta PHSP;
0.000010 K0 anti-K0 anti-K0 pi0 PHSP;
0.000005 K0 anti-K0 anti-K0 eta PHSP;
0.000003900 phi anti-K_0*0 PHSP; #[Reconstructed PDG2011]
0.000010 K- K+ anti-K_0*0 PHSP;
0.000010 K- K_0*+ anti-K0 PHSP;
0.000010 K_0*- K+ anti-K0 PHSP;
0.000002 phi anti-K'_10 PHSP;
0.000002 K- K+ anti-K'_10 PHSP;
#--- 4-body phi-phi, phi-K-K, K-K-K-K --------------------------------
0.0000001 phi phi SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
0.0000001 phi K+ K- PHSP;
0.0000001 phi K0 anti-K0 PHSP;
0.0000001 K+ K- K+ K- PHSP;
0.0000001 K0 anti-K0 K+ K- PHSP;
0.0000001 K0 anti-K0 K0 anti-K0 PHSP;
#5-body modes from Jim Smith - 1/30/05
#--- 5-body rho-rho-pi, rho-pi-pi-pi, pi-pi-pi-pi-pi -------------------
0.000010 rho0 rho0 pi0 PHSP;
0.000005 rho0 rho0 eta PHSP;
0.000002 rho0 rho0 eta' PHSP;
0.000010 omega rho0 pi0 PHSP;
0.000005 omega rho0 eta PHSP;
0.000002 omega rho0 eta' PHSP;
0.000010 rho+ rho- pi0 PHSP;
0.000005 rho+ rho- eta PHSP;
0.000002 rho+ rho- eta' PHSP;
0.000010 rho0 rho+ pi- PHSP;
0.000010 omega rho+ pi- PHSP;
0.000002 rho0 pi+ pi- pi0 PHSP;
0.000001 rho0 pi+ pi- eta PHSP;
0.000001 rho0 pi+ pi- eta' PHSP;
0.000002 rho0 pi0 pi0 pi0 PHSP;
0.000001 rho0 pi0 eta pi0 PHSP;
0.000001 rho0 pi0 eta' pi0 PHSP;
0.000002 omega pi+ pi- pi0 PHSP;
0.000001 omega pi+ pi- eta PHSP;
0.000001 omega pi+ pi- eta' PHSP;
0.000002 omega pi0 pi0 pi0 PHSP;
0.000001 omega pi0 eta pi0 PHSP;
0.000001 omega pi0 eta' pi0 PHSP;
0.000002 rho+ pi- pi+ pi- PHSP;
0.000002 rho+ pi- pi0 pi0 PHSP;
0.000001 rho+ pi- pi0 eta PHSP;
0.000001 rho+ pi- pi0 eta' PHSP;
0.000002 rho- pi+ pi+ pi- PHSP;
0.000002 rho- pi+ pi0 pi0 PHSP;
0.000001 rho- pi+ pi0 eta PHSP;
0.000001 rho- pi+ pi0 eta' PHSP;
0.000001 pi+ pi- pi+ pi- pi0 PHSP;
0.0000005 pi+ pi- pi+ pi- eta PHSP;
0.0000002 pi+ pi- pi+ pi- eta' PHSP;
0.000001 pi+ pi- pi0 pi0 pi0 PHSP;
0.0000005 pi+ pi- pi0 pi0 eta PHSP;
0.0000002 pi+ pi- pi0 pi0 eta' PHSP;
0.0000005 pi0 pi0 pi0 pi0 pi0 PHSP;
0.0000002 pi0 pi0 pi0 pi0 eta PHSP;
0.0000001 pi0 pi0 pi0 pi0 eta' PHSP;
0.000010 a_1+ rho- PHSP;
0.000010 a_1- rho+ PHSP;
0.000010 a_10 rho0 PHSP;
0.000010 a_10 omega PHSP;
0.000010 a_1+ pi- pi0 PHSP;
0.000005 a_1+ pi- eta PHSP;
0.000002 a_1+ pi- eta' PHSP;
0.000010 a_10 pi+ pi- PHSP;
0.000010 a_10 pi0 pi0 PHSP;
0.000005 a_10 pi0 eta PHSP;
0.000002 a_10 pi0 eta' PHSP;
0.000010 a_1- pi+ pi0 PHSP;
0.000005 a_1- pi+ eta PHSP;
0.000002 a_1- pi+ eta' PHSP;
0.000010 rho+ f_0 pi- PHSP;
0.000010 rho- f_0 pi+ PHSP;
0.000010 rho0 f_0 pi0 PHSP;
0.000005 rho0 f_0 eta PHSP;
0.000002 rho0 f_0 eta' PHSP;
0.000010 omega f_0 pi0 PHSP;
0.000010 rho+ a_00 pi- PHSP;
0.000010 rho+ a_0- pi0 PHSP;
0.000005 rho+ a_0- eta PHSP;
0.000002 rho+ a_0- eta' PHSP;
0.000010 rho0 a_0+ pi- PHSP;
0.000010 rho0 a_0- pi+ PHSP;
0.000010 rho0 a_00 pi0 PHSP;
0.000005 rho0 a_00 eta PHSP;
0.000002 rho0 a_00 eta' PHSP;
0.000010 rho- a_00 pi+ PHSP;
0.000010 rho- a_0+ pi0 PHSP;
0.000005 rho- a_0+ eta PHSP;
0.000002 rho- a_0+ eta' PHSP;
0.000002 f_0 pi+ pi- pi0 PHSP;
0.000001 f_0 pi+ pi- eta PHSP;
0.000001 f_0 pi+ pi- eta' PHSP;
0.000002 f_0 pi0 pi0 pi0 PHSP;
0.000001 f_0 pi0 pi0 eta PHSP;
0.000001 f_0 pi0 pi0 eta' PHSP;
0.000001 a_00 pi+ pi- pi0 PHSP;
0.000001 a_00 pi0 pi0 pi0 PHSP;
0.000001 a_0+ pi- pi+ pi- PHSP;
0.000001 a_0+ pi- pi0 pi0 PHSP;
0.000001 a_0- pi+ pi+ pi- PHSP;
0.000001 a_0- pi+ pi0 pi0 PHSP;
#--- 5-body rho-K*pi, rho-2pi-K, K*-3pi, 4pi-K -------------------
0.000010 rho0 K*- pi+ PHSP;
0.000010 rho+ K*- pi0 PHSP;
0.000005 rho+ K*- eta PHSP;
0.000002 rho+ K*- eta' PHSP;
0.000010 rho0 anti-K*0 pi0 PHSP;
0.000005 rho0 anti-K*0 eta PHSP;
0.000002 rho0 anti-K*0 eta' PHSP;
0.000010 rho- anti-K*0 pi+ PHSP;
0.000010 rho+ anti-K*0 pi- PHSP;
0.000010 omega K*- pi+ PHSP;
0.000010 omega anti-K*0 pi0 PHSP;
0.000005 omega anti-K*0 eta PHSP;
0.000002 omega anti-K*0 eta' PHSP;
0.000010 pi+ pi- K*- pi+ PHSP;
0.000010 pi+ pi0 K*- pi0 PHSP;
0.000005 pi+ pi0 K*- eta PHSP;
0.000002 pi+ pi0 K*- eta' PHSP;
0.000010 pi+ pi- anti-K*0 pi0 PHSP;
0.000005 pi+ pi- anti-K*0 eta PHSP;
0.000002 pi+ pi- anti-K*0 eta' PHSP;
0.000010 pi0 pi0 anti-K*0 pi0 PHSP;
0.000005 pi0 pi0 anti-K*0 eta PHSP;
0.000002 pi0 pi0 anti-K*0 eta' PHSP;
0.000010 rho- K- pi+ pi+ PHSP;
0.000010 rho0 K- pi+ pi0 PHSP;
0.000005 rho0 K- pi+ eta PHSP;
0.000002 rho0 K- pi+ eta' PHSP;
0.000010 rho+ K- pi+ pi- PHSP;
0.000010 rho+ K- pi0 pi0 PHSP;
0.000005 rho+ K- pi0 eta PHSP;
0.000002 rho+ K- pi0 eta' PHSP;
0.000010 rho0 anti-K0 pi+ pi- PHSP;
0.000010 rho0 anti-K0 pi0 pi0 PHSP;
0.000005 rho0 anti-K0 pi0 eta PHSP;
0.000002 rho0 anti-K0 pi0 eta' PHSP;
0.000010 rho- anti-K0 pi+ pi0 PHSP;
0.000005 rho- anti-K0 pi+ eta PHSP;
0.000002 rho- anti-K0 pi+ eta' PHSP;
0.000010 rho+ anti-K0 pi- pi0 PHSP;
0.000005 rho+ anti-K0 pi- eta PHSP;
0.000002 rho+ anti-K0 pi- eta' PHSP;
0.000010 omega K- pi+ pi0 PHSP;
0.000005 omega K- pi+ eta PHSP;
0.000002 omega K- pi+ eta' PHSP;
0.000010 omega anti-K0 pi+ pi- PHSP;
0.000010 omega anti-K0 pi0 pi0 PHSP;
0.000005 omega anti-K0 pi0 eta PHSP;
0.000002 omega anti-K0 pi0 eta' PHSP;
0.000002 pi+ pi- pi+ K- pi0 PHSP;
0.000001 pi+ pi- pi+ K- eta PHSP;
0.000001 pi+ pi- pi+ K- eta' PHSP;
0.000002 pi0 pi0 pi+ K- pi0 PHSP;
0.000001 pi0 pi0 pi+ K- eta PHSP;
0.000001 pi0 pi0 pi+ K- eta' PHSP;
0.000002 pi+ pi- pi+ anti-K0 pi- PHSP;
0.000002 pi+ pi- pi0 anti-K0 pi0 PHSP;
0.000001 pi+ pi- pi0 anti-K0 eta PHSP;
0.000001 pi+ pi- pi0 anti-K0 eta' PHSP;
0.000002 pi0 pi0 pi0 anti-K0 pi0 PHSP;
0.000001 pi0 pi0 pi0 anti-K0 eta PHSP;
0.000001 pi0 pi0 pi0 anti-K0 eta' PHSP;
0.000010 rho+ anti-K_0*0 pi- PHSP;
0.000010 rho0 anti-K_0*0 pi0 PHSP;
0.000005 rho0 anti-K_0*0 eta PHSP;
0.000002 rho0 anti-K_0*0 eta' PHSP;
0.000010 rho- anti-K_0*0 pi+ PHSP;
0.000010 rho+ K_0*- pi0 PHSP;
0.000005 rho+ K_0*- eta PHSP;
0.000002 rho+ K_0*- eta' PHSP;
0.000010 rho0 K_0*- pi+ PHSP;
0.000010 pi+ pi- K_0*- pi+ PHSP;
0.000010 pi+ pi0 K_0*- pi0 PHSP;
0.000005 pi+ pi0 K_0*- eta PHSP;
0.000002 pi+ pi0 K_0*- eta' PHSP;
0.000010 pi+ pi- anti-K_0*0 pi0 PHSP;
0.000005 pi+ pi- anti-K_0*0 eta PHSP;
0.000002 pi+ pi- anti-K_0*0 eta' PHSP;
0.000010 pi0 pi0 anti-K_0*0 pi0 PHSP;
0.000005 pi0 pi0 anti-K_0*0 eta PHSP;
0.000005 pi0 pi0 anti-K_0*0 eta' PHSP;
0.000010 K*- f_0 pi+ PHSP;
0.000010 anti-K*0 f_0 pi0 PHSP;
0.000005 anti-K*0 f_0 eta PHSP;
0.000002 anti-K*0 f_0 eta' PHSP;
0.000010 a_1+ anti-K0 pi- PHSP;
0.000020 a_10 anti-K*0 PHSP;
0.000010 a_10 K- pi+ PHSP;
0.000010 a_10 anti-K0 pi0 PHSP;
0.000020 a_1+ K*- PHSP;
0.000010 a_1+ K- pi0 PHSP;
0.000010 a_1- anti-K0 pi+ PHSP;
0.000005 K*- a_00 pi+ PHSP;
0.000005 K*- a_0+ pi0 PHSP;
0.000002 K*- a_0+ eta PHSP;
0.000001 K*- a_0+ eta' PHSP;
0.000005 anti-K*0 a_0+ pi- PHSP;
0.000005 anti-K*0 a_00 pi0 PHSP;
0.000002 anti-K*0 a_00 eta PHSP;
0.000001 anti-K*0 a_00 eta' PHSP;
#--- 5-body K*-K*pi,rho-K-Kpi, K*-K-2pi, phi-3pi, 3pi-K-K ------------
0.000001 K*- K*0 pi+ PHSP;
0.000001 K*+ anti-K*0 pi- PHSP;
0.000001 K*- K*+ pi0 PHSP;
0.0000005 K*- K*+ eta PHSP;
0.0000002 K*- K*+ eta' PHSP;
0.000001 K*0 anti-K*0 pi0 PHSP;
0.0000005 K*0 anti-K*0 eta PHSP;
0.0000002 K*0 anti-K*0 eta' PHSP;
0.000001 phi rho+ pi- PHSP;
0.000001 phi rho- pi+ PHSP;
0.000001 phi rho0 pi0 PHSP;
0.0000005 phi rho0 eta PHSP;
0.0000002 phi rho0 eta' PHSP;
0.000001 phi omega pi0 PHSP;
0.000001 phi f_0 pi0 PHSP;
0.000001 phi pi+ pi- pi0 PHSP;
0.0000005 phi pi+ pi- eta PHSP;
0.0000002 phi pi+ pi- eta' PHSP;
0.000001 phi pi0 pi0 pi0 PHSP;
0.0000005 phi pi0 pi0 eta PHSP;
0.0000002 phi pi0 pi0 eta' PHSP;
0.000001 K0 K- rho0 pi+ PHSP;
0.000001 K0 K- rho+ pi0 PHSP;
0.0000005 K0 K- rho+ eta PHSP;
0.0000002 K0 K- rho+ eta' PHSP;
0.000001 K0 K- omega pi+ PHSP;
0.000001 K0 K- f_0 pi+ PHSP;
0.000001 K+ K- rho+ pi- PHSP;
0.000001 K+ K- rho- pi+ PHSP;
0.000001 K+ K- rho0 pi0 PHSP;
0.0000005 K+ K- rho0 eta PHSP;
0.0000002 K+ K- rho0 eta' PHSP;
0.000001 K+ K- omega pi0 PHSP;
0.0000005 K+ K- omega eta PHSP;
0.0000002 K+ K- omega eta' PHSP;
0.000001 K+ K- f_0 pi0 PHSP;
0.000001 K0 anti-K0 rho+ pi- PHSP;
0.000001 K0 anti-K0 rho- pi+ PHSP;
0.000001 K0 anti-K0 rho0 pi0 PHSP;
0.0000005 K0 anti-K0 rho0 eta PHSP;
0.0000002 K0 anti-K0 rho0 eta' PHSP;
0.000001 K0 anti-K0 omega pi0 PHSP;
0.000001 K0 anti-K0 f_0 pi0 PHSP;
0.000001 anti-K0 K+ rho0 pi- PHSP;
0.000001 anti-K0 K+ rho- pi0 PHSP;
0.0000005 anti-K0 K+ rho- eta PHSP;
0.0000002 anti-K0 K+ rho- eta' PHSP;
0.000001 anti-K0 K+ omega pi- PHSP;
0.000001 anti-K0 K+ f_0 pi- PHSP;
0.0000002 K0 K- pi+ pi+ pi- PHSP;
0.0000002 K0 K- pi+ pi0 pi0 PHSP;
0.0000001 K0 K- pi+ pi0 eta PHSP;
0.0000001 K0 K- pi+ pi0 eta' PHSP;
0.0000002 K+ K- pi+ pi- pi0 PHSP;
0.0000001 K+ K- pi+ pi- eta PHSP;
0.0000001 K+ K- pi+ pi- eta' PHSP;
0.0000002 K+ K- pi0 pi0 pi0 PHSP;
0.0000001 K+ K- pi0 pi0 eta PHSP;
0.0000001 K+ K- pi0 pi0 eta' PHSP;
0.0000002 K0 anti-K0 pi+ pi- pi0 PHSP;
0.0000001 K0 anti-K0 pi+ pi- eta PHSP;
0.0000001 K0 anti-K0 pi+ pi- eta' PHSP;
0.0000002 K0 anti-K0 pi0 pi0 pi0 PHSP;
0.0000001 K0 anti-K0 pi0 pi0 eta PHSP;
0.0000001 K0 anti-K0 pi0 pi0 eta' PHSP;
0.0000002 anti-K0 K+ pi- pi+ pi- PHSP;
0.0000002 anti-K0 K+ pi- pi0 pi0 PHSP;
0.0000001 anti-K0 K+ pi- pi0 eta PHSP;
0.0000001 anti-K0 K+ pi- pi0 eta' PHSP;
0.000001 K*0 K- pi+ pi0 PHSP;
0.0000005 K*0 K- pi+ eta PHSP;
0.0000002 K*0 K- pi+ eta' PHSP;
0.000001 K*- K0 pi+ pi0 PHSP;
0.0000005 K*- K0 pi+ eta PHSP;
0.0000002 K*- K0 pi+ eta' PHSP;
0.000001 K*+ K- pi+ pi- PHSP;
0.000001 K*+ K- pi0 pi0 PHSP;
0.0000005 K*+ K- pi0 eta PHSP;
0.0000002 K*+ K- pi0 eta' PHSP;
0.000001 K*- K+ pi+ pi- PHSP;
0.000001 K*- K+ pi0 pi0 PHSP;
0.0000005 K*- K+ pi0 eta PHSP;
0.0000002 K*- K+ pi0 eta' PHSP;
0.000001 K*0 anti-K0 pi+ pi- PHSP;
0.000001 K*0 anti-K0 pi0 pi0 PHSP;
0.0000005 K*0 anti-K0 pi0 eta PHSP;
0.0000002 K*0 anti-K0 pi0 eta' PHSP;
0.000001 K0 anti-K*0 pi+ pi- PHSP;
0.000001 K0 anti-K*0 pi0 pi0 PHSP;