source: trunk/DataCheck/QuickLook/Step2a.sh@ 17050

#!/bin/bash
#
source `dirname $0`/../Sourcefile.sh
printprocesslog "INFO starting $0"

if [ "$1" = "" ]
then 
   hour=`date +%k`
   # get date (before 18h there is no new data to be processed)
   if [ $hour -le 7 ] || [ $hour -ge 19 ]
   then
      datepaths=( `date --date="-19HOUR" +%Y/%m/%d` )
   else
      datepaths=( `date +%Y/%m/%d --date="-12hour"` `date +%Y/%m/%d --date="-36hour"` `date +%Y/%m/%d --date="-60hour"` \
              `date +%Y/%m/%d --date="-84hour"` `date +%Y/%m/%d --date="-108hour"` `date +%Y/%m/%d --date="-132hour"` \
#              `date +%Y/%m/%d --date="-156hour"` `date +%Y/%m/%d --date="-180hour"` `date +%Y/%m/%d --date="-204hour"` \
              )
   fi
   #datepaths=( `find -L $datapath/star -mindepth 3 -type d | sort -r | sed "s/\${datapath_for_sed}\/star//g" | sed -e 's/^\///'` ) #all available dates in star-directory
else
   datepaths=( $1 )
fi

for datepath in ${datepaths[@]}
do 
   printprocesslog "INFO processing "$datepath
   #echo "INFO processing "$datepath
   night=`echo $datepath | sed -e 's/\///g'`
   
   numganymeds=0
   
   # getting all image files for this night
   printprocesslog "INFO get file list for night "$datepath
   files=( `find $anapath/star/$datepath -type f -name '*_I.root' 2>/dev/null | sort` )
   if [ ${#files[@]} -eq 0 ]
   then
      printprocesslog "INFO no image files available for night "$datepath
      continue
   fi
   # write data set file
   #dspath=$anapath/$source/datasets`printf %03d $min`min/$datepath
   #dspath=$anapath/$source/datasets_run/$datepath
   dspath=$anapath/datasets_run/$datepath
   makedir $dspath

   for (( i=1 ; i<=${#files[@]} ; i++ ))
   do
      night=`basename ${files[$i-1]} | cut -d_ -f1`
      run=`basename ${files[$i-1]} | cut -d_ -f2`
      printprocesslog "DEBUG night: "$night" run: "$run
      dsfile=$dspath"/dataset"$night"_"$run".txt"
      #echo $dsfile
      printprocesslog "INFO writing dataset file "$dsfile
      # write dataset only if merpp-log is available to make sure that star-file is complete
      if echo $resulttable1 | grep ISDC >/dev/null 2>&1
      then 
         query="SELECT LPAD(fSequenceID, 3, 0) FROM RunInfo WHERE fNight="$night" AND fRunID="$run
         seq=`sendquery`
         merpplog=`dirname ${files[$i-1]}`"/"$night"_"$seq"-merpp.log"
      else 
         merpplog=`echo ${files[$i-1]} | sed -e 's/_I.root/-merpp.log/'`
      fi
      stillrunning=`find $merpplog -mmin -1 2>/dev/null`
      if [ -e $merpplog ] && [ "$stillrunning" == "" ]
      then
         echo `dirname ${files[$i-1]}`" "`basename ${files[$i-1]}` > $dsfile
      else
         printprocesslog "INFO merpp still running - wait with file "${files[$i-1]}
         continue
      fi
      # process only if ganymedlogfiles is not yet there
      logfile=`echo $dsfile | sed -e 's/datasets/ganymed/' -e 's/dataset//' -e 's/[.]txt/-ganymed[.]log/'`
      if ! ls $logfile >/dev/null 2>&1
      then
         # fixme: why is this double?
         night=`basename ${files[$i-1]} | cut -d_ -f1`
         runid=`basename ${files[$i-1]} | cut -d_ -f2 | cut -d. -f1`
         drivefile=$auxdata"/"$datepath"/"`basename ${files[$i-1]} | cut -d_ -f1`".DRIVE_CONTROL_SOURCE_POSITION.fits"
         rawfilepath=$rawdata"/"$datepath"/"
         rawfilename=`basename ${files[$i-1]} | sed -e 's/_I.root/.fits/'`
         rawfile=`find $rawfilepath -name $rawfilename* 2>/dev/null`
         if ! [ "$rawfile" == "" ]
         then 
            printprocesslog "INFO processing "$rawfile
            tstopi=`$factpath/fitsdump -h $rawfile  2>/dev/null | grep TSTOPI | grep -E -o '[0-9]+'`
            tstopf=`$factpath/fitsdump -h $rawfile  2>/dev/null | grep TSTOPF | grep -E -o '[.][0-9]+'`
            tstop=${tstopi}${tstopf}
         fi
         #echo "tstop: "$tstop
         if [ ls $drivefile >/dev/null 2>&1 ] && ! [ "$tstop" == "" ]
         then
            printprocesslog "INFO getting coordinates from "$drivefile
            coordinates=( `${factpath}/fitsdump ${drivefile} -c Ra_src Dec_src -r --filter='Time<'${tstop}  2>/dev/null | tail -1 2>&1` )
            if [ "${coordinates[0]}" == "" ] || [ "${coordinates[1]}" == "" ]
            then
               printprocesslog "WARN couldn't get coordinates ("${coordinates[@]}") from "$drivefile
               #echo "WARN couldn't get coordinates ("${coordinates[@]}") from "$drivefile
               continue
            fi
            if [ "${coordinates[0]}" == "0" ] || [ "${coordinates[1]}" == "0" ]
            then
               printprocesslog "WARN coordinates "${coordinates[@]}
               #echo "WARN coordinates "${coordinates[@]}
               continue
            fi
         else
            printprocesslog "INFO getting coordinates from DB."
            query="SELECT fSourceKEY FROM RunInfo WHERE fNight="$night" AND fRunID="$runid
            sourcekey=`sendquery`
            if [ "$sourcekey" == "" ]
            then
               printprocesslog "WARN sourcekey empty - coordinates"${coordinates[@]}
               #echo "WARN sourcekey empty - coordinates"${coordinates[@]}" for "$rawfile
               continue
            fi
            query="SELECT Round(fRightAscension,6), Round(fDeclination,6) from source WHERE fSourceKey="$sourcekey
            coordinates=( `sendquery` )
            if [ "${coordinates[0]}" == "" ] || [ "${coordinates[1]}" == "" ]
            then
               printprocesslog "WARN couldn't get coordinates ("${coordinates[@]}") from DB."
               #echo "WARN couldn't get coordinates ("${coordinates[@]}") from "$drivefile
               continue
            fi
            if [ "${coordinates[0]}" == "0" ] || [ "${coordinates[1]}" == "0" ]
            then
               printprocesslog "WARN coordinates "${coordinates[@]}
               #echo "WARN coordinates "${coordinates[@]}
               continue
            fi
         fi

         printprocesslog "DEBUG "$logfile" "${files[$i-1]}" -"${coordinates[0]}"-"${coordinates[1]}
         printprocesslog "DEBUG `dirname $0`/RunGanymed.sh $dsfile "${coordinates[0]}" "${coordinates[1]}
         if [ "$qlasge" = "yes" ]
         then
            printprocesslog "$sgepath/qsub -b y -q fact_short -v AUTOMATIONSETUP=$AUTOMATIONSETUP -v dsfile=$dsfile -v ra=${coordinates[0]} -v dec=${coordinates[1]} -e bla.txt -o bla.txt `dirname $0`/RunGanymed.sh"
            $sgepath/qsub -b y -q fact_short -v AUTOMATIONSETUP=$AUTOMATIONSETUP -v dsfile=$dsfile -v ra=${coordinates[0]} -v dec=${coordinates[1]} -e bla.txt -o bla.txt `dirname $0`/RunGanymed.sh
         else
            `dirname $0`/RunGanymed.sh $dsfile ${coordinates[0]} ${coordinates[1]}
            #echo "ganymed for "$dsfile" "${coordinates[0]}" "${coordinates[1]}
         fi
         numganymeds=`echo " $numganymeds + 1 " | bc -l`
      fi
   done
   printprocesslog "INFO found "$numganymeds" ganymeds."
done

finish