]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ITS/ITSSPDSCANda.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
IMPORTANT BUG FIXES svn diff RAW Causing a segv in AliMDC::ClassName method when...
[u/mrichter/AliRoot.git] / ITS / ITSSPDSCANda.cxx
CommitLineData
6727e2db 1/*
6e2aae15 2- Contact: - henrik.tydesjo@cern.ch
3- Link: - http://tydes.home.cern.ch/tydes/doc/CalibrationOverview/CalibrationAlgorithms/
4- Run Type: - DAQ_MIN_TH_SCAN,DAQ_MEAN_TH_SCAN,DAQ_UNIFORMITY_SCAN,DAQ_NOISY_PIX_SCAN,DAQ_PIX_DELAY_SCAN,DAQ_FO_UNIF_SCAN
5- DA Type: - LDC
6- Number of events needed:
7- Input Files: - daq db config files: spd_standal_params,spd_perm_noisy , previous noisy lists: ./calibResults/ScanNoisy/* , raw data
8- Output Files: - Ref Data: ./calibResults/ScanReference/* , DCS update: ./calibResults/ScanDCSconfigToFXS/* , OCDB update: ./calibResults/ScanNoisyToFXS/* , persistent files: ./calibResults/ScanNoisy/*
9- Trigger types used:
6727e2db 10*/
11
1fd93b67 12////////////////////////////////////////////////////////////////////////////////
6727e2db 13// This program can be compiled in two modes. //
1fd93b67 14// //
15// 1. With the DAQ DA framework on. This is the default operating mode. //
16// Call this program with the name of the executable followed by the //
17// data files to process. //
18// //
19// 2. Without the DAQ DA framework on. Define the SPD_DA_OFF environment var. //
20// Call this program with the name of the executable followed by the //
21// runNr and the data files to process. //
22// //
23////////////////////////////////////////////////////////////////////////////////
3ec5024c 24
3f917533 25#ifndef SPD_DA_OFF
3ec5024c 26extern "C" {
27#include "daqDA.h"
28}
3f917533 29#endif
3ec5024c 30#include "event.h"
31#include "monitor.h"
32#include "AliRawReaderDate.h"
33#include "AliITSRawStreamSPD.h"
34#include "AliITSOnlineSPDscan.h"
35#include "AliITSOnlineSPDscanSingle.h"
36#include "AliITSOnlineSPDscanMultiple.h"
37#include "AliITSOnlineSPDscanMeanTh.h"
38#include "AliITSOnlineSPDscanAnalyzer.h"
53ae21ce 39#include "AliITSOnlineCalibrationSPDhandler.h"
3ec5024c 40#include "AliLog.h"
41#include <iostream>
42#include <fstream>
43#include <TROOT.h>
44#include <TPluginManager.h>
53ae21ce 45#include <TObjArray.h>
46#include <TString.h>
3ec5024c 47
48int main(int argc, char **argv) {
49 if (argc<2) {
50 printf("Wrong number of arguments\n");
51 return -1;
52 }
53
6727e2db 54 // directory structure, hard coded
55 char *saveDirNoisy = "./calibResults/ScanNoisy"; // may NOT delete content
56 char *saveDirNoisyToFXS = "./calibResults/ScanNoisyToFXS"; // may delete content
57 char *saveDirDCSconfigToFXS= "./calibResults/ScanDCSconfigToFXS"; // may delete content
58 char *saveDirRef = "./calibResults/ScanReference"; // may delete content
4f88491f 59 char *saveDirIdsToFXS = "./calibResults/IdsToFXS"; // may delete content
6727e2db 60 char *configFilesDir = "./configFiles"; // may delete content
4f88491f 61
6727e2db 62 // make sure the directory structure is correct:
1fd93b67 63 system("mkdir ./calibResults >& /dev/null");
6727e2db 64 system("mkdir ./calibResults/ScanNoisy >& /dev/null");
65 system("mkdir ./calibResults/ScanNoisyToFXS >& /dev/null");
66 system("mkdir ./calibResults/ScanDCSconfigToFXS >& /dev/null");
67 system("mkdir ./calibResults/ScanReference >& /dev/null");
4f88491f 68 system("mkdir ./calibResults/IdsToFXS >& /dev/null");
6727e2db 69 system("mkdir ./configFiles >& /dev/null");
70 // prameters config files
71 TString paramsFileName = Form("%s/standal_params.txt",configFilesDir);
72 TString permNoisyFileName = Form("%s/perm_noisy.txt",configFilesDir);
53ae21ce 73
3ec5024c 74 // This line is needed in case of a stand-alone application w/o
75 // $ROOTSYS/etc/system.rootrc file
76 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
77 "*",
78 "TStreamerInfo",
79 "RIO",
80 "TStreamerInfo()");
81
82 // turn off annoying warning messages
83 new AliLog;
1fd93b67 84 AliLog::Instance()->SetGlobalDebugLevel(-20);
3ec5024c 85
86 // calib scan types
87 enum calib_types{MINTH,MEANTH,DAC,UNIMA,NOISE,DELAY};
88
89
4ee23d3d 90 // ********* STEP 0: Get configuration files from db (if there are any) , then read parameters*********
53ae21ce 91 UInt_t nrTuningParams = 0;
92 TObjArray paramNames; paramNames.SetOwner(kTRUE);
93 TObjArray paramVals; paramVals.SetOwner(kTRUE);
94
95 // tuning parameters:
96 Int_t status = 0;
97#ifndef SPD_DA_OFF
98 TString idp = "spd_standal_params";
6727e2db 99 status=daqDA_DB_getFile(idp.Data(),paramsFileName.Data());
53ae21ce 100 if (status) {
101 printf("Failed to get config file %s: status=%d. Using default tuning parameters.\n",idp.Data(),status);
102 }
103#endif
104 if (status==0) {
105 ifstream paramsFile;
6727e2db 106 paramsFile.open(paramsFileName.Data(), ifstream::in);
53ae21ce 107 if (paramsFile.fail()) {
6727e2db 108 printf("No config file (%s) present. Using default tuning parameters.\n",paramsFileName.Data());
53ae21ce 109 }
110 else {
111 while(1) {
112 Char_t paramN[50];
113 Char_t paramV[50];
114 paramsFile >> paramN;
115 if (paramsFile.eof()) break;
116 paramsFile >> paramV;
117 TString* paramNS = new TString(paramN);
118 TString* paramVS = new TString(paramV);
119 paramNames.AddAtAndExpand((TObject*)paramNS,nrTuningParams);
120 paramVals.AddAtAndExpand((TObject*)paramVS,nrTuningParams);
121 nrTuningParams++;
122 if (paramsFile.eof()) break;
123 }
124 paramsFile.close();
125 }
126 }
127 // for (UInt_t i=0; i<nrTuningParams; i++) {
128 // printf("Entry %d: N=%s , V=%s\n",i,((TString*)paramNames.At(i))->Data(),((TString*)paramVals.At(i))->Data());
129 // }
130
131 // perm noisy list:
132 Int_t permstatus = 0;
133#ifndef SPD_DA_OFF
134 TString idn = "spd_perm_noisy";
6727e2db 135 permstatus=daqDA_DB_getFile(idn.Data(),permNoisyFileName.Data());
53ae21ce 136 if (permstatus) {
6727e2db 137 printf("Failed to get config file %s: status=%d. No permanently noisy pixels will be added.\n",idn.Data(),permstatus);
53ae21ce 138 }
139#endif
140
141
142
143
144 // ********* STEP 1: Produce scan container files (Reference Data). ***********************************
3ec5024c 145 int startSeg = 1;
146
1fd93b67 147#ifndef SPD_DA_OFF
3ec5024c 148 int runNr = atoi(getenv("DATE_RUN_NUMBER"));
149#else
150 int runNr = atoi(argv[1]);
151 startSeg = 2;
152#endif
153
154 // container objects
155 AliITSOnlineSPDscan *scanObj[20];
156 Bool_t bScanInit[20];
157 for (UInt_t eqId=0; eqId<20; eqId++) {
158 scanObj[eqId]=NULL;
159 bScanInit[eqId]=kFALSE;
160 }
161 // header data variables
162 UInt_t routerNr[20];
163 Bool_t halfStaveScanned[20][6];
164 UInt_t type[20];
165 Bool_t dataFormat[20];
166 UInt_t triggers[20];
167 Bool_t chipPresent[20][6][10];
168 UInt_t dacStart[20];
169 UInt_t dacEnd[20];
170 UInt_t dacStep[20];
171 UInt_t dacId[20];
172 UInt_t rowStart[20];
173 UInt_t rowEnd[20];
174 UInt_t rowValue[20];
175 UInt_t dacValue[20];
176 UInt_t dacHigh[20][6];
177 UInt_t dacLow[20][6];
178 UInt_t TPAmp[20][6];
1fd93b67 179 Bool_t minTHchipPresent[20][10];
3ec5024c 180 // current scan step flag
181 Int_t currentStep[20];
182 for (UInt_t eqId=0; eqId<20; eqId++) currentStep[eqId] = 9999;
183
184 // loop over run segments
185 for (int segNr=startSeg; segNr<argc; segNr++) {
186
3ec5024c 187 int status;
188
189 /* define data source : this is argument 1 */
190 status=monitorSetDataSource( argv[segNr] );
191 if (status!=0) {
192 printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
193 return -1;
194 }
195 /* declare monitoring program */
196 status=monitorDeclareMp("ITS_SPD_CAL");
197 if (status!=0) {
198 printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
199 return -1;
200 }
201 /* define wait event timeout - 1s max */
202 monitorSetNowait();
203 monitorSetNoWaitNetworkTimeout(1000);
6727e2db 204
205
3ec5024c 206 UInt_t eventNr=0;
207
3ec5024c 208 /* main loop (infinite) */
209 for(;;) {
6727e2db 210
211 struct eventHeaderStruct *event;
212 eventTypeType eventT;
213
3ec5024c 214 /* check shutdown condition */
1fd93b67 215#ifndef SPD_DA_OFF
3ec5024c 216 if (daqDA_checkShutdown()) {break;}
1fd93b67 217#endif
3ec5024c 218 /* get next event (blocking call until timeout) */
219 status=monitorGetEventDynamic((void **)&event);
220 if (status==MON_ERR_EOF) {
221 printf ("End of File detected\n");
222 break; /* end of monitoring file has been reached */
223 }
224 if (status!=0) {
225 printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
226 break;
227 }
228 /* retry if got no event */
229 if (event==NULL) {
3ec5024c 230 continue;
231 }
232
233 eventT=event->eventType;
6727e2db 234 if (eventT == PHYSICS_EVENT) {
235
3ec5024c 236 eventNr++;
1fd93b67 237 // printf("eventNr %d\n",eventNr);
6727e2db 238
3ec5024c 239 AliRawReader *reader = new AliRawReaderDate((void*)event);
240 AliITSRawStreamSPD *str = new AliITSRawStreamSPD(reader);
6727e2db 241
3ec5024c 242 for (UInt_t eqId=0; eqId<20; eqId++) {
6727e2db 243
3ec5024c 244 reader->Reset();
245 reader->Select("ITSSPD",eqId,eqId);
6727e2db 246
3ec5024c 247 // Hit Event flags, specific for one event
248 Bool_t hitEventHSIncremented[6];
249 Bool_t hitEventChipIncremented[6][10];
250 for (UInt_t hs=0; hs<6; hs++) {
251 hitEventHSIncremented[hs] = kFALSE;
252 for (UInt_t chip=0; chip<10; chip++) {
253 hitEventChipIncremented[hs][chip] = kFALSE;
254 }
255 }
6727e2db 256
3f917533 257 if (str->ReadCalibHeader()>0) {
3ec5024c 258 // first check the type:
259 if (bScanInit[eqId] && type[eqId]!=str->GetHtype()) {
260 printf("Calib header problem. Type changed (%d -> %d)!\n",type[eqId],str->GetHtype());
261 }
262
263 // read calib values
1fd93b67 264 routerNr[eqId] = str->GetHrouterNr();
265 type[eqId] = str->GetHtype();
266 dataFormat[eqId] = str->GetHdataFormat();
267 triggers[eqId] = str->GetHtriggers();
268 dacStart[eqId] = str->GetHdacStart();
269 dacEnd[eqId] = str->GetHdacEnd();
270 dacStep[eqId] = str->GetHdacStep();
271 dacId[eqId] = str->GetHdacId();
272 rowStart[eqId] = str->GetHrowStart();
273 rowEnd[eqId] = str->GetHrowEnd();
274 rowValue[eqId] = str->GetHrowValue();
275 dacValue[eqId] = str->GetHdacValue();
3ec5024c 276 for (UInt_t hs=0; hs<6; hs++) {
1fd93b67 277 halfStaveScanned[eqId][hs] = str->GetHhalfStaveScanned(hs);
278 dacHigh[eqId][hs] = str->GetHdacHigh(hs);
279 dacLow[eqId][hs] = str->GetHdacLow(hs);
280 TPAmp[eqId][hs] = str->GetHTPAmp(hs);
3ec5024c 281 for (UInt_t chip=0; chip<10; chip++) {
1fd93b67 282 chipPresent[eqId][hs][chip] = str->GetHchipPresent(hs,chip);
3ec5024c 283 }
284 }
1fd93b67 285 for (UInt_t chip=0; chip<10; chip++) {
286 minTHchipPresent[eqId][chip] = str->GetHminTHchipPresent(chip);
287 }
3ec5024c 288
289 currentStep[eqId] = (dacValue[eqId]-dacStart[eqId])/dacStep[eqId];
290 if (type[eqId]==DELAY) {
291 currentStep[eqId]=currentStep[eqId]*2;
292 dacValue[eqId]=dacValue[eqId]*2;
293 if (dacHigh[eqId][0]==128) { // misc_ctrl value
294 currentStep[eqId]=currentStep[eqId]+1;
295 dacValue[eqId]=dacValue[eqId]+1;
296 }
297 }
298
299 // router nr check:
300 if (routerNr[eqId]!=eqId) {
301 printf("Router nr problem? Router nr %d != EqID %d\n",routerNr[eqId],eqId);
302 }
303
304 if (!bScanInit[eqId]) {
305 // initialize container object
53ae21ce 306 TString fileName = Form("%s/SPDcal_run_%d_eq_%d.root",saveDirRef,runNr,eqId);
3ec5024c 307 switch (type[eqId]) {
308 case NOISE:
309 case UNIMA:
53ae21ce 310 scanObj[eqId] = new AliITSOnlineSPDscanSingle(fileName.Data());
3ec5024c 311 ((AliITSOnlineSPDscanSingle*)scanObj[eqId])->ClearThis();
312 bScanInit[eqId]=kTRUE;
313 break;
314 case MINTH:
315 case DAC:
316 case DELAY:
53ae21ce 317 scanObj[eqId] = new AliITSOnlineSPDscanMultiple(fileName.Data());
3ec5024c 318 scanObj[eqId]->ClearThis();
319 bScanInit[eqId]=kTRUE;
320 break;
321 case MEANTH:
53ae21ce 322 scanObj[eqId] = new AliITSOnlineSPDscanMeanTh(fileName.Data());
3ec5024c 323 scanObj[eqId]->ClearThis();
324 bScanInit[eqId]=kTRUE;
325 break;
326 default:
327 printf("Unknown scan type: %d.\n",type[eqId]);
328 }
329 // some multiple scan data
330 if (type[eqId]==MINTH || type[eqId]==MEANTH || type[eqId]==DAC || type[eqId]==DELAY) {
331 ((AliITSOnlineSPDscanMultiple*)scanObj[eqId])->SetDacId(dacId[eqId]);
332 }
333 // some common data
334 scanObj[eqId]->SetRunNr((UInt_t)runNr);
335 scanObj[eqId]->SetRouterNr(routerNr[eqId]);
336 for (UInt_t hs=0; hs<6; hs++) {
337 scanObj[eqId]->SetHalfStaveScanned(hs,halfStaveScanned[eqId][hs]);
338 }
339 scanObj[eqId]->SetType(type[eqId]);
340 scanObj[eqId]->SetDataFormat(dataFormat[eqId]);
341 for (Int_t hs=0; hs<6; hs++) {
342
343 // remove later when the chip present is set correctly !!!!!!!!!!!!!!!!!!!!!!!!!!!
344 Bool_t halfStavePresent = str->GetHalfStavePresent(hs);
345 // remove later when the chip present is set correctly !!!!!!!!!!!!!!!!!!!!!!!!!!!
346
347 for (UInt_t chip=0; chip<10; chip++) {
348 scanObj[eqId]->SetChipPresent(hs,chip,chipPresent[eqId][hs][chip]);
349
350 // remove later when the chip present is set correctly !!!!!!!!!!!!!!!!!!!!!!!!!!!
351 if (halfStavePresent) scanObj[eqId]->SetChipPresent(hs,chip,kTRUE);
352 // remove later when the chip present is set correctly !!!!!!!!!!!!!!!!!!!!!!!!!!!
353
354 }
355 }
356 scanObj[eqId]->SetRowStart(rowStart[eqId]);
357 scanObj[eqId]->SetRowEnd(rowEnd[eqId]);
358 scanObj[eqId]->SetDacStart(dacStart[eqId]);
359 scanObj[eqId]->SetDacEnd(dacEnd[eqId]);
360 scanObj[eqId]->SetDacStep(dacStep[eqId]);
361 }
362
363 if (type[eqId]==MINTH) {
364 scanObj[eqId]->SetTriggers(currentStep[eqId],triggers[eqId]);
365 }
366 if (type[eqId]==UNIMA || type[eqId]==NOISE) {
367 if (currentStep[eqId]==9999) printf("SPDcalibratorStep1 (eq %d): single step\n",eqId);
368 currentStep[eqId]=0;
369 }
370 if (type[eqId]==MINTH || type[eqId]==MEANTH || type[eqId]==DAC || type[eqId]==DELAY) {
371 ((AliITSOnlineSPDscanMultiple*)scanObj[eqId])->SetDacValue(currentStep[eqId],dacValue[eqId]);
372 if (type[eqId]==DELAY) {
373 printf("SPDcalibratorStep1 (eq %d): DAC %d/%d , step %d\n",eqId,dacValue[eqId]/2,dacHigh[eqId][0],currentStep[eqId]);
374 }
375 else {
376 printf("SPDcalibratorStep1 (eq %d): DAC %d , step %d\n",eqId,dacValue[eqId],currentStep[eqId]);
377 }
378 }
379 if (type[eqId]==MEANTH) {
380 for (Int_t hs=0; hs<6; hs++) {
381 ((AliITSOnlineSPDscanMeanTh*)scanObj[eqId])->SetDacLow(currentStep[eqId],hs,dacLow[eqId][hs]);
382 ((AliITSOnlineSPDscanMeanTh*)scanObj[eqId])->SetDacHigh(currentStep[eqId],hs,dacHigh[eqId][hs]);
383 ((AliITSOnlineSPDscanMeanTh*)scanObj[eqId])->SetTPAmp(currentStep[eqId],hs,TPAmp[eqId][hs]);
384 }
385 }
386
387
388 }
53ae21ce 389
3ec5024c 390 if (bScanInit[eqId]) {
391 while (str->Next()) {
392 UInt_t hs = str->GetHalfStaveNr();
393 UInt_t chip = str->GetChipAddr();
53ae21ce 394 //***remove last condition when minthpresent put correctly in calib header?
1fd93b67 395#ifndef SPD_DA_OFF
53ae21ce 396 if (type[eqId]!=MINTH || minTHchipPresent[eqId][chip] || runNr<=416900) {
1fd93b67 397#else
53ae21ce 398 if (type[eqId]!=MINTH || minTHchipPresent[eqId][chip] || runNr<=416900) {
1fd93b67 399#endif
6727e2db 400 //*************************************************************************
3ec5024c 401 scanObj[eqId]->IncrementHits(currentStep[eqId],hs,chip,str->GetChipCol(),str->GetChipRow());
6727e2db 402
3ec5024c 403 if (!hitEventHSIncremented[hs]) {
404 scanObj[eqId]->IncrementHitEventsTot(currentStep[eqId],hs);
405 hitEventHSIncremented[hs]=kTRUE;
406 }
6727e2db 407
3ec5024c 408 if (!hitEventChipIncremented[hs][chip]) {
409 scanObj[eqId]->IncrementHitEvents(currentStep[eqId],hs,chip);
410 hitEventChipIncremented[hs][chip]=kTRUE;
411 }
412 }
413
414 }
415
416 if (type[eqId]!=MINTH) { // for minth, triggers are set from header info
417 scanObj[eqId]->IncrementTriggers(currentStep[eqId]);
418 }
419
420 }
421
422 }
423
424 delete str;
425 delete reader;
426
427 }
428
429 /* free resources */
430 free(event);
431
432 }
433
1fd93b67 434#ifndef SPD_DA_OFF
3ec5024c 435 daqDA_progressReport((unsigned int)( ((Float_t)(segNr-startSeg+1))/(argc-startSeg)*50 ));
436#else
437 printf("progress: %d\n",(unsigned int)( ((Float_t)(segNr-startSeg+1))/(argc-startSeg)*50 ));
438#endif
439
440 }
441
6727e2db 442
443 // clean up scan objects (which also saves them)
3ec5024c 444 for (UInt_t eqId=0; eqId<20; eqId++) {
445 if (scanObj[eqId]!=NULL) delete scanObj[eqId];
446 }
447
448
449
450
451
53ae21ce 452 // ********* STEP 2: Analyze scan container files. ************************************************
453
454 // clear noisyToFXS and DCSconfigToFXS dirs:
455 TString command = Form("cd %s; rm -f *",saveDirNoisyToFXS);
456 system(command.Data());
457 TString command2 = Form("cd %s; rm -f *",saveDirDCSconfigToFXS);
458 system(command2.Data());
459 UInt_t nrNoisyFilesProduced=0;
460 UInt_t nrDCSconfigFilesProduced=0;
3ec5024c 461
6727e2db 462 AliITSOnlineCalibrationSPDhandler* handler = new AliITSOnlineCalibrationSPDhandler();
463 AliITSOnlineSPDscanAnalyzer *analyzer = NULL;
464 AliITSOnlineCalibrationSPDhandler* handlerPermNoisy = NULL;
465 // fill permanent noisy list to add later...
466 if (permstatus==0) {
467 handlerPermNoisy = new AliITSOnlineCalibrationSPDhandler();
468 UInt_t permNoisy = handlerPermNoisy->ReadNoisyFromText(permNoisyFileName.Data(),240); // 240 = read for all modules
469 if (permNoisy>0) {
470 printf("%d noisy pixels read from permanent list.\n",permNoisy);
471 }
472 }
473
4f88491f 474 Bool_t reset_made = kFALSE;
3ec5024c 475
476 // *** *** *** start loop over equipments (eq_id)
477 for (int eqId=0; eqId<20; eqId++) {
478
53ae21ce 479 // create analyzer for this eq
480 TString fileName = Form("%s/SPDcal_run_%d_eq_%d.root",saveDirRef,runNr,eqId);
6727e2db 481 analyzer = new AliITSOnlineSPDscanAnalyzer(fileName.Data(),handler);
53ae21ce 482
483 // configure analyzer with tuning parameters etc:
484 for (UInt_t i=0; i<nrTuningParams; i++) {
485 analyzer->SetParam(((TString*)paramNames.At(i))->Data(),((TString*)paramVals.At(i))->Data());
486 }
3ec5024c 487
4ee23d3d 488 UInt_t type = analyzer->GetType();
489 UInt_t dacId = analyzer->GetDacId();
490 UInt_t routerNr = analyzer->GetRouterNr();
3ec5024c 491 if (type!=99) {
492 if (type==DAC) {
6727e2db 493 printf("SPD scan calibrator Step2: eqId %d, type %d, dacId %d\n",eqId,type,dacId);
3ec5024c 494 }
6727e2db 495 else printf("SPD scan calibrator Step2: eqId %d type %d\n",eqId,type);
3ec5024c 496 }
4ee23d3d 497
498
3ec5024c 499
53ae21ce 500 // algorithms for the different types of scans:
3ec5024c 501
502 if (type==UNIMA) {
53ae21ce 503
3ec5024c 504 }
53ae21ce 505
3ec5024c 506 else if (type==NOISE) {
4f88491f 507 // read previous noisy list (clear if overwriting)
6727e2db 508 handler->SetFileLocation(saveDirNoisy);
4f88491f 509 if (analyzer->IsOverWriteSet() && !reset_made) {
510 handler->ResetNoisy();
511 handler->WriteToFilesAlways();
512 reset_made=kTRUE;
513 }
514 else {
515 handler->ReadFromFiles();
516 }
6727e2db 517 if (analyzer->ProcessNoisyPixels(/*saveDirNoisy*/)) {
518 if (permstatus==0) {
519 handler->AddNoisyFrom(handlerPermNoisy);
520 }
53ae21ce 521 // init dcs config text file
522 TString dcsConfigFileName = Form("%s/dcsConfig_run_%d_eq_%d.txt",saveDirDCSconfigToFXS,runNr,eqId);
523 ofstream dcsfile;
524 dcsfile.open(dcsConfigFileName.Data());
525 dcsfile << "[SPD SCAN]\n";
526 dcsfile << "RunNumber=" << runNr << "\n";
527 dcsfile << "Type=" << type << "\n";
4ee23d3d 528 dcsfile << "Router=" << routerNr << "\n";
6727e2db 529 dcsfile << "ActualDetConfiguration=" << "0,-1,-1\n"; // dummy values for now
53ae21ce 530 dcsfile << "[NOISY]\n";
531 nrDCSconfigFilesProduced++;
4f88491f 532
533 for (UInt_t hs=0; hs<6; hs++) {
534 for (UInt_t chip=0; chip<10; chip++) {
535 if (analyzer->IsChipPresent(hs,chip) || analyzer->IsOverWriteSet()) {
536 dcsfile << "-" << eqId << "," << hs << "," << chip << "\n";
537 UInt_t nrNoisy = handler->GetNrNoisyC(eqId,hs,chip);
538 for (UInt_t ind=0; ind<nrNoisy; ind++) {
539 UInt_t col = handler->GetNoisyColAtC(eqId,hs,chip,ind);
540 UInt_t row = handler->GetNoisyRowAtC(eqId,hs,chip,ind);
541 dcsfile << col << "," << row << "\n";
53ae21ce 542 }
53ae21ce 543 }
3ec5024c 544 }
545 }
4f88491f 546 handler->SetFileLocation(saveDirNoisy);
547 handler->WriteNoisyToFile(eqId);
548 handler->SetFileLocation(saveDirNoisyToFXS);
549 handler->WriteNoisyToFile(eqId);
550 nrNoisyFilesProduced++;
551
53ae21ce 552 dcsfile.close();
3ec5024c 553 }
554 }
53ae21ce 555
3ec5024c 556 else if (type==MINTH || (type==DAC && dacId==39)) {
53ae21ce 557 // init dcs config text file
558 TString dcsConfigFileName = Form("%s/dcsConfig_run_%d_eq_%d.txt",saveDirDCSconfigToFXS,runNr,eqId);
559 ofstream dcsfile;
560 dcsfile.open(dcsConfigFileName.Data());
561 dcsfile << "[SPD SCAN]\n";
562 dcsfile << "RunNumber=" << runNr << "\n";
563 dcsfile << "Type=" << type << "\n";
4ee23d3d 564 dcsfile << "Router=" << routerNr << "\n";
4f88491f 565 dcsfile << "ActualDetConfiguration=" << "0,-1,-1\n"; // dummy values for now
53ae21ce 566 dcsfile << "[DACvalues]\n";
567 nrDCSconfigFilesProduced++;
3ec5024c 568 for (UInt_t hs=0; hs<6; hs++) {
569 for (UInt_t chipNr=0; chipNr<10; chipNr++) {
53ae21ce 570 Int_t minTh = -1;
571 if (analyzer->GetOnlineScan()->GetChipPresent(hs,chipNr)) {
572 minTh = analyzer->GetMinTh(hs,chipNr);
573 if (minTh!=-1) {
574 dcsfile << "39," << eqId << "," << hs << "," << chipNr << "=" << minTh << "\n";
575 }
576 else {
577 printf("MinTh failed for Eq %d , HS %d , Chip %d\n",eqId,hs,chipNr);
578 }
579 }
3ec5024c 580 }
3ec5024c 581 }
53ae21ce 582 dcsfile.close();
3ec5024c 583 }
53ae21ce 584
3ec5024c 585 else if (type==DELAY) {
53ae21ce 586 // init dcs config text file
587 TString dcsConfigFileName = Form("%s/dcsConfig_run_%d_eq_%d.txt",saveDirDCSconfigToFXS,runNr,eqId);
588 ofstream dcsfile;
589 dcsfile.open(dcsConfigFileName.Data());
590 dcsfile << "[SPD SCAN]\n";
591 dcsfile << "RunNumber=" << runNr << "\n";
592 dcsfile << "Type=" << type << "\n";
4ee23d3d 593 dcsfile << "Router=" << routerNr << "\n";
53ae21ce 594 dcsfile << "ActualDetCoonfiguration=" << "0,-1,-1\n"; // dummy values for now
595 dcsfile << "[DACvalues]\n";
596 nrDCSconfigFilesProduced++;
3ec5024c 597 for (UInt_t hs=0; hs<6; hs++) {
598 for (UInt_t chipNr=0; chipNr<10; chipNr++) {
53ae21ce 599 Int_t clockCycle = -1;
600 Int_t delayCtrl = -1;
601 Int_t miscCtrl = -1;
602 if (analyzer->GetOnlineScan()->GetChipPresent(hs,chipNr)) {
603 clockCycle = analyzer->GetDelay(hs,chipNr);
604 delayCtrl = clockCycle/2;
605 miscCtrl = 192;
606 if (clockCycle!=-1) {
607 if (clockCycle%2==1) miscCtrl = 128;
608 dcsfile << "42," << eqId << "," << hs << "," << chipNr << "=" << delayCtrl << "\n";
609 dcsfile << "43," << eqId << "," << hs << "," << chipNr << "=" << miscCtrl << "\n";
610 }
611 else {
612 printf("Delay failed for Eq %d , HS %d , Chip %d\n",eqId,hs,chipNr);
613 }
614 }
3ec5024c 615 }
3ec5024c 616 }
53ae21ce 617 dcsfile.close();
3ec5024c 618 }
619
3ec5024c 620 delete analyzer;
621
1fd93b67 622#ifndef SPD_DA_OFF
3ec5024c 623 daqDA_progressReport((unsigned int)(50+(eqId+1)*2.5));
624#else
625 printf("progress: %d\n",(unsigned int)(50+(eqId+1)*2.5));
626#endif
627
628 }
629 // *** *** *** end loop over equipments (eq_id)
630
4f88491f 631
6727e2db 632 delete handler;
633 if (handlerPermNoisy!=NULL) {
634 delete handlerPermNoisy;
635 }
636
637
4f88491f 638 TString idsFXSFileName = Form("%s/FXSids_run_%d.txt",saveDirIdsToFXS,runNr);
639 ofstream idsFXSfile;
640 idsFXSfile.open(idsFXSFileName.Data());
641
3ec5024c 642
6727e2db 643
644 // send noisy data to FXS
1fd93b67 645 if (nrNoisyFilesProduced>0) {
646 // send a tared file of all new noisy maps
6727e2db 647 TString command = Form("cd %s; tar -cf noisy_scan.tar *",saveDirNoisyToFXS);
648 // printf("\n\n%s\n\n",command.Data());
53ae21ce 649 system(command.Data());
6727e2db 650 TString fileName = Form("%s/noisy_scan.tar",saveDirNoisyToFXS);
651 TString id = "SPD_scan_noisy";
4f88491f 652#ifndef SPD_DA_OFF
53ae21ce 653 Int_t status = daqDA_FES_storeFile(fileName.Data(),id.Data());
1fd93b67 654 if (status!=0) {
53ae21ce 655 printf("Failed to export file %s , status %d\n",fileName.Data(),status);
1fd93b67 656 return -1;
657 }
53ae21ce 658#endif
4f88491f 659 idsFXSfile << Form("%s\n",id.Data());
3ec5024c 660 }
53ae21ce 661
6727e2db 662 // send dcs config files to FXS
53ae21ce 663 if (nrDCSconfigFilesProduced>0) {
664 // send a tared file of all the dcsConfig text files
665 TString command = Form("cd %s; tar -cf dcsConfig.tar *",saveDirDCSconfigToFXS);
666 // printf("\n\n%s\n\n",command.Data());
667 system(command.Data());
53ae21ce 668 TString fileName = Form("%s/dcsConfig.tar",saveDirDCSconfigToFXS);
669 TString id = "SPD_dcsConfig";
4f88491f 670#ifndef SPD_DA_OFF
53ae21ce 671 Int_t status = daqDA_FES_storeFile(fileName.Data(),id.Data());
672 if (status!=0) {
673 printf("Failed to export file %s , status %d\n",fileName.Data(),status);
674 return -1;
675 }
3ec5024c 676#endif
4f88491f 677 //idsFXSfile << Form("%s\n",id.Data()); // do NOT write this id (this is not for preprocessor)
53ae21ce 678 }
679
6727e2db 680 // send reference data to FXS
6727e2db 681 for (UInt_t eqId=0; eqId<20; eqId++) {
682 if (bScanInit[eqId]) {
683 TString fileName = Form("%s/SPDcal_run_%d_eq_%d.root",saveDirRef,runNr,eqId);
684 TString id = Form("SPD_ref_scan_%d",eqId);
4f88491f 685#ifndef SPD_DA_OFF
6727e2db 686 Int_t status = daqDA_FES_storeFile(fileName.Data(),id.Data());
687 if (status!=0) {
688 printf("Failed to export file %s , status %d\n",fileName.Data(),status);
689 return -1;
690 }
4f88491f 691#endif
692 idsFXSfile << Form("%s\n",id.Data());
6727e2db 693 }
694 }
6727e2db 695
696
4f88491f 697 idsFXSfile.close();
6727e2db 698
3ec5024c 699
700
701 return 0;
702}