|
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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
|
#!/bin/bash
set -o verbose
echo " ______________________________________ "
echo " Running Madgraph5.... "
echo " ______________________________________ "
repo=${1}
echo "%MSG-MG5 repository = $repo"
name=${2}
echo "%MSG-MG5 gridpack = $name"
decay=${3}
echo "%MSG-MG5 run decay = $decay"
replace=${4}
echo "%MSG-MG5 replace = $replace"
process=${5}
echo "%MSG-MG5 process = $process"
maxjetflavor=${6}
echo "%MSG-MG5 maxjetflavor = $maxjetflavor"
qcut=${7}
echo "%MSG-MG5 qcut = $qcut"
minmax_jet=${8}
echo "%MSG-MG5 minmax_jet = $minmax_jet"
min_jets=${9}
max_jets=${10}
echo "%MSG-MG5 min/max jet multiplicity = $min_jets / $max_jets"
nevt=${11}
echo "%MSG-MG5 number of events requested = $nevt"
rnum=${12}
echo "%MSG-MG5 random seed used for the run = $rnum"
ncpu=${13}
echo "%MSG-MG5 thread count requested = $ncpu (ignored)"
# retrieve the wanted gridpack from the official repository
fn-fileget -c `cmsGetFnConnect frontier://smallfiles` ${repo}/${name}_gridpack.tar.gz
# force the f77 compiler to be the CMS defined one
ln -sf `which gfortran` f77
ln -sf `which gfortran` g77
PATH=`pwd`:${PATH}
tar xzf ${name}_gridpack.tar.gz ; rm -f ${name}_gridpack.tar.gz ; cd madevent
## rename addmasses.py to addmasses.py.no
find . -name addmasses.py -exec mv {} {}.no \;
########### BEGIN - REPLACE process ################
# REPLACE script is runned automatically by run.sh if REPLACE dir is found ###
# REPLACE will replace el with el/mu/taus by default, if you need something else you need to edit the replace_card1.dat
if [ ${replace} == true ] ; then
cd ..
mkdir REPLACE
cat > replace_card1.dat <<EOF
# Enter here any particles you want replaced in the event file after ME run
# In the syntax PID : PID1 PID2 PID3 ...
# End with "done" or <newline>
11:11 13 15
-12: -12 -14 -16
-11:-11 -13 -15
12: 12 14 16
done
EOF
cp ./madevent/bin/replace.pl ./replace_card1.dat ./REPLACE/
cd -
fi
########## END - REPLACE #########################
###### BEGIN - DECAY process #################
# Decay is runned automatically by run.sh if DECAY dir is found
# To avoid compilation problem DECAY dir has been created and compiled when gridpack is created
# if not needed DECAY dir is deleted
if [ "${decay}" != true ] ; then
rm -rf DECAY
fi
######END - DECAY #####
# run the production stage
./run.sh ${nevt} ${rnum}
ls -al
file="events"
if [ ! -f ${file}.lhe.gz ]; then
echo "%MSG-MG5 events.lhe.gz file is not in the same folder with run.sh script, abort !!! "
exit
fi
cp ${file}.lhe.gz ${file}_orig.lhe.gz
gzip -d ${file}.lhe.gz
#_______________________________________________________________________________________
# check the seed number in LHE file.
echo " ______________________________________ "
echo " post processing started "
echo " ______________________________________ "
echo
if [ -f ${file}.lhe ] ; then
seed=`awk 'BEGIN{FS=" = gseed "}/gseed/{print $1}' ${file}.lhe`
number_event=`grep -c "</event>" ${file}.lhe`
fi
if [ $seed -eq $rnum ] ;then
echo "GSEED :$seed"
if [ $number_event -eq $nevt ] ;then
echo "NEVENT : $nevt "
else
echo "%MSG-MG5 Error: The are less events ( $number_event ) Post Production is cancelled."
# TO-DO You might want to save the events in case of inspection the events.
exit 1
fi
else
echo "%MSG-MG5 Error: Seed numbers doesnt match ( $seed )"
exit 1
fi
#_______________________________________________________________________________________
# post-process the LHE file.
#__________________________________________
# wjets/zjets
if [[ ${process} == wjets || ${process} == zjets ]] ; then
echo "%MSG-MG5 process V+jets"
python madevent/bin/mgPostProcv2.py -o ${file}_qcut${qcut}_mgPostv2.lhe -m -w -j ${maxjetflavor} -q ${qcut} -e 5 -s ${file}.lhe
fi
# qcd
if [ ${process} == qcd ] ; then
echo "%MSG-MG5 process QCD"
python madevent/bin/mgPostProcv2.py -o ${file}_qcut${qcut}_mgPostv2.lhe -q ${qcut} -j ${maxjetflavor} -e 5 -s ${file}.lhe
fi
# ttbar
if [ ${process} == ttbar ] ; then
echo "%MSG-MG5 process ttbar"
python madevent/bin/mgPostProcv2.py -o ${file}_qcut${qcut}_mgPostv2.lhe -m -w -t -j ${maxjetflavor} -q ${qcut} -e 5 -s ${file}.lhe
sed -i -e '/Rnd seed/d' -e '/MC partial width/d' -e '/Number of Events/d' -e '/Max wgt/d' -e '/Average wgt/d' -e '/Integrated weight/d' ${file}_qcut${qcut}_mgPostv2.lhe
fi
#__________________________________________
# If you have HT binned samples min/max jets might be different from file to file.
# So you can override the min/max jets decision and put by hand these from the command line
if [ $minmax_jet == true ] ;then
sed -i "s/ [0-9]* = minjets ! Smallest number of additional light flavour jets/ $min_jets = minjets ! Smallest number of additional light flavour jets/g" \
${file}_qcut${qcut}_mgPostv2.lhe
sed -i "s/ [0-9]* = maxjets ! Largest number (inclusive ktMLM matching multipl.)/ $max_jets = maxjets ! Largest number (inclusive ktMLM matching multipl.)/g" \
${file}_qcut${qcut}_mgPostv2.lhe
fi
mv ${file}_qcut${qcut}_mgPostv2.lhe ${name}_final.lhe
ls -l
echo
exit 0
|