Commit 93518829 authored by Joany Manjarres's avatar Joany Manjarres
Browse files

See ChangeLog

parent 57419629
2020-04-15 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-03-24/
* add JO 451706-451707 451700-451703
* add common MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py
2020-04-09 Cyril Becot <cbecot@cern.ch>
* tag MC15JobOptions-01-03-23/
* up share/evgeninputfiles.csv
* add 313352-313371
2020-04-06 Joany Manjarres <joany@cern.ch>
2020-04-06 Joany Manjarres <joany@cern.ch>
* tag MC15JobOptions-01-03-22
* add MadGraphControl_SimplifiedModel_BB_direct_upperMET.py
......
from MadGraphControl.MadGraphUtils import *
import math
fcard = open('proc_card_mg5.dat','w')
# generate ... QED=0 QCD=3
fcard.write("""
import model DMsimp_s_spin1 -modelname
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define all = 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+ h w- ta- ta+
""")
if "ee" in runArgs.jobConfig[0]:
fcard.write("""
generate p p > Y1 > e+ e- j, Y1 > e+ e-
""")
elif "mumu" in runArgs.jobConfig[0]:
#generate p p > Y1 > mu+ mu-
fcard.write("""
generate p p > Y1 > mu+ mu-
add process p p > Y1 all, (Y1 > mu+ mu-)
display diagrams
""")
else:
raise RuntimeError("No dilepton channel specified.")
fcard.write("""
output -f
""")
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()
#Fetch default LO run_card.dat and set parameters
extras = {'lhe_version':'2.0',
'cut_decays' :'F',
'pdlabel' : "'lhapdf'",
'lhaid' : 263000,
'ickkw' : 1,
'ptj' : 23,
'pta' : 0,
'ptl' : 4,
'ptl1min' : 11,
#'xptl' : 13,
'etaj' : -1,
'etaa' :-1,
'etal' :3.0,
'drll' :0.0,
'drjj' :0.0,
'draa' :0.0,
'draj' :0.0,
'drja' :0.0,
'dral' :0.0,
'drll' : 0
}
#'xptj' : 23,
#'etaj' : 5 }
#from https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MadGraph5aMCatNLOForAtlas#Problems_with_run_card_dat_in_ne
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat', xqcut=10,
nevts=runArgs.maxEvents*2,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=extras)
print_cards()
paramcard = subprocess.Popen(['get_files','-data','MadGraph_param_card_DMsimp_s_spin1.dat'])
paramcard.wait()
if not os.access('MadGraph_param_card_DMsimp_s_spin1.dat',os.R_OK):
print 'ERROR: Could not get param card'
elif os.access('param_card.dat',os.R_OK):
print 'ERROR: Old param card in the current directory. Dont want to clobber it. Please move it first.'
else:
oldcard = open('MadGraph_param_card_DMsimp_s_spin1.dat','r')
newcard = open('param_card.dat','w')
for line in oldcard:
if '# gVXc' in line:
newcard.write(' 1 %e # gVXc \n'%(gVXc))
elif '# gVXd' in line:
newcard.write(' 2 %e # gVXd \n'%(gVXd))
elif '# gAXd' in line:
newcard.write(' 3 %e # gAXd \n'%(gAXd))
elif '# gVd11' in line:
newcard.write(' 4 %e # gVd11 \n'%(gVd11))
elif '# gVu11' in line:
newcard.write(' 5 %e # gVu11 \n'%(gVu11))
elif '# gVd22' in line:
newcard.write(' 6 %e # gVd22 \n'%(gVd22))
elif '# gVu22' in line:
newcard.write(' 7 %e # gVu22 \n'%(gVu22))
elif '# gVd33' in line:
newcard.write(' 8 %e # gVd33 \n'%(gVd33))
elif '# gVu33' in line:
newcard.write(' 9 %e # gVu33 \n'%(gVu33))
elif '# gVl11' in line:
newcard.write(' 10 %e # gVl11 \n'%(gVl11))
elif '# gVl22' in line:
newcard.write(' 11 %e # gVl22 \n'%(gVl22))
elif '# gAd11' in line:
newcard.write(' 13 %e # gAd11 \n'%(gAd11))
elif '# gAu11' in line:
newcard.write(' 14 %e # gAu11 \n'%(gAu11))
elif '# gAd22' in line:
newcard.write(' 15 %e # gAd22 \n'%(gAd22))
elif '# gAu22' in line:
newcard.write(' 16 %e # gAu22 \n'%(gAu22))
elif '# gAd33' in line:
newcard.write(' 17 %e # gAd33 \n'%(gAd33))
elif '# gAu33' in line:
newcard.write(' 18 %e # gAu33 \n'%(gAu33))
elif '# gAl11' in line:
newcard.write(' 19 %e # gAl11 \n'%(gAl11))
elif '# gAl22' in line:
newcard.write(' 20 %e # gAl22 \n'%(gAl22))
elif '# gAl33' in line:
newcard.write(' 21 %e # gAl22 \n'%(gAl33))
elif ' MY1 ' in line:
newcard.write(' 5000001 %e # MY1 \n'%(MY1))
elif 'DECAY 5000001' in line :
# do I need to change this to the value calculated from WY1
newcard.write('DECAY 5000001 auto #WY1 \n')
elif ' xd : MXd ' in line:
newcard.write(' 1000022 %e # xd : MXd \n'%(MXd))
elif ' # MXd' in line:
newcard.write(' 1000022 %e # MXd \n'%(MXd))
elif ' # MXr' in line:
newcard.write(' 5000511 %e # MXr \n'%(MXr))
elif ' # MXc' in line:
newcard.write(' 5000512 %e # MXc \n'%(MXc))
else:
newcard.write(line)
oldcard.close()
newcard.close()
runName='run_01'
generate(run_card_loc='run_card.dat',param_card_loc='param_card.dat',mode=0,njobs=1,run_name=runName,proc_dir=process_dir)
arrange_output(run_name=runName,proc_dir=process_dir,outputDS=runName+'._00001.events.tar.gz')
#### Shower
#evgenConfig.description = "Wimp dmA mediator from DMSimp, ptj>"+str(xptj)+" GeV"
evgenConfig.description = "Wimp dmA mediator from DMSimp with ISR"
evgenConfig.keywords = ["exotic","BSM"]
evgenConfig.process = "pp > Y1 > ll "
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
evgenConfig.contact = ["Marie-Helene Genest <mgenest@cern.ch>"]
evgenConfig.minevents= 10000
include("MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("MC15JobOptions/Pythia8_MadGraph.py")
#include("MC15JobOptions/Pythia8_aMcAtNlo.py")
#particle data = name antiname spin=2s+1 3xcharge colour mass width (left out, so set to 0: mMin mMax tau0)
genSeq.Pythia8.Commands += ["1000022:all = xd xd~ 2 0 0 %d 0" %(MXd),
"1000022:isVisible = false"]
MXd = 10000
MY1 = 15
gVl33 = 0.0
gVu33 = 0.0
gAXd = 0.0
gVd33 = 0.0
gVl11 = 0.0
gAl22 = 0.1
gVu11 = 0.0
gAu22 = 0.1
gAd22 = 0.1
gVl22 = 0.0
gVu22 = 0.0
gVd22 = 0.0
MXr = 10000
gAl11 = 0.1
MXc = 10000
gAl33 = 0.1
gAd33 = 0.1
gAu11 = 0.1
gAu33 = 0.1
gVXd = 0.0
gVXc = 0.0
gVd11 = 0.0
gAd11 = 0.1
import os
include("MC15JobOptions/MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py")
MXd = 10000
MY1 = 20
gVl33 = 0.0
gVu33 = 0.0
gAXd = 0.0
gVd33 = 0.0
gVl11 = 0.0
gAl22 = 0.1
gVu11 = 0.0
gAu22 = 0.1
gAd22 = 0.1
gVl22 = 0.0
gVu22 = 0.0
gVd22 = 0.0
MXr = 10000
gAl11 = 0.1
MXc = 10000
gAl33 = 0.1
gAd33 = 0.1
gAu11 = 0.1
gAu33 = 0.1
gVXd = 0.0
gVXc = 0.0
gVd11 = 0.0
gAd11 = 0.1
import os
include("MC15JobOptions/MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py")
MXd = 10000
MY1 = 25
gVl33 = 0.0
gVu33 = 0.0
gAXd = 0.0
gVd33 = 0.0
gVl11 = 0.0
gAl22 = 0.1
gVu11 = 0.0
gAu22 = 0.1
gAd22 = 0.1
gVl22 = 0.0
gVu22 = 0.0
gVd22 = 0.0
MXr = 10000
gAl11 = 0.1
MXc = 10000
gAl33 = 0.1
gAd33 = 0.1
gAu11 = 0.1
gAu33 = 0.1
gVXd = 0.0
gVXc = 0.0
gVd11 = 0.0
gAd11 = 0.1
import os
include("MC15JobOptions/MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py")
MXd = 10000
MY1 = 30
gVl33 = 0.0
gVu33 = 0.0
gAXd = 0.0
gVd33 = 0.0
gVl11 = 0.0
gAl22 = 0.1
gVu11 = 0.0
gAu22 = 0.1
gAd22 = 0.1
gVl22 = 0.0
gVu22 = 0.0
gVd22 = 0.0
MXr = 10000
gAl11 = 0.1
MXc = 10000
gAl33 = 0.1
gAd33 = 0.1
gAu11 = 0.1
gAu33 = 0.1
gVXd = 0.0
gVXc = 0.0
gVd11 = 0.0
gAd11 = 0.1
import os
include("MC15JobOptions/MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py")
MXd = 10000
MY1 = 35
gVl33 = 0.0
gVu33 = 0.0
gAXd = 0.0
gVd33 = 0.0
gVl11 = 0.0
gAl22 = 0.1
gVu11 = 0.0
gAu22 = 0.1
gAd22 = 0.1
gVl22 = 0.0
gVu22 = 0.0
gVd22 = 0.0
MXr = 10000
gAl11 = 0.1
MXc = 10000
gAl33 = 0.1
gAd33 = 0.1
gAu11 = 0.1
gAu33 = 0.1
gVXd = 0.0
gVXc = 0.0
gVd11 = 0.0
gAd11 = 0.1
import os
include("MC15JobOptions/MadGraphControl_MGPy8EG_N30LO_A14N23LO_dmA_oic.py")
include("MC15JobOptions/Sherpa_NNPDF30NNLO_Common.py")
evgenConfig.description = "5FS PP-> bA+X with A-> b ~b with up to 4 jets in ME+PS"
evgenConfig.keywords = [ "BSMHiggs", "bHiggs","bbbar"]
evgenConfig.contact = [ "timb@slac.stanford.edu" ]
evgenConfig.inputconfcheck = "bhbb"
evgenConfig.process="""
(run){
MASS[25]=0950.; WIDTH[25]=0.;
YUKAWA[5] 4.8;
PARTICLE_CONTAINER 98 B 5 -5;
CSS_IS_AS_FAC 1;
HARD_DECAYS On;
HDH_STATUS[25,5,-5]=2;
HDH_BR_WEIGHTS=1;
}(run);
(processes){
Process 93 93 -> 25 98 93{1};
CKKW sqr(20.0/E_CMS);
Order (*,1);
End process;
}(processes);
(selector){
PT 98 1 E_CMS;
}(selector);
"""
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,5,-5]=7.830" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,15,-15]=0.696" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,6,-6]=0.174" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,24,-24]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,23,23]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,13,-13]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,4,-4]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,3,-3]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,21,21]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,22,22]=0." ]
include("MC15JobOptions/Sherpa_NNPDF30NNLO_Common.py")
evgenConfig.description = "5FS PP-> bA+X with A-> b ~b with up to 4 jets in ME+PS"
evgenConfig.keywords = [ "BSMHiggs", "bHiggs","bbbar"]
evgenConfig.contact = [ "timb@slac.stanford.edu" ]
evgenConfig.inputconfcheck = "bhbb"
evgenConfig.process="""
(run){
MASS[25]=1350.; WIDTH[25]=0.;
YUKAWA[5] 4.8;
PARTICLE_CONTAINER 98 B 5 -5;
CSS_IS_AS_FAC 1;
HARD_DECAYS On;
HDH_STATUS[25,5,-5]=2;
HDH_BR_WEIGHTS=1;
}(run);
(processes){
Process 93 93 -> 25 98 93{1};
CKKW sqr(20.0/E_CMS);
Order (*,1);
End process;
}(processes);
(selector){
PT 98 1 E_CMS;
}(selector);
"""
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,5,-5]=10.421" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,15,-15]=0.964" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,6,-6]=0.241" ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,24,-24]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,23,23]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,13,-13]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,4,-4]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,3,-3]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,21,21]=0." ]
genSeq.Sherpa_i.Parameters += [ "HDH_WIDTH[25,22,22]=0." ]
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