[ViewVC] Annotation of: cvsroot/UserCode/auterman/SusyScan/GeneratorScan/fill

#!/bin/bash

function get () {
        RESULT=`grep -m $5 $1 $4 | cut -b $2-$3`
        if [ "$RESULT" = '' ]
        then
           RESULT=0.0
        fi
        return
}  

OPTIONS=1
#OUTPUT="scan.dat"
#rm -f $OUTPUT
#FILE=isajet_200_200_50_0_1_170.9.dat
FILE=$1
OUTPUT=$2
touch $OUTPUT

#MSuGra parameters
MZERO=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b -10`
MHALF=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b 11-20`
AZERO=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b 21-30`
TANB=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b 31-40`
SGNMU=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b 41-48`
MTOP=`grep -m 1 -A 1 "M_0,  M_(1/2)" $FILE | grep -v "M_0,  M_(1/2)" | cut -b 49-58`

#Model parameters
get "1/alpha_em" 14 26 $FILE $OPTIONS
iALPHAEM=$RESULT

get "M_1" 14 26 $FILE $OPTIONS
M1=$RESULT

get "mu(Q)" 14 26 $FILE $OPTIONS
muQ=$RESULT

get "M_H1" 14 26 $FILE $OPTIONS
MH1sq=$RESULT

get "2(thetaw)" 43 53 $FILE $OPTIONS
THETA_W=$RESULT

get "M_2" 43 53 $FILE $OPTIONS
M2=$RESULT

get "B(Q)" 43 53 $FILE $OPTIONS
BQ=$RESULT

get "M_H2" 43 53 $FILE $OPTIONS
MH2sq=$RESULT

get "a_s" 65 90 $FILE $OPTIONS
ALPHAS=$RESULT

get "M_3" 65 90 $FILE $OPTIONS
M3=$RESULT

get "mu(Q)" 65 90 $FILE $OPTIONS
Q=$RESULT

get "TANBQ" 65 90 $FILE $OPTIONS
TANBQ=$RESULT


#SUSY masses
get "M(GL)" 10 18 $FILE $OPTIONS
MGL=$RESULT

get "M(UL)" 10 18 $FILE $OPTIONS
MUL=$RESULT

get "M(B1)" 10 18 $FILE $OPTIONS
MB1=$RESULT

get "M(SN)" 10 18 $FILE $OPTIONS
MSN=$RESULT

get "M(NTAU)" 10 18 $FILE $OPTIONS
MNTAU=$RESULT

get "M(Z1)" 10 18 $FILE $OPTIONS
MZ1=$RESULT

get "M(W1)" 10 18 $FILE $OPTIONS
MW1=$RESULT

get "M(HL)" 10 18 $FILE $OPTIONS
MHL=$RESULT

get "M(UR)" 30 38 $FILE $OPTIONS
MUR=$RESULT

#get "M(UR)" 30 38 $FILE $OPTIONS
#MUR=$RESULT

get "M(B2)" 30 38 $FILE $OPTIONS
MB2=$RESULT

get "M(EL)" 30 38 $FILE $OPTIONS
MEL=$RESULT

get "M(TAU1)" 30 38 $FILE $OPTIONS
MTAU1=$RESULT

get "M(Z2)" 30 38 $FILE $OPTIONS
MZ2=$RESULT

get "M(W2)" 30 38 $FILE $OPTIONS
MW2=$RESULT

get "M(HH)" 30 38 $FILE $OPTIONS
MHH=$RESULT

get "M(DL)" 50 58 $FILE $OPTIONS
MDL=$RESULT

get "M(T1)" 50 58 $FILE $OPTIONS
MT1=$RESULT

get "M(ER)" 50 58 $FILE $OPTIONS
MER=$RESULT

get "M(TAU2)" 50 58 $FILE $OPTIONS
MTAU2=$RESULT

get "M(Z3)" 50 58 $FILE $OPTIONS
MZ3=$RESULT

get "M(HA)" 50 58 $FILE $OPTIONS
MHA=$RESULT

get "M(DR)" 70 78 $FILE $OPTIONS
MDR=$RESULT

get "M(T2)" 70 78 $FILE $OPTIONS
MT2=$RESULT

get "M(Z4)" 70 78 $FILE $OPTIONS
MZ4=$RESULT

get "M(H+)" 70 78 $FILE $OPTIONS
MHp=$RESULT

if [ "$SGNMU" != '' ]
then
        echo "$MZERO $MHALF $TANB $SGNMU $AZERO $MTOP $muQ $Q $M1 $M2 $M3 $MGL $MUL $MB1 $MSN $MNTAU $MZ1 $MW1 $MHL $MUR $MB2 $MEL $MTAU1 $MZ2 $MW2 $MHH $MDL $MT1 $MER $MTAU2 $MZ3 $MHA $MDR $MT2 $MZ4 $MHp" >> $OUTPUT
else
        echo "no solution for m0=$MZERO and m1/2=$MHALF"
fi
Revision:	1.1.1.1 (vendor branch)
Committed:	Fri Nov 30 13:10:50 2007 UTC (17 years, 5 months ago) by auterman
Branch:	tex, Demo, SusyScan, MAIN
CVS Tags:	start, HEAD
Changes since 1.1:	+0 -0 lines
Log Message:	Scripts for a SUSY parameter scan
#	User	Rev	Content
1	auterman	1.1	#!/bin/bash
2
3			function get () {
4			RESULT=`grep -m $5 $1 $4 \| cut -b $2-$3`
5			if [ "$RESULT" = '' ]
6			then
7			RESULT=0.0
8			fi
9			return
10			}
11
12			OPTIONS=1
13			#OUTPUT="scan.dat"
14			#rm -f $OUTPUT
15			#FILE=isajet_200_200_50_0_1_170.9.dat
16			FILE=$1
17			OUTPUT=$2
18			touch $OUTPUT
19
20			#MSuGra parameters
21			MZERO=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b -10`
22			MHALF=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b 11-20`
23			AZERO=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b 21-30`
24			TANB=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b 31-40`
25			SGNMU=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b 41-48`
26			MTOP=`grep -m 1 -A 1 "M_0, M_(1/2)" $FILE \| grep -v "M_0, M_(1/2)" \| cut -b 49-58`
27
28			#Model parameters
29			get "1/alpha_em" 14 26 $FILE $OPTIONS
30			iALPHAEM=$RESULT
31
32			get "M_1" 14 26 $FILE $OPTIONS
33			M1=$RESULT
34
35			get "mu(Q)" 14 26 $FILE $OPTIONS
36			muQ=$RESULT
37
38			get "M_H1" 14 26 $FILE $OPTIONS
39			MH1sq=$RESULT
40
41			get "2(thetaw)" 43 53 $FILE $OPTIONS
42			THETA_W=$RESULT
43
44			get "M_2" 43 53 $FILE $OPTIONS
45			M2=$RESULT
46
47			get "B(Q)" 43 53 $FILE $OPTIONS
48			BQ=$RESULT
49
50			get "M_H2" 43 53 $FILE $OPTIONS
51			MH2sq=$RESULT
52
53			get "a_s" 65 90 $FILE $OPTIONS
54			ALPHAS=$RESULT
55
56			get "M_3" 65 90 $FILE $OPTIONS
57			M3=$RESULT
58
59			get "mu(Q)" 65 90 $FILE $OPTIONS
60			Q=$RESULT
61
62			get "TANBQ" 65 90 $FILE $OPTIONS
63			TANBQ=$RESULT
64
65
66			#SUSY masses
67			get "M(GL)" 10 18 $FILE $OPTIONS
68			MGL=$RESULT
69
70			get "M(UL)" 10 18 $FILE $OPTIONS
71			MUL=$RESULT
72
73			get "M(B1)" 10 18 $FILE $OPTIONS
74			MB1=$RESULT
75
76			get "M(SN)" 10 18 $FILE $OPTIONS
77			MSN=$RESULT
78
79			get "M(NTAU)" 10 18 $FILE $OPTIONS
80			MNTAU=$RESULT
81
82			get "M(Z1)" 10 18 $FILE $OPTIONS
83			MZ1=$RESULT
84
85			get "M(W1)" 10 18 $FILE $OPTIONS
86			MW1=$RESULT
87
88			get "M(HL)" 10 18 $FILE $OPTIONS
89			MHL=$RESULT
90
91			get "M(UR)" 30 38 $FILE $OPTIONS
92			MUR=$RESULT
93
94			#get "M(UR)" 30 38 $FILE $OPTIONS
95			#MUR=$RESULT
96
97			get "M(B2)" 30 38 $FILE $OPTIONS
98			MB2=$RESULT
99
100			get "M(EL)" 30 38 $FILE $OPTIONS
101			MEL=$RESULT
102
103			get "M(TAU1)" 30 38 $FILE $OPTIONS
104			MTAU1=$RESULT
105
106			get "M(Z2)" 30 38 $FILE $OPTIONS
107			MZ2=$RESULT
108
109			get "M(W2)" 30 38 $FILE $OPTIONS
110			MW2=$RESULT
111
112			get "M(HH)" 30 38 $FILE $OPTIONS
113			MHH=$RESULT
114
115			get "M(DL)" 50 58 $FILE $OPTIONS
116			MDL=$RESULT
117
118			get "M(T1)" 50 58 $FILE $OPTIONS
119			MT1=$RESULT
120
121			get "M(ER)" 50 58 $FILE $OPTIONS
122			MER=$RESULT
123
124			get "M(TAU2)" 50 58 $FILE $OPTIONS
125			MTAU2=$RESULT
126
127			get "M(Z3)" 50 58 $FILE $OPTIONS
128			MZ3=$RESULT
129
130			get "M(HA)" 50 58 $FILE $OPTIONS
131			MHA=$RESULT
132
133			get "M(DR)" 70 78 $FILE $OPTIONS
134			MDR=$RESULT
135
136			get "M(T2)" 70 78 $FILE $OPTIONS
137			MT2=$RESULT
138
139			get "M(Z4)" 70 78 $FILE $OPTIONS
140			MZ4=$RESULT
141
142			get "M(H+)" 70 78 $FILE $OPTIONS
143			MHp=$RESULT
144
145			if [ "$SGNMU" != '' ]
146			then
147			echo "$MZERO $MHALF $TANB $SGNMU $AZERO $MTOP $muQ $Q $M1 $M2 $M3 $MGL $MUL $MB1 $MSN $MNTAU $MZ1 $MW1 $MHL $MUR $MB2 $MEL $MTAU1 $MZ2 $MW2 $MHH $MDL $MT1 $MER $MTAU2 $MZ3 $MHA $MDR $MT2 $MZ4 $MHp" >> $OUTPUT
148			else
149			echo "no solution for m0=$MZERO and m1/2=$MHALF"
150			fi