Page MenuHomeHEPForge
Files F8309584

BatchSubmit_NSO.py
No OneTemporary
Actions

BatchSubmit_NSO.py
View Options

#!/usr/bin/python
import sys
import os
import string
import tarfile
import random
import errno
import time
import subprocess
import numpy as np
def mkdir_p(path):
try:
os.makedirs(path)
except OSError as exc: # Python >2.5
if exc.errno == errno.EEXIST and os.path.isdir(path):
pass
else:
raise
import shutil
def copytree(src, dst, symlinks=False, ignore=None):
if not os.path.exists(dst):
os.makedirs(dst)
for item in os.listdir(src):
s = os.path.join(src, item)
d = os.path.join(dst, item)
if os.path.isdir(s):
copytree(s, d, symlinks, ignore)
else:
if not os.path.exists(d) or os.stat(s).st_mtime - os.stat(d).st_mtime > 1:
shutil.copy2(s, d)
pwd = os.getcwd()
i=0
k=0
HerwigSetup="source /unix/cedar/software/sl6/Herwig-7.1.0/setupEnv.sh"
ConturSetup="source /home/jmb/svn-managed/contur/setupContur.sh"
for i in range(10,100,10):
for j in range(1000,10000,500):
i = int(i)
j = int(j)
modelpath = 'mP_'+ str(i) + '_f0_' + str(j)
copytree(pwd + '/GridPack',modelpath)
#mkdir_p(str(modelpath))
HerwigString = ''
HC=open('GridPack/Neutral_scalar_CPodd_UFO/nscpo_8.in', 'r')
HerwigString += 'read FRModel.model \n'
HerwigString += 'set /Herwig/FRModel/Particles/phiNP:NominalMass ' +str(i) + '.*GeV \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0B '+str(j)+' \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0W '+str(j)+' \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0G 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0H 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0u 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0d 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0l 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0gam 1000000 \n'
HerwigString += 'set /Herwig/FRModel/FRModel:f0Z 1000000 \n'
HerwigString += str(HC.read())
HC.close()
RunCard = open(str(modelpath + '/LHC.in'), 'w')
RunCard.write(str(HerwigString))
RunCard.close()
subprocess.call([HerwigSetup], shell=True)
subprocess.call([ConturSetup], shell=True)
os.chdir(modelpath)
subprocess.call(['Herwig read LHC.in'], shell=True)
batch_command = ''
batch_command += HerwigSetup + '; '
batch_command += 'cd ' + pwd + '/' + modelpath +'; '
batch_command += ConturSetup + "; "
numEv=30000
batch_command += 'Herwig run --seed='+str(i)+str(j)+' --tag='+str(modelpath)+' --jobs=5 --numevents='+ str(numEv) +' LHC.run;'
batch_filename = str(modelpath)+'.sh'
batch_submit = open(batch_filename, 'w')
batch_submit.write(batch_command)
batch_submit.close()
subprocess.call([ "qsub -q medium " + batch_filename],shell=True )
os.chdir(pwd)

File Metadata

Mime Type
text/x-python
Expires
Sat, Dec 21, 3:50 PM (1 d, 20 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023349
Default Alt Text
BatchSubmit_NSO.py (2 KB)

Event Timeline