Commit 5e2c04ed authored by Cyril Pascal Becot's avatar Cyril Pascal Becot
Browse files

See ChangeLog


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@816233 4525493e-7705-40b1-a816-d608a930855b
parent 10513e84
2018-12-13 Cyril Becot <cyril.becot@cern.ch>
* tag MC15JobOptions-00-12-05
* add 450190,450200,450210,450247
* up MadGraphControl_SimplifiedModel_C1C1_Stau_Filter.py MadGraphControl_SimplifiedModel_C1N2_Stau_Filter.py
* up MadGraphPythia8EvtGenControl_Hplus_HighMass_NLO.py MadGraphPythia8EvtGenControl_Hplus_IntMass_LO.py MadGraphPythia8EvtGenControl_Hplus_LowMass_LO.py
* up 346213
* up MadGraphControl_SimplifiedModel_GG_mixedC1LLP.py
* up 346189
* add 413002,413003
* add 436200-436205
* up 345919-345923, 345928
* up 450028,450029
* add SLHA files for 448390-448499
* add 450185-4501859
* up 450005-450013
* add 450082-450111
* add 425400,425401
* up MadGraphPythia8EvtGenControl_Hplus_HighMass_NLO.py MadGraphPythia8EvtGenControl_Hplus_IntMass_LO.py MadGraphPythia8EvtGenControl_Hplus_LowMass_LO.py
* add 450112-450115
* up MadGraphControl_SimplifiedModel_higgsino.py
* up MadGraphControl_NUHM2_higgsino_MadSpin.py
* add PhPy8EG_A14_Wt_inclusive_filt_control.py
* add 411181-411204
* add 436072,436206-436235
* up 436200-436205
2018-12-13 E. M. Lobodzinska <ewelina@mail.desy.de>
* add 363749-61 and two H7 commons
* tag MC15JobOptions-00-12-04
......
......@@ -168,13 +168,15 @@ if '2L2MET75' in splitConfig[-2]:
filtSeq.Expression = "MultiElecMuTauFilter and MissingEtFilter"
evt_multiplier = 100
evt_multiplier = 10
if gentype=="C1C1" and int(m12)>=500:
evt_multiplier = 150
evt_multiplier = 15
if int(m12)>=6000:
evt_multiplier = 500
evt_multiplier = 50
if gentype=="N2N1":
evgenConfig.minevents = 1000
else:
evgenConfig.minevents = 2000
#--------------------------------------------------------------
# Pass everything to the PostInclude JO file
......
......@@ -11,8 +11,8 @@ masses['1000024'] = MassToFloat(splits[4])
masses['1000023'] = MassToFloat(splits[4])
masses['1000022'] = MassToFloat(splits[5].split('.')[0])
splitting = 0.5
if len(splits) > 6:
splitting = MassToFloat(splits[6].split('.')[0])
if len(splits) > 7:
splitting = MassToFloat(splits[-2].split('.')[0])
evgenLog.info("Using mass splitting of %f for mslep" % splitting)
mslep = ( masses['1000024'] - masses['1000022'] )*splitting + masses['1000022']
evgenLog.info(" Stau / sneutrino mass is %f GeV" % mslep)
......@@ -38,6 +38,9 @@ evgenConfig.description = '~chi1+/~chi1- production, decay via stau in simplifie
# evt_multiplier = 3 is only sufficient for large mass splittings
evt_multiplier = 20
extras['event_norm']='sum'
extras['use_syst']='F'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
......@@ -48,22 +51,24 @@ if njets>0:
#--------------------------------------------------------------
# Lepton filter
#--------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilter")
filtSeq += MultiElecMuTauFilter("TauFilter")
MultiElecMuTauFilter1 = filtSeq.DileptonFilter
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 5000.
MultiElecMuTauFilter1.MaxEta = 2.7
MultiElecMuTauFilter1.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter1.IncludeHadTaus = 1 # include hadronic taus
MultiElecMuTauFilter2 = filtSeq.TauFilter
MultiElecMuTauFilter2.NLeptons = 1
MultiElecMuTauFilter2.MinPt = 1e10
MultiElecMuTauFilter2.MaxEta = 2.7
MultiElecMuTauFilter2.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter2.IncludeHadTaus = 1 # include hadronic taus
filtSeq.Expression = "(DileptonFilter and TauFilter)"
if '2L1T' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('2L1T filter is applied')
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilter")
filtSeq += MultiElecMuTauFilter("TauFilter")
MultiElecMuTauFilter1 = filtSeq.DileptonFilter
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 5000.
MultiElecMuTauFilter1.MaxEta = 2.7
MultiElecMuTauFilter1.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter1.IncludeHadTaus = 1 # include hadronic taus
MultiElecMuTauFilter2 = filtSeq.TauFilter
MultiElecMuTauFilter2.NLeptons = 1
MultiElecMuTauFilter2.MinPt = 1e10
MultiElecMuTauFilter2.MaxEta = 2.7
MultiElecMuTauFilter2.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter2.IncludeHadTaus = 1 # include hadronic taus
filtSeq.Expression = "(DileptonFilter and TauFilter)"
......@@ -11,8 +11,8 @@ masses['1000024'] = MassToFloat(splits[4])
masses['1000023'] = MassToFloat(splits[4])
masses['1000022'] = MassToFloat(splits[5].split('.')[0])
splitting = 0.5
if len(splits) > 6:
splitting = MassToFloat(splits[6].split('.')[0])
if len(splits) > 7:
splitting = MassToFloat(splits[-2].split('.')[0])
evgenLog.info("Using mass splitting of %f for mslep" % splitting)
mslep = ( masses['1000024'] - masses['1000022'] )*splitting + masses['1000022']
evgenLog.info(" Stau / sneutrino mass is %f GeV" % mslep)
......@@ -41,31 +41,43 @@ evgenConfig.description = '~chi1+/- ~chi20 production, decay via stau in simplif
# evt_multiplier = 3 is only sufficient for large mass splittings
evt_multiplier = 18 # 10 was too small for MC15.999002.MadGraphPythia8EvtGen_A14NNPDF23LO_C1N2_Stau_217p5_182p5.py
extras['event_norm']='sum'
extras['use_syst']='F'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
genSeq.Pythia8.Commands += [ "Merging:Process = pp>{x1+,1000024}{x1-,-1000024}{n2,1000023}",
"1000024:spinType = 1",
"1000023:spinType = 1" ]
# New feature in Pythia 8.230+
genSeq.Pythia8.Commands += ["Merging:Process = guess"]
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
genSeq.Pythia8.UserHooks += ['JetMergingaMCatNLO']
else:
genSeq.Pythia8.UserHook = 'JetMergingaMCatNLO'
#
genSeq.Pythia8.Commands += [ "1000024:spinType = 1",
"1000023:spinType = 1" ]
#--------------------------------------------------------------
# Lepton filter
#--------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilter")
filtSeq += MultiElecMuTauFilter("TauFilter")
MultiElecMuTauFilter1 = filtSeq.DileptonFilter
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 5000.
MultiElecMuTauFilter1.MaxEta = 2.7
MultiElecMuTauFilter1.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter1.IncludeHadTaus = 1 # include hadronic taus
MultiElecMuTauFilter2 = filtSeq.TauFilter
MultiElecMuTauFilter2.NLeptons = 1
MultiElecMuTauFilter2.MinPt = 1e10
MultiElecMuTauFilter2.MaxEta = 2.7
MultiElecMuTauFilter2.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter2.IncludeHadTaus = 1 # include hadronic taus
filtSeq.Expression = "(DileptonFilter and TauFilter)"
if '2L1T' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('2L1T filter is applied')
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("DileptonFilter")
filtSeq += MultiElecMuTauFilter("TauFilter")
MultiElecMuTauFilter1 = filtSeq.DileptonFilter
MultiElecMuTauFilter1.NLeptons = 2
MultiElecMuTauFilter1.MinPt = 5000.
MultiElecMuTauFilter1.MaxEta = 2.7
MultiElecMuTauFilter1.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter1.IncludeHadTaus = 1 # include hadronic taus
MultiElecMuTauFilter2 = filtSeq.TauFilter
MultiElecMuTauFilter2.NLeptons = 1
MultiElecMuTauFilter2.MinPt = 1e10
MultiElecMuTauFilter2.MaxEta = 2.7
MultiElecMuTauFilter2.MinVisPtHadTau = 15000. # pt-cut on the visible hadronic tau
MultiElecMuTauFilter2.IncludeHadTaus = 1 # include hadronic taus
filtSeq.Expression = "(DileptonFilter and TauFilter)"
......@@ -38,6 +38,10 @@ evgenConfig.contact = ["chihiro.kozakai@cern.ch"]
evgenConfig.keywords += ['simplifiedModel','gluino', 'chargino', 'longLived']
evgenConfig.description = 'gluino production, glu->qq+wino in simplified model, m_glu = %s GeV, m_wino = %s GeV, lifetime = %s ns' % (masses['1000021'], masses['1000022'], lifetime)
# AGENE-1542: disable syscalc and set event_norm explicitly
extras['event_norm']='sum'
extras['use_syst']='F'
include ('MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py')
print "---"
......
......@@ -316,11 +316,11 @@ if '1L' in config_after_UEtune:
evt_multiplier *= 1
evgenConfig.minevents = 5000
elif int(dM)<=2:
evt_multiplier *= 1.
evt_multiplier *= 1.5
evgenConfig.minevents = 5000
if abs(masses['1000023'])>=300:
evt_multiplier *= 1.3
evt_multiplier *= 1.25
from GeneratorFilters.GeneratorFiltersConf import MultiElecMuTauFilter
filtSeq += MultiElecMuTauFilter("MultiElecMuTauFilter")
ElecMuTauFilter = filtSeq.MultiElecMuTauFilter
......
......@@ -187,7 +187,7 @@ else:
generate(run_card_loc='run_card.dat',param_card_loc='param_card.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)
outputDS=arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=3,saveProcDir=True)
#### Shower
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -105,7 +105,7 @@ runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='param_card.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)
outputDS=arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=3,saveProcDir=True)
#### Shower
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -151,8 +151,8 @@ print_cards()
runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='param_card.dat',mode=mode,proc_dir=process_dir,run_name=runName)
arrange_output(run_name=runName,lhe_version=3,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
outputDS=arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz',lhe_version=3,saveProcDir=True)
#### Shower
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
#--------------------------------------------------------------
# EVGEN configuration
#--------------------------------------------------------------
evgenConfig.description = 'POWHEG+Pythia8+EvtGen Wt production, HT/MET Filter, DR/DS scheme, inclusive, with Powheg hdamp equal 1.5*top mass, A14 tune, ME NNPDF30 NLO, A14 NNPDF23 LO'
evgenConfig.keywords = [ 'SM', 'top', 'singleTop', 'Wt', 'inclusive']
evgenConfig.contact = [ 'hass.abouzeid@cern.ch' ]
#--------------------------------------------------------------
# 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' ]
#--------------------------------------------------------------
# Event filter
#--------------------------------------------------------------
# Depending on the DSID, set up MET/HT filters
thisDSID = runArgs.runNumber
HT1k5_filter = [411183,411186,411189,411192]
HT1k_1k5_filter = [411182,411185,411188,411191]
HT6c_1k_filter = [411181,411184,411187,411190]
MET200_filter = [411193,411196,411199,411202]
MET300_filter = [411194,411197,411200,411203]
MET400_filter = [411195,411198,411201,411204]
if thisDSID in MET200_filter:
include('MC15JobOptions/MissingEtFilter.py')
filtSeq.MissingEtFilter.METCut = 200*GeV
filtSeq.MissingEtFilterUpperCut.METCut = 300*GeV
if thisDSID in MET300_filter:
include('MC15JobOptions/MissingEtFilter.py')
filtSeq.MissingEtFilter.METCut = 300*GeV
filtSeq.MissingEtFilterUpperCut.METCut = 400*GeV
if thisDSID in MET400_filter:
include('MC15JobOptions/MissingEtFilter.py')
filtSeq.MissingEtFilter.METCut = 400*GeV
if thisDSID in HT6c_1k_filter:
include('MC15JobOptions/HTFilter.py')
filtSeq.HTFilter.MinJetPt = 35.*GeV # Min pT to consider jet in HT
filtSeq.HTFilter.MaxJetEta = 2.5 # Max eta to consider jet in HT
filtSeq.HTFilter.MinHT = 600.*GeV # Min HT to keep event
filtSeq.HTFilter.MaxHT = 1000.*GeV # Max HT to keep event
filtSeq.HTFilter.UseLeptonsFromWZTau = True # Include e/mu from the MC event in the HT
filtSeq.HTFilter.MinLeptonPt = 25.*GeV # Min pT to consider muon in HT
filtSeq.HTFilter.MaxLeptonEta = 2.5 # Max eta to consider muon in HT
if thisDSID in HT1k_1k5_filter:
include('MC15JobOptions/HTFilter.py')
filtSeq.HTFilter.MinJetPt = 35.*GeV # Min pT to consider jet in HT
filtSeq.HTFilter.MaxJetEta = 2.5 # Max eta to consider jet in HT
filtSeq.HTFilter.MinHT = 1000.*GeV # Min HT to keep event
filtSeq.HTFilter.MaxHT = 1500.*GeV # Max HT to keep event
filtSeq.HTFilter.UseLeptonsFromWZTau = True # Include e/mu from the MC event in the HT
filtSeq.HTFilter.MinLeptonPt = 25.*GeV # Min pT to consider muon in HT
filtSeq.HTFilter.MaxLeptonEta = 2.5 # Max eta to consider muon in HT
if thisDSID in HT1k5_filter:
include('MC15JobOptions/HTFilter.py')
filtSeq.HTFilter.MinJetPt = 35.*GeV # Min pT to consider jet in HT
filtSeq.HTFilter.MaxJetEta = 2.5 # Max eta to consider jet in HT
filtSeq.HTFilter.MinHT = 1500.*GeV # Min HT to keep event
filtSeq.HTFilter.MaxHT = 20000.*GeV # Max HT to keep event
filtSeq.HTFilter.UseLeptonsFromWZTau = True # Include e/mu from the MC event in the HT
filtSeq.HTFilter.MinLeptonPt = 25.*GeV # Min pT to consider muon in HT
filtSeq.HTFilter.MaxLeptonEta = 2.5 # Max eta to consider muon in HT
......@@ -99,8 +99,9 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#### Shower
evgenConfig.description = 'MG5 ttbarX to ta+ ta-, mX = 20 GeV'
evgenConfig.keywords+=['ttbar','jets','tau']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#evgenConfig.inputfilecheck = runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
runArgs.inputGeneratorFile=outputDS
evgenConfig.contact = ['Zirui Wang <zirui.wang@cern.ch>']
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -99,8 +99,9 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#### Shower
evgenConfig.description = 'MG5 ttbarX to ta+ ta-, mX = 40 GeV'
evgenConfig.keywords+=['ttbar','jets','tau']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#evgenConfig.inputfilecheck = runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
runArgs.inputGeneratorFile=outputDS
evgenConfig.contact = ['Zirui Wang <zirui.wang@cern.ch>']
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -99,8 +99,9 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#### Shower
evgenConfig.description = 'MG5 ttbarX to ta+ ta-, mX = 80 GeV'
evgenConfig.keywords+=['ttbar','jets','tau']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#evgenConfig.inputfilecheck = runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
runArgs.inputGeneratorFile=outputDS
evgenConfig.contact = ['Zirui Wang <zirui.wang@cern.ch>']
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -99,8 +99,9 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#### Shower
evgenConfig.description = 'MG5 ttbarX to ta+ ta-, mX = 100 GeV'
evgenConfig.keywords+=['ttbar','jets','tau']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#evgenConfig.inputfilecheck = runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
runArgs.inputGeneratorFile=outputDS
evgenConfig.contact = ['Zirui Wang <zirui.wang@cern.ch>']
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
......@@ -99,8 +99,9 @@ arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.e
#### Shower
evgenConfig.description = 'MG5 ttbarX to ta+ ta-, mX = 60 GeV'
evgenConfig.keywords+=['ttbar','jets','tau']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
#evgenConfig.inputfilecheck = runName
#runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
runArgs.inputGeneratorFile=outputDS
evgenConfig.contact = ['Zirui Wang <zirui.wang@cern.ch>']
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
......
include('MadGraphControl_ttH_gamgam_NLO.py')
include('MC15JobOptions/MadGraphControl_ttH_gamgam_NLO.py')
genSeq.Pythia8.Commands += [
'25:onMode = off', # switch OFF all Higgs decay channels
'25:onIfMatch = 22 22' # H -> gamma gamma
]
evgenConfig.inputconfcheck = "TXT"
evgenConfig.minevents = 2000
evgenConfig.inputFilesPerJob = 14
include('MC15JobOptions/PhPy8EG_A14_Wt_inclusive_filt_control.py')
evgenConfig.minevents = 200
evgenConfig.inputFilesPerJob = 6
include('MC15JobOptions/PhPy8EG_A14_Wt_inclusive_filt_control.py')
evgenConfig.minevents = 100
evgenConfig.inputFilesPerJob = 13
include('MC15JobOptions/PhPy8EG_A14_Wt_inclusive_filt_control.py')
evgenConfig.minevents = 2000
evgenConfig.inputFilesPerJob = 14
include('MC15JobOptions/PhPy8EG_A14_Wt_inclusive_filt_control.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