Commit a31938dd authored by mcjo's avatar mcjo
Browse files

add 410511-2, 403665-9, 410550,60


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@786712 4525493e-7705-40b1-a816-d608a930855b
parent 15dd3ca5
2016-26-11 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-05-02
* add 410511-2, 403665-9, 410550,60
* update MadGraphControl_tt012FxFx.py
2016-24-11 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-05-01
......
########################################################
# This file is to generate SM tZq events using LO calculations
# There are two processes: trilepton including tau and inclusive decay excluding all Had decays
# contact alhroob@cern.ch
########################################################
from MadGraphControl.MadGraphUtils import *
import fileinput
import shutil
import subprocess
import os
import sys
# General settings
minevents=5000
nevents=5500
# MG merging settings
maxjetflavor=4
ickkw=0
dyn_scale = '0'
### DSID lists (extensions can include systematics samples)
Nominal=[410550,410560]
DSID = runArgs.runNumber
if DSID == 410550:
# the following to generate events with three leptons including tau
mgproc="""generate p p > t b~ j Z $$ W+ W-, (t > w+ b, w+ > l+ vl) , Z > l+ l-\n"""
mgproc +="""add process p p > t~ b j Z $$ W+ W- , (t~ > w- b~, w- > l- vl~) , Z > l+ l-"""
name='tZq'
keyword=['SM','singleTop','tZ','lepton']
elif DSID == 410560:
# the following to generate events with at least one lepton
mgproc="""generate p p > t b~ j Z $$ W+ W- , (t > w+ b , w+ > l+ vl) , Z > all all\n"""
mgproc +="""add process p p > t b~ j Z $$ W+ W- , (t > w+ b , w+ > j j) , Z > l+ l-\n"""
mgproc +="""add process p p > t~ b j Z $$ W+ W- , (t~ > w- b~ , w- > l- vl~) , Z > all all\n"""
mgproc +="""add process p p > t~ b j Z $$ W+ W- , (t~ > w- b~ , w- > j j) , Z > l+ l-"""
name='tZq'
keyword=['SM','singleTop','tZ','lepton']
else:
raise RuntimeError("runNumber %i not recognised in these jobOptions."%runArgs.runNumber)
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(name)
get_param_card = subprocess.Popen(['get_files','-data','aMcAtNlo_param_card_tt.dat'])
get_param_card.wait()
fcard = open('proc_card_mg5.dat','w')
mode=1
cluster_type='condor'
cluster_queue='None'
fcard.write("""
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define q = u c d s b
define q~ = u~ c~ d~ s~ b~
define all = q q~ l- l+ vl vl~
"""+mgproc+"""
output -f
""")
fcard.close()
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
pdflabel="lhapdf"
lhaid=262000
#Fetch default LO run_card.dat and set parameters
extras = { 'lhe_version' : '3.0',
'cut_decays' : 'F',
'pdlabel' : "'"+pdflabel+"'",
'lhaid' : lhaid,
'maxjetflavor' : maxjetflavor,
'asrwgtflavor' : maxjetflavor,
'ickkw' : 0,
'ptj' : 1,
'ptb' : 0,
'mmjj' : 0,
'drjj' : 0.4,
'drll' : 0.4,
'drjl' : 0.4,
'ptl' : 0,
'etal' : -1,
'etab' : -1,
'etaj' : 5,
'dynamical_scale_choice' : dyn_scale
}
scalefact = 1.0
#if runArgs.runNumber in systScalefactUp:
# name=name+'_facsc2_rensc2'
# scalefact = 2.0
#elif runArgs.runNumber in systScalefactDn:
# name=name+'_facsc0p5_rensc0p5'
# scalefact = 0.5
process_dir = new_process()
# Cook the setscales file for the user defined dynamical scale
fileN = process_dir+'/SubProcesses/setscales.f'
mark = ' elseif(dynamical_scale_choice.eq.0) then'
rmLines = ['ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc',
'cc USER-DEFINED SCALE: ENTER YOUR CODE HERE cc',
'cc to use this code you must set cc',
'cc dynamical_scale_choice = 0 cc',
'cc in the run_card (run_card.dat) cc',
'write(*,*) "User-defined scale not set"',
'stop 1',
'temp_scale_id=\'User-defined dynamical scale\' ! use a meaningful string',
'tmp = 0',
'cc USER-DEFINED SCALE: END OF USER CODE cc',
'stop 21',
'rscale = 0'
]
flag=0
for line in fileinput.input(fileN, inplace=1):
toKeep = True
for rmLine in rmLines:
if line.find(rmLine) >= 0:
toKeep = False
break
if toKeep:
print line,
if line.startswith(mark) and flag==0:
flag +=1
print """
c Q^2= 4 * sqrt(mb^2 + pt(b)^2), b quark is the one coming from the gluon splitting. The b quark is always the last in the list, so for tZq at LO(only) and when the decay is done on-the-fly,
c it is always the 6th particle in the list, otherwise we have to change the number or use a loop to find the last particle in the list
rscale=0d0 !factorization scale**2 for pdf1
!factorization scale**2 for pdf2
xmtc=dot(P(0,6),P(0,6))
rscale = 4d0 * sqrt(pt(P(0,6))**2+xmtc)
cc An example how to find the last particle in the list
c do i=3,nexternal
c xmtc=dot(P(0,i),P(0,i))
c if ( xmtc < 30 .and. xmtc > 10 ) then
c rscale = 4d0 * sqrt(pt(P(0,i))**2+xmtc)
c endif
c enddo
"""
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,xqcut=0.,
extras=extras,scalefact=scalefact)
print_cards()
generate(run_card_loc='run_card.dat',
param_card_loc= 'aMcAtNlo_param_card_tt.dat',
mode=mode,
cluster_type=cluster_type,
cluster_queue=cluster_queue,
proc_dir=process_dir)
arrange_output(proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz')
## Shower
evgenConfig.description = 'MadGraph_'+str(name)
evgenConfig.keywords+=keyword
evgenConfig.inputfilecheck = stringy
evgenConfig.minevents = minevents
runArgs.inputGeneratorFile=stringy+'._00001.events.tar.gz'
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[4].split('.')[0])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t1 t1~ $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @0
generate p p > t1 t1~ j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
generate p p > t1 t1~ j j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
'''
njets = 0
evgenLog.info('Registered generation of RPV stop pair production with UDD coupling, stop to two b-quarks; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
evgenConfig.contact = [ "simone.amoroso@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop','RPV']
evgenConfig.description = 'stop direct pair production, RPV decay stop->b+b in simplified model, mstop = '+str(runArgs.jobConfig[0].split('_')[4].split('.')[0])
evgenConfig.minevents = 5000
include ('MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py')
if njets>0:
genSeq.Pythia8.Commands += ["Merging:Process = pp>{t1,1000006}{t1~,-1000006}"]
......@@ -49,7 +49,7 @@ if runArgs.runNumber in inclusive:
msdecay = "decay t > w+ b, w+ > all all \ndecay t~ > w- b~, w- > all all"
elif runArgs.runNumber in nonallhad:
nleptons = -1
nevents = 3.0*runArgs.maxEvents # Filter efficiency
nevents = 3.0*nevents # Filter efficiency
msdecay = "decay t > w+ b, w+ > all all \ndecay t~ > w- b~, w- > all all"
elif runArgs.runNumber in dilepton:
nleptons = 2
......@@ -119,18 +119,23 @@ build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_n
print_cards()
madspin_card_loc='madspin_card.dat'
# Location is dependent on the method
madspin_card_loc = ''
mscard = open(madspin_card_loc,'w')
if not hasattr(runArgs, 'inputGenConfFile'):
madspin_card_loc = 'madspin_card_loc'
fMadSpinCard = open('madspin_card.dat','w')
fMadSpinCard.write('import Events/'+runName+'/events.lhe.gz\n')
fMadSpinCard.write('set ms_dir MadSpin\n')
else:
madspin_card_loc = gridpack_dir+'Cards/madspin_card.dat'
os.unlink(gridpack_dir+'Cards/madspin_card.dat')
fMadSpinCard = open(gridpack_dir+'Cards/madspin_card.dat','w')
fMadSpinCard.write('import '+gridpack_dir+'Events/'+runName+'/events.lhe.gz\n')
fMadSpinCard.write('set ms_dir '+gridpack_dir+'MadSpin\n')
# We cannot use this directory if we wish to use different ME for the decay
# Removing this line will ensure that the MS decay is correct for the process
# and is uncorrelated to the actual gridpack production in this case
# fMadSpinCard.write('set ms_dir '+gridpack_dir+'MadSpin\n')
fMadSpinCard.write('set seed '+str(10000000+int(runArgs.randomSeed))+'\n')
fMadSpinCard.write('''set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
......
######################################################################
## 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.103989E+04 # H (derived)
253 0.102874E+04 # H3p (derived)
255 0.100000E+04 # H5pp (derived)
254 0.102874E+04 # H3z (derived)
256 0.100000E+04 # H5p (derived)
257 0.100000E+04 # H5z (derived)
###################################
## INFORMATION FOR POTENTIALPARAM
###################################
Block potentialparam
1 0.400000E+00 # lam2
2 -.100000E+00 # lam3
3 0.200000E+00 # lam4
4 -.103844E+00 # lam5
5 0.152418E+04 # M1coeff
6 0.540305E+02 # 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.258199E+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 0.100000E+01 # WT
DECAY 23 0.100000E+01 # WZ
DECAY 24 0.100000E+01 # 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)
......@@ -21,12 +21,12 @@ Block mass
1 0.504000E-02 # MD
2 0.255000E-02 # MU
3 0.000000E+00 # MS
4 0.164801E+01 # MC
5 0.491263E+01 # MB
6 0.172000E+03 # MT
4 0.164384E+01 # MC
5 0.457936E+01 # MB
6 0.172500E+03 # MT
11 0.511000E-03 # Me
13 0.105660E+00 # MM
15 0.177700E+01 # MTA
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. ##
......@@ -38,11 +38,11 @@ Block mass
16 0.000000E+00 # vt
21 0.000000E+00 # g
22 0.000000E+00 # a
24 0.803980E+02 # w+ (derived)
252 0.964192E+03 # H (derived)
253 0.969546E+03 # H3p (derived)
24 0.803850E+02 # w+ (derived)
252 0.103989E+04 # H (derived)
253 0.102874E+04 # H3p (derived)
255 0.100000E+04 # H5pp (derived)
254 0.969546E+03 # H3z (derived)
254 0.102874E+04 # H3z (derived)
256 0.100000E+04 # H5p (derived)
257 0.100000E+04 # H5z (derived)
......@@ -50,26 +50,26 @@ Block mass
## INFORMATION FOR POTENTIALPARAM
###################################
Block potentialparam
1 0.000000E+00 # lam2
2 -.120000E+01 # lam3
3 0.120000E+01 # lam4
4 -.345094E+01 # lam5
5 0.300000E+04 # M1coeff
6 0.900000E+03 # M2coeff
1 0.400000E+00 # lam2
2 -.100000E+00 # lam3
3 0.200000E+00 # lam4
4 -.103844E+00 # lam5
5 0.152418E+04 # M1coeff
6 0.540305E+02 # M2coeff
###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
1 0.132341E+03 # aEWM1
2 0.116637E-04 # Gf
3 0.118400E+00 # aS
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
1 0.258199E+00 # tanth
###################################
## INFORMATION FOR YUKAWA
......@@ -78,26 +78,26 @@ Block yukawa
1 0.504000E-02 # ymdo
2 0.255000E-02 # ymup
3 0.000000E+00 # yms
4 0.164801E+01 # ymc
5 0.491263E+01 # ymb
6 0.172000E+03 # ymt
4 0.164384E+01 # ymc
5 0.457936E+01 # ymb
6 0.172500E+03 # ymt
11 0.511000E-03 # yme
13 0.105660E+00 # ymm
15 0.177700E+01 # ymtau
13 0.105658E+00 # ymm
15 0.177682E+01 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.51013490E+00 # WT
DECAY 23 2.44639985E+00 # WZ
DECAY 24 2.03535570E+00 # WW
DECAY 25 4.27608700E-03 # 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
DECAY 6 0.100000E+01 # WT
DECAY 23 0.100000E+01 # WZ
DECAY 24 0.100000E+01 # 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 ##
......
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_directTT4B.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_directTT4B.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_directTT4B.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_directTT4B.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_directTT4B.py' )
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 3.0*top mass, A14 tune var 3c up, scale=0.5, at least one lepton filter'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'james.robinson@cern.ch','andrea.helen.knue@cern.ch','onofrio@liverpool.ac.uk','ian.connelly@cern.ch' ]
include('PowhegControl/PowhegControl_tt_Common.py')
PowhegConfig.topdecaymode = 22222
PowhegConfig.hdamp = 517.5
PowhegConfig.PDF = 260000
PowhegConfig.mu_F = 0.5000
PowhegConfig.mu_R = 0.5000
#Information on how to run with multiple weights: https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/PowhegForATLAS#Running_with_multiple_scale_PDF
#PDFs - you can see a listing here: https://lhapdf.hepforge.org/pdfsets.html; picked these three as they are the inputs to the PDF4LHC2015 prescription (http://arxiv.org/pdf/1510.03865v2.pdf).
# Weight group configuration
PowhegConfig.define_event_weight_group( group_name='scales', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, NNPDF
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_NNPDF', parameter_values=[ 2.0, 1.0, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_NNPDF', parameter_values=[ 0.5, 1.0, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muR_NNPDF', parameter_values=[ 1.0, 2.0, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muR_NNPDF', parameter_values=[ 1.0, 0.5, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_0p5muR_NNPDF', parameter_values=[ 0.5, 0.5, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_2muR_NNPDF', parameter_values=[ 2.0, 2.0, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_2muR_NNPDF', parameter_values=[ 0.5, 2.0, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_0p5muR_NNPDF', parameter_values=[ 2.0, 0.5, 260000] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='1muF_1muR_NNPDF', parameter_values=[ 1.0, 1.0, 260000] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='1muF_1muR_MMHT', parameter_values=[ 1.0, 1.0, 25200] )
# Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales', weight_name='1muF_1muR_CT14', parameter_values=[ 1.0, 1.0, 13165] )