runset.sh

#!/bin/bash

. $SRT_PRIVATE_CONTEXT/ligo/env.sh

if ! [ $1 ]; then
  echo "$(basename $0) [type] [time] [files]"
  echo "  type:  The analysis type, that is, the part of"
  echo "         the name of the fcl between ligojob_ and .fcl"
  echo "  time:  Unix time stamp of the GW event or whatever"
  echo "  files: Any number of art files"
  exit 1
fi

# Should only be needed temporarily
CET_PLUGIN_PATH+=:$SRT_PRIVATE_CONTEXT/lib/$SRT_SUBDIR
export CET_PLUGIN_PATH

type=$1
shift

if ! [ $1 ]; then
  echo specify the unix time of the event
  exit 1
fi
unixtime=$1
shift

export TZ=UTC
fracsec=$(cut -d. -f 2 -s <<< $unixtime)
rfctime=$(date "+%Y-%m-%dT%H:%M:%S" -d @$unixtime).${fracsec}Z

if [ $? -ne 0 ]; then
  exit 1
fi

if ! $SRT_PRIVATE_CONTEXT/ligo/makefcl.sh $type $rfctime "$realgweventtime" $skymap; then
#if ! $SRT_PRIVATE_CONTEXT/ligo/makefcl.sh $type $rfctime "2019-03-24T13:29:01Z" $skymap; then
#if ! $SRT_PRIVATE_CONTEXT/ligo/makefcl.sh $type $rfctime "2019-03-24T13:29:01Z" ""; then
  exit 1
fi

fcl=$SRT_PRIVATE_CONTEXT/job/$(printf ligojob_$type.$rfctime.$realgweventtime.$(basename $skymap).fcl | sed s/:/-/g)

if ! [ -e $fcl ]; then
  echo I thought the fcl was $fcl
  echo but that does not exist
  exit 1
fi

if ! [ $1 ]; then
  echo You did not specify any files to processes.  Technical success\!
  exit 0
fi

out=$type.ligohists.root

for f in $@; do
  if ! [ -e "$f" ]; then
    echo $f does not exist > /dev/stderr
    exit 1
  fi
done

for f in $@; do
  base=$(basename $(basename $(basename $f .root) .raw) .artdaq)
  reco=$base.ligo${type}.root
  if echo $type | grep -q noreco; then
    reco=/dev/null
  fi
  log=$base.ligo.$type.log
  hist=$base.hists.root
  hists="$hists $hist"

  if [ -e $hist ]; then
    echo $hist already exists. Will overwrite.
  fi

  if [ -e $log ]; then
    echo $log already exists. Will overwrite.
  fi

  echo $SRT_PRIVATE_CONTEXT/ligo/runone.sh "$f" "$fcl" "$reco" "$hist" "$log"
done | ~mstrait/bin/parallel -j 3 --line-buffer

if [ $? -eq 0 ]; then
  hadd -f $out $hists
  echo Full histograms in $out
fi