Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:03:33

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 from Configuration.Generator.PythiaUESettings_cfi import *
0004 generator = cms.EDFilter("Pythia6GeneratorFilter",
0005     #untracked int32 maxEvents = 30000
0006     pythiaPylistVerbosity = cms.untracked.int32(0),
0007     filterEfficiency = cms.untracked.double(0.00013),
0008     pythiaHepMCVerbosity = cms.untracked.bool(False),
0009     crossSection = cms.untracked.double(54710000000.0),
0010     maxEventsToPrint = cms.untracked.int32(0),
0011     comEnergy = cms.double(10000.0),
0012     PythiaParameters = cms.PSet(
0013         pythiaUESettingsBlock,
0014         processParameters = cms.vstring('MSEL        = 1        ! user defined subprocess', 
0015             'MSTJ(26)    = 0        ! Mixing off', 
0016             'MDME(863,1) = 3 ! 0.020000    nu_e     e+      D-   ', 
0017             'MDME(864,1) = 3 ! 0.055000    nu_e     e+      D*-  ', 
0018             'MDME(865,1) = 3 ! 0.005000    nu_e     e+      D_1- ', 
0019             'MDME(866,1) = 3 ! 0.005000    nu_e     e+      D*_0-', 
0020             'MDME(867,1) = 3 ! 0.008000    nu_e     e+      D*_1-', 
0021             'MDME(868,1) = 3 ! 0.012000    nu_e     e+      D*_2-', 
0022             'MDME(869,1) = 3 ! 0.020000    nu_mu    mu+     D-   ', 
0023             'MDME(870,1) = 3 ! 0.055000    nu_mu    mu+     D*-  ', 
0024             'MDME(871,1) = 3 ! 0.005000    nu_mu    mu+     D_1- ', 
0025             'MDME(872,1) = 3 ! 0.005000    nu_mu    mu+     D*_0-', 
0026             'MDME(873,1) = 3 ! 0.008000    nu_mu    mu+     D*_1-', 
0027             'MDME(874,1) = 3 ! 0.012000    nu_mu    mu+     D*_2-', 
0028             'MDME(875,1) = 3 ! 0.010000    nu_tau   tau+    D-   ', 
0029             'MDME(876,1) = 3 ! 0.030000    nu_tau   tau+    D*-  ', 
0030             'MDME(877,1) = 3 ! 0.003500    D-       pi+          ', 
0031             'MDME(878,1) = 3 ! 0.011000    D-       rho+         ', 
0032             'MDME(879,1) = 3 ! 0.005500    D-       a_1+         ', 
0033             'MDME(880,1) = 3 ! 0.004200    D*-      pi+          ', 
0034             'MDME(881,1) = 3 ! 0.009000    D*-      rho+         ', 
0035             'MDME(882,1) = 3 ! 0.018000    D*-      a_1+         ', 
0036             'MDME(883,1) = 3 ! 0.015000    D-       D_s+         ', 
0037             'MDME(884,1) = 3 ! 0.018500    D-       D*_s+        ', 
0038             'MDME(885,1) = 3 ! 0.013500    D*-      D_s+         ', 
0039             'MDME(886,1) = 3 ! 0.025000    D*-      D*_s+        ', 
0040             'MDME(887,1) = 3 ! 0.000400    eta_c    K0           ', 
0041             'MDME(888,1) = 3 ! 0.000700    eta_c    K*0          ', 
0042             'MDME(889,1) = 2 ! 0.000800    J/psi    K0           ', 
0043             'MDME(890,1) = 2 ! 0.001400    J/psi    K*0          ', 
0044             'MDME(891,1) = 2 ! 0.001900    chi_1c   K0           ', 
0045             'MDME(892,1) = 2 ! 0.002500    chi_1c   K*0          ', 
0046             'MDME(893,1) = 3 ! 0.429100    u     dbar    cbar   d', 
0047             'MDME(894,1) = 3 ! 0.080000    u     cbar    dbar   d', 
0048             'MDME(895,1) = 2 ! 0.070000    c     sbar    cbar   d', 
0049             'MDME(896,1) = 2 ! 0.020000    c     cbar    sbar   d', 
0050             'MDME(897,1) = 3 ! 0.015000    u     dbar    ubar   d', 
0051             'MDME(898,1) = 3 ! 0.005000    c     sbar    ubar   d', 
0052             'MDME(908,1) = 3 ! 0.020000    nu_e     e+    Dbar0   ', 
0053             'MDME(909,1) = 3 ! 0.055000    nu_e     e+    D*bar0  ', 
0054             'MDME(910,1) = 3 ! 0.005000    nu_e     e+    D_1bar0 ', 
0055             'MDME(911,1) = 3 ! 0.005000    nu_e     e+    D*_0bar0', 
0056             'MDME(912,1) = 3 ! 0.008000    nu_e     e+    D*_1bar0', 
0057             'MDME(913,1) = 3 ! 0.012000    nu_e     e+    D*_2bar0', 
0058             'MDME(914,1) = 3 ! 0.020000    nu_mu    mu+   Dbar0   ', 
0059             'MDME(915,1) = 3 ! 0.055000    nu_mu    mu+   D*bar0  ', 
0060             'MDME(916,1) = 3 ! 0.005000    nu_mu    mu+   D_1bar0 ', 
0061             'MDME(917,1) = 3 ! 0.005000    nu_mu    mu+   D*_0bar0', 
0062             'MDME(918,1) = 3 ! 0.008000    nu_mu    mu+   D*_1bar0', 
0063             'MDME(919,1) = 3 ! 0.012000    nu_mu    mu+   D*_2bar0', 
0064             'MDME(920,1) = 3 ! 0.010000    nu_tau   tau+  Dbar0   ', 
0065             'MDME(921,1) = 3 ! 0.030000    nu_tau   tau+  D*bar0  ', 
0066             'MDME(922,1) = 3 ! 0.003500    Dbar0      pi+         ', 
0067             'MDME(923,1) = 3 ! 0.011000    Dbar0      rho+        ', 
0068             'MDME(924,1) = 3 ! 0.005500    Dbar0      a_1+        ', 
0069             'MDME(925,1) = 3 ! 0.004200    D*bar0     pi+         ', 
0070             'MDME(926,1) = 3 ! 0.009000    D*bar0     rho+        ', 
0071             'MDME(927,1) = 3 ! 0.018000    D*bar0     a_1+        ', 
0072             'MDME(928,1) = 3 ! 0.015000    Dbar0      D_s+        ', 
0073             'MDME(929,1) = 3 ! 0.018500    Dbar0      D*_s+       ', 
0074             'MDME(930,1) = 3 ! 0.013500    D*bar0     D_s+        ', 
0075             'MDME(931,1) = 3 ! 0.025000    D*bar0     D*_s+       ', 
0076             'MDME(932,1) = 3 ! 0.000400    eta_c      K+          ', 
0077             'MDME(933,1) = 3 ! 0.000700    eta_c      K*+         ', 
0078             'MDME(934,1) = 2 ! 0.000800    J/psi      K+          ', 
0079             'MDME(935,1) = 2 ! 0.001400    J/psi      K*+         ', 
0080             'MDME(936,1) = 2 ! 0.001900    chi_1c     K+          ', 
0081             'MDME(937,1) = 2 ! 0.002500    chi_1c     K*+         ', 
0082             'MDME(938,1) = 3 ! 0.429100    u     dbar    cbar    u', 
0083             'MDME(939,1) = 3 ! 0.080000    u     cbar    dbar    u', 
0084             'MDME(940,1) = 2 ! 0.070000    c     sbar    cbar    u', 
0085             'MDME(941,1) = 2 ! 0.020000    c     cbar    sbar    u', 
0086             'MDME(942,1) = 3 ! 0.015000    u     dbar    ubar    u', 
0087             'MDME(943,1) = 3 ! 0.005000    c     sbar    ubar    u', 
0088             'MDME(953,1) = 3 ! 0.020000    nu_e     e+      D_s-  ', 
0089             'MDME(954,1) = 3 ! 0.055000    nu_e     e+      D*_s- ', 
0090             'MDME(955,1) = 3 ! 0.005000    nu_e     e+      D_1s- ', 
0091             'MDME(956,1) = 3 ! 0.005000    nu_e     e+      D*_0s-', 
0092             'MDME(957,1) = 3 ! 0.008000    nu_e     e+      D*_1s-', 
0093             'MDME(958,1) = 3 ! 0.012000    nu_e     e+      D*_2s-', 
0094             'MDME(959,1) = 3 ! 0.020000    nu_mu    mu+     D_s-  ', 
0095             'MDME(960,1) = 3 ! 0.055000    nu_mu    mu+     D*_s- ', 
0096             'MDME(961,1) = 3 ! 0.005000    nu_mu    mu+     D_1s- ', 
0097             'MDME(962,1) = 3 ! 0.005000    nu_mu    mu+     D*_0s-', 
0098             'MDME(963,1) = 3 ! 0.008000    nu_mu    mu+     D*_1s-', 
0099             'MDME(964,1) = 3 ! 0.012000    nu_mu    mu+     D*_2s-', 
0100             'MDME(965,1) = 3 ! 0.010000    nu_tau   tau+    D_s-  ', 
0101             'MDME(966,1) = 3 ! 0.030000    nu_tau   tau+    D*_s- ', 
0102             'MDME(967,1) = 3 ! 0.003500    D_s-        pi+        ', 
0103             'MDME(968,1) = 3 ! 0.011000    D_s-        rho+       ', 
0104             'MDME(969,1) = 3 ! 0.005500    D_s-        a_1+       ', 
0105             'MDME(970,1) = 3 ! 0.004200    D*_s-       pi+        ', 
0106             'MDME(971,1) = 3 ! 0.009000    D*_s-       rho+       ', 
0107             'MDME(972,1) = 3 ! 0.018000    D*_s-       a_1+       ', 
0108             'MDME(973,1) = 3 ! 0.015000    D_s-        D_s+       ', 
0109             'MDME(974,1) = 3 ! 0.018500    D_s-        D*_s+      ', 
0110             'MDME(975,1) = 3 ! 0.013500    D*_s-       D_s+       ', 
0111             'MDME(976,1) = 3 ! 0.025000    D*_s-       D*_s+      ', 
0112             'MDME(977,1) = 3 ! 0.000200    eta_c       eta        ', 
0113             'MDME(978,1) = 3 ! 0.000200    eta_c       eta_       ', 
0114             'MDME(979,1) = 3 ! 0.000700    eta_c       phi        ', 
0115             'MDME(980,1) = 2 ! 0.000400    J/psi       eta        ', 
0116             'MDME(981,1) = 2 ! 0.000400    J/psi       eta_       ', 
0117             'MDME(982,1) = 2 ! 0.001400    J/psi       phi        ', 
0118             'MDME(983,1) = 2 ! 0.001000    chi_1c      eta        ', 
0119             'MDME(984,1) = 2 ! 0.000900    chi_1c      eta_       ', 
0120             'MDME(985,1) = 2 ! 0.002500    chi_1c      phi        ', 
0121             'MDME(986,1) = 3 ! 0.429100    u     dbar    cbar    s', 
0122             'MDME(987,1) = 3 ! 0.080000    u     cbar    dba     s', 
0123             'MDME(988,1) = 2 ! 0.070000    c     sbar    cbar    s', 
0124             'MDME(989,1) = 2 ! 0.020000    c     cbar    sbar    s', 
0125             'MDME(990,1) = 3 ! 0.015000    u     dbar    ubar    s', 
0126             'MDME(991,1) = 3 ! 0.005000    c     sbar    ubar    s', 
0127             'MDME(997,1) = 3 ! 0.047000    nu_tau      tau+       ', 
0128             'MDME(998,1) = 3 ! 0.122000    c           sbar       ', 
0129             'MDME(999,1) = 3 ! 0.006000    c           dbar       ', 
0130             'MDME(1000,1) = 3 ! 0.012000    nu_e      e+     eta_c', 
0131             'MDME(1001,1) = 2 ! 0.035000    nu_e      e+     J/psi', 
0132             'MDME(1002,1) = 3 ! 0.012000    nu_mu     mu+    eta_c', 
0133             'MDME(1003,1) = 2 ! 0.035000    nu_mu     mu+    J/psi', 
0134             'MDME(1004,1) = 3 ! 0.003000    nu_tau    tau+   eta_c', 
0135             'MDME(1005,1) = 2 ! 0.007000    nu_tau    tau+   J/psi', 
0136             'MDME(1006,1) = 3 ! 0.150000    u    dbar   cbar     c', 
0137             'MDME(1007,1) = 2 ! 0.037000    u    cbar   dbar     c', 
0138             'MDME(1008,1) = 2 ! 0.008000    u    sbar   cbar     c', 
0139             'MDME(1009,1) = 2 ! 0.002000    u    cbar   sbar     c', 
0140             'MDME(1010,1) = 2 ! 0.050000    c    sbar   cbar     c', 
0141             'MDME(1011,1) = 2 ! 0.015000    c    cbar   sbar     c', 
0142             'MDME(1012,1) = 2 ! 0.003000    c    dbar   cbar     c', 
0143             'MDME(1013,1) = 2 ! 0.001000    c    cbar   dbar     c', 
0144             'MDME(1014,1) = 1 ! 0.014000    e+      nu_e     B_s0 ', 
0145             'MDME(1015,1) = 1 ! 0.042000    e+      nu_e     B*_s0', 
0146             'MDME(1016,1) = 1 ! 0.014000    mu+     nu_mu    B_s0 ', 
0147             'MDME(1017,1) = 1 ! 0.042000    mu+     nu_mu    B*_s0', 
0148             'MDME(1018,1) = 3 ! 0.240000    dbar    u    s    bbar', 
0149             'MDME(1019,1) = 3 ! 0.065000    dbar    s    u    bbar', 
0150             'MDME(1020,1) = 3 ! 0.012000    sbar    u    s    bbar', 
0151             'MDME(1021,1) = 3 ! 0.003000    sbar    s    u    bbar', 
0152             'MDME(1022,1) = 1 ! 0.001000    e+      nu_e      B0  ', 
0153             'MDME(1023,1) = 1 ! 0.002000    e+      nu_e      B*0 ', 
0154             'MDME(1024,1) = 1 ! 0.001000    mu+     nu_mu     B0  ', 
0155             'MDME(1025,1) = 1 ! 0.002000    mu+     nu_mu     B*0 ', 
0156             'MDME(1026,1) = 3 ! 0.014000    dbar    u    d    bbar', 
0157             'MDME(1027,1) = 3 ! 0.003000    dbar    d    u    bbar', 
0158             'MDME(1219,1) = 3 ! 0.105000    nu_ebar     e-     Lambda_c+', 
0159             'MDME(1220,1) = 3 ! 0.105000    nu_mubar    mu-    Lambda_c+', 
0160             'MDME(1221,1) = 3 ! 0.040000    nu_taubar   tau-   Lambda_c+', 
0161             'MDME(1222,1) = 3 ! 0.007700    Lambda_c+       pi-         ', 
0162             'MDME(1223,1) = 3 ! 0.020000    Lambda_c+       rho-        ', 
0163             'MDME(1224,1) = 3 ! 0.023500    Lambda_c+       a_1-        ', 
0164             'MDME(1225,1) = 3 ! 0.028500    Lambda_c+       D_s-        ', 
0165             'MDME(1226,1) = 3 ! 0.043500    Lambda_c+       D*_s-       ', 
0166             'MDME(1227,1) = 3 ! 0.001100    eta_c           Lambda0     ', 
0167             'MDME(1228,1) = 2 ! 0.002200    J/psi           Lambda0     ', 
0168             'MDME(1229,1) = 2  ! 0.004400    chi_1c          Lambda0     ', 
0169             'MDME(1230,1) = 3 ! 0.429100    ubar    d     c    ud_0     ', 
0170             'MDME(1231,1) = 3 ! 0.080000    ubar    c     d    ud_0     ', 
0171             'MDME(1232,1) = 2 ! 0.070000    cbar    s     c    ud_0     ', 
0172             'MDME(1233,1) = 2 ! 0.020000    cbar    c     s    ud_0     ', 
0173             'MDME(1234,1) = 3 ! 0.015000    ubar    d     u    ud_0     ', 
0174             'MDME(1235,1) = 3 ! 0.005000    cbar    s     u    ud_0     ', 
0175             'BRAT(891) = 0.000519 ! B0->Chi_1c 0.001900*0.273000', 
0176             'BRAT(892) = 0.000683 ! B0->Chi_1c 0.002500*0.273000', 
0177             'BRAT(936) = 0.000519 ! B+->Chi_1c 0.001900*0.273000', 
0178             'BRAT(936) = 0.000683 ! B+->Chi_1c 0.002500*0.273000', 
0179             'BRAT(983) = 0.000273 ! Bs->Chi_1c 0.001000*0.273000', 
0180             'BRAT(984) = 0.000246 ! Bs->Chi_1c 0.000900*0.273000', 
0181             'BRAT(985) = 0.000683 ! Bs->Chi_1c 0.002500*0.273000', 
0182             'BRAT(1229)= 0.001201 ! Lambda_b0->Chi_1c+lambda0 0.004400*0.273000', 
0183             'MDME(1555,1) = 1 ! 0.273000    J/psi           gamma ', 
0184             'MDME(1556,1) = 0 ! 0.727000    rndmflav        rndmflavbar ', 
0185             'MDME(858,1) = 0  ! 0.060200    e-    e+', 
0186             'MDME(859,1) = 1  ! 0.060100    mu-  mu+', 
0187             'MDME(860,1) = 0  ! 0.879700    rndmflav        rndmflavbar'),
0188         # This is a vector of ParameterSet names to be read, in this order
0189         parameterSets = cms.vstring('pythiaUESettings', 
0190             'processParameters')
0191     )
0192 )
0193 
0194 bfilter = cms.EDFilter("PythiaFilter",
0195     ParticleID = cms.untracked.int32(5)
0196 )
0197 
0198 jpsifilter = cms.EDFilter("PythiaFilter",
0199     Status = cms.untracked.int32(2),
0200     MaxEta = cms.untracked.double(20.0),
0201     MinEta = cms.untracked.double(-20.0),
0202     MinPt = cms.untracked.double(0.0),
0203     ParticleID = cms.untracked.int32(443)
0204 )
0205 
0206 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
0207     Status = cms.untracked.vint32(1, 1),
0208     MinPt = cms.untracked.vdouble(2.5, 2.5),
0209     MaxEta = cms.untracked.vdouble(2.5, 2.5),
0210     MinEta = cms.untracked.vdouble(-2.5, -2.5),
0211     ParticleCharge = cms.untracked.int32(-1),
0212     MaxInvMass = cms.untracked.double(4.0),
0213     MinInvMass = cms.untracked.double(2.0),
0214     ParticleID1 = cms.untracked.vint32(13),
0215     ParticleID2 = cms.untracked.vint32(13)
0216 )
0217 
0218 ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*mumugenfilter)