/
githubmirror
/
cmssw
Обзор
Документация
Войти
/
githubmirror
/
cmssw
Код
Запросы
0
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
GeneratorInterface/Pythia8Interface/test/randomizedParametersSLHA.py
565 строк
21 KB
Josh Bendavid
update examples
24 фев 2016, 11:06
24 фев 2016, 11:06
e4f1593
Код
Авторство
О чём код?
#example gen fragment where a randomized parameter scan is performed over parameters #in the SLHA table #since the randomization occurs only at lumi section boundaries, test with command like #cmsDriver.py randomizedParametersSLHA.py -s GEN --conditions auto:mc -n -1 --filein file:cmsgrid_final.lhe --eventcontent AODSIM --no_exec --customise_command "process.source.numberEventsInLuminosityBlock = cms.untracked.uint32(200)" import FWCore.ParameterSet.Config as cms from Configuration.Generator.Pythia8CommonSettings_cfi import * from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import * basePythiaParameters = cms.PSet( pythia8CommonSettingsBlock, pythia8CUEP8M1SettingsBlock, 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 = 58', #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 ), parameterSets = cms.vstring('pythia8CommonSettings', 'pythia8CUEP8M1Settings', 'JetMatchingParameters' ) ) baseSLHATable=""" ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ################################### ## INFORMATION FOR HMIX ################################### BLOCK HMIX Q= 1.000000e+00 # 1 3.576810e+02 # rmuh 2 9.748624e+00 # tb 4 1.664391e+05 # ma2 3 2.449353e+02 # higgs vev(Q) MSSM DRb ################################### ## INFORMATION FOR MASS ################################### BLOCK MASS # 5 0.000000e+00 # mb 6 1.730000e+02 # mt 15 1.777000e+00 # mta 23 9.118760e+01 # mz 24 80.419002 # mw 25 1.25e02 # mh01 35 1.00e05 # mh02 36 1.00e05 # ma0 37 1.00e05 # mh 1000001 1.00e05 # set of param :1*msd1, 1*msd2 1000002 1.00e05 # set of param :1*msu1, 1*msu2 1000005 4.70e02 # msd3 1000006 1.00e05 # msu3 1000011 1.00e05 # set of param :1*msl1, 1*msl2 1000012 1.00e05 # set of param :1*msn1, 1*msn2 1000015 1.00e05 # msl3 1000016 1.00e05 # msn3 1000021 1.00e05 # mgo 1000022 %MLSP% # mneu1 1000023 %MNLSP% # mneu2 1000024 1.00e05 # mch1 1000025 -1.00e05 # mneu3 1000035 1.00e05 # mneu4 1000037 1.00e05 # mch2 2000001 1.00e05 # set of param :1*msd4, 1*msd5 2000002 1.00e05 # set of param :1*msu4, 1*msu5 2000005 1.00e05 # msd6 2000006 1.00e05 # msu6 2000011 1.00e05 # set of param :1*msl4, 1*msl5 2000015 1.0005 # msl6 1 0.000000e00 # d : 0.0 2 0.000000e00 # u : 0.0 3 0.000000e00 # s : 0.0 4 0.000000e00 # c : 0.0 11 0.000000e00 # e- : 0.0 12 0.000000e00 # ve : 0.0 13 0.000000e00 # mu- : 0.0 14 0.000000e00 # vm : 0.0 16 0.000000e00 # vt : 0.0 21 0.000000e00 # g : 0.0 22 0.000000e00 # a : 0.0 1000014 1.00e05 # svm : msn1 1000013 1.00e05 # mul- : msl1 2000013 1.00e05 # mur- : msl4 1000004 1.00e05 # cl : msu1 2000004 1.00e05 # cr : msu4 1000003 1.00e05 # sl : msd1 2000003 1.00e05 # sr : msd4 ################################### ## INFORMATION FOR MSOFT ################################### BLOCK MSOFT Q= 1.000000e+00 # 1 1.013965e+02 # rmx1 2 1.915042e+02 # rmx2 3 5.882630e+02 # rmx3 21 3.233749e+04 # mhd2 22 -1.288001e+05 # mhu2 31 1.953348e+02 # set of param :1*rml21x1, 1*rml22x2 32 1.953348e+02 # mg5 will not use this value use instead 1*mdl_rml21x1 33 1.944960e+02 # rml23x3 34 1.364941e+02 # set of param :1*rme21x1, 1*rme22x2 35 1.364941e+02 # mg5 will not use this value use instead 1*mdl_rme21x1 36 1.340434e+02 # rme23x3 41 5.475735e+02 # set of param :1*rmq21x1, 1*rmq22x2 42 5.475735e+02 # mg5 will not use this value use instead 1*mdl_rmq21x1 43 4.987639e+02 # rmq23x3 44 5.295112e+02 # set of param :1*rmu21x1, 1*rmu22x2 45 5.295112e+02 # mg5 will not use this value use instead 1*mdl_rmu21x1 46 4.232459e+02 # rmu23x3 47 5.231488e+02 # set of param :1*rmd21x1, 1*rmd22x2 48 5.231488e+02 # mg5 will not use this value use instead 1*mdl_rmd21x1 49 5.198673e+02 # rmd23x3 ################################### ## INFORMATION FOR NMIX ################################### BLOCK NMIX # 1 1 9.863644e-01 # rnn1x1 1 2 -5.311036e-02 # rnn1x2 1 3 1.464340e-01 # rnn1x3 1 4 -5.311861e-02 # rnn1x4 2 1 9.935054e-02 # rnn2x1 2 2 9.449493e-01 # rnn2x2 2 3 -2.698467e-01 # rnn2x3 2 4 1.561507e-01 # rnn2x4 3 1 -6.033880e-02 # rnn3x1 3 2 8.770049e-02 # rnn3x2 3 3 6.958775e-01 # rnn3x3 3 4 7.102270e-01 # rnn3x4 4 1 -1.165071e-01 # rnn4x1 4 2 3.107390e-01 # rnn4x2 4 3 6.492260e-01 # rnn4x3 4 4 -6.843778e-01 # rnn4x4 ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.279340e+02 # aewm1 3 1.180000e-01 # as 6 1.750000e+02 # mt 7 1.777000e+00 # mta 4 9.118760e+01 # mz 2 1.166370e-05 # G_F [GeV^-2] ################################### ## INFORMATION FOR UMIX ################################### BLOCK UMIX # 1 1 9.168349e-01 # ruu1x1 1 2 -3.992666e-01 # ruu1x2 2 1 3.992666e-01 # ruu2x1 2 2 9.168349e-01 # ruu2x2 ################################### ## INFORMATION FOR VMIX ################################### BLOCK VMIX # 1 1 9.725578e-01 # rvv1x1 1 2 -2.326612e-01 # rvv1x2 2 1 2.326612e-01 # rvv2x1 2 2 9.725578e-01 # rvv2x2 ################################### ## INFORMATION FOR YD ################################### BLOCK YD Q= 1.000000e+00 # 3 3 0.000000e+00 # ryd3x3 ################################### ## INFORMATION FOR YE ################################### BLOCK YE Q= 1.000000e+00 # 3 3 1.008908e-01 # rye3x3 ################################### ## INFORMATION FOR YU ################################### BLOCK YU Q= 1.000000e+00 # 3 3 8.928445e-01 # ryu3x3 ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.491500E+00 # wt DECAY 23 2.441404e+00 # wz DECAY 24 2.047600e+00 # ww DECAY 25 6.382339e-03 # wh01 DECAY 35 5.748014e-01 # wh02 DECAY 36 6.321785e-01 # wa0 DECAY 37 5.469628e-01 # wh DECAY 1000001 5.312788e+00 # wsd1 DECAY 1000002 5.477195e+00 # wsu1 DECAY 1000003 5.312788e+00 # wsd2 DECAY 1000004 5.477195e+00 # wsu2 DECAY 1000005 1.000000e+00 # wsd3 1.00000000E+00 2 1000023 5 DECAY 1000006 2.021596e+00 # wsu3 DECAY 1000011 2.136822e-01 # wsl1 DECAY 1000012 1.498816e-01 # wsn1 DECAY 1000013 2.136822e-01 # wsl2 DECAY 1000014 1.498816e-01 # wsn2 DECAY 1000015 1.483273e-01 # wsl3 DECAY 1000016 1.475190e-01 # wsn3 DECAY 1000021 5.506754e+00 # wgo DECAY 1000023 1.000000e-01 # wneu2 1.00000000E+00 2 1000022 25 DECAY 1000024 1.704145e-02 # wch1 DECAY 1000025 1.915985e+00 # wneu3 DECAY 1000035 2.585851e+00 # wneu4 DECAY 1000037 2.486895e+00 # wch2 DECAY 2000001 2.858123e-01 # wsd4 DECAY 2000002 1.152973e+00 # wsu4 DECAY 2000003 2.858123e-01 # wsd5 DECAY 2000004 1.152973e+00 # wsu5 DECAY 2000005 8.015663e-01 # wsd6 DECAY 2000006 7.373133e+00 # wsu6 DECAY 2000011 2.161216e-01 # wsl4 DECAY 2000013 2.161216e-01 # wsl5 DECAY 2000015 2.699061e-01 # wsl6 DECAY 1 0.000000e+00 # d : 0.0 DECAY 2 0.000000e+00 # u : 0.0 DECAY 3 0.000000e+00 # s : 0.0 DECAY 4 0.000000e+00 # c : 0.0 DECAY 5 0.000000e+00 # b : 0.0 DECAY 11 0.000000e+00 # e- : 0.0 DECAY 12 0.000000e+00 # ve : 0.0 DECAY 13 0.000000e+00 # mu- : 0.0 DECAY 14 0.000000e+00 # vm : 0.0 DECAY 15 0.000000e+00 # ta- : 0.0 DECAY 16 0.000000e+00 # vt : 0.0 DECAY 21 0.000000e+00 # g : 0.0 DECAY 22 0.000000e+00 # a : 0.0 DECAY 1000022 0.000000e+00 # n1 : 0.0 ################################### ## INFORMATION FOR QNUMBERS 1000022 ################################### BLOCK QNUMBERS 1000022 # n1 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) ################################### ## INFORMATION FOR QNUMBERS 1000023 ################################### BLOCK QNUMBERS 1000023 # n2 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) ################################### ## INFORMATION FOR QNUMBERS 1000025 ################################### BLOCK QNUMBERS 1000025 # n3 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) ################################### ## INFORMATION FOR QNUMBERS 1000035 ################################### BLOCK QNUMBERS 1000035 # n4 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) ################################### ## INFORMATION FOR QNUMBERS 1000024 ################################### BLOCK QNUMBERS 1000024 # x1+ 1 3 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000037 ################################### BLOCK QNUMBERS 1000037 # x2+ 1 3 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000021 ################################### BLOCK QNUMBERS 1000021 # go 1 0 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 35 ################################### BLOCK QNUMBERS 35 # h2 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) ################################### ## INFORMATION FOR QNUMBERS 36 ################################### BLOCK QNUMBERS 36 # h3 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) ################################### ## INFORMATION FOR QNUMBERS 37 ################################### BLOCK QNUMBERS 37 # h+ 1 3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000012 ################################### BLOCK QNUMBERS 1000012 # sve 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 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000014 ################################### BLOCK QNUMBERS 1000014 # svm 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 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000016 ################################### BLOCK QNUMBERS 1000016 # svt 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 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000011 ################################### BLOCK QNUMBERS 1000011 # el- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000013 ################################### BLOCK QNUMBERS 1000013 # mul- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000015 ################################### BLOCK QNUMBERS 1000015 # ta1- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000011 ################################### BLOCK QNUMBERS 2000011 # er- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000013 ################################### BLOCK QNUMBERS 2000013 # mur- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000015 ################################### BLOCK QNUMBERS 2000015 # ta2- 1 -3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000002 ################################### BLOCK QNUMBERS 1000002 # ul 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000004 ################################### BLOCK QNUMBERS 1000004 # cl 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000006 ################################### BLOCK QNUMBERS 1000006 # t1 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000002 ################################### BLOCK QNUMBERS 2000002 # ur 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000004 ################################### BLOCK QNUMBERS 2000004 # cr 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000006 ################################### BLOCK QNUMBERS 2000006 # t2 1 2 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000001 ################################### BLOCK QNUMBERS 1000001 # dl 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000003 ################################### BLOCK QNUMBERS 1000003 # sl 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 1000005 ################################### BLOCK QNUMBERS 1000005 # b1 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000001 ################################### BLOCK QNUMBERS 2000001 # dr 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000003 ################################### BLOCK QNUMBERS 2000003 # sr 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 2000005 ################################### BLOCK QNUMBERS 2000005 # b2 1 -1 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 3 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR MODSEL ################################### BLOCK MODSEL # 1 1 # ################################### ## INFORMATION FOR STOPMIX ################################### BLOCK STOPMIX # 1 1 5.536450e-01 # rru3x3 1 2 8.327528e-01 # rru3x6 2 1 8.327528e-01 # rru6x3 2 2 -5.536450e-01 # rru6x6 ################################### ## INFORMATION FOR SBOTMIX ################################### BLOCK SBOTMIX # 1 1 9.387379e-01 # rrd3x3 1 2 3.446319e-01 # rrd3x6 2 1 -3.446319e-01 # rrd6x3 2 2 9.387379e-01 # rrd6x6 ################################### ## INFORMATION FOR STAUMIX ################################### BLOCK STAUMIX # 1 1 2.824872e-01 # rrl3x3 1 2 9.592711e-01 # rrl3x6 2 1 9.592711e-01 # rrl6x3 2 2 -2.824872e-01 # rrl6x6 ################################### ## INFORMATION FOR ALPHA ################################### BLOCK ALPHA # -1.138252e-01 # alp ################################### ## INFORMATION FOR AE ################################### BLOCK AE Q= 1.000000e+00 # 3 3 -2.517769e+02 # A_tau(Q) DRbar 1 1 0.000000e+00 # A_e(Q) DRbar 2 2 0.000000e+00 # A_mu(Q) DRbar ################################### ## INFORMATION FOR AU ################################### BLOCK AU Q= 1.000000e+00 # 3 3 -4.981299e+02 # A_t(Q) DRbar 1 1 0.000000e+00 # A_u(Q) DRbar 2 2 0.000000e+00 # A_c(Q) DRbar ################################### ## INFORMATION FOR AD ################################### BLOCK AD Q= 1.000000e+00 # 3 3 0.000000e+00 # A_b(Q) DRbar 1 1 0.000000e+00 # A_d(Q) DRbar 2 2 0.000000e+00 # A_s(Q) DRbar """ 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.), RandomizedParameters = cms.VPSet(), ) for i in range(10): mlsp = 100. + i*5.0 mnlsp = 230. + i*10.0 slhatable = baseSLHATable.replace('%MLSP%','%e' % mlsp) slhatable = slhatable.replace('%MNLSP%','%e' % mnlsp) generator.RandomizedParameters.append( cms.PSet( ConfigWeight = cms.double(1.0), ConfigDescription = cms.string('MLSP = %e, MNLSP = %e' % (mlsp,mnlsp)), SLHATableForPythia8 = cms.string('%s' % slhatable), PythiaParameters = basePythiaParameters, ), )