Added some ignores and moved script

[u/mrichter/AliRoot.git] / PWG2 / FORWARD / analysis2 / Run.sh

#!/bin/bash 

nev=10000
noanal=0
nodraw=0
rebin=1
vzmin=-10
vzmax=10
batch=0
gdb=0
proof=0
type=
cms=900
hhd=1
comp=1
anal=0
term=0
tit=

usage()
{
cat<<EOF
Usage: $0 [OPTIONS]

Do Pass1 and Pass2 on ESD files in current directory.  

Options:
	-h,--help		This help                  
	-n,--events N		Number of events           ($nev)
	-1,--pass1 		Run only pass 1, no draw   ($nodraw)
	-2,--pass2		Run only pass 2, just draw ($noanal)
	-r,--rebin N		Rebin by N                 ($rebin)
	-v,--vz-min CM          Minimum value of vz        ($vzmin)
	-V,--vz-max CM          Maximum value of vz        ($vzmax)
	-t,--trigger TYPE       Select trigger TYPE        ($type)
	-e,--energy CMS         Center of mass energy      ($cms)
	-b,--batch              Do batch processing        ($batch)
	-P,--proof		Run in PROOF(Lite) mode    ($proof)
	-A,--analyse-only       Run only analysis          ($anal)
	-T,--terminate-only     Run only terminate         ($term)
	-S,--title STRING       Set the title string       ($tit)
	-g,--gdb		Run in GDB mode    	   ($gdb)
	-H,--hhd		Do comparison to HHD	   ($hhd)
	-O,--other		Do comparison to other	   ($comp)

TYPE is a comma or space separated list of 
 
  INEL	      Inelastic triggers (V0A|V0C|SPD)
  INEL>0      As above + N_ch > 0 in -0.5<eta<+0.5
  NSD         Non-single diffractive ((VOA&VOC)|N_ch > 5 -1.9<eta<+1.9)

EOF
}

toggle()
{
    echo $((($1+1)%2))
}


while test $# -gt 0 ; do
    case $1 in 
	-h|--help)           usage            ; exit 0;; 
	-n|--events)         nev=$2           ; shift ;; 
	-2|--pass2)          noanal=`toggle $noanal`   ;; 
	-1|--pass1)          nodraw=`toggle $nodraw`   ;; 
	-b|--batch)          batch=`toggle $batch`   ;; 
	-P|--proof)          proof=`toggle $proof`   ;; 
	-A|--analyse-only)   anal=`toggle $anal`   ;; 
	-T|--terminate-only) term=`toggle $term`   ;; 
	-g|--gdb)            gdb=`toggle $gdb`   ;; 
	-H|--hhd)            hhd=`toggle $hhd`   ;; 
	-O|--other)          other=`toggle $other`   ;; 
	-r|--rebin)          rebin=$2         ; shift ;; 
	-v|--vz-min)         vzmin=$2         ; shift ;; 
	-V|--vz-max)         vzmax=$2         ; shift ;; 
	-e|--energy)         cms=$2           ; shift ;;
	-S|--title)          tit="$2"         ; shift ;;
	-t|--type)           
	    if test "x$type" = "x" ; then type=$2 ; else type="$type|$2"; fi
	    shift ;;
	*) echo "$0: Unknown option '$1'" >> /dev/stderr ; exit 1 ;;
    esac
    shift
done 

base=`printf dndeta_%07d $nev`
opts="-l -x"
redir=
if test $batch -gt 0 ; then 
    opts="-l -b -q -x" 
    redir="2>&1 | tee ${base}.log"
fi 
if test $noanal -lt 1 ; then 
    rm -f AnalysisResult.root AliAODs.root
    rm -f fmdana.png
    
    # Setup analysis flags
    af=0
    if test $proof -gt 0 ; then 
	af=2 
    else 
	if test $anal -gt 0 ; then 
	    af=8
	elif test $term -gt 0 ; then 
	    af=16
	else 
	    af=4
	fi
    fi

    if test $gdb -gt 0 ; then 
	export PROOF_WRAPPERCMD="gdb -batch -x $ALICE_ROOT/PWG2/FORWARD/analysis2/gdb_cmds --args"
    fi
    aliroot $opts $ALICE_ROOT/PWG2/FORWARD/analysis2/Pass1.C\(\".\",$nev,$af\) \
	$redir 
    rm -f event_stat.root EventStat_temp.root outputs_valid
    if test ! -f AnalysisResults.root || test ! -f AliAODs.root ; then 
	echo "Analysis failed" 
	exit 1
    fi
    echo "Analysis done"
fi

if test $nodraw -lt 1 ; then
    rm -f result.root 
    if test "x$tit" = "x" ; then 
	tit="$nev\ events,\ v_{z}#in[$vzmin,$vzmax],\ $type"
    else 
	tit=`echo $tit | tr ' ' '\ '` 
    fi
    aliroot ${opts} $ALICE_ROOT/PWG2/FORWARD/analysis2/Pass2.C\(\"AliAODs.root\",\"$type\",$cms,$vzmin,$vzmax,$rebin,\"$tit\",$hhd,$comp\)
fi


#
# EOF
#