import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import * import math #import os,sys,math #sys.path.append(os.path.dirname(os.path.dirname(os.path.abspath(__file__)))) baseSLHATable=""" BLOCK MODSEL # Model selection 1 1 sugra # BLOCK QNUMBERS 3000001 # f 1 0 # 3 times electric charge 2 2 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # Particle/Antiparticle distinction (0=own anti) # BLOCK QNUMBERS 3000002 # n 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # Particle/Antiparticle distinction (0=own anti) BLOCK MASS 2000001 1.00000000E+10 2000002 1.00000000E+10 2000003 1.00000000E+10 2000004 1.00000000E+10 2000005 1.00000000E+10 2000006 1.00000000E+10 2000011 1.00000000E+10 2000013 1.00000000E+10 2000015 1.00000000E+10 1000001 1.00000000E+10 1000002 1.00000000E+10 1000003 1.00000000E+10 1000004 1.00000000E+10 1000005 1.00000000E+10 1000006 1.00000000E+10 1000011 1.00000000E+10 1000012 1.00000000E+10 1000013 1.00000000E+10 1000014 1.00000000E+10 1000015 1.00000000E+10 1000016 1.00000000E+10 1000021 2000 1000022 1000 1000023 1.00000000E+10 1000024 1000 1000025 1.00000000E+10 1000035 1.00000000E+10 1000037 1.00000000E+10 1000039 0.0 # Gravitino 3000001 1.00000000E+02 # f: stealth sector 3000002 9.00000000E+01 # n: stealth sector # DECAY 2000001 0.00000000E+00 DECAY 2000002 0.00000000E+00 DECAY 2000003 0.00000000E+00 DECAY 2000004 0.00000000E+00 DECAY 2000005 0.00000000E+00 DECAY 2000006 0.00000000E+00 DECAY 2000011 0.00000000E+00 DECAY 2000013 0.00000000E+00 DECAY 2000015 0.00000000E+00 DECAY 1000001 0.00000000E+00 DECAY 1000002 0.00000000E+00 DECAY 1000003 0.00000000E+00 DECAY 1000004 0.00000000E+00 DECAY 1000005 0.00000000E+00 DECAY 1000006 0.00000000E+00 DECAY 1000011 0.00000000E+00 DECAY 1000012 0.00000000E+00 DECAY 1000013 0.00000000E+00 DECAY 1000014 0.00000000E+00 DECAY 1000015 0.00000000E+00 DECAY 1000016 0.00000000E+00 DECAY 1000021 1.00000000E+00 0.00000000E+00 3 -6 6 1000022 1.00000000E-01 3 -1 1 1000022 1.00000000E-01 3 -2 2 1000022 1.00000000E-01 3 -3 3 1000022 1.00000000E-01 3 -4 4 1000022 1.00000000E-01 3 -5 5 1000022 1.25000000E-01 3 -2 1 1000024 1.25000000E-01 3 2 -1 -1000024 1.25000000E-01 3 -4 3 1000024 1.25000000E-01 3 4 -3 -1000024 DECAY 1000022 1.00000000E-03 0.0000000 2 3000001 22 1.0000000 2 3000001 22 DECAY 1000024 1.00000000E-03 0.0000000 3 3000001 -1 2 1.0000000 2 3000001 24 DECAY 3000001 1.00000000E-06 # f decays 1.00000000e+00 2 1000039 3000002 # BR(f -> gs n ) DECAY 3000002 1.00000000E-06 # n decays 1.00000000e+00 2 21 21 # BR(n -> g g ) DECAY 1000039 0.00000000E+00 """ externalLHEProducer = cms.EDProducer("ExternalLHEProducer", args = cms.vstring("/cvmfs/cms.cern.ch/phys_generator/gridpacks/UL/13TeV/madgraph/V5_2.6.5/sus_sms/SMS-GlGl_v3/SMS-GlGl_mGl-2000_slc7_amd64_gcc700_CMSSW_10_6_19_tarball.tar.xz"), nEvents = cms.untracked.uint32(500), numberOfParameters = cms.uint32(1), outputFile = cms.string('cmsgrid_final.lhe'), scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh'), generateConcurrently = cms.untracked.bool(False) ) generator = cms.EDFilter("Pythia8HadronizerFilter", 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.), PythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CP5SettingsBlock, JetMatchingParameters = cms.vstring( 'JetMatching:setMad = off', 'JetMatching:scheme = 1', 'JetMatching:merge = on', 'JetMatching:jetAlgorithm = 2', 'JetMatching:etaJetMax = 5.', 'JetMatching:coneRadius = 1.', 'JetMatching:slowJetPower = 1', 'JetMatching:qCut = 156', #this is the actual merging scale 'JetMatching:nQmatch = 5', #4 corresponds to 4-flavour scheme (no matching of b-quarks), 5 for 5-flavour scheme 'JetMatching:nJetMax = 2', #number of partons in born matrix element for highest multiplicity 'JetMatching:doShowerKt = off', #off for MLM matching, turn on for shower-kT matching '6:m0 = 172.5', '24:mMin = 0.1', # allows W offshelness to go down to 7 GeV (very low LSP points in the scan) 'Check:abortIfVeto = on', ), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CP5Settings', 'JetMatchingParameters' ) ), SLHATableForPythia8 = cms.string(baseSLHATable), )