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.75) # 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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS125.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.00259608107509008,'acceptance':0.412},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS125.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.00259608107509008,'acceptance':0.613},
{'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 3.000 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.75,'name':'SUEP_mS125.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.00259608107509008,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS200.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.0048232274473043275,'acceptance':0.423},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS200.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.0048232274473043275,'acceptance':0.6244},
{'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 3.000 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.75,'name':'SUEP_mS200.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.0048232274473043275,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS300.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.009035833057989004,'acceptance':0.5},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS300.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.009035833057989004,'acceptance':0.6757},
{'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 3.000 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.75,'name':'SUEP_mS300.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.009035833057989004,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS400.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.01466909822798684,'acceptance':0.518},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS400.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.01466909822798684,'acceptance':0.652},
{'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 3.000 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.75,'name':'SUEP_mS400.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.01466909822798684,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS500.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.021723022957297833,'acceptance':0.5146},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS500.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.021723022957297833,'acceptance':0.6577},
{'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 3.000 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.75,'name':'SUEP_mS500.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.021723022957297833,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS600.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.030197607245921985,'acceptance':0.5111},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS600.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.030197607245921985,'acceptance':0.6634},
{'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 3.000 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.75,'name':'SUEP_mS600.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.030197607245921985,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS700.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.0400928510938593,'acceptance':0.5077},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS700.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.0400928510938593,'acceptance':0.6691},
{'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 3.000 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.75,'name':'SUEP_mS700.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.0400928510938593,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS800.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.051408754501109766,'acceptance':0.5068},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS800.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.051408754501109766,'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 = 800.000','999999:all = GeneralResonance void 0 0 0 3.000 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.75,'name':'SUEP_mS800.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.051408754501109766,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS900.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.0641453174676734,'acceptance':0.5084},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS900.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.0641453174676734,'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 = 900.000','999999:all = GeneralResonance void 0 0 0 3.000 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.75,'name':'SUEP_mS900.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.0641453174676734,'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 1.0 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 1.0 101 211 -211'],'T':0.75,'name':'SUEP_mS1000.000_mPhi3.000_T0.750_modegeneric','filterEfficiency':0.07830253999355018,'acceptance':0.51},
{'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 3.000 0.001 0.0 0.0 0.0','999998:all = GeneralResonance void 1 0 0 0.5 0.001 0.0 0.0 0.0','999999:oneChannel = 1 1.0 101 999998 999998','999998:oneChannel = 1 0.4 101 11 -11','999998:addChannel = 1 0.4 101 13 -13','999998:addChannel = 1 0.2 101 211 -211'],'T':0.75,'name':'SUEP_mS1000.000_mPhi3.000_T0.750_modeleptonic','filterEfficiency':0.07830253999355018,'acceptance':0.668},
{'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 3.000 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.75,'name':'SUEP_mS1000.000_mPhi3.000_T0.750_modehadronic','filterEfficiency':0.07830253999355018,'acceptance':0.672}
]  # These are the parameters for all generated points
mcm_eff = 0.005828 # 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.75) # 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(1000.) #genHT cut
)

ProductionFilterSequence = cms.Sequence(generator*
                                        tmpGenParticles * tmpGenParticlesForJetsNoNu *
                                        tmpAk4GenJetsNoNu * genHTFilter
)