]>
Commit | Line | Data |
---|---|---|
803e6db5 | 1 | /* |
2 | Contact: henrik.tydesjo@cern.ch | |
3 | Link: tydes.home.cern.ch/tydes/doc/CalibrationOverview/CalibrationAlgorithms/ | |
4 | Run Type: PHYSICS | |
5 | DA Type: MON | |
6 | Number of events needed: Depending on muliplicity per event | |
7 | Input Files: spd_physics_params , ./calibResults/DeadReferenceTmp/* , ./calibResults/DeadToFXS/* | |
8 | Output Files: ./calibResults/NoisyReference/* , ./calibResults/DeadReference/* , ./calibResults/NoisyToFXS/* , ./calibResults/DeadReferenceTmp/*,./calibResults/DeadToFXS/* | |
9 | Trigger types used: PHYSICS | |
10 | */ | |
11 | ||
12 | //////////////////////////////////////////////////////////////////////////////// | |
13 | // This program can be compiled in two modes. // | |
14 | // // | |
15 | // 1. Online. With the DAQ DA framework. This is the default operating mode. // | |
16 | // // | |
17 | // 2. Offline. Without the DAQ DA framework. Define the SPD_DA_OFF // | |
18 | // environment var. Call this program with the name of the executable // | |
19 | // followed by the runNr and the data files to process. // | |
20 | // // | |
21 | //////////////////////////////////////////////////////////////////////////////// | |
22 | ||
23 | #ifndef SPD_DA_OFF | |
24 | extern "C" { | |
25 | #include "daqDA.h" | |
26 | } | |
27 | #endif | |
28 | #include "event.h" | |
29 | #include "monitor.h" | |
30 | #include "AliRawReaderDate.h" | |
31 | #include "AliITSRawStreamSPD.h" | |
32 | #include "AliITSOnlineSPDphys.h" | |
33 | #include "AliITSOnlineSPDphysAnalyzer.h" | |
34 | #include "AliITSOnlineCalibrationSPDhandler.h" | |
35 | #include "AliLog.h" | |
36 | #include <iostream> | |
37 | #include <fstream> | |
38 | #include <TROOT.h> | |
39 | #include <TPluginManager.h> | |
40 | #include <TObjArray.h> | |
41 | #include <TString.h> | |
42 | ||
43 | int main(int argc, char **argv) { | |
44 | if (argc<2) { | |
45 | printf("Wrong number of arguments\n"); | |
46 | return -1; | |
47 | } | |
48 | ||
49 | // directory structure, hard coded | |
50 | char *saveDirDead = "./calibResults/Dead"; // may NOT delete content | |
51 | char *saveDirDeadToFXS = "./calibResults/DeadToFXS"; // may delete content | |
52 | char *saveDirDeadRef = "./calibResults/DeadReference"; // may delete content | |
53 | char *saveDirDeadRefTmp = "./calibResults/DeadReferenceTmp"; // may NOT delete content | |
54 | char *saveDirNoisyToFXS = "./calibResults/NoisyToFXS"; // may delete content | |
55 | char *saveDirNoisyRef = "./calibResults/NoisyReference"; // may delete content | |
56 | char *saveDirIdsToFXS = "./calibResults/IdsToFXS"; // may delete content | |
57 | char *configFilesDir = "./configFiles"; // may delete content | |
58 | // make sure the directory structure is put up correctly: | |
59 | system("mkdir ./calibResults >& /dev/null"); | |
60 | system("mkdir ./calibResults/Dead >& /dev/null"); | |
61 | system("mkdir ./calibResults/DeadToFXS >& /dev/null"); | |
62 | system("mkdir ./calibResults/DeadReference >& /dev/null"); | |
63 | system("mkdir ./calibResults/DeadReferenceTmp >& /dev/null"); | |
64 | system("mkdir ./calibResults/NoisyToFXS >& /dev/null"); | |
65 | system("mkdir ./calibResults/NoisyReference >& /dev/null"); | |
66 | system("mkdir ./calibResults/IdsToFXS >& /dev/null"); | |
67 | system("mkdir ./configFiles >& /dev/null"); | |
68 | // parameters config file | |
69 | TString paramsFileName = Form("%s/physics_params.txt",configFilesDir); | |
70 | ||
71 | // This line is needed in case of a stand-alone application w/o | |
72 | // $ROOTSYS/etc/system.rootrc file | |
73 | gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", | |
74 | "*", | |
75 | "TStreamerInfo", | |
76 | "RIO", | |
77 | "TStreamerInfo()"); | |
78 | ||
79 | // turn off annoying warning messages | |
80 | new AliLog; | |
81 | AliLog::Instance()->SetGlobalDebugLevel(-20); | |
82 | ||
83 | ||
84 | // ********* STEP 0: Get configuration files from db (if there are any) , then read parameters********* | |
85 | UInt_t nrTuningParams = 0; | |
86 | TObjArray paramNames; paramNames.SetOwner(kTRUE); | |
87 | TObjArray paramVals; paramVals.SetOwner(kTRUE); | |
88 | ||
89 | // tuning parameters: | |
90 | Int_t par_status = 0; | |
91 | #ifndef SPD_DA_OFF | |
92 | TString idp = "spd_physics_params"; | |
93 | par_status=daqDA_DB_getFile(idp.Data(),paramsFileName.Data()); | |
94 | if (par_status) { | |
95 | printf("Failed to get config file %s: status=%d. Using default tuning parameters.\n",idp.Data(),par_status); | |
96 | } | |
97 | #endif | |
98 | if (par_status==0) { | |
99 | ifstream paramsFile; | |
100 | paramsFile.open(paramsFileName.Data(), ifstream::in); | |
101 | if (paramsFile.fail()) { | |
102 | printf("No config file (%s) present. Using default tuning parameters.\n",paramsFileName.Data()); | |
103 | } | |
104 | else { | |
105 | while(1) { | |
106 | Char_t paramN[50]; | |
107 | Char_t paramV[50]; | |
108 | paramsFile >> paramN; | |
109 | if (paramsFile.eof()) break; | |
110 | paramsFile >> paramV; | |
111 | TString* paramNS = new TString(paramN); | |
112 | TString* paramVS = new TString(paramV); | |
113 | paramNames.AddAtAndExpand((TObject*)paramNS,nrTuningParams); | |
114 | paramVals.AddAtAndExpand((TObject*)paramVS,nrTuningParams); | |
115 | nrTuningParams++; | |
116 | if (paramsFile.eof()) break; | |
117 | } | |
118 | paramsFile.close(); | |
119 | } | |
120 | } | |
121 | // for (UInt_t i=0; i<nrTuningParams; i++) { | |
122 | // printf("Entry %d: N=%s , V=%s\n",i,((TString*)paramNames.At(i))->Data(),((TString*)paramVals.At(i))->Data()); | |
123 | // } | |
124 | ||
125 | ||
126 | ||
127 | ||
128 | ||
129 | // ********* STEP 1: Produce phys container files (Reference Data). *********************************** | |
130 | ||
131 | #ifndef SPD_DA_OFF | |
132 | if (getenv("DATE_RUN_NUMBER")==0) { | |
133 | printf("DATE_RUN_NUMBER not properly set.\n"); | |
134 | return -1; | |
135 | } | |
136 | int runNr = atoi(getenv("DATE_RUN_NUMBER")); | |
137 | #else | |
138 | int runNr = atoi(argv[1]); | |
139 | int startSeg = 2; | |
140 | #endif | |
141 | ||
142 | ||
143 | // container objects | |
144 | AliITSOnlineSPDphys *physObj[20]; | |
145 | Bool_t bPhysInit[20]; | |
146 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
147 | physObj[eqId]=NULL; | |
148 | bPhysInit[eqId]=kFALSE; | |
149 | } | |
150 | ||
151 | ||
152 | // loop over run segments in case of offline mode | |
153 | #ifdef SPD_DA_OFF | |
154 | for (int segNr=startSeg; segNr<argc; segNr++) { | |
155 | #endif | |
156 | ||
157 | int status; | |
158 | ||
159 | /* define data source : */ | |
160 | #ifndef SPD_DA_OFF | |
161 | status=monitorSetDataSource( argv[1] ); // should be "^SPD" in order to get full detector online | |
162 | #else | |
163 | status=monitorSetDataSource( argv[segNr] ); | |
164 | #endif | |
165 | if (status!=0) { | |
166 | printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); | |
167 | return -1; | |
168 | } | |
169 | /* declare monitoring program */ | |
170 | status=monitorDeclareMp("ITS_SPD_PHYS"); | |
171 | if (status!=0) { | |
172 | printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status)); | |
173 | return -1; | |
174 | } | |
175 | /* define wait event timeout - 1s max */ | |
176 | monitorSetNowait(); | |
177 | monitorSetNoWaitNetworkTimeout(1000); | |
178 | ||
179 | ||
180 | UInt_t eventNr=0; | |
181 | ||
182 | /* main loop (infinite) */ | |
183 | for(;;) { | |
184 | ||
185 | struct eventHeaderStruct *event; | |
186 | eventTypeType eventT; | |
187 | ||
188 | /* check shutdown condition */ | |
189 | #ifndef SPD_DA_OFF | |
190 | if (daqDA_checkShutdown()) {break;} | |
191 | #endif | |
192 | ||
193 | /* get next event (blocking call until timeout) */ | |
194 | status=monitorGetEventDynamic((void **)&event); | |
195 | if (status==MON_ERR_EOF) { | |
196 | printf ("End of File detected\n"); | |
197 | break; /* end of monitoring file has been reached */ | |
198 | } | |
199 | ||
200 | if (status!=0) { | |
201 | printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); | |
202 | break; | |
203 | } | |
204 | ||
205 | /* retry if got no event */ | |
206 | if (event==NULL) { | |
207 | continue; | |
208 | } | |
209 | ||
210 | eventT=event->eventType; | |
211 | if (eventT == PHYSICS_EVENT){ | |
212 | ||
213 | eventNr++; | |
214 | // printf("eventNr %d\n",eventNr); | |
215 | ||
216 | AliRawReader *reader = new AliRawReaderDate((void*)event); | |
217 | AliITSRawStreamSPD *str = new AliITSRawStreamSPD(reader); | |
218 | ||
219 | // for (UInt_t eqId=0; eqId<20; eqId++) { | |
220 | // reader->Reset(); | |
221 | // reader->Select("ITSSPD",eqId,eqId); | |
222 | ||
223 | while (str->Next()) { | |
224 | ||
225 | Int_t eqId = reader->GetDDLID(); | |
226 | if (eqId>=0 && eqId<20) { | |
227 | if (!bPhysInit[eqId]) { | |
228 | TString fileName = Form("%s/SPDphys_run_%d_eq_%d.root",saveDirNoisyRef,runNr,eqId); | |
229 | physObj[eqId] = new AliITSOnlineSPDphys(fileName.Data()); | |
230 | physObj[eqId]->AddRunNr(runNr); | |
231 | physObj[eqId]->SetEqNr(eqId); | |
232 | bPhysInit[eqId]=kTRUE; | |
233 | } | |
234 | ||
235 | UInt_t hs = str->GetHalfStaveNr(); | |
236 | UInt_t chip = str->GetChipAddr(); | |
237 | physObj[eqId]->IncrementHits(hs,chip,str->GetChipCol(),str->GetChipRow()); | |
238 | ||
239 | } | |
240 | } | |
241 | ||
242 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
243 | if (bPhysInit[eqId]) { | |
244 | physObj[eqId]->IncrementNrEvents(); | |
245 | } | |
246 | } | |
247 | ||
248 | // } | |
249 | ||
250 | delete str; | |
251 | delete reader; | |
252 | ||
253 | } | |
254 | ||
255 | /* free resources */ | |
256 | free(event); | |
257 | ||
258 | } | |
259 | ||
260 | ||
261 | #ifdef SPD_DA_OFF | |
262 | printf("progress: %d\n",(unsigned int)( ((Float_t)(segNr-startSeg+1))/(argc-startSeg)*50 )); | |
263 | } | |
264 | #endif | |
265 | ||
266 | // clean up phys objects (also saves them) | |
267 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
268 | if (physObj[eqId]!=NULL) delete physObj[eqId]; | |
269 | } | |
270 | ||
271 | ||
272 | ||
273 | ||
274 | ||
275 | // ********* STEP 2: Analyze phys container files. ************************************************ | |
276 | ||
277 | // clear noisyToFXS dir: | |
278 | TString command; | |
279 | command = Form("cd %s; rm -f *",saveDirNoisyToFXS); | |
280 | system(command.Data()); | |
281 | // clear deadToFXS dir: | |
282 | command = Form("cd %s; rm -f *",saveDirDeadToFXS); | |
283 | system(command.Data()); | |
284 | ||
285 | ||
286 | // create calibration handler and read dead from previous calibrations | |
287 | AliITSOnlineCalibrationSPDhandler* handler = new AliITSOnlineCalibrationSPDhandler(); | |
288 | handler->SetFileLocation(saveDirDead); | |
289 | handler->ReadDeadFromFiles(); | |
290 | ||
291 | ||
292 | UInt_t firstRunNrDead = runNr; | |
293 | ||
294 | ||
295 | UInt_t nrEnoughStatNoisy = 0; | |
296 | UInt_t nrEqActiveNoisy = 0; | |
297 | Bool_t eqActiveNoisy[20]; | |
298 | ||
299 | // *** *** *** start loop over equipments (eq_id) | |
300 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
301 | eqActiveNoisy[eqId] = kFALSE; | |
302 | ||
303 | // create analyzer for this eq | |
304 | TString fileName = Form("%s/SPDphys_run_%d_eq_%d.root",saveDirNoisyRef,runNr,eqId); | |
305 | AliITSOnlineSPDphysAnalyzer *noisyAnalyzer = new AliITSOnlineSPDphysAnalyzer(fileName.Data(),handler); | |
306 | ||
307 | // check data in container | |
308 | if (noisyAnalyzer->GetEqNr() != eqId) { | |
309 | if (noisyAnalyzer->GetEqNr() != 999) { | |
310 | printf("Error: Mismatching EqId in Container data and filename (%d!=%d). Skipping.\n", | |
311 | noisyAnalyzer->GetEqNr(),eqId); | |
312 | } | |
313 | delete noisyAnalyzer; | |
314 | continue; | |
315 | } | |
316 | ||
317 | nrEqActiveNoisy++; | |
318 | eqActiveNoisy[eqId] = kTRUE; | |
319 | ||
320 | // configure analyzer with tuning parameters etc: | |
321 | for (UInt_t i=0; i<nrTuningParams; i++) { | |
322 | noisyAnalyzer->SetParam(((TString*)paramNames.At(i))->Data(),((TString*)paramVals.At(i))->Data()); | |
323 | } | |
324 | ||
325 | printf("SPD phys STEP 2: Noisy search for eq %d\n",eqId); | |
326 | ||
327 | // search for noisy pixels: | |
328 | nrEnoughStatNoisy += noisyAnalyzer->ProcessNoisyPixels(); | |
329 | ||
330 | // copy this phys obj to temporary dead reference dir to process after noisy search | |
331 | TString fileNameDead = Form("%s/SPDphys_dead_run_0_0_eq_%d.root",saveDirDeadRefTmp,eqId); | |
332 | AliITSOnlineSPDphys* physObj = new AliITSOnlineSPDphys(fileNameDead.Data()); | |
333 | physObj->AddPhys(noisyAnalyzer->GetOnlinePhys()); | |
334 | if (physObj->GetNrRuns()>0) { | |
335 | UInt_t firstRunNr = physObj->GetRunNr(0); | |
336 | if (firstRunNrDead>firstRunNr) { | |
337 | firstRunNrDead=firstRunNr; | |
338 | } | |
339 | } | |
340 | // remove noisy pixels from dead hitmap | |
341 | for (UInt_t hs=0; hs<6; hs++) { | |
342 | for (UInt_t chip=0; chip<10; chip++) { | |
343 | for (UInt_t ind=0; ind<handler->GetNrNoisyC(eqId,hs,chip); ind++) { | |
344 | UInt_t col = handler->GetNoisyColAtC(eqId,hs,chip,ind); | |
345 | UInt_t row = handler->GetNoisyRowAtC(eqId,hs,chip,ind); | |
346 | physObj->AddHits(hs,chip,col,row,-noisyAnalyzer->GetOnlinePhys()->GetHits(hs,chip,col,row)); | |
347 | } | |
348 | } | |
349 | } | |
350 | ||
351 | delete physObj; | |
352 | delete noisyAnalyzer; | |
353 | ||
354 | #ifndef SPD_DA_OFF | |
355 | daqDA_progressReport((unsigned int)((eqId+1)*2.5)); | |
356 | #else | |
357 | printf("progress: %d\n",(unsigned int)(50+(eqId+1)*1.25)); | |
358 | #endif | |
359 | } | |
360 | // *** *** *** end loop over equipments (eq_id) | |
361 | ||
362 | printf("Noisy search finished. %d noisy pixels found. %d chips (%d) had enough statistics.\n", | |
363 | handler->GetNrNoisy(),nrEnoughStatNoisy,nrEqActiveNoisy*60); | |
364 | handler->SetFileLocation(saveDirNoisyToFXS); | |
365 | handler->WriteNoisyToFiles(); | |
366 | ||
367 | ||
368 | ||
369 | ||
370 | ||
371 | ||
372 | ||
373 | ||
374 | UInt_t nrEnoughStatChips = 0; | |
375 | UInt_t nrDeadChips = 0; | |
376 | UInt_t nrInefficientChips = 0; | |
377 | UInt_t nrEqActiveDead = 0; | |
378 | Bool_t eqActiveDead[20]; | |
379 | ||
380 | // *** *** *** start loop over equipments (eq_id) | |
381 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
382 | eqActiveDead[eqId] = kFALSE; | |
383 | ||
384 | // setup analyzer for dead search | |
385 | TString fileNameDead = Form("%s/SPDphys_dead_run_0_0_eq_%d.root",saveDirDeadRefTmp,eqId); | |
386 | AliITSOnlineSPDphys* physObj = new AliITSOnlineSPDphys(fileNameDead.Data()); | |
387 | AliITSOnlineSPDphysAnalyzer* deadAnalyzer = new AliITSOnlineSPDphysAnalyzer(physObj,handler); | |
388 | // check data in container | |
389 | if (deadAnalyzer->GetEqNr() != eqId) { | |
390 | if (deadAnalyzer->GetEqNr() != 999) { | |
391 | printf("Error: Mismatching EqId in Dead Container data and filename (%d!=%d). Skipping.\n", | |
392 | deadAnalyzer->GetEqNr(),eqId); | |
393 | } | |
394 | delete deadAnalyzer; | |
395 | continue; | |
396 | } | |
397 | ||
398 | nrEqActiveDead++; | |
399 | eqActiveDead[eqId] = kTRUE; | |
400 | ||
401 | // configure analyzer with tuning parameters etc: | |
402 | for (UInt_t i=0; i<nrTuningParams; i++) { | |
403 | deadAnalyzer->SetParam(((TString*)paramNames.At(i))->Data(),((TString*)paramVals.At(i))->Data()); | |
404 | } | |
405 | ||
406 | printf("SPD phys STEP 2: Dead search for eq %d\n",eqId); | |
407 | ||
408 | // search for dead pixels: | |
409 | nrEnoughStatChips += deadAnalyzer->ProcessDeadPixels(); | |
410 | nrDeadChips += deadAnalyzer->GetNrDeadChips(); | |
411 | nrInefficientChips += deadAnalyzer->GetNrInefficientChips(); | |
412 | ||
413 | delete deadAnalyzer; | |
414 | ||
415 | #ifndef SPD_DA_OFF | |
416 | daqDA_progressReport((unsigned int)(50+(eqId+1)*2.5)); | |
417 | #else | |
418 | printf("progress: %d\n",(unsigned int)(75+(eqId+1)*1.25)); | |
419 | #endif | |
420 | } | |
421 | // *** *** *** end loop over equipments (eq_id) | |
422 | ||
423 | ||
424 | printf("Dead search finished. %d dead pixels in total.\n%d chips (%d) had enough statistics. %d chips were dead. %d chips were inefficient.\n",handler->GetNrDead(),nrEnoughStatChips,nrEqActiveDead*60,nrDeadChips,nrInefficientChips); | |
425 | handler->SetFileLocation(saveDirDead); | |
426 | handler->WriteDeadToFilesAlways(); | |
427 | ||
428 | handler->SetFileLocation(saveDirDeadToFXS); | |
429 | handler->WriteDeadToFilesAlways(); | |
430 | // *** old code (used if not all dead data should be uploaded) | |
431 | // UInt_t nrDeadFilesToTar = 0; | |
432 | // handler->SetFileLocation(saveDirDeadToFXS); | |
433 | // for (UInt_t eqId=0; eqId<20; eqId++) { | |
434 | // if (eqActiveDead[eqId]) { | |
435 | // handler->WriteDeadToFile(eqId); | |
436 | // nrDeadFilesToTar++; | |
437 | // } | |
438 | // } | |
439 | //*** | |
440 | ||
441 | printf("Opening id list file\n"); | |
442 | TString idsFXSFileName = Form("%s/FXSids_run_%d.txt",saveDirIdsToFXS,runNr); | |
443 | ofstream idsFXSfile; | |
444 | idsFXSfile.open(idsFXSFileName.Data()); | |
445 | ||
446 | ||
447 | // if there is no chip in category "needsMoreStat" | |
448 | if (nrEnoughStatChips+nrDeadChips+nrInefficientChips == nrEqActiveDead*60) { | |
449 | // calibration is complete | |
450 | printf("Dead calibration is complete.\n"); | |
451 | ||
452 | ||
453 | ||
454 | printf("Preparing dead reference data\n"); | |
455 | // send reference data for dead pixels to FXS | |
456 | TString tarFiles = ""; | |
457 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
458 | if (eqActiveDead[eqId]) { | |
459 | printf("Preparing dead pixels for eq %d\n",eqId); | |
460 | // move file to ref dir | |
461 | TString fileName = Form("%s/SPDphys_dead_run_0_0_eq_%d.root",saveDirDeadRefTmp,eqId); | |
462 | TString newFileName = Form("%s/SPDphys_dead_run_%d_%d_eq_%d.root",saveDirDeadRef,firstRunNrDead,runNr,eqId); | |
463 | TString command = Form("mv -f %s %s",fileName.Data(),newFileName.Data()); | |
464 | system(command.Data()); | |
465 | ||
466 | tarFiles.Append(Form("SPDphys_dead_run_%d_%d_eq_%d.root ",firstRunNrDead,runNr,eqId)); | |
467 | } | |
468 | } | |
469 | TString send_command = Form("cd %s; tar -cf ref_phys_dead.tar %s",saveDirDeadRef,tarFiles.Data()); | |
470 | system(send_command.Data()); | |
471 | TString fileName = Form("%s/ref_phys_dead.tar",saveDirDeadRef); | |
472 | TString id = "SPD_ref_phys_dead"; | |
473 | #ifndef SPD_DA_OFF | |
474 | status = daqDA_FES_storeFile(fileName.Data(),id.Data()); | |
475 | if (status!=0) { | |
476 | printf("Failed to export file %s , status %d\n",fileName.Data(),status); | |
477 | return -1; | |
478 | } | |
479 | #endif | |
480 | idsFXSfile << Form("%s\n",id.Data()); | |
481 | // for (UInt_t eqId=0; eqId<20; eqId++) { // OLD CODE NOT TARED | |
482 | // if (eqActiveDead[eqId]) { | |
483 | // TString fileName = Form("%s/SPDphys_dead_run_0_0_eq_%d.root",saveDirDeadRefTmp,eqId); | |
484 | // // move file to ref dir | |
485 | // TString newFileName = Form("%s/SPDphys_dead_run_%d_%d_eq_%d.root",saveDirDeadRef,firstRunNrDead,runNr,eqId); | |
486 | // TString command = Form("mv -f %s %s",fileName.Data(),newFileName.Data()); | |
487 | // system(command.Data()); | |
488 | // // send ref data to FXS | |
489 | // TString id = Form("SPD_ref_phys_dead_%d",eqId); | |
490 | //#ifndef SPD_DA_OFF | |
491 | // Int_t status = daqDA_FES_storeFile(newFileName.Data(),id.Data()); | |
492 | // if (status!=0) { | |
493 | // printf("Failed to export file %s , status %d\n",newFileName.Data(),status); | |
494 | // return -1; | |
495 | // } | |
496 | //#endif | |
497 | // idsFXSfile << Form("%s\n",id.Data()); | |
498 | // } | |
499 | // else { | |
500 | // TString command = Form("rm -f %s/SPDphys_dead_run_0_0_eq_%d.root",saveDirDeadRefTmp,eqId); | |
501 | // system(command.Data()); | |
502 | // } | |
503 | // } | |
504 | } | |
505 | ||
506 | ||
507 | printf("Preparing noisy reference data\n"); | |
508 | // send reference data for this run to FXS | |
509 | TString tarFiles = ""; | |
510 | for (UInt_t eqId=0; eqId<20; eqId++) { | |
511 | if (eqActiveNoisy[eqId]) { | |
512 | tarFiles.Append(Form("SPDphys_run_%d_eq_%d.root ",runNr,eqId)); | |
513 | } | |
514 | } | |
515 | TString send_command = Form("cd %s; tar -cf ref_phys.tar %s",saveDirNoisyRef,tarFiles.Data()); | |
516 | system(send_command.Data()); | |
517 | TString fileName = Form("%s/ref_phys.tar",saveDirNoisyRef); | |
518 | TString id = "SPD_ref_phys"; | |
519 | #ifndef SPD_DA_OFF | |
520 | status = daqDA_FES_storeFile(fileName.Data(),id.Data()); | |
521 | if (status!=0) { | |
522 | printf("Failed to export file %s , status %d\n",fileName.Data(),status); | |
523 | return -1; | |
524 | } | |
525 | #endif | |
526 | idsFXSfile << Form("%s\n",id.Data()); | |
527 | ||
528 | ||
529 | // // send reference data for this run to FXS - OLD CODE NOT TARED | |
530 | // for (UInt_t eqId=0; eqId<20; eqId++) { | |
531 | // if (eqActiveNoisy[eqId]) { | |
532 | // TString fileName = Form("%s/SPDphys_run_%d_eq_%d.root",saveDirNoisyRef,runNr,eqId); | |
533 | // TString id = Form("SPD_ref_phys_%d",eqId); | |
534 | //#ifndef SPD_DA_OFF | |
535 | // Int_t status = daqDA_FES_storeFile(fileName.Data(),id.Data()); | |
536 | // if (status!=0) { | |
537 | // printf("Failed to export file %s , status %d\n",fileName.Data(),status); | |
538 | // return -1; | |
539 | // } | |
540 | // idsFXSfile << Form("%s\n",id.Data()); | |
541 | // } | |
542 | // } | |
543 | //#endif | |
544 | ||
545 | ||
546 | ||
547 | ||
548 | printf("Preparing dead files\n"); | |
549 | // send dead pixels to FXS | |
550 | // if (nrDeadFilesToTar>0) { //*** old code (used if not all dead data should be uploaded) | |
551 | // send a tared file of all the dead files | |
552 | send_command = Form("cd %s; tar -cf dead_phys.tar *",saveDirDeadToFXS); | |
553 | // printf("\n\n%s\n\n",command.Data()); | |
554 | system(send_command.Data()); | |
555 | fileName = Form("%s/dead_phys.tar",saveDirDeadToFXS); | |
556 | id = "SPD_phys_dead"; | |
557 | #ifndef SPD_DA_OFF | |
558 | Int_t send_status = daqDA_FES_storeFile(fileName.Data(),id.Data()); | |
559 | if (send_status!=0) { | |
560 | printf("Failed to export file %s , status %d\n",fileName.Data(),send_status); | |
561 | return -1; | |
562 | } | |
563 | #endif | |
564 | idsFXSfile << Form("%s\n",id.Data()); | |
565 | // } //*** old code (used if not all dead data should be uploaded) | |
566 | ||
567 | ||
568 | // send noisy pixels to FXS | |
569 | if (handler->GetNrNoisy()>0) { // there must be at least one file created | |
570 | printf("Preparing noisy files\n"); | |
571 | // send a tared file of all the noisy files | |
572 | TString command = Form("cd %s; tar -cf noisy_phys.tar *",saveDirNoisyToFXS); | |
573 | // printf("\n\n%s\n\n",command.Data()); | |
574 | system(command.Data()); | |
575 | TString fileName = Form("%s/noisy_phys.tar",saveDirNoisyToFXS); | |
576 | TString id = "SPD_phys_noisy"; | |
577 | #ifndef SPD_DA_OFF | |
578 | status = daqDA_FES_storeFile(fileName.Data(),id.Data()); | |
579 | if (status!=0) { | |
580 | printf("Failed to export file %s , status %d\n",fileName.Data(),status); | |
581 | return -1; | |
582 | } | |
583 | #endif | |
584 | idsFXSfile << Form("%s\n",id.Data()); | |
585 | } | |
586 | ||
587 | ||
588 | // send ids file to FXS | |
589 | idsFXSfile.close(); | |
590 | id = "SPD_id_list"; | |
591 | #ifndef SPD_DA_OFF | |
592 | status = daqDA_FES_storeFile(idsFXSFileName.Data(),id.Data()); | |
593 | if (status!=0) { | |
594 | printf("Failed to export file %s , status %d\n",idsFXSFileName.Data(),status); | |
595 | return -1; | |
596 | } | |
597 | #endif | |
598 | ||
599 | ||
600 | ||
601 | ||
602 | delete handler; | |
603 | ||
604 | ||
605 | return 0; | |
606 | } |