Commit cb3efcaf authored by Cyril Pascal Becot's avatar Cyril Pascal Becot
Browse files

See changelog


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@811174 4525493e-7705-40b1-a816-d608a930855b
parent 1f7caf78
2017-10-17 Cyril Becot <cyril.becot@cern.ch>
* add MadGraphControl_GM_VBS_add_lvll.py
* up MadGraph_param_card_GMm250.dat MadGraph_param_card_GMm350.dat
* add MadGraph_param_card_GMm425.dat MadGraph_param_card_GMm450.dat
* add 309497,309498,309500,309501
* up 410392,410393
* tag MC15JobOptions-00-08-06
2017-10-16 Franziska Iltzsche <franziska.speiser@cern.ch>
* tagging MC15JobOptions-00-08-05
* update common/MadGraph/MadGraphControl_SimplifiedModel_RPV_mixedDecay.py
......
from MadGraphControl.MadGraphUtils import *
nevents=40000
mode=0
### DSID lists
DSIDDic = {309500:'250',309501:'350',309497:'425',309498:'450'}
mass = DSIDDic[runArgs.runNumber]
### Cards
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm
define vl~ = ve~ vm~
import model GM_UFO
generate p p > H5p > w+ z j j QCD=0, w+ > l+ vl, z > l+ l- @1
add process p p > H5p~ > w- z j j QCD=0, w- > l- vl~, z > l+ l- @2
add process p p > H5p > w+ z j j QCD=0, w+ > ta+ vt , z > l+ l- @3
add process p p > H5p~ > w- z j j QCD=0, w- > ta- vt~, z > l+ l- @4
output -f
""")
fcard.close()
extras = {
'pdlabel':"'lhapdf'",
'lhaid':"247000",
'ptj':"15",
'pta':"0",
'ptb':"0",
'etal':"2.7",
'drjj':"0",
'drll':"0",
'draa':"0",
'draj':"0",
'drjl':"0",
'dral':"0",
'sys_pdf':"NNPDF23_lo_as_0130_qed.LHgrid" }
build_run_card(run_card_old='MadGraph_run_card_GM.dat',run_card_new='run_card.dat',nevts=2*nevents,rand_seed=runArgs.randomSeed,beamEnergy=runArgs.ecmEnergy/2.,extras=extras)
print_cards()
runName='GM_VBS_H5p_lvll_'+mass+'_ISHdot5_qcd0'
process_dir = new_process(card_loc='proc_card_mg5.dat')
evgenConfig.generators = ["MadGraph"]
generate(run_card_loc='run_card.dat',param_card_loc='MadGraph_param_card_GMm'+mass+'.dat',mode=mode,proc_dir=process_dir,run_name=runName)
############################
# Shower JOs will go here
import os
if 'ATHENA_PROC_NUMBER' in os.environ:
njobs = os.environ.pop('ATHENA_PROC_NUMBER')
# Try to modify the opts underfoot
if not hasattr(opts,'nprocs'): print 'Warning: Did not see option!'
else: opts.nprocs = 0
print opts
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
include("MC15JobOptions/nonStandard/Pythia8_A14_NNPDF23LO_Common.py")
include("MC15JobOptions/Pythia8_LHEF.py")
include("MC15JobOptions/Pythia8_EvtGen.py")
evgenConfig.generators += ["MadGraph", "Pythia8","EvtGen"]
evgenConfig.contact = ["Benjamin Freund <Benjamin.Freund@cern.ch>"]
genSeq.Pythia8.Commands += [
"Init:showAllParticleData = on"]
evgenConfig.description = 'MadGraph_GM_VBS_H5p'
evgenConfig.keywords+=['BSM','VBS','Higgs','WZ','leptonic']
#stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_EffDM'
#runArgs.inputGeneratorFile=stringy+'._00001.events.tar.gz'
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
include('MC15JobOptions/MultiLeptonFilter.py')
MultiLeptonFilter = filtSeq.MultiLeptonFilter
MultiLeptonFilter.Ptcut = 3500.
MultiLeptonFilter.Etacut = 2.7
MultiLeptonFilter.NLeptons = 3
......@@ -88,9 +88,9 @@ Block yukawa
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 0.100000E+01 # WT
DECAY 23 0.100000E+01 # WZ
DECAY 24 0.100000E+01 # WW
DECAY 6 Auto # WT
DECAY 23 Auto # WZ
DECAY 24 Auto # WW
DECAY 25 Auto # Wh
DECAY 252 Auto # Wh__2
DECAY 253 Auto # WH3p
......
......@@ -88,9 +88,9 @@ Block yukawa
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 0.100000E+01 # WT
DECAY 23 0.100000E+01 # WZ
DECAY 24 0.100000E+01 # WW
DECAY 6 Auto # WT
DECAY 23 Auto # WZ
DECAY 24 Auto # WW
DECAY 25 Auto # Wh
DECAY 252 Auto # Wh__2
DECAY 253 Auto # WH3p
......
######################################################################
## PARAM_CARD GENERATED BY GMCALC ##
######################################################################
## ##
## Width set on Auto will be computed following the ##
## information present in the decay.py files of the ##
## model. By default, this is only 1->2 decay modes. ##
## ##
######################################################################
###################################
## INFORMATION FOR CKMBLOCK
###################################
Block ckmblock
1 0.227736E+00 # cabi
###################################
## INFORMATION FOR MASS
###################################
Block mass
1 0.504000E-02 # MD
2 0.255000E-02 # MU
3 0.000000E+00 # MS
4 0.164384E+01 # MC
5 0.457936E+01 # MB
6 0.172500E+03 # MT
11 0.511000E-03 # Me
13 0.105658E+00 # MM
15 0.177682E+01 # MTA
23 0.911876E+02 # MZ
25 0.125000E+03 # Mh
## Dependent parameters, given by model restrictions. ##
## MG5 ignores these values but they are important ##
## for interfacing the output of MG5 to external ##
## programs such as Pythia. ##
12 0.000000E+00 # ve
14 0.000000E+00 # vm
16 0.000000E+00 # vt
21 0.000000E+00 # g
22 0.000000E+00 # a
24 0.803850E+02 # w+ (derived)
252 0.462693E+03 # H (derived)
253 0.463722E+03 # H3p (derived)
255 0.425000E+03 # H5pp (derived)
254 0.463722E+03 # H3z (derived)
256 0.425000E+03 # H5p (derived)
257 0.425000E+03 # H5z (derived)
###################################
## INFORMATION FOR POTENTIALPARAM
###################################
Block potentialparam
1 0.170000E+00 # lam2
2 -.100000E+00 # lam3
3 0.200000E+00 # lam4
4 -.868238E+00 # lam5
5 0.693140E+03 # M1coeff
6 0.115523E+03 # M2coeff
###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
1 0.132233E+03 # aEWM1
2 0.116638E-04 # Gf
3 0.118000E+00 # aS
###################################
## INFORMATION FOR VEV
###################################
Block vev
1 0.577350E+00 # tanth
###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
1 0.504000E-02 # ymdo
2 0.255000E-02 # ymup
3 0.000000E+00 # yms
4 0.164384E+01 # ymc
5 0.457936E+01 # ymb
6 0.172500E+03 # ymt
11 0.511000E-03 # yme
13 0.105658E+00 # ymm
15 0.177682E+01 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 Auto # WT
DECAY 23 Auto # WZ
DECAY 24 Auto # WW
DECAY 25 Auto # Wh
DECAY 252 Auto # Wh__2
DECAY 253 Auto # WH3p
DECAY 254 Auto # WH3z
DECAY 255 Auto # WH5pp
DECAY 256 Auto # WH5p
DECAY 257 Auto # WH5z
## Dependent parameters, given by model restrictions. ##
## MG5 ignores these values but they are important ##
## for interfacing the output of MG5 to external ##
## programs such as Pythia. ##
DECAY 1 0.000000E+00 # d
DECAY 2 0.000000E+00 # u
DECAY 3 0.000000E+00 # s
DECAY 4 0.000000E+00 # c
DECAY 5 0.000000E+00 # b
DECAY 11 0.000000E+00 # e-
DECAY 12 0.000000E+00 # ve
DECAY 13 0.000000E+00 # mu-
DECAY 14 0.000000E+00 # vm
DECAY 15 0.000000E+00 # ta-
DECAY 16 0.000000E+00 # vt
DECAY 21 0.000000E+00 # g
DECAY 22 0.000000E+00 # a
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================
Block QNUMBERS 252 # H
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 QNUMBERS 253 # H3p
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)
Block QNUMBERS 254 # H3z
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 QNUMBERS 255 # H5++
1 6 # 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)
Block QNUMBERS 256 # H5p
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)
Block QNUMBERS 257 # H5z
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)
######################################################################
## PARAM_CARD GENERATED BY GMCALC ##
######################################################################
## ##
## Width set on Auto will be computed following the ##
## information present in the decay.py files of the ##
## model. By default, this is only 1->2 decay modes. ##
## ##
######################################################################
###################################
## INFORMATION FOR CKMBLOCK
###################################
Block ckmblock
1 0.227736E+00 # cabi
###################################
## INFORMATION FOR MASS
###################################
Block mass
1 0.504000E-02 # MD
2 0.255000E-02 # MU
3 0.000000E+00 # MS
4 0.164384E+01 # MC
5 0.457936E+01 # MB
6 0.172500E+03 # MT
11 0.511000E-03 # Me
13 0.105658E+00 # MM
15 0.177682E+01 # MTA
23 0.911876E+02 # MZ
25 0.125000E+03 # Mh
## Dependent parameters, given by model restrictions. ##
## MG5 ignores these values but they are important ##
## for interfacing the output of MG5 to external ##
## programs such as Pythia. ##
12 0.000000E+00 # ve
14 0.000000E+00 # vm
16 0.000000E+00 # vt
21 0.000000E+00 # g
22 0.000000E+00 # a
24 0.803850E+02 # w+ (derived)
252 0.485960E+03 # H (derived)
253 0.486748E+03 # H3p (derived)
255 0.450000E+03 # H5pp (derived)
254 0.486748E+03 # H3z (derived)
256 0.450000E+03 # H5p (derived)
257 0.450000E+03 # H5z (derived)
###################################
## INFORMATION FOR POTENTIALPARAM
###################################
Block potentialparam
1 0.180000E+00 # lam2
2 -.100000E+00 # lam3
3 0.200000E+00 # lam4
4 -.868525E+00 # lam5
5 0.756018E+03 # M1coeff
6 0.126003E+03 # M2coeff
###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
1 0.132233E+03 # aEWM1
2 0.116638E-04 # Gf
3 0.118000E+00 # aS
###################################
## INFORMATION FOR VEV
###################################
Block vev
1 0.577350E+00 # tanth
###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
1 0.504000E-02 # ymdo
2 0.255000E-02 # ymup
3 0.000000E+00 # yms
4 0.164384E+01 # ymc
5 0.457936E+01 # ymb
6 0.172500E+03 # ymt
11 0.511000E-03 # yme
13 0.105658E+00 # ymm
15 0.177682E+01 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 Auto # WT
DECAY 23 Auto # WZ
DECAY 24 Auto # WW
DECAY 25 Auto # Wh
DECAY 252 Auto # Wh__2
DECAY 253 Auto # WH3p
DECAY 254 Auto # WH3z
DECAY 255 Auto # WH5pp
DECAY 256 Auto # WH5p
DECAY 257 Auto # WH5z
## Dependent parameters, given by model restrictions. ##
## MG5 ignores these values but they are important ##
## for interfacing the output of MG5 to external ##
## programs such as Pythia. ##
DECAY 1 0.000000E+00 # d
DECAY 2 0.000000E+00 # u
DECAY 3 0.000000E+00 # s
DECAY 4 0.000000E+00 # c
DECAY 5 0.000000E+00 # b
DECAY 11 0.000000E+00 # e-
DECAY 12 0.000000E+00 # ve
DECAY 13 0.000000E+00 # mu-
DECAY 14 0.000000E+00 # vm
DECAY 15 0.000000E+00 # ta-
DECAY 16 0.000000E+00 # vt
DECAY 21 0.000000E+00 # g
DECAY 22 0.000000E+00 # a
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================
Block QNUMBERS 252 # H
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 QNUMBERS 253 # H3p
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)
Block QNUMBERS 254 # H3z
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 QNUMBERS 255 # H5++
1 6 # 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)
Block QNUMBERS 256 # H5p
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)
Block QNUMBERS 257 # H5z
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)
include('MC15JobOptions/MadGraphControl_GM_VBS_add_lvll.py')
evgenConfig.description = 'MadGraph_GM_VBS_H5p'
evgenConfig.keywords+=['BSM','VBS','Higgs','WZ','leptonic']
include('MC15JobOptions/MadGraphControl_GM_VBS_add_lvll.py')
evgenConfig.description = 'MadGraph_GM_VBS_H5p'
evgenConfig.keywords+=['BSM','VBS','Higgs','WZ','leptonic']
include('MC15JobOptions/MadGraphControl_GM_VBS_add_lvll.py')
evgenConfig.description = 'MadGraph_GM_VBS_H5p'
evgenConfig.keywords+=['BSM','VBS','Higgs','WZ','leptonic']
include('MC15JobOptions/MadGraphControl_GM_VBS_add_lvll.py')
evgenConfig.description = 'MadGraph_GM_VBS_H5p'
evgenConfig.keywords+=['BSM','VBS','Higgs','WZ','leptonic']
......@@ -5,14 +5,12 @@ evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equ
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'sonja.verena.bartkowski@cern.ch' ]
include('PowhegControl/PowhegControl_tt_Common.py')
PowhegConfig.decay_mode = "t t~ > none"
PowhegConfig.MadSpin_decays= ["decay t > w+ b, w+ > l+ vl", "decay t~ > w- s~, w- > l- vl~","decay t~ > w- d~, w- > l- vl~","decay t > w+ b, w+ > j j", "decay t~ > w- s~, w- > j j","decay t~ > w- d~, w- > j j"]
PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]\n add process b g > t t~ b [QCD]\n add process b~ g > t t~ b~ [QCD]\n add process b b~ > t t~[QCD]"
PowhegConfig.hdamp = 172.5
PowhegConfig.hdamp = 258.75
# compensate filter efficiency
PowhegConfig.nEvents *= 3.
......
......@@ -11,7 +11,7 @@ PowhegConfig.decay_mode = "t t~ > none"
PowhegConfig.MadSpin_decays= ["decay t > w+ s, w+ > l+ vl","decay t > w+ d, w+ > l+ vl","decay t~ > w- b~, w- > l- vl~","decay t > w+ s, w+ > j j","decay t > w+ d, w+ > j j","decay t~ > w- b~, w- > j j"]
PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]\n add process b g > t t~ b [QCD]\n add process b~ g > t t~ b~ [QCD] \n add process b b~ > t t~[QCD]"
PowhegConfig.hdamp = 172.5
PowhegConfig.hdamp = 258.75
# compensate filter efficiency
PowhegConfig.nEvents *= 3.
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment