beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found")
#make the proc card
fcard = open('proc_card_mg5.dat','w')
fcard.write(proc_card + """
output %s -f
""" % proc_name)
fcard.close()
#does mg5 to determine all the diagrams etc ... generates the process directory corresponding to proc_name
process_dir = MadGraphUtils.new_process()
if proc_name != "" and process_dir != proc_name:
raise RuntimeError("Unexpected process_dir %s when expected %s" % (process_dir,proc_name))
#create the param_card
if os.access('param_card.dat',os.R_OK):
print("Deleting old param_card.dat")
os.remove('param_card.dat')
if param_card == "": param_card = "%s/Cards/param_card.dat" % process_dir
if( MadGraphUtils.build_param_card(param_card_old=param_card,param_card_new='param_card.dat',params=param_card_extras) == -1):
raise RuntimeError("Could not create param_card.dat")
#create the run card: FIXME: Should check for success
if os.access('run_card.dat',os.R_OK):
elif '# WEtaU' in line:
newcard.write('DECAY 9000007 {0:0.6e} # WEtaU \n'.format(wEtau))
elif '# WEtaS' in line:
newcard.write('DECAY 9000008 {0:0.6e} # WEtaS \n'.format(wEtas))
elif '# WEtaC' in line:
newcard.write('DECAY 9000009 {0:0.6e} # WEtaC \n'.format(wEtac))
elif '# WEtaT' in line:
newcard.write('DECAY 9000011 {0:0.6e} # WEtaT \n'.format(wEtat))
elif '# WEtaB' in line:
newcard.write('DECAY 9000014 {0:0.6e} # WEtaB \n'.format(wEtab))
else:
newcard.write(line)
oldcard.close()
newcard.close()
process_dir = MG.new_process()
runName='run_01'
MG.generate(run_card_loc='run_card.dat', param_card_loc='param_card.dat',
mode=0, njobs=1, run_name=runName, proc_dir=process_dir)
MG.arrange_output(run_name=runName, proc_dir=process_dir,
outputDS=runName+'._00001.events.tar.gz')
evgenConfig.description = ('Wimp pair mono-WZ with simplified model, dmST, '
'M_WIMP={0:d} GeV, M_MED={1:d} GeV').format(int(mDM),
int(mMed))
evgenConfig.keywords = ['WZ','exotic','BSM','WIMP']
evgenConfig.process = 'pp -> Z chi chibar'
elif '# gxiB' in line:
newcard.write(' 6 {0:0.6e} # gxiB \n'.format(gxiU))
elif '# gxiT' in line:
newcard.write(' 7 {0:0.6e} # gxiT \n'.format(gxiU))
elif 'Mchi' in line:
newcard.write(' 1000022 {0:0.6e} # Mchi\n'.format(mDM))
elif 'Mxi' in line:
newcard.write(' 101 {0:0.6e} # Mxi\n'.format(mMed))
elif 'Wxi' in line:
newcard.write('DECAY 101 {0:0.6e} # Wxi\n'.format(wMed))
else:
newcard.write(line)
oldcard.close()
newcard.close()
process_dir = MG.new_process()
runName='run_01'
MG.generate(run_card_loc='run_card.dat', param_card_loc='param_card.dat',
mode=0, njobs=1, run_name=runName, proc_dir=process_dir)
MG.arrange_output(run_name=runName, proc_dir=process_dir,
outputDS=runName+'._00001.events.tar.gz')
evgenConfig.description = ('Wimp pair mono-WZ with simplified model, dmV, '
'M_WIMP={0:d} GeV, M_MED={1:d} GeV').format(int(mDM),
int(mMed))
evgenConfig.keywords = ['WZ','exotic','BSM','WIMP']
evgenConfig.process = 'pp -> W chi chibar'
beamEnergy = -999
if hasattr(runArgs, 'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found")
#make the proc card
fcard = open('proc_card_mg5.dat', 'w')
fcard.write(proc_card + """
output %s -f
""" % proc_name)
fcard.close()
#does mg5 to determine all the diagrams etc ... generates the process directory corresponding to proc_name
process_dir = MadGraphUtils.new_process()
if proc_name != "" and process_dir != proc_name:
raise RuntimeError("Unexpected process_dir %s when expected %s" %
(process_dir, proc_name))
#create the param_card
if os.access('param_card.dat', os.R_OK):
print("Deleting old param_card.dat")
os.remove('param_card.dat')
if param_card == "": param_card = "%s/Cards/param_card.dat" % process_dir
if (MadGraphUtils.build_param_card(param_card_old=param_card,
param_card_new='param_card.dat',
params=param_card_extras) == -1):
raise RuntimeError("Could not create param_card.dat")
