Pythia8Interface/test/randomizedParametersSLHA.py

0001 #example gen fragment where a randomized parameter scan is performed over parameters
0002 #in the SLHA table
0003
0004 #since the randomization occurs only at lumi section boundaries, test with command like
0005 #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)"
0006
0007
0008 import FWCore.ParameterSet.Config as cms
0009
0010 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0011 from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
0012
0013 basePythiaParameters = cms.PSet(
0014   pythia8CommonSettingsBlock,
0015   pythia8CUEP8M1SettingsBlock,
0016   JetMatchingParameters = cms.vstring(
0017     'JetMatching:setMad = off',
0018     'JetMatching:scheme = 1',
0019     'JetMatching:merge = on',
0020     'JetMatching:jetAlgorithm = 2',
0021     'JetMatching:etaJetMax = 5.',
0022     'JetMatching:coneRadius = 1.',
0023     'JetMatching:slowJetPower = 1',
0024     'JetMatching:qCut = 58', #this is the actual merging scale
0025     'JetMatching:nQmatch = 5', #4 corresponds to 4-flavour scheme (no matching of b-quarks), 5 for 5-flavour scheme
0026     'JetMatching:nJetMax = 2', #number of partons in born matrix element for highest multiplicity
0027     'JetMatching:doShowerKt = off', #off for MLM matching, turn on for shower-kT matching
0028   ),
0029   parameterSets = cms.vstring('pythia8CommonSettings',
0030     'pythia8CUEP8M1Settings',
0031     'JetMatchingParameters'
0032   )
0033 )
0034
0035 baseSLHATable="""
0036 ######################################################################
0037 ## PARAM_CARD AUTOMATICALY GENERATED BY MG5                       ####
0038 ######################################################################
0039 ###################################
0040 ## INFORMATION FOR HMIX
0041 ###################################
0042 BLOCK HMIX Q= 1.000000e+00 #
0043       1 3.576810e+02 #  rmuh
0044       2 9.748624e+00 #  tb
0045       4 1.664391e+05 #  ma2
0046       3 2.449353e+02 # higgs vev(Q) MSSM DRb
0047 ###################################
0048 ## INFORMATION FOR MASS
0049 ###################################
0050 BLOCK MASS #
0051       5 0.000000e+00 #  mb
0052       6 1.730000e+02 #  mt
0053       15 1.777000e+00 #  mta
0054       23 9.118760e+01 #  mz
0055       24 80.419002 #  mw
0056       25 1.25e02 #  mh01
0057       35 1.00e05 #  mh02
0058       36 1.00e05 #  ma0
0059       37 1.00e05 #  mh
0060       1000001 1.00e05 #  set of param :1*msd1, 1*msd2
0061       1000002 1.00e05 #  set of param :1*msu1, 1*msu2
0062       1000005 4.70e02 #  msd3
0063       1000006 1.00e05 #  msu3
0064       1000011 1.00e05 #  set of param :1*msl1, 1*msl2
0065       1000012 1.00e05 #  set of param :1*msn1, 1*msn2
0066       1000015 1.00e05 #  msl3
0067       1000016 1.00e05 #  msn3
0068       1000021 1.00e05 #  mgo
0069       1000022 %MLSP% #  mneu1
0070       1000023 %MNLSP% #  mneu2
0071       1000024 1.00e05 #  mch1
0072       1000025 -1.00e05 #  mneu3
0073       1000035 1.00e05 #  mneu4
0074       1000037 1.00e05 #  mch2
0075       2000001 1.00e05 #  set of param :1*msd4, 1*msd5
0076       2000002 1.00e05 #  set of param :1*msu4, 1*msu5
0077       2000005 1.00e05 #  msd6
0078       2000006 1.00e05 #  msu6
0079       2000011 1.00e05 #  set of param :1*msl4, 1*msl5
0080       2000015 1.0005 #  msl6
0081       1 0.000000e00 #  d : 0.0
0082       2 0.000000e00 #  u : 0.0
0083       3 0.000000e00 #  s : 0.0
0084       4 0.000000e00 #  c : 0.0
0085       11 0.000000e00 #  e- : 0.0
0086       12 0.000000e00 #  ve : 0.0
0087       13 0.000000e00 #  mu- : 0.0
0088       14 0.000000e00 #  vm : 0.0
0089       16 0.000000e00 #  vt : 0.0
0090       21 0.000000e00 #  g : 0.0
0091       22 0.000000e00 #  a : 0.0
0092       1000014 1.00e05 #  svm : msn1
0093       1000013 1.00e05 #  mul- : msl1
0094       2000013 1.00e05 #  mur- : msl4
0095       1000004 1.00e05 #  cl : msu1
0096       2000004 1.00e05 #  cr : msu4
0097       1000003 1.00e05 #  sl : msd1
0098       2000003 1.00e05 #  sr : msd4
0099 ###################################
0100 ## INFORMATION FOR MSOFT
0101 ###################################
0102 BLOCK MSOFT Q= 1.000000e+00 #
0103       1 1.013965e+02 #  rmx1
0104       2 1.915042e+02 #  rmx2
0105       3 5.882630e+02 #  rmx3
0106       21 3.233749e+04 #  mhd2
0107       22 -1.288001e+05 #  mhu2
0108       31 1.953348e+02 #  set of param :1*rml21x1, 1*rml22x2
0109       32 1.953348e+02 #  mg5 will not use this value use instead 1*mdl_rml21x1
0110       33 1.944960e+02 #  rml23x3
0111       34 1.364941e+02 #  set of param :1*rme21x1, 1*rme22x2
0112       35 1.364941e+02 #  mg5 will not use this value use instead 1*mdl_rme21x1
0113       36 1.340434e+02 #  rme23x3
0114       41 5.475735e+02 #  set of param :1*rmq21x1, 1*rmq22x2
0115       42 5.475735e+02 #  mg5 will not use this value use instead 1*mdl_rmq21x1
0116       43 4.987639e+02 #  rmq23x3
0117       44 5.295112e+02 #  set of param :1*rmu21x1, 1*rmu22x2
0118       45 5.295112e+02 #  mg5 will not use this value use instead 1*mdl_rmu21x1
0119       46 4.232459e+02 #  rmu23x3
0120       47 5.231488e+02 #  set of param :1*rmd21x1, 1*rmd22x2
0121       48 5.231488e+02 #  mg5 will not use this value use instead 1*mdl_rmd21x1
0122       49 5.198673e+02 #  rmd23x3
0123 ###################################
0124 ## INFORMATION FOR NMIX
0125 ###################################
0126 BLOCK NMIX #
0127       1 1 9.863644e-01 #  rnn1x1
0128       1 2 -5.311036e-02 #  rnn1x2
0129       1 3 1.464340e-01 #  rnn1x3
0130       1 4 -5.311861e-02 #  rnn1x4
0131       2 1 9.935054e-02 #  rnn2x1
0132       2 2 9.449493e-01 #  rnn2x2
0133       2 3 -2.698467e-01 #  rnn2x3
0134       2 4 1.561507e-01 #  rnn2x4
0135       3 1 -6.033880e-02 #  rnn3x1
0136       3 2 8.770049e-02 #  rnn3x2
0137       3 3 6.958775e-01 #  rnn3x3
0138       3 4 7.102270e-01 #  rnn3x4
0139       4 1 -1.165071e-01 #  rnn4x1
0140       4 2 3.107390e-01 #  rnn4x2
0141       4 3 6.492260e-01 #  rnn4x3
0142       4 4 -6.843778e-01 #  rnn4x4
0143 ###################################
0144 ## INFORMATION FOR SMINPUTS
0145 ###################################
0146 BLOCK SMINPUTS #
0147       1 1.279340e+02 #  aewm1
0148       3 1.180000e-01 #  as
0149       6 1.750000e+02 #  mt
0150       7 1.777000e+00 #  mta
0151       4 9.118760e+01 #  mz
0152       2 1.166370e-05 # G_F [GeV^-2]
0153 ###################################
0154 ## INFORMATION FOR UMIX
0155 ###################################
0156 BLOCK UMIX #
0157       1 1 9.168349e-01 #  ruu1x1
0158       1 2 -3.992666e-01 #  ruu1x2
0159       2 1 3.992666e-01 #  ruu2x1
0160       2 2 9.168349e-01 #  ruu2x2
0161 ###################################
0162 ## INFORMATION FOR VMIX
0163 ###################################
0164 BLOCK VMIX #
0165       1 1 9.725578e-01 #  rvv1x1
0166       1 2 -2.326612e-01 #  rvv1x2
0167       2 1 2.326612e-01 #  rvv2x1
0168       2 2 9.725578e-01 #  rvv2x2
0169 ###################################
0170 ## INFORMATION FOR YD
0171 ###################################
0172 BLOCK YD Q= 1.000000e+00 #
0173       3 3 0.000000e+00 #  ryd3x3
0174 ###################################
0175 ## INFORMATION FOR YE
0176 ###################################
0177 BLOCK YE Q= 1.000000e+00 #
0178       3 3 1.008908e-01 #  rye3x3
0179 ###################################
0180 ## INFORMATION FOR YU
0181 ###################################
0182 BLOCK YU Q= 1.000000e+00 #
0183       3 3 8.928445e-01 #  ryu3x3
0184 ###################################
0185 ## INFORMATION FOR DECAY
0186 ###################################
0187 DECAY 6 1.491500E+00 #  wt
0188 DECAY 23 2.441404e+00 #  wz
0189 DECAY 24 2.047600e+00 #  ww
0190 DECAY 25 6.382339e-03 #  wh01
0191 DECAY 35 5.748014e-01 #  wh02
0192 DECAY 36 6.321785e-01 #  wa0
0193 DECAY 37 5.469628e-01 #  wh
0194 DECAY 1000001 5.312788e+00 #  wsd1
0195 DECAY 1000002 5.477195e+00 #  wsu1
0196 DECAY 1000003 5.312788e+00 #  wsd2
0197 DECAY 1000004 5.477195e+00 #  wsu2
0198 DECAY 1000005 1.000000e+00 #  wsd3
0199      1.00000000E+00    2     1000023       5
0200 DECAY 1000006 2.021596e+00 #  wsu3
0201 DECAY 1000011 2.136822e-01 #  wsl1
0202 DECAY 1000012 1.498816e-01 #  wsn1
0203 DECAY 1000013 2.136822e-01 #  wsl2
0204 DECAY 1000014 1.498816e-01 #  wsn2
0205 DECAY 1000015 1.483273e-01 #  wsl3
0206 DECAY 1000016 1.475190e-01 #  wsn3
0207 DECAY 1000021 5.506754e+00 #  wgo
0208 DECAY 1000023 1.000000e-01 #  wneu2
0209      1.00000000E+00    2     1000022       25
0210 DECAY 1000024 1.704145e-02 #  wch1
0211 DECAY 1000025 1.915985e+00 #  wneu3
0212 DECAY 1000035 2.585851e+00 #  wneu4
0213 DECAY 1000037 2.486895e+00 #  wch2
0214 DECAY 2000001 2.858123e-01 #  wsd4
0215 DECAY 2000002 1.152973e+00 #  wsu4
0216 DECAY 2000003 2.858123e-01 #  wsd5
0217 DECAY 2000004 1.152973e+00 #  wsu5
0218 DECAY 2000005 8.015663e-01 #  wsd6
0219 DECAY 2000006 7.373133e+00 #  wsu6
0220 DECAY 2000011 2.161216e-01 #  wsl4
0221 DECAY 2000013 2.161216e-01 #  wsl5
0222 DECAY 2000015 2.699061e-01 #  wsl6
0223 DECAY 1 0.000000e+00 #  d : 0.0
0224 DECAY 2 0.000000e+00 #  u : 0.0
0225 DECAY 3 0.000000e+00 #  s : 0.0
0226 DECAY 4 0.000000e+00 #  c : 0.0
0227 DECAY 5 0.000000e+00 #  b : 0.0
0228 DECAY 11 0.000000e+00 #  e- : 0.0
0229 DECAY 12 0.000000e+00 #  ve : 0.0
0230 DECAY 13 0.000000e+00 #  mu- : 0.0
0231 DECAY 14 0.000000e+00 #  vm : 0.0
0232 DECAY 15 0.000000e+00 #  ta- : 0.0
0233 DECAY 16 0.000000e+00 #  vt : 0.0
0234 DECAY 21 0.000000e+00 #  g : 0.0
0235 DECAY 22 0.000000e+00 #  a : 0.0
0236 DECAY 1000022 0.000000e+00 #  n1 : 0.0
0237 ###################################
0238 ## INFORMATION FOR QNUMBERS 1000022
0239 ###################################
0240 BLOCK QNUMBERS 1000022 #  n1
0241       1 0 #  3 times electric charge
0242       2 2 #  number of spin states (2s+1)
0243       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0244       4 0 #  particle/antiparticle distinction (0=own anti)
0245 ###################################
0246 ## INFORMATION FOR QNUMBERS 1000023
0247 ###################################
0248 BLOCK QNUMBERS 1000023 #  n2
0249       1 0 #  3 times electric charge
0250       2 2 #  number of spin states (2s+1)
0251       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0252       4 0 #  particle/antiparticle distinction (0=own anti)
0253 ###################################
0254 ## INFORMATION FOR QNUMBERS 1000025
0255 ###################################
0256 BLOCK QNUMBERS 1000025 #  n3
0257       1 0 #  3 times electric charge
0258       2 2 #  number of spin states (2s+1)
0259       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0260       4 0 #  particle/antiparticle distinction (0=own anti)
0261 ###################################
0262 ## INFORMATION FOR QNUMBERS 1000035
0263 ###################################
0264 BLOCK QNUMBERS 1000035 #  n4
0265       1 0 #  3 times electric charge
0266       2 2 #  number of spin states (2s+1)
0267       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0268       4 0 #  particle/antiparticle distinction (0=own anti)
0269 ###################################
0270 ## INFORMATION FOR QNUMBERS 1000024
0271 ###################################
0272 BLOCK QNUMBERS 1000024 #  x1+
0273       1 3 #  3 times electric charge
0274       2 2 #  number of spin states (2s+1)
0275       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0276       4 1 #  particle/antiparticle distinction (0=own anti)
0277 ###################################
0278 ## INFORMATION FOR QNUMBERS 1000037
0279 ###################################
0280 BLOCK QNUMBERS 1000037 #  x2+
0281       1 3 #  3 times electric charge
0282       2 2 #  number of spin states (2s+1)
0283       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0284       4 1 #  particle/antiparticle distinction (0=own anti)
0285 ###################################
0286 ## INFORMATION FOR QNUMBERS 1000021
0287 ###################################
0288 BLOCK QNUMBERS 1000021 #  go
0289       1 0 #  3 times electric charge
0290       2 2 #  number of spin states (2s+1)
0291       3 8 #  colour rep (1: singlet, 3: triplet, 8: octet)
0292       4 0 #  particle/antiparticle distinction (0=own anti)
0293 ###################################
0294 ## INFORMATION FOR QNUMBERS 35
0295 ###################################
0296 BLOCK QNUMBERS 35 #  h2
0297       1 0 #  3 times electric charge
0298       2 1 #  number of spin states (2s+1)
0299       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0300       4 0 #  particle/antiparticle distinction (0=own anti)
0301 ###################################
0302 ## INFORMATION FOR QNUMBERS 36
0303 ###################################
0304 BLOCK QNUMBERS 36 #  h3
0305       1 0 #  3 times electric charge
0306       2 1 #  number of spin states (2s+1)
0307       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0308       4 0 #  particle/antiparticle distinction (0=own anti)
0309 ###################################
0310 ## INFORMATION FOR QNUMBERS 37
0311 ###################################
0312 BLOCK QNUMBERS 37 #  h+
0313       1 3 #  3 times electric charge
0314       2 1 #  number of spin states (2s+1)
0315       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0316       4 1 #  particle/antiparticle distinction (0=own anti)
0317 ###################################
0318 ## INFORMATION FOR QNUMBERS 1000012
0319 ###################################
0320 BLOCK QNUMBERS 1000012 #  sve
0321       1 0 #  3 times electric charge
0322       2 1 #  number of spin states (2s+1)
0323       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0324       4 1 #  particle/antiparticle distinction (0=own anti)
0325 ###################################
0326 ## INFORMATION FOR QNUMBERS 1000014
0327 ###################################
0328 BLOCK QNUMBERS 1000014 #  svm
0329       1 0 #  3 times electric charge
0330       2 1 #  number of spin states (2s+1)
0331       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0332       4 1 #  particle/antiparticle distinction (0=own anti)
0333 ###################################
0334 ## INFORMATION FOR QNUMBERS 1000016
0335 ###################################
0336 BLOCK QNUMBERS 1000016 #  svt
0337       1 0 #  3 times electric charge
0338       2 1 #  number of spin states (2s+1)
0339       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0340       4 1 #  particle/antiparticle distinction (0=own anti)
0341 ###################################
0342 ## INFORMATION FOR QNUMBERS 1000011
0343 ###################################
0344 BLOCK QNUMBERS 1000011 #  el-
0345       1 -3 #  3 times electric charge
0346       2 1 #  number of spin states (2s+1)
0347       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0348       4 1 #  particle/antiparticle distinction (0=own anti)
0349 ###################################
0350 ## INFORMATION FOR QNUMBERS 1000013
0351 ###################################
0352 BLOCK QNUMBERS 1000013 #  mul-
0353       1 -3 #  3 times electric charge
0354       2 1 #  number of spin states (2s+1)
0355       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0356       4 1 #  particle/antiparticle distinction (0=own anti)
0357 ###################################
0358 ## INFORMATION FOR QNUMBERS 1000015
0359 ###################################
0360 BLOCK QNUMBERS 1000015 #  ta1-
0361       1 -3 #  3 times electric charge
0362       2 1 #  number of spin states (2s+1)
0363       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0364       4 1 #  particle/antiparticle distinction (0=own anti)
0365 ###################################
0366 ## INFORMATION FOR QNUMBERS 2000011
0367 ###################################
0368 BLOCK QNUMBERS 2000011 #  er-
0369       1 -3 #  3 times electric charge
0370       2 1 #  number of spin states (2s+1)
0371       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0372       4 1 #  particle/antiparticle distinction (0=own anti)
0373 ###################################
0374 ## INFORMATION FOR QNUMBERS 2000013
0375 ###################################
0376 BLOCK QNUMBERS 2000013 #  mur-
0377       1 -3 #  3 times electric charge
0378       2 1 #  number of spin states (2s+1)
0379       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0380       4 1 #  particle/antiparticle distinction (0=own anti)
0381 ###################################
0382 ## INFORMATION FOR QNUMBERS 2000015
0383 ###################################
0384 BLOCK QNUMBERS 2000015 #  ta2-
0385       1 -3 #  3 times electric charge
0386       2 1 #  number of spin states (2s+1)
0387       3 1 #  colour rep (1: singlet, 3: triplet, 8: octet)
0388       4 1 #  particle/antiparticle distinction (0=own anti)
0389 ###################################
0390 ## INFORMATION FOR QNUMBERS 1000002
0391 ###################################
0392 BLOCK QNUMBERS 1000002 #  ul
0393       1 2 #  3 times electric charge
0394       2 1 #  number of spin states (2s+1)
0395       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0396       4 1 #  particle/antiparticle distinction (0=own anti)
0397 ###################################
0398 ## INFORMATION FOR QNUMBERS 1000004
0399 ###################################
0400 BLOCK QNUMBERS 1000004 #  cl
0401       1 2 #  3 times electric charge
0402       2 1 #  number of spin states (2s+1)
0403       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0404       4 1 #  particle/antiparticle distinction (0=own anti)
0405 ###################################
0406 ## INFORMATION FOR QNUMBERS 1000006
0407 ###################################
0408 BLOCK QNUMBERS 1000006 #  t1
0409       1 2 #  3 times electric charge
0410       2 1 #  number of spin states (2s+1)
0411       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0412       4 1 #  particle/antiparticle distinction (0=own anti)
0413 ###################################
0414 ## INFORMATION FOR QNUMBERS 2000002
0415 ###################################
0416 BLOCK QNUMBERS 2000002 #  ur
0417       1 2 #  3 times electric charge
0418       2 1 #  number of spin states (2s+1)
0419       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0420       4 1 #  particle/antiparticle distinction (0=own anti)
0421 ###################################
0422 ## INFORMATION FOR QNUMBERS 2000004
0423 ###################################
0424 BLOCK QNUMBERS 2000004 #  cr
0425       1 2 #  3 times electric charge
0426       2 1 #  number of spin states (2s+1)
0427       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0428       4 1 #  particle/antiparticle distinction (0=own anti)
0429 ###################################
0430 ## INFORMATION FOR QNUMBERS 2000006
0431 ###################################
0432 BLOCK QNUMBERS 2000006 #  t2
0433       1 2 #  3 times electric charge
0434       2 1 #  number of spin states (2s+1)
0435       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0436       4 1 #  particle/antiparticle distinction (0=own anti)
0437 ###################################
0438 ## INFORMATION FOR QNUMBERS 1000001
0439 ###################################
0440 BLOCK QNUMBERS 1000001 #  dl
0441       1 -1 #  3 times electric charge
0442       2 1 #  number of spin states (2s+1)
0443       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0444       4 1 #  particle/antiparticle distinction (0=own anti)
0445 ###################################
0446 ## INFORMATION FOR QNUMBERS 1000003
0447 ###################################
0448 BLOCK QNUMBERS 1000003 #  sl
0449       1 -1 #  3 times electric charge
0450       2 1 #  number of spin states (2s+1)
0451       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0452       4 1 #  particle/antiparticle distinction (0=own anti)
0453 ###################################
0454 ## INFORMATION FOR QNUMBERS 1000005
0455 ###################################
0456 BLOCK QNUMBERS 1000005 #  b1
0457       1 -1 #  3 times electric charge
0458       2 1 #  number of spin states (2s+1)
0459       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0460       4 1 #  particle/antiparticle distinction (0=own anti)
0461 ###################################
0462 ## INFORMATION FOR QNUMBERS 2000001
0463 ###################################
0464 BLOCK QNUMBERS 2000001 #  dr
0465       1 -1 #  3 times electric charge
0466       2 1 #  number of spin states (2s+1)
0467       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0468       4 1 #  particle/antiparticle distinction (0=own anti)
0469 ###################################
0470 ## INFORMATION FOR QNUMBERS 2000003
0471 ###################################
0472 BLOCK QNUMBERS 2000003 #  sr
0473       1 -1 #  3 times electric charge
0474       2 1 #  number of spin states (2s+1)
0475       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0476       4 1 #  particle/antiparticle distinction (0=own anti)
0477 ###################################
0478 ## INFORMATION FOR QNUMBERS 2000005
0479 ###################################
0480 BLOCK QNUMBERS 2000005 #  b2
0481       1 -1 #  3 times electric charge
0482       2 1 #  number of spin states (2s+1)
0483       3 3 #  colour rep (1: singlet, 3: triplet, 8: octet)
0484       4 1 #  particle/antiparticle distinction (0=own anti)
0485 ###################################
0486 ## INFORMATION FOR MODSEL
0487 ###################################
0488 BLOCK MODSEL #
0489       1 1 #
0490 ###################################
0491 ## INFORMATION FOR STOPMIX
0492 ###################################
0493 BLOCK STOPMIX #
0494       1 1 5.536450e-01 #  rru3x3
0495       1 2 8.327528e-01 #  rru3x6
0496       2 1 8.327528e-01 #  rru6x3
0497       2 2 -5.536450e-01 #  rru6x6
0498 ###################################
0499 ## INFORMATION FOR SBOTMIX
0500 ###################################
0501 BLOCK SBOTMIX #
0502       1 1 9.387379e-01 #  rrd3x3
0503       1 2 3.446319e-01 #  rrd3x6
0504       2 1 -3.446319e-01 #  rrd6x3
0505       2 2 9.387379e-01 #  rrd6x6
0506 ###################################
0507 ## INFORMATION FOR STAUMIX
0508 ###################################
0509 BLOCK STAUMIX #
0510       1 1 2.824872e-01 #  rrl3x3
0511       1 2 9.592711e-01 #  rrl3x6
0512       2 1 9.592711e-01 #  rrl6x3
0513       2 2 -2.824872e-01 #  rrl6x6
0514 ###################################
0515 ## INFORMATION FOR ALPHA
0516 ###################################
0517 BLOCK ALPHA #
0518         -1.138252e-01 #  alp
0519 ###################################
0520 ## INFORMATION FOR AE
0521 ###################################
0522 BLOCK AE Q= 1.000000e+00 #
0523       3 3 -2.517769e+02 # A_tau(Q) DRbar
0524       1 1 0.000000e+00 # A_e(Q) DRbar
0525       2 2 0.000000e+00 # A_mu(Q) DRbar
0526 ###################################
0527 ## INFORMATION FOR AU
0528 ###################################
0529 BLOCK AU Q= 1.000000e+00 #
0530       3 3 -4.981299e+02 # A_t(Q) DRbar
0531       1 1 0.000000e+00 # A_u(Q) DRbar
0532       2 2 0.000000e+00 # A_c(Q) DRbar
0533 ###################################
0534 ## INFORMATION FOR AD
0535 ###################################
0536 BLOCK AD Q= 1.000000e+00 #
0537       3 3 0.000000e+00 # A_b(Q) DRbar
0538       1 1 0.000000e+00 # A_d(Q) DRbar
0539       2 2 0.000000e+00 # A_s(Q) DRbar
0540 """
0541
0542 generator = cms.EDFilter("Pythia8HadronizerFilter",
0543   maxEventsToPrint = cms.untracked.int32(1),
0544   pythiaPylistVerbosity = cms.untracked.int32(1),
0545   filterEfficiency = cms.untracked.double(1.0),
0546   pythiaHepMCVerbosity = cms.untracked.bool(False),
0547   comEnergy = cms.double(13000.),
0548   RandomizedParameters = cms.VPSet(),
0549 )
0550
0551 for i in range(10):
0552   mlsp = 100. + i*5.0
0553   mnlsp = 230. + i*10.0
0554
0555   slhatable = baseSLHATable.replace('%MLSP%','%e' % mlsp)
0556   slhatable = slhatable.replace('%MNLSP%','%e' % mnlsp)
0557
0558   generator.RandomizedParameters.append(
0559     cms.PSet(
0560       ConfigWeight = cms.double(1.0),
0561       ConfigDescription = cms.string('MLSP = %e, MNLSP = %e' % (mlsp,mnlsp)),
0562       SLHATableForPythia8 = cms.string('%s' % slhatable),
0563       PythiaParameters = basePythiaParameters,
0564     ),
0565   )