← Back to branch summary

~maddevelopers/mg5amcnlo/WWW_tg

« back to all changes in this revision

Viewing changes to users/mardelcourt/PROC_407857/PROC_407857/bin/internal/Gridpack/run.sh

  • Committer: Martin Delcourt
  • Date: 2013-03-25 20:56:51 UTC
  • Revision ID: mdel@perso.be-20130325205651-f3bisj4z2jlfyjh3
added test et test2

Show diffs side-by-side

added added

removed removed

Lines of Context:
users/mardelcourt/PROC_407857/PROC_407857/bin/internal/Gridpack/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
 
# For Linux
23
 
export LD_LIBRARY_PATH=${LD_LIBRARY_PATH}:${PWD}/madevent/lib:${PWD}/HELAS/lib
24
 
# For Mac OS X
25
 
export DYLD_LIBRARY_PATH=${DYLD_LIBRARY_PATH}:${PWD}/madevent/lib:${PWD}/HELAS/lib
26
 
 
27
 
card=./madevent/Cards/grid_card.dat
28
 
 
29
 
 
30
 
if [[  ($1 != "") && ("$2" != "") && ("$3" == "") ]]; then
31
 
   num_events=$1
32
 
   seed=$2
33
 
   echo "Updating grid_card.dat..."
34
 
   sed -i.bak "s/\s*\d*.*gevents/  $num_events = gevents/g" $card
35
 
   sed -i.bak "s/\s*\d*.*gseed/  $seed = gseed/g" $card
36
 
   gran=`awk '/^[^#].*=.*ngran/{print $1}' $card`
37
 
elif [[  ($1 != "") && ("$2" != "") && ("$3" != "") ]]; then
38
 
   num_events=$1
39
 
   seed=$2
40
 
   gran=$3
41
 
   echo "Updating grid_card.dat..."
42
 
   sed -i.bak "s/\s*\d*.*gevents/  $num_events = gevents/g" $card
43
 
   sed -i.bak "s/\s*\d*.*gseed/  $seed = gseed/g" $card
44
 
   sed -i.bak "s/\s*\d*.*ngran/  $gran = ngran/g" $card
45
 
else
46
 
   echo "Warning: input is not correct, using values from the grid_card.dat."
47
 
   if [[ ! -e $card ]]; then
48
 
        echo "Error: $card not found !"
49
 
        exit
50
 
   fi
51
 
   num_events=`awk '/^[^#].*=.*gevents/{print $1}' $card`
52
 
   seed=`awk '/^[^#].*=.*gseed/{print $1}' $card`
53
 
   gran=`awk '/^[^#].*=.*ngran/{print $1}' $card`
54
 
fi
55
 
 
56
 
 
57
 
echo "Now generating $num_events events with random seed $seed and granularity $gran"
58
 
 
59
 
 
60
 
if [[ ! -x ./madevent/bin/gridrun ]]; then
61
 
    echo "Error: gridrun script not found !"
62
 
    exit
63
 
else
64
 
    cd ./madevent
65
 
    ./bin/gridrun $num_events $seed
66
 
fi
67
 
 
68
 
if [[ -e ./Events/GridRun_${seed}/unweighted_events.lhe.gz ]]; then
69
 
        gunzip ./Events/GridRun_${seed}/unweighted_events.lhe.gz
70
 
fi
71
 
 
72
 
if [[ ! -e  ./Events/GridRun_${seed}/unweighted_events.lhe ]]; then
73
 
    echo "Error: event file not found !"
74
 
    exit
75
 
else
76
 
    echo "Moving events from  events.lhe"
77
 
    mv ./Events/GridRun_${seed}/unweighted_events.lhe ../events.lhe
78
 
    cd ..
79
 
fi
80
 
 
81
 
if [[ -e ./DECAY/decay ]]; then
82
 
    cd DECAY
83
 
    echo -$seed > iseed.dat
84
 
    for ((i = 1 ;  i <= 20;  i++)) ; do
85
 
        if [[ -e decay_$i\.in ]]; then
86
 
            echo "Decaying events..."
87
 
            mv ../events.lhe ../events_in.lhe
88
 
            ./decay < decay_$i\.in
89
 
        fi
90
 
    done
91
 
    cd ..
92
 
fi
93
 
 
94
 
if [[ -e ./REPLACE/replace.pl ]]; then
95
 
    for ((i = 1 ;  i <= 20;  i++)) ; do
96
 
        if [[ -e ./REPLACE/replace_card$i\.dat ]];then
97
 
            echo "Adding flavors..."
98
 
            mv ./events.lhe ./events_in.lhe
99
 
            cd ./REPLACE
100
 
            ./replace.pl ../events_in.lhe ../events.lhe < replace_card$i\.dat
101
 
            cd ..
102
 
        fi
103
 
    done
104
 
fi
105
 
 
106
 
# part added by Stephen Mrenna to correct the kinematics of the replaced
107
 
#  particles
108
 
if [[ -e ./madevent/bin/internal/addmasses.py ]]; then
109
 
  mv ./events.lhe ./events.lhe.0
110
 
  python ./madevent/bin/internal/addmasses.py ./events.lhe.0 ./events.lhe
111
 
fi  
112
 
 
113
 
gzip -f events.lhe
114
 
exit