test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
bJpsiX_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
4 generator = cms.EDFilter("Pythia6GeneratorFilter",
5  #untracked int32 maxEvents = 30000
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(10000.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'),
188  # This is a vector of ParameterSet names to be read, in this order
189  parameterSets = cms.vstring('pythiaUESettings',
190  'processParameters')
191  )
192 )
193 
194 bfilter = cms.EDFilter("PythiaFilter",
195  ParticleID = cms.untracked.int32(5)
196 )
197 
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)
204 )
205 
206 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
207  Status = cms.untracked.vint32(1, 1),
208  MinPt = cms.untracked.vdouble(2.5, 2.5),
209  MaxEta = cms.untracked.vdouble(2.5, 2.5),
210  MinEta = cms.untracked.vdouble(-2.5, -2.5),
211  ParticleCharge = cms.untracked.int32(-1),
212  MaxInvMass = cms.untracked.double(4.0),
213  MinInvMass = cms.untracked.double(2.0),
214  ParticleID1 = cms.untracked.vint32(13),
215  ParticleID2 = cms.untracked.vint32(13)
216 )
217 
218 ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*mumugenfilter)