CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Configuration/Generator/python/bJpsiX_cfi.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 from Configuration.Generator.PythiaUESettings_cfi import *
00004 Generator = cms.EDFilter("Pythia6GeneratorFilter",
00005     #untracked int32 maxEvents = 30000
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         # This is a vector of ParameterSet names to be read, in this order
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)