Commit d86fdd72 authored by Giovanni Siragusa's avatar Giovanni Siragusa
Browse files

Add 410274-410275, Update 410503, 410505


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@780497 4525493e-7705-40b1-a816-d608a930855b
parent b6b5cc6f
2016-10-26 Giovanni Siragusa <siragusa@cern.ch>
* Update common/MadGraph/MadGraphControl_ttbar.py
* Update 410503, 410505
* Add 410274-410275
2016-10-24 Christian Gutschow <chris.g@cern.ch>
* Add DiLeptonMassSliceFilter.py
* tagging MC15JobOptions-00-04-74
......
......@@ -5,9 +5,10 @@ import subprocess
import os
# DSID list of ttbar samples
tt_nonallhad = [410225]
tt_dilepton = [410226]
tt_nonallhad = [410225, 410274]
tt_dilepton = [410226, 410275]
tt_allhad = [410227]
b_filter = [410274, 410275]
thisDSID = runArgs.runNumber
# --------------------------------------------------------------
......@@ -62,6 +63,10 @@ elif thisDSID in tt_dilepton:
else:
raise RuntimeError("DSID %i is not recoginsed for this control files."%(thisDSID))
# B-filtering
if thisDSID in b_filter:
nevents = nevents*25
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(name)
# Proc card writing
......@@ -218,7 +223,13 @@ 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
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
# --------------------------------------------------------------
# If required by JO, apply b-filter
# --------------------------------------------------------------
if thisDSID in b_filter:
include('MC15JobOptions/TTbarPlusHFFilter.py')
filtSeq.TTbarPlusBFilter.SelectB = True
......@@ -2,12 +2,12 @@
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, dilepton filter, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', '2lepton']
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')
# Initial settings
PowhegConfig.topdecaymode = 22222 # Inclusive
PowhegConfig.topdecaymode = 22200 # Dilepton
PowhegConfig.hdamp = 258.75 # 1.5 * mtop
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
......@@ -52,7 +52,7 @@ PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2mu
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 90900] )
PowhegConfig.nEvents *= 25. # compensate filter efficiency
PowhegConfig.nEvents *= 1.1 # Safety factor
PowhegConfig.generate()
#--------------------------------------------------------------
......
......@@ -7,7 +7,7 @@ evgenConfig.contact = [ 'james.robinson@cern.ch','andrea.helen.knue@cern.ch'
include('PowhegControl/PowhegControl_tt_Common.py')
# Initial settings
PowhegConfig.topdecaymode = 22222 # Inclusive
PowhegConfig.topdecaymode = 22200 # Dilepton
PowhegConfig.hdamp = 258.75 # 1.5 * mtop
PowhegConfig.mu_F = [1.0, 2.0, 0.5, 1.0, 1.0, 0.5, 2.0, 0.5, 2.0] # List of factorisation scales which pairs with renormalisation scale below
PowhegConfig.mu_R = [1.0, 1.0, 1.0, 2.0, 0.5, 0.5, 2.0, 2.0, 0.5] # List of renormalisation scales
......@@ -52,7 +52,7 @@ PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2mu
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15', parameter_values=[ 2.0, 0.5, 90900] )
PowhegConfig.nEvents *= 250. # compensate filter efficiency
PowhegConfig.nEvents *= 25. # compensate filter efficiency
PowhegConfig.generate()
#--------------------------------------------------------------
......
......@@ -18542,6 +18542,8 @@ MC15.410249.Sherpa_221_NNPDF30NNLO_ttbar_AllHadronic_MEPS_NLO.py
MC15.410250.Sherpa_221_NNPDF30NNLO_ttbar_SingleLeptonP_MEPS_NLO.py
MC15.410251.Sherpa_221_NNPDF30NNLO_ttbar_SingleLeptonM_MEPS_NLO.py
MC15.410252.Sherpa_221_NNPDF30NNLO_ttbar_dilepton_MEPS_NLO.py
MC15.410274.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttbar_nonallhad_bfil.py
MC15.410275.aMcAtNloPythia8EvtGen_MEN30NLO_A14N23LO_ttbar_dil_bfil.py
MC15.410500.PowhegPythia8EvtGen_A14_ttbar_hdamp172p5_nonallhad.py
MC15.410501.PowhegPythia8EvtGen_A14_ttbar_hdamp258p75_nonallhad.py
MC15.410502.PowhegPythia8EvtGen_A14_ttbar_hdamp258p75_allhad.py
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