POWHEG-MC-Event-generation/create_scripts.py at production · HerrHorizontal/POWHEG-MC-Event-generation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87


import sys
import os
import stat
import string
from glob import glob

def create_scripts (nbatches, processes, decay = False):

    for process in processes:

        # check if given argument is a process directory in POWHEG-BOX-V2 or POWHEG-BOX-RES
        if not(os.path.isdir(process)) or  'POWHEG' not in os.path.dirname(process):
            print 'Error: Argument ' + str(process) + ' is not a POWHEG process directory' + '\nGeneration aborted!'
            continue

        # generate directory system
        work_dir = os.getcwd()
        process_name = os.path.basename(os.path.abspath(process))
        if not os.path.exists('./GenData'):
            os.mkdir('./GenData')
        if not os.path.exists('./GenData/' + process_name):
            os.mkdir('./GenData/' + process_name)

        # remove all already existing jobscripts in the GenData directory
        os.chdir(os.path.join(os.path.abspath("./GenData/"), process_name))
        for script in glob('*.sh'):
            os.remove(script)
        os.chdir(work_dir)

        # generate a jobscript for each batch
        for batch in range(nbatches):
            filename = os.path.abspath('./GenData/' + process_name) + '/jobscript_batch_' + str(batch) + '.sh'
            if os.path.isfile(filename):
                os.remove(filename)
            with open(filename, 'wb') as scriptfile:
                scriptfile.write('#!/bin/bash\n\n')
                scriptfile.writelines(['# set the CMSSW environment and load all needed modules\n',
                                       '#module use -a /afs/desy.de/group/cms/modulefiles/\n',
                                       'export VO_CMS_SW_DIR=/cvmfs/cms.cern.ch\n',
                                       'source $VO_CMS_SW_DIR/cmsset_default.sh\n',
                                       'export CMSSW_GIT_REFERENCE=/nfs/dust/cms/user/mhorzela/.cmsgit-cache\n',
                                       'alias cd=\'cd -P\'\n',
                                       'myvarcwd=$PWD\n',
                                       'cd /nfs/dust/cms/user/mhorzela/CMSSW/CMSSW_9_4_9/src\n',
                                       'eval `scramv1 runtime -sh`\n',
                                       'cd ~\n',
                                       'echo "setup CMSSW_949 and stuff"\n',
                                       'cd $myvarcwd\n\n',
                                       '#add the LHAPDF library path to PATH\n',
                                       'PATH=$PATH:/cvmfs/cms.cern.ch/slc6_amd64_gcc630/external/lhapdf/6.2.1-fmblme/bin/\n',
                                       'LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/cvmfs/cms.cern.ch/slc6_amd64_gcc630/external/lhapdf/6.2.1-fmblme/bin/\n',
                                       '#add the FASTJET library path to PATH\n',
                                       'PATH=$PATH:/cvmfs/cms.cern.ch/slc6_amd64_gcc630/external/fastjet/3.1.0/bin/\n',
                                       'echo "setup POWHEG"\n\n'])

                scriptfile.writelines(['# run POWHEG process ' + str(process_name) + ' batch number ' + str(batch) + '\n',
                                       'cd ' + str(os.path.abspath(process)) + '\n'])
                if decay == True:
                    scriptfile.writelines(['echo pwgevents-' + str(batch).zfill(4) + '.lhe | ./lhef_decay\n',
                                            'echo "</LesHouchesEvents>" | gzip - | cat - >> pwgevents-' + str(batch).zfill(4) + '-decayed.lhe \n'])

                else:
                    scriptfile.writelines(['echo ' + str(batch) + ' | ./pwhg_main' '\n'])

                scriptfile.write('cd ' + str(os.path.abspath(work_dir)))


        print 'created scripts for POWHEG process ' + process_name
        status = os.stat(filename)
        os.chmod(filename, status.st_mode | stat.S_IEXEC)

def main (args = sys.argv[1:]):
    if not(args[0].isdigit()):
        print 'Wrong usage! First argument has to be an integer, representing the number of batches, further arguments should be directors to a POWHEG process\nAbort!'
        exit(0)

    nbatches = int(args[0])
    processes = args[1:]

    create_scripts(nbatches = nbatches, processes = processes)


if __name__ == '__main__':
    main()