Removed as obsolete.
authorpolicheh <policheh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 Feb 2011 10:57:24 +0000 (10:57 +0000)
committerpolicheh <policheh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 Feb 2011 10:57:24 +0000 (10:57 +0000)
PHOS/PHOSPHYSda.cxx [deleted file]

diff --git a/PHOS/PHOSPHYSda.cxx b/PHOS/PHOSPHYSda.cxx
deleted file mode 100644 (file)
index abe7e20..0000000
+++ /dev/null
@@ -1,193 +0,0 @@
-/*
-
-PHOSPHYSda.cxx
-
-This program reads the DAQ data files passed as argument using the monitoring library.
-It stores the required event information into an output file.
-Messages on stdout are exported to DAQ log system.
-
-contact: Hisayuki.Torii@cern.ch
-
-*/
-
-#include "event.h"
-#include "monitor.h"
-#include "daqDA.h"
-
-#include <stdio.h>
-#include <stdlib.h>
-
-#include "AliRawReader.h"
-#include "AliRawReaderDate.h"
-#include "AliPHOSRawDecoder.h"
-#include "AliCaloAltroMapping.h"
-#include "AliPHOSDApi0mip.h"
-
-/* Main routine
-      Arguments: list of DATE raw data files
-*/
-int main(int argc, char **argv) {
-
-  int status;
-
-  /* log start of process */
-  printf("PHOSPHYSda program started\n");  
-
-  /* Time out counters*/
-  int timeout = 180; // 3mins
-  time_t start_time = time(0);
-
-  /* check that we got some arguments = list of files */
-  if (argc<2) {
-    printf("Wrong number of arguments\n");
-    return -1;
-  }
-
-
-  /* Retrieve mapping files from DAQ DB */ 
-  const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
-  for(Int_t iFile=0; iFile<4; iFile++) {
-    int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
-    if(failed) { 
-      printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]);
-      return -1;
-    }
-  }
-
-  /* Open mapping files */
-  AliAltroMapping *mapping[4];
-  TString path = "./";
-  path += "RCU";
-  TString path2;
-  for(Int_t i = 0; i < 4; i++) {
-    path2 = path;
-    path2 += i;
-    path2 += ".data";
-    mapping[i] = new AliCaloAltroMapping(path2.Data());
-  }
-
-  /* report progress */
-  daqDA_progressReport(10);
-
-  /* init some counters */
-  int nevents_physics=0;
-  int nevents_total=0;
-
-  /* init PHOS reader*/
-  AliRawReader *rawReader = NULL;
-
-  /* Prepare DA algorithm */
-  int iMod = 0;
-  AliPHOSDApi0mip* dapi0mip = new AliPHOSDApi0mip(iMod);
-
-  /* read the data files */
-  for(int n=1;n<argc;n++) {
-    status=monitorSetDataSource( argv[n] );
-    if (status!=0) {
-      printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
-      return -1;
-    }
-
-    /* report progress */
-    /* in this example, indexed on the number of files */
-    daqDA_progressReport(10+80*(n-1)/argc);
-
-    /* read the file */
-    for(;;) {
-      struct eventHeaderStruct *event;
-      eventTypeType eventT;
-
-      /* get next event */
-      status=monitorGetEventDynamic((void **)&event);
-      if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */
-      if (status!=0) {
-        printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
-        return -1;
-      }
-
-      /* retry if got no event */
-      if (event==NULL) {
-        break;
-      }
-
-      /* Status report */
-      // This is just for debugging
-      //if( nevents_total%1000 == 0 ){
-      //daqDA_progressReport(10+80*(n-1)/argc+nevents_total/10000);
-      //printf(" DEBUGDEBUG:: next event %d\n",nevents_total);
-      //}
-      
-      /* get time stamp */
-      time_t t;
-      t = event->eventTimestamp;
-      dapi0mip->SetTime(t);
-
-      /* Applying into DA algorithm */
-      int ix, iz, igain;
-      eventT=event->eventType;
-      if (eventT==PHYSICS_EVENT || eventT==CALIBRATION_EVENT) {
-       // ---------------------------------------------------------------
-       // User Defined Function
-       // ---------------------------------------------------------------
-       dapi0mip->NewEvent();
-       rawReader = new AliRawReaderDate((void*)event);
-       //AliPHOSRawDecoderv1 dc(rawReader,mapping);
-       AliPHOSRawDecoder dc(rawReader,mapping);
-       dc.SubtractPedestals(false);
-       while(dc.NextDigit()) {
-         //printf(".");
-         ix = dc.GetRow() - 1;
-         iz = dc.GetColumn() - 1;
-         if(dc.IsLowGain()) igain = 0; else igain = 1;
-         if( igain == 1 ){
-           dapi0mip->FillDigit(dc.GetEnergy(),ix,iz);
-         }
-         //dc.GetTime();
-       }
-       delete rawReader;
-       //dapi0mip->FillHist(); // This is not necesarry in this DA algorithm
-       dapi0mip->FillTree();
-       //if( nevents_total%1000 == 0 ){ // Dump for debugging
-       //dapi0mip->Print();
-       //}
-       // ---------------------------------------------------------------
-        nevents_physics++;
-      }
-      nevents_total++;
-
-      /* Check the time out */
-      if( nevents_total%1000 == 0 ){
-       time_t current_time = time(0);
-       if( current_time - start_time > timeout ){
-         free(event);
-         printf(" Warning: Exit due to the processing time exceed the limitation of %d sec\n",timeout);
-         n = argc;
-         break;
-       }
-      }
-      
-      /* free resources */
-      free(event);
-    }
-    
-  }
-
-  /* report progress */
-  daqDA_progressReport(90);
-
-  /* Delete DA algorithm for saving output into a file */
-  delete dapi0mip;
-
-  /* Store output files to the File Exchange Server */
-  char localfile[1024];
-  char remotefile[1024];
-  sprintf(localfile,"AliPHOSDApi0mip_mod%d.root",iMod);
-  sprintf(remotefile,"PHOSDApi0mip",iMod);
-  daqDA_FES_storeFile(localfile,remotefile);
-
-  /* report progress */
-  daqDA_progressReport(100);
-  
-  return status;
-}
-