import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import * from Configuration.Generator.PSweightsPythia.PythiaPSweightsSettings_cfi import * generator = cms.EDFilter("Pythia8GeneratorFilter", maxEventsToPrint = cms.untracked.int32(1), pythiaPylistVerbosity = cms.untracked.int32(1), filterEfficiency = cms.untracked.double(1.0), pythiaHepMCVerbosity = cms.untracked.bool(False), comEnergy = cms.double(13000.0), RandomizedParameters = cms.VPSet(), UserCustomization = cms.VPSet(cms.PSet( pluginName = cms.string("SuepDecay"), idDark = cms.int32(999999), # pdgId of the dark meson idMediator = cms.int32(25), # pdgId of the mediator temperature = cms.double(0.35) # Temperature of the thermal distribution ) ), ) points = [ {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 125.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS125.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.022561693317296093,'acceptance':0.582}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 200.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS200.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.045499775471722015,'acceptance':0.5531}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 300.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS300.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.08115244345939307,'acceptance':0.678}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 400.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS400.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.12259774967384872,'acceptance':0.679}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 500.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS500.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.16983569411508892,'acceptance':0.6793}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 600.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS600.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.22286627678311371,'acceptance':0.6796}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 700.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS700.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.2816894976779231,'acceptance':0.6799}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 800.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS800.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.3463053567995171,'acceptance':0.6784}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 900.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS900.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.4167138541478956,'acceptance':0.6752}, {'processParameters':['Higgs:useBSM = on','HiggsBSM:gg2H1 = on','HiggsH1:coup2d = 1','HiggsH1:coup2u = 0','HiggsH1:coup2Z = 0','HiggsH1:coup2W = 0','HiggsH1:coup2l = 0','25:m0 = 1000.000','999999:all = GeneralResonance void 0 0 0 1.400 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.7 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.15 101 11 -11','999998:addChannel = 1 0.15 101 13 -13','999998:addChannel = 1 0.7 101 211 -211'],'T':0.35,'name':'SUEP_mS1000.000_mPhi1.400_T0.350_modehadronic','filterEfficiency':0.4929149897230587,'acceptance':0.672} ] # These are the parameters for all generated points mcm_eff = 0.086422 # This is the overall mcm efficiency, averaged across all points for point in points: basePythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CP5SettingsBlock, pythia8PSweightsSettingsBlock, processParameters = cms.vstring(point['processParameters']), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CP5Settings', 'pythia8PSweightsSettings', 'processParameters', ) ) baseUserCustomization = cms.VPSet( cms.PSet( pluginName = cms.string("SuepDecay"), idDark = cms.int32(999999), # pdgId of the dark meson idMediator = cms.int32(25), # pdgId of the mediator temperature = cms.double(0.35) # Temperature of the thermal distribution ) ) generator.RandomizedParameters.append( cms.PSet( ConfigWeight = cms.double(1./(point['filterEfficiency']*point['acceptance'])), # I.e. bias events so the number of events in SR is independent of point ConfigDescription = cms.string(point['name']), PythiaParameters = basePythiaParameters, UserCustomization = baseUserCustomization, ), ) # Filter setup # ------------------------ # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/PhysicsTools/HepMCCandAlgos/python/genParticles_cfi.py tmpGenParticles = cms.EDProducer("GenParticleProducer", saveBarCodes = cms.untracked.bool(True), src = cms.InputTag("generator","unsmeared"), abortOnUnknownPDGCode = cms.untracked.bool(False) ) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/Configuration/python/GenJetParticles_cff.py # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoMET/Configuration/python/GenMETParticles_cff.py tmpGenParticlesForJetsNoNu = cms.EDProducer("InputGenJetsParticleSelector", src = cms.InputTag("tmpGenParticles"), ignoreParticleIDs = cms.vuint32( 1000022, 1000012, 1000014, 1000016, 2000012, 2000014, 2000016, 1000039, 5100039, 4000012, 4000014, 4000016, 9900012, 9900014, 9900016, 39,12,14,16), partonicFinalState = cms.bool(False), excludeResonances = cms.bool(False), excludeFromResonancePids = cms.vuint32(12, 13, 14, 16), tausAsJets = cms.bool(False) ) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/JetProducers/python/AnomalousCellParameters_cfi.py AnomalousCellParameters = cms.PSet( maxBadEcalCells = cms.uint32(9999999), maxRecoveredEcalCells = cms.uint32(9999999), maxProblematicEcalCells = cms.uint32(9999999), maxBadHcalCells = cms.uint32(9999999), maxRecoveredHcalCells = cms.uint32(9999999), maxProblematicHcalCells = cms.uint32(9999999) ) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/JetProducers/python/GenJetParameters_cfi.py GenJetParameters = cms.PSet( src = cms.InputTag("tmpGenParticlesForJetsNoNu"), srcPVs = cms.InputTag(''), jetType = cms.string('GenJet'), jetPtMin = cms.double(3.0), inputEtMin = cms.double(0.0), inputEMin = cms.double(0.0), doPVCorrection = cms.bool(False), # pileup with offset correction doPUOffsetCorr = cms.bool(False), # if pileup is false, these are not read: nSigmaPU = cms.double(1.0), radiusPU = cms.double(0.5), # fastjet-style pileup doAreaFastjet = cms.bool(False), doRhoFastjet = cms.bool(False), # if doPU is false, these are not read: Active_Area_Repeats = cms.int32(5), GhostArea = cms.double(0.01), Ghost_EtaMax = cms.double(6.0), Rho_EtaMax = cms.double(4.5), useDeterministicSeed= cms.bool( True ), minSeed = cms.uint32( 14327 ) ) # https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/RecoJets/JetProducers/python/ak4GenJets_cfi.py tmpAk4GenJetsNoNu = cms.EDProducer( "FastjetJetProducer", GenJetParameters, AnomalousCellParameters, jetAlgorithm = cms.string("AntiKt"), rParam = cms.double(0.4) ) genHTFilter = cms.EDFilter("GenHTFilter", src = cms.InputTag("tmpAk4GenJetsNoNu"), #GenJet collection as input jetPtCut = cms.double(30.0), #GenJet pT cut for HT jetEtaCut = cms.double(2.5), #GenJet eta cut for HT genHTcut = cms.double(400.) #genHT cut ) ProductionFilterSequence = cms.Sequence(generator* tmpGenParticles * tmpGenParticlesForJetsNoNu * tmpAk4GenJetsNoNu * genHTFilter )