00001 import FWCore.ParameterSet.Config as cms
00002
00003 from Configuration.Generator.PythiaUESettings_cfi import *
00004 source = cms.Source("PythiaSource",
00005 pythiaPylistVerbosity = cms.untracked.int32(0),
00006 filterEfficiency = cms.untracked.double(0.141),
00007 pythiaHepMCVerbosity = cms.untracked.bool(False),
00008 crossSection = cms.untracked.double(9090000.0),
00009 comEnergy = cms.untracked.double(10000.0),
00010 maxEventsToPrint = cms.untracked.int32(0),
00011 PythiaParameters = cms.PSet(
00012 pythiaUESettingsBlock,
00013 processParameters = cms.vstring(
00014 'MSEL=62 ! Quarkonia NRQCD bb',
00015 'CKIN(3)=20. ! Min pthard',
00016 'CKIN(4)=-1. ! Max pthard',
00017 'MDME(1034,1)=0 ! 0.025200 e- e+',
00018 'MDME(1035,1)=1 ! 0.024800 mu- mu+',
00019 'MDME(1036,1)=0 ! 0.026700 tau- tau+',
00020 'MDME(1037,1)=0 ! 0.015000 d dbar',
00021 'MDME(1038,1)=0 ! 0.045000 u ubar',
00022 'MDME(1039,1)=0 ! 0.015000 s sbar',
00023 'MDME(1040,1)=0 ! 0.045000 c cbar',
00024 'MDME(1041,1)=0 ! 0.774300 g g g',
00025 'MDME(1042,1)=0 ! 0.029000 gamma g',
00026 'MSTP(142)=2 ! turns on the PYEVWT Pt re-weighting routine',
00027 'PARJ(13)=0.750 ! probability that a c or b meson has S=1',
00028 'PARJ(14)=0.162 ! probability that a meson with S=0 is produced with L=1, J=1',
00029 'PARJ(15)=0.018 ! probability that a meson with S=1 is produced with L=1, J=0',
00030 'PARJ(16)=0.054 ! probability that a meson with S=1 is produced with L=1, J=1',
00031 'MSTP(145)=0 ! choice of polarization',
00032 'MSTP(146)=0 ! choice of polarization frame ONLY when mstp(145)=1',
00033 'MSTP(147)=0 ! particular helicity or density matrix component when mstp(145)=1',
00034 'MSTP(148)=1 ! possibility to allow for final-state shower evolution, extreme case!',
00035 'MSTP(149)=1 ! if mstp(148)=1, it determines the kinematics of the QQ~3S1(8)->QQ~3S1(8)+g branching',
00036 'PARP(141)=1.16 ! New values for COM matrix elements',
00037 'PARP(142)=0.0119 ! New values for COM matrix elements',
00038 'PARP(143)=0.01 ! New values for COM matrix elements',
00039 'PARP(144)=0.01 ! New values for COM matrix elements',
00040 'PARP(145)=0.05 ! New values for COM matrix elements',
00041 'PARP(146)=9.28 ! New values for COM matrix elements',
00042 'PARP(147)=0.15 ! New values for COM matrix elements',
00043 'PARP(148)=0.02 ! New values for COM matrix elements',
00044 'PARP(149)=0.02 ! New values for COM matrix elements',
00045 'PARP(150)=0.09 ! New values for COM matrix elements'
00046 ),
00047
00048 parameterSets = cms.vstring('pythiaUESettings',
00049 'processParameters',
00050 'CSAParameters'),
00051 CSAParameters = cms.vstring('CSAMODE = 6 ! cross-section reweighted quarkonia')
00052 )
00053 )
00054
00055 oniafilter = cms.EDFilter("PythiaFilter",
00056 Status = cms.untracked.int32(2),
00057 MaxEta = cms.untracked.double(1000.0),
00058 MinEta = cms.untracked.double(-1000.0),
00059 MinPt = cms.untracked.double(0.0),
00060 ParticleID = cms.untracked.int32(553)
00061 )
00062
00063 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
00064 Status = cms.untracked.vint32(1, 1),
00065 MinPt = cms.untracked.vdouble(2.5, 2.5),
00066 MaxEta = cms.untracked.vdouble(2.5, 2.5),
00067 MinEta = cms.untracked.vdouble(-2.5, -2.5),
00068 ParticleCharge = cms.untracked.int32(-1),
00069 ParticleID1 = cms.untracked.vint32(13),
00070 ParticleID2 = cms.untracked.vint32(13)
00071 )
00072
00073 ProductionFilterSequence = cms.Sequence(oniafilter*mumugenfilter)
1.5.4