820623431fada7f09af86d82008e205e30d723fd
[u/mrichter/AliRoot.git] / PHOS / PHOSDA1da.cxx
1 /*
2
3 DAcase2.c
4
5 This program connects to the DAQ data source passed as argument
6 and populates local "./result.txt" file with the ids of events received
7 during the run.
8
9 The program exits when being asked to shut down (daqDA_checkshutdown)
10 or End of Run event.
11
12 Messages on stdout are exported to DAQ log system.
13
14 contact: alice-datesupport@cern.ch
15
16 */
17
18
19 #include "event.h"
20 #include "monitor.h"
21 extern "C" {
22 #include "daqDA.h"
23 }
24
25 #include <stdio.h>
26 #include <stdlib.h>
27
28 #include <TSystem.h>
29
30 #include "AliRawReader.h"
31 #include "AliRawReaderDate.h"
32 #include "AliPHOSDA1.h"
33 #include "AliPHOSRawDecoderv1.h"
34 #include "AliCaloAltroMapping.h"
35
36
37 /* Main routine
38       Arguments: 
39       1- monitoring data source
40 */
41 int main(int argc, char **argv) {
42
43   int status;
44   
45   if (argc!=2) {
46     printf("Wrong number of arguments\n");
47     return -1;
48   }
49
50
51   /* open result file */
52   FILE *fp=NULL;
53   fp=fopen("./result.txt","a");
54   if (fp==NULL) {
55     printf("Failed to open file\n");
56     return -1;
57   }
58
59   /* Open mapping files */
60   AliAltroMapping *mapping[4];
61   TString path = gSystem->Getenv("ALICE_ROOT");
62   path += "/PHOS/mapping/RCU";
63   TString path2;
64   for(Int_t i = 0; i < 4; i++) {
65     path2 = path;
66     path2 += i;
67     path2 += ".data";
68     mapping[i] = new AliCaloAltroMapping(path2.Data());
69   }
70   
71
72   /* define data source : this is argument 1 */  
73   status=monitorSetDataSource( argv[1] );
74   if (status!=0) {
75     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
76     return -1;
77   }
78
79
80   /* declare monitoring program */
81   status=monitorDeclareMp( __FILE__ );
82   if (status!=0) {
83     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
84     return -1;
85   }
86
87
88   /* define wait event timeout - 1s max */
89   monitorSetNowait();
90   monitorSetNoWaitNetworkTimeout(1000);
91   
92
93   /* log start of process */
94   printf("DA example case2 monitoring program started\n");  
95
96
97   /* init some counters */
98   int nevents_physics=0;
99   int nevents_total=0;
100
101   AliRawReader *rawReader = NULL;
102
103   AliPHOSDA1 da1(2); // DA1 (Calibration DA) for module2
104   
105   Float_t e[64][56][2];
106   Float_t t[64][56][2];
107
108   Int_t gain = -1;
109   Int_t X = -1;
110   Int_t Z = -1;
111
112   /* main loop (infinite) */
113   for(;;) {
114     struct eventHeaderStruct *event;
115     eventTypeType eventT;
116   
117     /* check shutdown condition */
118     if (daqDA_checkShutdown()) {break;}
119     
120     /* get next event (blocking call until timeout) */
121     status=monitorGetEventDynamic((void **)&event);
122     if (status==MON_ERR_EOF) {
123       printf ("End of File detected\n");
124       break; /* end of monitoring file has been reached */
125     }
126     
127     if (status!=0) {
128       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
129       break;
130     }
131
132     /* retry if got no event */
133     if (event==NULL) {
134       continue;
135     }
136
137
138     /* use event - here, just write event id to result file */
139     eventT=event->eventType;
140     
141     if (eventT==PHYSICS_EVENT) {
142       fprintf(fp,"Run #%lu, event size: %lu, BC:%u, Orbit:%u, Period:%u\n",
143           (unsigned long)event->eventRunNb,
144         (unsigned long)event->eventSize,
145         EVENT_ID_GET_BUNCH_CROSSING(event->eventId),
146         EVENT_ID_GET_ORBIT(event->eventId),
147         EVENT_ID_GET_PERIOD(event->eventId)
148       );
149       
150       for(Int_t iX=0; iX<64; iX++) {
151         for(Int_t iZ=0; iZ<56; iZ++) {
152           for(Int_t iGain=0; iGain<2; iGain++) {
153             e[iX][iZ][iGain] = 0.;
154             t[iX][iZ][iGain] = 0.;
155           }
156         }
157       }
158
159       rawReader = new AliRawReaderDate((void*)event);
160       AliPHOSRawDecoderv1 dc(rawReader,mapping);
161       dc.SubtractPedestals(kTRUE);
162       
163       while(dc.NextDigit()) {
164
165         X = dc.GetRow() - 1;
166         Z = dc.GetColumn() - 1;
167
168         if(dc.IsLowGain()) gain = 0;
169         else
170           gain = 1;
171         
172         e[X][Z][gain] = dc.GetEnergy();
173         t[X][Z][gain] = dc.GetTime();
174         
175       }
176
177       da1.FillHistograms(e,t);
178       //da1.UpdateHistoFile();
179       
180       delete rawReader;     
181       nevents_physics++;
182     }
183     
184     nevents_total++;
185     
186     /* free resources */
187     free(event);
188     
189     /* exit when last event received, no need to wait for TERM signal */
190     if (eventT==END_OF_RUN) {
191       printf("EOR event detected\n");
192       break;
193     }
194   }
195
196   for(Int_t i = 0; i < 4; i++) delete mapping[i];  
197
198   /* write report */
199   fprintf(fp,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);
200
201   /* close result file */
202   fclose(fp);
203
204
205   return status;
206 }