Generator/python/BsToPhiPhi_KKKK_14TeV_TuneCP5_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0004 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0005 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0006
0007 generator = cms.EDFilter("Pythia8GeneratorFilter",
0008     pythiaHepMCVerbosity = cms.untracked.bool(False),
0009     maxEventsToPrint = cms.untracked.int32(0),
0010     pythiaPylistVerbosity = cms.untracked.int32(0),
0011     filterEfficiency = cms.untracked.double(1.38e-3),
0012     crossSection = cms.untracked.double(540000000.),
0013     comEnergy = cms.double(14000.0),
0014     ExternalDecays = cms.PSet(
0015         EvtGen130 = cms.untracked.PSet(
0016             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
0017             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0018             user_decay_embedded= cms.vstring(
0019                 'Define Hp 0.49',
0020                 'Define Hz 0.775',
0021                 'Define Hm 0.4',
0022                 'Define pHp 2.50',
0023                 'Define pHz 0.0',
0024                 'Define pHm -0.17',
0025                 '#',
0026                 'Alias      MyB_s0   B_s0',
0027                 'Alias      Myanti-B_s0   anti-B_s0',
0028                 'ChargeConj Myanti-B_s0   MyB_s0',
0029                 'Alias      MyPhi    phi',
0030                 'ChargeConj MyPhi    MyPhi',
0031                 '#',
0032                 'Decay MyB_s0',
0033                 '  1.000         MyPhi      MyPhi        PVV_CPLH 0.02 1 Hp pHp Hz pHz Hm pHm;',
0034                 '#',
0035                 'Enddecay',
0036                 'Decay Myanti-B_s0',
0037                  ' 1.000         MyPhi      MyPhi        PVV_CPLH 0.02 1 Hp pHp Hz pHz Hm pHm;',
0038                 'Enddecay',
0039                 '#',
0040                 'Decay MyPhi',
0041                  ' 1.000         K+          K-           VSS;',
0042                 'Enddecay',
0043                 'End'
0044             ),
0045             list_forced_decays = cms.vstring('MyB_s0','Myanti-B_s0'),
0046             operates_on_particles = cms.vint32(),
0047         ),
0048         parameterSets = cms.vstring('EvtGen130')
0049     ),
0050     PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
0051         pythia8CP5SettingsBlock,
0052         processParameters = cms.vstring(
0053             "SoftQCD:nonDiffractive = on",
0054             'PTFilter:filter = on', # this turn on the filter

0055             'PTFilter:quarkToFilter = 5', # PDG id of q quark

0056             'PTFilter:scaleToFilter = 1.0'
0057         ),
0058         parameterSets = cms.vstring('pythia8CommonSettings',
0059             'pythia8CP5Settings',
0060             'processParameters',
0061         )
0062     )
0063 )
0064
0065 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0066
0067 configurationMetadata = cms.untracked.PSet(
0068     version = cms.untracked.string('$Revision: 1.0 $'),
0069     name = cms.untracked.string('$Source: /Configuration/Generator/python/BsPhiPhi_KKKK_PhaseII_cfi.py $'),
0070     annotation = cms.untracked.string('PhaseII: Pythia8+EvtGen130 generation of Bs --> phi phi (4K), 14TeV, Tune CP5')
0071 )
0072
0073 bfilter = cms.EDFilter(
0074     "PythiaFilter",
0075     MaxEta = cms.untracked.double(9999.),
0076     MinEta = cms.untracked.double(-9999.),
0077     ParticleID = cms.untracked.int32(531)
0078 )
0079
0080 decayfilter = cms.EDFilter(
0081     "PythiaDauVFilter",
0082     verbose         = cms.untracked.int32(1),
0083     NumberDaughters = cms.untracked.int32(2),
0084     ParticleID      = cms.untracked.int32(531),
0085     DaughterIDs     = cms.untracked.vint32(333, 333),  ## phi phi

0086     MinPt           = cms.untracked.vdouble(3.0, 3.0),
0087     MinEta          = cms.untracked.vdouble(-3.0, -3.0),
0088     MaxEta          = cms.untracked.vdouble( 3.0,  3.0)
0089 )
0090
0091 phifilter = cms.EDFilter(
0092     "PythiaDauVFilter",
0093     verbose = cms.untracked.int32(0),
0094     NumberDaughters = cms.untracked.int32(2),
0095     MotherID = cms.untracked.int32(531),
0096     ParticleID = cms.untracked.int32(333),
0097     DaughterIDs = cms.untracked.vint32(321, -321),
0098     MinPt = cms.untracked.vdouble(1.95, 1.95),
0099     MinEta = cms.untracked.vdouble(-2.5, -2.5),
0100     MaxEta = cms.untracked.vdouble(2.5, 2.5)
0101 )
0102
0103 ProductionFilterSequence = cms.Sequence(generator*bfilter*decayfilter*phifilter)