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):
print("Deleting old run_card.dat")
os.remove('run_card.dat')
if run_card == "": run_card = MadGraphUtils.get_default_runcard(proc_dir=process_dir)
if MadGraphUtils.build_run_card(run_card_old=run_card,run_card_new='run_card.dat',nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,extras=run_card_extras)==-1:
raise RuntimeError("Could not create run_card.dat")
#ready to generate
MadGraphUtils.generate(run_card_loc='run_card.dat',param_card_loc='param_card.dat',mode=0,njobs=1,run_name='Test',proc_dir=process_dir)
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_'+str(proc_name)
if "post_lhe_hook" in dir():
print("Calling post_lhe_hook....")
post_lhe_hook()
MadGraphUtils.arrange_output(run_name='Test',proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz')
if 'ATHENA_PROC_NUMBER' in os.environ:
njobs = os.environ.pop('ATHENA_PROC_NUMBER')
if not hasattr(opts,'nprocs'): mglog.warning('Did not see option!')
else: opts.nprocs = 0
evgenConfig.inputfilecheck = stringy
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'
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
evgenConfig.contact = ['*****@*****.**']
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
print("Deleting old run_card.dat")
os.remove('run_card.dat')
if run_card == "":
run_card = MadGraphUtils.get_default_runcard(proc_dir=process_dir)
if MadGraphUtils.build_run_card(run_card_old=run_card,
run_card_new='run_card.dat',
nevts=nevents,
rand_seed=runArgs.randomSeed,
beamEnergy=beamEnergy,
extras=run_card_extras) == -1:
raise RuntimeError("Could not create run_card.dat")
#ready to generate
MadGraphUtils.generate(run_card_loc='run_card.dat',
param_card_loc='param_card.dat',
mode=0,
njobs=1,
run_name='Test',
proc_dir=process_dir)
stringy = 'madgraph.' + str(runArgs.runNumber) + '.MadGraph_' + str(proc_name)
if "post_lhe_hook" in dir():
print("Calling post_lhe_hook....")
post_lhe_hook()
MadGraphUtils.arrange_output(run_name='Test',
proc_dir=process_dir,
outputDS=stringy + '._00001.events.tar.gz')
if 'ATHENA_PROC_NUMBER' in os.environ:
njobs = os.environ.pop('ATHENA_PROC_NUMBER')
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'
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile=runName+'._00001.events.tar.gz'
evgenConfig.contact = ['*****@*****.**']
include('MC15JobOptions/Pythia8_A14_NNPDF23LO_EvtGen_Common.py')
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 -> W chi chibar"
evgenConfig.inputfilecheck = runName
runArgs.inputGeneratorFile = runName + "._00001.events.tar.gz"
evgenConfig.contact = ["*****@*****.**"]
