Forward declaration of TVector3 is replaced by #include

[u/mrichter/AliRoot.git] / PHOS / PHOSPHYSda.cxx

diff --git a/PHOS/PHOSPHYSda.cxx b/PHOS/PHOSPHYSda.cxx
deleted file mode 120000 (symlink)
index 97c59029520796a8c9e1f94ab1778a539c1e9b2d..0000000000000000000000000000000000000000
--- a/PHOS/PHOSPHYSda.cxx
+++ /dev/null
@@ -1 +0,0 @@
-/home/polishch/Shuttle/HisaDA/DAsrc/PHOSPHYSda.cxx
\ No newline at end of file

diff --git a/PHOS/PHOSPHYSda.cxx b/PHOS/PHOSPHYSda.cxx
new file mode 100644 (file)
index 0000000000000000000000000000000000000000..abe7e2055bdb883b73d24d4f101344aa94c7003e
--- /dev/null
+++ b/PHOS/PHOSPHYSda.cxx
@@ -0,0 +1,193 @@
+/*
+
+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;
+}
+