Commit edbe1a80 authored by Franziska Iltzsche's avatar Franziska Iltzsche
Browse files

add 410440


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@811098 4525493e-7705-40b1-a816-d608a930855b
parent 562708ae
2017-10-13 Franziska Iltzsche <franziska.speiser@cern.ch>
* add 410440
* update common/MadGraph/MadGraphControl_ttbar.py
2017-10-12 Franziska Iltzsche <franziska.speiser@cern.ch>
* add 404940 - 404966
* add common/MadGraph/MadGraphControl_SimplifiedModel_RPV_mixedDecay.py
......
......@@ -8,6 +8,7 @@ import os
tt_nonallhad = [410225, 410274, 407045, 407046, 407047, 407048]
tt_dilepton = [410226, 410275]
tt_allhad = [410227,410368]
tt_inclusive = [410440]
b_filter = [410274, 410275]
HT1k5_filter = [407045]
HT1k_1k5_filter = [407046]
......@@ -15,6 +16,7 @@ HT6c_1k_filter = [407047]
MET200_filter = [407048]
pT300_cut = [410368]
enhanced_bwcut = [410440]
thisDSID = runArgs.runNumber
# --------------------------------------------------------------
......@@ -36,8 +38,9 @@ parton_shower='PYTHIA8'
muR_over_ref = 1.0
muF1_over_ref = 1.0
muF2_over_ref = 1.0
dyn_scale = '0' # user-defined scale -> Dominic's definition of mt+1/2*(pt^2+ptx^2)
dyn_scale = '10' # user-defined scale -> Dominic's definition of mt+1/2*(pt^2+ptx^2)
lhe_version=3
bwcut = 15
# --------------------------------------------------------------
# Setting up the process
......@@ -66,9 +69,17 @@ elif thisDSID in tt_dilepton:
evgenConfig.description += ', dilepton'
name = "ttbar_NLO5FS_dilepton"
wdecay = "define lv = e+ mu+ ta+ ve vm vt e- mu- ta- ve~ vm~ vt~ \ndecay t > w+ b, w+ > lv lv \ndecay t~ > w- b~, w- > lv lv \n"
elif thisDSID in tt_inclusive:
evgenConfig.description += ', inclusive'
name = "ttbar_NLO5FS_inclusive"
wdecay = "decay t > w+ b, w+ > all all \ndecay t~ > w- b~, w- > all all \n"
else:
raise RuntimeError("DSID %i is not recoginsed for this control files."%(thisDSID))
if thisDSID in enhanced_bwcut:
bwcut = 50
# Adjust nEvents for the filters
if thisDSID in HT1k5_filter:
......@@ -115,7 +126,7 @@ 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'
mark = ' elseif(dynamical_scale_choice.eq.10) then'
rmLines = ['ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc',
'cc USER-DEFINED SCALE: ENTER YOUR CODE HERE cc',
'cc to use this code you must set cc',
......@@ -187,12 +198,14 @@ mscard.write("""#************************************************************
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set max_weight_ps_point 500 # number of PS to estimate the maximum for each event
set Nevents_for_max_weigth 500
set BW_cut %i
set seed %i
define j = g u c d s b u~ c~ d~ s~ b~
%s
launch
"""%(runArgs.randomSeed, wdecay))
"""%(bwcut, runArgs.randomSeed, wdecay))
mscard.close()
# --------------------------------------------------------------
......@@ -271,9 +284,11 @@ include("MC15JobOptions/Pythia8_aMcAtNlo.py")
# Apply TTbarWToLeptonFilter
# --------------------------------------------------------------
# Only truely needed for non-all hadronic decays
include("MC15JobOptions/TTbarWToLeptonFilter.py")
filtSeq.TTbarWToLeptonFilter.NumLeptons = nLeptons
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
if nLeptons != -99:
include("MC15JobOptions/TTbarWToLeptonFilter.py")
filtSeq.TTbarWToLeptonFilter.NumLeptons = nLeptons
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
# --------------------------------------------------------------
# If required by JO, apply b-filter
......
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