isolated cell search algo is added
[u/mrichter/AliRoot.git] / TOF / TOFnoiseda.cxx
CommitLineData
29c170aa 1/*
2
3TOF DA for online calibration from pulser data
4
5Contact: Chiara.Zampolli@bo.infn.it
6Link: www.bo.infn.it/~zampolli
7fffa85b 7Run Type: NOISE
29c170aa 8DA Type: LDC
9Number of events needed: 10000
10Input Files: TOF<nrun>.raw, where <nrun> is the run number
11Output Files: TOFoutPulserLDC_<LDCid>.root, where <LDCid> is the id of the LDC which is read (2 characters field, e.g. TOFoutPulserLDC_03.root),to be exported to the DAQ FXS
12Trigger types used: PHYSICS_EVENT (for the time being)
13
14*/
15
16// DATE
17#include "event.h"
18#include "monitor.h"
89b3ea8c 19#include "daqDA.h"
29c170aa 20
21#include <stdio.h>
22#include <stdlib.h>
23
24//AliRoot
25#include <AliTOFRawStream.h>
26#include <AliRawReaderDate.h>
27#include <AliRawReader.h>
28#include <AliTOFGeometry.h>
29#include <AliDAQ.h>
30#include <AliTOFHitData.h>
31#include <AliTOFHitDataBuffer.h>
32#include <AliTOFDecoder.h>
fa2852b6 33#include <AliTOFNoiseConfigHandler.h>
29c170aa 34
35//ROOT
36#include <TFile.h>
37#include <TKey.h>
38#include <TH1I.h>
39#include <TObject.h>
40#include <TMath.h>
41#include <TSystem.h>
3bc66139 42#include "TROOT.h"
43#include "TPluginManager.h"
fa2852b6 44#include "TSAXParser.h"
29c170aa 45
46/* Main routine
47 Arguments: list of DATE raw data files
48*/
49int main(int argc, char **argv) {
50
3bc66139 51/* magic line from Rene */
52 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
53 "*",
54 "TStreamerInfo",
55 "RIO",
56 "TStreamerInfo()");
57
29c170aa 58 AliTOFGeometry * geom = new AliTOFGeometry();
59
60 static const Int_t size = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
61 TH1F::AddDirectory(0);
62 TH1F * htofNoise = new TH1F("hTOFnoise","histo with signals on TOF during noise", size,-0.5,size-0.5);
63 for (Int_t ibin =1;ibin<=size;ibin++){
64 htofNoise->SetBinContent(ibin,-1);
65 }
66 UInt_t ldcId=0;
67
68 int status;
69
70 /* log start of process */
71 printf("TOF DA started\n");
72
73 /* check that we got some arguments = list of files */
74 if (argc<2) {
75 printf("Wrong number of arguments\n");
76 return -1;
77 }
78
79 /* open result file */
80 FILE *fp=NULL;
81 fp=fopen("./result.txt","a");
82 if (fp==NULL) {
83 printf("Failed to open file\n");
84 return -1;
85 }
86
fa2852b6 87 /* retrieve config file */
88 int getConfigFile = daqDA_DB_getFile("TOFNoiseConfig.xml","TOFNoiseConfig.xml");
89 if (getConfigFile != 0){
90 printf("Failed to retrieve config file from DB! returning...\n");
91 return -1;
92 }
93
94 AliTOFNoiseConfigHandler* tofHandler = new AliTOFNoiseConfigHandler();
95 TSAXParser *parser = new TSAXParser();
96 parser->ConnectToHandler("AliTOFNoiseConfigHandler", tofHandler);
46c74706 97 if (parser->ParseFile("./TOFNoiseConfig.xml") != 0) {
98 printf("Failed parsing config file! retunring... \n");
99 return -1;
fa2852b6 100 }
013ac2dd 101
fa2852b6 102 Int_t debugFlag = tofHandler->GetDebugFlag();
013ac2dd 103
fa2852b6 104 printf("the debug flag is %i\n",debugFlag);
fa2852b6 105
29c170aa 106 /* init some counters */
107 int nevents_physics=0;
108 int nevents_total=0;
109
46c74706 110
111 Int_t nPDBEntriesToT = 0;
112 Int_t nDBEntriesToT = 0;
113 AliTOFHitData *HitData;
114 Int_t dummy = -1;
115 Int_t Volume[5];
116 for (Int_t i=0;i<5;i++) Volume[i]=-1;
117 AliTOFRawStream *rawStreamTOF = new AliTOFRawStream();
118 const AliRawDataHeader *currentCDH;
119 AliTOFDecoder * decoderTOF = new AliTOFDecoder();
120 AliTOFHitDataBuffer *DataBuffer[AliDAQ::NumberOfDdls("TOF")];
121 AliTOFHitDataBuffer *PackedDataBuffer[AliDAQ::NumberOfDdls("TOF")];
122 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++) {
123 DataBuffer[i]=new AliTOFHitDataBuffer();
124 PackedDataBuffer[i]=new AliTOFHitDataBuffer();
125 }
126 Int_t currentEquipment;
127 Int_t currentDDL;
128 UChar_t *data = 0x0;
129
29c170aa 130 /* read the data files */
131 int n;
132 for (n=1;n<argc;n++) {
133
134 status=monitorSetDataSource( argv[n] );
135 if (status!=0) {
136 printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
137 return -1;
138 }
139
140 /* read the file */
141 for(;;) {
142 struct eventHeaderStruct *event;
143 eventTypeType eventT;
144
145 /* get next event */
146 status=monitorGetEventDynamic((void **)&event);
147 if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */
148 if (status!=0) {
149 printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
150 return -1;
151 }
152
153 /* retry if got no event */
154 if (event==NULL) {
155 break;
156 }
157
158 /* use event - here, just write event id to result file */
159 eventT=event->eventType;
160
161 if (eventT==PHYSICS_EVENT) {
162 //printf ("event %i \n", nevents_physics);
46c74706 163
164 nPDBEntriesToT = 0;
165 nDBEntriesToT = 0;
166 dummy = -1;
29c170aa 167 AliRawReader *rawReader = new AliRawReaderDate((void*)event);
46c74706 168 rawStreamTOF->SetRawReader(rawReader);
169 data = 0x0;
170
29c170aa 171 while(rawReader->ReadHeader()){
29c170aa 172 ldcId = rawReader->GetLDCId();
46c74706 173 //debugging printings
174 //if (debugFlag) printf ("ldcId = %i \n",ldcId);
29c170aa 175
176 //memory leak prevention (actually data should be always 0x0 here)
177 if (data != 0x0) delete [] data;
178
179 //get equipment infos
180 currentEquipment = rawReader->GetEquipmentId();
181 currentDDL = rawReader->GetDDLID();
9512a964 182 currentCDH = rawReader->GetDataHeader();
29c170aa 183 Int_t nchDDL = 0;
184 if (currentDDL%2==0) {
185 nchDDL = 2160;
186 }
187 else {
188 nchDDL = 2208;
189 }
46c74706 190
29c170aa 191 Int_t * array = new Int_t[nchDDL];
46c74706 192
193 for (Int_t ii=0; ii<nchDDL; ii++) array[ii] = 0;
29c170aa 194 decoderTOF->GetArrayDDL(array, currentDDL);
195
46c74706 196 for (Int_t i=0;i<nchDDL;i++) {
29c170aa 197 if (htofNoise->GetBinContent(array[i]+1)<0) htofNoise->SetBinContent(array[i]+1,0);
198 }
199
46c74706 200 //debugging printings
201 // if (debugFlag) printf(" Equipment = %i, and DDL = %i \n", currentEquipment,currentDDL);
29c170aa 202 const Int_t kDataSize = rawReader->GetDataSize();
203 const Int_t kDataWords = kDataSize / 4;
204 data = new UChar_t[kDataSize];
205 decoderTOF->SetDataBuffer(DataBuffer[currentDDL]);
206 decoderTOF->SetPackedDataBuffer(PackedDataBuffer[currentDDL]);
207 //start decoding
208 if (!rawReader->ReadNext(data, kDataSize)){
209 rawReader->AddMajorErrorLog(AliTOFRawStream::kDDLdataReading);
210 printf("Error while reading DDL data. Go to next equipment \n");
211 delete [] data;
212 data = 0x0;
213 continue;
214 }
46c74706 215 if (decoderTOF->Decode((UInt_t *)data, kDataWords, currentCDH) == kTRUE ) {
29c170aa 216 rawReader->AddMajorErrorLog(AliTOFRawStream::kDDLDecoder,Form("DDL # = %d",currentDDL));
217 printf("Error while decoding DDL # %d: decoder returned with errors \n", currentDDL);
218 }
219
220 Int_t nDBEntries = DataBuffer[currentDDL]->GetEntries();
221 Int_t nPDBEntries = PackedDataBuffer[currentDDL]->GetEntries();
222 nPDBEntriesToT+=nPDBEntries;
223 nDBEntriesToT+=nDBEntries;
224 /* reset buffer */
225 DataBuffer[currentDDL]->Reset();
226
227 /* read data buffer hits */
228 for (Int_t iHit = 0; iHit < nPDBEntries; iHit++){
229 HitData = PackedDataBuffer[currentDDL]->GetHit(iHit);
230 /* add volume information */
231 HitData->SetDDLID(currentDDL);
46c74706 232 for (Int_t i=0;i<5;i++) Volume[i]=-1;
29c170aa 233 rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
234 if (Volume[0]==-1 ||
235 Volume[1]==-1 ||
236 Volume[2]==-1 ||
237 Volume[3]==-1 ||
238 Volume[4]==-1) continue;
239 else {
240 dummy = Volume[3];
241 Volume[3] = Volume[4];
242 Volume[4] = dummy;
243 Int_t index = geom->GetIndex(Volume);
29c170aa 244 // to check array indexes
245 /*
013ac2dd 246 Bool_t found =kFALSE;
29c170aa 247 for (Int_t j=0;j<nchDDL;j++){
248 if (index==array[j]) {
249 found = kTRUE;
250 break;
251 }
252 }
253 printf ("index = %i, found = %i\n",index, (Int_t)found);
254 */
255 //printf ("index = %i \n",index);
256 htofNoise->Fill(index); //channel index start from 0, bin index from 1
257 //debugging printings
46c74706 258 // if (debugFlag) printf("sector %i, plate %i, strip %i, padz %i, padx %i \n",Volume[0],Volume[1],Volume[2],Volume[3],Volume[4]); // too verbose
29c170aa 259 }
260 }
261 /* reset buffer */
262 PackedDataBuffer[currentDDL]->Reset();
263 delete [] data;
264 data = 0x0;
265 delete [] array;
266 }
267
46c74706 268 //debugging printings
269 //if (debugFlag) {
270 // printf(" Packed Hit Buffer Entries = %i \n",nPDBEntriesToT);
271 // printf(" Hit Buffer Entries = %i \n",nDBEntriesToT);
272 // }
29c170aa 273
274 delete rawReader;
275 rawReader = 0x0;
276 nevents_physics++;
277 }
278 nevents_total++;
279
280 /* free resources */
281 free(event);
282 }
283
284 }
46c74706 285
286 delete decoderTOF;
287 decoderTOF=0x0;
288
289 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
290 delete DataBuffer[i];
291 delete PackedDataBuffer[i];
292 }
293
294 delete rawStreamTOF;
295 rawStreamTOF = 0x0;
296
29c170aa 297 delete geom;
298 geom = 0x0;
299
fa2852b6 300 Int_t noisyChannels = 0;
301 Int_t checkedChannels = 0;
013ac2dd 302 Float_t time = nevents_physics; // acquisition time in number of events
46c74706 303 //debugging printings
304 if (debugFlag) printf(" Noise run lasted %f s \n",time);
29c170aa 305 for (Int_t ibin =1;ibin<=size;ibin++){
306 Float_t cont = htofNoise->GetBinContent(ibin);
307 if (cont!=-1) {
fa2852b6 308 checkedChannels++;
46c74706 309 //debugging printings
310 if (debugFlag) printf(" content = %f \n", cont);
29c170aa 311 htofNoise->SetBinContent(ibin,cont/time);
46c74706 312 //debugging printings
313 if (debugFlag) printf(" scaled content = %f \n", cont/time);
fa2852b6 314 if (cont != 0){
315 noisyChannels++;
316 }
29c170aa 317 }
318 }
319
46c74706 320 if (debugFlag) {
321 printf("Number of checked channels = %i\n",checkedChannels);
322 printf("Number of noisy channels (before applying threshold) = %i\n",noisyChannels);
fa2852b6 323 }
324
29c170aa 325 //write the Run level file
326 char filename[100];
327 sprintf(filename,"TOFoutNoiseLDC_%02i.root",ldcId);
328 TFile * fileRun = new TFile (filename,"RECREATE");
329 htofNoise->Write();
330 fileRun->Close();
331
332 /* write report */
333 fprintf(fp,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);
334
335 /* close result file */
336 fclose(fp);
337
338 status =0;
339
340 /* store the result file on FES */
98a837bf 341 status=daqDA_FES_storeFile(filename,"NOISE");
29c170aa 342 if (status) {
343 status = -2;
344 }
345
346 return status;
347}