Commit 1fce184f authored by mcjo's avatar mcjo
Browse files

add 310480, 406456-503, 391934-67, 414000


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@815101 4525493e-7705-40b1-a816-d608a930855b
parent b0359d7d
2018-08-15 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-10-82
* add 310480, 406456-503, 391934-67, 414000
* update MadGraphControl_SimplifiedModel_T1T1_onestepN2N2.py
2018-08-13 Christian Gutschow <chris.g@cern.ch>
* change b-quark mass to 4.95 in 429755-8
* tagging MC15JobOptions-00-10-81
......
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000021'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t t1~ go $ susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
'''
njets = 2
evgenLog.info('Generation of anti-stop top gluino production, stop to bffN; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
#jobConfigParts = runArgs.jobConfig[0].split("_")
if 'MadSpin' in runArgs.jobConfig[0].split('_')[-2]:
evgenLog.info('Running w/ MadSpin option')
madspin_card='madspin_card_test.dat'
mscard = open(madspin_card,'w')
decay_chains = [ "t1~ > n1 fu~ fd b~" , "t > b W+, W+ > qu qd~, W+ > ll+ vl" ]
# evgenLog.info('Setting only leptonic decay chains')
# decay_chains = [ "t1 > n1 ll+ vl b" , "t1~ > n1 ll- vl~ b~" ]
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
set spinmode none
# specify the decay for the final state particles
define ll+ = e+ mu+ ta+
define vl = ve vm vt
define qu = u c
define qd~ = d~ s~
decay %s
decay %s
# running the actual code
launch"""%(runArgs.randomSeed,decay_chains[0],decay_chains[1]))
mscard.close()
if 'MET100' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('MET100 filter is applied')
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 100*GeV
evt_multiplier = 100.
elif 'M1501L20' in runArgs.jobConfig[0].split('_')[-1]:
evgenLog.info('1Lepton or MET100 filter is applied')
include ( 'MC15JobOptions/MultiElecMuTauFilter.py' )
filtSeq.MultiElecMuTauFilter.MinPt = 20000.
filtSeq.MultiElecMuTauFilter.MaxEta = 2.8
filtSeq.MultiElecMuTauFilter.NLeptons = 1
filtSeq.MultiElecMuTauFilter.IncludeHadTaus = 0
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 150000.
filtSeq.Expression = "MultiElecMuTauFilter or MissingEtFilter"
evt_multiplier = 5.
evgenConfig.contact = [ "paola.arrubarrena@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop']
evgenConfig.description = 'top stop gluino production, st->bffN in simplified model, m_guino = %s GeV, m_stop = %s GeV'%(masses['1000021'],masses['1000006'])
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
# genSeq.Pythia8.Commands += ["Merging:Process = pp>{t,6}{t1~,-1000006}{go,1000021}"]
evgenLog.info('Performing the Pythia guess CKKW-L merging feature')
# id:all = name antiName spinType chargeType colType m0 mWidth mMin mMax tau0
# genSeq.Pythia8.Commands += ["10000104:all = qb4p5 qb4p5bar 2 13.5 0 900"]
genSeq.Pythia8.Commands += ["Merging:Process = guess"]
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
genSeq.Pythia8.UserHooks += ['JetMergingaMCatNLO']
else:
genSeq.Pythia8.UserHook = 'JetMergingaMCatNLO'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t t1~ n1 $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
'''
njets = 2
evgenLog.info('Registered generation of top antistop neutralino production, stop to bffN; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
#jobConfigParts = runArgs.jobConfig[0].split("_")
if 'MadSpin' in runArgs.jobConfig[0].split('_')[-2]:
evgenLog.info('Running w/ MadSpin option')
madspin_card='madspin_card_test.dat'
mscard = open(madspin_card,'w')
decay_chains = [ "t1~ > n1 fu~ fd b~" , " t > b W+, W+ > ll+ vl, W+ > qu qd~"]
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
set spinmode none
# specify the decay for the final state particles
define ll+ = e+ mu+
define vl = ve vm
define qu = u c
define qd~ = d~ s~
decay %s
decay %s
# running the actual code
launch"""%(runArgs.randomSeed,decay_chains[0],decay_chains[1]))
mscard.close()
if 'MET100' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('MET100 filter is applied')
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 100*GeV
evt_multiplier = 50.
elif 'M1501L20' in runArgs.jobConfig[0].split('_')[-1]:
evgenLog.info('1Lepton or MET100 filter is applied')
include ( 'MC15JobOptions/MultiElecMuTauFilter.py' )
filtSeq.MultiElecMuTauFilter.MinPt = 20000.
filtSeq.MultiElecMuTauFilter.MaxEta = 2.8
filtSeq.MultiElecMuTauFilter.NLeptons = 1
filtSeq.MultiElecMuTauFilter.IncludeHadTaus = 0
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 150000.
filtSeq.Expression = "MultiElecMuTauFilter or MissingEtFilter"
evt_multiplier = 5.
evgenConfig.contact = [ "paola.arrubarrena@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop']
evgenConfig.description = 'top stop neutralino production, st->bffN in simplified model, m_stop = %s GeV, m_neutralino = %s GeV'%(masses['1000006'],masses['1000022'])
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
# genSeq.Pythia8.Commands += ["Merging:Process = pp>{t,6}{t1~,-1000006}{go,1000021}"]
evgenLog.info('Performing the Pythia guess CKKW-L merging feature')
# id:all = name antiName spinType chargeType colType m0 mWidth mMin mMax tau0
# genSeq.Pythia8.Commands += ["10000104:all = qb4p5 qb4p5bar 2 13.5 0 900"]
genSeq.Pythia8.Commands += ["Merging:Process = guess"]
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
genSeq.Pythia8.UserHooks += ['JetMergingaMCatNLO']
else:
genSeq.Pythia8.UserHook = 'JetMergingaMCatNLO'
......@@ -4,6 +4,12 @@ masses['1000006'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000023'] = float(runArgs.jobConfig[0].split('_')[5])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
extras['event_norm']='sum'
extras['use_syst']='F'
extras['pdgs_for_merging_cut']='1, 2, 3, 4, 21, 1000001, 1000002, 1000003, 1000004, 1000021, 2000001, 2000002, 2000003, 2000004'
process = '''
generate p p > t1 t1~ $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
add process p p > t1 t1~ j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
......
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000021'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t~ t1 go $ susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
'''
njets = 2
evgenLog.info('Generation of stop anti-top gluino production, stop to bffN; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
#jobConfigParts = runArgs.jobConfig[0].split("_")
if 'MadSpin' in runArgs.jobConfig[0].split('_')[-2]:
evgenLog.info('Running w/ MadSpin option')
madspin_card='madspin_card_test.dat'
mscard = open(madspin_card,'w')
decay_chains = ["t1 > n1 fu fd~ b", "t~ > b~ W-, W- > ll- vl~, W- > qu~ qd" ]
# if '2L15' in runArgs.jobConfig[0].split('_')[-1] or 'LepBRonly' in runArgs.jobConfig[0].split('_'):
# evgenLog.info('Setting only leptonic decay chains')
# decay_chains = [ "t1 > n1 ll+ vl b" , "t1~ > n1 ll- vl~ b~" ]
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
set spinmode none
# specify the decay for the final state particles
define ll- = e- mu- ta-
define vl~ = ve~ vm~ vt~
deifine qu~ = ~u c~
define qd = d s
decay %s
decay %s
# running the actual code
launch"""%(runArgs.randomSeed,decay_chains[0],decay_chains[1]))
mscard.close()
if 'MET100' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('MET100 filter is applied')
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 100*GeV
evt_multiplier = 100.
elif 'M1501L20' in runArgs.jobConfig[0].split('_')[-1]:
evgenLog.info('1Lepton or MET100 filter is applied')
include ( 'MC15JobOptions/MultiElecMuTauFilter.py' )
filtSeq.MultiElecMuTauFilter.MinPt = 20000.
filtSeq.MultiElecMuTauFilter.MaxEta = 2.8
filtSeq.MultiElecMuTauFilter.NLeptons = 1
filtSeq.MultiElecMuTauFilter.IncludeHadTaus = 0
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 150000.
filtSeq.Expression = "MultiElecMuTauFilter or MissingEtFilter"
evt_multiplier = 5.
evgenConfig.contact = [ "takashi.yamanaka@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop']
evgenConfig.description = 'top stop gluino production, st->bffN in simplified model, m_guino = %s GeV, m_stop = %s GeV'%(masses['1000021'],masses['1000006'])
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
# genSeq.Pythia8.Commands += ["Merging:Process = pp>{t,6}{t1~,-1000006}{go,1000021}"]
evgenLog.info('Performing the Pythia guess CKKW-L merging feature')
# id:all = name antiName spinType chargeType colType m0 mWidth mMin mMax tau0
# genSeq.Pythia8.Commands += ["10000104:all = qb4p5 qb4p5bar 2 13.5 0 900"]
genSeq.Pythia8.Commands += ["Merging:Process = guess"]
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
genSeq.Pythia8.UserHooks += ['JetMergingaMCatNLO']
else:
genSeq.Pythia8.UserHook = 'JetMergingaMCatNLO'
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
masses['1000006'] = float(runArgs.jobConfig[0].split('_')[4])
masses['1000022'] = float(runArgs.jobConfig[0].split('_')[5].split('.')[0])
gentype = str(runArgs.jobConfig[0].split('_')[2])
decaytype = str(runArgs.jobConfig[0].split('_')[3])
process = '''
generate p p > t~ t1 n1 $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
'''
njets = 2
evgenLog.info('Registered generation of antitop stop neutralino production, stop to bffN; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000006']))
#jobConfigParts = runArgs.jobConfig[0].split("_")
if 'MadSpin' in runArgs.jobConfig[0].split('_')[-2]:
evgenLog.info('Running w/ MadSpin option')
madspin_card='madspin_card_test.dat'
mscard = open(madspin_card,'w')
decay_chains = [ "t1 > n1 fu fd~ b" , " t~ > b~ W-, W- > ll- vl~, W- > qu~ qd"]
mscard.write("""#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
#
set seed %i
set spinmode none
# specify the decay for the final state particles
define ll+ = e+ mu+
define vl = ve vm
deifine qu~ = ~u c~
define qd = d s
decay %s
decay %s
# running the actual code
launch"""%(runArgs.randomSeed,decay_chains[0],decay_chains[1]))
mscard.close()
if 'MET100' in runArgs.jobConfig[0].split("_")[-1]:
evgenLog.info('MET100 filter is applied')
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 100*GeV
evt_multiplier = 50.
elif 'M1501L20' in runArgs.jobConfig[0].split('_')[-1]:
evgenLog.info('1Lepton or MET100 filter is applied')
include ( 'MC15JobOptions/MultiElecMuTauFilter.py' )
filtSeq.MultiElecMuTauFilter.MinPt = 20000.
filtSeq.MultiElecMuTauFilter.MaxEta = 2.8
filtSeq.MultiElecMuTauFilter.NLeptons = 1
filtSeq.MultiElecMuTauFilter.IncludeHadTaus = 0
include ( 'MC15JobOptions/MissingEtFilter.py' )
filtSeq.MissingEtFilter.METCut = 150000.
filtSeq.Expression = "MultiElecMuTauFilter or MissingEtFilter"
evt_multiplier = 5.
evgenConfig.contact = [ "paola.arrubarrena@cern.ch" ]
evgenConfig.keywords += ['simplifiedModel', 'stop']
evgenConfig.description = 'top stop neutralino production, st->bffN in simplified model, m_stop = %s GeV, m_neutralino = %s GeV'%(masses['1000006'],masses['1000022'])
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
if njets>0:
# genSeq.Pythia8.Commands += ["Merging:Process = pp>{t,6}{t1~,-1000006}{go,1000021}"]
evgenLog.info('Performing the Pythia guess CKKW-L merging feature')
# id:all = name antiName spinType chargeType colType m0 mWidth mMin mMax tau0
# genSeq.Pythia8.Commands += ["10000104:all = qb4p5 qb4p5bar 2 13.5 0 900"]
genSeq.Pythia8.Commands += ["Merging:Process = guess"]
if "UserHooks" in genSeq.Pythia8.__slots__.keys():
genSeq.Pythia8.UserHooks += ['JetMergingaMCatNLO']
else:
genSeq.Pythia8.UserHook = 'JetMergingaMCatNLO'
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
# JO for Pythia 8 jet jet with large-R jet filtering
evgenConfig.description = "Dijet truth jet filtered with the A14 NNPDF23 LO tune"
evgenConfig.keywords = ["QCD", "jets", "SM"]
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
genSeq.Pythia8.Commands += ["HardQCD:all = on",
"PhaseSpace:pTHatMin = 250."]
## Truth jet filter common config for all JZx and JZxW
include("JetFilterAkt10.py")
filtSeq.QCDTruthJetFilter.MinPt = 400*GeV
filtSeq.QCDTruthJetFilter.MaxPt = 7000*GeV
evgenConfig.minevents = 500
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
include( 'MC15JobOptions/MadGraphControl_SimplifiedModel_ATGT_bffN.py' )
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