import FWCore.ParameterSet.Config as cms

# link to card:
# https://github.com/cms-sw/genproductions/tree/master/bin/Powheg/production/2017/13TeV/Higgs/VBF_H_NNPDF31_13TeV

externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
                                     args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/UL/13TeV/powheg/V2/VBF_H_slc7_amd64_gcc820_CMSSW_10_6_20_VBF_H_M150/v1/VBF_H_slc7_amd64_gcc820_CMSSW_10_6_20_VBF_H_M150.tgz'),
    nEvents = cms.untracked.uint32(5000),
    numberOfParameters = cms.uint32(1),
    outputFile = cms.string('cmsgrid_final.lhe'),
    generateConcurrently = cms.untracked.bool(True),
    scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh')
)

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from Configuration.Generator.PSweightsPythia.PythiaPSweightsSettings_cfi import *
from Configuration.Generator.Pythia8PowhegEmissionVetoSettings_cfi import *

generator = cms.EDFilter("Pythia8ConcurrentHadronizerFilter",
                                 comEnergy = cms.double(13000.0),
                       	         filterEfficiency = cms.untracked.double(1),
                       	         maxEventsToPrint = cms.untracked.int32(1),
                       	         pythiaHepMCVerbosity = cms.untracked.bool(False),
                       	         pythiaPylistVerbosity = cms.untracked.int32(1),
                       	         PythiaParameters = cms.PSet(
      					 pythia8CommonSettingsBlock,
      					 pythia8CP5SettingsBlock,
                         pythia8PSweightsSettingsBlock,
					 	 pythia8PowhegEmissionVetoSettingsBlock,
                                         processParameters = cms.vstring('POWHEG:nFinal = 3',
                                                                         'SpaceShower:dipoleRecoil = on',
                                                                         'TauDecays:mode = 2', 
                                                                         'TauDecays:tauPolarization = 0', 
                                                                         'TauDecays:tauMother = 25', 
                                                                         '25:m0 = 150.0', 
                                                                         '25:mWidth = 0.0175', 
                                                                         '25:addChannel 1 0.1 100 11 -15', 
                                                                         '25:addChannel 1 0.1 100 15 -11', 
                                                                         '25:onMode = off', 
                                                                         '25:onIfMatch 11 15'),
                                         
      					parameterSets = cms.vstring(
                           			    'pythia8PowhegEmissionVetoSettings',
                            			'pythia8CommonSettings',
                            			'pythia8CP5Settings',
                                        'pythia8PSweightsSettings',
                            			'processParameters',
                            			)
      					)
			)

ProductionFilterSequence = cms.Sequence(generator)