Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-08-18 22:30:03

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0003 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0004 from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *
0005 
0006 generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
0007   maxEventsToPrint = cms.untracked.int32(1),
0008   pythiaPylistVerbosity = cms.untracked.int32(1),
0009   filterEfficiency = cms.untracked.double(0.4045),
0010   pythiaHepMCVerbosity = cms.untracked.bool(False),
0011   comEnergy = cms.double(14000.),
0012   PythiaParameters = cms.PSet(
0013     pythia8CommonSettingsBlock,
0014     pythia8CP5SettingsBlock,
0015     processParameters = cms.vstring(
0016       'Higgs:useBSM = on',
0017       'HiggsBSM:gg2H2 = on',
0018       'HiggsH2:coup2d = 10.0',
0019       'HiggsH2:coup2u = 10.0',
0020       'HiggsH2:coup2Z = 0.0',
0021       'HiggsH2:coup2W = 0.0',
0022       'HiggsA3:coup2H2Z = 0.0',
0023       'HiggsH2:coup2A3A3 = 0.0',
0024       'HiggsH2:coup2H1H1 = 0.0',
0025       '443:onMode = off',
0026       '443:onIfMatch 13 -13',
0027       ############# Fix to mass of etaB (9.4 GeV) #############
0028       '35:mMin = 0',
0029       '35:mMax = 50.0',
0030       '35:m0 = 9.4',
0031       '35:mWidth = 0.00',
0032       '35:addChannel 1 1.00 100 443 443',
0033       '35:onMode = off',
0034       '35:onIfMatch 443 443' #Jpsi Jpsi decay
0035     ),
0036       # This is a vector of ParameterSet names to be read, in this order
0037     parameterSets = cms.vstring(
0038       'pythia8CommonSettings',
0039       'pythia8CP5Settings',
0040       'processParameters'
0041     )
0042   )
0043 )
0044 
0045 generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)
0046 
0047 etafilter = cms.EDFilter("PythiaFilter",
0048   MaxEta = cms.untracked.double(9999.0),
0049   MinEta = cms.untracked.double(-9999.0),
0050   ParticleID = cms.untracked.int32(35)
0051 )
0052 
0053 etatojpsipairfilter = cms.EDFilter("PythiaDauVFilter",
0054   MotherID = cms.untracked.int32(0),
0055   verbose = cms.untracked.int32(0),
0056   ParticleID = cms.untracked.int32(35),
0057   MaxEta = cms.untracked.vdouble(2.6, 2.6),
0058   MinEta = cms.untracked.vdouble(-2.6, -2.6),
0059   DaughterIDs = cms.untracked.vint32(443, 443),
0060   MinPt = cms.untracked.vdouble(0., 0.),
0061   NumberDaughters = cms.untracked.int32(2)
0062 )
0063 
0064 jpsifilter = cms.EDFilter("PythiaDauVFilter",
0065   MotherID = cms.untracked.int32(35),
0066   verbose = cms.untracked.int32(0),
0067   ParticleID = cms.untracked.int32(443),
0068   MaxEta = cms.untracked.vdouble(2.5, 2.5),
0069   MinEta = cms.untracked.vdouble(-2.5, -2.5),
0070   DaughterIDs = cms.untracked.vint32(13, -13),
0071   MinPt = cms.untracked.vdouble(1.8, 1.8),
0072   NumberDaughters = cms.untracked.int32(2)
0073 )
0074 
0075 ProductionFilterSequence = cms.Sequence(generator*etafilter*etatojpsipairfilter*jpsifilter)