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