CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JpsiMM_14TeV_cfi.py
Go to the documentation of this file.
2 
4 generator = cms.EDFilter("Pythia6GeneratorFilter",
5  pythiaPylistVerbosity = cms.untracked.int32(0),
6  filterEfficiency = cms.untracked.double(0.138),
7  pythiaHepMCVerbosity = cms.untracked.bool(False),
8  crossSection = cms.untracked.double(1872000.0),
9  comEnergy = cms.double(14000.0),
10  maxEventsToPrint = cms.untracked.int32(0),
11  PythiaParameters = cms.PSet(
12  pythiaUESettingsBlock,
13  processParameters = cms.vstring(
14  'MSEL=61 ! Quarkonia',
15  'CKIN(3)=10. ! Min pthard',
16  'CKIN(4)=-1. ! Max pthard',
17  'MDME(858,1) = 0 ! 0.060200 e- e+',
18  'MDME(859,1) = 1 ! 0.060100 mu- mu+',
19  'MDME(860,1) = 0 ! 0.879700 rndmflav rndmflavbar',
20  'MSTP(142)=2 ! turns on the PYEVWT Pt re-weighting routine',
21  'PARJ(13)=0.750 ! probability that a c or b meson has S=1',
22  'PARJ(14)=0.162 ! probability that a meson with S=0 is produced with L=1, J=1',
23  'PARJ(15)=0.018 ! probability that a meson with S=1 is produced with L=1, J=0',
24  'PARJ(16)=0.054 ! probability that a meson with S=1 is produced with L=1, J=1',
25  'MSTP(145)=0 ! choice of polarization',
26  'MSTP(146)=0 ! choice of polarization frame ONLY when mstp(145)=1',
27  'MSTP(147)=0 ! particular helicity or density matrix component when mstp(145)=1',
28  'MSTP(148)=1 ! possibility to allow for final-state shower evolution, extreme case !',
29  'MSTP(149)=1 ! if mstp(148)=1, it determines the kinematics of the QQ~3S1(8)->QQ~3S1(8)+g branching',
30  'PARP(141)=1.16 ! New values for COM matrix elements',
31  'PARP(142)=0.0119 ! New values for COM matrix elements',
32  'PARP(143)=0.01 ! New values for COM matrix elements',
33  'PARP(144)=0.01 ! New values for COM matrix elements',
34  'PARP(145)=0.05 ! New values for COM matrix elements',
35  'PARP(146)=9.28 ! New values for COM matrix elements',
36  'PARP(147)=0.15 ! New values for COM matrix elements',
37  'PARP(148)=0.02 ! New values for COM matrix elements',
38  'PARP(149)=0.02 ! New values for COM matrix elements',
39  'PARP(150)=0.085 ! New values for COM matrix elements',
40  'BRAT(859)=1.000 ! J/psi->mu+mu-',
41  'BRAT(861)=0.000 ! chi_2c->J/psi gamma',
42  'BRAT(862)=0.798 ! chi_2c->rndmflav rndmflavbar',
43  'BRAT(1501)=0.000 ! chi_0c->J/psi gamma',
44  'BRAT(1502)=0.987 ! chi_0c->rndmflav rndmflavbar',
45  'BRAT(1555)=0.000 ! chi_1c->J/psi gamma',
46  'BRAT(1556)=0.644 ! chi_1c->rndmflav rndmflavbar',
47  'BRAT(1569)=0.186600 ! psi(2S) -> rndmflav rndmflavbar',
48  'BRAT(1570)=0.000 ! psi(2S) ->J/psi pi+ pi-',
49  'BRAT(1571)=0.000 ! psi(2S) ->J/psi pi0 pi0',
50  'BRAT(1572)=0.000 ! psi(2S) ->J/psi eta',
51  'BRAT(1573)=0.000 ! psi(2S) ->J/psi pi0'
52  ),
53  # This is a vector of ParameterSet names to be read, in this order
54  parameterSets = cms.vstring('pythiaUESettings',
55  'processParameters',
56  'CSAParameters'),
57  CSAParameters = cms.vstring('CSAMODE = 6 ! cross-section reweighted quarkonia')
58  )
59 )
60 
61 oniafilter = cms.EDFilter("PythiaFilter",
62  Status = cms.untracked.int32(2),
63  MaxEta = cms.untracked.double(1000.0),
64  MinEta = cms.untracked.double(-1000.0),
65  MinPt = cms.untracked.double(0.0),
66  ParticleID = cms.untracked.int32(443)
67 )
68 
69 mumugenfilter = cms.EDFilter("MCParticlePairFilter",
70  Status = cms.untracked.vint32(1, 1),
71  MinPt = cms.untracked.vdouble(0.5, 0.5),
72  MinP = cms.untracked.vdouble(2.7, 2.7),
73  MaxEta = cms.untracked.vdouble(2.5, 2.5),
74  MinEta = cms.untracked.vdouble(-2.5, -2.5),
75  ParticleCharge = cms.untracked.int32(-1),
76  ParticleID1 = cms.untracked.vint32(13),
77  ParticleID2 = cms.untracked.vint32(13)
78 )
79 
80 mugenfilter = cms.EDFilter("MCSingleParticleFilter",
81  Status = cms.untracked.vint32(1,1),
82  MinPt = cms.untracked.vdouble(10.0,10.0),
83  ParticleID = cms.untracked.vint32(13,-13),
84 )
85 
86 ProductionFilterSequence = cms.Sequence(generator*oniafilter*mumugenfilter*~mugenfilter)