import FWCore.ParameterSet.Config as cms externalLHEProducer = cms.EDProducer("ExternalLHEProducer", args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/UL/13TeV/madgraph/V5_2.6.5/sus_sms/StopStop/SMS-StopStop_mStop-1400_slc7_amd64_gcc700_CMSSW_10_6_19_tarball.tar.xz'), 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') ) import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import * # Some paremeters # Width_4Body = 2.000000000E-17 # BR_4 = 80.0E-02 # BR_2 = 100.0E-02 - BR_4 # Width_tot = Width_4Body/BR_4 hcut = 6.5821E-25#GeV c = 2.9979E11 #mm/s #ctau = (hcut/Width_tot) * c ctau = 0.1 Width = (hcut/ctau) * c qcut, tru_eff = 70, 0.243 mstop = 1400 baseSLHATable=""" BLOCK MASS # Mass Spectrum # PDG code mass particle 1000001 1.00000000E+05 # ~d_L 2000001 1.00000000E+05 # ~d_R 1000002 1.00000000E+05 # ~u_L 2000002 1.00000000E+05 # ~u_R 1000003 1.00000000E+05 # ~s_L 2000003 1.00000000E+05 # ~s_R 1000004 1.00000000E+05 # ~c_L 2000004 1.00000000E+05 # ~c_R 1000005 1.00000000E+05 # ~b_1 2000005 1.00000000E+05 # ~b_2 1000006 %MSTOP% # ~t_1 2000006 1.00000000E+05 # ~t_2 1000011 1.00000000E+05 # ~e_L 2000011 1.00000000E+05 # ~e_R 1000012 1.00000000E+05 # ~nu_eL 1000013 1.00000000E+05 # ~mu_L 2000013 1.00000000E+05 # ~mu_R 1000014 1.00000000E+05 # ~nu_muL 1000015 1.00000000E+05 # ~tau_1 2000015 1.00000000E+05 # ~tau_2 1000016 1.00000000E+05 # ~nu_tauL 1000021 1.00000000E+05 # ~g 1000022 1.00000000E+05 # ~chi_10 1000023 1.00000000E+05 # ~chi_20 1000025 1.00000000E+05 # ~chi_30 1000035 1.00000000E+05 # ~chi_40 1000024 1.00000000E+05 # ~chi_1+ 1000037 1.00000000E+05 # ~chi_2+ # DECAY TABLE # PDG Width DECAY 1000001 0.00000000E+00 # sdown_L decays DECAY 2000001 0.00000000E+00 # sdown_R decays DECAY 1000002 0.00000000E+00 # sup_L decays DECAY 2000002 0.00000000E+00 # sup_R decays DECAY 1000003 0.00000000E+00 # sstrange_L decays DECAY 2000003 0.00000000E+00 # sstrange_R decays DECAY 1000004 0.00000000E+00 # scharm_L decays DECAY 2000004 0.00000000E+00 # scharm_R decays DECAY 1000005 0.00000000E+00 # sbottom1 decays DECAY 2000005 0.00000000E+00 # sbottom2 decays DECAY 1000006 %Width% # stop1 decays # 00.00000000E+00 4 1000022 5 -1 2 # dummy allowed decay, in order to turn on off-shell decays (currently not using) # BR NDA ID1 ID2 3.33333333E-01 2 1 -11 3.33333333E-01 2 1 -13 3.33333333E-01 2 1 -15 DECAY 2000006 0.00000000E+00 # stop2 decays DECAY 1000011 0.00000000E+00 # selectron_L decays DECAY 2000011 0.00000000E+00 # selectron_R decays DECAY 1000012 0.00000000E+00 # snu_elL decays DECAY 1000013 0.00000000E+00 # smuon_L decays DECAY 2000013 0.00000000E+00 # smuon_R decays DECAY 1000014 0.00000000E+00 # snu_muL decays DECAY 1000015 0.00000000E+00 # stau_1 decays DECAY 2000015 0.00000000E+00 # stau_2 decays DECAY 1000016 0.00000000E+00 # snu_tauL decays DECAY 1000021 0.00000000E+00 # gluino decays DECAY 1000022 0.00000000E+00 # neutralino1 decays DECAY 1000023 0.00000000E+00 # neutralino2 decays DECAY 1000024 0.00000000E+00 # chargino1+ decays DECAY 1000025 0.00000000E+00 # neutralino3 decays DECAY 1000035 0.00000000E+00 # neutralino4 decays DECAY 1000037 0.00000000E+00 # chargino2+ decays """ slhatable = baseSLHATable.replace('%MSTOP%','%e' % mstop) slhatable = slhatable.replace('%Width%','%e' % Width) basePythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CP5SettingsBlock, processParameters = cms.vstring( 'LesHouches:setLifetime = 2', '1000006:tau0 = %.1f' % ctau, 'RHadrons:allow = on', 'RHadrons:allowDecay = on', 'RHadrons:setMasses = on', ), 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 = %.0f' % qcut, #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', 'Check:abortIfVeto = on', ), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CP5Settings', 'JetMatchingParameters', 'processParameters' ) ) generator = cms.EDFilter("Pythia8ConcurrentHadronizerFilter", maxEventsToPrint = cms.untracked.int32(1), pythiaPylistVerbosity = cms.untracked.int32(1), pythiaHepMCVerbosity = cms.untracked.bool(False), comEnergy = cms.double(13000.), SLHATableForPythia8 = cms.string('%s' % slhatable), PythiaParameters = basePythiaParameters, ) # Filter setup # ------------------------ # any gen level filters go here ProductionFilterSequence = cms.Sequence(generator) # Link to generator fragment: # genFragments/Hadronizer/13TeV/DisplacedSUSY/StopToLD/DisplacedSUSY_stopToLD_M_1400_0p1mm_TuneCP5_13TeV_pythia8_cff.py