ROOTPWA
doInt.sge
Go to the documentation of this file.
1 #!/bin/bash
2 # Iospin symmetrizing
3 # Integration
4 # Select queue
5 #$ -l short=TRUE,h_vmem=100M
6 # Divert output
7 #$ -o /nfs/nas/user/sneubert/LOGS
8 # Write error in same file
9 #$ -j y
10 
11 cd /nfs/nas/user/sneubert/scratch/COMPASS/rootpwa/trunk
12 #cd $ROOTPWA
13 source setup.sh
14 cd scripts/
15 
16 DATADIR=/nfs/hicran/project/compass/analysis/sneubert/5PiLTData3/
17 #DATADIR=/nfs/nas/user/sneubert/2004/5PiLTData3/
18 #DATADIR=/afs/e18/compass/analysis/sneubert/5PiModels/ImpTest/pieces2/
19 #DATADIR=$HOME/disk2/analysis/Kp_Kppipi/pwa/WA03_like/
20 
21 WAVESETDIR=$ROOTPWA/keyfiles/key5pi/
22 #WAVESETDIR=$ROOTPWA/keyfiles/keyKpipi/
23 
24 # choose which amplitudes to do
25 DOAMPS=1
26 DOPSP=1
27 DOACC=1
28 
29 
30 
31 cd $DATADIR
32 
33 # get bin for this TaskID
34 BINS=*
35 echo $BINS
36 # get the correct entry for this task
37 echo SGE_TASK_ID=$SGE_TASK_ID
38 if [ "$SGE_TASK_ID" == "undefined" ]; then
40 elif [ ! -n "$SGE_TASK_ID" ]; then
41  SGE_TASK_ID=1;
42 fi
43 echo SGE_TASK_ID=$SGE_TASK_ID
44 # start counting at 0 !!!
45 let "INDEX = SGE_TASK_ID"
46 echo INDEX=$INDEX
47 
48 MYBIN=`echo $BINS | gawk -F " " '{ print $'${INDEX}"}";`
49 
50 echo MYBIN=$MYBIN
51 
52 if [ ! -d $DATADIR/$MYBIN ]; then
53  echo "$BINDIR not found";
54  exit;
55 fi
56 
57 #WEIGHTS="-w ../$MYBIN.genbod.wht"
58 
59 if [ -n "$DOPSP" ]; then
60  echo "--- Doing PSPAMPS ----------------";
61  cd $DATADIR/$MYBIN/PSPAMPS;
62 
63  for i in $WAVESETDIR/*list.dat; do
64  echo "Processing $i";
65  ls -l $i;
66  addamp $i ./; #OLD/;
67  done;
68 
69  int $WEIGHTS *.amp > norm.int;
70  cp norm.int $DATADIR/$MYBIN/AMPS;
71 fi;
72 
73 if [ -n "$DOACC" ]; then
74  echo "--- Doing ACCAMPS ----------------";
75  cd $DATADIR/$MYBIN/ACCAMPS;
76 
77  for i in $WAVESETDIR/*list.dat; do
78  echo "Processing $i";
79  ls -l $i;
80  addamp $i ./; #OLD/;
81  done;
82 
83  int *.amp > accnorm.int;
84  cp accnorm.int $DATADIR/$MYBIN/AMPS;
85 fi;
86 
87 
88 if [ -n "$DOAMPS" ]; then
89  echo "--- Doing AMPS ----------------";
90  cd $DATADIR/$MYBIN/AMPS;
91 
92  for i in $WAVESETDIR/*list.dat; do
93  echo "Processing $i";
94  ls -l $i;
95  addamp $i ./; #OLD/;
96  done;
97 
98 fi;
99 
100 echo "------------------"
101 
102