14 |
#include <TSQLRow.h> |
#include <TSQLRow.h> |
15 |
#include <TTree.h> |
#include <TTree.h> |
16 |
#include <TGraph.h> |
#include <TGraph.h> |
|
#include <TDatime.h> |
|
17 |
#include <TTimeStamp.h> |
#include <TTimeStamp.h> |
18 |
#include <TF1.h> |
#include <TF1.h> |
19 |
// |
// |
157 |
}; |
}; |
158 |
|
|
159 |
/** |
/** |
160 |
|
* Set the autoboot flag |
161 |
|
* |
162 |
|
*/ |
163 |
|
void PamelaDBOperations::SetAutoBoot(Bool_t dbg){ |
164 |
|
AUTOBOOT = dbg; |
165 |
|
}; |
166 |
|
|
167 |
|
/** |
168 |
* Set the nofrag flag |
* Set the nofrag flag |
169 |
* |
* |
170 |
*/ |
*/ |
303 |
UInt_t nevent = 0; |
UInt_t nevent = 0; |
304 |
UInt_t pktlast = 0; |
UInt_t pktlast = 0; |
305 |
UInt_t obtlast = 0; |
UInt_t obtlast = 0; |
306 |
UInt_t t_pktlast = 0; |
Long64_t t_pktlast = 0LL; |
307 |
UInt_t t_obtlast = 0; |
// UInt_t t_obtlast = 0; |
308 |
UInt_t upperpkt2 = 0; |
Long64_t t_obtlast = 0LL; |
309 |
ULong64_t upperobt2 = 0; |
Long64_t upperpkt2 = 0LL; |
310 |
|
Long64_t upperobt2 = 0LL; |
311 |
UInt_t zomp = 0; |
UInt_t zomp = 0; |
312 |
UInt_t jump = 50000; // was 5000 |
UInt_t jump = 50000; // was 5000 |
313 |
EventCounter *code=0; |
EventCounter *code=0; |
314 |
// |
// |
315 |
UInt_t deltapkt = 5000; |
Long64_t deltapkt = 5000LL; |
316 |
ULong64_t deltaobt = 50000; |
Long64_t deltaobt = 50000LL; |
317 |
// |
// |
318 |
// pcksList packetsNames; |
// pcksList packetsNames; |
319 |
// pcksList::iterator Iter; |
// pcksList::iterator Iter; |
351 |
upperobt = OBT(obtlast); |
upperobt = OBT(obtlast); |
352 |
upperentry = nevent-1; |
upperentry = nevent-1; |
353 |
// |
// |
354 |
if ( IsDebug() ) printf(" First entries are: OBT %llu pkt_num %i \n",obtfirst,pktfirst); |
if ( IsDebug() ) printf(" First entries are: OBT %i pkt_num %i \n",obtfirst,pktfirst); |
355 |
// |
// |
356 |
if ( IsDebug() ) printf(" Last entries are: OBT %llu pkt_num %i entry %i\n",upperobt,upperpkt,upperentry); |
if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld entry %i\n",upperobt,upperpkt,upperentry); |
357 |
// |
// |
358 |
if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(1); |
if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(1); |
359 |
// |
// |
384 |
upperpkt2 = PKT(ph->GetCounter()); |
upperpkt2 = PKT(ph->GetCounter()); |
385 |
upperobt2 = OBT(ph->GetOrbitalTime()); |
upperobt2 = OBT(ph->GetOrbitalTime()); |
386 |
// |
// |
387 |
if ( (t_pktlast < upperpkt && t_obtlast > upperobt) || (t_pktlast < upperpkt2 && t_obtlast > upperobt2) ) throw -13; |
if ( (t_pktlast < upperpkt && t_obtlast > upperobt) || (t_pktlast < upperpkt2 && t_obtlast > upperobt2) ){ |
388 |
|
if ( IsDebug() ) printf(" .-. upperpkt2 %lld upperobt2 %lld \n",upperpkt2,upperobt2); |
389 |
|
if ( IsDebug() ) printf(" .-. upperpkt %lld t_pktlast %lld upperobt %lld t_obtlast %lld \n",upperpkt,t_pktlast,upperobt,t_obtlast); |
390 |
|
if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %u upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i); |
391 |
|
throw -13; |
392 |
|
}; |
393 |
// |
// |
394 |
if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){ |
if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){ |
395 |
zomp = i + jump + 1; |
zomp = i + jump + 1; |
396 |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
397 |
if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %i pktlast %i upperobt %llu obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i); |
if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i); |
398 |
break; |
break; |
399 |
}; |
}; |
400 |
// |
// |
429 |
// |
// |
430 |
rhev = rh->GetEntries(); |
rhev = rh->GetEntries(); |
431 |
rtev = rt->GetEntries(); |
rtev = rt->GetEntries(); |
432 |
UInt_t sobtt = 0; |
Long64_t sobtt = 0LL; |
433 |
UInt_t sobth = 0; |
Long64_t sobth = 0LL; |
434 |
UInt_t spktt = 0; |
Long64_t spktt = 0LL; |
435 |
UInt_t spkth = 0; |
Long64_t spkth = 0LL; |
436 |
UInt_t pktt = 0; |
Long64_t pktt = 0LL; |
437 |
ULong64_t obtt = 0; |
Long64_t obtt = 0LL; |
438 |
UInt_t pkth = 0; |
Long64_t pkth = 0LL; |
439 |
ULong64_t obth = 0; |
Long64_t obth = 0LL; |
440 |
// |
// |
441 |
T->GetEntry(upperentry); |
T->GetEntry(upperentry); |
442 |
code = eh->GetCounter(); |
code = eh->GetCounter(); |
456 |
obth = OBT(phh->GetOrbitalTime()); |
obth = OBT(phh->GetOrbitalTime()); |
457 |
}; |
}; |
458 |
// |
// |
459 |
if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
460 |
if ( pkth > upperpkt && obth > upperobt ){ |
if ( pkth > upperpkt && obth > upperobt ){ |
461 |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt); |
462 |
upperpkt = pkth; |
upperpkt = pkth; |
463 |
upperobt = obth; |
upperobt = obth; |
464 |
rhev = lasthead+1; |
rhev = lasthead+1; |
465 |
} else { |
} else { |
466 |
rhev = lasthead; |
rhev = lasthead; |
467 |
}; |
}; |
468 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
469 |
// |
// |
470 |
if ( IsDebug() ) printf(" rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rtev beforev %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
471 |
if ( pktt > upperpkt && obtt > upperobt ){ |
if ( pktt > upperpkt && obtt > upperobt ){ |
472 |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt); |
473 |
upperpkt = pktt; |
upperpkt = pktt; |
474 |
upperobt = obtt; |
upperobt = obtt; |
475 |
rtev = lasttrail+1; |
rtev = lasttrail+1; |
476 |
} else { |
} else { |
477 |
rtev = lasttrail; |
rtev = lasttrail; |
478 |
}; |
}; |
479 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
480 |
// goto kikko; |
// goto kikko; |
481 |
// |
// |
482 |
// |
// |
483 |
// Check if runtrailer/runheader are within lower limits |
// Check if runtrailer/runheader are within lower limits |
484 |
// |
// |
485 |
// |
// |
486 |
pkth = 0; |
pkth = 0LL; |
487 |
obth = 0; |
obth = 0LL; |
488 |
spkth = 0; |
spkth = 0LL; |
489 |
sobth = 0; |
sobth = 0LL; |
490 |
for (Int_t k=0; k<rhev; k++){ |
for (Int_t k=0; k<rhev; k++){ |
491 |
if ( k > 0 ){ |
if ( k > 0 ){ |
492 |
spkth = pkth; |
spkth = pkth; |
534 |
upperentry = evbefh-1; |
upperentry = evbefh-1; |
535 |
}; |
}; |
536 |
}; |
}; |
537 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
538 |
goto kikko0; |
goto kikko0; |
539 |
}; |
}; |
540 |
}; |
}; |
542 |
// |
// |
543 |
// |
// |
544 |
// |
// |
545 |
pktt = 0; |
pktt = 0LL; |
546 |
obtt = 0; |
obtt = 0LL; |
547 |
spktt = 0; |
spktt = 0LL; |
548 |
sobtt = 0; |
sobtt = 0LL; |
549 |
for (Int_t k=0; k<rtev; k++){ |
for (Int_t k=0; k<rtev; k++){ |
550 |
if ( k > 0 ){ |
if ( k > 0 ){ |
551 |
spktt = pktt; |
spktt = pktt; |
565 |
rt->GetEntry(rtev); |
rt->GetEntry(rtev); |
566 |
pktt = spktt; |
pktt = spktt; |
567 |
obtt = sobtt; |
obtt = sobtt; |
568 |
if ( IsDebug() ) printf(" lasttrail %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
569 |
// |
// |
570 |
UInt_t evbeft = 0; |
UInt_t evbeft = 0; |
571 |
code = eht->GetCounter(); |
code = eht->GetCounter(); |
594 |
upperentry = evbeft-1; |
upperentry = evbeft-1; |
595 |
}; |
}; |
596 |
}; |
}; |
597 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
598 |
goto kikko; |
goto kikko; |
599 |
// break; |
// break; |
600 |
// |
// |
622 |
obth = OBT(phh->GetOrbitalTime()); |
obth = OBT(phh->GetOrbitalTime()); |
623 |
}; |
}; |
624 |
// |
// |
625 |
if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
626 |
if ( pkth > upperpkt && obth > upperobt ){ |
if ( pkth > upperpkt && obth > upperobt ){ |
627 |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt); |
628 |
upperpkt = pkth; |
upperpkt = pkth; |
629 |
upperobt = obth; |
upperobt = obth; |
630 |
rhev = lasthead+1; |
rhev = lasthead+1; |
631 |
} else { |
} else { |
632 |
rhev = lasthead; |
rhev = lasthead; |
633 |
}; |
}; |
634 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
635 |
// |
// |
636 |
if ( IsDebug() ) printf(" rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rtev beforev %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
637 |
if ( pktt > upperpkt && obtt > upperobt ){ |
if ( pktt > upperpkt && obtt > upperobt ){ |
638 |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt); |
639 |
upperpkt = pktt; |
upperpkt = pktt; |
640 |
upperobt = obtt; |
upperobt = obtt; |
641 |
rtev = lasttrail+1; |
rtev = lasttrail+1; |
642 |
} else { |
} else { |
643 |
rtev = lasttrail; |
rtev = lasttrail; |
644 |
}; |
}; |
645 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
646 |
// |
// |
647 |
if ( IsDebug() ) printf(" Upper limits are: OBT %llu pkt_num %i upper entry %i \n",upperobt,upperpkt,upperentry); |
if ( IsDebug() ) printf(" Upper limits are: OBT %lld pkt_num %lld upper entry %i \n",upperobt,upperpkt,upperentry); |
648 |
// |
// |
649 |
return(0); |
return(0); |
650 |
} |
} |
782 |
// |
// |
783 |
// set DB timezone to UTC |
// set DB timezone to UTC |
784 |
// |
// |
785 |
stringstream oss; |
stringstream oss; |
786 |
// |
// |
787 |
oss.str(""); |
oss.str(""); |
788 |
oss << "SET time_zone='+0:00';"; |
oss << "SET time_zone='+0:00';"; |
789 |
conn->Query(oss.str().c_str()); |
TSQLResult *result = 0; |
790 |
|
result = conn->Query(oss.str().c_str()); |
791 |
|
if ( !result ) throw -10; |
792 |
// |
// |
793 |
}; |
}; |
794 |
|
|
795 |
/** |
/** |
796 |
* Return the correct packet number if we went back to zero |
* Return the correct packet number if we went back to zero |
797 |
*/ |
*/ |
798 |
UInt_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
799 |
// |
// |
800 |
// if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u (pktfirst - (UInt_t)(16777212/2)) is %u \n",pkt_num,pktfirst,(pktfirst - (UInt_t)(16777212/2))); |
if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u (UInt_t)(16777214/2)) is %u \n",pkt_num,pktfirst,(UInt_t)(16777214/2)); |
801 |
// |
// |
802 |
if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2) ) return((pkt_num+16777215)); |
if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2) ){ |
803 |
|
if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL); |
804 |
|
return((Long64_t)pkt_num+16777215LL); |
805 |
|
}; |
806 |
// |
// |
807 |
if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){ |
if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){ |
808 |
if ( (pkt_num-16777215) < 0 ){ |
if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL); |
809 |
return((16777215-pkt_num)); |
return((Long64_t)pkt_num-16777215LL); |
|
} else { |
|
|
return((pkt_num-16777215)); |
|
|
}; |
|
810 |
}; |
}; |
811 |
// |
// |
812 |
return(pkt_num); |
if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)pkt_num); |
813 |
|
return((Long64_t)pkt_num); |
814 |
// |
// |
815 |
}; |
}; |
816 |
|
|
817 |
/** |
/** |
818 |
* Return the correct On Board Time if we went back to zero |
* Return the correct On Board Time if we went back to zero |
819 |
*/ |
*/ |
820 |
ULong64_t PamelaDBOperations::OBT(UInt_t obt){ |
Long64_t PamelaDBOperations::OBT(UInt_t obt){ |
821 |
// |
// |
822 |
if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((ULong64_t)(obt+numeric_limits<UInt_t>::max())); |
if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((Long64_t)(obt+numeric_limits<UInt_t>::max())); |
823 |
// |
// |
824 |
if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
825 |
if ( (obt-numeric_limits<UInt_t>::max()) < 0 ){ |
return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
|
return((ULong64_t)(numeric_limits<UInt_t>::max()-obt)); |
|
|
} else { |
|
|
return((ULong64_t)(obt-numeric_limits<UInt_t>::max())); |
|
|
}; |
|
826 |
}; |
}; |
827 |
// |
// |
828 |
return((ULong64_t)obt); |
return((Long64_t)obt); |
829 |
}; |
}; |
830 |
|
|
831 |
/** |
/** |
966 |
if ( this->GetID_RAW() == 0 ) throw -11; |
if ( this->GetID_RAW() == 0 ) throw -11; |
967 |
// |
// |
968 |
oss.str(""); |
oss.str(""); |
969 |
oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='" |
oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='" |
970 |
<< this->GetRawFile().Data() << "';"; |
<< this->GetRawFile().Data() << "';"; |
971 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
972 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
975 |
// |
// |
976 |
if ( !row ){ |
if ( !row ){ |
977 |
oss.str(""); |
oss.str(""); |
978 |
oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< " |
oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< " |
979 |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
980 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
981 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
984 |
if ( !row ) throw -10; |
if ( !row ) throw -10; |
985 |
}; |
}; |
986 |
// |
// |
987 |
TDatime ti = TDatime(row->GetField(0)); |
TTimeStamp tu = TTimeStamp((UInt_t)atoi(row->GetField(0)),(UInt_t)atoi(row->GetField(1)),(UInt_t)atoi(row->GetField(2)),(UInt_t)atoi(row->GetField(3)),(UInt_t)atoi(row->GetField(4)),(UInt_t)atoi(row->GetField(5)),0,true,0); |
|
TTimeStamp tu = TTimeStamp((UInt_t)ti.GetYear(),(UInt_t)ti.GetMonth(),(UInt_t)ti.GetDay(),(UInt_t)ti.GetHour(),(UInt_t)ti.GetMinute(),(UInt_t)ti.GetSecond(),0,true,0); |
|
988 |
t0 = (UInt_t)tu.GetSec(); |
t0 = (UInt_t)tu.GetSec(); |
989 |
|
if ( IsDebug() ) printf(" t0 is %u ti is %s %s %s %s %s %s\n",t0,row->GetField(0),row->GetField(1),row->GetField(2),row->GetField(3),row->GetField(4),row->GetField(5)); |
990 |
// |
// |
|
// t0 = (UInt_t)TDatime(row->GetField(0)).Convert(); |
|
991 |
/* |
/* |
992 |
* Verify that the TIMESYNC have been not already processed |
* Verify that the TIMESYNC have been not already processed |
993 |
*/ |
*/ |
1002 |
if (result == NULL) throw -10; |
if (result == NULL) throw -10; |
1003 |
row = result->Next(); |
row = result->Next(); |
1004 |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
1005 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
1006 |
toffset = (UInt_t)atoll(row->GetField(2)) - (UInt_t)(this->OBT((UInt_t)atoll(row->GetField(1)))/1000) + t0; |
toffset = (UInt_t)atoll(row->GetField(2)) - (UInt_t)(this->OBT((UInt_t)atoll(row->GetField(1)))/1000) + t0; |
1007 |
|
// |
1008 |
|
tsync = (UInt_t)atoll(row->GetField(2)); |
1009 |
|
obt0 = (UInt_t)atoll(row->GetField(1)); |
1010 |
|
// |
1011 |
return(1); |
return(1); |
1012 |
}; |
}; |
1013 |
// |
// |
1201 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
1202 |
if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str()); |
1203 |
// |
// |
1204 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
1205 |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
1206 |
// |
// |
1207 |
|
tsync = TSYNC; |
1208 |
|
obt0 = OBT; |
1209 |
|
// |
1210 |
delete result; |
delete result; |
1211 |
return(signal); |
return(signal); |
1212 |
} |
} |
1306 |
VarDumpEvent *vde = 0; |
VarDumpEvent *vde = 0; |
1307 |
VarDumpRecord *vdr = 0; |
VarDumpRecord *vdr = 0; |
1308 |
// |
// |
1309 |
|
Bool_t found = false; |
1310 |
trDumpEv->SetBranchAddress("VarDump", &vde); |
trDumpEv->SetBranchAddress("VarDump", &vde); |
1311 |
if ( trDumpEv->GetEntries() > 0 ){ |
if ( trDumpEv->GetEntries() > 0 ){ |
1312 |
Bool_t found = false; |
found = false; |
1313 |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
1314 |
trDumpEv->GetEntry(i); |
trDumpEv->GetEntry(i); |
1315 |
vde->Records->GetEntries(); |
// vde->Records->GetEntries(); |
1316 |
if ( vde->Records->GetEntries()>0 ){ |
if ( vde->Records->GetEntries()>5 ){ |
1317 |
found = true; |
found = true; |
1318 |
goto fill; |
goto fill; |
1319 |
}; |
}; |
1326 |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
1327 |
// |
// |
1328 |
} else { |
} else { |
1329 |
if ( !this->GetBOOTnumber() ) return(4); |
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4); |
1330 |
|
}; |
1331 |
|
} else { |
1332 |
|
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2); |
1333 |
|
}; |
1334 |
|
// |
1335 |
|
UInt_t bn = 0; |
1336 |
|
Bool_t afound = false; |
1337 |
|
if ( !found && this->AutoBoot()){ |
1338 |
|
afound = true; |
1339 |
|
// |
1340 |
|
// Search for other files with similar timesync |
1341 |
|
// |
1342 |
|
if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0); |
1343 |
|
UInt_t upperts = tsync-(obt0/1000)+5; |
1344 |
|
UInt_t lowerts = tsync-(obt0/1000)-5; |
1345 |
|
oss.str(""); |
1346 |
|
oss << "select GL_RAW.BOOT_NUMBER from GL_TIMESYNC LEFT JOIN GL_RAW ON GL_RAW.ID = GL_TIMESYNC.ID_RAW where TIMESYNC-(OBT0/1000)<" |
1347 |
|
<< upperts |
1348 |
|
<< " AND TIMESYNC-(OBT0/1000)>" |
1349 |
|
<< lowerts |
1350 |
|
<< " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;"; |
1351 |
|
result = conn->Query(oss.str().c_str()); |
1352 |
|
if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str()); |
1353 |
|
// |
1354 |
|
if ( !result ) throw -4;; |
1355 |
|
found = true; |
1356 |
|
if ( result->GetRowCount()<3 ){ |
1357 |
|
if ( IsDebug() ) printf(" AGH! no results!\n"); |
1358 |
|
found = false; |
1359 |
|
} else { |
1360 |
|
row = result->Next(); |
1361 |
|
bn = (UInt_t)atoll(row->GetField(0)); |
1362 |
|
for ( Int_t r=1; r<result->GetRowCount() ;r++){ |
1363 |
|
if ( !row ) throw -4; |
1364 |
|
if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0)); |
1365 |
|
if ( bn != (UInt_t)atoll(row->GetField(0)) ){ |
1366 |
|
if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0))); |
1367 |
|
found = false; |
1368 |
|
}; |
1369 |
|
row = result->Next(); |
1370 |
|
}; |
1371 |
}; |
}; |
1372 |
|
}; |
1373 |
|
// |
1374 |
|
Int_t sgn = 0; |
1375 |
|
// |
1376 |
|
if ( !found ){ |
1377 |
|
throw -29; |
1378 |
} else { |
} else { |
1379 |
if ( !this->GetBOOTnumber() ) return(2); |
if ( afound ){ |
1380 |
|
this->SetBOOTnumber(bn); |
1381 |
|
sgn = 8; |
1382 |
|
}; |
1383 |
}; |
}; |
1384 |
// |
// |
1385 |
oss.str(""); |
oss.str(""); |
1389 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
1390 |
// |
// |
1391 |
delete result; |
delete result; |
1392 |
return(0); |
return(sgn); |
1393 |
}; |
}; |
1394 |
|
|
1395 |
/** |
/** |
1707 |
delete result; |
delete result; |
1708 |
// |
// |
1709 |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
1710 |
|
if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n"); |
1711 |
return(signal); |
return(signal); |
1712 |
}; |
}; |
1713 |
|
|
1821 |
// |
// |
1822 |
} else { |
} else { |
1823 |
if ( IsDebug() ) printf(" The run is already present in the fragment table \n"); |
if ( IsDebug() ) printf(" The run is already present in the fragment table \n"); |
1824 |
|
return; |
1825 |
}; |
}; |
1826 |
// |
// |
1827 |
// |
// |
2820 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
2821 |
// |
// |
2822 |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section); |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section); |
2823 |
if ( fromtime < 1150867000 ){ //1150866904 |
if ( fromtime < 1150871000 ){ //1150866904 |
2824 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
2825 |
fromtime = 0;// the first flight calibration was taken at about 1156429100 s, this line allow to analyze first runs in raw mode |
fromtime = 0;// the first flight calibration was taken at about 1156429100 s, this line allow to analyze first runs in raw mode |
2826 |
}; |
}; |
2965 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
2966 |
// |
// |
2967 |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n"); |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n"); |
2968 |
if ( fromtime < 1150867000 ) fromtime = 0; // the first flight calibration was taken at about 1150863300 s, this line allows to analyze first runs in raw mode |
if ( fromtime < 1150871000 ) fromtime = 0; // the first flight calibration was taken at about 1150863300 s, this line allows to analyze first runs in raw mode |
2969 |
// |
// |
2970 |
oss.str(""); |
oss.str(""); |
2971 |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
3094 |
pkt1 = ph1->GetCounter(); |
pkt1 = ph1->GetCounter(); |
3095 |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
3096 |
// |
// |
3097 |
valid = 1; |
// valid = 1; |
3098 |
// |
// // |
3099 |
if ( caltrk1->unpackError != 0 && caltrk1->good0 == 0 ) valid = 0;// CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT1 |
// if ( caltrk1->unpackError != 0 && caltrk1->good0 == 0 ) valid = 0;// CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT1 |
3100 |
// |
// |
3101 |
// |
// |
3102 |
if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){ |
if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){ |
3103 |
// |
// |
3104 |
if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1); |
if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1); |
3105 |
|
// |
3106 |
|
valid = ValidateTrkCalib( caltrk1, eh1 ); |
3107 |
|
if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl; |
3108 |
// |
// |
3109 |
// Do we have the second calibration packet? |
// Do we have the second calibration packet? |
3110 |
// |
// |
3121 |
obt2 = ph2->GetOrbitalTime(); |
obt2 = ph2->GetOrbitalTime(); |
3122 |
pkt2 = ph2->GetCounter(); |
pkt2 = ph2->GetCounter(); |
3123 |
// |
// |
3124 |
if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2 |
// if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2 |
3125 |
// |
// |
3126 |
} else { |
} else { |
3127 |
// |
// |
3152 |
// |
// |
3153 |
if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2); |
if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2); |
3154 |
// |
// |
3155 |
|
UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 ); |
3156 |
|
if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl; |
3157 |
|
valid = valid & valid2; |
3158 |
|
// |
3159 |
// Handle good calib |
// Handle good calib |
3160 |
// |
// |
3161 |
this->HandleTRK_CALIB(true,true); |
this->HandleTRK_CALIB(true,true); |
3328 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
3329 |
// |
// |
3330 |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n"); |
if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n"); |
3331 |
if ( fromtime < 1150867000 ){ |
if ( fromtime < 1150871000 ){ |
3332 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
3333 |
fromtime = 0;// the first flight calibration was taken at about 1156429100 s, this line allow to analyze first runs in raw mode |
fromtime = 0;// the first flight calibration was taken at about 1156429100 s, this line allow to analyze first runs in raw mode |
3334 |
}; |
}; |
4424 |
// |
// |
4425 |
}; |
}; |
4426 |
}; |
}; |
4427 |
|
|
4428 |
|
/** |
4429 |
|
* |
4430 |
|
* Rearrange calibration tables |
4431 |
|
* |
4432 |
|
**/ |
4433 |
|
UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){ |
4434 |
|
|
4435 |
|
Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0}; |
4436 |
|
UInt_t timeaftercalib=120000; //2000; |
4437 |
|
// ---------- |
4438 |
|
// Check CRCs |
4439 |
|
// ---------- |
4440 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
4441 |
|
if( caltrk->crc_hcal[ipkt] )return 0; // :-( |
4442 |
|
for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-( |
4443 |
|
} |
4444 |
|
// ----------------------- |
4445 |
|
// Check missing packets: |
4446 |
|
// ----------------------- |
4447 |
|
// Readout order: |
4448 |
|
// ------------------ |
4449 |
|
// DSP packet board |
4450 |
|
// ------------------ |
4451 |
|
// 12 0 1 |
4452 |
|
// 10 1 1 |
4453 |
|
// 8 2 1 |
4454 |
|
// 4 3 1 |
4455 |
|
// 6 4 1 |
4456 |
|
// 2 5 1 |
4457 |
|
// ------------------ |
4458 |
|
// 11 0 2 |
4459 |
|
// 9 1 2 |
4460 |
|
// 7 2 2 |
4461 |
|
// 3 3 2 |
4462 |
|
// 5 4 2 |
4463 |
|
// 1 5 2 |
4464 |
|
// ------------------ |
4465 |
|
// ------------------------------------------------- |
4466 |
|
// Check if it is first or second calibration packet |
4467 |
|
// ------------------------------------------------- |
4468 |
|
UInt_t build=0; |
4469 |
|
TString classname = caltrk->GetName(); |
4470 |
|
UInt_t base=0; |
4471 |
|
UInt_t mask=0; |
4472 |
|
if(classname.Contains("CalibTrk1Event")){ |
4473 |
|
base=12; |
4474 |
|
mask=0x03F000; |
4475 |
|
} |
4476 |
|
if(classname.Contains("CalibTrk2Event")){ |
4477 |
|
base=18; |
4478 |
|
mask=0xFC0000; |
4479 |
|
} |
4480 |
|
// ------------------------------------------------- |
4481 |
|
// Count number of packets and set build variable |
4482 |
|
// ------------------------------------------------- |
4483 |
|
Int_t npkts=0; |
4484 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
4485 |
|
if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){ |
4486 |
|
npkts++; |
4487 |
|
build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) ); |
4488 |
|
} |
4489 |
|
} |
4490 |
|
// if( npkts==6 )return 1; // :-) |
4491 |
|
|
4492 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl; |
4493 |
|
|
4494 |
|
// ----------------------------------------------- |
4495 |
|
// If missing packets: check the acq configuration |
4496 |
|
// (some DSPs might be excluded from acquisition) |
4497 |
|
// ----------------------------------------------- |
4498 |
|
|
4499 |
|
// ----------------------------------------------- |
4500 |
|
// retrieve the first run header after calib |
4501 |
|
// ----------------------------------------------- |
4502 |
|
PacketType *pctp; |
4503 |
|
EventCounter *cod; |
4504 |
|
cod = eh->GetCounter(); |
4505 |
|
Int_t irun = cod->Get(pctp->RunHeader); |
4506 |
|
TTree *rh=(TTree*)file->Get("RunHeader"); |
4507 |
|
if ( !rh || rh->IsZombie() ) throw -17; |
4508 |
|
if( rh->GetEntries() == irun ){ |
4509 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1) -- cannot validate :-( "<<endl; |
4510 |
|
return 0; // :-( |
4511 |
|
} |
4512 |
|
|
4513 |
|
RunHeaderEvent *run = 0; |
4514 |
|
EventHeader *hrun = 0; |
4515 |
|
rh->SetBranchAddress("RunHeader", &run); |
4516 |
|
rh->SetBranchAddress("Header", &hrun); |
4517 |
|
rh->GetEntry(irun); |
4518 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl; |
4519 |
|
|
4520 |
|
if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){ |
4521 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl; |
4522 |
|
return 0; // :-( |
4523 |
|
} |
4524 |
|
|
4525 |
|
if( !run->RM_ACQ_AFTER_CALIB ){ |
4526 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0 -- cannot validate :-( "<<endl; |
4527 |
|
return 0; // :-( |
4528 |
|
} |
4529 |
|
|
4530 |
|
UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime()); |
4531 |
|
if( dtime > timeaftercalib ){ |
4532 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl; |
4533 |
|
return 0; // :-( |
4534 |
|
} |
4535 |
|
|
4536 |
|
|
4537 |
|
|
4538 |
|
if( (run->ACQ_BUILD_INFO & mask) != build ){ |
4539 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl; |
4540 |
|
return 0; // :-( |
4541 |
|
} |
4542 |
|
return 1; // :-) |
4543 |
|
|
4544 |
|
} |
4545 |
|
|
4546 |
|
/** |
4547 |
|
* |
4548 |
|
* Check the DB (only for overlapping runs at the moment) |
4549 |
|
* |
4550 |
|
**/ |
4551 |
|
UInt_t PamelaDBOperations::Check(){ |
4552 |
|
// |
4553 |
|
UInt_t test = 0; |
4554 |
|
// |
4555 |
|
UInt_t thisrht = 0; |
4556 |
|
UInt_t thisrtt = 0; |
4557 |
|
UInt_t thisid = 0; |
4558 |
|
UInt_t prevrht = 0; |
4559 |
|
UInt_t prevrtt = 0; |
4560 |
|
UInt_t previd = 0; |
4561 |
|
// |
4562 |
|
UInt_t prevl0id = 0; |
4563 |
|
UInt_t thisl0id = 0; |
4564 |
|
// |
4565 |
|
stringstream oss; |
4566 |
|
TSQLResult *result = 0; |
4567 |
|
TSQLRow *row = 0; |
4568 |
|
TSQLResult *result2 = 0; |
4569 |
|
TSQLRow *row2 = 0; |
4570 |
|
oss.str(""); |
4571 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN order by RUNHEADER_TIME asc;"; |
4572 |
|
// oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;"; |
4573 |
|
result = conn->Query(oss.str().c_str()); |
4574 |
|
// |
4575 |
|
if ( !result ) throw -4;; |
4576 |
|
// |
4577 |
|
row = result->Next(); |
4578 |
|
// |
4579 |
|
while ( row ){ |
4580 |
|
thisid = (UInt_t)atoll(row->GetField(0)); |
4581 |
|
thisl0id = (UInt_t)atoll(row->GetField(1)); |
4582 |
|
thisrht = (UInt_t)atoll(row->GetField(2)); |
4583 |
|
thisrtt = (UInt_t)atoll(row->GetField(3)); |
4584 |
|
// |
4585 |
|
// if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){ |
4586 |
|
// if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){ |
4587 |
|
if ( (thisrht < prevrtt) && (thisrht != prevrht) ){ |
4588 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4589 |
|
printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid); |
4590 |
|
TString prevf = ""; |
4591 |
|
TString thisf = ""; |
4592 |
|
oss.str(""); |
4593 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4594 |
|
result2 = conn->Query(oss.str().c_str()); |
4595 |
|
if ( !result2 ) throw -4;; |
4596 |
|
row2 = result2->Next(); |
4597 |
|
prevf = (TString)row2->GetField(0); |
4598 |
|
oss.str(""); |
4599 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4600 |
|
result2 = conn->Query(oss.str().c_str()); |
4601 |
|
if ( !result2 ) throw -4;; |
4602 |
|
row2 = result2->Next(); |
4603 |
|
thisf = (TString)row2->GetField(0); |
4604 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4605 |
|
test = 1; |
4606 |
|
}; |
4607 |
|
// |
4608 |
|
if ( (thisrtt < prevrht) && (thisrht != prevrht) ){ |
4609 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4610 |
|
printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid); |
4611 |
|
TString prevf = ""; |
4612 |
|
TString thisf = ""; |
4613 |
|
oss.str(""); |
4614 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4615 |
|
result2 = conn->Query(oss.str().c_str()); |
4616 |
|
if ( !result2 ) throw -4;; |
4617 |
|
row2 = result2->Next(); |
4618 |
|
prevf = (TString)row2->GetField(0); |
4619 |
|
oss.str(""); |
4620 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4621 |
|
result2 = conn->Query(oss.str().c_str()); |
4622 |
|
if ( !result2 ) throw -4;; |
4623 |
|
row2 = result2->Next(); |
4624 |
|
thisf = (TString)row2->GetField(0); |
4625 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4626 |
|
test = 1; |
4627 |
|
}; |
4628 |
|
// |
4629 |
|
if ( (thisrht > thisrtt) && (thisrht != prevrht) ){ |
4630 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4631 |
|
printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid); |
4632 |
|
TString prevf = ""; |
4633 |
|
TString thisf = ""; |
4634 |
|
oss.str(""); |
4635 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4636 |
|
result2 = conn->Query(oss.str().c_str()); |
4637 |
|
if ( !result2 ) throw -4;; |
4638 |
|
row2 = result2->Next(); |
4639 |
|
prevf = (TString)row2->GetField(0); |
4640 |
|
oss.str(""); |
4641 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4642 |
|
result2 = conn->Query(oss.str().c_str()); |
4643 |
|
if ( !result2 ) throw -4;; |
4644 |
|
row2 = result2->Next(); |
4645 |
|
thisf = (TString)row2->GetField(0); |
4646 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4647 |
|
test = 1; |
4648 |
|
}; |
4649 |
|
|
4650 |
|
// |
4651 |
|
prevrht = thisrht; |
4652 |
|
prevrtt = thisrtt; |
4653 |
|
previd = thisid; |
4654 |
|
prevl0id = thisl0id; |
4655 |
|
row = result->Next(); |
4656 |
|
}; |
4657 |
|
// |
4658 |
|
return(test); |
4659 |
|
// |
4660 |
|
}; |