main.c

/* Copyright (c) Anglo-Australian Telescope Board, 1999.
   Not to be used for commercial purposes without AATB permission. */


#include "DitsFix.h"
#include "Sdp.h"
#include "Git.h"
#include "dul_err.h"
#include "drs_err.h" 

/* globals */
/* string containing drexec version */
char versionStr[128];

/* externals */
extern void combine_image_(DitsArgType*,StatusType*);
extern void combine_spectra_(DitsArgType*,StatusType*);
/* extern void flx_calib_(DitsArgType*,StatusType*); */
/* extern void gen_calib_(DitsArgType*,StatusType*); */
extern void cmfim_gfibpos_(DitsArgType*,StatusType*);
extern void plot_file_(DitsArgType*,StatusType*);
extern void reduce_arc_(DitsArgType*,StatusType*);
extern void reduce_bias_(DitsArgType*,StatusType*);
extern void reduce_dark_(DitsArgType*,StatusType*);
extern void reduce_fflat_(DitsArgType*,StatusType*);
extern void reduce_lflat_(DitsArgType*,StatusType*);
extern void reduce_fflx_(DitsArgType*,StatusType*);
extern void reduce_object_(DitsArgType*,StatusType*);
extern void reduce_sky_(DitsArgType*,StatusType*);
extern void splice_aaomega_(DitsArgType*,StatusType*);
extern void detect_cosmic_streaks_(DitsArgType*,StatusType*);

/* forward declarations */
static void combineImageAction(StatusType*);
static void combineSpectraAction(StatusType*);
/* static void flxCalibAction(StatusType*); */
/* static void genCalibAction(StatusType*); */
static void find2dfFibresAction(StatusType*);
static void plotAction(StatusType*);
static void plotKick(StatusType*);
static void reduceArcAction(StatusType*);
static void reduceBiasAction(StatusType*);
static void reduceDarkAction(StatusType*);
static void reduceFFlatAction(StatusType*);
static void reduceLFlatAction(StatusType*);
static void reduceFFluxAction(StatusType*);
static void reduceObjAction(StatusType*);
static void reduceSkyAction(StatusType*);
static void spliceAction(StatusType*);
static void detectCosmicStreaksAction(StatusType*);
static void exitAction(StatusType*);
static void testAction(StatusType*);
static void testKick(StatusType*);
void SetPgplotRequest(char key, float x, float y);
void GetPgplotRequest(char* key, float* x, float* y);


/*---------------------------------------------------------------------------*/

int main(int argc,char** argv)
/*
 * Execution of drexec server DRAMA task begins here.
 *
 * This task acts as an action server for 2dfdr. drcontrol acts as the client
 * DRAMA task sending this task requests for 2dfdr processing.
 *
 * N.B. This task BREAKS one of the cardinal rules of DRAMA -
 *      1. Tasks accepting and servicing actions shall reschedule immediately
 *         thereby being ready for any new actions (or kicks).
 *      Notice the DRAMA actions served by this task are quite CPU intensive,
 *      and the program does NOT reschedule until they complete.
 */
{
    StatusType status = DITS__OK;
    char taskName[40];

    /* this task executes following DRAMA actions */
    DitsActionDetailsType actions[] = {
	{ combineImageAction,  plotKick,0,0,0,0,"COMBINE_IMAGE"},
	{ combineSpectraAction,plotKick,0,0,0,0,"COMBINE_SPECTRA"},
/*	{ flxCalibAction,      0,0,0,0,0,"FLX_CALIB"}, */
/*	{ genCalibAction,      0,0,0,0,0,"GEN_CALIB"}, */
	{ find2dfFibresAction, plotKick,0,0,0,0,"GFIBPOS"},
	{ plotAction,          plotKick,0,0,0,0,"PLOT"},
	{ reduceArcAction,     plotKick,0,0,0,0,"REDUCE_ARC"},
	{ reduceBiasAction,    plotKick,0,0,0,0,"REDUCE_BIAS"},
	{ reduceDarkAction,    plotKick,0,0,0,0,"REDUCE_DARK"},
	{ reduceFFlatAction,   plotKick,0,0,0,0,"REDUCE_FFLAT"},
	{ reduceLFlatAction,   plotKick,0,0,0,0,"REDUCE_LFLAT"},
	{ reduceFFluxAction,   plotKick,0,0,0,0,"REDUCE_FFLX"},
	{ reduceObjAction,     plotKick,0,0,0,0,"REDUCE_OBJECT"},
	{ reduceSkyAction,     plotKick,0,0,0,0,"REDUCE_SKY"},
	{ spliceAction,        plotKick,0,0,0,0,"SPLICE"},
	{ detectCosmicStreaksAction,plotKick,0,0,0,0,"DETECT_CSTRKS"},
	{ exitAction,          0,0,0,0,0,"EXIT"},
        { testAction,          plotKick, 0, 0, 0, 0, "TEST"}
    };
/*	{ testAction,          0,0,0,0,0,"TEST"}*/

    /* The task name is normally given in argv[1]. Multiple copies of the
       task can  be started using different task names. A default name is
       used if one is not given. */
    strcpy(taskName,(argc > 1) ? argv[1] : "DREXEC");
 
    /* initialise DRAMA
       DITS_M_NOEXHAND => we will install exit handler */
    DitsAppInit(taskName,300000,DITS_M_NOEXHAND,0,&status);

    /* define supported actions */
    DitsPutActions(DitsNumber(actions),actions,&status);  

    /* enter DRAMA main loop
       This returns when action returns DITS_REQ_EXIT */
    DitsMainLoop(&status);

     /* here returned from DRAMA main loop, shutdown DRAMA and exit */
    return DitsStop(taskName,&status);
}

/*---------------------------------------------------------------------------*/

void combineImageAction(
    StatusType* pStatus)
/*
 * Handler for "COMBINE_IMAGE" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* call Fortran COMBINE_IMAGE() */
    combine_image_(&args,pStatus);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void combineSpectraAction(
    StatusType* pStatus)
/*
 * Handler for "COMBINE_SPECTRA" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout 
    SdsList(args,pStatus);*/

    /* call Fortran COMBINE_SPECTRA() */
    combine_spectra_(&args,pStatus);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void flxCalibAction(
    StatusType* pStatus)
/*
 * Handler for "FLX_CALIB" DRAMA action.
 *
 * This action is not currently working. See flux_cal.f.
 */
{
    DitsArgType args;
/*    SdsIdType localArgs; */

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
 /*    SdsCopy(args,&localArgs,pStatus); */

    /* call Fortran and ???() */
/*    flx_calib_(&args,pStatus); */

    /* delete local SDS copy */
/*    SdsDelete(localArgs,pStatus); */

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void genCalibAction(
    StatusType* pStatus)
/*
 * Handler for "GEN_CALIB" DRAMA action.
 *
 * This action is not currently working. See flux_cal.f.
 */
{
    DitsArgType args;
/*    SdsIdType localArgs; */

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout 
    SdsList(args,pStatus);*/

    /* make copy so it can be changed (external SDS -> internal) */
 /*    SdsCopy(args,&localArgs,pStatus); */

    /* call Fortran and ???() */
/*    gen_calib_(&args,pStatus); */

    /* delete local SDS copy */
/*    SdsDelete(localArgs,pStatus); */

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void find2dfFibresAction(
    StatusType* pStatus)
/*
 * Handler for "GFIBPOS" (find 2df fibre positions) DRAMA action.
 *
 * This action is not currently working. See find_fib_2df.f.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout 
    SdsList(args,pStatus);*/


    /* call Fortran CMFIM_GFIBPOS() */
/*    cmfim_gfibpos_(&args,pStatus); */

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void plotAction(
    StatusType* pStatus)
/*
 * Handler for "PLOT" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType inpArgs, outArgs;
    SdsIdType id2;


    if (*pStatus != DITS__OK) return;

    /* Get our input SDS argument structure */
    /* Note: This is an external SDS object and cannot be modified*/
    args = DitsGetArgument();

    SdsCopy(args,&inpArgs,pStatus);
    
    /* Create a return sds object that can be written to*/
    SdsNew(0,"retArgs",0,NULL,SDS_STRUCT,0,NULL,&outArgs,pStatus);

    /* call the Fortran subroutine PLOT2() */
    plot2_(&inpArgs,&outArgs,pStatus);

    /* Set the return SDS object to be outArgs*/
    DitsPutArgument (outArgs,DITS_ARG_DELETE, pStatus);
 
    DitsPutRequest(DITS_REQ_END,pStatus);

}

/*---------------------------------------------------------------------------*/

void plotKick(
    StatusType *pStatus)
/*
 * Handler for kicking the "Plot" DRAMA action.
 */
{
    DitsArgType args;
    char key;
    char Buffer[10];
    float x,y;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();
    
    /* Obtain the KEY character and the X Y values*/    
       /* Note: ArgGetc(args,"KEY",&key,pStatus) doesn't appear to work so
                we use ArgGetString instead                                */
    ArgGetString(args,"KEY",10,Buffer,pStatus);
    key=Buffer[0];
    ArgGetf(args,"X",&x,pStatus);
    ArgGetf(args,"Y",&y,pStatus);  

    /* Store these values with Pgplot Request Set Functions*/
    SetPgplotRequest(key,x,y);

    /* Now this is the bit I am not sure about but is supposed to unblock
       the DitsActionWait command */
    DitsPutRequest(DITS_REQ_STAGE,pStatus); 

}

/*---------------------------------------------------------------------------*/

void reduceArcAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_ARC" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran REDUCE_ARC() */
    reduce_arc_(&localArgs,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce arc error");
    }

    /* delete local SDS copy and free id */
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceBiasAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_BIAS" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* call Fortran REDUCE_BIAS() */
    reduce_bias_(&args,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce bias error");
    }

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceDarkAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_DARK" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* call Fortran REDUCE_DARK() */
    reduce_dark_(&args,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce dark error");
    }

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceFFlatAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_FFLAT" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran REDUCE_FFLAT() */
    reduce_fflat_(&localArgs,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce flat error");
    }

    /* delete local SDS copy and free id */
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceLFlatAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_LFLAT" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* call Fortran REDUCE_LFLAT() */
    reduce_lflat_(&args,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce long flat error");
    }

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceFFluxAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_FFLX" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran REDUCE_OBJECT() */
    reduce_fflx_(&localArgs,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce object error");
    }

    /* delete local SDS copy and free its id*/
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}


/*---------------------------------------------------------------------------*/

void reduceObjAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_OBJECT" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran REDUCE_OBJECT() */
    reduce_object_(&localArgs,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce object error");
    }

    /* delete local SDS copy and free its id*/
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void reduceSkyAction(
    StatusType* pStatus)
/*
 * Handler for "REDUCE_SKY" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran REDUCE_SKY() */
    reduce_sky_(&localArgs,pStatus);
    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","reduce sky error");
    }

    /* delete local SDS copy and free its id*/
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void spliceAction(
    StatusType* pStatus)
/*
 * Handler for "SPLICE" DRAMA action.
 */
{
    DitsArgType args;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* call Fortran SPLICE_AAOMEGA() */
    splice_aaomega_(&args,pStatus);
    if (*pStatus != DITS__OK) {
        ErsRep(0,pStatus,"%s","splice error");
    }

    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void detectCosmicStreaksAction(
    StatusType* pStatus)
/*
 * Handler for "DETECT_CSTRKS" DRAMA action.
 */
{
    DitsArgType args;
    SdsIdType localArgs = 0;
    StatusType ignore = STATUS__OK;

    if (*pStatus != DITS__OK) return;

    /* get our SDS argument structure */
    args = DitsGetArgument();

    /* for debug print our argument structure on stdout
    SdsList(args,pStatus); */

    /* make copy so it can be changed (external SDS -> internal) */
    SdsCopy(args,&localArgs,pStatus);

    /* call Fortran DETECT_COSMIC_STREAKS()  */
    detect_cosmic_streaks_(&args,pStatus);

    if (*pStatus!=DITS__OK) {
        ErsRep(0,pStatus,"%s","cosmic streak detection error error");
    }

    /* delete local SDS copy and free its id*/
    SdsDelete(localArgs,&ignore);
    SdsFreeId(localArgs,&ignore);

    DitsPutRequest(DITS_REQ_END,pStatus);
}


/*---------------------------------------------------------------------------*/

void exitAction(
    StatusType *pStatus)
/*
 * Handler for "EXIT" DRAMA action.
 */
{
    DitsPutRequest(DITS_REQ_EXIT,pStatus);
}

/*---------------------------------------------------------------------------*/

void testAction(
    StatusType *pStatus)
/*
 * Handler for "TEST" DRAMA action.
 */
{
    printf("Test Drama Action Called\n");
    DitsPutRequest(DITS_REQ_MESSAGE,pStatus);
    DitsPutRequest(DITS_REQ_END,pStatus);
}

/*---------------------------------------------------------------------------*/

void testKick(
    StatusType *pStatus)
/*
 * Handler for "TESTKICK" DRAMA action.
 */
{
    printf("Test Drama Kick Called\n");
    DitsPutRequest(DITS_REQ_END,pStatus);
}