1 |
// |
// |
2 |
// Check for possible errors and shows raw distribution of variables coming from the DSP (no calibrations needed) - Emiliano Mocchiutti |
// Check for possible errors and shows raw distribution of variables coming from the DSP (no calibrations needed) - Emiliano Mocchiutti |
3 |
// |
// |
4 |
// FCaloQLOOK.c version 1.08 (2006-05-29) |
// FCaloQLOOK.c version 1.09 (2006-05-29) |
5 |
// |
// |
6 |
// The only input needed is the path to the directory created by YODA for the data file you want to analyze. |
// The only input needed is the path to the directory created by YODA for the data file you want to analyze. |
7 |
// |
// |
8 |
// Changelog: |
// Changelog: |
9 |
// |
// |
10 |
// 1.07 - 1.08 (2006-05-29): Fixed bug in output filename when input is not in the form DW_YYMMDD_NNN. |
// 1.07 - 1.09 (2006-05-29): Fixed bug in output filename when input is not in the form DW_YYMMDD_NNN. Bug in latchup alarm not coming from the status word, fixed. |
11 |
// |
// |
12 |
// 1.06 - 1.07 (2006-05-23): Don't print warning in case of latchup not recognized from the status word, fixed. Added "expert" figures from 21 to 28 which appears only in the |
// 1.06 - 1.07 (2006-05-23): Don't print warning in case of latchup not recognized from the status word, fixed. Added "expert" figures from 21 to 28 which appears only in the |
13 |
// case of latchup alarm. |
// case of latchup alarm. |
281 |
bool isRAW = 0; |
bool isRAW = 0; |
282 |
Int_t alldexy=0; |
Int_t alldexy=0; |
283 |
Int_t alldexy2=0; |
Int_t alldexy2=0; |
284 |
Int_t stri=0; |
Int_t planebases=0; |
285 |
|
// Int_t stri=0; |
286 |
Int_t fcheck = 0; |
Int_t fcheck = 0; |
287 |
Int_t cestw=0; |
Int_t cestw=0; |
288 |
Int_t cmask = 127 ; |
Int_t cmask = 127 ; |
433 |
bdone = 0; |
bdone = 0; |
434 |
bl++; |
bl++; |
435 |
allbase = ce->base[l][ii][bl]; |
allbase = ce->base[l][ii][bl]; |
436 |
|
// alldexy=0; |
437 |
|
// alldexy2=0; |
438 |
|
// stri=0; |
439 |
|
// for (Int_t e = 0; e < 16 ; e++ ){ |
440 |
|
// stri = e + 16 * bl; |
441 |
|
// alldexy += (int)ce->dexyc[l][ii][stri]; |
442 |
|
// alldexy2 += (int)ce->dexy[l][ii][stri]; |
443 |
|
// }; |
444 |
|
}; |
445 |
|
// |
446 |
|
if ( kk == 0 ){ |
447 |
|
planebases = 0; |
448 |
alldexy=0; |
alldexy=0; |
449 |
alldexy2=0; |
alldexy2=0; |
450 |
stri=0; |
for (Int_t e = 0; e < 96 ; e++ ){ |
451 |
for (Int_t e = 0; e < 16 ; e++ ){ |
if ( e < 6 ) planebases += (int)ce->base[l][ii][e]; |
452 |
stri = e + 16 * bl; |
alldexy += (int)ce->dexyc[l][ii][e]; |
453 |
alldexy += (int)ce->dexyc[l][ii][stri]; |
alldexy2 += (int)ce->dexy[l][ii][e]; |
|
alldexy2 += (int)ce->dexy[l][ii][stri]; |
|
454 |
}; |
}; |
455 |
}; |
}; |
456 |
|
// |
457 |
if ( !isRAW ) { |
if ( !isRAW ) { |
458 |
// |
// |
459 |
if ( !pdone ){ |
if ( !pdone ){ |
460 |
if ( (ce->base[l][ii][bl]>32000 || ce->base[l][ii][bl] == 0 ) && ( alldexy > 512000 || alldexy == 0) && ce->perror[se] == 0 ) { |
// if ( (ce->base[l][ii][bl]>32000 || ce->base[l][ii][bl] == 0 ) && ( alldexy > 512000 || alldexy == 0) && ce->perror[se] == 0 ) { |
461 |
|
if ( (planebases>192000 || planebases == 0 ) && ( alldexy > 3072000 || alldexy == 0) && ce->perror[se] == 0 ) { |
462 |
pdone = 1; |
pdone = 1; |
463 |
pshit[se][pl]++ ; |
pshit[se][pl]++ ; |
464 |
if ( (ce->stwerr[se] & (1 << 4)) == 0 ) { |
if ( (ce->stwerr[se] & (1 << 4)) == 0 ) { |
526 |
// |
// |
527 |
} else { |
} else { |
528 |
if ( !pdone ){ |
if ( !pdone ){ |
529 |
if ( (alldexy2>512000 || alldexy2 == 0) && ce->perror[se] == 0 ) { |
if ( ( alldexy2 > 3072000 || alldexy2 == 0) && ce->perror[se] == 0 ) { |
530 |
|
// if ( (alldexy2>512000 || alldexy2 == 0) && ce->perror[se] == 0 ) { |
531 |
pdone = 1; |
pdone = 1; |
532 |
pshit[se][pl]++ ; |
pshit[se][pl]++ ; |
533 |
if ( (ce->stwerr[se] & (1 << 4)) == 0 ){ |
if ( (ce->stwerr[se] & (1 << 4)) == 0 ){ |