|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
comEnergy = cms.double(14000.0),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'Charmonium:states(3S1) = 100443',
'Charmonium:O(3S1)[3S1(1)] = 0.76',
'Charmonium:O(3S1)[3S1(8)] = 0.0050',
'Charmonium:O(3S1)[1S0(8)] = 0.004',
'Charmonium:O(3S1)[3P0(8)] = 0.004',
'Charmonium:gg2ccbar(3S1)[3S1(1)]g = on',
'Charmonium:gg2ccbar(3S1)[3S1(1)]gm = on',
'Charmonium:gg2ccbar(3S1)[3S1(8)]g = on',
'Charmonium:qg2ccbar(3S1)[3S1(8)]q = on',
'Charmonium:qqbar2ccbar(3S1)[3S1(8)]g = on',
'Charmonium:gg2ccbar(3S1)[1S0(8)]g = on',
'Charmonium:qg2ccbar(3S1)[1S0(8)]q = on',
'Charmonium:qqbar2ccbar(3S1)[1S0(8)]g = on',
'Charmonium:gg2ccbar(3S1)[3PJ(8)]g = on',
'Charmonium:qg2ccbar(3S1)[3PJ(8)]q = on',
'Charmonium:qqbar2ccbar(3S1)[3PJ(8)]g = on',
'100443:onMode = off',
'100443:onIfMatch = 443 211 -211',
'443:onMode = off',
'443:onIfMatch = 13 -13',
'PhaseSpace:pTHatMin = 10.'
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)
)
)
# Filter with high pT cut on dimuon, trying to accomodate trigger requirements.
psi2SIDfilter = cms.EDFilter("PythiaFilter",
ParticleID = cms.untracked.int32(100443),
MinPt = cms.untracked.double(0.0),
MinEta = cms.untracked.double(-2.4),
MaxEta = cms.untracked.double(2.4),
Status = cms.untracked.int32(2)
)
jpsifilter = cms.EDFilter("PythiaFilter",
ParticleID = cms.untracked.int32(443),
MinPt = cms.untracked.double(0.0),
MinEta = cms.untracked.double(-2.4),
MaxEta = cms.untracked.double(2.4),
Status = cms.untracked.int32(2)
)
# Next two muon filter are derived from muon reconstruction
muminusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.),
ParticleID = cms.untracked.int32(443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(-2.4),
MaxEta = cms.untracked.vdouble(2.4),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(-13)
)
muplusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(2.),
ParticleID = cms.untracked.int32(443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(-2.4),
MaxEta = cms.untracked.vdouble(2.4),
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(13)
)
# two pion filter
piminusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(0.0),
ParticleID = cms.untracked.int32(100443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(-2.4), # or 3.0 ?
MaxEta = cms.untracked.vdouble(2.4), # or 3.0 ?
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(-211)
)
piplusfilter = cms.EDFilter("PythiaDauVFilter",
MotherID = cms.untracked.int32(0),
MinPt = cms.untracked.vdouble(0.0),
ParticleID = cms.untracked.int32(100443),
ChargeConjugation = cms.untracked.bool(False),
MinEta = cms.untracked.vdouble(-2.4), # or 3.0 ?
MaxEta = cms.untracked.vdouble(2.4), # or 3.0 ?
NumberDaughters = cms.untracked.int32(1),
DaughterIDs = cms.untracked.vint32(211)
)
ProductionFilterSequence = cms.Sequence(generator*psi2SIDfilter*jpsifilter*muminusfilter*muplusfilter*piminusfilter*piplusfilter)
|