PWG2/FORWARD/analysis2/Run.sh

   1 #!/bin/bash
   2
   3 ana=$ALICE_ROOT/PWG2/FORWARD/analysis2
   4 nev=-1
   5 noanal=0
   6 nodraw=0
   7 rebin=1
   8 vzmin=-10
   9 vzmax=10
  10 batch=0
  11 gdb=0
  12 proof=0
  13 mc=0
  14 type=
  15 cms=900
  16 hhd=1
  17 comp=1
  18 tit=
  19 pass1=Pass1.C
  20
  21 usage()
  22 {
  23 cat<<EOF
  24 Usage: $0 [OPTIONS]
  25
  26 Do Pass1 and Pass2 on ESD files in current directory.
  27
  28 Options:
  29         -h,--help               This help
  30         -n,--events N           Number of events           ($nev)
  31         -1,--pass1              Run only pass 1, no draw   ($nodraw)
  32         -2,--pass2              Run only pass 2, just draw ($noanal)
  33         -r,--rebin N            Rebin by N                 ($rebin)
  34         -v,--vz-min CM          Minimum value of vz        ($vzmin)
  35         -V,--vz-max CM          Maximum value of vz        ($vzmax)
  36         -t,--trigger TYPE       Select trigger TYPE        ($type)
  37         -e,--energy CMS         Center of mass energy      ($cms)
  38         -b,--batch              Do batch processing        ($batch)
  39         -P,--proof NWORKERS     Run in PROOF(Lite) mode    ($proof)
  40         -M,--mc                 Run over MC data           ($mc)
  41         -S,--title STRING       Set the title string       ($tit)
  42         -g,--gdb                Run in GDB mode            ($gdb)
  43         -H,--hhd                Do comparison to HHD       ($hhd)
  44         -O,--other              Do comparison to other     ($comp)
  45         -E,--eloss              Run energy loss script
  46
  47 TYPE is a comma or space separated list of
  48
  49   INEL        Inelastic triggers (V0A|V0C|SPD)
  50   INEL>0      As above + N_ch > 0 in -0.5<eta<+0.5
  51   NSD         Non-single diffractive ((VOA&VOC)|N_ch > 5 -1.9<eta<+1.9)
  52
  53 If NWORKERS is 0, then the analysis will be run in local mode.
  54 EOF
  55 }
  56
  57 toggle()
  58 {
  59     echo $((($1+1)%2))
  60 }
  61
  62
  63 while test $# -gt 0 ; do
  64     case $1 in
  65         -h|--help)           usage            ; exit 0;;
  66         -n|--events)         nev=$2           ; shift ;;
  67         -2|--pass2)          noanal=`toggle $noanal`   ;;
  68         -1|--pass1)          nodraw=`toggle $nodraw`   ;;
  69         -b|--batch)          batch=`toggle $batch`   ;;
  70         -P|--proof)          proof=$2         ; shift ;;
  71         -M|--mc)             mc=`toggle $mc`   ;;
  72         -g|--gdb)            gdb=`toggle $gdb`   ;;
  73         -H|--hhd)            hhd=`toggle $hhd`   ;;
  74         -O|--other)          comp=`toggle $comp`   ;;
  75         -r|--rebin)          rebin=$2         ; shift ;;
  76         -v|--vz-min)         vzmin=$2         ; shift ;;
  77         -V|--vz-max)         vzmax=$2         ; shift ;;
  78         -e|--energy)         cms=$2           ; shift ;;
  79         -S|--title)          tit="$2"         ; shift ;;
  80         -E|--eloss)          pass1=MakeELossFits.C ; nodraw=1 ;;
  81         -t|--type)
  82             if test "x$type" = "x" ; then type=$2 ; else type="$type|$2"; fi
  83             shift ;;
  84         *) echo "$0: Unknown option '$1'" >> /dev/stderr ; exit 1 ;;
  85     esac
  86     shift
  87 done
  88
  89 if test $nev -lt 0 ; then
  90     base=dndeta_xxxxxxx
  91 else
  92     base=`printf dndeta_%07d $nev`
  93 fi
  94 opts="-l -x"
  95 opts1=""
  96 redir=
  97 if test $nodraw -lt 1 ; then
  98     opts1="-q"
  99 fi
 100 if test $batch -gt 0 ; then
 101     opts="-b -q $opts"
 102     redir="2>&1 | tee ${base}.log"
 103     echo "redir=$redir"
 104 fi
 105 if test $noanal -lt 1 ; then
 106     rm -f AnalysisResult.root AliAODs.root
 107     rm -f fmdana.png
 108
 109     if test $gdb -gt 0 ; then
 110         export PROOF_WRAPPERCMD="gdb -batch -x $ALICE_ROOT/PWG2/FORWARD/analysis2/gdb_cmds --args"
 111     fi
 112     echo "Running aliroot ${opts} ${opts1} ${ana}/${pass1}\(\".\",$nev,$proof,$mc\) $redir"
 113     if test $batch -gt 0 ; then
 114         aliroot $opts $opts1 ${ana}/${pass1}\(\".\",$nev,$proof,$mc\) 2>&1 | tee ${base}.log
 115     else
 116         aliroot $opts $opts1 ${ana}/${pass1}\(\".\",$nev,$proof,$mc\)
 117     fi
 118     fail=$?
 119     rm -f event_stat.root \
 120         EventStat_temp.root \
 121         outputs_valid \
 122         `printf %09d.stat $nev`
 123     if  test $fail -gt 0               || \
 124         test ! -f AnalysisResults.root || \
 125         test ! -f AliAODs.root ; then
 126         echo "Analysis failed"
 127         exit 1
 128     fi
 129     echo "Analysis done"
 130 fi
 131
 132 if test $nodraw -lt 1 ; then
 133     rm -f result.root
 134     tit=`echo $tit | tr ' ' '@'`
 135     echo "Running aliroot ${opts} ${opts1} ${ana}/Pass2.C\(\".\",\"$type\",$cms,$vzmin,$vzmax,$rebin,$hhd,$comp\)"
 136     aliroot ${opts} ${ana}/Pass2.C\(\".\",\"$type\",$cms,$vzmin,$vzmax,$rebin,\"$tit\",$hhd,$comp\)
 137 fi
 138
 139
 140 #
 141 # EOF
 142 #