source: trunk/MagicSoft/Mars/datacenter/scripts/runreflector@ 9054

#!/bin/sh
#
# ========================================================================
#
# *
# * This file is part of MARS, the MAGIC Analysis and Reconstruction
# * Software. It is distributed to you in the hope that it can be a useful
# * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
# * It is distributed WITHOUT ANY WARRANTY.
# *
# * Permission to use, copy, modify and distribute this software and its
# * documentation for any purpose is hereby granted without fee,
# * provided that the above copyright notice appear in all copies and
# * that both that copyright notice and this permission notice appear
# * in supporting documentation. It is provided "as is" without express
# * or implied warranty.
# *
#
#
#   Author(s): Daniela Dorner  05/2007 <mailto:dorner@astro.uni-wuerzburg.de>
#
#   Copyright: MAGIC Software Development, 2000-2007
#
#
# ========================================================================
#
# This script runs the reflector.
#

source `dirname $0`/sourcefile
printprocesslog "INFO starting $0"
program=reflector
column=fReflectorFileAvail
detectordir=/home/msmeyer/MC/MagicSoft/Simulation/Detector

set -C

# printing input card to stdout
function printinputcard()
{
   echo $reflectorversion
   echo "verbose_level 4"
   echo "seeds "$seeds
   echo "fixed_target 0. 0. "
   echo "telescopes_layout 1 1"
   echo "telescope_position 0.0 0.0"
   if [ $wobble -eq 2 ]
   then 
      echo "wobble_mode 1"
   else
      echo "wobble_mode 0"
   fi
   echo "ct_file $detectordir/Data/magic.def"
   echo "reflectivity_file $detectordir/Data/reflectivity.dat"
   echo "axisdev_file $detectordir/Data/axisdev.dat"
   echo "output_file "$outfile
   echo "atm_model ATM_MagicWinter"
   echo "cer_files "
   echo $infile
}


scriptlog=$runlogpath/run$program-$datetime.log
date >> $scriptlog 2>&1

# get sequence # 
gettodo >> $scriptlog 2>&1
run=$process
# get reflector runnumber
query="SELECT fReflectorRunNumber FROM MCRunData where fMCRunNumber="$run
reflectorrunno=`sendquery`

# lock sequ
lockfile=$lockpath/lock-$table-$column-$reflectorrunno.txt
checklock >> $scriptlog 2>&1

echo "run $program for run $run..." >> $scriptlog 2>&1
printprocesslog "INFO starting $program for run $run"

primvar=$run
setstatus "start" >> $scriptlog 2>&1

echo "run $program for run $run "
echo "create input card... "

# get values for inputcard
query="SELECT fReflectorSeed, fReflectorSeed+1 FROM MCReflectorRunData where fReflectorRunNumber="$reflectorrunno
seeds=`sendquery`
query="SELECT fObservationModeKEY FROM MCReflectorRunData where fReflectorRunNumber="$reflectorrunno
wobble=`sendquery`
query="SELECT DATE_FORMAT(fRunStart, '%Y/%m/%d') FROM MCReflectorRunData where fReflectorRunNumber="$reflectorrunno
date=`sendquery`
query="SELECT fParticleTypeName FROM MCCorsikaRunData "
query=$query" LEFT JOIN ParticleType ON MCCorsikaRunData.fParticleTypeKEY=ParticleType.fParticleTypeKEY "
query=$query" LEFT JOIN MCReflectorRunData ON MCCorsikaRunData.fCorsikaRunNumber=MCReflectorRunData.fCorsikaRunNumber "
query=$query" WHERE fReflectorRunNumber="$reflectorrunno
particle=`sendquery`
outpath=$mcpath"/"$program"/"$date
date2=`echo $date | sed -e 's/\///g'`
outfile=$outpath"/"$date2"_"`printf %06d $reflectorrunno`"_"$particle".rfl"
query="SELECT fCorsikaRunNumber FROM MCReflectorRunData where fReflectorRunNumber="$reflectorrunno
infile=$mcpath"/corsika/"$date"/cer"`printf %06d \`sendquery\``

echo "" >> $scriptlog 2>&1
echo "INPUTCARD:" >> $scriptlog 2>&1
echo "----------" >> $scriptlog 2>&1
echo "" >> $scriptlog 2>&1
echo "" >> $scriptlog 2>&1
printinputcard >> $scriptlog 2>&1

makedir $outpath >> $scriptlog 2>&1
if ls $outfile >/dev/null 2>&1
then
   rm -v $outfile
fi
logfile=$outpath/$program"-"`printf %06d $reflectorrunno`".log"

cd $detectordir/ReflectorII

echo "INPUTCARD:" >> $logfile 2>&1
printinputcard >> $logfile 2>&1
echo "" >> $logfile 2>&1
echo "" >> $logfile 2>&1
printinputcard | ./reflector

check1=$?

echo "the reflector returned: "$check1

case $check1 in
   22)   echo " check1=$check1 -> everything ok " >> $scriptlog 2>&1
         printprocesslog "INFO $program finished successfully for sequence $sequence"
         ;;
    *)   echo " check1=$check1 -> ERROR -> step has to be repeated" >> $scriptlog 2>&1
         printprocesslog "ERROR $program failed for sequence $sequence"
         com=$Freflector
         check=$check1
         ;;
esac

setstatus "stop" >> $scriptlog 2>&1

finish >> $scriptlog 2>&1