Reverted to old name to comply the offline agreements (requirement of C.Cheshkov).
authorpolicheh <policheh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 3 Jun 2008 18:39:40 +0000 (18:39 +0000)
committerpolicheh <policheh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 3 Jun 2008 18:39:40 +0000 (18:39 +0000)
PHOS/PHOSGAINda.cxx [new file with mode: 0644]

diff --git a/PHOS/PHOSGAINda.cxx b/PHOS/PHOSGAINda.cxx
new file mode 100644 (file)
index 0000000..bc0c994
--- /dev/null
@@ -0,0 +1,174 @@
+/*
+contact: Boris.Polishchuk@cern.ch
+link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx
+reference run: /castor/cern.ch/alice/phos/2007/10/04/18/07000008249001.1000.root
+run type: PHYSICS, STANDALONE
+DA type: MON 
+number of events needed: 1000
+input files: RCU0.data  RCU1.data  RCU2.data  RCU3.data
+Output files: PHOS_Module2_Calib.root
+trigger types used: L0
+*/
+
+
+#include "event.h"
+#include "monitor.h"
+extern "C" {
+#include "daqDA.h"
+}
+
+#include <stdio.h>
+#include <stdlib.h>
+
+#include <TSystem.h>
+
+#include "AliRawReader.h"
+#include "AliRawReaderDate.h"
+#include "AliPHOSRcuDA1.h"
+#include "AliPHOSRawDecoder.h"
+#include "AliCaloAltroMapping.h"
+
+
+/* Main routine
+      Arguments: 
+      1- monitoring data source
+*/
+int main(int argc, char **argv) {
+
+  int status;
+  
+  if (argc!=2) {
+    printf("Wrong number of arguments\n");
+    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());
+  }
+  
+
+  /* define data source : this is argument 1 */  
+  status=monitorSetDataSource( argv[1] );
+  if (status!=0) {
+    printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
+    return -1;
+  }
+
+
+  /* declare monitoring program */
+  status=monitorDeclareMp( __FILE__ );
+  if (status!=0) {
+    printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
+    return -1;
+  }
+
+
+  /* define wait event timeout - 1s max */
+  monitorSetNowait();
+  monitorSetNoWaitNetworkTimeout(1000);
+  
+   /* init some counters */
+  int nevents_physics=0;
+  int nevents_total=0;
+
+  AliRawReader *rawReader = NULL;
+
+  AliPHOSRcuDA1 da1(2,-1); // DA1 (Calibration DA) for module2
+  
+  Float_t e[64][56][2];
+  Float_t t[64][56][2];
+
+  Int_t gain = -1;
+  Int_t X = -1;
+  Int_t Z = -1;
+
+  /* main loop (infinite) */
+  for(;;) {
+    struct eventHeaderStruct *event;
+    eventTypeType eventT;
+  
+    /* check shutdown condition */
+    if (daqDA_checkShutdown()) {break;}
+    
+    /* get next event (blocking call until timeout) */
+    status=monitorGetEventDynamic((void **)&event);
+    if (status==MON_ERR_EOF) {
+      printf ("End of File detected\n");
+      break; /* end of monitoring file has been reached */
+    }
+    
+    if (status!=0) {
+      printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
+      break;
+    }
+
+    /* retry if got no event */
+    if (event==NULL) {
+      continue;
+    }
+
+
+    /* use event - here, just write event id to result file */
+    eventT=event->eventType;
+    
+    if (eventT==PHYSICS_EVENT) {
+      
+      for(Int_t iX=0; iX<64; iX++) {
+       for(Int_t iZ=0; iZ<56; iZ++) {
+         for(Int_t iGain=0; iGain<2; iGain++) {
+           e[iX][iZ][iGain] = 0.;
+           t[iX][iZ][iGain] = 0.;
+         }
+       }
+      }
+
+      rawReader = new AliRawReaderDate((void*)event);
+//       AliPHOSRawDecoderv1 dc(rawReader,mapping);
+      AliPHOSRawDecoder dc(rawReader,mapping);
+      dc.SubtractPedestals(kTRUE);
+      
+      while(dc.NextDigit()) {
+
+       X = dc.GetRow() - 1;
+       Z = dc.GetColumn() - 1;
+
+       if(dc.IsLowGain()) gain = 0;
+       else
+         gain = 1;
+       
+       e[X][Z][gain] = dc.GetEnergy();
+       t[X][Z][gain] = dc.GetTime();
+       
+      }
+
+      da1.FillHistograms(e,t);
+      //da1.UpdateHistoFile();
+      
+      delete rawReader;     
+      nevents_physics++;
+    }
+    
+    nevents_total++;
+    
+    /* free resources */
+    free(event);
+    
+    /* exit when last event received, no need to wait for TERM signal */
+    if (eventT==END_OF_RUN) {
+      printf("EOR event detected\n");
+      break;
+    }
+  }
+  
+  for(Int_t i = 0; i < 4; i++) delete mapping[i];  
+  
+  return status;
+}