Commit 13ea5518 authored by Joany Manjarres's avatar Joany Manjarres
Browse files

see ChangeLog

parent c8da2957
2019-12-06 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-02-21
* update MadGraphControl_tHjb_CP_NLO.py MadGraphControl_tWH_CP_NLO.py
* add JO 346768-775
2019-12-06 Joany Manjarres <joany@cern.ch> 2019-12-06 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-02-20 * tag MC15JobOptions-01-02-20
* add param_SMlimit_CR_2p5.dat param_SMlimit_CR_4p9.dat param_SMlimit_CR_6.dat param_SMlimit_CR_m2p5.dat param_SMlimit_CR_m4p9.dat param_SMlimit_CR_m6.dat * add param_SMlimit_CR_2p5.dat param_SMlimit_CR_4p9.dat param_SMlimit_CR_6.dat param_SMlimit_CR_m2p5.dat param_SMlimit_CR_m4p9.dat param_SMlimit_CR_m6.dat
......
...@@ -14,7 +14,7 @@ mode=2 ...@@ -14,7 +14,7 @@ mode=2
njobs=132 njobs=132
runName='run_01' runName='run_01'
tx0jb_4fl = [346533, 346534, 346535, 346536, 346537, 346538, 346539, 346540, 346541, 346542, 346543, 346544, 346676, 346677, 346764, 346765, 346766, 346767] tx0jb_4fl = [346533, 346534, 346535, 346536, 346537, 346538, 346539, 346540, 346541, 346542, 346543, 346544, 346676, 346677, 346764, 346765, 346766, 346767, 346772, 346773, 346774, 346775]
tx0jb_4fl_CPalpha_0 = [346676, 346677] tx0jb_4fl_CPalpha_0 = [346676, 346677]
tx0jb_4fl_CPalpha_15 = [346533, 346539] tx0jb_4fl_CPalpha_15 = [346533, 346539]
...@@ -24,10 +24,10 @@ tx0jb_4fl_CPalpha_60 = [346536, 346542] ...@@ -24,10 +24,10 @@ tx0jb_4fl_CPalpha_60 = [346536, 346542]
tx0jb_4fl_CPalpha_75 = [346537, 346543] tx0jb_4fl_CPalpha_75 = [346537, 346543]
tx0jb_4fl_CPalpha_90 = [346538, 346544] tx0jb_4fl_CPalpha_90 = [346538, 346544]
tx0jb_yt_minus1_CPalpha_0 = [346764] tx0jb_yt_minus1_CPalpha_0 = [346764, 346772]
tx0jb_yt_plus0p5_CPalpha_0 = [346765] tx0jb_yt_plus0p5_CPalpha_0 = [346765, 346773]
tx0jb_yt_plus2_CPalpha_0 = [346766] tx0jb_yt_plus2_CPalpha_0 = [346766, 346774]
tx0jb_yt_plus2_CPalpha_45 = [346767] tx0jb_yt_plus2_CPalpha_45 = [346767, 346775]
if runArgs.runNumber in tx0jb_4fl: if runArgs.runNumber in tx0jb_4fl:
......
...@@ -16,10 +16,10 @@ cluster_queue='8nh' ...@@ -16,10 +16,10 @@ cluster_queue='8nh'
## map DSID to process settings ## map DSID to process settings
# select diagram removal scheme: DR1 or DR2? # select diagram removal scheme: DR1 or DR2?
tWH_DR1 = [] tWH_DR1 = []
tWH_DR2 = [ 346486, 346511, 346545, 346546, 346547, 346548, 346549, 346550, 346589, 346590, 346591, 346592, 346593, 346594, 346678, 346759, 346760, 346761, 346762, 346763, 346883] tWH_DR2 = [ 346486, 346511, 346545, 346546, 346547, 346548, 346549, 346550, 346589, 346590, 346591, 346592, 346593, 346594, 346678, 346759, 346760, 346761, 346762, 346763, 346768, 346769, 346770, 346771 ]
# select any BSM top Yukawa couplings (default: SM): # select any BSM top Yukawa couplings (default: SM):
tWH_CPalpha_0 = [ 346678, 346759, 346883] tWH_CPalpha_0 = [ 346678, 346759]
tWH_CPalpha_15 = [ 346545, 346589] tWH_CPalpha_15 = [ 346545, 346589]
tWH_CPalpha_30 = [ 346546, 346590] tWH_CPalpha_30 = [ 346546, 346590]
tWH_CPalpha_45 = [ 346547, 346591] tWH_CPalpha_45 = [ 346547, 346591]
...@@ -30,20 +30,10 @@ tWH_CPalpha_90 = [ 346550, 346594] ...@@ -30,20 +30,10 @@ tWH_CPalpha_90 = [ 346550, 346594]
tWH_yt_minus1 = [] tWH_yt_minus1 = []
tWH_yt_plus2 = [] tWH_yt_plus2 = []
tWH_yt_minus1_CPalpha_0 = [346760] tWH_yt_minus1_CPalpha_0 = [346760, 346768]
tWH_yt_plus0p5_CPalpha_0 = [346761] tWH_yt_plus0p5_CPalpha_0 = [346761, 346769]
tWH_yt_plus2_CPalpha_0 = [346762] tWH_yt_plus2_CPalpha_0 = [346762, 346770]
tWH_yt_plus2_CPalpha_45 = [346763] tWH_yt_plus2_CPalpha_45 = [346763, 346771]
# set shower subtraction terms: defaul: Pythia8, this array H7
tWH_H7 = [346883]
# in case of H7, we also set the Higgs decay channel to H7...
# diphotons:
tWH_H7_yy = [346883]
parton_shower='PYTHIA8'
if runArgs.runNumber in tWH_H7:
parton_shower='HERWIGPP'
tWH_BSM = tWH_yt_minus1+tWH_yt_plus2+tWH_CPalpha_0+tWH_CPalpha_15+tWH_CPalpha_30+tWH_CPalpha_45+tWH_CPalpha_60+tWH_CPalpha_75+tWH_CPalpha_90+tWH_yt_minus1_CPalpha_0+tWH_yt_plus0p5_CPalpha_0+tWH_yt_plus2_CPalpha_0+tWH_yt_plus2_CPalpha_45 tWH_BSM = tWH_yt_minus1+tWH_yt_plus2+tWH_CPalpha_0+tWH_CPalpha_15+tWH_CPalpha_30+tWH_CPalpha_45+tWH_CPalpha_60+tWH_CPalpha_75+tWH_CPalpha_90+tWH_yt_minus1_CPalpha_0+tWH_yt_plus0p5_CPalpha_0+tWH_yt_plus2_CPalpha_0+tWH_yt_plus2_CPalpha_45
...@@ -81,7 +71,7 @@ fcard.close() ...@@ -81,7 +71,7 @@ fcard.close()
extras = {'pdlabel' : "'lhapdf'", extras = {'pdlabel' : "'lhapdf'",
'lhaid' : 260000, 'lhaid' : 260000,
'parton_shower' : parton_shower, 'parton_shower' :'PYTHIA8',
'reweight_scale' : 'True', 'reweight_scale' : 'True',
'reweight_PDF' : 'True', 'reweight_PDF' : 'True',
'PDF_set_min' : 260001, 'PDF_set_min' : 260001,
...@@ -251,41 +241,17 @@ outputDS=arrange_output(run_name=runName,proc_dir=process_dir,lhe_version=3,save ...@@ -251,41 +241,17 @@ outputDS=arrange_output(run_name=runName,proc_dir=process_dir,lhe_version=3,save
evgenConfig.description = 'aMcAtNlo tHW' evgenConfig.description = 'aMcAtNlo tHW'
evgenConfig.keywords+=['Higgs', 'tHiggs'] evgenConfig.keywords+=['Higgs', 'tHiggs']
evgenConfig.contact = ['maria.moreno.llacer@cern.ch'] evgenConfig.contact = ['maria.moreno.llacer@cern.ch']
# running on the fly, name depends on MGUtils, so no check for particular pattern here: evgenConfig.inputfilecheck = outputDS
evgenConfig.inputfilecheck = ''
runArgs.inputGeneratorFile=outputDS runArgs.inputGeneratorFile=outputDS
###### ######
## shower settings: ## shower settings: if you change these to another shower eg Hw++,
## make sure the 'parton_shower' argument to MG is consistent with shower you use ## make sure you update subtraction term ('parton_shower' :'PYTHIA8')
if not runArgs.runNumber in tWH_H7: include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py") include("MC15JobOptions/Pythia8_aMcAtNlo.py")
include("MC15JobOptions/Pythia8_aMcAtNlo.py") include("MC15JobOptions/Pythia8_SMHiggs125_inc.py")
include("MC15JobOptions/Pythia8_SMHiggs125_inc.py") ## don't include shower weights, see ATLMCPROD-6135
## don't include shower weights, see ATLMCPROD-6135 #include("MC15JobOptions/Pythia8_ShowerWeights.py")
#include("MC15JobOptions/Pythia8_ShowerWeights.py")
elif runArgs.runNumber in tWH_H7:
evgenConfig.generators += ["aMcAtNlo", "Herwig7"]
evgenConfig.description = 'MG5_aMC@NLO+Herwig7+EvtGen '+name+' OTF, H7p1 default tune, ME NNPDF 3.0 NLO'
evgenConfig.tune = "H7.1-Default"
evgenConfig.contact += ['liza.mijovic@cern.ch']
include("MC15JobOptions/Herwig7_LHEF.py")
Herwig7Config.me_pdf_commands(order="NLO", name="NNPDF30_nlo_as_0118")
Herwig7Config.tune_commands()
Herwig7Config.lhef_mg5amc_commands(lhe_filename=runArgs.inputGeneratorFile, me_pdf_order="NLO")
include("MC15JobOptions/Herwig71_EvtGen.py")
if runArgs.runNumber in tWH_H7_yy:
# only consider H->gamgam decays
Herwig7Config.add_commands("""
# force H->gamgam decays
do /Herwig/Particles/h0:SelectDecayModes h0->gamma,gamma;
# print out Higgs decays modes and branching ratios to log.generate
do /Herwig/Particles/h0:PrintDecayModes
""")
# common to all tWH_H7
Herwig7Config.run()
else:
theApp.finalize()
theApp.exit()
include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')
include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')
include('MC15JobOptions/MadGraphControl_tWH_CP_NLO.py')
include('MC15JobOptions/MadGraphControl_tHjb_CP_NLO.py')
evgenConfig.inputconfcheck = "tHjb"
include('MC15JobOptions/MadGraphControl_tHjb_CP_NLO.py')
evgenConfig.inputconfcheck = "tHjb"
include('MC15JobOptions/MadGraphControl_tHjb_CP_NLO.py')
evgenConfig.inputconfcheck = "tHjb"
include('MC15JobOptions/MadGraphControl_tHjb_CP_NLO.py')
evgenConfig.inputconfcheck = "tHjb"
MC15JobOptions-01-02-20/ MC15JobOptions-01-02-21/
\ No newline at end of file \ No newline at end of file
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