| 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(""); |
| 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(); |
| 346 |
Long64_t deltapkt = 5000LL; |
Long64_t deltapkt = 5000LL; |
| 347 |
Long64_t deltaobt = 50000LL; |
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 %i pkt_num %i \n",obtfirst,pktfirst); |
if ( IsDebug() ) printf(" First entries are: OBT %u pkt_num %u \n",obtfirst,pktfirst); |
| 404 |
// |
// |
| 405 |
if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld 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(32); |
| 408 |
// |
// |
| 409 |
if ( !nevent ) return(2); |
if ( !nevent ) return(64); |
| 410 |
// |
// |
| 411 |
if ( nevent < 2 ) return(4); |
if ( nevent < 2 ) return(128); |
| 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 |
// |
// |
| 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 %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i); |
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); |
| 451 |
break; |
break; |
| 452 |
}; |
}; |
| 453 |
// |
// |
| 491 |
Long64_t pkth = 0LL; |
Long64_t pkth = 0LL; |
| 492 |
Long64_t obth = 0LL; |
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 %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
}; |
| 513 |
if ( pkth > upperpkt && obth > upperobt ){ |
// |
| 514 |
if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \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 %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
rhev = lasthead; |
| 522 |
// |
}; |
| 523 |
if ( IsDebug() ) printf(" rtev beforev %i pt %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \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 = 0LL; |
// |
| 540 |
obth = 0LL; |
// |
| 541 |
spkth = 0LL; |
pkth = 0LL; |
| 542 |
sobth = 0LL; |
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 %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
|
|
goto kikko0; |
|
| 595 |
}; |
}; |
| 596 |
}; |
kikko0: |
| 597 |
kikko0: |
// |
|
// |
|
|
// |
|
|
// |
|
|
pktt = 0LL; |
|
|
obtt = 0LL; |
|
|
spktt = 0LL; |
|
|
sobtt = 0LL; |
|
|
for (Int_t k=0; k<rtev; k++){ |
|
|
if ( k > 0 ){ |
|
|
spktt = pktt; |
|
|
sobtt = obtt; |
|
|
}; |
|
|
rt->GetEntry(k); |
|
|
pht = eht->GetPscuHeader(); |
|
|
pktt = PKT(pht->GetCounter()); |
|
|
obtt = OBT(pht->GetOrbitalTime()); |
|
| 598 |
// |
// |
|
// if ( IsDebug() ) printf(" k %i rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt); |
|
| 599 |
// |
// |
| 600 |
if ( pktt < spktt && obtt < sobtt ){ |
pktt = 0LL; |
| 601 |
if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev); |
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 %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld oh %lld uppero %lld \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 %lld upperp %lld oh %lld uppero %lld \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 %lld upperp %lld oh %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \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 %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
|
| 703 |
// |
// |
| 704 |
if ( IsDebug() ) printf(" Upper limits are: OBT %lld pkt_num %lld 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 |
// |
// |
| 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 |
// set DB timezone to UTC |
| 849 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 850 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 851 |
if ( !result ) throw -10; |
if ( !result ) throw -10; |
| 852 |
|
oss.str(""); |
| 853 |
|
oss << "SET wait_timeout=173000;"; |
| 854 |
|
conn->Query(oss.str().c_str()); |
| 855 |
// |
// |
| 856 |
}; |
}; |
| 857 |
|
|
| 860 |
*/ |
*/ |
| 861 |
Long64_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 (UInt_t)(16777214/2)) is %u \n",pkt_num,pktfirst,(UInt_t)(16777214/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) ){ |
if ( pkt_num < (ppktfirst/2) && ppktfirst > (16777214/2) ){ |
| 866 |
if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL); |
if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL); |
| 867 |
return((Long64_t)pkt_num+16777215LL); |
return((Long64_t)pkt_num+16777215LL); |
| 868 |
}; |
}; |
| 869 |
// |
// |
| 870 |
if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){ |
if ( pkt_num > ((Long64_t)ppktfirst*2) && pkt_num > (16777214/2) ){ |
| 871 |
if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL); |
if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL); |
| 872 |
return((Long64_t)pkt_num-16777215LL); |
return((Long64_t)pkt_num-16777215LL); |
| 873 |
}; |
}; |
| 882 |
*/ |
*/ |
| 883 |
Long64_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((Long64_t)(obt+numeric_limits<UInt_t>::max())); |
if ( IsDebug() ) printf(" obt conversion: obt 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 < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){ |
| 888 |
|
if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)(obt+numeric_limits<UInt_t>::max())); |
| 889 |
|
return((Long64_t)(obt+numeric_limits<UInt_t>::max())); |
| 890 |
|
}; |
| 891 |
// |
// |
| 892 |
if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
if ( obt > ((Long64_t)pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
| 893 |
|
if ( IsDebug() ) printf(" pobtfirst*2 %lld \n",((Long64_t)pobtfirst*2)); |
| 894 |
|
if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
| 895 |
return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
| 896 |
}; |
}; |
| 897 |
// |
// |
| 898 |
|
if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)obt); |
| 899 |
return((Long64_t)obt); |
return((Long64_t)obt); |
| 900 |
}; |
}; |
| 901 |
|
|
| 1011 |
if ( idr ) return(1); |
if ( idr ) return(1); |
| 1012 |
// |
// |
| 1013 |
oss.str(""); |
oss.str(""); |
| 1014 |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
if ( STATIC ){ |
| 1015 |
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
| 1016 |
|
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
| 1017 |
|
} else { |
| 1018 |
|
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')"; |
| 1019 |
|
}; |
| 1020 |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
| 1021 |
// |
// |
| 1022 |
idr = this->SetID_RAW(); |
idr = this->SetID_RAW(); |
| 1035 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 1036 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 1037 |
UInt_t t0 = 0; |
UInt_t t0 = 0; |
| 1038 |
|
Int_t signal = 0; |
| 1039 |
|
// |
| 1040 |
|
signal = this->SetUpperLimits(); |
| 1041 |
// |
// |
| 1042 |
stringstream oss; |
stringstream oss; |
| 1043 |
// |
// |
| 1090 |
}; |
}; |
| 1091 |
// |
// |
| 1092 |
TTree *T = 0; |
TTree *T = 0; |
|
Int_t signal = 0; |
|
| 1093 |
// |
// |
| 1094 |
UInt_t nevent = 0; |
UInt_t nevent = 0; |
| 1095 |
UInt_t recEntries = 0; |
UInt_t recEntries = 0; |
| 1299 |
UInt_t idtimesync = 0; |
UInt_t idtimesync = 0; |
| 1300 |
// |
// |
| 1301 |
oss.str(""); |
oss.str(""); |
| 1302 |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
if ( STATIC ){ |
| 1303 |
<< " LEFT JOIN GL_ROOT " |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
| 1304 |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
<< " LEFT JOIN GL_ROOT " |
| 1305 |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
| 1306 |
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
| 1307 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
| 1308 |
|
} else { |
| 1309 |
|
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
| 1310 |
|
<< " LEFT JOIN GL_ROOT " |
| 1311 |
|
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
| 1312 |
|
<< " WHERE GL_RAW.PATH = '$PAM_RAW' AND " |
| 1313 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
| 1314 |
|
}; |
| 1315 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1316 |
// |
// |
| 1317 |
if ( !result ) throw -12; |
if ( !result ) throw -12; |
| 1338 |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
| 1339 |
// |
// |
| 1340 |
oss.str(""); |
oss.str(""); |
| 1341 |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
if ( STATIC ){ |
| 1342 |
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
| 1343 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
| 1344 |
|
} else { |
| 1345 |
|
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
| 1346 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')"; |
| 1347 |
|
}; |
| 1348 |
// |
// |
| 1349 |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
| 1350 |
// |
// |
| 1371 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 1372 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 1373 |
oss.str(""); |
oss.str(""); |
| 1374 |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
if ( STATIC ){ |
| 1375 |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
| 1376 |
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
| 1377 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 1378 |
|
} else { |
| 1379 |
|
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
| 1380 |
|
<< " PATH = '$PAM_RAW' AND " |
| 1381 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 1382 |
|
}; |
| 1383 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1384 |
// |
// |
| 1385 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
| 1469 |
// |
// |
| 1470 |
Int_t sgn = 0; |
Int_t sgn = 0; |
| 1471 |
// |
// |
| 1472 |
if ( !found ){ |
if ( !found && !BOOTNO ){ |
| 1473 |
throw -29; |
throw -29; |
| 1474 |
} else { |
} else { |
| 1475 |
if ( afound ){ |
if ( afound ){ |
| 1498 |
stringstream oss; |
stringstream oss; |
| 1499 |
oss.str(""); |
oss.str(""); |
| 1500 |
// |
// |
| 1501 |
signal = this->SetUpperLimits(); |
// signal = this->SetUpperLimits(); |
| 1502 |
// |
// |
| 1503 |
// loop on runheader and runtrailer events |
// loop on runheader and runtrailer events |
| 1504 |
// |
// |
| 1516 |
rt->SetBranchAddress("RunTrailer", &runt); |
rt->SetBranchAddress("RunTrailer", &runt); |
| 1517 |
rt->SetBranchAddress("Header", &eht); |
rt->SetBranchAddress("Header", &eht); |
| 1518 |
// |
// |
| 1519 |
|
TTree *T = (TTree*)file->Get("Physics"); |
| 1520 |
|
if ( !T || T->IsZombie() ) throw -16; |
| 1521 |
|
EventHeader *eh = 0; |
| 1522 |
|
T->SetBranchAddress("Header", &eh); |
| 1523 |
|
// |
| 1524 |
|
if ( !(rh->GetEntries()) && !(rt->GetEntries()) && !(T->GetEntries()) ) return(16); |
| 1525 |
|
// |
| 1526 |
UInt_t obtt = 0; |
UInt_t obtt = 0; |
| 1527 |
UInt_t obth = 0; |
UInt_t obth = 0; |
| 1528 |
UInt_t pktt = 0; |
UInt_t pktt = 0; |
| 1556 |
// |
// |
| 1557 |
if ( !ptt && !(ptht+1) ){ |
if ( !ptt && !(ptht+1) ){ |
| 1558 |
// |
// |
| 1559 |
if ( IsDebug() ) printf(" Piece of run at the beginning of the file %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Piece of run at the beginning of the file %u %u %u \n",ptht,pth,ptt); |
| 1560 |
// |
// |
| 1561 |
this->HandleRunFragments(true,false,0,(evbeft-1)); |
this->HandleRunFragments(true,false,0,(evbeft-1)); |
| 1562 |
// |
// |
| 1563 |
// |
// |
| 1564 |
} else if ( pth == ptht ){ |
} else if ( pth == ptht ){ |
| 1565 |
// |
// |
| 1566 |
if ( IsDebug() ) printf(" Missing header %i %i %i\n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Missing header %u %u %u\n",ptht,pth,ptt); |
| 1567 |
// |
// |
| 1568 |
if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
| 1569 |
rt->GetEntry(ptt-1); |
rt->GetEntry(ptt-1); |
| 1572 |
rt->GetEntry(ptt); |
rt->GetEntry(ptt); |
| 1573 |
pht = eht->GetPscuHeader(); |
pht = eht->GetPscuHeader(); |
| 1574 |
// |
// |
| 1575 |
if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer %u %u %u \n",ptht,pth,ptt); |
| 1576 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1577 |
// |
// |
| 1578 |
this->HandleMissingHoT(true,false,evbefh,evbeft-1); |
this->HandleMissingHoT(true,false,evbefh,evbeft-1); |
| 1579 |
// |
// |
| 1580 |
} else { |
} else { |
| 1581 |
// |
// |
| 1582 |
if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %i %i %i\n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %u %u %u\n",ptht,pth,ptt); |
| 1583 |
// |
// |
| 1584 |
rh->GetEntry(ptht); |
rh->GetEntry(ptht); |
| 1585 |
phh = ehh->GetPscuHeader(); |
phh = ehh->GetPscuHeader(); |
| 1587 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1588 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
| 1589 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1590 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1591 |
// |
// |
| 1592 |
// handle this run |
// handle this run |
| 1593 |
// |
// |
| 1597 |
// |
// |
| 1598 |
if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){ |
if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){ |
| 1599 |
// |
// |
| 1600 |
if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER \n"); |
if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER evbefh = %u \n",evbefh); |
|
// |
|
|
this->HandleRunFragments(true,true,0,(evbefh-1)); |
|
| 1601 |
// |
// |
| 1602 |
|
if ( evbefh == 0 ) { |
| 1603 |
|
// |
| 1604 |
|
signal = 8; |
| 1605 |
|
if ( IsDebug() ) printf(" Not supported yet: run with no events, no runtrailer, no runheader \n"); |
| 1606 |
|
// |
| 1607 |
|
} else { |
| 1608 |
|
// |
| 1609 |
|
this->HandleRunFragments(true,true,0,(evbefh-1)); |
| 1610 |
|
// |
| 1611 |
|
}; |
| 1612 |
}; |
}; |
| 1613 |
// |
// |
| 1614 |
// |
// |
| 1615 |
if ( (ptht - pth) > 1 ){ |
if ( (ptht - pth) > 1 ){ |
| 1616 |
// |
// |
| 1617 |
if ( IsDebug() ) printf(" Missing runtrailers! \n"); |
if ( IsDebug() ) printf(" Missing runtrailers! \n"); |
| 1618 |
if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %u %u %u \n",ptht,pth,ptt); |
| 1619 |
// is not the consecutive header |
// is not the consecutive header |
| 1620 |
while ( pth != ptht ){ |
while ( pth != ptht ){ |
| 1621 |
// |
// |
| 1636 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1637 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1638 |
// |
// |
| 1639 |
if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader %u %u %u \n",ptht,pth,ptt); |
| 1640 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1641 |
// |
// |
| 1642 |
this->HandleMissingHoT(false,true,evbefh,evbeft-1); |
this->HandleMissingHoT(false,true,evbefh,evbeft-1); |
| 1643 |
// |
// |
| 1646 |
} else if ( !(ptht - pth) ){ |
} else if ( !(ptht - pth) ){ |
| 1647 |
// |
// |
| 1648 |
if ( IsDebug() ) printf(" Missing runheader! \n"); |
if ( IsDebug() ) printf(" Missing runheader! \n"); |
| 1649 |
if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %u %u %u \n",ptht,pth,ptt); |
| 1650 |
if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n"); |
if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n"); |
| 1651 |
// |
// |
| 1652 |
} else { |
} else { |
| 1667 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1668 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
| 1669 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1670 |
if ( IsDebug() ) printf(" Piece of run at the end of file %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" Piece of run at the end of file %u %u %u \n",pkth,obth,obtt); |
| 1671 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u %u %u \n",ptht,pth,ptt); |
| 1672 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i \n",rhev); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u \n",rhev); |
| 1673 |
// |
// |
| 1674 |
this->HandleRunFragments(false,true,evbefh,upperentry); |
this->HandleRunFragments(false,true,evbefh,upperentry); |
| 1675 |
} else { |
} else { |
| 2069 |
ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT()); |
ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT()); |
| 2070 |
UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT()); |
UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT()); |
| 2071 |
ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT()); |
ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT()); |
| 2072 |
if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev); |
if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev); |
| 2073 |
TTree *T= 0; |
TTree *T= 0; |
| 2074 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
| 2075 |
if ( !T || T->IsZombie() ) throw -16; |
if ( !T || T->IsZombie() ) throw -16; |
| 2083 |
bobt = OBT(ph->GetOrbitalTime()); |
bobt = OBT(ph->GetOrbitalTime()); |
| 2084 |
firstev++; |
firstev++; |
| 2085 |
}; |
}; |
| 2086 |
if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev); |
if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev); |
| 2087 |
// |
// |
| 2088 |
glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER()); |
glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER()); |
| 2089 |
glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER()); |
glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER()); |
| 2266 |
ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT()); |
ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT()); |
| 2267 |
UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT()); |
UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT()); |
| 2268 |
ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT()); |
ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT()); |
| 2269 |
if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
| 2270 |
TTree *T= 0; |
TTree *T= 0; |
| 2271 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
| 2272 |
if ( !T || T->IsZombie() ) throw -16; |
if ( !T || T->IsZombie() ) throw -16; |
| 2280 |
aobt = OBT(ph->GetOrbitalTime()); |
aobt = OBT(ph->GetOrbitalTime()); |
| 2281 |
lastev--; |
lastev--; |
| 2282 |
}; |
}; |
| 2283 |
if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
| 2284 |
// |
// |
| 2285 |
glrun->SetEV_TO(lastev); |
glrun->SetEV_TO(lastev); |
| 2286 |
glrun->SetNEVENTS(lastev-firstev+1); |
glrun->SetNEVENTS(lastev-firstev+1); |
| 2887 |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
| 2888 |
if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){ |
if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){ |
| 2889 |
// |
// |
| 2890 |
if ( IsDebug() ) printf(" Calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt); |
if ( IsDebug() ) printf(" Calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt); |
| 2891 |
// |
// |
| 2892 |
// check if the calibration has already been inserted |
// check if the calibration has already been inserted |
| 2893 |
// |
// |
| 2932 |
// |
// |
| 2933 |
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); |
| 2934 |
if ( fromtime < 1150871000 ){ //1150866904 |
if ( fromtime < 1150871000 ){ //1150866904 |
| 2935 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime); |
| 2936 |
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 |
| 2937 |
}; |
}; |
| 2938 |
// |
// |
| 2994 |
// |
// |
| 2995 |
} else { |
} else { |
| 2996 |
// |
// |
| 2997 |
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 %u obt %u pkt %u \n",section,fromtime,obt,pkt); |
| 2998 |
// |
// |
| 2999 |
}; |
}; |
| 3000 |
// |
// |
| 3212 |
// |
// |
| 3213 |
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) ){ |
| 3214 |
// |
// |
| 3215 |
if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1); |
if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1); |
| 3216 |
// |
// |
| 3217 |
valid = ValidateTrkCalib( caltrk1, eh1 ); |
valid = ValidateTrkCalib( caltrk1, eh1 ); |
| 3218 |
if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl; |
if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl; |
| 3255 |
// |
// |
| 3256 |
}; |
}; |
| 3257 |
// |
// |
| 3258 |
if ( IsDebug() ) printf(" Found trk calibration2 at obt %i pkt %i t2 is %i \n",obt2,pkt2,t2); |
if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2); |
| 3259 |
// |
// |
| 3260 |
// The calibration is good |
// The calibration is good |
| 3261 |
// |
// |
| 3262 |
if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){ |
if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){ |
| 3263 |
// |
// |
| 3264 |
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 %u pkt %u t2 is %u is good \n",obt2,pkt2,t2); |
| 3265 |
// |
// |
| 3266 |
UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 ); |
UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 ); |
| 3267 |
if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl; |
if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl; |
| 3275 |
// |
// |
| 3276 |
if ( t2 != pret2+1 ){ |
if ( t2 != pret2+1 ){ |
| 3277 |
// |
// |
| 3278 |
if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i pret2 is %i \n",obt2,pkt2,t2,pret2); |
if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u pret2 is %u \n",obt2,pkt2,t2,pret2); |
| 3279 |
// |
// |
| 3280 |
while ( t2 > pret2+1 ){ |
while ( t2 > pret2+1 ){ |
| 3281 |
// |
// |
| 3304 |
// |
// |
| 3305 |
// Check for missing calibtrk2 |
// Check for missing calibtrk2 |
| 3306 |
// |
// |
| 3307 |
if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2); |
if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2); |
| 3308 |
t2 = pret2; |
t2 = pret2; |
| 3309 |
// |
// |
| 3310 |
// handle missing calib2 |
// handle missing calib2 |
| 3318 |
// |
// |
| 3319 |
} else { |
} else { |
| 3320 |
// |
// |
| 3321 |
if ( IsDebug() ) printf(" Repetead trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1); |
if ( IsDebug() ) printf(" Repetead trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1); |
| 3322 |
// |
// |
| 3323 |
}; |
}; |
| 3324 |
// |
// |
| 3343 |
valid = 0; |
valid = 0; |
| 3344 |
if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){ |
if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){ |
| 3345 |
// |
// |
| 3346 |
if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2); |
if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2); |
| 3347 |
// |
// |
| 3348 |
this->HandleTRK_CALIB(false,true); |
this->HandleTRK_CALIB(false,true); |
| 3349 |
// |
// |
| 3397 |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
| 3398 |
if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){ |
if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){ |
| 3399 |
// |
// |
| 3400 |
if ( IsDebug() ) printf(" S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt); |
if ( IsDebug() ) printf(" S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt); |
| 3401 |
// |
// |
| 3402 |
// check if the calibration has already been inserted |
// check if the calibration has already been inserted |
| 3403 |
// |
// |
| 3440 |
// |
// |
| 3441 |
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"); |
| 3442 |
if ( fromtime < 1150871000 ){ |
if ( fromtime < 1150871000 ){ |
| 3443 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime); |
| 3444 |
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 |
| 3445 |
}; |
}; |
| 3446 |
// |
// |
| 3499 |
// |
// |
| 3500 |
} else { |
} else { |
| 3501 |
// |
// |
| 3502 |
if ( IsDebug() ) printf(" Repetead S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt); |
if ( IsDebug() ) printf(" Repeated S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt); |
| 3503 |
// |
// |
| 3504 |
}; |
}; |
| 3505 |
// |
// |
| 3710 |
row = result->Next(); |
row = result->Next(); |
| 3711 |
t_stop = (UInt_t)atoll(row->GetField(4)); |
t_stop = (UInt_t)atoll(row->GetField(4)); |
| 3712 |
}; |
}; |
| 3713 |
if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start); |
if ( IsDebug() ) printf("Validation interval: from time %u - to time %u \n\n",t_stop,t_start); |
| 3714 |
// -------------------------------------------------------------- |
// -------------------------------------------------------------- |
| 3715 |
// now retrieves runs to be validated |
// now retrieves runs to be validated |
| 3716 |
// -------------------------------------------------------------- |
// -------------------------------------------------------------- |
| 3881 |
}else printf("ValidateRuns ***WARNING*** : run sequence exceed assumed size (%i) \n",nseq_max); |
}else printf("ValidateRuns ***WARNING*** : run sequence exceed assumed size (%i) \n",nseq_max); |
| 3882 |
}; |
}; |
| 3883 |
|
|
| 3884 |
if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID); |
if ( IsDebug() ) printf("%i Run %u \n",nrow,this_run->ID); |
| 3885 |
nrow++; |
nrow++; |
| 3886 |
|
|
| 3887 |
}; |
}; |
| 3910 |
// which should be equal to the time between ascending-nodes. |
// which should be equal to the time between ascending-nodes. |
| 3911 |
//============================================================== |
//============================================================== |
| 3912 |
if ( t2 - trkcalib->FROM_TIME > 5700) { |
if ( t2 - trkcalib->FROM_TIME > 5700) { |
| 3913 |
if ( IsDebug() )printf("Long time between calib and run start %i :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME); |
if ( IsDebug() )printf("Long time between calib and run start %u :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME); |
| 3914 |
//============================================================== |
//============================================================== |
| 3915 |
// there might be a missing calibration, due to: |
// there might be a missing calibration, due to: |
| 3916 |
// - MM full |
// - MM full |
| 3931 |
// it is enough to say that there are no missing calibrations |
// it is enough to say that there are no missing calibrations |
| 3932 |
//============================================================== |
//============================================================== |
| 3933 |
// the long time interval bewteen runs might be due to download |
// the long time interval bewteen runs might be due to download |
| 3934 |
if ( IsDebug() )printf("Short time between calib and run start %i :-) ==> OK! \n",t2 - trkcalib->FROM_TIME); |
if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME); |
| 3935 |
return(false); |
return(false); |
| 3936 |
|
|
| 3937 |
}; |
}; |
| 4247 |
}; |
}; |
| 4248 |
}; |
}; |
| 4249 |
}; |
}; |
| 4250 |
if ( IsDebug() ) printf(" %i runs have been unvalidated \n",unv); |
if ( IsDebug() ) printf(" %u runs have been unvalidated \n",unv); |
| 4251 |
}; |
}; |
| 4252 |
|
|
| 4253 |
/** |
/** |
| 4278 |
if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n"); |
if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n"); |
| 4279 |
while ( Row ){ |
while ( Row ){ |
| 4280 |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
| 4281 |
if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0))); |
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4282 |
drun++; |
drun++; |
| 4283 |
Row = pResult->Next(); |
Row = pResult->Next(); |
| 4284 |
}; |
}; |
| 4299 |
} else { |
} else { |
| 4300 |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n"); |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n"); |
| 4301 |
while ( Row ){ |
while ( Row ){ |
| 4302 |
if ( IsDebug() ) printf(" restore run %i \n",(UInt_t)atoll(Row->GetField(1))); |
if ( IsDebug() ) printf(" restore run %u \n",(UInt_t)atoll(Row->GetField(1))); |
| 4303 |
delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS"); |
delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS"); |
| 4304 |
if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(1))); |
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(1))); |
| 4305 |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN"); |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN"); |
| 4306 |
if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){ |
if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){ |
| 4307 |
if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0))); |
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4308 |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
| 4309 |
}; |
}; |
| 4310 |
drun++; |
drun++; |
| 4330 |
} else { |
} else { |
| 4331 |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n"); |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n"); |
| 4332 |
while ( Row ){ |
while ( Row ){ |
| 4333 |
if ( IsDebug() ) printf(" del run idtrash %i \n",(UInt_t)atoll(Row->GetField(0))); |
if ( IsDebug() ) printf(" del run idtrash %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4334 |
myquery.str(""); |
myquery.str(""); |
| 4335 |
myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";"; |
myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";"; |
| 4336 |
conn->Query(myquery.str().c_str()); |
conn->Query(myquery.str().c_str()); |
| 4339 |
}; |
}; |
| 4340 |
}; |
}; |
| 4341 |
// |
// |
| 4342 |
if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_TRASH table \n",drun); |
if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_TRASH table \n",drun); |
| 4343 |
// |
// |
| 4344 |
// |
// |
| 4345 |
// |
// |
| 4357 |
} else { |
} else { |
| 4358 |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n"); |
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n"); |
| 4359 |
while ( Row ){ |
while ( Row ){ |
| 4360 |
if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0))); |
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4361 |
myquery.str(""); |
myquery.str(""); |
| 4362 |
myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";"; |
myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";"; |
| 4363 |
conn->Query(myquery.str().c_str()); |
conn->Query(myquery.str().c_str()); |
| 4366 |
}; |
}; |
| 4367 |
}; |
}; |
| 4368 |
// |
// |
| 4369 |
if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_FRAGMENTS table \n",drun); |
if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_FRAGMENTS table \n",drun); |
| 4370 |
// |
// |
| 4371 |
// |
// |
| 4372 |
// |
// |