Commit 29ba5849 authored by Joany Manjarres's avatar Joany Manjarres
Browse files

See ChangeLog

parent 23cfeb80
2020-01-25 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-02-55
* add JO 346877 451496
* update cvs 346877
* update JO 451324-338
2020-01-24 E. M. Lobodzinska
* tag MC15JobOptions-01-02-54
* correct MadGraphControl_NNPDF30NLO_lvlvjjEW6_EFTdim8.py
......
#--------------------------------------------------------------
# Herwig showering
#--------------------------------------------------------------
# initialize Herwig7 generator configuration for showering of LHE files
include("MC15JobOptions/Herwig7_LHEF.py")
# configure Herwig7
Herwig7Config.me_pdf_commands(order="NLO", name="NNPDF30_nnlo_as_0118")
Herwig7Config.lhef_powhegbox_commands(lhe_filename=runArgs.inputGeneratorFile, me_pdf_order="NLO")
Herwig7Config.tune_commands()
# add EvtGen
include("MC15JobOptions/Herwig71_EvtGen.py")
# only consider H->WW devays
Herwig7Config.add_commands("""
# force H->WW decays
do /Herwig/Particles/h0:SelectDecayModes h0->W+,W-;
# print out Higgs decays modes and branching ratios to the terminal/log.generate
do /Herwig/Particles/h0:PrintDecayModes
# force W->enu/munu decays
do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+; W+->nu_tau,tau+;
do /Herwig/Particles/W-:SelectDecayModes W-->nu_ebar,e-; W-->nu_mubar,mu-; W-->nu_taubar,tau-;
# print out W decays modes and branching ratios to the terminal/log.generate
do /Herwig/Particles/W+:PrintDecayModes
do /Herwig/Particles/W-:PrintDecayModes
""")
# run Herwig7
Herwig7Config.run()
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'Powheg+Herwig7 VBF Higgs production with H7UE tune'
evgenConfig.keywords = [ "SM", "Higgs", "VBF", "WW", "mH125" ]
evgenConfig.generators += [ 'Powheg', 'Herwig7' ]
evgenConfig.contact = [ 'spyridon.argyropoulos@cern.ch' ]
# 6600 events per file - filtering efficiency ~50% => 3k events per file
evgenConfig.minevents = 10000
evgenConfig.inputFilesPerJob = 4
evgenConfig.tune = "H7.1-Default"
#--------------------------------------------------------------
# Dilepton filter
#--------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import MultiLeptonFilter
filtSeq += MultiLeptonFilter("Multi1TLeptonFilter")
filtSeq += MultiLeptonFilter("Multi2LLeptonFilter")
Multi1TLeptonFilter = filtSeq.Multi1TLeptonFilter
Multi1TLeptonFilter.Ptcut = 15000.
Multi1TLeptonFilter.Etacut = 5.0
Multi1TLeptonFilter.NLeptons = 1
Multi2LLeptonFilter = filtSeq.Multi2LLeptonFilter
Multi2LLeptonFilter.Ptcut = 5000.
Multi2LLeptonFilter.Etacut = 5.0
Multi2LLeptonFilter.NLeptons = 2
filtSeq.Expression = "(Multi1TLeptonFilter) and (Multi2LLeptonFilter)"
......@@ -89,7 +89,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -87,7 +87,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2e)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -90,7 +90,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2e)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -88,7 +88,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2mu)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
......@@ -89,7 +89,8 @@ evgenConfig.keywords = [ "BSM", "Higgs", "BSMHiggs", "mH125"]
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.minevents = 10000
evgenConfig.minevents = 40000
evgenConfig.inputFilesPerJob = 5
evgenConfig.description = "POWHEG+PYTHIA8, ggH H->aZ -> (Z -> ll, a ->2e)"
evgenConfig.keywords = [ "SM", "Higgs", "SMHiggs","Z", "mH125" ]
evgenConfig.contact = [ 'trevor.willard.nelson@cern.ch' ] #Add my email to this thing
......
from MadGraphControl.MadGraphUtils import *
from os import path
safefactor=1.1
mode=0
### DSID lists (extensions can include systematics samples)
test=[451496]
fcard = open('proc_card_mg5.dat','w')
if runArgs.runNumber in test:
fcard.write("""
import model 2HDMCPF_NLO
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+
define l- = e- mu-
define vl = ve vm
define vl~ = ve~ vm~
generate p p > h3 > ta+ ta- j [QCD]
output -f""")
fcard.close()
else:
raise RuntimeError("runNumber %i not recognised in these jobOptions."%runArgs.runNumber)
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
#Fetch default LO run_card.dat and set parameters
extras = { 'lhe_version':'2.0',
'cut_decays':'F',
'pdlabel':"'lhapdf'",
'lhaid' :"263000",
'use_syst':"True",
'sys_pdf' : 'NNPDF30_lo_as_0130',
'ptj1min' : '120'
}
XMass = { '25':'1.2500e+02',
'35':'1.2500e+10',
'36':'20',
'37':'1.2500e+10',
}
Params={}
decays={'WH3':'1.801e-02',
}
Params['DECAY']=decays
runName='run_01'
process_dir = new_process()
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=runArgs.maxEvents/0.53*safefactor,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
build_param_card(param_card_old=path.join(process_dir,'Cards/param_card.dat'),param_card_new='param_card_new.dat',masses=XMass,params=Params)
print_cards()
generate(run_card_loc='run_card.dat',param_card_loc='param_card_new.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=3,saveProcDir=True)
#### Shower
evgenConfig.description = 'MG5 ggF X to ta+ ta-, mX = 20 GeV'
evgenConfig.keywords+=['jets','tau']
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
......@@ -3069,6 +3069,7 @@ Dsid,energy,inputeventfile,(inputconffile)
346874 , 13000, , group.phys-gener.sherpa020201.363356.Sherpa_221_NNPDF30NNLO_ZqqZll_13TeV.TXT.mc15_v1
346875 , 13000, , group.phys-gener.sherpa020201.363356.Sherpa_221_NNPDF30NNLO_ZqqZll_13TeV.TXT.mc15_v1
346876 , 13000, , group.phys-gener.sherpa020201.363357.Sherpa_221_NNPDF30NNLO_WqqZvv_13TeV.TXT.mc15_v1
346877 , 13000, group.phys-gener.powheg.344235.Powheg_NNPDF30_VBFH125_13TeV.TXT.mc15_v1
346878, 13000, mc15_13TeV.345916.Powheg_NNPDF30_VBFH125_LHE.evgen.TXT.e6901
346879, 13000, mc15_13TeV.345038.PowhegPythia8EvtGen_NNPDF30_AZNLO_ZH125J_Zincl_MINLO.evgen.TXT.e5590
346880, 13000, mc15_13TeV.345040.PowhegPythia8EvtGen_NNPDF30_AZNLO_WmH125J_Wincl_MINLO.evgen.TXT.e5590
......
MC15JobOptions-01-02-54/
\ No newline at end of file
MC15JobOptions-01-02-55/
\ 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