bug fix in GetRefFilePrefix
[u/mrichter/AliRoot.git] / PHOS / PHOSGAINda.cxx
CommitLineData
62adef6c 1/*
2contact: Boris.Polishchuk@cern.ch
3link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx
4reference run: /castor/cern.ch/alice/phos/2007/10/04/18/07000008249001.1000.root
82127abf 5run type: PHYSICS
62adef6c 6DA type: MON
7number of events needed: 1000
8input files: RCU0.data RCU1.data RCU2.data RCU3.data
9Output files: PHOS_Module2_Calib.root
82127abf 10Trigger types used: PHYSICS
62adef6c 11*/
12
13
14#include "event.h"
15#include "monitor.h"
16extern "C" {
17#include "daqDA.h"
18}
19
20#include <stdio.h>
21#include <stdlib.h>
22
23#include <TSystem.h>
c59c9912 24#include <TROOT.h>
25#include <TPluginManager.h>
62adef6c 26
27#include "AliRawReader.h"
28#include "AliRawReaderDate.h"
29#include "AliPHOSRcuDA1.h"
8083e819 30#include "AliPHOSRawFitterv0.h"
62adef6c 31#include "AliCaloAltroMapping.h"
8083e819 32#include "AliCaloRawStreamV3.h"
62adef6c 33
34
35/* Main routine
36 Arguments:
37 1- monitoring data source
38*/
39int main(int argc, char **argv) {
40
c59c9912 41 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
42 "*",
43 "TStreamerInfo",
44 "RIO",
45 "TStreamerInfo()");
46
62adef6c 47 int status;
48
49 if (argc!=2) {
50 printf("Wrong number of arguments\n");
51 return -1;
52 }
a9b2b02d 53
54 /* Retrieve ZS parameters from DAQ DB */
55 const char* zsfile = "zs.txt";
56 int failZS = daqDA_DB_getFile(zsfile, zsfile);
57
58 Int_t offset,threshold;
59
60 if(!failZS) {
61 FILE *f = fopen(zsfile,"r");
62 int scan = fscanf(f,"%d %d",&offset,&threshold);
63 }
64
c59c9912 65 /* Retrieve mapping files from DAQ DB */
66 const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
a9b2b02d 67
c59c9912 68 for(Int_t iFile=0; iFile<4; iFile++) {
69 int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
70 if(failed) {
71 printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]);
72 return -1;
73 }
74 }
62adef6c 75
76 /* Open mapping files */
77 AliAltroMapping *mapping[4];
78 TString path = "./";
79 path += "RCU";
80 TString path2;
81 for(Int_t i = 0; i < 4; i++) {
82 path2 = path;
83 path2 += i;
84 path2 += ".data";
85 mapping[i] = new AliCaloAltroMapping(path2.Data());
86 }
87
62adef6c 88 /* define data source : this is argument 1 */
89 status=monitorSetDataSource( argv[1] );
90 if (status!=0) {
91 printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
92 return -1;
93 }
94
95
96 /* declare monitoring program */
97 status=monitorDeclareMp( __FILE__ );
98 if (status!=0) {
99 printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
100 return -1;
101 }
102
103
104 /* define wait event timeout - 1s max */
105 monitorSetNowait();
106 monitorSetNoWaitNetworkTimeout(1000);
107
108 /* init some counters */
109 int nevents_physics=0;
110 int nevents_total=0;
111
112 AliRawReader *rawReader = NULL;
a9b2b02d 113 AliPHOSRcuDA1* dAs[5];
62adef6c 114
a9b2b02d 115 for(Int_t iMod=0; iMod<5; iMod++) {
116 dAs[iMod] = 0;
117 }
62adef6c 118
119 Float_t e[64][56][2];
120 Float_t t[64][56][2];
121
a9b2b02d 122 for(Int_t iX=0; iX<64; iX++) {
123 for(Int_t iZ=0; iZ<56; iZ++) {
124 for(Int_t iGain=0; iGain<2; iGain++) {
125 e[iX][iZ][iGain] = 0.;
126 t[iX][iZ][iGain] = 0.;
127 }
128 }
129 }
130
8083e819 131 Int_t cellX = -1;
132 Int_t cellZ = -1;
133 Int_t nBunches = 0;
134 Int_t sigStart, sigLength;
a9b2b02d 135 Int_t caloFlag;
62adef6c 136
137 /* main loop (infinite) */
138 for(;;) {
139 struct eventHeaderStruct *event;
140 eventTypeType eventT;
141
142 /* check shutdown condition */
143 if (daqDA_checkShutdown()) {break;}
144
145 /* get next event (blocking call until timeout) */
146 status=monitorGetEventDynamic((void **)&event);
147 if (status==MON_ERR_EOF) {
148 printf ("End of File detected\n");
149 break; /* end of monitoring file has been reached */
150 }
151
152 if (status!=0) {
153 printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
154 break;
155 }
156
157 /* retry if got no event */
158 if (event==NULL) {
159 continue;
160 }
161
162
163 /* use event - here, just write event id to result file */
164 eventT=event->eventType;
165
166 if (eventT==PHYSICS_EVENT) {
167
62adef6c 168 rawReader = new AliRawReaderDate((void*)event);
8083e819 169 AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping);
a9b2b02d 170 AliPHOSRawFitterv0 fitter;
8083e819 171 fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS
62adef6c 172
a9b2b02d 173 if(!failZS) {
174 fitter.SubtractPedestals(kFALSE);
175 fitter.SetAmpOffset(offset);
176 fitter.SetAmpThreshold(threshold);
177 }
178
8083e819 179 while (stream.NextDDL()) {
180 while (stream.NextChannel()) {
181
182 cellX = stream.GetCellX();
183 cellZ = stream.GetCellZ();
184 caloFlag = stream.GetCaloFlag(); // 0=LG, 1=HG, 2=TRU
2905f08d 185
186 if(caloFlag!=0 && caloFlag!=1) continue; //TRU data!
187
a9b2b02d 188 // In case of oscillating signals with ZS,
189 //a channel can have several bunches.
190
8083e819 191 nBunches = 0;
192 while (stream.NextBunch()) {
193 nBunches++;
194 if (nBunches > 1) continue;
a9b2b02d 195 sigStart = stream.GetStartTimeBin();
196 sigLength = stream.GetBunchLength();
92236b27 197 fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag);
198 fitter.Eval(stream.GetSignals(),sigStart,sigLength);
8083e819 199 } // End of NextBunch()
200
a9b2b02d 201 if (nBunches>1) continue;
8083e819 202
203 e[cellX][cellZ][caloFlag] = fitter.GetEnergy();
204 t[cellX][cellZ][caloFlag] = fitter.GetTime();
205 }
2905f08d 206
207 if(stream.GetModule()<0 || stream.GetModule()>4) continue;
208
a9b2b02d 209 if(dAs[stream.GetModule()])
210 dAs[stream.GetModule()]->FillHistograms(e,t);
2905f08d 211 else {
a9b2b02d 212 dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1);
2905f08d 213 dAs[stream.GetModule()]->FillHistograms(e,t);
214 }
a9b2b02d 215
216 for(Int_t iX=0; iX<64; iX++) {
217 for(Int_t iZ=0; iZ<56; iZ++) {
218 for(Int_t iGain=0; iGain<2; iGain++) {
219 e[iX][iZ][iGain] = 0.;
220 t[iX][iZ][iGain] = 0.;
221 }
222 }
223 }
224
62adef6c 225 }
226
a9b2b02d 227// da1.FillHistograms(e,t);
228// //da1.UpdateHistoFile();
62adef6c 229
a9b2b02d 230 delete rawReader;
231 nevents_physics++;
62adef6c 232 }
233
234 nevents_total++;
235
236 /* free resources */
237 free(event);
238
239 /* exit when last event received, no need to wait for TERM signal */
240 if (eventT==END_OF_RUN) {
241 printf("EOR event detected\n");
242 break;
243 }
244 }
245
246 for(Int_t i = 0; i < 4; i++) delete mapping[i];
247
c59c9912 248 /* Be sure that all histograms are saved */
c59c9912 249
c59c9912 250 char localfile[128];
251
252 for(Int_t iMod=0; iMod<5; iMod++) {
a9b2b02d 253 if(!dAs[iMod]) continue;
254
255 dAs[iMod]->UpdateHistoFile();
256 dAs[iMod]->SetWriteToFile(kFALSE);
257
258 /* Store output files to the File Exchange Server */
c59c9912 259 sprintf(localfile,"PHOS_Module%d_Calib.root",iMod);
604fe52e 260 daqDA_FES_storeFile(localfile,"AMPLITUDES");
c59c9912 261 }
262
62adef6c 263 return status;
264}