remove exec
[u/mrichter/AliRoot.git] / EMCAL / EMCALLEDda.cxx
1 /*
2   EMCAL DA for online calibration: for LED studies
3   
4   Contact: silvermy@ornl.gov
5   Run Type: PHYSICS or STANDALONE
6   DA Type: MON 
7   Number of events needed: continously accumulating for all runs, rate ~0.1-1 Hz
8   Input Files: argument list
9   Output Files: RESULT_FILE=EMCALLED.root, to be exported to the DAQ FXS
10   fileId:  FILE_ID=EMCALLED    
11   Trigger types used: CALIBRATION_EVENT 
12 */
13 /*
14   This process reads RAW data from the files provided as command line arguments
15   and save results (class itself) in a file (named from RESULT_FILE define - 
16   see below).
17 */
18
19 #define RESULT_FILE  "EMCALLED.root"
20 #define FILE_ID "signal"
21 #define AliDebugLevel() -1
22 #define FILE_SIGClassName "emcCalibSignal"
23 const int kNRCU = 4;
24 /* LOCAL_DEBUG is used to bypass daq* calls that do not work locally */
25 //#define LOCAL_DEBUG 1 // comment out to run normally                                                            
26 extern "C" {
27 #include <daqDA.h>
28 }
29 #include "event.h" /* in $DATE_COMMON_DEFS/; includes definition of event types */
30 #include "monitor.h" /* in $DATE_MONITOR_DIR/; monitor* interfaces */
31
32 #include "stdio.h"
33 #include "stdlib.h"
34
35 // ROOT includes
36 #include <TFile.h> 
37 #include <TROOT.h> 
38 #include <TPluginManager.h> 
39 #include <TSystem.h> 
40
41 //
42 //AliRoot includes
43 //
44 #include "AliRawReader.h"
45 #include "AliRawReaderDate.h"
46 #include "AliRawEventHeaderBase.h"
47 #include "AliCaloRawStreamV3.h"
48 #include "AliCaloAltroMapping.h"
49 #include "AliLog.h"
50 #include "AliDAQ.h"
51
52 //
53 // EMC calibration-helper algorithm includes
54 //
55 #include "AliCaloCalibSignal.h"
56
57 /*
58   Main routine, EMC signal detector algorithm 
59   Arguments: list of DATE raw data files
60 */
61
62 int main(int argc, char **argv) { // Main routine, EMC signal detector algorithm 
63
64   AliLog::SetClassDebugLevel("AliCaloRawStreamV3",-5);
65   AliLog::SetClassDebugLevel("AliRawReaderDate",-5);
66   AliLog::SetModuleDebugLevel("RAW",-5);
67
68   if (argc<2) {
69     printf("Wrong number of arguments\n");
70     return -1;  
71   }
72
73   /* magic line - for TStreamerInfo */
74   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
75                                         "*",
76                                         "TStreamerInfo",
77                                         "RIO",
78                                         "TStreamerInfo()"); 
79   
80   int i, status;
81   
82   /* log start of process */
83   printf("EMCAL DA started - %s\n",__FILE__);
84   
85   Int_t emcID = AliDAQ::DetectorID("EMCAL"); // bit 18..
86
87   /* declare monitoring program */
88   status=monitorDeclareMp( __FILE__ );
89   if (status!=0) {
90     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
91     return -1;
92   }
93   /* define wait event timeout - 1s max */
94   monitorSetNowait();
95   monitorSetNoWaitNetworkTimeout(1000);
96   
97   /* Retrieve mapping files from DAQ DB */ 
98   const char* mapFiles[kNRCU] = {"RCU0A.data","RCU1A.data","RCU0C.data","RCU1C.data"};
99   
100   for(Int_t iFile=0; iFile<kNRCU; iFile++) {
101     int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
102     if(failed) { 
103       printf("Cannot retrieve file %d : %s from DAQ DB. Exit now..\n",
104              iFile, mapFiles[iFile]);
105 #ifdef LOCAL_DEBUG
106 #else
107       return -1;
108 #endif
109     }
110   }
111   
112   /* Open mapping files */
113   AliCaloAltroMapping *mapping[kNRCU];
114   TString path = "./";
115   path += "RCU";
116   TString path2;
117   TString side[] = {"A","C"};//+ and - pseudorapidity supermodules
118   for(Int_t j = 0; j < 2; j++){
119     for(Int_t i = 0; i < 2; i++) {
120       path2 = path;
121       path2 += i;
122       path2 += side[j]; 
123       path2 += ".data";
124       mapping[j*2 + i] = new AliCaloAltroMapping(path2.Data());
125     }
126   }
127   
128   /* Retrieve cut=parameter file from DAQ DB */ 
129   const char* parameterFile = {"EMCALLEDda.dat"};
130   
131   int failed = daqDA_DB_getFile(parameterFile, parameterFile);
132   if (failed) { 
133     printf("Cannot retrieve file : %s from DAQ DB. Exit now..\n",
134            parameterFile);
135 #ifdef LOCAL_DEBUG
136 #else
137     return -1;
138 #endif
139   }
140   
141   /* set up our analysis classes */  
142   AliCaloCalibSignal * calibSignal = new AliCaloCalibSignal(AliCaloCalibSignal::kEmCal); 
143   calibSignal->SetAltroMapping( mapping );
144   calibSignal->SetParametersFromFile( parameterFile );
145   
146   AliRawReader *rawReader = NULL;
147   int nevents=0;
148   
149   /* loop over RAW data files */
150   for ( i=1; i<argc; i++ ) {
151     
152     /* define data source : this is argument i */
153     printf("Processing file %s\n", argv[i]);
154     status=monitorSetDataSource( argv[i] );
155     if (status!=0) {
156       printf("monitorSetDataSource() failed. Error=%s. Exiting ...\n", monitorDecodeError(status));
157       return -1;
158     }
159     
160     /* read until EOF */
161     struct eventHeaderStruct *event;
162     eventTypeType eventT;
163
164     for ( ; ; ) { // infinite loop
165       
166       /* check shutdown condition */
167       if (daqDA_checkShutdown()) {break;}
168       
169       /* get next event (blocking call until timeout) */
170       status=monitorGetEventDynamic((void **)&event);
171       if (status==MON_ERR_EOF) {
172         printf ("End of File %d (%s) detected\n", i, argv[i]);
173         break; /* end of monitoring file has been reached */
174       }
175       if (status!=0) {
176         printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
177         break;
178       }
179       
180       /* retry if got no event */
181       if (event==NULL) {
182         continue;
183       }
184       eventT = event->eventType; /* just convenient shorthand */
185       /* only look at calibration events */
186       if ( eventT != calibrationEvent ) {
187         free(event);    
188         continue;
189       }
190       
191       /* only look at events where EMCAL was included */
192       if (! TEST_DETECTOR_IN_PATTERN(event->eventDetectorPattern, emcID) ) {
193         free(event);    
194         continue;
195       }
196
197       nevents++; // count how many acceptable events we have
198       
199       //  Signal calibration
200       rawReader = new AliRawReaderDate((void*)event);
201       calibSignal->SetRunNumber(event->eventRunNb);
202       calibSignal->ProcessEvent(rawReader);
203       delete rawReader;
204       
205       /* free resources */
206       free(event);    
207       
208     } //until EOF
209   } // loop over files
210   
211   // calculate average values also, for the LED info
212   calibSignal->SetUseAverage(kTRUE);
213   calibSignal->Analyze();
214   
215   // by default, we only save the full info in debug mode  
216 #ifdef LOCAL_DEBUG
217 #else
218   // reset the full trees, when we are not in debug mode
219   calibSignal->GetTreeAmpVsTime()->Reset();
220   calibSignal->GetTreeLEDAmpVsTime()->Reset();
221 #endif
222   
223   //
224   // write class to rootfile
225   //
226   
227   printf ("%d calibration events processed.\n",nevents);
228   
229   TFile f(RESULT_FILE, "recreate");
230   if (!f.IsZombie()) { 
231     f.cd();
232     calibSignal->Write(FILE_SIGClassName);
233     f.Close();
234     printf("Objects saved to file \"%s\" as \"%s\".\n", 
235            RESULT_FILE, FILE_SIGClassName); 
236   } 
237   else {
238     printf("Could not save the object to file \"%s\".\n", 
239            RESULT_FILE);
240   }
241   
242   //
243   // closing down; see if we can delete our analysis helper(s) also
244   //
245   delete calibSignal;
246   for(Int_t iFile=0; iFile<kNRCU; iFile++) {
247     if (mapping[iFile]) delete mapping[iFile];
248   }
249   
250   /* store the result file on FES */
251 #ifdef LOCAL_DEBUG
252 #else
253   status = daqDA_FES_storeFile(RESULT_FILE, FILE_ID);
254 #endif
255   
256   return status;
257 }