Commit bab7184b authored by Joany Manjarres's avatar Joany Manjarres
Browse files

see ChangeLog

parent 94b1cdf6
2019-10-29 Joany Manjarres <joany@cern.ch>
* tagging MC15JobOptions-01-01-81
* update JO 412123
* update cvs 414005-018
* add JO 411343-345
2019-10-27 Cyril Becot <cyril.becot@cern.ch>
* add ../common/MadGraph/MadGraphControl_Pythia8EvtGen_A14NNPDF30LO_Pseudoscalar_2HDMa_monoH.py
* add 312662-312751
......@@ -14,7 +20,7 @@
2019-10-24 E. M. Lobodzinska <ewelina@mail.desy.de>
* add 364819
* tagging MC15JobOptions-01-01-78
2019-10-23 Cyril Becot <cyril.becot@cern.ch>
* up share/evgeninputfiles.csv
* tagging MC15JobOptions-01-01-77
......
......@@ -170,7 +170,7 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#---------------------------------------------------------------------------------------------------
#### Showering with Pythia 8
#---------------------------------------------------------------------------------------------------
evgenConfig.description = "ttW+j production with EW corrections, LO in QCD"
evgenConfig.description = "ttW+j EW LO production"
evgenConfig.keywords = ["ttW", "SM","top"]
evgenConfig.process = "generate p p > t t~ w j"
evgenConfig.generators = ['MadGraph', 'Pythia8', 'EvtGen']
......
#--------------------------------------------------------------
# EVGEN configuration
# Note: These JO are designed to run Powheg and make an LHE file and to not run a showering
# generator afterwards. Because the current Generate_tf.py requires an output file, we
# need to fake it a bit. We therefore will run Pythia8 on the first event in the LHE file.
# Note, since we do not intend to keep the EVNT file, the JO name doesn't include Pythia8
# The Powheg configuration is identical to DSID 410501
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch']
# This is to fake the system into not spending all its time running Pythia8, since we intend
# to throw out the EVNT file
evgenConfig.minevents = 1
include('PowhegControl/PowhegControl_tt_Common.py')
# Initial settings
PowhegConfig.decay_mode = "t t~ > undecayed" #don't decay the top quarks in Powheg
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
PowhegConfig.PDF = [260000, 25200, 13165, 90900, 265000, 266000, 303400] # NNPDF30_nlo_as_0118, MMHT2014nlo68clas118, CT14nlo_as_0118, PDF4LHC15_nlo_30, NNPDF30_nlo_as_0117, NNPDF30_nlo_as_0119, NNPDF31_nlo_as_0118 - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260001, 260101)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(90901, 90931)) # Include the PDF4LHC15_nlo_30 error set
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25200] )
## Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13165] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 0.5, 90900] )
# Scale variations, NNPDF30_nlo_as_0117
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0117', parameter_values=[ 2.0, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0117', parameter_values=[ 0.5, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0117', parameter_values=[ 1.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_0117', parameter_values=[ 1.0, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 0.5, 265000] )
# Scale variations, NNPDF30_nlo_as_0119
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0119', parameter_values=[ 2.0, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0119', parameter_values=[ 0.5, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 0.5, 266000] )
# Scale variations, NNPDF31_nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF31_NLO_0118', parameter_values=[ 2.0, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF31_NLO_0118', parameter_values=[ 0.5, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 0.5, 303400] )
PowhegConfig.nEvents *= 16500. # we'll request 1 event from Pythia8 just to get the metadata
### MadSpin configuration
# dilepton
PowhegConfig.MadSpin_decays = ["decay t > w+ b, w+ > l+ vl", "decay t~ > w- b~, w- > l- vl~"]
PowhegConfig.MadSpin_decays += ["set ms_dir MadSpin"] # hack to prevent MadSpin from crashing on some file systems
PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]"
PowhegConfig.MadSpin_nFlavours = 5
PowhegConfig.MadSpin_model = "loop_sm-no_b_mass"
# Parameters are hardcoded to avoid madspin picking up it's own defaults.
# Lines taken from t-channel single-top (4fs) Powheg joboption
PowhegConfig.alphaem_inv = 1.323489e+02
PowhegConfig.G_F = 1.166370e-05
PowhegConfig.alphaqcd = 1.184000e-01
PowhegConfig.mass_W = 80.399
PowhegConfig.mass_Z = 9.118760e+01
PowhegConfig.mass_H = 125.0
PowhegConfig.bwcutoff = 50
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#Override the generator names since we won't keep the Pythia8
evgenConfig.generators = ["Powheg"]
#--------------------------------------------------------------
# EVGEN configuration
# Note: These JO are designed to run Powheg and make an LHE file and to not run a showering
# generator afterwards. Because the current Generate_tf.py requires an output file, we
# need to fake it a bit. We therefore will run Pythia8 on the first event in the LHE file.
# Note, since we do not intend to keep the EVNT file, the JO name doesn't include Pythia8
# The Powheg configuration is identical to DSID 410501
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch']
# This is to fake the system into not spending all its time running Pythia8, since we intend
# to throw out the EVNT file
evgenConfig.minevents = 1
include('PowhegControl/PowhegControl_tt_Common.py')
# Initial settings
PowhegConfig.decay_mode = "t t~ > undecayed" #don't decay the top quarks in Powheg
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
PowhegConfig.PDF = [260000, 25200, 13165, 90900, 265000, 266000, 303400] # NNPDF30_nlo_as_0118, MMHT2014nlo68clas118, CT14nlo_as_0118, PDF4LHC15_nlo_30, NNPDF30_nlo_as_0117, NNPDF30_nlo_as_0119, NNPDF31_nlo_as_0118 - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260001, 260101)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(90901, 90931)) # Include the PDF4LHC15_nlo_30 error set
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25200] )
## Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13165] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 0.5, 90900] )
# Scale variations, NNPDF30_nlo_as_0117
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0117', parameter_values=[ 2.0, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0117', parameter_values=[ 0.5, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0117', parameter_values=[ 1.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_0117', parameter_values=[ 1.0, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 0.5, 265000] )
# Scale variations, NNPDF30_nlo_as_0119
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0119', parameter_values=[ 2.0, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0119', parameter_values=[ 0.5, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 0.5, 266000] )
# Scale variations, NNPDF31_nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF31_NLO_0118', parameter_values=[ 2.0, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF31_NLO_0118', parameter_values=[ 0.5, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 0.5, 303400] )
PowhegConfig.nEvents *= 16500. # we'll request 1 event from Pythia8 just to get the metadata
### MadSpin configuration
PowhegConfig.MadSpin_decays = ["decay t > w+ b, w+ > l+ vl", "decay t~ > w- b~, w- > j j"]
PowhegConfig.MadSpin_decays += ["set ms_dir MadSpin"] # hack to prevent MadSpin from crashing on some file systems
PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]"
PowhegConfig.MadSpin_nFlavours = 5
PowhegConfig.MadSpin_model = "loop_sm-no_b_mass"
# Parameters are hardcoded to avoid madspin picking up it's own defaults.
# Lines taken from t-channel single-top (4fs) Powheg joboption
PowhegConfig.alphaem_inv = 1.323489e+02
PowhegConfig.G_F = 1.166370e-05
PowhegConfig.alphaqcd = 1.184000e-01
PowhegConfig.mass_W = 80.399
PowhegConfig.mass_Z = 9.118760e+01
PowhegConfig.mass_H = 125.0
PowhegConfig.bwcutoff = 50
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#Override the generator names since we won't keep the Pythia8
evgenConfig.generators = ["Powheg"]
#--------------------------------------------------------------
# EVGEN configuration
# Note: These JO are designed to run Powheg and make an LHE file and to not run a showering
# generator afterwards. Because the current Generate_tf.py requires an output file, we
# need to fake it a bit. We therefore will run Pythia8 on the first event in the LHE file.
# Note, since we do not intend to keep the EVNT file, the JO name doesn't include Pythia8
# The Powheg configuration is identical to DSID 410501
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8 ttbar production with Powheg hdamp equal 1.5*top mass, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'dominic.hirschbuehl@cern.ch']
# This is to fake the system into not spending all its time running Pythia8, since we intend
# to throw out the EVNT file
evgenConfig.minevents = 1
include('PowhegControl/PowhegControl_tt_Common.py')
# Initial settings
PowhegConfig.decay_mode = "t t~ > undecayed" #don't decay the top quarks in Powheg
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
PowhegConfig.PDF = [260000, 25200, 13165, 90900, 265000, 266000, 303400] # NNPDF30_nlo_as_0118, MMHT2014nlo68clas118, CT14nlo_as_0118, PDF4LHC15_nlo_30, NNPDF30_nlo_as_0117, NNPDF30_nlo_as_0119, NNPDF31_nlo_as_0118 - PDF variations with nominal scale variation
PowhegConfig.PDF.extend(range(260001, 260101)) # Include the NNPDF error set
PowhegConfig.PDF.extend(range(90901, 90931)) # Include the PDF4LHC15_nlo_30 error set
# Define a weight group configuration for scale variations with different PDFs
# Nominal mu_F = mu_R = 1.0 is not required as this is captured by the PDF variation above
PowhegConfig.define_event_weight_group( group_name='scales_pdf', parameters_to_vary=['mu_F','mu_R','PDF'] )
# Scale variations, MMHT2014nlo68clas118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_MMHT', parameter_values=[ 2.0, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_MMHT', parameter_values=[ 0.5, 1.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_MMHT', parameter_values=[ 1.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_MMHT', parameter_values=[ 1.0, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_MMHT', parameter_values=[ 0.5, 0.5, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_MMHT', parameter_values=[ 2.0, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_MMHT', parameter_values=[ 0.5, 2.0, 25200] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_MMHT', parameter_values=[ 2.0, 0.5, 25200] )
## Scale variations, CT14nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_CT14', parameter_values=[ 2.0, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_CT14', parameter_values=[ 0.5, 1.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_CT14', parameter_values=[ 1.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_CT14', parameter_values=[ 1.0, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_CT14', parameter_values=[ 0.5, 0.5, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_CT14', parameter_values=[ 2.0, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_CT14', parameter_values=[ 0.5, 2.0, 13165] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_CT14', parameter_values=[ 2.0, 0.5, 13165] )
# Scale variations, PDF4LHC15_nlo_30
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 1.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 1.0, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 0.5, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_PDF4LHC15_NLO_30', parameter_values=[ 0.5, 2.0, 90900] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_PDF4LHC15_NLO_30', parameter_values=[ 2.0, 0.5, 90900] )
# Scale variations, NNPDF30_nlo_as_0117
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0117', parameter_values=[ 2.0, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0117', parameter_values=[ 0.5, 1.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0117', parameter_values=[ 1.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_0117', parameter_values=[ 1.0, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 0.5, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0117', parameter_values=[ 0.5, 2.0, 265000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0117', parameter_values=[ 2.0, 0.5, 265000] )
# Scale variations, NNPDF30_nlo_as_0119
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF_NLO_0119', parameter_values=[ 2.0, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF_NLO_0119', parameter_values=[ 0.5, 1.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF_NLO_0119', parameter_values=[ 1.0, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 0.5, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF_NLO_0119', parameter_values=[ 0.5, 2.0, 266000] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF_NLO_0119', parameter_values=[ 2.0, 0.5, 266000] )
# Scale variations, NNPDF31_nlo_as_0118
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_NNPDF31_NLO_0118', parameter_values=[ 2.0, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_NNPDF31_NLO_0118', parameter_values=[ 0.5, 1.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muR_NNPDF31_NLO_0118', parameter_values=[ 1.0, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 0.5, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='0p5muF_2muR_NNPDF31_NLO_0118', parameter_values=[ 0.5, 2.0, 303400] )
PowhegConfig.add_weight_to_group( group_name='scales_pdf', weight_name='2muF_0p5muR_NNPDF31_NLO_0118', parameter_values=[ 2.0, 0.5, 303400] )
PowhegConfig.nEvents *= 16500. # we'll request 1 event from Pythia8 just to get the metadata
### MadSpin configuration
PowhegConfig.MadSpin_decays = ["decay t > w+ b, w+ > j j", "decay t~ > w- b~, w- > l- vl~"]
PowhegConfig.MadSpin_decays += ["set ms_dir MadSpin"] # hack to prevent MadSpin from crashing on some file systems
PowhegConfig.MadSpin_process = "generate p p > t t~ [QCD]"
PowhegConfig.MadSpin_nFlavours = 5
PowhegConfig.MadSpin_model = "loop_sm-no_b_mass"
# Parameters are hardcoded to avoid madspin picking up it's own defaults.
# Lines taken from t-channel single-top (4fs) Powheg joboption
PowhegConfig.alphaem_inv = 1.323489e+02
PowhegConfig.G_F = 1.166370e-05
PowhegConfig.alphaqcd = 1.184000e-01
PowhegConfig.mass_W = 80.399
PowhegConfig.mass_Z = 9.118760e+01
PowhegConfig.mass_H = 125.0
PowhegConfig.bwcutoff = 50
PowhegConfig.generate()
#--------------------------------------------------------------
# Pythia8 showering
#--------------------------------------------------------------
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
include("MC15JobOptions/Pythia8_Powheg_Main31.py")
genSeq.Pythia8.Commands += [ 'Powheg:pTHard = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:NFinal = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTdef = 2' ]
genSeq.Pythia8.Commands += [ 'Powheg:veto = 1' ]
genSeq.Pythia8.Commands += [ 'Powheg:vetoCount = 3' ]
genSeq.Pythia8.Commands += [ 'Powheg:pTemt = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:emitted = 0' ]
genSeq.Pythia8.Commands += [ 'Powheg:MPIveto = 0' ]
#Override the generator names since we won't keep the Pythia8
evgenConfig.generators = ["Powheg"]
......@@ -8078,7 +8078,7 @@ Dsid,energy,inputeventfile,(inputconffile)
414006, 13000 ,, group.phys_gener.h7_matchbox.414006.tt_MCatNLO_dipSh_SingleLep.TXT.mc15_v0
414007, 13000 ,, group.phys-gener.h7_matchbox.414007.tt_POW_defOpt_angSh_SingleLep.TXT.mc15_v0
414008, 13000 ,, group.phys-gener.h7_matchbox.414008.tt_POW_defOpt_dipSh_SingleLep.TXT.mc15_v0
414009, 13000 ,, group.phys-gener.h7_matchbox.414009.tt_POW-group.phys-gener.h7_matchbox.414008.tt_POW_defOpt_dipSh_SingleLep.TXT.mc15_v0
414009, 13000 ,, group.phys-gener.h7_matchbox.414009.tt_POW_PWopt_angSh_SingleLep.TXT.mc15_v0
414010, 13000 ,, group.phys-gener.h7_matchbox.414010.tt_POW_PWopt_dipSh_SingleLep.TXT.mc15_v0
414011, 13000 ,, group.phys-gener.h7_matchbox.414011.tt_POW_PWopt_dipSh_NLOdec_SingleLep.TXT.mc15_v0
414012, 13000 ,, group.phys-gener.h7_matchbox.414012.tt_MCatNLO_angSh_dil.TXT.mc15_v0
......
MC15JobOptions-01-01-80/
\ No newline at end of file
MC15JobOptions-01-01-81/
\ 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