14 |
#include <TSQLRow.h> |
#include <TSQLRow.h> |
15 |
#include <TTree.h> |
#include <TTree.h> |
16 |
#include <TGraph.h> |
#include <TGraph.h> |
17 |
#include <TDatime.h> |
#include <TTimeStamp.h> |
18 |
#include <TF1.h> |
#include <TF1.h> |
19 |
// |
// |
20 |
#include <EventHeader.h> |
#include <EventHeader.h> |
56 |
* @param debug debug flag. |
* @param debug debug flag. |
57 |
* @param tlefilename ascii file with TLE 3 line elements. |
* @param tlefilename ascii file with TLE 3 line elements. |
58 |
*/ |
*/ |
59 |
PamelaDBOperations::PamelaDBOperations(TString host, TString user, TString password, TString filerawname, TString filerootname, UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t debug, TString tlefilename){ |
PamelaDBOperations::PamelaDBOperations(TString host, TString user, TString password, TString filerawname, TString filerootname, UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t debug, TString tlefilename, UInt_t dwinput, Bool_t staticp){ |
60 |
// |
// |
61 |
|
STATIC=false; |
62 |
|
if ( staticp ) STATIC=true; |
63 |
// |
// |
64 |
SetConnection(host,user,password); |
SetConnection(host,user,password); |
65 |
// |
// |
81 |
INSERT_ROOT = !filerootname.IsNull(); |
INSERT_ROOT = !filerootname.IsNull(); |
82 |
if( INSERT_ROOT ){ |
if( INSERT_ROOT ){ |
83 |
this->SetRootName(filerootname); |
this->SetRootName(filerootname); |
84 |
this->SetOrbitNo(); |
this->SetOrbitNo(dwinput); |
85 |
file = TFile::Open(this->GetRootName().Data()); |
file = TFile::Open(this->GetRootName().Data()); |
86 |
} else { |
} else { |
87 |
this->SetRootName(""); |
this->SetRootName(""); |
114 |
// |
// |
115 |
void PamelaDBOperations::CheckValidate(Long64_t olderthan){ |
void PamelaDBOperations::CheckValidate(Long64_t olderthan){ |
116 |
clean_time = new TDatime(); |
clean_time = new TDatime(); |
117 |
|
// |
118 |
if(olderthan >= 0){ |
if(olderthan >= 0){ |
119 |
VALIDATE = true; |
VALIDATE = true; |
120 |
UInt_t timelim = 0; |
UInt_t timelim = 0; |
121 |
timelim = (UInt_t)clean_time->Convert() - olderthan; |
timelim = (UInt_t)clean_time->Convert(true) - olderthan; |
122 |
clean_time->Set(timelim,false); |
clean_time->Set(timelim,false); |
123 |
}; |
}; |
124 |
}; |
}; |
159 |
}; |
}; |
160 |
|
|
161 |
/** |
/** |
162 |
|
* Set the autoboot flag |
163 |
|
* |
164 |
|
*/ |
165 |
|
void PamelaDBOperations::SetAutoBoot(Bool_t dbg){ |
166 |
|
AUTOBOOT = dbg; |
167 |
|
}; |
168 |
|
|
169 |
|
/** |
170 |
* Set the nofrag flag |
* Set the nofrag flag |
171 |
* |
* |
172 |
*/ |
*/ |
217 |
/** |
/** |
218 |
* Store the downlink orbit number from filename. |
* Store the downlink orbit number from filename. |
219 |
*/ |
*/ |
220 |
void PamelaDBOperations::SetOrbitNo(){ |
void PamelaDBOperations::SetOrbitNo(UInt_t dwinput){ |
221 |
dworbit = 0; |
dworbit = 0; |
222 |
|
// |
223 |
|
if ( dwinput ){ |
224 |
|
dworbit = dwinput; |
225 |
|
if ( IsDebug() ) printf(" Downlink orbit given by hand: %i \n",dworbit); |
226 |
|
return; |
227 |
|
}; |
228 |
|
// |
229 |
TString name = this->GetRootFile(); |
TString name = this->GetRootFile(); |
230 |
Int_t nlength = name.Length(); |
Int_t nlength = name.Length(); |
231 |
if ( nlength < 5 ) return; |
if ( nlength < 5 ) return; |
263 |
tlefilename = str; |
tlefilename = str; |
264 |
}; |
}; |
265 |
|
|
266 |
|
TString PamelaDBOperations::GetRawPath(){ |
267 |
|
if ( STATIC ){ |
268 |
|
return((TString)gSystem->DirName(filerawname.Data())+'/'); |
269 |
|
} else { |
270 |
|
return((TString)gSystem->ExpandPathName("$PAM_RAW")+'/'); |
271 |
|
}; |
272 |
|
}; |
273 |
|
|
274 |
|
TString PamelaDBOperations::GetRootPath(){ |
275 |
|
if ( STATIC ){ |
276 |
|
return((TString)gSystem->DirName(filerootname.Data())+'/'); |
277 |
|
} else { |
278 |
|
return((TString)gSystem->ExpandPathName("$PAM_L0")+'/'); |
279 |
|
}; |
280 |
|
}; |
281 |
|
|
282 |
/** |
/** |
283 |
* Store the olderthan variable |
* Store the olderthan variable |
284 |
* @param olderthan |
* @param olderthan |
295 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
296 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
297 |
oss.str(""); |
oss.str(""); |
298 |
oss << "SELECT ID FROM GL_RAW WHERE " |
if ( STATIC ){ |
299 |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
oss << "SELECT ID FROM GL_RAW WHERE " |
300 |
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
301 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
302 |
|
} else { |
303 |
|
oss << "SELECT ID FROM GL_RAW WHERE " |
304 |
|
<< " PATH = '$PAM_RAW' AND " |
305 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
306 |
|
} |
307 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
308 |
if ( result == NULL ) throw -4; |
if ( result == NULL ) throw -4; |
309 |
row = result->Next(); |
row = result->Next(); |
334 |
UInt_t nevent = 0; |
UInt_t nevent = 0; |
335 |
UInt_t pktlast = 0; |
UInt_t pktlast = 0; |
336 |
UInt_t obtlast = 0; |
UInt_t obtlast = 0; |
337 |
UInt_t t_pktlast = 0; |
Long64_t t_pktlast = 0LL; |
338 |
UInt_t t_obtlast = 0; |
// UInt_t t_obtlast = 0; |
339 |
UInt_t upperpkt2 = 0; |
Long64_t t_obtlast = 0LL; |
340 |
ULong64_t upperobt2 = 0; |
Long64_t upperpkt2 = 0LL; |
341 |
|
Long64_t upperobt2 = 0LL; |
342 |
UInt_t zomp = 0; |
UInt_t zomp = 0; |
343 |
UInt_t jump = 50000; // was 5000 |
UInt_t jump = 50000; // was 5000 |
344 |
EventCounter *code=0; |
EventCounter *code=0; |
345 |
// |
// |
346 |
UInt_t deltapkt = 5000; |
Long64_t deltapkt = 5000LL; |
347 |
ULong64_t deltaobt = 50000; |
Long64_t deltaobt = 50000LL; |
348 |
// |
// |
349 |
// pcksList packetsNames; |
pcksList packetsNames; |
350 |
// pcksList::iterator Iter; |
pcksList::iterator Iter; |
351 |
// getPacketsNames(packetsNames); |
getPacketsNames(packetsNames); |
352 |
// |
// |
353 |
pktfirst = 0; |
pktfirst = 0; |
354 |
obtfirst = 0; |
obtfirst = 0; |
355 |
|
ppktfirst = 0; |
356 |
|
pobtfirst = 0; |
357 |
// |
// |
358 |
TTree *T = 0; |
TTree *T = 0; |
359 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
367 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
368 |
pktfirst = ph->GetCounter(); |
pktfirst = ph->GetCounter(); |
369 |
obtfirst = ph->GetOrbitalTime(); |
obtfirst = ph->GetOrbitalTime(); |
370 |
|
ppktfirst = pktfirst; |
371 |
|
pobtfirst = obtfirst; |
372 |
// |
// |
373 |
// code = eh->GetCounter(); |
code = eh->GetCounter(); |
374 |
// UInt_t en = 0; |
UInt_t en = 0; |
375 |
// for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
376 |
// en = code->Get(GetPacketType(*Iter)); |
en = code->Get(GetPacketType(*Iter)); |
377 |
// if ( en ) printf(" Packet type is %s, entries: %i \n",*Iter,en); |
if ( !strcmp("CalibCalPed",*Iter) || !strcmp("CalibTrk1",*Iter) || !strcmp("CalibTrk2",*Iter) || !strcmp("CalibS4",*Iter) ){ |
378 |
//}; |
if ( IsDebug() ) printf(" We have a calibration before the first physic packet: packet type is %s, entries: %i \n",*Iter,en); |
379 |
|
// |
380 |
|
TTree *TC = 0; |
381 |
|
TC = (TTree*)file->Get("CalibCalPed"); |
382 |
|
if ( !TC || TC->IsZombie() ) throw -16; |
383 |
|
EventHeader *ehc = 0; |
384 |
|
PscuHeader *phc = 0; |
385 |
|
TC->SetBranchAddress("Header", &ehc); |
386 |
|
TC->GetEntry(0); |
387 |
|
phc = ehc->GetPscuHeader(); |
388 |
|
pktfirst = phc->GetCounter(); |
389 |
|
obtfirst = phc->GetOrbitalTime(); |
390 |
|
// |
391 |
|
}; |
392 |
|
}; |
393 |
// |
// |
394 |
T->GetEntry(nevent-1); |
T->GetEntry(nevent-1); |
395 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
400 |
upperobt = OBT(obtlast); |
upperobt = OBT(obtlast); |
401 |
upperentry = nevent-1; |
upperentry = nevent-1; |
402 |
// |
// |
403 |
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); |
404 |
// |
// |
405 |
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); |
406 |
// |
// |
407 |
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); |
408 |
// |
// |
409 |
if ( !nevent ) return(2); |
if ( !nevent ) return(2); |
410 |
// |
// |
411 |
if ( nevent < 2 ) return(4); |
if ( nevent < 2 ) return(4); |
412 |
|
if ( nevent < jump ) jump = 1; |
413 |
|
// if ( nevent < jump ) jump = int(nevent/10); |
414 |
|
// if ( !jump ) jump = 1; |
415 |
// |
// |
416 |
if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) ){ |
if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) && nevent > deltapkt ){ |
417 |
// |
// |
418 |
|
if ( IsDebug() ) printf(" starting jump %i \n",jump); |
419 |
// go back |
// go back |
420 |
zomp = nevent - 2; |
zomp = nevent - 2; |
421 |
// |
// |
437 |
upperpkt2 = PKT(ph->GetCounter()); |
upperpkt2 = PKT(ph->GetCounter()); |
438 |
upperobt2 = OBT(ph->GetOrbitalTime()); |
upperobt2 = OBT(ph->GetOrbitalTime()); |
439 |
// |
// |
440 |
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) ){ |
441 |
|
if ( IsDebug() ) printf(" .-. upperpkt2 %lld upperobt2 %lld \n",upperpkt2,upperobt2); |
442 |
|
if ( IsDebug() ) printf(" .-. upperpkt %lld t_pktlast %lld upperobt %lld t_obtlast %lld \n",upperpkt,t_pktlast,upperobt,t_obtlast); |
443 |
|
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); |
444 |
|
throw -13; |
445 |
|
}; |
446 |
// |
// |
447 |
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 ){ |
448 |
zomp = i + jump + 1; |
zomp = i + jump + 1; |
449 |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
450 |
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); |
451 |
break; |
break; |
452 |
}; |
}; |
453 |
// |
// |
482 |
// |
// |
483 |
rhev = rh->GetEntries(); |
rhev = rh->GetEntries(); |
484 |
rtev = rt->GetEntries(); |
rtev = rt->GetEntries(); |
485 |
UInt_t sobtt = 0; |
Long64_t sobtt = 0LL; |
486 |
UInt_t sobth = 0; |
Long64_t sobth = 0LL; |
487 |
UInt_t spktt = 0; |
Long64_t spktt = 0LL; |
488 |
UInt_t spkth = 0; |
Long64_t spkth = 0LL; |
489 |
UInt_t pktt = 0; |
Long64_t pktt = 0LL; |
490 |
ULong64_t obtt = 0; |
Long64_t obtt = 0LL; |
491 |
UInt_t pkth = 0; |
Long64_t pkth = 0LL; |
492 |
ULong64_t obth = 0; |
Long64_t obth = 0LL; |
493 |
// |
// |
494 |
T->GetEntry(upperentry); |
if ( rhev || rtev ){ |
495 |
code = eh->GetCounter(); |
|
496 |
Int_t lasttrail = code->Get(pctp->RunTrailer); |
T->GetEntry(upperentry); |
497 |
Int_t lasthead = code->Get(pctp->RunHeader); |
code = eh->GetCounter(); |
498 |
if ( lasttrail < rtev ){ |
Int_t lasttrail = code->Get(pctp->RunTrailer); |
499 |
rt->GetEntry(lasttrail); |
Int_t lasthead = code->Get(pctp->RunHeader); |
500 |
pht = eht->GetPscuHeader(); |
if ( lasttrail < rtev ){ |
501 |
pktt = PKT(pht->GetCounter()); |
rt->GetEntry(lasttrail); |
502 |
obtt = OBT(pht->GetOrbitalTime()); |
pht = eht->GetPscuHeader(); |
503 |
}; |
pktt = PKT(pht->GetCounter()); |
504 |
// |
obtt = OBT(pht->GetOrbitalTime()); |
505 |
if ( lasthead < rhev ){ |
}; |
506 |
rh->GetEntry(lasthead); |
// |
507 |
phh = ehh->GetPscuHeader(); |
if ( lasthead < rhev ){ |
508 |
pkth = PKT(phh->GetCounter()); |
rh->GetEntry(lasthead); |
509 |
obth = OBT(phh->GetOrbitalTime()); |
phh = ehh->GetPscuHeader(); |
510 |
}; |
pkth = PKT(phh->GetCounter()); |
511 |
// |
obth = OBT(phh->GetOrbitalTime()); |
512 |
if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
}; |
513 |
if ( pkth > upperpkt && obth > upperobt ){ |
// |
514 |
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(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
515 |
upperpkt = pkth; |
if ( pkth > upperpkt && obth > upperobt ){ |
516 |
upperobt = obth; |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt); |
517 |
rhev = lasthead+1; |
upperpkt = pkth; |
518 |
} else { |
upperobt = obth; |
519 |
rhev = lasthead; |
rhev = lasthead+1; |
520 |
}; |
} else { |
521 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
rhev = lasthead; |
522 |
// |
}; |
523 |
if ( IsDebug() ) printf(" rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
524 |
if ( pktt > upperpkt && obtt > upperobt ){ |
// |
525 |
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(" rtev beforev %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
526 |
upperpkt = pktt; |
if ( pktt > upperpkt && obtt > upperobt ){ |
527 |
upperobt = obtt; |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt); |
528 |
rtev = lasttrail+1; |
upperpkt = pktt; |
529 |
} else { |
upperobt = obtt; |
530 |
rtev = lasttrail; |
rtev = lasttrail+1; |
531 |
}; |
} else { |
532 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
rtev = lasttrail; |
533 |
// goto kikko; |
}; |
534 |
// |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
535 |
// |
// goto kikko; |
536 |
// Check if runtrailer/runheader are within lower limits |
// |
537 |
// |
// |
538 |
// |
// Check if runtrailer/runheader are within lower limits |
539 |
pkth = 0; |
// |
540 |
obth = 0; |
// |
541 |
spkth = 0; |
pkth = 0LL; |
542 |
sobth = 0; |
obth = 0LL; |
543 |
for (Int_t k=0; k<rhev; k++){ |
spkth = 0LL; |
544 |
if ( k > 0 ){ |
sobth = 0LL; |
545 |
spkth = pkth; |
for (Int_t k=0; k<rhev; k++){ |
546 |
sobth = obth; |
if ( k > 0 ){ |
547 |
}; |
spkth = pkth; |
548 |
rh->GetEntry(k); |
sobth = obth; |
549 |
phh = ehh->GetPscuHeader(); |
}; |
550 |
pkth = PKT(phh->GetCounter()); |
rh->GetEntry(k); |
551 |
obth = OBT(phh->GetOrbitalTime()); |
phh = ehh->GetPscuHeader(); |
552 |
// |
pkth = PKT(phh->GetCounter()); |
553 |
// if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth); |
obth = OBT(phh->GetOrbitalTime()); |
554 |
// |
// |
555 |
if ( pkth < spkth && obth < sobth ){ |
// if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth); |
556 |
if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i \n",rhev); |
// |
557 |
// |
if ( pkth < spkth && obth < sobth ){ |
558 |
rhev = k-1; |
if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i \n",rhev); |
559 |
rh->GetEntry(rhev); |
// |
560 |
pkth = spkth; |
rhev = k-1; |
561 |
obth = sobth; |
rh->GetEntry(rhev); |
562 |
// |
pkth = spkth; |
563 |
UInt_t evbefh = 0; |
obth = sobth; |
564 |
code = ehh->GetCounter(); |
// |
565 |
evbefh = code->Get(pctp->Physics); |
UInt_t evbefh = 0; |
566 |
if ( evbefh >= 0 ){ |
code = ehh->GetCounter(); |
567 |
T->GetEntry(evbefh); |
evbefh = code->Get(pctp->Physics); |
568 |
ph = eh->GetPscuHeader(); |
if ( evbefh >= 0 ){ |
569 |
t_pktlast = PKT(ph->GetCounter()); |
T->GetEntry(evbefh); |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
|
|
if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump |
|
|
upperpkt = pkth; |
|
|
upperobt = obth; |
|
|
upperentry = evbefh-1; |
|
|
} else { |
|
|
while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){ |
|
|
evbefh++; |
|
|
T->GetEntry(evbefh); |
|
|
ph = eh->GetPscuHeader(); |
|
|
t_pktlast = PKT(ph->GetCounter()); |
|
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
|
|
}; |
|
|
T->GetEntry(evbefh-1); |
|
570 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
571 |
upperpkt = PKT(ph->GetCounter()); |
t_pktlast = PKT(ph->GetCounter()); |
572 |
upperobt = OBT(ph->GetOrbitalTime()); |
t_obtlast = OBT(ph->GetOrbitalTime()); |
573 |
upperentry = evbefh-1; |
if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump |
574 |
}; |
upperpkt = pkth; |
575 |
|
upperobt = obth; |
576 |
|
upperentry = evbefh-1; |
577 |
|
} else { |
578 |
|
while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){ |
579 |
|
evbefh++; |
580 |
|
T->GetEntry(evbefh); |
581 |
|
ph = eh->GetPscuHeader(); |
582 |
|
t_pktlast = PKT(ph->GetCounter()); |
583 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
584 |
|
}; |
585 |
|
T->GetEntry(evbefh-1); |
586 |
|
ph = eh->GetPscuHeader(); |
587 |
|
upperpkt = PKT(ph->GetCounter()); |
588 |
|
upperobt = OBT(ph->GetOrbitalTime()); |
589 |
|
upperentry = evbefh-1; |
590 |
|
}; |
591 |
|
}; |
592 |
|
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
593 |
|
goto kikko0; |
594 |
}; |
}; |
|
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
|
|
goto kikko0; |
|
|
}; |
|
|
}; |
|
|
kikko0: |
|
|
// |
|
|
// |
|
|
// |
|
|
pktt = 0; |
|
|
obtt = 0; |
|
|
spktt = 0; |
|
|
sobtt = 0; |
|
|
for (Int_t k=0; k<rtev; k++){ |
|
|
if ( k > 0 ){ |
|
|
spktt = pktt; |
|
|
sobtt = obtt; |
|
595 |
}; |
}; |
596 |
rt->GetEntry(k); |
kikko0: |
|
pht = eht->GetPscuHeader(); |
|
|
pktt = PKT(pht->GetCounter()); |
|
|
obtt = OBT(pht->GetOrbitalTime()); |
|
597 |
// |
// |
|
// if ( IsDebug() ) printf(" k %i rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt); |
|
598 |
// |
// |
599 |
if ( pktt < spktt && obtt < sobtt ){ |
// |
600 |
if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev); |
pktt = 0LL; |
601 |
|
obtt = 0LL; |
602 |
|
spktt = 0LL; |
603 |
|
sobtt = 0LL; |
604 |
|
for (Int_t k=0; k<rtev; k++){ |
605 |
|
if ( k > 0 ){ |
606 |
|
spktt = pktt; |
607 |
|
sobtt = obtt; |
608 |
|
}; |
609 |
|
rt->GetEntry(k); |
610 |
|
pht = eht->GetPscuHeader(); |
611 |
|
pktt = PKT(pht->GetCounter()); |
612 |
|
obtt = OBT(pht->GetOrbitalTime()); |
613 |
// |
// |
614 |
rtev = k-1; |
// if ( IsDebug() ) printf(" k %i rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt); |
|
rt->GetEntry(rtev); |
|
|
pktt = spktt; |
|
|
obtt = sobtt; |
|
|
if ( IsDebug() ) printf(" lasttrail %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
|
615 |
// |
// |
616 |
UInt_t evbeft = 0; |
if ( pktt < spktt && obtt < sobtt ){ |
617 |
code = eht->GetCounter(); |
if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev); |
618 |
evbeft = code->Get(pctp->Physics); |
// |
619 |
if ( evbeft >= 0 ){ |
rtev = k-1; |
620 |
T->GetEntry(evbeft); |
rt->GetEntry(rtev); |
621 |
ph = eh->GetPscuHeader(); |
pktt = spktt; |
622 |
t_pktlast = PKT(ph->GetCounter()); |
obtt = sobtt; |
623 |
t_obtlast = OBT(ph->GetOrbitalTime()); |
if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
624 |
if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump |
// |
625 |
upperpkt = pktt; |
UInt_t evbeft = 0; |
626 |
upperobt = obtt; |
code = eht->GetCounter(); |
627 |
upperentry = evbeft-1; |
evbeft = code->Get(pctp->Physics); |
628 |
} else { |
if ( evbeft >= 0 ){ |
629 |
while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){ |
T->GetEntry(evbeft); |
630 |
evbeft++; |
ph = eh->GetPscuHeader(); |
631 |
T->GetEntry(evbeft); |
t_pktlast = PKT(ph->GetCounter()); |
632 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
633 |
|
if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump |
634 |
|
upperpkt = pktt; |
635 |
|
upperobt = obtt; |
636 |
|
upperentry = evbeft-1; |
637 |
|
} else { |
638 |
|
while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){ |
639 |
|
evbeft++; |
640 |
|
T->GetEntry(evbeft); |
641 |
|
ph = eh->GetPscuHeader(); |
642 |
|
t_pktlast = PKT(ph->GetCounter()); |
643 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
644 |
|
}; |
645 |
|
T->GetEntry(evbeft-1); |
646 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
647 |
t_pktlast = PKT(ph->GetCounter()); |
upperpkt = PKT(ph->GetCounter()); |
648 |
t_obtlast = OBT(ph->GetOrbitalTime()); |
upperobt = OBT(ph->GetOrbitalTime()); |
649 |
|
upperentry = evbeft-1; |
650 |
}; |
}; |
|
T->GetEntry(evbeft-1); |
|
|
ph = eh->GetPscuHeader(); |
|
|
upperpkt = PKT(ph->GetCounter()); |
|
|
upperobt = OBT(ph->GetOrbitalTime()); |
|
|
upperentry = evbeft-1; |
|
651 |
}; |
}; |
652 |
|
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
653 |
|
goto kikko; |
654 |
|
// break; |
655 |
|
// |
656 |
}; |
}; |
657 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
// |
|
goto kikko; |
|
|
// break; |
|
|
// |
|
658 |
}; |
}; |
659 |
// |
// |
660 |
}; |
kikko: |
661 |
// |
// |
662 |
kikko: |
T->GetEntry(upperentry); |
663 |
// |
code = eh->GetCounter(); |
664 |
T->GetEntry(upperentry); |
lasttrail = code->Get(pctp->RunTrailer); |
665 |
code = eh->GetCounter(); |
lasthead = code->Get(pctp->RunHeader); |
666 |
lasttrail = code->Get(pctp->RunTrailer); |
if ( lasttrail < rtev ){ |
667 |
lasthead = code->Get(pctp->RunHeader); |
rt->GetEntry(lasttrail); |
668 |
if ( lasttrail < rtev ){ |
pht = eht->GetPscuHeader(); |
669 |
rt->GetEntry(lasttrail); |
pktt = PKT(pht->GetCounter()); |
670 |
pht = eht->GetPscuHeader(); |
obtt = OBT(pht->GetOrbitalTime()); |
671 |
pktt = PKT(pht->GetCounter()); |
}; |
672 |
obtt = OBT(pht->GetOrbitalTime()); |
// |
673 |
}; |
if ( lasthead < rhev ){ |
674 |
// |
rh->GetEntry(lasthead); |
675 |
if ( lasthead < rhev ){ |
phh = ehh->GetPscuHeader(); |
676 |
rh->GetEntry(lasthead); |
pkth = PKT(phh->GetCounter()); |
677 |
phh = ehh->GetPscuHeader(); |
obth = OBT(phh->GetOrbitalTime()); |
678 |
pkth = PKT(phh->GetCounter()); |
}; |
679 |
obth = OBT(phh->GetOrbitalTime()); |
// |
680 |
}; |
if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
681 |
// |
if ( pkth > upperpkt && obth > upperobt ){ |
682 |
if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,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); |
683 |
if ( pkth > upperpkt && obth > upperobt ){ |
upperpkt = pkth; |
684 |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt); |
upperobt = obth; |
685 |
upperpkt = pkth; |
rhev = lasthead+1; |
686 |
upperobt = obth; |
} else { |
687 |
rhev = lasthead+1; |
rhev = lasthead; |
688 |
} else { |
}; |
689 |
rhev = lasthead; |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
690 |
}; |
// |
691 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
if ( IsDebug() ) printf(" rtev beforev %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
692 |
// |
if ( pktt > upperpkt && obtt > upperobt ){ |
693 |
if ( IsDebug() ) printf(" rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",rtev,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); |
694 |
if ( pktt > upperpkt && obtt > upperobt ){ |
upperpkt = pktt; |
695 |
if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt); |
upperobt = obtt; |
696 |
upperpkt = pktt; |
rtev = lasttrail+1; |
697 |
upperobt = obtt; |
} else { |
698 |
rtev = lasttrail+1; |
rtev = lasttrail; |
699 |
} else { |
}; |
700 |
rtev = lasttrail; |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
701 |
|
// |
702 |
}; |
}; |
|
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
|
703 |
// |
// |
704 |
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); |
705 |
// |
// |
706 |
return(0); |
return(0); |
707 |
} |
} |
835 |
if( !conn ) throw -1; |
if( !conn ) throw -1; |
836 |
bool connect = conn->IsConnected(); |
bool connect = conn->IsConnected(); |
837 |
if( !connect ) throw -1; |
if( !connect ) throw -1; |
838 |
|
// |
839 |
|
if ( IsDebug() ) printf("\n DB INFORMATIONS:\n SQL: %s Version: %s Host %s Port %i \n\n",conn->GetDBMS(),conn->ServerInfo(),conn->GetHost(),conn->GetPort()); |
840 |
|
// |
841 |
if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27; |
if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27; |
842 |
|
// |
843 |
|
// set DB timezone to UTC |
844 |
|
// |
845 |
|
stringstream oss; |
846 |
|
// |
847 |
|
oss.str(""); |
848 |
|
oss << "SET time_zone='+0:00';"; |
849 |
|
TSQLResult *result = 0; |
850 |
|
result = conn->Query(oss.str().c_str()); |
851 |
|
if ( !result ) throw -10; |
852 |
|
oss.str(""); |
853 |
|
oss << "SET wait_timeout=173000;"; |
854 |
|
conn->Query(oss.str().c_str()); |
855 |
|
// |
856 |
}; |
}; |
857 |
|
|
858 |
/** |
/** |
859 |
* Return the correct packet number if we went back to zero |
* Return the correct packet number if we went back to zero |
860 |
*/ |
*/ |
861 |
UInt_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
862 |
// |
// |
863 |
// 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,ppktfirst,(UInt_t)(16777214/2)); |
864 |
// |
// |
865 |
if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2) ) return((pkt_num+16777215)); |
if ( pkt_num < (ppktfirst/2) && ppktfirst > (16777214/2) ){ |
866 |
|
if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL); |
867 |
|
return((Long64_t)pkt_num+16777215LL); |
868 |
|
}; |
869 |
// |
// |
870 |
if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){ |
if ( pkt_num > ppktfirst*2 && pkt_num > (16777214/2) ){ |
871 |
if ( (pkt_num-16777215) < 0 ){ |
if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL); |
872 |
return((16777215-pkt_num)); |
return((Long64_t)pkt_num-16777215LL); |
|
} else { |
|
|
return((pkt_num-16777215)); |
|
|
}; |
|
873 |
}; |
}; |
874 |
// |
// |
875 |
return(pkt_num); |
if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)pkt_num); |
876 |
|
return((Long64_t)pkt_num); |
877 |
// |
// |
878 |
}; |
}; |
879 |
|
|
880 |
/** |
/** |
881 |
* Return the correct On Board Time if we went back to zero |
* Return the correct On Board Time if we went back to zero |
882 |
*/ |
*/ |
883 |
ULong64_t PamelaDBOperations::OBT(UInt_t obt){ |
Long64_t PamelaDBOperations::OBT(UInt_t obt){ |
884 |
// |
// |
885 |
if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((ULong64_t)(obt+numeric_limits<UInt_t>::max())); |
if ( IsDebug() ) printf(" obt conversion: ob is %u obtfirst is %u (numeric_limits<UInt_t>::max()/2) is %u \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2)); |
886 |
// |
// |
887 |
if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ) return((Long64_t)(obt+numeric_limits<UInt_t>::max())); |
888 |
if ( (obt-numeric_limits<UInt_t>::max()) < 0 ){ |
// |
889 |
return((ULong64_t)(numeric_limits<UInt_t>::max()-obt)); |
if ( obt > (pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
890 |
} else { |
return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
|
return((ULong64_t)(obt-numeric_limits<UInt_t>::max())); |
|
|
}; |
|
891 |
}; |
}; |
892 |
// |
// |
893 |
return((ULong64_t)obt); |
return((Long64_t)obt); |
894 |
}; |
}; |
895 |
|
|
896 |
/** |
/** |
1005 |
if ( idr ) return(1); |
if ( idr ) return(1); |
1006 |
// |
// |
1007 |
oss.str(""); |
oss.str(""); |
1008 |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
if ( STATIC ){ |
1009 |
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
1010 |
|
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
1011 |
|
} else { |
1012 |
|
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')"; |
1013 |
|
}; |
1014 |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
1015 |
// |
// |
1016 |
idr = this->SetID_RAW(); |
idr = this->SetID_RAW(); |
1035 |
if ( this->GetID_RAW() == 0 ) throw -11; |
if ( this->GetID_RAW() == 0 ) throw -11; |
1036 |
// |
// |
1037 |
oss.str(""); |
oss.str(""); |
1038 |
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='" |
1039 |
<< this->GetRawFile().Data() << "';"; |
<< this->GetRawFile().Data() << "';"; |
1040 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
1041 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1044 |
// |
// |
1045 |
if ( !row ){ |
if ( !row ){ |
1046 |
oss.str(""); |
oss.str(""); |
1047 |
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< " |
1048 |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
1049 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
1050 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1053 |
if ( !row ) throw -10; |
if ( !row ) throw -10; |
1054 |
}; |
}; |
1055 |
// |
// |
1056 |
t0 = (UInt_t)TDatime(row->GetField(0)).Convert(); |
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); |
1057 |
|
t0 = (UInt_t)tu.GetSec(); |
1058 |
|
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)); |
1059 |
|
// |
1060 |
/* |
/* |
1061 |
* Verify that the TIMESYNC have been not already processed |
* Verify that the TIMESYNC have been not already processed |
1062 |
*/ |
*/ |
1071 |
if (result == NULL) throw -10; |
if (result == NULL) throw -10; |
1072 |
row = result->Next(); |
row = result->Next(); |
1073 |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
1074 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
1075 |
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; |
1076 |
|
// |
1077 |
|
tsync = (UInt_t)atoll(row->GetField(2)); |
1078 |
|
obt0 = (UInt_t)atoll(row->GetField(1)); |
1079 |
|
// |
1080 |
return(1); |
return(1); |
1081 |
}; |
}; |
1082 |
// |
// |
1270 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
1271 |
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()); |
1272 |
// |
// |
1273 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
1274 |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
1275 |
// |
// |
1276 |
|
tsync = TSYNC; |
1277 |
|
obt0 = OBT; |
1278 |
|
// |
1279 |
delete result; |
delete result; |
1280 |
return(signal); |
return(signal); |
1281 |
} |
} |
1291 |
UInt_t idtimesync = 0; |
UInt_t idtimesync = 0; |
1292 |
// |
// |
1293 |
oss.str(""); |
oss.str(""); |
1294 |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
if ( STATIC ){ |
1295 |
<< " LEFT JOIN GL_ROOT " |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
1296 |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
<< " LEFT JOIN GL_ROOT " |
1297 |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
1298 |
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
1299 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
1300 |
|
} else { |
1301 |
|
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
1302 |
|
<< " LEFT JOIN GL_ROOT " |
1303 |
|
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
1304 |
|
<< " WHERE GL_RAW.PATH = '$PAM_RAW' AND " |
1305 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
1306 |
|
}; |
1307 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1308 |
// |
// |
1309 |
if ( !result ) throw -12; |
if ( !result ) throw -12; |
1330 |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
1331 |
// |
// |
1332 |
oss.str(""); |
oss.str(""); |
1333 |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
if ( STATIC ){ |
1334 |
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
1335 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
1336 |
|
} else { |
1337 |
|
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
1338 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')"; |
1339 |
|
}; |
1340 |
// |
// |
1341 |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
1342 |
// |
// |
1363 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
1364 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
1365 |
oss.str(""); |
oss.str(""); |
1366 |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
if ( STATIC ){ |
1367 |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
1368 |
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
1369 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
1370 |
|
} else { |
1371 |
|
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
1372 |
|
<< " PATH = '$PAM_RAW' AND " |
1373 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
1374 |
|
}; |
1375 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1376 |
// |
// |
1377 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
1394 |
VarDumpEvent *vde = 0; |
VarDumpEvent *vde = 0; |
1395 |
VarDumpRecord *vdr = 0; |
VarDumpRecord *vdr = 0; |
1396 |
// |
// |
1397 |
|
Bool_t found = false; |
1398 |
trDumpEv->SetBranchAddress("VarDump", &vde); |
trDumpEv->SetBranchAddress("VarDump", &vde); |
1399 |
if ( trDumpEv->GetEntries() > 0 ){ |
if ( trDumpEv->GetEntries() > 0 ){ |
1400 |
Bool_t found = false; |
found = false; |
1401 |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
1402 |
trDumpEv->GetEntry(i); |
trDumpEv->GetEntry(i); |
1403 |
vde->Records->GetEntries(); |
// vde->Records->GetEntries(); |
1404 |
if ( vde->Records->GetEntries()>0 ){ |
if ( vde->Records->GetEntries()>5 ){ |
1405 |
found = true; |
found = true; |
1406 |
goto fill; |
goto fill; |
1407 |
}; |
}; |
1414 |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
1415 |
// |
// |
1416 |
} else { |
} else { |
1417 |
if ( !this->GetBOOTnumber() ) return(4); |
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4); |
1418 |
}; |
}; |
1419 |
} else { |
} else { |
1420 |
if ( !this->GetBOOTnumber() ) return(2); |
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2); |
1421 |
|
}; |
1422 |
|
// |
1423 |
|
UInt_t bn = 0; |
1424 |
|
Bool_t afound = false; |
1425 |
|
if ( !found && this->AutoBoot()){ |
1426 |
|
afound = true; |
1427 |
|
// |
1428 |
|
// Search for other files with similar timesync |
1429 |
|
// |
1430 |
|
if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0); |
1431 |
|
UInt_t upperts = tsync-(obt0/1000)+5; |
1432 |
|
UInt_t lowerts = tsync-(obt0/1000)-5; |
1433 |
|
oss.str(""); |
1434 |
|
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)<" |
1435 |
|
<< upperts |
1436 |
|
<< " AND TIMESYNC-(OBT0/1000)>" |
1437 |
|
<< lowerts |
1438 |
|
<< " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;"; |
1439 |
|
result = conn->Query(oss.str().c_str()); |
1440 |
|
if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str()); |
1441 |
|
// |
1442 |
|
if ( !result ) throw -4;; |
1443 |
|
found = true; |
1444 |
|
if ( result->GetRowCount()<3 ){ |
1445 |
|
if ( IsDebug() ) printf(" AGH! no results!\n"); |
1446 |
|
found = false; |
1447 |
|
} else { |
1448 |
|
row = result->Next(); |
1449 |
|
bn = (UInt_t)atoll(row->GetField(0)); |
1450 |
|
for ( Int_t r=1; r<result->GetRowCount() ;r++){ |
1451 |
|
if ( !row ) throw -4; |
1452 |
|
if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0)); |
1453 |
|
if ( bn != (UInt_t)atoll(row->GetField(0)) ){ |
1454 |
|
if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0))); |
1455 |
|
found = false; |
1456 |
|
}; |
1457 |
|
row = result->Next(); |
1458 |
|
}; |
1459 |
|
}; |
1460 |
|
}; |
1461 |
|
// |
1462 |
|
Int_t sgn = 0; |
1463 |
|
// |
1464 |
|
if ( !found && !BOOTNO ){ |
1465 |
|
throw -29; |
1466 |
|
} else { |
1467 |
|
if ( afound ){ |
1468 |
|
this->SetBOOTnumber(bn); |
1469 |
|
sgn = 8; |
1470 |
|
}; |
1471 |
}; |
}; |
1472 |
// |
// |
1473 |
oss.str(""); |
oss.str(""); |
1477 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
1478 |
// |
// |
1479 |
delete result; |
delete result; |
1480 |
return(0); |
return(sgn); |
1481 |
}; |
}; |
1482 |
|
|
1483 |
/** |
/** |
1795 |
delete result; |
delete result; |
1796 |
// |
// |
1797 |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
1798 |
|
if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n"); |
1799 |
return(signal); |
return(signal); |
1800 |
}; |
}; |
1801 |
|
|
1909 |
// |
// |
1910 |
} else { |
} else { |
1911 |
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"); |
1912 |
|
return; |
1913 |
}; |
}; |
1914 |
// |
// |
1915 |
// |
// |
2908 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
2909 |
// |
// |
2910 |
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); |
2911 |
if ( fromtime < 1150863400 ){ |
if ( fromtime < 1150871000 ){ //1150866904 |
2912 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
2913 |
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 |
2914 |
}; |
}; |
2971 |
// |
// |
2972 |
} else { |
} else { |
2973 |
// |
// |
2974 |
if ( IsDebug() ) printf(" Repetead calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt); |
if ( IsDebug() ) printf(" Repeated calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt); |
2975 |
// |
// |
2976 |
}; |
}; |
2977 |
// |
// |
3053 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
3054 |
// |
// |
3055 |
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"); |
3056 |
if ( fromtime < 1150863400 ) 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 |
3057 |
// |
// |
3058 |
oss.str(""); |
oss.str(""); |
3059 |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
3182 |
pkt1 = ph1->GetCounter(); |
pkt1 = ph1->GetCounter(); |
3183 |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
3184 |
// |
// |
3185 |
valid = 1; |
// valid = 1; |
3186 |
// |
// // |
3187 |
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 |
3188 |
// |
// |
3189 |
// |
// |
3190 |
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) ){ |
3191 |
// |
// |
3192 |
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); |
3193 |
|
// |
3194 |
|
valid = ValidateTrkCalib( caltrk1, eh1 ); |
3195 |
|
if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl; |
3196 |
// |
// |
3197 |
// Do we have the second calibration packet? |
// Do we have the second calibration packet? |
3198 |
// |
// |
3209 |
obt2 = ph2->GetOrbitalTime(); |
obt2 = ph2->GetOrbitalTime(); |
3210 |
pkt2 = ph2->GetCounter(); |
pkt2 = ph2->GetCounter(); |
3211 |
// |
// |
3212 |
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 |
3213 |
// |
// |
3214 |
} else { |
} else { |
3215 |
// |
// |
3240 |
// |
// |
3241 |
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); |
3242 |
// |
// |
3243 |
|
UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 ); |
3244 |
|
if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl; |
3245 |
|
valid = valid & valid2; |
3246 |
|
// |
3247 |
// Handle good calib |
// Handle good calib |
3248 |
// |
// |
3249 |
this->HandleTRK_CALIB(true,true); |
this->HandleTRK_CALIB(true,true); |
3416 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
3417 |
// |
// |
3418 |
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"); |
3419 |
if ( fromtime < 1150863400 ){ |
if ( fromtime < 1150871000 ){ |
3420 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
3421 |
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 |
3422 |
}; |
}; |
4512 |
// |
// |
4513 |
}; |
}; |
4514 |
}; |
}; |
4515 |
|
|
4516 |
|
/** |
4517 |
|
* |
4518 |
|
* Rearrange calibration tables |
4519 |
|
* |
4520 |
|
**/ |
4521 |
|
UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){ |
4522 |
|
|
4523 |
|
Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0}; |
4524 |
|
UInt_t timeaftercalib=120000; //2000; |
4525 |
|
// ---------- |
4526 |
|
// Check CRCs |
4527 |
|
// ---------- |
4528 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
4529 |
|
if( caltrk->crc_hcal[ipkt] )return 0; // :-( |
4530 |
|
for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-( |
4531 |
|
} |
4532 |
|
// ----------------------- |
4533 |
|
// Check missing packets: |
4534 |
|
// ----------------------- |
4535 |
|
// Readout order: |
4536 |
|
// ------------------ |
4537 |
|
// DSP packet board |
4538 |
|
// ------------------ |
4539 |
|
// 12 0 1 |
4540 |
|
// 10 1 1 |
4541 |
|
// 8 2 1 |
4542 |
|
// 4 3 1 |
4543 |
|
// 6 4 1 |
4544 |
|
// 2 5 1 |
4545 |
|
// ------------------ |
4546 |
|
// 11 0 2 |
4547 |
|
// 9 1 2 |
4548 |
|
// 7 2 2 |
4549 |
|
// 3 3 2 |
4550 |
|
// 5 4 2 |
4551 |
|
// 1 5 2 |
4552 |
|
// ------------------ |
4553 |
|
// ------------------------------------------------- |
4554 |
|
// Check if it is first or second calibration packet |
4555 |
|
// ------------------------------------------------- |
4556 |
|
UInt_t build=0; |
4557 |
|
TString classname = caltrk->GetName(); |
4558 |
|
UInt_t base=0; |
4559 |
|
UInt_t mask=0; |
4560 |
|
if(classname.Contains("CalibTrk1Event")){ |
4561 |
|
base=12; |
4562 |
|
mask=0x03F000; |
4563 |
|
} |
4564 |
|
if(classname.Contains("CalibTrk2Event")){ |
4565 |
|
base=18; |
4566 |
|
mask=0xFC0000; |
4567 |
|
} |
4568 |
|
// ------------------------------------------------- |
4569 |
|
// Count number of packets and set build variable |
4570 |
|
// ------------------------------------------------- |
4571 |
|
Int_t npkts=0; |
4572 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
4573 |
|
if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){ |
4574 |
|
npkts++; |
4575 |
|
build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) ); |
4576 |
|
} |
4577 |
|
} |
4578 |
|
// if( npkts==6 )return 1; // :-) |
4579 |
|
|
4580 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl; |
4581 |
|
|
4582 |
|
// ----------------------------------------------- |
4583 |
|
// If missing packets: check the acq configuration |
4584 |
|
// (some DSPs might be excluded from acquisition) |
4585 |
|
// ----------------------------------------------- |
4586 |
|
|
4587 |
|
// ----------------------------------------------- |
4588 |
|
// retrieve the first run header after calib |
4589 |
|
// ----------------------------------------------- |
4590 |
|
PacketType *pctp; |
4591 |
|
EventCounter *cod; |
4592 |
|
cod = eh->GetCounter(); |
4593 |
|
Int_t irun = cod->Get(pctp->RunHeader); |
4594 |
|
TTree *rh=(TTree*)file->Get("RunHeader"); |
4595 |
|
if ( !rh || rh->IsZombie() ) throw -17; |
4596 |
|
if( rh->GetEntries() == irun ){ |
4597 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1) -- cannot validate :-( "<<endl; |
4598 |
|
return 0; // :-( |
4599 |
|
} |
4600 |
|
|
4601 |
|
RunHeaderEvent *run = 0; |
4602 |
|
EventHeader *hrun = 0; |
4603 |
|
rh->SetBranchAddress("RunHeader", &run); |
4604 |
|
rh->SetBranchAddress("Header", &hrun); |
4605 |
|
rh->GetEntry(irun); |
4606 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl; |
4607 |
|
|
4608 |
|
if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){ |
4609 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl; |
4610 |
|
return 0; // :-( |
4611 |
|
} |
4612 |
|
|
4613 |
|
if( !run->RM_ACQ_AFTER_CALIB ){ |
4614 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0 -- cannot validate :-( "<<endl; |
4615 |
|
return 0; // :-( |
4616 |
|
} |
4617 |
|
|
4618 |
|
UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime()); |
4619 |
|
if( dtime > timeaftercalib ){ |
4620 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl; |
4621 |
|
return 0; // :-( |
4622 |
|
} |
4623 |
|
|
4624 |
|
|
4625 |
|
|
4626 |
|
if( (run->ACQ_BUILD_INFO & mask) != build ){ |
4627 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl; |
4628 |
|
return 0; // :-( |
4629 |
|
} |
4630 |
|
return 1; // :-) |
4631 |
|
|
4632 |
|
} |
4633 |
|
|
4634 |
|
/** |
4635 |
|
* |
4636 |
|
* Check the DB (only for overlapping runs at the moment) |
4637 |
|
* |
4638 |
|
**/ |
4639 |
|
UInt_t PamelaDBOperations::Check(){ |
4640 |
|
// |
4641 |
|
UInt_t test = 0; |
4642 |
|
// |
4643 |
|
UInt_t thisrht = 0; |
4644 |
|
UInt_t thisrtt = 0; |
4645 |
|
UInt_t thisid = 0; |
4646 |
|
UInt_t prevrht = 0; |
4647 |
|
UInt_t prevrtt = 0; |
4648 |
|
UInt_t previd = 0; |
4649 |
|
// |
4650 |
|
UInt_t prevl0id = 0; |
4651 |
|
UInt_t thisl0id = 0; |
4652 |
|
// |
4653 |
|
stringstream oss; |
4654 |
|
TSQLResult *result = 0; |
4655 |
|
TSQLRow *row = 0; |
4656 |
|
TSQLResult *result2 = 0; |
4657 |
|
TSQLRow *row2 = 0; |
4658 |
|
oss.str(""); |
4659 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN order by RUNHEADER_TIME asc;"; |
4660 |
|
// oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;"; |
4661 |
|
result = conn->Query(oss.str().c_str()); |
4662 |
|
// |
4663 |
|
if ( !result ) throw -4;; |
4664 |
|
// |
4665 |
|
row = result->Next(); |
4666 |
|
// |
4667 |
|
while ( row ){ |
4668 |
|
thisid = (UInt_t)atoll(row->GetField(0)); |
4669 |
|
thisl0id = (UInt_t)atoll(row->GetField(1)); |
4670 |
|
thisrht = (UInt_t)atoll(row->GetField(2)); |
4671 |
|
thisrtt = (UInt_t)atoll(row->GetField(3)); |
4672 |
|
// |
4673 |
|
// if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){ |
4674 |
|
// if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){ |
4675 |
|
if ( (thisrht < prevrtt) && (thisrht != prevrht) ){ |
4676 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4677 |
|
printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid); |
4678 |
|
TString prevf = ""; |
4679 |
|
TString thisf = ""; |
4680 |
|
oss.str(""); |
4681 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4682 |
|
result2 = conn->Query(oss.str().c_str()); |
4683 |
|
if ( !result2 ) throw -4;; |
4684 |
|
row2 = result2->Next(); |
4685 |
|
prevf = (TString)row2->GetField(0); |
4686 |
|
oss.str(""); |
4687 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4688 |
|
result2 = conn->Query(oss.str().c_str()); |
4689 |
|
if ( !result2 ) throw -4;; |
4690 |
|
row2 = result2->Next(); |
4691 |
|
thisf = (TString)row2->GetField(0); |
4692 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4693 |
|
test = 1; |
4694 |
|
}; |
4695 |
|
// |
4696 |
|
if ( (thisrtt < prevrht) && (thisrht != prevrht) ){ |
4697 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4698 |
|
printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid); |
4699 |
|
TString prevf = ""; |
4700 |
|
TString thisf = ""; |
4701 |
|
oss.str(""); |
4702 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4703 |
|
result2 = conn->Query(oss.str().c_str()); |
4704 |
|
if ( !result2 ) throw -4;; |
4705 |
|
row2 = result2->Next(); |
4706 |
|
prevf = (TString)row2->GetField(0); |
4707 |
|
oss.str(""); |
4708 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4709 |
|
result2 = conn->Query(oss.str().c_str()); |
4710 |
|
if ( !result2 ) throw -4;; |
4711 |
|
row2 = result2->Next(); |
4712 |
|
thisf = (TString)row2->GetField(0); |
4713 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4714 |
|
test = 1; |
4715 |
|
}; |
4716 |
|
// |
4717 |
|
if ( (thisrht > thisrtt) && (thisrht != prevrht) ){ |
4718 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
4719 |
|
printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid); |
4720 |
|
TString prevf = ""; |
4721 |
|
TString thisf = ""; |
4722 |
|
oss.str(""); |
4723 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
4724 |
|
result2 = conn->Query(oss.str().c_str()); |
4725 |
|
if ( !result2 ) throw -4;; |
4726 |
|
row2 = result2->Next(); |
4727 |
|
prevf = (TString)row2->GetField(0); |
4728 |
|
oss.str(""); |
4729 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
4730 |
|
result2 = conn->Query(oss.str().c_str()); |
4731 |
|
if ( !result2 ) throw -4;; |
4732 |
|
row2 = result2->Next(); |
4733 |
|
thisf = (TString)row2->GetField(0); |
4734 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
4735 |
|
test = 1; |
4736 |
|
}; |
4737 |
|
|
4738 |
|
// |
4739 |
|
prevrht = thisrht; |
4740 |
|
prevrtt = thisrtt; |
4741 |
|
previd = thisid; |
4742 |
|
prevl0id = thisl0id; |
4743 |
|
row = result->Next(); |
4744 |
|
}; |
4745 |
|
// |
4746 |
|
return(test); |
4747 |
|
// |
4748 |
|
}; |