#! /bin/tcsh 
# Modifcation history:
# 1.2 06/06/02
# 1.3 10/09/08
# 2.0 25/08/11 - use fparkey to explicitly set 
#                keywords RA_PNT, DEC_PNT and 
#                PA_PNT, JAC, University of Leicester
# 2.1 30/08/11 - Small correction to the format, carriage return 
#                one line. More output to screen.
#

set version=2.0

nice +19

cat <<EOF 

--------------- SKYCAST VERSION $version -----------------

EOF

if ( $#argv != 3 ) then 
  cat <<EOF 

    Script to cast an XMM EPIC background dataset (or indeed any EPIC
    event dataset) onto the sky, at the position given by an input
    template event dataset (e.g. the event file you are interested in
    producing a background for). It is required that the XMM-SAS has been
    initiated, and that SAS variables, e.g. SAS_CCF, SAS_ODF, are
    correctly set.
    
    ( If using this script with official XMM "refilled" blank sky files, a
    subsequent usage of the SAS task backscale may result in the BACKSCAL
    keyword being set to 0. If this occurs, please see the web help pages at :
    http://xmm2.esac.esa.int/external/xmm_sw_cal/background/blank_sky.shtml )
    
    Use: skycast background template outfile
    
    background - Blank-sky background (or in fact any) event list
    template   - Template file (e.g. event file/image) with attitude info  
    outfile    - Background file with sky co-ords relevant to template file
    
    e.g. skycast /BGfiles/mos1BG.ds /mydata/myEvfile.fits mos1BGcorrectsky.fits
    
EOF
   exit
endif
    
set backgr=$1
set template=$2
set outfile=$3

echo "Skycast: input arguments, backgr, template, outfile: " $backgr $template $outfile

# SAS environment variables
echo "SAS environment variables are..."
printenv | grep SAS
echo ""

# XMM sas started?
if (! -e $SAS_DIR) then 
    echo "XMM SAS has not been started... exiting..."
    exit
endif

# Check input
if ( ! -e $backgr ) then
    echo "Error: User input file $backgr does not exist... exiting..."
    exit
else 
    echo "Input file... $backgr"
    if ( -e $outfile) then
	echo "Overwriting output file..."
	rm -f $outfile
    endif
    cp $backgr $outfile
endif

# CCF OK?
if (! -e $SAS_CCF) then 
    echo "Environment variable SAS_CCF does not exist...needs to be set... " 
    exit
endif

# Get pnt values from template dataset
if (! -e $template ) then 
    echo "Error: Template file $template does not exist... exiting..."
    exit
else 
    if (-e $template ) echo "Getting attitude information from template file $template..."
    rm -f temp0001
    fdump $template+0 temp0001 1 1-1
    set ra_pnt=`grep RA_PNT < temp0001 | head -1 | awk '{print $3}'`
    set dec_pnt=`grep DEC_PNT < temp0001 | head -1 | awk '{print $3}'`
    set pa_pnt=`grep PA_PNT < temp0001 | head -1 | awk '{print $3}'`
    set ranom=`grep REFXCRVL < temp0001 | head -1 | awk '{print $2}'`
    set decnom=`grep REFYCRVL < temp0001 | head -1 | awk '{print $2}'`
    rm -f temp0001
endif

# Using PNT attitude - do attcalc
echo "Calculating sky coordinates..."
attcalc -w 0 -V 0 eventset=$outfile attitudelabel=fixed fixedra=$ra_pnt fixeddec=$dec_pnt fixedposangle=$pa_pnt withatthkset=N refpointlabel=user nominalra=$ranom nominaldec=$decnom

# Add keywords to primary and events extensions
echo "Add keywords to primary and event extensions .."
fparkey $ra_pnt $outfile+0 RA_PNT add=yes
fparkey $ra_pnt $outfile+1 RA_PNT add=yes
fparkey $dec_pnt $outfile+0 DEC_PNT add=yes
fparkey $dec_pnt $outfile+1 DEC_PNT add=yes
fparkey $pa_pnt $outfile+0 PA_PNT add=yes
fparkey $pa_pnt $outfile+1 PA_PNT add=yes

echo "Skycast Completed."

exit