[ViewVC] Contents 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
Error occurred while calculating annotation data.
Log Message:	Scripts for a SUSY parameter scan
#	Content
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