1 import FWCore.ParameterSet.Config
as cms
4 generator = cms.EDFilter(
"Pythia6GeneratorFilter",
6 pythiaPylistVerbosity = cms.untracked.int32(0),
7 filterEfficiency = cms.untracked.double(0.00013),
8 pythiaHepMCVerbosity = cms.untracked.bool(
False),
9 crossSection = cms.untracked.double(54710000000.0),
10 maxEventsToPrint = cms.untracked.int32(0),
11 comEnergy = cms.double(7000.0),
12 PythiaParameters = cms.PSet(
13 pythiaUESettingsBlock,
14 processParameters = cms.vstring(
'MSEL = 1 ! user defined subprocess',
15 'MSTJ(26) = 0 ! Mixing off',
16 'MDME(863,1) = 3 ! 0.020000 nu_e e+ D- ',
17 'MDME(864,1) = 3 ! 0.055000 nu_e e+ D*- ',
18 'MDME(865,1) = 3 ! 0.005000 nu_e e+ D_1- ',
19 'MDME(866,1) = 3 ! 0.005000 nu_e e+ D*_0-',
20 'MDME(867,1) = 3 ! 0.008000 nu_e e+ D*_1-',
21 'MDME(868,1) = 3 ! 0.012000 nu_e e+ D*_2-',
22 'MDME(869,1) = 3 ! 0.020000 nu_mu mu+ D- ',
23 'MDME(870,1) = 3 ! 0.055000 nu_mu mu+ D*- ',
24 'MDME(871,1) = 3 ! 0.005000 nu_mu mu+ D_1- ',
25 'MDME(872,1) = 3 ! 0.005000 nu_mu mu+ D*_0-',
26 'MDME(873,1) = 3 ! 0.008000 nu_mu mu+ D*_1-',
27 'MDME(874,1) = 3 ! 0.012000 nu_mu mu+ D*_2-',
28 'MDME(875,1) = 3 ! 0.010000 nu_tau tau+ D- ',
29 'MDME(876,1) = 3 ! 0.030000 nu_tau tau+ D*- ',
30 'MDME(877,1) = 3 ! 0.003500 D- pi+ ',
31 'MDME(878,1) = 3 ! 0.011000 D- rho+ ',
32 'MDME(879,1) = 3 ! 0.005500 D- a_1+ ',
33 'MDME(880,1) = 3 ! 0.004200 D*- pi+ ',
34 'MDME(881,1) = 3 ! 0.009000 D*- rho+ ',
35 'MDME(882,1) = 3 ! 0.018000 D*- a_1+ ',
36 'MDME(883,1) = 3 ! 0.015000 D- D_s+ ',
37 'MDME(884,1) = 3 ! 0.018500 D- D*_s+ ',
38 'MDME(885,1) = 3 ! 0.013500 D*- D_s+ ',
39 'MDME(886,1) = 3 ! 0.025000 D*- D*_s+ ',
40 'MDME(887,1) = 3 ! 0.000400 eta_c K0 ',
41 'MDME(888,1) = 3 ! 0.000700 eta_c K*0 ',
42 'MDME(889,1) = 2 ! 0.000800 J/psi K0 ',
43 'MDME(890,1) = 2 ! 0.001400 J/psi K*0 ',
44 'MDME(891,1) = 2 ! 0.001900 chi_1c K0 ',
45 'MDME(892,1) = 2 ! 0.002500 chi_1c K*0 ',
46 'MDME(893,1) = 3 ! 0.429100 u dbar cbar d',
47 'MDME(894,1) = 3 ! 0.080000 u cbar dbar d',
48 'MDME(895,1) = 2 ! 0.070000 c sbar cbar d',
49 'MDME(896,1) = 2 ! 0.020000 c cbar sbar d',
50 'MDME(897,1) = 3 ! 0.015000 u dbar ubar d',
51 'MDME(898,1) = 3 ! 0.005000 c sbar ubar d',
52 'MDME(908,1) = 3 ! 0.020000 nu_e e+ Dbar0 ',
53 'MDME(909,1) = 3 ! 0.055000 nu_e e+ D*bar0 ',
54 'MDME(910,1) = 3 ! 0.005000 nu_e e+ D_1bar0 ',
55 'MDME(911,1) = 3 ! 0.005000 nu_e e+ D*_0bar0',
56 'MDME(912,1) = 3 ! 0.008000 nu_e e+ D*_1bar0',
57 'MDME(913,1) = 3 ! 0.012000 nu_e e+ D*_2bar0',
58 'MDME(914,1) = 3 ! 0.020000 nu_mu mu+ Dbar0 ',
59 'MDME(915,1) = 3 ! 0.055000 nu_mu mu+ D*bar0 ',
60 'MDME(916,1) = 3 ! 0.005000 nu_mu mu+ D_1bar0 ',
61 'MDME(917,1) = 3 ! 0.005000 nu_mu mu+ D*_0bar0',
62 'MDME(918,1) = 3 ! 0.008000 nu_mu mu+ D*_1bar0',
63 'MDME(919,1) = 3 ! 0.012000 nu_mu mu+ D*_2bar0',
64 'MDME(920,1) = 3 ! 0.010000 nu_tau tau+ Dbar0 ',
65 'MDME(921,1) = 3 ! 0.030000 nu_tau tau+ D*bar0 ',
66 'MDME(922,1) = 3 ! 0.003500 Dbar0 pi+ ',
67 'MDME(923,1) = 3 ! 0.011000 Dbar0 rho+ ',
68 'MDME(924,1) = 3 ! 0.005500 Dbar0 a_1+ ',
69 'MDME(925,1) = 3 ! 0.004200 D*bar0 pi+ ',
70 'MDME(926,1) = 3 ! 0.009000 D*bar0 rho+ ',
71 'MDME(927,1) = 3 ! 0.018000 D*bar0 a_1+ ',
72 'MDME(928,1) = 3 ! 0.015000 Dbar0 D_s+ ',
73 'MDME(929,1) = 3 ! 0.018500 Dbar0 D*_s+ ',
74 'MDME(930,1) = 3 ! 0.013500 D*bar0 D_s+ ',
75 'MDME(931,1) = 3 ! 0.025000 D*bar0 D*_s+ ',
76 'MDME(932,1) = 3 ! 0.000400 eta_c K+ ',
77 'MDME(933,1) = 3 ! 0.000700 eta_c K*+ ',
78 'MDME(934,1) = 2 ! 0.000800 J/psi K+ ',
79 'MDME(935,1) = 2 ! 0.001400 J/psi K*+ ',
80 'MDME(936,1) = 2 ! 0.001900 chi_1c K+ ',
81 'MDME(937,1) = 2 ! 0.002500 chi_1c K*+ ',
82 'MDME(938,1) = 3 ! 0.429100 u dbar cbar u',
83 'MDME(939,1) = 3 ! 0.080000 u cbar dbar u',
84 'MDME(940,1) = 2 ! 0.070000 c sbar cbar u',
85 'MDME(941,1) = 2 ! 0.020000 c cbar sbar u',
86 'MDME(942,1) = 3 ! 0.015000 u dbar ubar u',
87 'MDME(943,1) = 3 ! 0.005000 c sbar ubar u',
88 'MDME(953,1) = 3 ! 0.020000 nu_e e+ D_s- ',
89 'MDME(954,1) = 3 ! 0.055000 nu_e e+ D*_s- ',
90 'MDME(955,1) = 3 ! 0.005000 nu_e e+ D_1s- ',
91 'MDME(956,1) = 3 ! 0.005000 nu_e e+ D*_0s-',
92 'MDME(957,1) = 3 ! 0.008000 nu_e e+ D*_1s-',
93 'MDME(958,1) = 3 ! 0.012000 nu_e e+ D*_2s-',
94 'MDME(959,1) = 3 ! 0.020000 nu_mu mu+ D_s- ',
95 'MDME(960,1) = 3 ! 0.055000 nu_mu mu+ D*_s- ',
96 'MDME(961,1) = 3 ! 0.005000 nu_mu mu+ D_1s- ',
97 'MDME(962,1) = 3 ! 0.005000 nu_mu mu+ D*_0s-',
98 'MDME(963,1) = 3 ! 0.008000 nu_mu mu+ D*_1s-',
99 'MDME(964,1) = 3 ! 0.012000 nu_mu mu+ D*_2s-',
100 'MDME(965,1) = 3 ! 0.010000 nu_tau tau+ D_s- ',
101 'MDME(966,1) = 3 ! 0.030000 nu_tau tau+ D*_s- ',
102 'MDME(967,1) = 3 ! 0.003500 D_s- pi+ ',
103 'MDME(968,1) = 3 ! 0.011000 D_s- rho+ ',
104 'MDME(969,1) = 3 ! 0.005500 D_s- a_1+ ',
105 'MDME(970,1) = 3 ! 0.004200 D*_s- pi+ ',
106 'MDME(971,1) = 3 ! 0.009000 D*_s- rho+ ',
107 'MDME(972,1) = 3 ! 0.018000 D*_s- a_1+ ',
108 'MDME(973,1) = 3 ! 0.015000 D_s- D_s+ ',
109 'MDME(974,1) = 3 ! 0.018500 D_s- D*_s+ ',
110 'MDME(975,1) = 3 ! 0.013500 D*_s- D_s+ ',
111 'MDME(976,1) = 3 ! 0.025000 D*_s- D*_s+ ',
112 'MDME(977,1) = 3 ! 0.000200 eta_c eta ',
113 'MDME(978,1) = 3 ! 0.000200 eta_c eta_ ',
114 'MDME(979,1) = 3 ! 0.000700 eta_c phi ',
115 'MDME(980,1) = 2 ! 0.000400 J/psi eta ',
116 'MDME(981,1) = 2 ! 0.000400 J/psi eta_ ',
117 'MDME(982,1) = 2 ! 0.001400 J/psi phi ',
118 'MDME(983,1) = 2 ! 0.001000 chi_1c eta ',
119 'MDME(984,1) = 2 ! 0.000900 chi_1c eta_ ',
120 'MDME(985,1) = 2 ! 0.002500 chi_1c phi ',
121 'MDME(986,1) = 3 ! 0.429100 u dbar cbar s',
122 'MDME(987,1) = 3 ! 0.080000 u cbar dba s',
123 'MDME(988,1) = 2 ! 0.070000 c sbar cbar s',
124 'MDME(989,1) = 2 ! 0.020000 c cbar sbar s',
125 'MDME(990,1) = 3 ! 0.015000 u dbar ubar s',
126 'MDME(991,1) = 3 ! 0.005000 c sbar ubar s',
127 'MDME(997,1) = 3 ! 0.047000 nu_tau tau+ ',
128 'MDME(998,1) = 3 ! 0.122000 c sbar ',
129 'MDME(999,1) = 3 ! 0.006000 c dbar ',
130 'MDME(1000,1) = 3 ! 0.012000 nu_e e+ eta_c',
131 'MDME(1001,1) = 2 ! 0.035000 nu_e e+ J/psi',
132 'MDME(1002,1) = 3 ! 0.012000 nu_mu mu+ eta_c',
133 'MDME(1003,1) = 2 ! 0.035000 nu_mu mu+ J/psi',
134 'MDME(1004,1) = 3 ! 0.003000 nu_tau tau+ eta_c',
135 'MDME(1005,1) = 2 ! 0.007000 nu_tau tau+ J/psi',
136 'MDME(1006,1) = 3 ! 0.150000 u dbar cbar c',
137 'MDME(1007,1) = 2 ! 0.037000 u cbar dbar c',
138 'MDME(1008,1) = 2 ! 0.008000 u sbar cbar c',
139 'MDME(1009,1) = 2 ! 0.002000 u cbar sbar c',
140 'MDME(1010,1) = 2 ! 0.050000 c sbar cbar c',
141 'MDME(1011,1) = 2 ! 0.015000 c cbar sbar c',
142 'MDME(1012,1) = 2 ! 0.003000 c dbar cbar c',
143 'MDME(1013,1) = 2 ! 0.001000 c cbar dbar c',
144 'MDME(1014,1) = 1 ! 0.014000 e+ nu_e B_s0 ',
145 'MDME(1015,1) = 1 ! 0.042000 e+ nu_e B*_s0',
146 'MDME(1016,1) = 1 ! 0.014000 mu+ nu_mu B_s0 ',
147 'MDME(1017,1) = 1 ! 0.042000 mu+ nu_mu B*_s0',
148 'MDME(1018,1) = 3 ! 0.240000 dbar u s bbar',
149 'MDME(1019,1) = 3 ! 0.065000 dbar s u bbar',
150 'MDME(1020,1) = 3 ! 0.012000 sbar u s bbar',
151 'MDME(1021,1) = 3 ! 0.003000 sbar s u bbar',
152 'MDME(1022,1) = 1 ! 0.001000 e+ nu_e B0 ',
153 'MDME(1023,1) = 1 ! 0.002000 e+ nu_e B*0 ',
154 'MDME(1024,1) = 1 ! 0.001000 mu+ nu_mu B0 ',
155 'MDME(1025,1) = 1 ! 0.002000 mu+ nu_mu B*0 ',
156 'MDME(1026,1) = 3 ! 0.014000 dbar u d bbar',
157 'MDME(1027,1) = 3 ! 0.003000 dbar d u bbar',
158 'MDME(1219,1) = 3 ! 0.105000 nu_ebar e- Lambda_c+',
159 'MDME(1220,1) = 3 ! 0.105000 nu_mubar mu- Lambda_c+',
160 'MDME(1221,1) = 3 ! 0.040000 nu_taubar tau- Lambda_c+',
161 'MDME(1222,1) = 3 ! 0.007700 Lambda_c+ pi- ',
162 'MDME(1223,1) = 3 ! 0.020000 Lambda_c+ rho- ',
163 'MDME(1224,1) = 3 ! 0.023500 Lambda_c+ a_1- ',
164 'MDME(1225,1) = 3 ! 0.028500 Lambda_c+ D_s- ',
165 'MDME(1226,1) = 3 ! 0.043500 Lambda_c+ D*_s- ',
166 'MDME(1227,1) = 3 ! 0.001100 eta_c Lambda0 ',
167 'MDME(1228,1) = 2 ! 0.002200 J/psi Lambda0 ',
168 'MDME(1229,1) = 2 ! 0.004400 chi_1c Lambda0 ',
169 'MDME(1230,1) = 3 ! 0.429100 ubar d c ud_0 ',
170 'MDME(1231,1) = 3 ! 0.080000 ubar c d ud_0 ',
171 'MDME(1232,1) = 2 ! 0.070000 cbar s c ud_0 ',
172 'MDME(1233,1) = 2 ! 0.020000 cbar c s ud_0 ',
173 'MDME(1234,1) = 3 ! 0.015000 ubar d u ud_0 ',
174 'MDME(1235,1) = 3 ! 0.005000 cbar s u ud_0 ',
175 'BRAT(891) = 0.000519 ! B0->Chi_1c 0.001900*0.273000',
176 'BRAT(892) = 0.000683 ! B0->Chi_1c 0.002500*0.273000',
177 'BRAT(936) = 0.000519 ! B+->Chi_1c 0.001900*0.273000',
178 'BRAT(936) = 0.000683 ! B+->Chi_1c 0.002500*0.273000',
179 'BRAT(983) = 0.000273 ! Bs->Chi_1c 0.001000*0.273000',
180 'BRAT(984) = 0.000246 ! Bs->Chi_1c 0.000900*0.273000',
181 'BRAT(985) = 0.000683 ! Bs->Chi_1c 0.002500*0.273000',
182 'BRAT(1229)= 0.001201 ! Lambda_b0->Chi_1c+lambda0 0.004400*0.273000',
183 'MDME(1555,1) = 1 ! 0.273000 J/psi gamma ',
184 'MDME(1556,1) = 0 ! 0.727000 rndmflav rndmflavbar ',
185 'MDME(858,1) = 0 ! 0.060200 e- e+',
186 'MDME(859,1) = 1 ! 0.060100 mu- mu+',
187 'MDME(860,1) = 0 ! 0.879700 rndmflav rndmflavbar'),
189 parameterSets = cms.vstring(
'pythiaUESettings',
194 bfilter = cms.EDFilter(
"PythiaFilter",
195 ParticleID = cms.untracked.int32(5)
198 jpsifilter = cms.EDFilter(
"PythiaFilter",
199 Status = cms.untracked.int32(2),
200 MaxEta = cms.untracked.double(20.0),
201 MinEta = cms.untracked.double(-20.0),
202 MinPt = cms.untracked.double(0.0),
203 ParticleID = cms.untracked.int32(443)
206 mumugenfilter = cms.EDFilter(
"MCParticlePairFilter",
207 Status = cms.untracked.vint32(1, 1),
208 MinPt = cms.untracked.vdouble(0.5, 0.5),
209 MinP = cms.untracked.vdouble(2.7, 2.7),
210 MaxEta = cms.untracked.vdouble(2.5, 2.5),
211 MinEta = cms.untracked.vdouble(-2.5, -2.5),
212 ParticleCharge = cms.untracked.int32(-1),
213 MaxInvMass = cms.untracked.double(4.0),
214 MinInvMass = cms.untracked.double(2.0),
215 ParticleID1 = cms.untracked.vint32(13),
216 ParticleID2 = cms.untracked.vint32(13)
219 ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*mumugenfilter)