3e652f194d77bb5a0ff8f7b1a13c2db45ae852f1
[u/mrichter/AliRoot.git] / PHOS / PHOSLEDda.cxx
1 /*
2 contact: Boris.Polishchuk@cern.ch
3 link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx
4 reference run: /castor/cern.ch/alice/phos/2007/10/04/18/07000008249001.1000.root
5 run type: STANDALONE
6 DA type: MON 
7 number of events needed: 1000
8 input files: RCU0.data  RCU1.data  RCU2.data  RCU3.data
9 Output files: PHOS_Module2_LED.root
10 Trigger types used: CALIBRATION_EVENT
11 */
12
13
14 #include "event.h"
15 #include "monitor.h"
16 extern "C" {
17 #include "daqDA.h"
18 }
19
20 #include <stdio.h>
21 #include <stdlib.h>
22
23 #include <TSystem.h>
24 #include <TROOT.h>
25 #include <TPluginManager.h>
26
27 #include "AliRawReader.h"
28 #include "AliRawReaderDate.h"
29 #include "AliPHOSRcuDA1.h"
30 #include "AliPHOSRawDecoder.h"
31 #include "AliCaloAltroMapping.h"
32
33
34 /* Main routine
35       Arguments: 
36       1- monitoring data source
37 */
38 int main(int argc, char **argv) {
39
40   gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
41                                         "*",
42                                         "TStreamerInfo",
43                                         "RIO",
44                                         "TStreamerInfo()");
45
46   int status;
47   
48   if (argc!=2) {
49     printf("Wrong number of arguments\n");
50     return -1;
51   }
52
53   /* Retrieve mapping files from DAQ DB */ 
54   const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
55
56   for(Int_t iFile=0; iFile<4; iFile++) {
57     int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
58     if(failed) { 
59       printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]);
60       return -1;
61     }
62   }
63   
64   /* Open mapping files */
65   AliAltroMapping *mapping[4];
66   TString path = "./";
67   path += "RCU";
68   TString path2;
69   for(Int_t i = 0; i < 4; i++) {
70     path2 = path;
71     path2 += i;
72     path2 += ".data";
73     mapping[i] = new AliCaloAltroMapping(path2.Data());
74   }
75   
76
77   /* define data source : this is argument 1 */  
78   status=monitorSetDataSource( argv[1] );
79   if (status!=0) {
80     printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
81     return -1;
82   }
83
84
85   /* declare monitoring program */
86   status=monitorDeclareMp( __FILE__ );
87   if (status!=0) {
88     printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
89     return -1;
90   }
91
92
93   /* define wait event timeout - 1s max */
94   monitorSetNowait();
95   monitorSetNoWaitNetworkTimeout(1000);
96   
97    /* init some counters */
98   int nevents_physics=0;
99   int nevents_total=0;
100
101   Int_t fMod = 2; // module 2!
102   AliRawReader *rawReader = NULL;
103
104   AliPHOSRcuDA1 da1(fMod,-1,0); // DA1 for module2, no input/output file
105   
106   Float_t e[64][56][2];
107   Float_t t[64][56][2];
108
109   Int_t gain = -1;
110   Int_t X = -1;
111   Int_t Z = -1;
112   Int_t nFired = -1;
113
114   TH1I fFiredCells("fFiredCells","Number of fired cells per event",100,0,1000);
115
116   /* main loop (infinite) */
117   for(;;) {
118     struct eventHeaderStruct *event;
119     eventTypeType eventT;
120   
121     /* check shutdown condition */
122     if (daqDA_checkShutdown()) {break;}
123     
124     /* get next event (blocking call until timeout) */
125     status=monitorGetEventDynamic((void **)&event);
126     if (status==MON_ERR_EOF) {
127       printf ("End of File detected\n");
128       break; /* end of monitoring file has been reached */
129     }
130     
131     if (status!=0) {
132       printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
133       break;
134     }
135
136     /* retry if got no event */
137     if (event==NULL) {
138       continue;
139     }
140
141
142     /* use event - here, just write event id to result file */
143     eventT=event->eventType;
144     
145     if (eventT==PHYSICS_EVENT) {
146       
147       for(Int_t iX=0; iX<64; iX++) {
148         for(Int_t iZ=0; iZ<56; iZ++) {
149           for(Int_t iGain=0; iGain<2; iGain++) {
150             e[iX][iZ][iGain] = 0.;
151             t[iX][iZ][iGain] = 0.;
152           }
153         }
154       }
155
156       nFired = 0;
157
158       rawReader = new AliRawReaderDate((void*)event);
159 //       AliPHOSRawDecoderv1 dc(rawReader,mapping);
160       AliPHOSRawDecoder 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         if(gain && dc.GetEnergy()>40)
176           nFired++;
177
178       }
179       
180       da1.FillHistograms(e,t);
181       fFiredCells.Fill(nFired);
182
183       delete rawReader;     
184       nevents_physics++;
185     }
186     
187     nevents_total++;
188     
189     /* free resources */
190     free(event);
191     
192     /* exit when last event received, no need to wait for TERM signal */
193     if (eventT==END_OF_RUN) {
194       printf("EOR event detected\n");
195       break;
196     }
197   }
198   
199   for(Int_t i = 0; i < 4; i++) delete mapping[i];  
200   
201   /* Be sure that all histograms are saved */
202
203   char localfile[128];
204   sprintf(localfile,"PHOS_Module%d_LED.root",fMod);
205   TFile* f = new TFile(localfile,"recreate");
206   
207   const TH2F* h2=0;
208   const TH1F* h1=0;
209
210   Int_t nGood=0;    // >10 entries in peak
211   Int_t nMax=-111;  // max. number of entries in peak
212   Int_t iXmax=-1;
213   Int_t iZmax=-1;
214
215   for(Int_t iX=0; iX<64; iX++) {
216     for(Int_t iZ=0; iZ<56; iZ++) {
217       
218       h1 = da1.GetHgLgRatioHistogram(iX,iZ); // High Gain/Low Gain ratio
219       if(h1) {
220         if(h1->GetMaximum()>10.) nGood++;
221         if(h1->GetMaximum()>nMax) {nMax = (Int_t)h1->GetMaximum(); iXmax=iX; iZmax=iZ;}
222         h1->Write(); 
223       }
224       
225       for(Int_t iGain=0; iGain<2; iGain++) {
226         h2 = da1.GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy
227         if(h2) h2->Write();
228       }
229       
230     }
231   }
232   
233   fFiredCells.Write();
234   f->Close();
235   
236   /* Store output files to the File Exchange Server */
237   daqDA_FES_storeFile("PHOS_Module2_LED.root","LED");
238   
239   printf("%d physics events of %d total processed.\n",nevents_physics,nevents_total);
240   printf("%d histograms has >10 entries in maximum, max. is %d entries ",nGood,nMax);
241   printf("(module,iX,iZ)=(%d,%d,%d)",fMod,iXmax,iZmax);
242   printf("\n");
243
244   return status;
245 }