← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to Template/bin/run.sh

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/bin/run.sh
 
1
#!/bin/bash
 
2
 
 
3
#############################################################################
 
4
#                                                                          ##
 
5
#                    MadGraph/MadEvent                                     ##
 
6
#                                                                          ##
 
7
# FILE : run.sh                                                            ##
 
8
# VERSION : 1.0                                                            ##
 
9
# DATE : 29 January 2008                                                   ##
 
10
# AUTHOR : Michel Herquet (UCL-CP3)                                        ##
 
11
#                                                                          ##
 
12
# DESCRIPTION : script to save command line param in a grid card and       ##
 
13
#   call gridrun                                                           ##
 
14
# USAGE : run [num_events] [iseed]                                         ##
 
15
#############################################################################
 
16
 
 
17
if [[ ! -d ./madevent ]]; then
 
18
        echo "Error: no madevent directory found !"
 
19
        exit
 
20
fi
 
21
 
 
22
export LD_LIBRARY_PATH=${LD_LIBRARY_PATH}:${PWD}/madevent/lib
 
23
 
 
24
card=./madevent/Cards/grid_card.dat
 
25
 
 
26
 
 
27
if [[  ($1 != "") && ("$2" != "") && ("$3" == "") ]]; then
 
28
   num_events=$1
 
29
   seed=$2
 
30
   echo "Updating grid_card.dat..."
 
31
   sed -i.bak "s/\s*\d*.*gevents/  $num_events = gevents/g" $card
 
32
   sed -i.bak "s/\s*\d*.*gseed/  $seed = gseed/g" $card
 
33
   gran=`awk '/^[^#].*=.*ngran/{print $1}' $card`
 
34
elif [[  ($1 != "") && ("$2" != "") && ("$3" != "") ]]; then
 
35
   num_events=$1
 
36
   seed=$2
 
37
   gran=$3
 
38
   echo "Updating grid_card.dat..."
 
39
   sed -i.bak "s/\s*\d*.*gevents/  $num_events = gevents/g" $card
 
40
   sed -i.bak "s/\s*\d*.*gseed/  $seed = gseed/g" $card
 
41
   sed -i.bak "s/\s*\d*.*ngran/  $gran = ngran/g" $card
 
42
else
 
43
   echo "Warning: input is not correct, using values from the grid_card.dat."
 
44
   if [[ ! -e $card ]]; then
 
45
        echo "Error: $card not found !"
 
46
        exit
 
47
   fi
 
48
   num_events=`awk '/^[^#].*=.*gevents/{print $1}' $card`
 
49
   seed=`awk '/^[^#].*=.*gseed/{print $1}' $card`
 
50
   gran=`awk '/^[^#].*=.*ngran/{print $1}' $card`
 
51
fi
 
52
 
 
53
 
 
54
echo "Now generating $num_events events with random seed $seed and granularity $gran"
 
55
 
 
56
 
 
57
if [[ ! -x ./madevent/bin/gridrun ]]; then
 
58
    echo "Error: gridrun script not found !"
 
59
    exit
 
60
else
 
61
    cd ./madevent
 
62
    ./bin/gridrun
 
63
fi
 
64
 
 
65
if [[ ! -e ./Events/unweighted_events.lhe ]]; then
 
66
    echo "Error: event file not found !"
 
67
    exit
 
68
else
 
69
    echo "Moving events"
 
70
    mv ./Events/unweighted_events.lhe ../events.lhe
 
71
    cd ..
 
72
fi
 
73
 
 
74
if [[ -e ./DECAY/decay ]]; then
 
75
    cd DECAY
 
76
    echo -$seed > iseed.dat
 
77
    for ((i = 1 ;  i <= 20;  i++)) ; do
 
78
        if [[ -e decay_$i\.in ]]; then
 
79
            echo "Decaying events..."
 
80
            mv ../events.lhe ../events_in.lhe
 
81
            ./decay < decay_$i\.in
 
82
        fi
 
83
    done
 
84
    cd ..
 
85
fi
 
86
 
 
87
if [[ -e ./REPLACE/replace.pl ]]; then
 
88
    for ((i = 1 ;  i <= 20;  i++)) ; do
 
89
        if [[ -e ./REPLACE/replace_card$i\.dat ]];then
 
90
            echo "Adding flavors..."
 
91
            mv ./events.lhe ./events_in.lhe
 
92
            cd ./REPLACE
 
93
            ./replace.pl ../events_in.lhe ../events.lhe < replace_card$i\.dat
 
94
            cd ..
 
95
        fi
 
96
    done
 
97
fi
 
98
 
 
99
# part added by Stephen Mrenna to correct the kinematics of the replaced
 
100
#  particles
 
101
if [[ -e ./madevent/bin/addmasses.py ]]; then
 
102
  mv ./events.lhe ./events.lhe.0
 
103
  python ./madevent/bin/addmasses.py ./events.lhe.0 ./events.lhe
 
104
fi  
 
105
 
 
106
gzip -f events.lhe
 
107
exit