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

Add 410255-410256


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@780786 4525493e-7705-40b1-a816-d608a930855b
parent 14426fc3
2016-10-27 Giovanni Siragusa <siragusa@cern.ch>
* Add common/MadGraph/MadGraphControl_tt012FxFx.py
* Add 410255-410256
* Add HI_fragments/BoostFiltered_pPb8TeV.py
* Add HI_fragments/BoostFiltered_Pbp8TeV.py
......
from MadGraphControl.MadGraphUtils import *
from AthenaCommon import Logging
JOlog = Logging.logging.getLogger('FxFxJobOption')
JOlog.info(" Loading MadGraphControl file for tt012 FxFx ")
# General settings
nevents = 5.0*runArgs.maxEvents # 25% efficiency for FxFx matching plus a safety factor
gridpack_dir = None
gridpack_mode = False
runName = 'run_01'
# DSID dictionary - Will allow expansion into a control file in the future
# key - DSID
# arg1 - muQ (qCut) for Pythia8 FxFx
# arg2 - qCutME for Pythia8 FxFx (also MG ptj parameter value)
ttbar_5fl={410255:["70.0","20.0"],
410256:["30.0","15.0"]}
# We can define the madspin decay pattern based on DSID
inclusive = [410255, 410256]
nonallhad = []
dilepton = []
allhad = []
# Using gridpacks, so make sure there is a check on the config file
if hasattr(runArgs,'inputGenConfFile'):
evgenConfig.inputconfcheck = "tt012_FxFx"
if runArgs.runNumber in ttbar_5fl:
name = 'ttbarjets'
process = "pp>tt"
gridpack_mode = True
gridpack_dir = 'madevent/'
else:
raise RuntimeError("runNumber %i not recognised in these jobOptions."%runArgs.runNumber)
msdecay = ""
nleptons = -2
if runArgs.runNumber in inclusive:
nleptons = -2
msdecay = "decay t > w+ b, w+ > all all \ndecay t~ > w- b~, w- > all all"
elif runArgs.runNumber in nonallhad:
nleptons = -1
nevents = 3.0*runArgs.maxEvents # Filter efficiency
msdecay = "decay t > w+ b, w+ > all all \ndecay t~ > w- b~, w- > all all"
elif runArgs.runNumber in dilepton:
nleptons = 2
msdecay = "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"
elif runArgs.runNumber in allhad:
nleptons = 0
msdecay = "decay t > w+ b, w+ > j j \ndecay t~ > w- b~, w- > j j"
else:
raise RuntimeError("Do not know which decay is required for this DSID")
JOlog.info(" MadSpin decay defined as "+msdecay)
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(name)
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
import model loop_sm-no_b_mass
define p = p b b~
define j = p
define l+ = e+ mu+
define l- = e- mu-
import model loop_sm-no_b_mass
generate p p > t t~ [QCD] @0
add process p p > t t~ j [QCD] @1
add process p p > t t~ j j [QCD] @2
output ttbar_FxFx_2jets
""")
fcard.close()
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
process_dir = new_process(grid_pack=gridpack_dir)
#Fetch default LO run_card.dat and set parameters
extras = { 'lhe_version' :'3.0',
#'req_acc' :'0.005',
'req_acc_FO' :'0.001',
'pdlabel' :"'lhapdf'",
'lhaid' :260000,
'parton_shower' :'PYTHIA8',
'cut_decays' :'F',
'reweight_scale':'True',
'reweight_PDF' :'True',
'PDF_set_min' :260001,
'PDF_set_max' :260100,
'ickkw' :'3',
'jetalgo' :'1',
'jetradius' :'1.0',
'ptj' : ttbar_5fl[runArgs.runNumber][1],
'etaj' :'5',
'ptl' :'20',
'etal' :'3',
'mll' :'40',
'mll_sf' :'40',
'isoEM' :'True',
'maxjetflavor' :'5',
'use_syst' :"False"
}
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,xqcut=0.,
extras=extras)
print_cards()
madspin_card_loc='madspin_card.dat'
mscard = open(madspin_card_loc,'w')
if not hasattr(runArgs, 'inputGenConfFile'):
fMadSpinCard = open('madspin_card.dat','w')
fMadSpinCard.write('import Events/'+runName+'/events.lhe.gz\n')
fMadSpinCard.write('set ms_dir MadSpin\n')
else:
os.unlink(gridpack_dir+'Cards/madspin_card.dat')
fMadSpinCard = open(gridpack_dir+'Cards/madspin_card.dat','w')
fMadSpinCard.write('import '+gridpack_dir+'Events/'+runName+'/events.lhe.gz\n')
fMadSpinCard.write('set ms_dir '+gridpack_dir+'MadSpin\n')
fMadSpinCard.write('set seed '+str(10000000+int(runArgs.randomSeed))+'\n')
fMadSpinCard.write('''set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
%s
launch'''%(msdecay))
fMadSpinCard.close()
get_param_card = subprocess.Popen(['get_files','-data','aMcAtNlo_param_card_loopsmnobmass.dat'])
get_param_card.wait()
print_cards()
generate(required_accuracy = 0.005,
run_card_loc = 'run_card.dat',
param_card_loc = 'aMcAtNlo_param_card_loopsmnobmass.dat',
njobs = 21,
proc_dir = process_dir,
run_name = runName,
madspin_card_loc = madspin_card_loc,
grid_pack = gridpack_mode,
gridpack_dir = gridpack_dir,
nevents = nevents,
random_seed = runArgs.randomSeed,
mode = 2,
cluster_type = 'pbs')
arrange_output(run_name = runName,
proc_dir = process_dir,
outputDS = runName+'._00001.events.tar.gz',
lhe_version = 3)
# Shower
evgenConfig.description = 'aMcAtNlo_ttbar_fxfx inclusive'
evgenConfig.contact = ['francesco.giuli@cern.ch','mcfayden@cern.ch','ian.connelly@cern.ch']
evgenConfig.keywords += ['ttbar','jets','top']
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile = runName+'._00001.events.tar.gz'
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_aMcAtNlo.py")
testSeq.TestHepMC.MaxVtxDisp = 100000000
# Matching settings
genSeq.Pythia8.Commands += ["JetMatching:merge = on",
"JetMatching:scheme = 1",
"JetMatching:setMad = off",
"JetMatching:qCut = "+ttbar_5fl[runArgs.runNumber][0], # Dependent on DSID
"JetMatching:coneRadius = 1.0",
"JetMatching:etaJetMax = 10.0",
"JetMatching:doFxFx = on",
"JetMatching:qCutME = "+ttbar_5fl[runArgs.runNumber][1], # Dependent on DSID
"JetMatching:nJetMax = 2" ]
# Bit more logging info
JOlog.info("JetMatching:qCut -> "+ttbar_5fl[runArgs.runNumber][0])
JOlog.info("JetMatching:qCutME -> "+ttbar_5fl[runArgs.runNumber][1])
JOlog.info("aMcAtNlo:ptj -> "+ttbar_5fl[runArgs.runNumber][1])
# Pythia8_i settings
genSeq.Pythia8.UserHook = 'JetMatchingMadgraph'
genSeq.Pythia8.FxFxXS = True
# --------------------------------------------------------------
# Apply TTbarWToLeptonFilter
# --------------------------------------------------------------
# Primarily needed for non-all hadronic decays
# nleptons == -2 means inclusive decay ergo no filtering on leptons
if nleptons != -2:
include("MC15JobOptions/TTbarWToLeptonFilter.py")
filtSeq.TTbarWToLeptonFilter.NumLeptons = nleptons
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
......@@ -18544,6 +18544,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.410255.aMcAtNloPythia8EvtGen_A14N23LO_MEN30NLO_tt012FxFx_70_inc.py
MC15.410256.aMcAtNloPythia8EvtGen_A14N23LO_MEN30NLO_tt012FxFx_30_inc.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
......@@ -4106,6 +4106,8 @@ DSID,energy,inputeventfile,(inputconffile)
410250, 13000,, group.phys-gener.sherpa020201.410249.Sherpa_221_NNPDF30NNLO_ttbar_AllHadronic_MEPS_NLO_13TeV.TXT.mc15_v1
410251, 13000,, group.phys-gener.sherpa020201.410249.Sherpa_221_NNPDF30NNLO_ttbar_AllHadronic_MEPS_NLO_13TeV.TXT.mc15_v1
410252, 13000,, group.phys-gener.sherpa020201.410249.Sherpa_221_NNPDF30NNLO_ttbar_AllHadronic_MEPS_NLO_13TeV.TXT.mc15_v1
410255, 13000,, group.phys-gener.MG5_aMcAtNlo.410255.tt012_FxFx_inclusive.TXT.mc15_v4
410256, 13000,, group.phys-gener.MG5_aMcAtNlo.410256.tt012_FxFx_inclusive.TXT.mc15_v1
420021, 5020,, group.phys-gener.powheg.420021.jetjet_JZ1_5TeV.TXT.mc15_v1
420022, 5020,, group.phys-gener.powheg.420022.jetjet_JZ2_5TeV.TXT.mc15_v1
420023, 5020,, group.phys-gener.powheg.420023.jetjet_JZ3_5TeV.TXT.mc15_v1
......
Supports Markdown
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