00001 import FWCore.ParameterSet.Config as cms
00002
00003 from Configuration.Generator.PythiaUESettings_cfi import *
00004 Generator = cms.EDFilter("Pythia6GeneratorFilter",
00005
00006 pythiaPylistVerbosity = cms.untracked.int32(0),
00007 filterEfficiency = cms.untracked.double(0.00013),
00008 pythiaHepMCVerbosity = cms.untracked.bool(False),
00009 crossSection = cms.untracked.double(54710000000.0),
00010 maxEventsToPrint = cms.untracked.int32(0),
00011 comEnergy = cms.double(10000.0),
00012 PythiaParameters = cms.PSet(
00013 pythiaUESettingsBlock,
00014 processParameters = cms.vstring('MSEL = 1 ! user defined subprocess',
00015 'MSTJ(26) = 0 ! Mixing off',
00016 'MDME(863,1) = 3 ! 0.020000 nu_e e+ D- ',
00017 'MDME(864,1) = 3 ! 0.055000 nu_e e+ D*- ',
00018 'MDME(865,1) = 3 ! 0.005000 nu_e e+ D_1- ',
00019 'MDME(866,1) = 3 ! 0.005000 nu_e e+ D*_0-',
00020 'MDME(867,1) = 3 ! 0.008000 nu_e e+ D*_1-',
00021 'MDME(868,1) = 3 ! 0.012000 nu_e e+ D*_2-',
00022 'MDME(869,1) = 3 ! 0.020000 nu_mu mu+ D- ',
00023 'MDME(870,1) = 3 ! 0.055000 nu_mu mu+ D*- ',
00024 'MDME(871,1) = 3 ! 0.005000 nu_mu mu+ D_1- ',
00025 'MDME(872,1) = 3 ! 0.005000 nu_mu mu+ D*_0-',
00026 'MDME(873,1) = 3 ! 0.008000 nu_mu mu+ D*_1-',
00027 'MDME(874,1) = 3 ! 0.012000 nu_mu mu+ D*_2-',
00028 'MDME(875,1) = 3 ! 0.010000 nu_tau tau+ D- ',
00029 'MDME(876,1) = 3 ! 0.030000 nu_tau tau+ D*- ',
00030 'MDME(877,1) = 3 ! 0.003500 D- pi+ ',
00031 'MDME(878,1) = 3 ! 0.011000 D- rho+ ',
00032 'MDME(879,1) = 3 ! 0.005500 D- a_1+ ',
00033 'MDME(880,1) = 3 ! 0.004200 D*- pi+ ',
00034 'MDME(881,1) = 3 ! 0.009000 D*- rho+ ',
00035 'MDME(882,1) = 3 ! 0.018000 D*- a_1+ ',
00036 'MDME(883,1) = 3 ! 0.015000 D- D_s+ ',
00037 'MDME(884,1) = 3 ! 0.018500 D- D*_s+ ',
00038 'MDME(885,1) = 3 ! 0.013500 D*- D_s+ ',
00039 'MDME(886,1) = 3 ! 0.025000 D*- D*_s+ ',
00040 'MDME(887,1) = 3 ! 0.000400 eta_c K0 ',
00041 'MDME(888,1) = 3 ! 0.000700 eta_c K*0 ',
00042 'MDME(889,1) = 2 ! 0.000800 J/psi K0 ',
00043 'MDME(890,1) = 2 ! 0.001400 J/psi K*0 ',
00044 'MDME(891,1) = 2 ! 0.001900 chi_1c K0 ',
00045 'MDME(892,1) = 2 ! 0.002500 chi_1c K*0 ',
00046 'MDME(893,1) = 3 ! 0.429100 u dbar cbar d',
00047 'MDME(894,1) = 3 ! 0.080000 u cbar dbar d',
00048 'MDME(895,1) = 2 ! 0.070000 c sbar cbar d',
00049 'MDME(896,1) = 2 ! 0.020000 c cbar sbar d',
00050 'MDME(897,1) = 3 ! 0.015000 u dbar ubar d',
00051 'MDME(898,1) = 3 ! 0.005000 c sbar ubar d',
00052 'MDME(908,1) = 3 ! 0.020000 nu_e e+ Dbar0 ',
00053 'MDME(909,1) = 3 ! 0.055000 nu_e e+ D*bar0 ',
00054 'MDME(910,1) = 3 ! 0.005000 nu_e e+ D_1bar0 ',
00055 'MDME(911,1) = 3 ! 0.005000 nu_e e+ D*_0bar0',
00056 'MDME(912,1) = 3 ! 0.008000 nu_e e+ D*_1bar0',
00057 'MDME(913,1) = 3 ! 0.012000 nu_e e+ D*_2bar0',
00058 'MDME(914,1) = 3 ! 0.020000 nu_mu mu+ Dbar0 ',
00059 'MDME(915,1) = 3 ! 0.055000 nu_mu mu+ D*bar0 ',
00060 'MDME(916,1) = 3 ! 0.005000 nu_mu mu+ D_1bar0 ',
00061 'MDME(917,1) = 3 ! 0.005000 nu_mu mu+ D*_0bar0',
00062 'MDME(918,1) = 3 ! 0.008000 nu_mu mu+ D*_1bar0',
00063 'MDME(919,1) = 3 ! 0.012000 nu_mu mu+ D*_2bar0',
00064 'MDME(920,1) = 3 ! 0.010000 nu_tau tau+ Dbar0 ',
00065 'MDME(921,1) = 3 ! 0.030000 nu_tau tau+ D*bar0 ',
00066 'MDME(922,1) = 3 ! 0.003500 Dbar0 pi+ ',
00067 'MDME(923,1) = 3 ! 0.011000 Dbar0 rho+ ',
00068 'MDME(924,1) = 3 ! 0.005500 Dbar0 a_1+ ',
00069 'MDME(925,1) = 3 ! 0.004200 D*bar0 pi+ ',
00070 'MDME(926,1) = 3 ! 0.009000 D*bar0 rho+ ',
00071 'MDME(927,1) = 3 ! 0.018000 D*bar0 a_1+ ',
00072 'MDME(928,1) = 3 ! 0.015000 Dbar0 D_s+ ',
00073 'MDME(929,1) = 3 ! 0.018500 Dbar0 D*_s+ ',
00074 'MDME(930,1) = 3 ! 0.013500 D*bar0 D_s+ ',
00075 'MDME(931,1) = 3 ! 0.025000 D*bar0 D*_s+ ',
00076 'MDME(932,1) = 3 ! 0.000400 eta_c K+ ',
00077 'MDME(933,1) = 3 ! 0.000700 eta_c K*+ ',
00078 'MDME(934,1) = 2 ! 0.000800 J/psi K+ ',
00079 'MDME(935,1) = 2 ! 0.001400 J/psi K*+ ',
00080 'MDME(936,1) = 2 ! 0.001900 chi_1c K+ ',
00081 'MDME(937,1) = 2 ! 0.002500 chi_1c K*+ ',
00082 'MDME(938,1) = 3 ! 0.429100 u dbar cbar u',
00083 'MDME(939,1) = 3 ! 0.080000 u cbar dbar u',
00084 'MDME(940,1) = 2 ! 0.070000 c sbar cbar u',
00085 'MDME(941,1) = 2 ! 0.020000 c cbar sbar u',
00086 'MDME(942,1) = 3 ! 0.015000 u dbar ubar u',
00087 'MDME(943,1) = 3 ! 0.005000 c sbar ubar u',
00088 'MDME(953,1) = 3 ! 0.020000 nu_e e+ D_s- ',
00089 'MDME(954,1) = 3 ! 0.055000 nu_e e+ D*_s- ',
00090 'MDME(955,1) = 3 ! 0.005000 nu_e e+ D_1s- ',
00091 'MDME(956,1) = 3 ! 0.005000 nu_e e+ D*_0s-',
00092 'MDME(957,1) = 3 ! 0.008000 nu_e e+ D*_1s-',
00093 'MDME(958,1) = 3 ! 0.012000 nu_e e+ D*_2s-',
00094 'MDME(959,1) = 3 ! 0.020000 nu_mu mu+ D_s- ',
00095 'MDME(960,1) = 3 ! 0.055000 nu_mu mu+ D*_s- ',
00096 'MDME(961,1) = 3 ! 0.005000 nu_mu mu+ D_1s- ',
00097 'MDME(962,1) = 3 ! 0.005000 nu_mu mu+ D*_0s-',
00098 'MDME(963,1) = 3 ! 0.008000 nu_mu mu+ D*_1s-',
00099 'MDME(964,1) = 3 ! 0.012000 nu_mu mu+ D*_2s-',
00100 'MDME(965,1) = 3 ! 0.010000 nu_tau tau+ D_s- ',
00101 'MDME(966,1) = 3 ! 0.030000 nu_tau tau+ D*_s- ',
00102 'MDME(967,1) = 3 ! 0.003500 D_s- pi+ ',
00103 'MDME(968,1) = 3 ! 0.011000 D_s- rho+ ',
00104 'MDME(969,1) = 3 ! 0.005500 D_s- a_1+ ',
00105 'MDME(970,1) = 3 ! 0.004200 D*_s- pi+ ',
00106 'MDME(971,1) = 3 ! 0.009000 D*_s- rho+ ',
00107 'MDME(972,1) = 3 ! 0.018000 D*_s- a_1+ ',
00108 'MDME(973,1) = 3 ! 0.015000 D_s- D_s+ ',
00109 'MDME(974,1) = 3 ! 0.018500 D_s- D*_s+ ',
00110 'MDME(975,1) = 3 ! 0.013500 D*_s- D_s+ ',
00111 'MDME(976,1) = 3 ! 0.025000 D*_s- D*_s+ ',
00112 'MDME(977,1) = 3 ! 0.000200 eta_c eta ',
00113 'MDME(978,1) = 3 ! 0.000200 eta_c eta_ ',
00114 'MDME(979,1) = 3 ! 0.000700 eta_c phi ',
00115 'MDME(980,1) = 2 ! 0.000400 J/psi eta ',
00116 'MDME(981,1) = 2 ! 0.000400 J/psi eta_ ',
00117 'MDME(982,1) = 2 ! 0.001400 J/psi phi ',
00118 'MDME(983,1) = 2 ! 0.001000 chi_1c eta ',
00119 'MDME(984,1) = 2 ! 0.000900 chi_1c eta_ ',
00120 'MDME(985,1) = 2 ! 0.002500 chi_1c phi ',
00121 'MDME(986,1) = 3 ! 0.429100 u dbar cbar s',
00122 'MDME(987,1) = 3 ! 0.080000 u cbar dba s',
00123 'MDME(988,1) = 2 ! 0.070000 c sbar cbar s',
00124 'MDME(989,1) = 2 ! 0.020000 c cbar sbar s',
00125 'MDME(990,1) = 3 ! 0.015000 u dbar ubar s',
00126 'MDME(991,1) = 3 ! 0.005000 c sbar ubar s',
00127 'MDME(997,1) = 3 ! 0.047000 nu_tau tau+ ',
00128 'MDME(998,1) = 3 ! 0.122000 c sbar ',
00129 'MDME(999,1) = 3 ! 0.006000 c dbar ',
00130 'MDME(1000,1) = 3 ! 0.012000 nu_e e+ eta_c',
00131 'MDME(1001,1) = 2 ! 0.035000 nu_e e+ J/psi',
00132 'MDME(1002,1) = 3 ! 0.012000 nu_mu mu+ eta_c',
00133 'MDME(1003,1) = 2 ! 0.035000 nu_mu mu+ J/psi',
00134 'MDME(1004,1) = 3 ! 0.003000 nu_tau tau+ eta_c',
00135 'MDME(1005,1) = 2 ! 0.007000 nu_tau tau+ J/psi',
00136 'MDME(1006,1) = 3 ! 0.150000 u dbar cbar c',
00137 'MDME(1007,1) = 2 ! 0.037000 u cbar dbar c',
00138 'MDME(1008,1) = 2 ! 0.008000 u sbar cbar c',
00139 'MDME(1009,1) = 2 ! 0.002000 u cbar sbar c',
00140 'MDME(1010,1) = 2 ! 0.050000 c sbar cbar c',
00141 'MDME(1011,1) = 2 ! 0.015000 c cbar sbar c',
00142 'MDME(1012,1) = 2 ! 0.003000 c dbar cbar c',
00143 'MDME(1013,1) = 2 ! 0.001000 c cbar dbar c',
00144 'MDME(1014,1) = 1 ! 0.014000 e+ nu_e B_s0 ',
00145 'MDME(1015,1) = 1 ! 0.042000 e+ nu_e B*_s0',
00146 'MDME(1016,1) = 1 ! 0.014000 mu+ nu_mu B_s0 ',
00147 'MDME(1017,1) = 1 ! 0.042000 mu+ nu_mu B*_s0',
00148 'MDME(1018,1) = 3 ! 0.240000 dbar u s bbar',
00149 'MDME(1019,1) = 3 ! 0.065000 dbar s u bbar',
00150 'MDME(1020,1) = 3 ! 0.012000 sbar u s bbar',
00151 'MDME(1021,1) = 3 ! 0.003000 sbar s u bbar',
00152 'MDME(1022,1) = 1 ! 0.001000 e+ nu_e B0 ',
00153 'MDME(1023,1) = 1 ! 0.002000 e+ nu_e B*0 ',
00154 'MDME(1024,1) = 1 ! 0.001000 mu+ nu_mu B0 ',
00155 'MDME(1025,1) = 1 ! 0.002000 mu+ nu_mu B*0 ',
00156 'MDME(1026,1) = 3 ! 0.014000 dbar u d bbar',
00157 'MDME(1027,1) = 3 ! 0.003000 dbar d u bbar',
00158 'MDME(1219,1) = 3 ! 0.105000 nu_ebar e- Lambda_c+',
00159 'MDME(1220,1) = 3 ! 0.105000 nu_mubar mu- Lambda_c+',
00160 'MDME(1221,1) = 3 ! 0.040000 nu_taubar tau- Lambda_c+',
00161 'MDME(1222,1) = 3 ! 0.007700 Lambda_c+ pi- ',
00162 'MDME(1223,1) = 3 ! 0.020000 Lambda_c+ rho- ',
00163 'MDME(1224,1) = 3 ! 0.023500 Lambda_c+ a_1- ',
00164 'MDME(1225,1) = 3 ! 0.028500 Lambda_c+ D_s- ',
00165 'MDME(1226,1) = 3 ! 0.043500 Lambda_c+ D*_s- ',
00166 'MDME(1227,1) = 3 ! 0.001100 eta_c Lambda0 ',
00167 'MDME(1228,1) = 2 ! 0.002200 J/psi Lambda0 ',
00168 'MDME(1229,1) = 2 ! 0.004400 chi_1c Lambda0 ',
00169 'MDME(1230,1) = 3 ! 0.429100 ubar d c ud_0 ',
00170 'MDME(1231,1) = 3 ! 0.080000 ubar c d ud_0 ',
00171 'MDME(1232,1) = 2 ! 0.070000 cbar s c ud_0 ',
00172 'MDME(1233,1) = 2 ! 0.020000 cbar c s ud_0 ',
00173 'MDME(1234,1) = 3 ! 0.015000 ubar d u ud_0 ',
00174 'MDME(1235,1) = 3 ! 0.005000 cbar s u ud_0 ',
00175 'BRAT(891) = 0.000519 ! B0->Chi_1c 0.001900*0.273000',
00176 'BRAT(892) = 0.000683 ! B0->Chi_1c 0.002500*0.273000',
00177 'BRAT(936) = 0.000519 ! B+->Chi_1c 0.001900*0.273000',
00178 'BRAT(936) = 0.000683 ! B+->Chi_1c 0.002500*0.273000',
00179 'BRAT(983) = 0.000273 ! Bs->Chi_1c 0.001000*0.273000',
00180 'BRAT(984) = 0.000246 ! Bs->Chi_1c 0.000900*0.273000',
00181 'BRAT(985) = 0.000683 ! Bs->Chi_1c 0.002500*0.273000',
00182 'BRAT(1229)= 0.001201 ! Lambda_b0->Chi_1c+lambda0 0.004400*0.273000',
00183 'MDME(1555,1) = 1 ! 0.273000 J/psi gamma ',
00184 'MDME(1556,1) = 0 ! 0.727000 rndmflav rndmflavbar ',
00185 'MDME(858,1) = 0 ! 0.060200 e- e+',
00186 'MDME(859,1) = 1 ! 0.060100 mu- mu+',
00187 'MDME(860,1) = 0 ! 0.879700 rndmflav rndmflavbar'),
00188
00189 parameterSets = cms.vstring('pythiaUESettings',
00190 'processParameters')
00191 )
00192 )
00193
00194 bfilter = cms.EDFilter("PythiaFilter",
00195 ParticleID = cms.untracked.int32(5)
00196 )
00197
00198 jpsifilter = cms.EDFilter("PythiaFilter",
00199 Status = cms.untracked.int32(2),
00200 MaxEta = cms.untracked.double(20.0),
00201 MinEta = cms.untracked.double(-20.0),
00202 MinPt = cms.untracked.double(0.0),
00203 ParticleID = cms.untracked.int32(443)
00204 )
00205
00206 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
00207 Status = cms.untracked.vint32(1, 1),
00208 MinPt = cms.untracked.vdouble(2.5, 2.5),
00209 MaxEta = cms.untracked.vdouble(2.5, 2.5),
00210 MinEta = cms.untracked.vdouble(-2.5, -2.5),
00211 ParticleCharge = cms.untracked.int32(-1),
00212 MaxInvMass = cms.untracked.double(4.0),
00213 MinInvMass = cms.untracked.double(2.0),
00214 ParticleID1 = cms.untracked.vint32(13),
00215 ParticleID2 = cms.untracked.vint32(13)
00216 )
00217
00218 ProductionFilterSequence = cms.Sequence(Generator*bfilter*jpsifilter*mumugenfilter)