TPC/TPCCEda.cxx

   1 /*
   2 TPC DA for online calibration
   3
   4 Contact: Haavard.Helstrup@cern.ch
   5 Link:
   6 Run Type: PHYSICS_RUN STANDALONE DAQ
   7 DA Type: MON
   8 Number of events needed: 500
   9 Input Files:
  10 Output Files: tpcCE.root, to be exported to the DAQ FXS
  11 Trigger types used: PHYSICS_EVENT
  12
  13 */
  14
  15 /*
  16
  17 TPCCEda.cxx - calibration algorithm for TPC Central Electrode events
  18
  19 10/06/2007  sylvain.chapeland@cern.ch :  first version - clean skeleton based on DAQ DA case1
  20 06/12/2007  haavard.helstrup@cern.ch  :  created CE DA based on pulser code
  21
  22 contact: marian.ivanov@cern.ch
  23
  24
  25 This process reads RAW data from the files provided as command line arguments
  26 and save results in a file (named from RESULT_FILE define - see below).
  27
  28 */
  29
  30 #define RESULT_FILE "tpcCE.root"
  31
  32
  33 #include <daqDA.h>
  34 #include "event.h"
  35 #include "monitor.h"
  36 #include <stdio.h>
  37 #include <stdlib.h>
  38
  39 //
  40 //Root includes
  41 //
  42 #include <TFile.h>
  43 #include "TROOT.h"
  44 #include "TPluginManager.h"
  45 //
  46 //AliRoot includes
  47 //
  48 #include "AliRawReader.h"
  49 #include "AliRawReaderDate.h"
  50 #include "AliTPCRawStream.h"
  51 #include "AliTPCROC.h"
  52 #include "AliTPCCalROC.h"
  53 #include "AliTPCCalPad.h"
  54 #include "AliMathBase.h"
  55 #include "TTreeStream.h"
  56
  57 //
  58 // TPC calibration algorithm includes
  59 //
  60 #include "AliTPCCalibCE.h"
  61
  62
  63
  64
  65 /* Main routine
  66       Arguments: list of DATE raw data files
  67 */
  68 int main(int argc, char **argv) {
  69
  70  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
  71                                          "*",
  72                                          "TStreamerInfo",
  73                                          "RIO",
  74                                          "TStreamerInfo()");
  75
  76
  77   int i,status;
  78   AliTPCCalibCE calibCE;   // pedestal and noise calibration
  79
  80   if (argc<2) {
  81     printf("Wrong number of arguments\n");
  82     return -1;
  83   }
  84
  85
  86   /* log start of process */
  87   printf("TPC CE DA started - %s\n",__FILE__);
  88
  89
  90   /* set time bin range */
  91   calibCE.SetRangeTime(800,940);
  92
  93
  94   /* declare monitoring program */
  95   status=monitorDeclareMp( __FILE__ );
  96   if (status!=0) {
  97     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
  98     return -1;
  99   }
 100
 101
 102   /* loop over RAW data files */
 103   int nevents=0;
 104   for(i=1;i<argc;i++) {
 105
 106     /* define data source : this is argument i */
 107     printf("Processing file %s\n", argv[i]);
 108     status=monitorSetDataSource( argv[i] );
 109     if (status!=0) {
 110       printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
 111       return -1;
 112     }
 113
 114     /* read until EOF */
 115     while (true) {
 116       struct eventHeaderStruct *event;
 117
 118       /* check shutdown condition */
 119       if (daqDA_checkShutdown()) {break;}
 120
 121       /* get next event (blocking call until timeout) */
 122       status=monitorGetEventDynamic((void **)&event);
 123       if (status==MON_ERR_EOF) {
 124         printf ("End of File %d detected\n",i);
 125         break; /* end of monitoring file has been reached */
 126       }
 127
 128       if (status!=0) {
 129         printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
 130         break;
 131       }
 132
 133       /* retry if got no event */
 134       if (event==NULL) {
 135         continue;
 136       }
 137       nevents++;
 138
 139       //  Pulser calibration
 140
 141       AliRawReader *rawReader = new AliRawReaderDate((void*)event);
 142       calibCE.ProcessEvent(rawReader);
 143       delete rawReader;
 144
 145       /* free resources */
 146       free(event);
 147     }
 148   }
 149
 150   calibCE.Analyse();
 151   printf ("%d events processed\n",nevents);
 152
 153   TFile * fileTPC = new TFile (RESULT_FILE,"recreate");
 154   calibCE.Write("calibCE");
 155   delete fileTPC;
 156   printf("Wrote %s\n",RESULT_FILE);
 157
 158   /* store the result file on FES */
 159
 160   status=daqDA_FES_storeFile(RESULT_FILE,RESULT_FILE);
 161   if (status) {
 162     status = -2;
 163   }
 164
 165   return status;
 166 }