Generator/python/Psi2SToJPsiPiPi_14TeV_TuneCP5_pythia8_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0003 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0004
0005 generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
0006     pythiaPylistVerbosity = cms.untracked.int32(0),
0007     pythiaHepMCVerbosity = cms.untracked.bool(False),
0008     maxEventsToPrint = cms.untracked.int32(0),
0009     comEnergy = cms.double(14000.0),
0010     PythiaParameters = cms.PSet(
0011         pythia8CommonSettingsBlock,
0012         pythia8CP5SettingsBlock,
0013         processParameters = cms.vstring(
0014                 'Charmonium:states(3S1) = 100443',
0015                 'Charmonium:O(3S1)[3S1(1)] = 0.76',
0016                 'Charmonium:O(3S1)[3S1(8)] = 0.0050',
0017                 'Charmonium:O(3S1)[1S0(8)] = 0.004',
0018                 'Charmonium:O(3S1)[3P0(8)] = 0.004',
0019                 'Charmonium:gg2ccbar(3S1)[3S1(1)]g = on',
0020                 'Charmonium:gg2ccbar(3S1)[3S1(1)]gm = on',
0021                 'Charmonium:gg2ccbar(3S1)[3S1(8)]g = on',
0022                 'Charmonium:qg2ccbar(3S1)[3S1(8)]q = on',
0023                 'Charmonium:qqbar2ccbar(3S1)[3S1(8)]g = on',
0024                 'Charmonium:gg2ccbar(3S1)[1S0(8)]g = on',
0025                 'Charmonium:qg2ccbar(3S1)[1S0(8)]q = on',
0026                 'Charmonium:qqbar2ccbar(3S1)[1S0(8)]g = on',
0027                 'Charmonium:gg2ccbar(3S1)[3PJ(8)]g = on',
0028                 'Charmonium:qg2ccbar(3S1)[3PJ(8)]q = on',
0029                 'Charmonium:qqbar2ccbar(3S1)[3PJ(8)]g = on',
0030                 '100443:onMode = off',
0031                 '100443:onIfMatch = 443 211 -211',
0032                 '443:onMode = off',
0033                 '443:onIfMatch = 13 -13',
0034                 'PhaseSpace:pTHatMin = 10.'
0035             ),
0036         parameterSets = cms.vstring('pythia8CommonSettings',
0037                                     'pythia8CP5Settings',
0038                                     'processParameters',
0039                                     )
0040     )
0041 )
0042
0043 # Filter with high pT cut on dimuon, trying to accomodate trigger requirements.
0044
0045 psi2SIDfilter = cms.EDFilter("PythiaFilter",
0046     ParticleID = cms.untracked.int32(100443),
0047     MinPt = cms.untracked.double(0.0),
0048     MinEta = cms.untracked.double(-2.4),
0049     MaxEta = cms.untracked.double(2.4),
0050     Status = cms.untracked.int32(2)
0051 )
0052
0053 jpsifilter = cms.EDFilter("PythiaFilter",
0054     ParticleID = cms.untracked.int32(443),
0055     MinPt = cms.untracked.double(0.0),
0056     MinEta = cms.untracked.double(-2.4),
0057     MaxEta = cms.untracked.double(2.4),
0058     Status = cms.untracked.int32(2)
0059 )
0060
0061 # Next two muon filter are derived from muon reconstruction
0062
0063 muminusfilter = cms.EDFilter("PythiaDauVFilter",
0064     MotherID = cms.untracked.int32(0),
0065     MinPt = cms.untracked.vdouble(2.),
0066     ParticleID = cms.untracked.int32(443),
0067     ChargeConjugation = cms.untracked.bool(False),
0068     MinEta = cms.untracked.vdouble(-2.4),
0069     MaxEta = cms.untracked.vdouble(2.4),
0070     NumberDaughters = cms.untracked.int32(1),
0071     DaughterIDs = cms.untracked.vint32(-13)
0072 )
0073
0074 muplusfilter = cms.EDFilter("PythiaDauVFilter",
0075     MotherID = cms.untracked.int32(0),
0076     MinPt = cms.untracked.vdouble(2.),
0077     ParticleID = cms.untracked.int32(443),
0078     ChargeConjugation = cms.untracked.bool(False),
0079     MinEta = cms.untracked.vdouble(-2.4),
0080     MaxEta = cms.untracked.vdouble(2.4),
0081     NumberDaughters = cms.untracked.int32(1),
0082     DaughterIDs = cms.untracked.vint32(13)
0083 )
0084
0085 #  two pion filter
0086 piminusfilter = cms.EDFilter("PythiaDauVFilter",
0087     MotherID = cms.untracked.int32(0),
0088     MinPt = cms.untracked.vdouble(0.0),
0089     ParticleID = cms.untracked.int32(100443),
0090     ChargeConjugation = cms.untracked.bool(False),
0091     MinEta = cms.untracked.vdouble(-2.4), # or 3.0 ?
0092     MaxEta = cms.untracked.vdouble(2.4), # or 3.0 ?
0093     NumberDaughters = cms.untracked.int32(1),
0094     DaughterIDs = cms.untracked.vint32(-211)
0095 )
0096
0097 piplusfilter = cms.EDFilter("PythiaDauVFilter",
0098     MotherID = cms.untracked.int32(0),
0099     MinPt = cms.untracked.vdouble(0.0),
0100     ParticleID = cms.untracked.int32(100443),
0101     ChargeConjugation = cms.untracked.bool(False),
0102     MinEta = cms.untracked.vdouble(-2.4), # or 3.0 ?
0103     MaxEta = cms.untracked.vdouble(2.4), # or 3.0 ?
0104     NumberDaughters = cms.untracked.int32(1),
0105     DaughterIDs = cms.untracked.vint32(211)
0106 )
0107
0108 ProductionFilterSequence = cms.Sequence(generator*psi2SIDfilter*jpsifilter*muminusfilter*muplusfilter*piminusfilter*piplusfilter)