Commit 7ba1d18b authored by mcjo's avatar mcjo
Browse files

add 410545-6, update 345672,4, add 375075-115, remove 385075-115


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@810910 4525493e-7705-40b1-a816-d608a930855b
parent 7565d137
2017-10-06 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-07-96
* add 410545-6, update 345672,4, add 375075-115, remove 385075-115
* update MadGraphControl_aMcAtNlo_Wt_DR2.py
* update MadGraphControl_SimplifiedModel_HiggsinoLLP.py
* update param_card.SM.C1N2.HiggsinoLLP.dat
2017-10-05 Ewelina Lobodzinska <ewelina.maria.lobodzinska@cern.ch>
* add 361281-9
* add 361281-9
2017-10-03 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-00-07-95
......
......@@ -49,7 +49,7 @@ else:
raise RunTimeError("ERROR: deltaM and lifetime are not correctly set")
masses['1000024'] = masses['1000022'] + deltaM
masses['1000023'] = masses['1000022']
masses['1000023'] = masses['1000022'] * (-1)
if str(splitConfig[5]) == "0p1":
lifetime = 0.1
......@@ -62,7 +62,6 @@ elif str(splitConfig[5]) == "0p3":
# process lines
# ----------------------------------------------
mergeproc = ""
if gentype == "C1pN1":
process = '''
generate p p > x1+ n1 $ susystrong @1
......@@ -70,7 +69,7 @@ add process p p > x1+ n1 j $ susystrong @2
add process p p > x1+ n1 j j $ susystrong @3
'''
mergeproc = "{x1+,1000024}{n1,1000022}"
if gentype == "C1mN1":
elif gentype == "C1mN1":
process = '''
generate p p > x1- n1 $ susystrong @1
add process p p > x1- n1 j $ susystrong @2
......@@ -83,7 +82,7 @@ generate p p > x1+ x1- $ susystrong @1
add process p p > x1+ x1- j $ susystrong @2
add process p p > x1+ x1- j j $ susystrong @3
'''
mergeproc = "{x1+,1000024}{x1-,1000024}"
mergeproc = "{x1+,1000024}{x1-,-1000024}"
njets = 2
# ----------------------------------------------
......@@ -108,25 +107,26 @@ EleBR = 0.030
MuBR = 0.015
evgenLog.info('Registered generation of ~chi1+/- ~chi20 production (long-lived higgsinos)'+str(runArgs.runNumber))
evgenConfig.contact = ["keisuke.yoshihara@cern.ch"]
evgenConfig.keywords += ['simplifiedModel','gaugino', 'neutralino', 'chargino', 'longLived']
evgenConfig.description = '~chi1+/- ~chi20 production (long-lived higgsinos), m_C1 = %s GeV, m_N1 = %s GeV, lifetime = %s ns' %(masses['1000024'],masses['1000022'], lifetime)
gentype = "C1N2" # dummy gentype to set xqcut correctly
include ('MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py')
if lifetime != 0:
evgenConfig.specialConfig = 'AMSBC1Mass=%s*GeV;AMSBN1Mass=%s*GeV;AMSBC1ToEleBR=%s;AMSBC1ToMuBR=%s;AMSBC1Lifetime=%s*ns;preInclude=SimulationJobOptions/preInclude.AMSB.py' % (masses['1000024'],masses['1000022'],EleBR,MuBR,lifetime)
gentype = "C1N2" # dummy gentype to set xqcut correctly
include ('MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py')
# ----------------------------------------------
# pythia8 commands etc
# ----------------------------------------------
genSeq.Pythia8.Commands += ["1000022:all = Chi Chi~ 1 0 0 %s" % (masses['1000022'])]
if njets>0:
genSeq.Pythia8.Commands += [ "Merging:Process = pp>%s" % mergeproc,
"1000022:spinType = 1",
"1000023:spinType = 1",
"1000024:spinType = 1" ]
bonus_file = open('pdg_extras.dat','w')
......
......@@ -255,19 +255,6 @@ if parton_shower=='PYTHIA8':
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_aMcAtNlo.py")
genSeq.Pythia8.Commands += ["SpaceShower:pTmaxMatch = 1",
"SpaceShower:pTmaxFudge = 1",
"SpaceShower:MEcorrections = off",
"TimeShower:pTmaxMatch = 1",
"TimeShower:pTmaxFudge = 1",
"TimeShower:MEcorrections = off",
"TimeShower:globalRecoil = on",
"TimeShower:limitPTmaxGlobal = on",
"TimeShower:nMaxGlobalRecoil = 1",
"TimeShower:globalRecoilMode = 2",
"TimeShower:nMaxGlobalBranch = 1.",
"Check:epTolErr = 1e-2" ]
elif parton_shower=="HERWIGPP":
## run Herwig7 on-the-fly -----------------------------------------------------
## Provide config information
......
......@@ -2,7 +2,7 @@
# Powheg ttH setup starting from ATLAS defaults
#--------------------------------------------------------------
include('PowhegControl/PowhegControl_ttH_Common.py')
PowhegConfig.topdecaymode = 22222 #ttbar we use 22222 which turns on all decay modes
PowhegConfig.topdecaymode = 00022 #ttbar we use 22222 which turns on all decay modes
PowhegConfig.hdamp = 352.5
PowhegConfig.runningscales = 1 ## dynamic scale
......@@ -10,7 +10,7 @@ PowhegConfig.storeinfo_rwgt = 1
PowhegConfig.PDF = range(260000, 260101) + range(90400, 90433) + [11068] + [25200] + [13165] + [25300] # PDF variations
PowhegConfig.mu_F = [ 1.0, 1.0, 1.0, 0.5, 0.5, 0.5, 2.0, 2.0, 2.0 ] # scale variations: first pair is the nominal setting
PowhegConfig.mu_R = [ 1.0, 0.5, 2.0, 1.0, 0.5, 2.0, 1.0, 0.5, 2.0 ]
PowhegConfig.nEvents *= 3. #to compensate for efficiency filter
PowhegConfig.nEvents *= 1.1 #safety factor
PowhegConfig.generate()
#--------------------------------------------------------------
......@@ -27,9 +27,6 @@ genSeq.Pythia8.UserModes += [ 'Main31:pTemt = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:emitted = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:MPIveto = 0' ]
include("MC15JobOptions/Pythia8_SMHiggs125_inc.py")
include("MC15JobOptions/TTbarWToLeptonFilter.py")
filtSeq.TTbarWToLeptonFilter.NumLeptons = 0 #(-1: non-all had, 0: all had, 1: l+jets, 2: dilepton)
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
#--------------------------------------------------------------
# EVGEN configuration
......
......@@ -2,7 +2,7 @@
# Powheg ttH setup starting from ATLAS defaults
#--------------------------------------------------------------
include('PowhegControl/PowhegControl_ttH_Common.py')
PowhegConfig.topdecaymode = 22222 #ttbar we use 22222 which turns on all decay modes
PowhegConfig.topdecaymode = 22200 #ttbar we use 22222 which turns on all decay modes
PowhegConfig.hdamp = 352.5
PowhegConfig.runningscales = 1 ## dynamic scale
......@@ -10,7 +10,7 @@ PowhegConfig.storeinfo_rwgt = 1
PowhegConfig.PDF = range(260000, 260101) + range(90400, 90433) + [11068] + [25200] + [13165] + [25300] # PDF variations
PowhegConfig.mu_F = [ 1.0, 1.0, 1.0, 0.5, 0.5, 0.5, 2.0, 2.0, 2.0 ] # scale variations: first pair is the nominal setting
PowhegConfig.mu_R = [ 1.0, 0.5, 2.0, 1.0, 0.5, 2.0, 1.0, 0.5, 2.0 ]
PowhegConfig.nEvents *= 3. #to compensate for efficiency filter
PowhegConfig.nEvents *= 1.1 #to compensate for efficiency filter
PowhegConfig.generate()
#--------------------------------------------------------------
......@@ -27,9 +27,6 @@ genSeq.Pythia8.UserModes += [ 'Main31:pTemt = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:emitted = 0' ]
genSeq.Pythia8.UserModes += [ 'Main31:MPIveto = 0' ]
include("MC15JobOptions/Pythia8_SMHiggs125_inc.py")
include("MC15JobOptions/TTbarWToLeptonFilter.py")
filtSeq.TTbarWToLeptonFilter.NumLeptons = 2 #(-1: non-all had, 0: all had, 1: l+jets, 2: dilepton)
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.0
#--------------------------------------------------------------
# EVGEN configuration
......
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