Commit 680a400f authored by MC JO registration's avatar MC JO registration
Browse files

add 450694-6, 450697-8, 312243, 449616-772, 346587, update 346551

parent 911167cb
2019-05-16 Evelina Bouhova-Thacker <e.bouhova@cern.ch>
* tagging MC15JobOptions-01-00-37
* add 450694-6, 450697-8, 312243, 449616-772, 346587, update 346551
2019-06-03 E. M. Lobodzinska <ewelina@mail.desy.de>
* tag MC15JobOptions-01-00-36
* add 361269
* tag MC15JobOptions-01-00-36
* add 361269
2019-05-31 Cyril Becot <cyril.becot@cern.ch>
* up MadGraphControl_SimplifiedModel_C1N2_Wh
......@@ -20,8 +24,8 @@
2019-05-28 E. M. Lobodzinska <ewelina@mail.desy.de>
* tag MC15JobOptions-01-00-33
* correct 363252-3
* tag MC15JobOptions-01-00-33
* correct 363252-3
2019-05-27 Cyril Becot <cyril.becot@cern.ch>
* tag MC15JobOptions-01-00-32/
......@@ -47,8 +51,8 @@
* tag MC15JobOptions-01-00-29
2019-05-23 E. M. Lobodzinska <ewelina@mail.desy.de>
* correct 361292
* tag MC15JobOptions-01-00-28
* correct 361292
* tag MC15JobOptions-01-00-28
2019-05-23 Joany Manjarres <joany@cern.ch>
* tagging MC15JobOptions-01-00-27
......@@ -66,8 +70,8 @@
* update cvs line 413024
2019-05-22 E. M. Lobodzinska <ewelina@mail.desy.de>
* correct 363240-3, 361281-3
* tag MC15JobOptions-01-00-24
* correct 363240-3, 361281-3
* tag MC15JobOptions-01-00-24
2019-05-21 Joany Manjarres <joany@cern.ch>
* tagging MC15JobOptions-01-00-23
......
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
splitConfig = runArgs.jobConfig[0].rstrip('.py').split('_')
#G/N2 degenerate
masses['1000022'] = float(splitConfig[5]) #--chi_10 LSP
masses['1000023'] = float(splitConfig[4]) #--chi_20 NLSP
if masses['1000022']<0.5:
masses['1000022']=0.5
lft = float(splitConfig[6]) #--chi_20 NLSP lifetime
if lft<1:
lft=0.5
width = 6.58212196173e-16/lft
#will be WH125
gentype = splitConfig[2]
#will be N2N2
decaytype = splitConfig[3]
decays['1000023'] = '''
# PDG Width
DECAY 1000023 ''' +str(width)+ ''' # neutralino2 decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 22 # BR(~chi_2 -> gamma ~cHi-_1 )'''
process = '''
generate p p > h01 w+ , (w+ > l+ vl) @1
add process p p > h01 w- , (w- > l- vl~) @1
'''
njets = 0
evgenLog.info('Registered generation of h01 w production, decay via h > n2 n2, w > l vl; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000023']) + ' ' + str(masses['1000022']) )
evgenConfig.contact = [ "haichen.wang@cern.ch" ]
evgenConfig.keywords += ['W', 'neutralino', 'lepton','longLived']
evgenConfig.description = 'h01 w production, decay via n2, n2 in simplified model, m_N1 = %s GeV, m_N2 = %s GeV, tau = %s ns'%(masses['1000022'],masses['1000023'],lft)
#--------------------------------------------------------------
# add some filter here
#--------------------------------------------------------------
# need more events from MG due to filter - this needs to be set before MadGraphControl_SimplifiedModelPostInclude.py is run (it's set at 2 there)
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
testSeq.TestHepMC.MaxVtxDisp = 10000*10000 #in mm
testSeq.TestHepMC.MaxTransVtxDisp = 10000*10000
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
splitConfig = runArgs.jobConfig[0].rstrip('.py').split('_')
#G/N2 degenerate
masses['1000022'] = float(splitConfig[5]) #--chi_10 LSP
masses['1000023'] = float(splitConfig[4]) #--chi_20 NLSP
lft = float(splitConfig[6]) #--chi_20 NLSP lifetime
if lft<1:
lft=0.5
width = 6.58212196173e-16/lft
if masses['1000022']<0.5:
masses['1000022']=0.5
#will be ZH125
gentype = splitConfig[2]
#will be N2N2
decaytype = splitConfig[3]
decays['1000023'] = '''
# PDG Width
DECAY 1000023 ''' +str(width)+ ''' # neutralino2 decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 22 # BR(~chi_2 -> gamma ~cHi-_1 )'''
process = '''
generate p p > h01 z , (z > l+ l-) @1
'''
njets = 0
evgenLog.info('Registered generation of h01 z production, decay via h01 > n2 n2, z > l l; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000023']) + ' ' + str(masses['1000022']) )
evgenConfig.contact = [ "sai.neha.santpur@cern.ch" ]
evgenConfig.keywords += ['Z', 'neutralino', 'lepton','longLived']
evgenConfig.description = 'h01 z production, decay via n2, n2 in simplified model, m_N1 = %s GeV, m_N2 = %s GeV, tau = %s ns'%(masses['1000022'],masses['1000023'],lft)
#--------------------------------------------------------------
# add some filter here
#--------------------------------------------------------------
# need more events from MG due to filter - this needs to be set before MadGraphControl_SimplifiedModelPostInclude.py is run (it's set at 2 there)
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
testSeq.TestHepMC.MaxVtxDisp = 10000*10000 #in mm
testSeq.TestHepMC.MaxTransVtxDisp = 10000*10000
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPreInclude.py' )
splitConfig = runArgs.jobConfig[0].rstrip('.py').split('_')
#G/N2 degenerate
masses['1000022'] = float(splitConfig[5]) #--chi_10 LSP
masses['1000023'] = float(splitConfig[4]) #--chi_20 NLSP
lft = float(splitConfig[6]) #--chi_20 NLSP lifetime
if lft<1:
lft=0.5
width = 6.58212196173e-16/lft
if masses['1000022']<0.5:
masses['1000022']=0.5
#will be ttH125
gentype = splitConfig[2]
#will be N2N2
decaytype = splitConfig[3]
decays['1000023'] = '''
# PDG Width
DECAY 1000023 ''' +str(width)+ ''' # neutralino2 decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 22 # BR(~chi_2 -> gamma ~cHi-_1 )'''
process = '''
define allsm = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h01 w- ta- ta+
generate p p > h01 t t~ , ( t > w+ b, w+ > allsm allsm), ( t~ > w- b~, w- > allsm allsm) @1
'''
njets = 0
evgenLog.info('Registered generation of h01 t tbar production, decay via h01 > n2 n2; grid point '+str(runArgs.runNumber)+' decoded into mass point ' + str(masses['1000023']) + ' ' + str(masses['1000022']) )
evgenConfig.contact = [ "sai.neha.santpur@cern.ch" ]
evgenConfig.keywords += ['Z', 'neutralino', 'longLived']
evgenConfig.description = 'h01 t tbar production, decay via n2, n2 in simplified model, m_N1 = %s GeV, m_N2 = %s GeV, tau = %s ns'%(masses['1000022'],masses['1000023'],lft)
#--------------------------------------------------------------
# add some filter here
#--------------------------------------------------------------
include('MC15JobOptions/TTbarWToLeptonFilter.py')
filtSeq.TTbarWToLeptonFilter.NumLeptons = -1 # -1 for non-all had, 0 for all had, 1 for 1lep, 2 for 2lep
filtSeq.TTbarWToLeptonFilter.Ptcut = 0.
# need more events from MG due to filter - this needs to be set before MadGraphControl_SimplifiedModelPostInclude.py is run (it's set at 2 there)
include ( 'MC15JobOptions/MadGraphControl_SimplifiedModelPostInclude.py' )
testSeq.TestHepMC.MaxVtxDisp = 10000*10000 #in mm
testSeq.TestHepMC.MaxTransVtxDisp = 10000*10000
from MadGraphControl.MadGraphUtils import *
import fileinput
import shutil
import subprocess
import os
mode=0
nJobs=1
gridpack_mode=True
gridpack_dir='madevent/'
cluster_type=None
cluster_queue=None
cluster_nb_retry=5
if 'ATHENA_PROC_NUMBER' in os.environ:
print 'Noticed that you have run with an athena MP-like whole-node setup. Will re-configure now to make sure that the remainder of the job runs serially.'
njobs = os.environ.pop('ATHENA_PROC_NUMBER')
# Try to modify the opts underfoot
if not hasattr(opts,'nprocs'): print 'Did not see option!'
else: opts.nprocs = njobs
nJobs=int(njobs)
mode=2
print opts
name = runArgs.jobConfig[0]
thisDSID = int(name.split(".")[1])
runName = str(thisDSID)
isHInv = [312243]
isPythia = [312243]
# ----------------------------------------------
# Some global production settings
# ----------------------------------------------
# Make some excess events to allow for Pythia8 failures
nevents=1.2*runArgs.maxEvents if runArgs.maxEvents>0 else 5500
if(thisDSID in isHInv):
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
set complex_mass_scheme
import model loop_qcd_qed_sm_Gmu-atlas
define p = g u c d s u~ c~ d~ s~ b b~
define j = g u c d s u~ c~ d~ s~ b b~
generate p p > h a j j [QCD]
output -f""")
fcard.close()
else:
print "Unknown DSID"
raise RuntimeError("DSID not known")
name = "hajj"
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(name)
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="madevent/")
#os.system("cp FKS_params.dat "+process_dir+"/Cards/")
#os.system("cp setscales.f "+process_dir+"/SubProcesses/setscales.f")
maxjetflavor = 5
parton_shower='PYTHIA8'
extras = { 'lhe_version' :'3.0',
'pdlabel' : "'lhapdf'",
'lhaid' : 260000, # maybe we use 303400?
'maxjetflavor' : maxjetflavor,
'parton_shower' : "'"+parton_shower+"'",
'jetalgo' : -1,
'jetradius' : 0.4,
'ptj' : 15.0,
'etaj' : -1,
'ickkw' : 0,
'ptgmin' : 10.0,
'epsgamma' :'0.1',
'R0gamma' :'0.1',
'xn' :'2',
'isoEM' :'True',
'etagamma' :'3.0',
'bwcutoff' :'15',
'store_rwgt_info' : '.true.',
'event_norm' : 'sum',
'reweight_scale': '.true.',
'reweight_PDF' : '.true.',
'dynamical_scale_choice': 10
}
madspin_card_loc=None
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,extras=extras)
print_cards()
generate(run_card_loc='run_card.dat',param_card_loc=None,mode=mode,njobs=nJobs,proc_dir=process_dir,run_name=runName,madspin_card_loc=madspin_card_loc,grid_pack=gridpack_mode,gridpack_dir=gridpack_dir,random_seed=runArgs.randomSeed,required_accuracy=0.001,nevents=nevents)
outputDS=arrange_output(run_name=runName, proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz',lhe_version=3,saveProcDir=True)
runArgs.inputGeneratorFile=outputDS
if (thisDSID in isPythia):
#---------------------------------------------------------------------------------------------------
# Pythia8 Showering with A14_NNPDF23LO
#---------------------------------------------------------------------------------------------------
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
#include("MC15JobOptions/Pythia8_AZNLO_CTEQ6L1_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_aMcAtNlo.py")
#--------------------------------------------------------------
# Higgs at Pythia8
#--------------------------------------------------------------
#genSeq.Pythia8.UserModes += [ 'Main31:NFinal = 4' ]
genSeq.Pythia8.Commands += [ '25:onMode = off', # decay of Higgs
'25:onIfMatch = 23 23',
'23:onMode = off', # decay of Z
'23:mMin = 2.0',
'23:onIfMatch = 12 12',
'23:onIfMatch = 14 14',
'23:onIfMatch = 16 16',
]
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
from MadGraphControl.MadGraphUtils import *
evgenConfig.generators = ["MadGraph", "Pythia8", "EvtGen"]
evgenConfig.description = 'MadGraph_ewkhjj_EW6_NLO'
evgenConfig.keywords = ['SM','diboson','VBS','2jet','NLO']
evgenConfig.contact = ["schae@cern.ch","mdshapiro@lbl.gov"]
include("MC15JobOptions/MadGraph_NNPDF30_HInv_NLO.py")
evgenConfig.inputconfcheck="MGaMcAtNloPy8EG_HInv_NLO"
......@@ -6,7 +6,7 @@ evgenConfig.keywords = [ 'SM', 'top', 'ttbar', 'lepton']
evgenConfig.contact = [ 'Christian.Schmitt@cern.ch' ]
evgenConfig.minevents = 1000
filterFactor = 40
filterFactor = 45
if runArgs.trfSubstepName == 'generate' :
......
include("MC15JobOptions/Sherpa_2.2.5_NNPDF30NNLO_Common.py")
evgenConfig.description = "Sherpa bbb SM background"
evgenConfig.keywords = [ "bbbar", "jets"]
evgenConfig.contact = [ 'matthew.h.klein@cern.ch' ]
evgenConfig.inputconfcheck = "bb_MassiveCB_2Bjets_DiPt150"
evgenConfig.minevents = 500
evgenConfig.process="""
(run){
MASSIVE[5] 1;
MASSIVE[4] 1;
}(run);
(me){
EVENT_GENERATION_MODE = PartiallyUnweighted
}(me)
(processes){
# 1 the class anything -> b
Process 93 5 -> 93 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 5 93 -> 93 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 93 -5 -> 93 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 93 -> 93 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 5 93 -> 5 4 -4 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 93 5 -> 5 4 -4 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 93 -> -5 4 -4 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 93 -5 -> -5 4 -4 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
# 2 the class anything -> b b~
Process 93 93 -> 5 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 5 -5 -> 5 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 5 -> 5 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
# 3 the class anything -> b b
Process 5 5 -> 5 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 -5 -> -5 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
# 4 the class anything -> b b b~
Process 93 5 -> 5 -5 5 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 5 93 -> 5 -5 5 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 93 -5 -> 5 -5 -5 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -5 93 -> 5 -5 -5 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
# 5 the class anything -> b b b~ b~
Process 93 93 -> 5 -5 5 -5;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 5 -5 -> 5 -5 5 -5;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -5 5 -> 5 -5 5 -5;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
# 6 the class anything -> b b b b~
Process 5 5 -> 5 5 5 -5;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -5 -5 -> -5 -5 -5 5;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
# 7 the class c+x -> b
Process 4 5 -> 4 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 5 4 -> 4 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 4 5 -> 4 5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 5 4 -> 4 5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 4 -5 -> 4 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 4 -> 4 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 4 -5 -> 4 -5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -5 4 -> 4 -5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
# 8 the class c~+x -> b
Process -4 5 -> -4 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -4 5 -> -4 5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -4 -5 -> -4 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -4 -5 -> -4 -5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process 5 -4 -> -4 5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 5 -4 -> -4 5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
Process -5 -4 -> -4 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process -5 -4 -> -4 -5 4 -4;
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;
Order (*,0);
End process;
# 9 the class charm+x -> b b~
Process 4 -4 -> 5 -5 93{2};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS);
Order (*,0);
End process;
Process 4 93 -> 4 5 -5 93{1};
Max_Epsilon 0.1
CKKW sqr(20/E_CMS); Cut_Core 1;