| 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, Bool_t gpamela){ |
| 60 |
// |
// |
| 61 |
|
STATIC=false; |
| 62 |
|
if ( staticp ) STATIC=true; |
| 63 |
// |
// |
| 64 |
SetConnection(host,user,password); |
SetConnection(host,user,password); |
| 65 |
// |
// |
| 68 |
glrun = new GL_RUN(); |
glrun = new GL_RUN(); |
| 69 |
// |
// |
| 70 |
if ( !boot ) SetNOBOOT(false); |
if ( !boot ) SetNOBOOT(false); |
| 71 |
SetBOOTnumber(boot); |
SetTsync(tsync,gpamela); |
| 72 |
SetTsync(tsync); |
SetBOOTnumber(boot,gpamela); |
| 73 |
SetObt0(obt0); |
SetObt0(obt0); |
| 74 |
// |
// |
| 75 |
SetTLEPath(tlefilename); |
SetTLEPath(tlefilename); |
| 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 { |
| 87 |
|
this->SetRootName(""); |
| 88 |
}; |
}; |
| 89 |
// |
// |
| 90 |
this->SetID_RAW(0); |
this->SetID_RAW(0); |
| 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 |
| 171 |
|
* |
| 172 |
|
*/ |
| 173 |
|
void PamelaDBOperations::SetNoFrag(Bool_t nf){ |
| 174 |
|
NOFRAG = nf; |
| 175 |
|
}; |
| 176 |
|
|
| 177 |
|
/** |
| 178 |
* Store the BOOT number of the RAW file. |
* Store the BOOT number of the RAW file. |
| 179 |
* @param boot BOOT number of the RAW file |
* @param boot BOOT number of the RAW file |
| 180 |
*/ |
*/ |
| 181 |
void PamelaDBOperations::SetBOOTnumber(UInt_t boot){ |
void PamelaDBOperations::SetBOOTnumber(UInt_t boot){ |
| 182 |
|
this->SetBOOTnumber(boot,false); |
| 183 |
|
}; |
| 184 |
|
|
| 185 |
|
/** |
| 186 |
|
* Store the BOOT number of the RAW file. |
| 187 |
|
* @param boot BOOT number of the RAW file |
| 188 |
|
*/ |
| 189 |
|
void PamelaDBOperations::SetBOOTnumber(UInt_t boot, Bool_t gpamela){ |
| 190 |
BOOTNO=boot; |
BOOTNO=boot; |
| 191 |
|
if ( gpamela ){ |
| 192 |
|
stringstream oss; |
| 193 |
|
TSQLResult *result = 0; |
| 194 |
|
TSQLRow *row = 0; |
| 195 |
|
if ( !boot ){ |
| 196 |
|
// |
| 197 |
|
BOOTNO = 1; |
| 198 |
|
// |
| 199 |
|
// look in the DB for the last timesync and the last run |
| 200 |
|
// |
| 201 |
|
oss.str(""); |
| 202 |
|
oss << "SELECT BOOT_NUMBER FROM GL_RUN order by RUNHEADER_TIME desc limit 1;"; |
| 203 |
|
result = conn->Query(oss.str().c_str()); |
| 204 |
|
if ( result ){ |
| 205 |
|
row = result->Next(); |
| 206 |
|
if ( row ){ |
| 207 |
|
BOOTNO = (UInt_t)atoll(row->GetField(0)) + 1; |
| 208 |
|
}; |
| 209 |
|
}; |
| 210 |
|
}; |
| 211 |
|
}; |
| 212 |
}; |
}; |
| 213 |
|
|
| 214 |
/** |
/** |
| 216 |
* @param boot time sync |
* @param boot time sync |
| 217 |
*/ |
*/ |
| 218 |
void PamelaDBOperations::SetTsync(UInt_t ts){ |
void PamelaDBOperations::SetTsync(UInt_t ts){ |
| 219 |
|
this->SetTsync(ts,false); |
| 220 |
|
}; |
| 221 |
|
|
| 222 |
|
/** |
| 223 |
|
* Store the time sync of the RAW file. |
| 224 |
|
* @param boot time sync |
| 225 |
|
*/ |
| 226 |
|
void PamelaDBOperations::SetTsync(UInt_t ts, Bool_t gpamela){ |
| 227 |
|
// |
| 228 |
|
// if not gpamela or given tsync file set ts |
| 229 |
|
// |
| 230 |
tsync=ts; |
tsync=ts; |
| 231 |
|
if ( gpamela ){ |
| 232 |
|
stringstream oss; |
| 233 |
|
TSQLResult *result = 0; |
| 234 |
|
TSQLRow *row = 0; |
| 235 |
|
TSQLResult *result2 = 0; |
| 236 |
|
TSQLRow *row2 = 0; |
| 237 |
|
if ( !ts ){ |
| 238 |
|
// |
| 239 |
|
tsync = 1; |
| 240 |
|
// |
| 241 |
|
// look in the DB for the last timesync and the last run |
| 242 |
|
// |
| 243 |
|
oss.str(""); |
| 244 |
|
oss << "SELECT TIMESYNC FROM GL_TIMESYNC order by TIMESYNC desc limit 1;"; |
| 245 |
|
result = conn->Query(oss.str().c_str()); |
| 246 |
|
if ( result ){ |
| 247 |
|
row = result->Next(); |
| 248 |
|
if ( row ){ |
| 249 |
|
tsync = (UInt_t)atoll(row->GetField(0)) + 1; |
| 250 |
|
oss.str(""); |
| 251 |
|
oss << "SELECT (RUNTRAILER_TIME-RUNHEADER_TIME) FROM GL_RUN order by RUNHEADER_TIME desc limit 1;"; |
| 252 |
|
result2 = conn->Query(oss.str().c_str()); |
| 253 |
|
if ( result2 ){ |
| 254 |
|
row2 = result2->Next(); |
| 255 |
|
if ( row2 ){ |
| 256 |
|
tsync += (UInt_t)atoll(row2->GetField(0)); |
| 257 |
|
}; |
| 258 |
|
} |
| 259 |
|
}; |
| 260 |
|
}; |
| 261 |
|
}; |
| 262 |
|
}; |
| 263 |
}; |
}; |
| 264 |
|
|
| 265 |
/** |
/** |
| 289 |
/** |
/** |
| 290 |
* Store the downlink orbit number from filename. |
* Store the downlink orbit number from filename. |
| 291 |
*/ |
*/ |
| 292 |
void PamelaDBOperations::SetOrbitNo(){ |
void PamelaDBOperations::SetOrbitNo(UInt_t dwinput){ |
| 293 |
dworbit = 0; |
dworbit = 0; |
| 294 |
|
// |
| 295 |
|
if ( dwinput ){ |
| 296 |
|
dworbit = dwinput; |
| 297 |
|
if ( IsDebug() ) printf(" Downlink orbit given by hand: %i \n",dworbit); |
| 298 |
|
return; |
| 299 |
|
}; |
| 300 |
|
// |
| 301 |
TString name = this->GetRootFile(); |
TString name = this->GetRootFile(); |
| 302 |
Int_t nlength = name.Length(); |
Int_t nlength = name.Length(); |
| 303 |
if ( nlength < 5 ) return; |
if ( nlength < 5 ) return; |
| 335 |
tlefilename = str; |
tlefilename = str; |
| 336 |
}; |
}; |
| 337 |
|
|
| 338 |
|
TString PamelaDBOperations::GetRawPath(){ |
| 339 |
|
if ( STATIC ){ |
| 340 |
|
return((TString)gSystem->DirName(filerawname.Data())+'/'); |
| 341 |
|
} else { |
| 342 |
|
return((TString)gSystem->ExpandPathName("$PAM_RAW")+'/'); |
| 343 |
|
}; |
| 344 |
|
}; |
| 345 |
|
|
| 346 |
|
TString PamelaDBOperations::GetRootPath(){ |
| 347 |
|
if ( STATIC ){ |
| 348 |
|
return((TString)gSystem->DirName(filerootname.Data())+'/'); |
| 349 |
|
} else { |
| 350 |
|
return((TString)gSystem->ExpandPathName("$PAM_L0")+'/'); |
| 351 |
|
}; |
| 352 |
|
}; |
| 353 |
|
|
| 354 |
/** |
/** |
| 355 |
* Store the olderthan variable |
* Store the olderthan variable |
| 356 |
* @param olderthan |
* @param olderthan |
| 367 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 368 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 369 |
oss.str(""); |
oss.str(""); |
| 370 |
oss << "SELECT ID FROM GL_RAW WHERE " |
if ( STATIC ){ |
| 371 |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
oss << "SELECT ID FROM GL_RAW WHERE " |
| 372 |
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
| 373 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 374 |
|
} else { |
| 375 |
|
oss << "SELECT ID FROM GL_RAW WHERE " |
| 376 |
|
<< " PATH = '$PAM_RAW' AND " |
| 377 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 378 |
|
} |
| 379 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 380 |
if ( result == NULL ) throw -4; |
if ( result == NULL ) throw -4; |
| 381 |
row = result->Next(); |
row = result->Next(); |
| 406 |
UInt_t nevent = 0; |
UInt_t nevent = 0; |
| 407 |
UInt_t pktlast = 0; |
UInt_t pktlast = 0; |
| 408 |
UInt_t obtlast = 0; |
UInt_t obtlast = 0; |
| 409 |
UInt_t t_pktlast = 0; |
Long64_t t_pktlast = 0LL; |
| 410 |
UInt_t t_obtlast = 0; |
// UInt_t t_obtlast = 0; |
| 411 |
UInt_t upperpkt2 = 0; |
Long64_t t_obtlast = 0LL; |
| 412 |
ULong64_t upperobt2 = 0; |
Long64_t upperpkt2 = 0LL; |
| 413 |
|
Long64_t upperobt2 = 0LL; |
| 414 |
UInt_t zomp = 0; |
UInt_t zomp = 0; |
| 415 |
UInt_t jump = 50000; // was 5000 |
UInt_t jump = 50000; // was 5000 |
| 416 |
EventCounter *code=0; |
EventCounter *code=0; |
| 417 |
// |
// |
| 418 |
UInt_t deltapkt = 5000; |
Long64_t deltapkt = 5000LL; |
| 419 |
ULong64_t deltaobt = 50000; |
Long64_t deltaobt = 50000LL; |
| 420 |
// |
// |
| 421 |
// pcksList packetsNames; |
pcksList packetsNames; |
| 422 |
// pcksList::iterator Iter; |
pcksList::iterator Iter; |
| 423 |
// getPacketsNames(packetsNames); |
getPacketsNames(packetsNames); |
| 424 |
// |
// |
| 425 |
pktfirst = 0; |
pktfirst = 0; |
| 426 |
obtfirst = 0; |
obtfirst = 0; |
| 427 |
|
ppktfirst = 0; |
| 428 |
|
pobtfirst = 0; |
| 429 |
// |
// |
| 430 |
TTree *T = 0; |
TTree *T = 0; |
| 431 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
| 439 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
| 440 |
pktfirst = ph->GetCounter(); |
pktfirst = ph->GetCounter(); |
| 441 |
obtfirst = ph->GetOrbitalTime(); |
obtfirst = ph->GetOrbitalTime(); |
| 442 |
|
ppktfirst = pktfirst; |
| 443 |
|
pobtfirst = obtfirst; |
| 444 |
// |
// |
| 445 |
// code = eh->GetCounter(); |
code = eh->GetCounter(); |
| 446 |
// UInt_t en = 0; |
UInt_t en = 0; |
| 447 |
// for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
| 448 |
// en = code->Get(GetPacketType(*Iter)); |
en = code->Get(GetPacketType(*Iter)); |
| 449 |
// 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) ){ |
| 450 |
//}; |
if ( IsDebug() ) printf(" We have a calibration before the first physic packet: packet type is %s, entries: %i \n",*Iter,en); |
| 451 |
|
// |
| 452 |
|
TTree *TC = 0; |
| 453 |
|
TC = (TTree*)file->Get("CalibCalPed"); |
| 454 |
|
if ( !TC || TC->IsZombie() ) throw -16; |
| 455 |
|
EventHeader *ehc = 0; |
| 456 |
|
PscuHeader *phc = 0; |
| 457 |
|
TC->SetBranchAddress("Header", &ehc); |
| 458 |
|
TC->GetEntry(0); |
| 459 |
|
phc = ehc->GetPscuHeader(); |
| 460 |
|
pktfirst = phc->GetCounter(); |
| 461 |
|
obtfirst = phc->GetOrbitalTime(); |
| 462 |
|
// |
| 463 |
|
}; |
| 464 |
|
}; |
| 465 |
// |
// |
| 466 |
T->GetEntry(nevent-1); |
T->GetEntry(nevent-1); |
| 467 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
| 472 |
upperobt = OBT(obtlast); |
upperobt = OBT(obtlast); |
| 473 |
upperentry = nevent-1; |
upperentry = nevent-1; |
| 474 |
// |
// |
| 475 |
if ( IsDebug() ) printf(" First entries are: OBT %llu pkt_num %i \n",obtfirst,pktfirst); |
if ( IsDebug() ) printf(" First entries are: OBT %u pkt_num %u \n",obtfirst,pktfirst); |
| 476 |
// |
// |
| 477 |
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); |
| 478 |
// |
// |
| 479 |
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); |
| 480 |
// |
// |
| 481 |
if ( !nevent ) return(2); |
if ( !nevent ) return(64); |
| 482 |
// |
// |
| 483 |
if ( nevent < 2 ) return(4); |
if ( nevent < 2 ) return(128); |
| 484 |
|
if ( nevent < jump ) jump = 1; |
| 485 |
|
// if ( nevent < jump ) jump = int(nevent/10); |
| 486 |
|
// if ( !jump ) jump = 1; |
| 487 |
// |
// |
| 488 |
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 ){ |
| 489 |
// |
// |
| 490 |
|
if ( IsDebug() ) printf(" starting jump %i \n",jump); |
| 491 |
// go back |
// go back |
| 492 |
zomp = nevent - 2; |
zomp = nevent - 2; |
| 493 |
// |
// |
| 509 |
upperpkt2 = PKT(ph->GetCounter()); |
upperpkt2 = PKT(ph->GetCounter()); |
| 510 |
upperobt2 = OBT(ph->GetOrbitalTime()); |
upperobt2 = OBT(ph->GetOrbitalTime()); |
| 511 |
// |
// |
| 512 |
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) ){ |
| 513 |
|
if ( IsDebug() ) printf(" .-. upperpkt2 %lld upperobt2 %lld \n",upperpkt2,upperobt2); |
| 514 |
|
if ( IsDebug() ) printf(" .-. upperpkt %lld t_pktlast %lld upperobt %lld t_obtlast %lld \n",upperpkt,t_pktlast,upperobt,t_obtlast); |
| 515 |
|
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); |
| 516 |
|
throw -13; |
| 517 |
|
}; |
| 518 |
// |
// |
| 519 |
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 ){ |
| 520 |
zomp = i + jump + 1; |
zomp = i + jump + 1; |
| 521 |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
if ( zomp > nevent-2 ) zomp = nevent - 2; |
| 522 |
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 %u upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i); |
| 523 |
break; |
break; |
| 524 |
}; |
}; |
| 525 |
// |
// |
| 554 |
// |
// |
| 555 |
rhev = rh->GetEntries(); |
rhev = rh->GetEntries(); |
| 556 |
rtev = rt->GetEntries(); |
rtev = rt->GetEntries(); |
| 557 |
UInt_t sobtt = 0; |
Long64_t sobtt = 0LL; |
| 558 |
UInt_t sobth = 0; |
Long64_t sobth = 0LL; |
| 559 |
UInt_t spktt = 0; |
Long64_t spktt = 0LL; |
| 560 |
UInt_t spkth = 0; |
Long64_t spkth = 0LL; |
| 561 |
UInt_t pktt = 0; |
Long64_t pktt = 0LL; |
| 562 |
ULong64_t obtt = 0; |
Long64_t obtt = 0LL; |
| 563 |
UInt_t pkth = 0; |
Long64_t pkth = 0LL; |
| 564 |
ULong64_t obth = 0; |
Long64_t obth = 0LL; |
| 565 |
// |
// |
| 566 |
T->GetEntry(upperentry); |
if ( rhev || rtev ){ |
| 567 |
code = eh->GetCounter(); |
|
| 568 |
Int_t lasttrail = code->Get(pctp->RunTrailer); |
T->GetEntry(upperentry); |
| 569 |
Int_t lasthead = code->Get(pctp->RunHeader); |
code = eh->GetCounter(); |
| 570 |
if ( lasttrail < rtev ){ |
Int_t lasttrail = code->Get(pctp->RunTrailer); |
| 571 |
rt->GetEntry(lasttrail); |
Int_t lasthead = code->Get(pctp->RunHeader); |
| 572 |
pht = eht->GetPscuHeader(); |
if ( lasttrail < rtev ){ |
| 573 |
pktt = PKT(pht->GetCounter()); |
rt->GetEntry(lasttrail); |
| 574 |
obtt = OBT(pht->GetOrbitalTime()); |
pht = eht->GetPscuHeader(); |
| 575 |
}; |
pktt = PKT(pht->GetCounter()); |
| 576 |
// |
obtt = OBT(pht->GetOrbitalTime()); |
| 577 |
if ( lasthead < rhev ){ |
}; |
| 578 |
rh->GetEntry(lasthead); |
// |
| 579 |
phh = ehh->GetPscuHeader(); |
if ( lasthead < rhev ){ |
| 580 |
pkth = PKT(phh->GetCounter()); |
rh->GetEntry(lasthead); |
| 581 |
obth = OBT(phh->GetOrbitalTime()); |
phh = ehh->GetPscuHeader(); |
| 582 |
}; |
pkth = PKT(phh->GetCounter()); |
| 583 |
// |
obth = OBT(phh->GetOrbitalTime()); |
| 584 |
if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
}; |
| 585 |
if ( pkth > upperpkt && obth > upperobt ){ |
// |
| 586 |
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); |
| 587 |
upperpkt = pkth; |
if ( pkth > upperpkt && obth > upperobt ){ |
| 588 |
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); |
| 589 |
rhev = lasthead+1; |
upperpkt = pkth; |
| 590 |
} else { |
upperobt = obth; |
| 591 |
rhev = lasthead; |
rhev = lasthead+1; |
| 592 |
}; |
} else { |
| 593 |
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
rhev = lasthead; |
| 594 |
// |
}; |
| 595 |
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); |
| 596 |
if ( pktt > upperpkt && obtt > upperobt ){ |
// |
| 597 |
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); |
| 598 |
upperpkt = pktt; |
if ( pktt > upperpkt && obtt > upperobt ){ |
| 599 |
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); |
| 600 |
rtev = lasttrail+1; |
upperpkt = pktt; |
| 601 |
} else { |
upperobt = obtt; |
| 602 |
rtev = lasttrail; |
rtev = lasttrail+1; |
| 603 |
}; |
} else { |
| 604 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
rtev = lasttrail; |
| 605 |
// goto kikko; |
}; |
| 606 |
// |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
| 607 |
// |
// goto kikko; |
| 608 |
// Check if runtrailer/runheader are within lower limits |
// |
| 609 |
// |
// |
| 610 |
// |
// Check if runtrailer/runheader are within lower limits |
| 611 |
pkth = 0; |
// |
| 612 |
obth = 0; |
// |
| 613 |
spkth = 0; |
pkth = 0LL; |
| 614 |
sobth = 0; |
obth = 0LL; |
| 615 |
for (Int_t k=0; k<rhev; k++){ |
spkth = 0LL; |
| 616 |
if ( k > 0 ){ |
sobth = 0LL; |
| 617 |
spkth = pkth; |
for (Int_t k=0; k<rhev; k++){ |
| 618 |
sobth = obth; |
if ( k > 0 ){ |
| 619 |
}; |
spkth = pkth; |
| 620 |
rh->GetEntry(k); |
sobth = obth; |
| 621 |
phh = ehh->GetPscuHeader(); |
}; |
| 622 |
pkth = PKT(phh->GetCounter()); |
rh->GetEntry(k); |
| 623 |
obth = OBT(phh->GetOrbitalTime()); |
phh = ehh->GetPscuHeader(); |
| 624 |
// |
pkth = PKT(phh->GetCounter()); |
| 625 |
// 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()); |
| 626 |
// |
// |
| 627 |
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); |
| 628 |
if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i \n",rhev); |
// |
| 629 |
// |
if ( pkth < spkth && obth < sobth ){ |
| 630 |
rhev = k-1; |
if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i \n",rhev); |
| 631 |
rh->GetEntry(rhev); |
// |
| 632 |
pkth = spkth; |
rhev = k-1; |
| 633 |
obth = sobth; |
rh->GetEntry(rhev); |
| 634 |
// |
pkth = spkth; |
| 635 |
UInt_t evbefh = 0; |
obth = sobth; |
| 636 |
code = ehh->GetCounter(); |
// |
| 637 |
evbefh = code->Get(pctp->Physics); |
UInt_t evbefh = 0; |
| 638 |
if ( evbefh >= 0 ){ |
code = ehh->GetCounter(); |
| 639 |
T->GetEntry(evbefh); |
evbefh = code->Get(pctp->Physics); |
| 640 |
ph = eh->GetPscuHeader(); |
if ( evbefh >= 0 ){ |
| 641 |
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); |
|
| 642 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
| 643 |
upperpkt = PKT(ph->GetCounter()); |
t_pktlast = PKT(ph->GetCounter()); |
| 644 |
upperobt = OBT(ph->GetOrbitalTime()); |
t_obtlast = OBT(ph->GetOrbitalTime()); |
| 645 |
upperentry = evbefh-1; |
if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump |
| 646 |
}; |
upperpkt = pkth; |
| 647 |
|
upperobt = obth; |
| 648 |
|
upperentry = evbefh-1; |
| 649 |
|
} else { |
| 650 |
|
while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){ |
| 651 |
|
evbefh++; |
| 652 |
|
T->GetEntry(evbefh); |
| 653 |
|
ph = eh->GetPscuHeader(); |
| 654 |
|
t_pktlast = PKT(ph->GetCounter()); |
| 655 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
| 656 |
|
}; |
| 657 |
|
T->GetEntry(evbefh-1); |
| 658 |
|
ph = eh->GetPscuHeader(); |
| 659 |
|
upperpkt = PKT(ph->GetCounter()); |
| 660 |
|
upperobt = OBT(ph->GetOrbitalTime()); |
| 661 |
|
upperentry = evbefh-1; |
| 662 |
|
}; |
| 663 |
|
}; |
| 664 |
|
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
| 665 |
|
goto kikko0; |
| 666 |
}; |
}; |
|
if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt); |
|
|
goto kikko0; |
|
| 667 |
}; |
}; |
| 668 |
}; |
kikko0: |
| 669 |
kikko0: |
// |
|
// |
|
|
// |
|
|
// |
|
|
pktt = 0; |
|
|
obtt = 0; |
|
|
spktt = 0; |
|
|
sobtt = 0; |
|
|
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()); |
|
| 670 |
// |
// |
|
// if ( IsDebug() ) printf(" k %i rtev beforev %i pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt); |
|
| 671 |
// |
// |
| 672 |
if ( pktt < spktt && obtt < sobtt ){ |
pktt = 0LL; |
| 673 |
if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev); |
obtt = 0LL; |
| 674 |
|
spktt = 0LL; |
| 675 |
|
sobtt = 0LL; |
| 676 |
|
for (Int_t k=0; k<rtev; k++){ |
| 677 |
|
if ( k > 0 ){ |
| 678 |
|
spktt = pktt; |
| 679 |
|
sobtt = obtt; |
| 680 |
|
}; |
| 681 |
|
rt->GetEntry(k); |
| 682 |
|
pht = eht->GetPscuHeader(); |
| 683 |
|
pktt = PKT(pht->GetCounter()); |
| 684 |
|
obtt = OBT(pht->GetOrbitalTime()); |
| 685 |
// |
// |
| 686 |
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); |
|
| 687 |
// |
// |
| 688 |
UInt_t evbeft = 0; |
if ( pktt < spktt && obtt < sobtt ){ |
| 689 |
code = eht->GetCounter(); |
if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev); |
| 690 |
evbeft = code->Get(pctp->Physics); |
// |
| 691 |
if ( evbeft >= 0 ){ |
rtev = k-1; |
| 692 |
T->GetEntry(evbeft); |
rt->GetEntry(rtev); |
| 693 |
ph = eh->GetPscuHeader(); |
pktt = spktt; |
| 694 |
t_pktlast = PKT(ph->GetCounter()); |
obtt = sobtt; |
| 695 |
t_obtlast = OBT(ph->GetOrbitalTime()); |
if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
| 696 |
if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump |
// |
| 697 |
upperpkt = pktt; |
UInt_t evbeft = 0; |
| 698 |
upperobt = obtt; |
code = eht->GetCounter(); |
| 699 |
upperentry = evbeft-1; |
evbeft = code->Get(pctp->Physics); |
| 700 |
} else { |
if ( evbeft >= 0 ){ |
| 701 |
while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){ |
T->GetEntry(evbeft); |
| 702 |
evbeft++; |
ph = eh->GetPscuHeader(); |
| 703 |
T->GetEntry(evbeft); |
t_pktlast = PKT(ph->GetCounter()); |
| 704 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
| 705 |
|
if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump |
| 706 |
|
upperpkt = pktt; |
| 707 |
|
upperobt = obtt; |
| 708 |
|
upperentry = evbeft-1; |
| 709 |
|
} else { |
| 710 |
|
while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){ |
| 711 |
|
evbeft++; |
| 712 |
|
T->GetEntry(evbeft); |
| 713 |
|
ph = eh->GetPscuHeader(); |
| 714 |
|
t_pktlast = PKT(ph->GetCounter()); |
| 715 |
|
t_obtlast = OBT(ph->GetOrbitalTime()); |
| 716 |
|
}; |
| 717 |
|
T->GetEntry(evbeft-1); |
| 718 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
| 719 |
t_pktlast = PKT(ph->GetCounter()); |
upperpkt = PKT(ph->GetCounter()); |
| 720 |
t_obtlast = OBT(ph->GetOrbitalTime()); |
upperobt = OBT(ph->GetOrbitalTime()); |
| 721 |
|
upperentry = evbeft-1; |
| 722 |
}; |
}; |
|
T->GetEntry(evbeft-1); |
|
|
ph = eh->GetPscuHeader(); |
|
|
upperpkt = PKT(ph->GetCounter()); |
|
|
upperobt = OBT(ph->GetOrbitalTime()); |
|
|
upperentry = evbeft-1; |
|
| 723 |
}; |
}; |
| 724 |
|
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
| 725 |
|
goto kikko; |
| 726 |
|
// break; |
| 727 |
|
// |
| 728 |
}; |
}; |
| 729 |
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
// |
|
goto kikko; |
|
|
// break; |
|
|
// |
|
| 730 |
}; |
}; |
| 731 |
// |
// |
| 732 |
}; |
kikko: |
| 733 |
// |
// |
| 734 |
kikko: |
T->GetEntry(upperentry); |
| 735 |
// |
code = eh->GetCounter(); |
| 736 |
T->GetEntry(upperentry); |
lasttrail = code->Get(pctp->RunTrailer); |
| 737 |
code = eh->GetCounter(); |
lasthead = code->Get(pctp->RunHeader); |
| 738 |
lasttrail = code->Get(pctp->RunTrailer); |
if ( lasttrail < rtev ){ |
| 739 |
lasthead = code->Get(pctp->RunHeader); |
rt->GetEntry(lasttrail); |
| 740 |
if ( lasttrail < rtev ){ |
pht = eht->GetPscuHeader(); |
| 741 |
rt->GetEntry(lasttrail); |
pktt = PKT(pht->GetCounter()); |
| 742 |
pht = eht->GetPscuHeader(); |
obtt = OBT(pht->GetOrbitalTime()); |
| 743 |
pktt = PKT(pht->GetCounter()); |
}; |
| 744 |
obtt = OBT(pht->GetOrbitalTime()); |
// |
| 745 |
}; |
if ( lasthead < rhev ){ |
| 746 |
// |
rh->GetEntry(lasthead); |
| 747 |
if ( lasthead < rhev ){ |
phh = ehh->GetPscuHeader(); |
| 748 |
rh->GetEntry(lasthead); |
pkth = PKT(phh->GetCounter()); |
| 749 |
phh = ehh->GetPscuHeader(); |
obth = OBT(phh->GetOrbitalTime()); |
| 750 |
pkth = PKT(phh->GetCounter()); |
}; |
| 751 |
obth = OBT(phh->GetOrbitalTime()); |
// |
| 752 |
}; |
if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
| 753 |
// |
if ( pkth > upperpkt && obth > upperobt ){ |
| 754 |
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); |
| 755 |
if ( pkth > upperpkt && obth > upperobt ){ |
upperpkt = pkth; |
| 756 |
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; |
| 757 |
upperpkt = pkth; |
rhev = lasthead+1; |
| 758 |
upperobt = obth; |
} else { |
| 759 |
rhev = lasthead+1; |
rhev = lasthead; |
| 760 |
} else { |
}; |
| 761 |
rhev = lasthead; |
if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt); |
| 762 |
}; |
// |
| 763 |
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); |
| 764 |
// |
if ( pktt > upperpkt && obtt > upperobt ){ |
| 765 |
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); |
| 766 |
if ( pktt > upperpkt && obtt > upperobt ){ |
upperpkt = pktt; |
| 767 |
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; |
| 768 |
upperpkt = pktt; |
rtev = lasttrail+1; |
| 769 |
upperobt = obtt; |
} else { |
| 770 |
rtev = lasttrail+1; |
rtev = lasttrail; |
| 771 |
} else { |
}; |
| 772 |
rtev = lasttrail; |
if ( IsDebug() ) printf(" rtev after %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt); |
| 773 |
|
// |
| 774 |
}; |
}; |
|
if ( IsDebug() ) printf(" rtev after %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt); |
|
| 775 |
// |
// |
| 776 |
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); |
| 777 |
// |
// |
| 778 |
return(0); |
return(0); |
| 779 |
} |
} |
| 907 |
if( !conn ) throw -1; |
if( !conn ) throw -1; |
| 908 |
bool connect = conn->IsConnected(); |
bool connect = conn->IsConnected(); |
| 909 |
if( !connect ) throw -1; |
if( !connect ) throw -1; |
| 910 |
if ( !dworbit ) throw -27; |
// |
| 911 |
|
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()); |
| 912 |
|
// |
| 913 |
|
if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27; |
| 914 |
|
// |
| 915 |
|
// set DB timezone to UTC |
| 916 |
|
// |
| 917 |
|
stringstream oss; |
| 918 |
|
// |
| 919 |
|
oss.str(""); |
| 920 |
|
oss << "SET time_zone='+0:00';"; |
| 921 |
|
TSQLResult *result = 0; |
| 922 |
|
result = conn->Query(oss.str().c_str()); |
| 923 |
|
if ( !result ) throw -10; |
| 924 |
|
oss.str(""); |
| 925 |
|
oss << "SET wait_timeout=173000;"; |
| 926 |
|
conn->Query(oss.str().c_str()); |
| 927 |
|
// |
| 928 |
}; |
}; |
| 929 |
|
|
| 930 |
/** |
/** |
| 931 |
* Return the correct packet number if we went back to zero |
* Return the correct packet number if we went back to zero |
| 932 |
*/ |
*/ |
| 933 |
UInt_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){ |
| 934 |
// |
// |
| 935 |
// 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)); |
| 936 |
// |
// |
| 937 |
if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2) ) return((pkt_num+16777215)); |
if ( pkt_num < (ppktfirst/2) && ppktfirst > (16777214/2) ){ |
| 938 |
|
if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL); |
| 939 |
|
return((Long64_t)pkt_num+16777215LL); |
| 940 |
|
}; |
| 941 |
// |
// |
| 942 |
if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){ |
if ( pkt_num > ((Long64_t)ppktfirst*2) && pkt_num > (16777214/2) ){ |
| 943 |
if ( (pkt_num-16777215) < 0 ){ |
if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL); |
| 944 |
return((16777215-pkt_num)); |
return((Long64_t)pkt_num-16777215LL); |
|
} else { |
|
|
return((pkt_num-16777215)); |
|
|
}; |
|
| 945 |
}; |
}; |
| 946 |
// |
// |
| 947 |
return(pkt_num); |
if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)pkt_num); |
| 948 |
|
return((Long64_t)pkt_num); |
| 949 |
// |
// |
| 950 |
}; |
}; |
| 951 |
|
|
| 952 |
/** |
/** |
| 953 |
* Return the correct On Board Time if we went back to zero |
* Return the correct On Board Time if we went back to zero |
| 954 |
*/ |
*/ |
| 955 |
ULong64_t PamelaDBOperations::OBT(UInt_t obt){ |
Long64_t PamelaDBOperations::OBT(UInt_t obt){ |
| 956 |
// |
// |
| 957 |
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: 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)); |
| 958 |
// |
// |
| 959 |
if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){ |
| 960 |
if ( (obt-numeric_limits<UInt_t>::max()) < 0 ){ |
if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)(obt+numeric_limits<UInt_t>::max())); |
| 961 |
return((ULong64_t)(numeric_limits<UInt_t>::max()-obt)); |
return((Long64_t)(obt+numeric_limits<UInt_t>::max())); |
| 962 |
} else { |
}; |
| 963 |
return((ULong64_t)(obt-numeric_limits<UInt_t>::max())); |
// |
| 964 |
}; |
if ( obt > ((Long64_t)pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){ |
| 965 |
|
if ( IsDebug() ) printf(" pobtfirst*2 %lld \n",((Long64_t)pobtfirst*2)); |
| 966 |
|
if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
| 967 |
|
return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max()); |
| 968 |
}; |
}; |
| 969 |
// |
// |
| 970 |
return((ULong64_t)obt); |
if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)obt); |
| 971 |
|
return((Long64_t)obt); |
| 972 |
}; |
}; |
| 973 |
|
|
| 974 |
/** |
/** |
| 1083 |
if ( idr ) return(1); |
if ( idr ) return(1); |
| 1084 |
// |
// |
| 1085 |
oss.str(""); |
oss.str(""); |
| 1086 |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
if ( STATIC ){ |
| 1087 |
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
| 1088 |
|
<< this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
| 1089 |
|
} else { |
| 1090 |
|
oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')"; |
| 1091 |
|
}; |
| 1092 |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
| 1093 |
// |
// |
| 1094 |
idr = this->SetID_RAW(); |
idr = this->SetID_RAW(); |
| 1107 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 1108 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 1109 |
UInt_t t0 = 0; |
UInt_t t0 = 0; |
| 1110 |
|
Int_t signal = 0; |
| 1111 |
|
// |
| 1112 |
|
signal = this->SetUpperLimits(); |
| 1113 |
// |
// |
| 1114 |
stringstream oss; |
stringstream oss; |
| 1115 |
// |
// |
| 1116 |
if ( this->GetID_RAW() == 0 ) throw -11; |
if ( this->GetID_RAW() == 0 ) throw -11; |
| 1117 |
// |
// |
| 1118 |
oss.str(""); |
oss.str(""); |
| 1119 |
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='" |
| 1120 |
<< this->GetRawFile().Data() << "';"; |
<< this->GetRawFile().Data() << "';"; |
| 1121 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
| 1122 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1125 |
// |
// |
| 1126 |
if ( !row ){ |
if ( !row ){ |
| 1127 |
oss.str(""); |
oss.str(""); |
| 1128 |
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< " |
| 1129 |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
| 1130 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
| 1131 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1134 |
if ( !row ) throw -10; |
if ( !row ) throw -10; |
| 1135 |
}; |
}; |
| 1136 |
// |
// |
| 1137 |
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); |
| 1138 |
|
t0 = (UInt_t)tu.GetSec(); |
| 1139 |
|
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)); |
| 1140 |
|
// |
| 1141 |
/* |
/* |
| 1142 |
* Verify that the TIMESYNC have been not already processed |
* Verify that the TIMESYNC have been not already processed |
| 1143 |
*/ |
*/ |
| 1152 |
if (result == NULL) throw -10; |
if (result == NULL) throw -10; |
| 1153 |
row = result->Next(); |
row = result->Next(); |
| 1154 |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){ |
| 1155 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
| 1156 |
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; |
| 1157 |
|
// |
| 1158 |
|
tsync = (UInt_t)atoll(row->GetField(2)); |
| 1159 |
|
obt0 = (UInt_t)atoll(row->GetField(1)); |
| 1160 |
|
// |
| 1161 |
return(1); |
return(1); |
| 1162 |
}; |
}; |
| 1163 |
// |
// |
| 1164 |
TTree *T = 0; |
TTree *T = 0; |
|
Int_t signal = 0; |
|
| 1165 |
// |
// |
| 1166 |
UInt_t nevent = 0; |
UInt_t nevent = 0; |
| 1167 |
UInt_t recEntries = 0; |
UInt_t recEntries = 0; |
| 1350 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
| 1351 |
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()); |
| 1352 |
// |
// |
| 1353 |
|
if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0); |
| 1354 |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0; |
| 1355 |
// |
// |
| 1356 |
|
tsync = TSYNC; |
| 1357 |
|
obt0 = OBT; |
| 1358 |
|
// |
| 1359 |
delete result; |
delete result; |
| 1360 |
return(signal); |
return(signal); |
| 1361 |
} |
} |
| 1371 |
UInt_t idtimesync = 0; |
UInt_t idtimesync = 0; |
| 1372 |
// |
// |
| 1373 |
oss.str(""); |
oss.str(""); |
| 1374 |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
if ( STATIC ){ |
| 1375 |
<< " LEFT JOIN GL_ROOT " |
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
| 1376 |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
<< " LEFT JOIN GL_ROOT " |
| 1377 |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
| 1378 |
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
<< " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
| 1379 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
| 1380 |
|
} else { |
| 1381 |
|
oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
| 1382 |
|
<< " LEFT JOIN GL_ROOT " |
| 1383 |
|
<< " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
| 1384 |
|
<< " WHERE GL_RAW.PATH = '$PAM_RAW' AND " |
| 1385 |
|
<< " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
| 1386 |
|
}; |
| 1387 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1388 |
// |
// |
| 1389 |
if ( !result ) throw -12; |
if ( !result ) throw -12; |
| 1410 |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
| 1411 |
// |
// |
| 1412 |
oss.str(""); |
oss.str(""); |
| 1413 |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
if ( STATIC ){ |
| 1414 |
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
| 1415 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
| 1416 |
|
} else { |
| 1417 |
|
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
| 1418 |
|
<< this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')"; |
| 1419 |
|
}; |
| 1420 |
// |
// |
| 1421 |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
| 1422 |
// |
// |
| 1443 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 1444 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 1445 |
oss.str(""); |
oss.str(""); |
| 1446 |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
if ( STATIC ){ |
| 1447 |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
| 1448 |
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
<< " PATH = '" << this->GetRawPath().Data() << "' AND " |
| 1449 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 1450 |
|
} else { |
| 1451 |
|
oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE " |
| 1452 |
|
<< " PATH = '$PAM_RAW' AND " |
| 1453 |
|
<< " NAME = '" << this->GetRawFile().Data() << "' "; |
| 1454 |
|
}; |
| 1455 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1456 |
// |
// |
| 1457 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
| 1474 |
VarDumpEvent *vde = 0; |
VarDumpEvent *vde = 0; |
| 1475 |
VarDumpRecord *vdr = 0; |
VarDumpRecord *vdr = 0; |
| 1476 |
// |
// |
| 1477 |
|
Bool_t found = false; |
| 1478 |
trDumpEv->SetBranchAddress("VarDump", &vde); |
trDumpEv->SetBranchAddress("VarDump", &vde); |
| 1479 |
if ( trDumpEv->GetEntries() > 0 ){ |
if ( trDumpEv->GetEntries() > 0 ){ |
| 1480 |
Bool_t found = false; |
found = false; |
| 1481 |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){ |
| 1482 |
trDumpEv->GetEntry(i); |
trDumpEv->GetEntry(i); |
| 1483 |
vde->Records->GetEntries(); |
// vde->Records->GetEntries(); |
| 1484 |
if ( vde->Records->GetEntries()>0 ){ |
if ( vde->Records->GetEntries()>5 ){ |
| 1485 |
found = true; |
found = true; |
| 1486 |
goto fill; |
goto fill; |
| 1487 |
}; |
}; |
| 1494 |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
this->SetBOOTnumber((Int_t)vdr->VAR_VALUE); |
| 1495 |
// |
// |
| 1496 |
} else { |
} else { |
| 1497 |
if ( !this->GetBOOTnumber() ) return(4); |
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4); |
| 1498 |
|
}; |
| 1499 |
|
} else { |
| 1500 |
|
if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2); |
| 1501 |
|
}; |
| 1502 |
|
// |
| 1503 |
|
UInt_t bn = 0; |
| 1504 |
|
Bool_t afound = false; |
| 1505 |
|
if ( !found && this->AutoBoot()){ |
| 1506 |
|
afound = true; |
| 1507 |
|
// |
| 1508 |
|
// Search for other files with similar timesync |
| 1509 |
|
// |
| 1510 |
|
if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0); |
| 1511 |
|
UInt_t upperts = tsync-(obt0/1000)+5; |
| 1512 |
|
UInt_t lowerts = tsync-(obt0/1000)-5; |
| 1513 |
|
oss.str(""); |
| 1514 |
|
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)<" |
| 1515 |
|
<< upperts |
| 1516 |
|
<< " AND TIMESYNC-(OBT0/1000)>" |
| 1517 |
|
<< lowerts |
| 1518 |
|
<< " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;"; |
| 1519 |
|
result = conn->Query(oss.str().c_str()); |
| 1520 |
|
if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str()); |
| 1521 |
|
// |
| 1522 |
|
if ( !result ) throw -4;; |
| 1523 |
|
found = true; |
| 1524 |
|
if ( result->GetRowCount()<3 ){ |
| 1525 |
|
if ( IsDebug() ) printf(" AGH! no results!\n"); |
| 1526 |
|
found = false; |
| 1527 |
|
} else { |
| 1528 |
|
row = result->Next(); |
| 1529 |
|
bn = (UInt_t)atoll(row->GetField(0)); |
| 1530 |
|
for ( Int_t r=1; r<result->GetRowCount() ;r++){ |
| 1531 |
|
if ( !row ) throw -4; |
| 1532 |
|
if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0)); |
| 1533 |
|
if ( bn != (UInt_t)atoll(row->GetField(0)) ){ |
| 1534 |
|
if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0))); |
| 1535 |
|
found = false; |
| 1536 |
|
}; |
| 1537 |
|
row = result->Next(); |
| 1538 |
|
}; |
| 1539 |
}; |
}; |
| 1540 |
|
}; |
| 1541 |
|
// |
| 1542 |
|
Int_t sgn = 0; |
| 1543 |
|
// |
| 1544 |
|
if ( !found && !BOOTNO ){ |
| 1545 |
|
throw -29; |
| 1546 |
} else { |
} else { |
| 1547 |
if ( !this->GetBOOTnumber() ) return(2); |
if ( afound ){ |
| 1548 |
|
this->SetBOOTnumber(bn); |
| 1549 |
|
sgn = 8; |
| 1550 |
|
}; |
| 1551 |
}; |
}; |
| 1552 |
// |
// |
| 1553 |
oss.str(""); |
oss.str(""); |
| 1557 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
| 1558 |
// |
// |
| 1559 |
delete result; |
delete result; |
| 1560 |
return(0); |
return(sgn); |
| 1561 |
}; |
}; |
| 1562 |
|
|
| 1563 |
/** |
/** |
| 1570 |
stringstream oss; |
stringstream oss; |
| 1571 |
oss.str(""); |
oss.str(""); |
| 1572 |
// |
// |
| 1573 |
signal = this->SetUpperLimits(); |
// signal = this->SetUpperLimits(); |
| 1574 |
// |
// |
| 1575 |
// loop on runheader and runtrailer events |
// loop on runheader and runtrailer events |
| 1576 |
// |
// |
| 1588 |
rt->SetBranchAddress("RunTrailer", &runt); |
rt->SetBranchAddress("RunTrailer", &runt); |
| 1589 |
rt->SetBranchAddress("Header", &eht); |
rt->SetBranchAddress("Header", &eht); |
| 1590 |
// |
// |
| 1591 |
|
TTree *T = (TTree*)file->Get("Physics"); |
| 1592 |
|
if ( !T || T->IsZombie() ) throw -16; |
| 1593 |
|
EventHeader *eh = 0; |
| 1594 |
|
T->SetBranchAddress("Header", &eh); |
| 1595 |
|
// |
| 1596 |
|
if ( !(rh->GetEntries()) && !(rt->GetEntries()) && !(T->GetEntries()) ) return(16); |
| 1597 |
|
// |
| 1598 |
UInt_t obtt = 0; |
UInt_t obtt = 0; |
| 1599 |
UInt_t obth = 0; |
UInt_t obth = 0; |
| 1600 |
UInt_t pktt = 0; |
UInt_t pktt = 0; |
| 1628 |
// |
// |
| 1629 |
if ( !ptt && !(ptht+1) ){ |
if ( !ptt && !(ptht+1) ){ |
| 1630 |
// |
// |
| 1631 |
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); |
| 1632 |
// |
// |
| 1633 |
this->HandleRunFragments(true,false,0,(evbeft-1)); |
this->HandleRunFragments(true,false,0,(evbeft-1)); |
| 1634 |
// |
// |
| 1635 |
// |
// |
| 1636 |
} else if ( pth == ptht ){ |
} else if ( pth == ptht ){ |
| 1637 |
// |
// |
| 1638 |
if ( IsDebug() ) printf(" Missing header %i %i %i\n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" Missing header %u %u %u\n",ptht,pth,ptt); |
| 1639 |
// |
// |
| 1640 |
if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
| 1641 |
rt->GetEntry(ptt-1); |
rt->GetEntry(ptt-1); |
| 1644 |
rt->GetEntry(ptt); |
rt->GetEntry(ptt); |
| 1645 |
pht = eht->GetPscuHeader(); |
pht = eht->GetPscuHeader(); |
| 1646 |
// |
// |
| 1647 |
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); |
| 1648 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1649 |
// |
// |
| 1650 |
this->HandleMissingHoT(true,false,evbefh,evbeft-1); |
this->HandleMissingHoT(true,false,evbefh,evbeft-1); |
| 1651 |
// |
// |
| 1652 |
} else { |
} else { |
| 1653 |
// |
// |
| 1654 |
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); |
| 1655 |
// |
// |
| 1656 |
rh->GetEntry(ptht); |
rh->GetEntry(ptht); |
| 1657 |
phh = ehh->GetPscuHeader(); |
phh = ehh->GetPscuHeader(); |
| 1659 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1660 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
| 1661 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1662 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1663 |
// |
// |
| 1664 |
// handle this run |
// handle this run |
| 1665 |
// |
// |
| 1669 |
// |
// |
| 1670 |
if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){ |
if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){ |
| 1671 |
// |
// |
| 1672 |
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)); |
|
| 1673 |
// |
// |
| 1674 |
|
if ( evbefh == 0 ) { |
| 1675 |
|
// |
| 1676 |
|
signal = 8; |
| 1677 |
|
if ( IsDebug() ) printf(" Not supported yet: run with no events, no runtrailer, no runheader \n"); |
| 1678 |
|
// |
| 1679 |
|
} else { |
| 1680 |
|
// |
| 1681 |
|
this->HandleRunFragments(true,true,0,(evbefh-1)); |
| 1682 |
|
// |
| 1683 |
|
}; |
| 1684 |
}; |
}; |
| 1685 |
// |
// |
| 1686 |
// |
// |
| 1687 |
if ( (ptht - pth) > 1 ){ |
if ( (ptht - pth) > 1 ){ |
| 1688 |
// |
// |
| 1689 |
if ( IsDebug() ) printf(" Missing runtrailers! \n"); |
if ( IsDebug() ) printf(" Missing runtrailers! \n"); |
| 1690 |
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); |
| 1691 |
// is not the consecutive header |
// is not the consecutive header |
| 1692 |
while ( pth != ptht ){ |
while ( pth != ptht ){ |
| 1693 |
// |
// |
| 1708 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1709 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1710 |
// |
// |
| 1711 |
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); |
| 1712 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt); |
| 1713 |
// |
// |
| 1714 |
this->HandleMissingHoT(false,true,evbefh,evbeft-1); |
this->HandleMissingHoT(false,true,evbefh,evbeft-1); |
| 1715 |
// |
// |
| 1718 |
} else if ( !(ptht - pth) ){ |
} else if ( !(ptht - pth) ){ |
| 1719 |
// |
// |
| 1720 |
if ( IsDebug() ) printf(" Missing runheader! \n"); |
if ( IsDebug() ) printf(" Missing runheader! \n"); |
| 1721 |
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); |
| 1722 |
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"); |
| 1723 |
// |
// |
| 1724 |
} else { |
} else { |
| 1739 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
| 1740 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
| 1741 |
evbefh = cod->Get(pctp->Physics); |
evbefh = cod->Get(pctp->Physics); |
| 1742 |
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); |
| 1743 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i %i %i \n",ptht,pth,ptt); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u %u %u \n",ptht,pth,ptt); |
| 1744 |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i \n",rhev); |
if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u \n",rhev); |
| 1745 |
// |
// |
| 1746 |
this->HandleRunFragments(false,true,evbefh,upperentry); |
this->HandleRunFragments(false,true,evbefh,upperentry); |
| 1747 |
} else { |
} else { |
| 1789 |
// ( runhead_time <= _our_runhead_time && runtrail_time >= _our_runtrail_time) && |
// ( runhead_time <= _our_runhead_time && runtrail_time >= _our_runtrail_time) && |
| 1790 |
// ( runhead_obt <= _our_runheadobt || runhead_pkt <= _our_runheadpkt ) && |
// ( runhead_obt <= _our_runheadobt || runhead_pkt <= _our_runheadpkt ) && |
| 1791 |
// ( runtrail_obt <= _our_runtrailobt || runtrail_pkt <= _our_runtrailpkt ) ) |
// ( runtrail_obt <= _our_runtrailobt || runtrail_pkt <= _our_runtrailpkt ) ) |
| 1792 |
|
// || |
| 1793 |
|
// ( runhead_time = _our_runhead_time && runtrail_time = _our_runtrail_time && nevents > 100 ) |
| 1794 |
|
// ) |
| 1795 |
|
// || |
| 1796 |
|
// ( runhead_time = _our_runhead_time && runtrail_time > _our_runtrail_time && nevents > 100 ) |
| 1797 |
|
// ) |
| 1798 |
|
// || |
| 1799 |
|
// ( runhead_time < _our_runhead_time && runtrail_time = _our_runtrail_time && nevents > 100 ) |
| 1800 |
// ) |
// ) |
| 1801 |
// |
// |
| 1802 |
oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE " |
oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE " |
| 1812 |
<< " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR " |
<< " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR " |
| 1813 |
<< " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
<< " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
| 1814 |
<< " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR " |
<< " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR " |
| 1815 |
<< " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));"; |
<< " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") OR " |
| 1816 |
|
<< " (RUNHEADER_TIME=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " // these two lines in a certain way disable the patch below... |
| 1817 |
|
<< " RUNTRAILER_TIME=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100) OR" // |
| 1818 |
|
<< " (RUNHEADER_TIME=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " // |
| 1819 |
|
<< " RUNTRAILER_TIME>" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100) OR" // |
| 1820 |
|
<< " (RUNHEADER_TIME<" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " // |
| 1821 |
|
<< " RUNTRAILER_TIME=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100)" // |
| 1822 |
|
<< " ));"; |
| 1823 |
// |
// |
| 1824 |
if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str()); |
| 1825 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 1841 |
// |
// |
| 1842 |
// the run has already been inserted |
// the run has already been inserted |
| 1843 |
// |
// |
| 1844 |
// return(true); //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
| 1845 |
|
return(true); //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here |
| 1846 |
// |
// |
| 1847 |
// PATCH! |
// PATCH! |
| 1848 |
// we keep the processing run if (in order of growing importance) 1) we have the runtrailer while the old run doesn't have it 2) we have the runheader |
// we keep the processing run if (in order of growing importance) 1) we have the runtrailer while the old run doesn't have it 2) we have the runheader |
| 1851 |
if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){ |
if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){ |
| 1852 |
// |
// |
| 1853 |
if ( IsDebug() ) printf(" The new run has more events than the old one \n"); |
if ( IsDebug() ) printf(" The new run has more events than the old one \n"); |
| 1854 |
oss.str(""); |
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN"); |
| 1855 |
oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
// oss.str(""); |
| 1856 |
if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
// oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
| 1857 |
conn->Query(oss.str().c_str()); |
// if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
| 1858 |
|
// conn->Query(oss.str().c_str()); |
| 1859 |
if ( signal ) signal = false; |
if ( signal ) signal = false; |
| 1860 |
goto gonext; |
goto gonext; |
| 1861 |
// |
// |
| 1869 |
// |
// |
| 1870 |
if ( IsDebug() ) printf(" The new run has the same number of events and the runheader the old one miss the runheader \n"); |
if ( IsDebug() ) printf(" The new run has the same number of events and the runheader the old one miss the runheader \n"); |
| 1871 |
// |
// |
| 1872 |
oss.str(""); |
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN"); |
| 1873 |
oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
// oss.str(""); |
| 1874 |
if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
// oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
| 1875 |
conn->Query(oss.str().c_str()); |
// if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
| 1876 |
|
// conn->Query(oss.str().c_str()); |
| 1877 |
// |
// |
| 1878 |
if ( signal ) signal = false; |
if ( signal ) signal = false; |
| 1879 |
goto gonext; |
goto gonext; |
| 1887 |
// |
// |
| 1888 |
if ( IsDebug() ) printf(" The new run has the same number of events, the runheader and the runtrailer the old one miss the runtrailer \n"); |
if ( IsDebug() ) printf(" The new run has the same number of events, the runheader and the runtrailer the old one miss the runtrailer \n"); |
| 1889 |
// |
// |
| 1890 |
oss.str(""); |
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN"); |
| 1891 |
oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
// oss.str(""); |
| 1892 |
if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
// oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";"; |
| 1893 |
conn->Query(oss.str().c_str()); |
// if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str()); |
| 1894 |
|
// conn->Query(oss.str().c_str()); |
| 1895 |
if ( signal ) signal = false; |
if ( signal ) signal = false; |
| 1896 |
// |
// |
| 1897 |
}; |
}; |
| 1906 |
delete result; |
delete result; |
| 1907 |
// |
// |
| 1908 |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
if ( signal && IsDebug() ) printf(" The run has already been inserted \n"); |
| 1909 |
|
if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n"); |
| 1910 |
return(signal); |
return(signal); |
| 1911 |
}; |
}; |
| 1912 |
|
|
| 1955 |
UInt_t rhfirstev = firstev; |
UInt_t rhfirstev = firstev; |
| 1956 |
UInt_t rtlastev = lastev; |
UInt_t rtlastev = lastev; |
| 1957 |
Bool_t found = false; |
Bool_t found = false; |
| 1958 |
|
Bool_t foundinrun = false; |
| 1959 |
// |
// |
| 1960 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 1961 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 1982 |
// |
// |
| 1983 |
if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n"); |
if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n"); |
| 1984 |
// |
// |
| 1985 |
// can we find the other piece of the run in the GL_RUN_FRAGMENTS table? |
// First of all insert the run in the fragment table... |
| 1986 |
// |
// |
| 1987 |
if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is |
oss.str(""); |
| 1988 |
// missing it no way we can found a piece in the frag table |
oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE " |
| 1989 |
|
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND (" |
| 1990 |
|
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
| 1991 |
|
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
| 1992 |
|
<< " RUNHEADER_OBT>=" << glrun->GetRUNHEADER_OBT() << " OR " |
| 1993 |
|
<< " RUNHEADER_PKT>=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
| 1994 |
|
<< " RUNTRAILER_OBT<=" << glrun->GetRUNTRAILER_OBT() << " OR " |
| 1995 |
|
<< " RUNTRAILER_PKT<=" << glrun->GetRUNTRAILER_PKT() << ") ) OR " |
| 1996 |
|
<< " (RUNHEADER_TIME<=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " |
| 1997 |
|
<< " RUNTRAILER_TIME>=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND (" |
| 1998 |
|
<< " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR " |
| 1999 |
|
<< " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
| 2000 |
|
<< " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR " |
| 2001 |
|
<< " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));"; |
| 2002 |
|
// |
| 2003 |
|
if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str()); |
| 2004 |
|
result = conn->Query(oss.str().c_str()); |
| 2005 |
|
// |
| 2006 |
|
if ( !result ) throw -4; |
| 2007 |
|
// |
| 2008 |
|
row = result->Next(); |
| 2009 |
|
// |
| 2010 |
|
if ( !row ){ |
| 2011 |
// |
// |
| 2012 |
oss.str(""); |
// no, insert this run in the GL_RUN_FRAGMENTS table (check if exist before!) |
| 2013 |
oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE " |
// |
| 2014 |
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
if ( IsDebug() ) printf(" The run is new \n"); |
| 2015 |
<< " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() |
if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n"); |
| 2016 |
<< " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!! |
// |
| 2017 |
|
glrun->SetID(this->AssignRunID()); |
| 2018 |
|
glrun->SetID_RUN_FRAG(0); |
| 2019 |
|
glrun->Fill_GL_RUN_FRAGMENTS(conn); |
| 2020 |
|
// |
| 2021 |
|
} else { |
| 2022 |
|
if ( IsDebug() ) printf(" The run is already present in the fragment table \n"); |
| 2023 |
|
return; |
| 2024 |
|
}; |
| 2025 |
|
// |
| 2026 |
|
// |
| 2027 |
|
// can we find the other piece of the run in the GL_RUN_FRAGMENTS table? |
| 2028 |
|
// |
| 2029 |
|
if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is |
| 2030 |
|
// missing it no way we can found a piece in the frag table |
| 2031 |
|
// |
| 2032 |
|
oss.str(""); |
| 2033 |
|
oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE " |
| 2034 |
|
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
| 2035 |
|
<< " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " |
| 2036 |
|
<< " ID != " << glrun->ID |
| 2037 |
|
<< " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!! |
| 2038 |
// |
// |
| 2039 |
if ( IsDebug() ) printf(" look for runheader in the fragments table: query is \n %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" look for runheader in the fragments table: query is \n %s \n",oss.str().c_str()); |
| 2040 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
| 2043 |
// |
// |
| 2044 |
row = result->Next(); |
row = result->Next(); |
| 2045 |
// |
// |
| 2046 |
|
if ( !row && NoFrag() ){ |
| 2047 |
|
// |
| 2048 |
|
oss.str(""); |
| 2049 |
|
oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN WHERE " |
| 2050 |
|
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
| 2051 |
|
<< " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " |
| 2052 |
|
<< " ID != " << glrun->ID |
| 2053 |
|
<< " AND ID=ID_RUN_FRAG ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!! |
| 2054 |
|
// |
| 2055 |
|
if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str()); |
| 2056 |
|
result = conn->Query(oss.str().c_str()); |
| 2057 |
|
// |
| 2058 |
|
if ( !result ) throw -4; |
| 2059 |
|
// |
| 2060 |
|
foundinrun = true; |
| 2061 |
|
// |
| 2062 |
|
row = result->Next(); |
| 2063 |
|
// |
| 2064 |
|
}; |
| 2065 |
|
// |
| 2066 |
if ( !row ){ |
if ( !row ){ |
| 2067 |
if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n"); |
if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n"); |
| 2068 |
found = false; |
found = false; |
| 2137 |
// |
// |
| 2138 |
if ( IsDebug() ) printf(" now you can handle the piece of the run \n "); |
if ( IsDebug() ) printf(" now you can handle the piece of the run \n "); |
| 2139 |
// |
// |
| 2140 |
|
if ( foundinrun ){ |
| 2141 |
|
glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS"); |
| 2142 |
|
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN"); |
| 2143 |
|
}; |
| 2144 |
|
// |
| 2145 |
GL_RUN *glrun1 = new GL_RUN(); |
GL_RUN *glrun1 = new GL_RUN(); |
| 2146 |
// |
// |
| 2147 |
UInt_t idfrag = (UInt_t)atoll(row->GetField(0)); |
// UInt_t idfrag = (UInt_t)atoll(row->GetField(0)); |
| 2148 |
// |
// |
| 2149 |
oss.str(""); |
oss.str(""); |
| 2150 |
oss << " ID="<<row->GetField(0)<<";"; |
oss << " ID="<<row->GetField(0)<<";"; |
| 2157 |
ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT()); |
ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT()); |
| 2158 |
UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT()); |
UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT()); |
| 2159 |
ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT()); |
ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT()); |
| 2160 |
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); |
| 2161 |
TTree *T= 0; |
TTree *T= 0; |
| 2162 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
| 2163 |
if ( !T || T->IsZombie() ) throw -16; |
if ( !T || T->IsZombie() ) throw -16; |
| 2171 |
bobt = OBT(ph->GetOrbitalTime()); |
bobt = OBT(ph->GetOrbitalTime()); |
| 2172 |
firstev++; |
firstev++; |
| 2173 |
}; |
}; |
| 2174 |
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); |
| 2175 |
// |
// |
| 2176 |
glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER()); |
glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER()); |
| 2177 |
glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER()); |
glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER()); |
| 2181 |
// |
// |
| 2182 |
glrun->SetEV_FROM(firstev); |
glrun->SetEV_FROM(firstev); |
| 2183 |
glrun->SetNEVENTS(lastev-firstev+1); |
glrun->SetNEVENTS(lastev-firstev+1); |
| 2184 |
|
// |
| 2185 |
glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME()); |
glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME()); |
| 2186 |
glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT()); |
glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT()); |
| 2187 |
glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT()); |
glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT()); |
| 2201 |
// |
// |
| 2202 |
if ( !IsRunAlreadyInserted() ){ |
if ( !IsRunAlreadyInserted() ){ |
| 2203 |
// |
// |
| 2204 |
glrun->SetID(this->AssignRunID()); |
// glrun->SetID(this->AssignRunID()); |
| 2205 |
glrun->SetID_RUN_FRAG(glrun1->GetID()); |
glrun->SetID_RUN_FRAG(glrun1->GetID()); |
| 2206 |
glrun->Fill_GL_RUN(conn); |
glrun->Fill_GL_RUN(conn); |
| 2207 |
// |
// |
| 2208 |
// get id number |
// set id number |
|
// |
|
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND " |
|
|
// << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND " |
|
|
// << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str()); |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// UInt_t idrun0 = 0; |
|
|
// if ( !row ){ |
|
|
// throw -10; |
|
|
// } else { |
|
|
// idrun0 = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
| 2209 |
// |
// |
|
// glrun1->SetID_RUN_FRAG(idrun0); |
|
| 2210 |
glrun1->SetID_RUN_FRAG(glrun->GetID()); |
glrun1->SetID_RUN_FRAG(glrun->GetID()); |
| 2211 |
glrun1->Fill_GL_RUN(conn); |
glrun1->Fill_GL_RUN(conn); |
| 2212 |
// |
// |
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str()); |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// UInt_t idrun1 = 0; |
|
|
// if ( !row ){ |
|
|
// throw -10; |
|
|
// } else { |
|
|
// idrun1 = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;"; |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// |
|
| 2213 |
}; |
}; |
|
// |
|
|
delete glrun1; |
|
|
// |
|
| 2214 |
// delete old entry in fragment table |
// delete old entry in fragment table |
| 2215 |
// |
// |
| 2216 |
oss.str(""); |
glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
| 2217 |
|
glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
| 2218 |
// |
// |
| 2219 |
oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";"; |
delete glrun1; |
| 2220 |
// |
// |
|
if ( IsDebug() ) printf(" Delete from frag table the old run :\n query is \n %s \n",oss.str().c_str()); |
|
|
result = conn->Query(oss.str().c_str()); |
|
|
if ( !result ) throw -4; |
|
| 2221 |
// |
// |
| 2222 |
return; |
return; |
| 2223 |
// |
// |
| 2226 |
}; |
}; |
| 2227 |
// |
// |
| 2228 |
if ( mistrail && rtlastev == lastev ) { // look for runtrailer (only when at the end of the file, if at the beginning and the runh is |
if ( mistrail && rtlastev == lastev ) { // look for runtrailer (only when at the end of the file, if at the beginning and the runh is |
| 2229 |
// missing it no way we can found a piece in the frag table |
// missing it no way we can found a piece in the frag table |
| 2230 |
// |
// |
| 2231 |
oss.str(""); |
oss.str(""); |
| 2232 |
oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN_FRAGMENTS WHERE " |
oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN_FRAGMENTS WHERE " |
| 2233 |
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
| 2234 |
<< " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() |
<< " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND " |
| 2235 |
|
<< " ID != " << glrun->ID |
| 2236 |
<< " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;"; |
<< " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;"; |
| 2237 |
// |
// |
| 2238 |
if ( IsDebug() ) printf(" look for runtrailer in the fragments table: query is \n %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" look for runtrailer in the fragments table: query is \n %s \n",oss.str().c_str()); |
| 2242 |
// |
// |
| 2243 |
row = result->Next(); |
row = result->Next(); |
| 2244 |
// |
// |
| 2245 |
|
if ( !row && NoFrag() ){ |
| 2246 |
|
// |
| 2247 |
|
oss.str(""); |
| 2248 |
|
oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN WHERE " |
| 2249 |
|
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND " |
| 2250 |
|
<< " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND " |
| 2251 |
|
<< " ID != " << glrun->ID |
| 2252 |
|
<< " AND ID=ID_RUN_FRAG ORDER BY RUNTRAILER_TIME ASC LIMIT 1;"; |
| 2253 |
|
// |
| 2254 |
|
if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str()); |
| 2255 |
|
result = conn->Query(oss.str().c_str()); |
| 2256 |
|
// |
| 2257 |
|
if ( !result ) throw -4; |
| 2258 |
|
// |
| 2259 |
|
foundinrun = true; |
| 2260 |
|
row = result->Next(); |
| 2261 |
|
// |
| 2262 |
|
}; |
| 2263 |
|
// |
| 2264 |
if ( !row ){ |
if ( !row ){ |
| 2265 |
if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n"); |
if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n"); |
| 2266 |
found = false; |
found = false; |
| 2334 |
// |
// |
| 2335 |
if ( IsDebug() ) printf(" now you can handle the piece of the run \n "); |
if ( IsDebug() ) printf(" now you can handle the piece of the run \n "); |
| 2336 |
// |
// |
| 2337 |
|
if ( foundinrun ){ |
| 2338 |
|
glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS"); |
| 2339 |
|
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN"); |
| 2340 |
|
}; |
| 2341 |
|
// |
| 2342 |
GL_RUN *glrun1 = new GL_RUN(); |
GL_RUN *glrun1 = new GL_RUN(); |
| 2343 |
// |
// |
| 2344 |
UInt_t idfrag = (UInt_t)atoll(row->GetField(0)); |
// UInt_t idfrag = (UInt_t)atoll(row->GetField(0)); |
| 2345 |
// |
// |
| 2346 |
oss.str(""); |
oss.str(""); |
| 2347 |
oss << " ID="<<row->GetField(0)<<";"; |
oss << " ID="<<row->GetField(0)<<";"; |
| 2354 |
ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT()); |
ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT()); |
| 2355 |
UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT()); |
UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT()); |
| 2356 |
ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT()); |
ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT()); |
| 2357 |
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); |
| 2358 |
TTree *T= 0; |
TTree *T= 0; |
| 2359 |
T = (TTree*)file->Get("Physics"); |
T = (TTree*)file->Get("Physics"); |
| 2360 |
if ( !T || T->IsZombie() ) throw -16; |
if ( !T || T->IsZombie() ) throw -16; |
| 2368 |
aobt = OBT(ph->GetOrbitalTime()); |
aobt = OBT(ph->GetOrbitalTime()); |
| 2369 |
lastev--; |
lastev--; |
| 2370 |
}; |
}; |
| 2371 |
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); |
| 2372 |
// |
// |
| 2373 |
glrun->SetEV_TO(lastev); |
glrun->SetEV_TO(lastev); |
| 2374 |
glrun->SetNEVENTS(lastev-firstev+1); |
glrun->SetNEVENTS(lastev-firstev+1); |
| 2397 |
// |
// |
| 2398 |
if ( !IsRunAlreadyInserted() ){ |
if ( !IsRunAlreadyInserted() ){ |
| 2399 |
// |
// |
| 2400 |
glrun->SetID(this->AssignRunID()); |
// glrun->SetID(this->AssignRunID()); |
| 2401 |
// |
// |
| 2402 |
glrun->SetID_RUN_FRAG(glrun1->GetID()); |
glrun->SetID_RUN_FRAG(glrun1->GetID()); |
| 2403 |
glrun->Fill_GL_RUN(conn); |
glrun->Fill_GL_RUN(conn); |
| 2404 |
// |
// |
| 2405 |
// get id number |
// set id number |
| 2406 |
// |
// |
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND " |
|
|
// << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND " |
|
|
// << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str()); |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// UInt_t idrun0 = 0; |
|
|
// if ( !row ){ |
|
|
// throw -10; |
|
|
// } else { |
|
|
// idrun0 = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
|
// // |
|
|
// glrun1->SetID_RUN_FRAG(idrun0); |
|
|
// |
|
| 2407 |
glrun1->SetID_RUN_FRAG(glrun->GetID()); |
glrun1->SetID_RUN_FRAG(glrun->GetID()); |
| 2408 |
glrun1->Fill_GL_RUN(conn); |
glrun1->Fill_GL_RUN(conn); |
| 2409 |
// |
// |
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str()); |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// UInt_t idrun1 = 0; |
|
|
// if ( !row ){ |
|
|
// throw -10; |
|
|
// } else { |
|
|
// idrun1 = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;"; |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -4; |
|
|
// |
|
| 2410 |
}; |
}; |
| 2411 |
// |
// |
| 2412 |
delete glrun1; |
// delete old entries in fragment table |
|
// |
|
|
// delete old entry in fragment table |
|
|
// |
|
|
oss.str(""); |
|
| 2413 |
// |
// |
| 2414 |
oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";"; |
glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
| 2415 |
// |
glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
|
if ( IsDebug() ) printf(" Delete from frag table the old run :\n query is \n %s \n",oss.str().c_str()); |
|
|
result = conn->Query(oss.str().c_str()); |
|
|
if ( !result ) throw -4; |
|
| 2416 |
// |
// |
| 2417 |
|
delete glrun1; |
| 2418 |
// |
// |
| 2419 |
return; |
return; |
| 2420 |
// |
// |
| 2451 |
// |
// |
| 2452 |
row = result->Next(); |
row = result->Next(); |
| 2453 |
// |
// |
| 2454 |
if ( !row ){ |
if ( row ){ |
|
// |
|
|
oss.str(""); |
|
|
oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE " |
|
|
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND (" |
|
|
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
|
|
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
|
|
<< " RUNHEADER_OBT>=" << glrun->GetRUNHEADER_OBT() << " OR " |
|
|
<< " RUNHEADER_PKT>=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
|
|
<< " RUNTRAILER_OBT<=" << glrun->GetRUNTRAILER_OBT() << " OR " |
|
|
<< " RUNTRAILER_PKT<=" << glrun->GetRUNTRAILER_PKT() << ") ) OR " |
|
|
<< " (RUNHEADER_TIME<=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND " |
|
|
<< " RUNTRAILER_TIME>=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND (" |
|
|
<< " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR " |
|
|
<< " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND (" |
|
|
<< " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR " |
|
|
<< " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));"; |
|
|
// |
|
|
if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str()); |
|
|
result = conn->Query(oss.str().c_str()); |
|
|
// |
|
|
if ( !result ) throw -4; |
|
|
// |
|
|
row = result->Next(); |
|
|
// |
|
|
if ( !row ){ |
|
|
// |
|
|
// no, insert this run in the GL_RUN_FRAGMENTS table (check if exist before!) |
|
|
// |
|
|
if ( IsDebug() ) printf(" The run is new \n"); |
|
|
if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n"); |
|
|
// |
|
|
glrun->SetID(this->AssignRunID()); |
|
|
glrun->SetID_RUN_FRAG(0); |
|
|
glrun->Fill_GL_RUN_FRAGMENTS(conn); |
|
|
// |
|
|
} else { |
|
|
if ( IsDebug() ) printf(" The run is already present in the fragment table \n"); |
|
|
}; |
|
|
} else { |
|
| 2455 |
if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n"); |
if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n"); |
| 2456 |
|
} else { |
| 2457 |
|
if ( NoFrag() ){ |
| 2458 |
|
glrun->SetID_RUN_FRAG(glrun->GetID()); |
| 2459 |
|
glrun->Fill_GL_RUN(conn); |
| 2460 |
|
glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
| 2461 |
|
}; |
| 2462 |
}; |
}; |
| 2463 |
}; |
}; |
| 2464 |
// |
// |
| 2975 |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
| 2976 |
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) ){ |
| 2977 |
// |
// |
| 2978 |
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); |
| 2979 |
// |
// |
| 2980 |
// check if the calibration has already been inserted |
// check if the calibration has already been inserted |
| 2981 |
// |
// |
| 3019 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
| 3020 |
// |
// |
| 3021 |
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); |
| 3022 |
if ( fromtime < 1150863400 ){ |
if ( fromtime < 1150871000 ){ //1150866904 |
| 3023 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime); |
| 3024 |
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 |
| 3025 |
}; |
}; |
| 3026 |
// |
// |
| 3082 |
// |
// |
| 3083 |
} else { |
} else { |
| 3084 |
// |
// |
| 3085 |
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); |
| 3086 |
// |
// |
| 3087 |
}; |
}; |
| 3088 |
// |
// |
| 3164 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
| 3165 |
// |
// |
| 3166 |
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"); |
| 3167 |
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 |
| 3168 |
// |
// |
| 3169 |
oss.str(""); |
oss.str(""); |
| 3170 |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
| 3293 |
pkt1 = ph1->GetCounter(); |
pkt1 = ph1->GetCounter(); |
| 3294 |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph1->GetOrbitalTime()); |
| 3295 |
// |
// |
| 3296 |
valid = 1; |
// valid = 1; |
| 3297 |
// |
// // |
| 3298 |
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 |
| 3299 |
// |
// |
| 3300 |
// |
// |
| 3301 |
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) ){ |
| 3302 |
// |
// |
| 3303 |
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); |
| 3304 |
|
// |
| 3305 |
|
valid = ValidateTrkCalib( caltrk1, eh1 ); |
| 3306 |
|
if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl; |
| 3307 |
// |
// |
| 3308 |
// Do we have the second calibration packet? |
// Do we have the second calibration packet? |
| 3309 |
// |
// |
| 3320 |
obt2 = ph2->GetOrbitalTime(); |
obt2 = ph2->GetOrbitalTime(); |
| 3321 |
pkt2 = ph2->GetCounter(); |
pkt2 = ph2->GetCounter(); |
| 3322 |
// |
// |
| 3323 |
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 |
| 3324 |
// |
// |
| 3325 |
} else { |
} else { |
| 3326 |
// |
// |
| 3343 |
// |
// |
| 3344 |
}; |
}; |
| 3345 |
// |
// |
| 3346 |
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); |
| 3347 |
// |
// |
| 3348 |
// The calibration is good |
// The calibration is good |
| 3349 |
// |
// |
| 3350 |
if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){ |
if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){ |
| 3351 |
// |
// |
| 3352 |
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); |
| 3353 |
|
// |
| 3354 |
|
UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 ); |
| 3355 |
|
if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl; |
| 3356 |
|
valid = valid & valid2; |
| 3357 |
// |
// |
| 3358 |
// Handle good calib |
// Handle good calib |
| 3359 |
// |
// |
| 3363 |
// |
// |
| 3364 |
if ( t2 != pret2+1 ){ |
if ( t2 != pret2+1 ){ |
| 3365 |
// |
// |
| 3366 |
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); |
| 3367 |
// |
// |
| 3368 |
while ( t2 > pret2+1 ){ |
while ( t2 > pret2+1 ){ |
| 3369 |
// |
// |
| 3392 |
// |
// |
| 3393 |
// Check for missing calibtrk2 |
// Check for missing calibtrk2 |
| 3394 |
// |
// |
| 3395 |
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); |
| 3396 |
t2 = pret2; |
t2 = pret2; |
| 3397 |
// |
// |
| 3398 |
// handle missing calib2 |
// handle missing calib2 |
| 3406 |
// |
// |
| 3407 |
} else { |
} else { |
| 3408 |
// |
// |
| 3409 |
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); |
| 3410 |
// |
// |
| 3411 |
}; |
}; |
| 3412 |
// |
// |
| 3431 |
valid = 0; |
valid = 0; |
| 3432 |
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) ){ |
| 3433 |
// |
// |
| 3434 |
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); |
| 3435 |
// |
// |
| 3436 |
this->HandleTRK_CALIB(false,true); |
this->HandleTRK_CALIB(false,true); |
| 3437 |
// |
// |
| 3485 |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
fromtime = this->GetAbsTime(ph->GetOrbitalTime()); |
| 3486 |
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) ){ |
| 3487 |
// |
// |
| 3488 |
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); |
| 3489 |
// |
// |
| 3490 |
// check if the calibration has already been inserted |
// check if the calibration has already been inserted |
| 3491 |
// |
// |
| 3527 |
// no calibrations in the db contain our calibration |
// no calibrations in the db contain our calibration |
| 3528 |
// |
// |
| 3529 |
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"); |
| 3530 |
if ( fromtime < 1150863400 ){ |
if ( fromtime < 1150871000 ){ |
| 3531 |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime); |
if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime); |
| 3532 |
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 |
| 3533 |
}; |
}; |
| 3534 |
// |
// |
| 3587 |
// |
// |
| 3588 |
} else { |
} else { |
| 3589 |
// |
// |
| 3590 |
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); |
| 3591 |
// |
// |
| 3592 |
}; |
}; |
| 3593 |
// |
// |
| 3600 |
* Scan the fragment table and move old fragments to the GL_RUN table |
* Scan the fragment table and move old fragments to the GL_RUN table |
| 3601 |
*/ |
*/ |
| 3602 |
Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){ |
Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){ |
| 3603 |
|
return(this->CleanGL_RUN_FRAGMENTS("")); |
| 3604 |
|
}; |
| 3605 |
|
|
| 3606 |
|
/** |
| 3607 |
|
* Scan the fragment table and move old fragments to the GL_RUN table |
| 3608 |
|
*/ |
| 3609 |
|
Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(TString fcleanfile){ |
| 3610 |
// |
// |
| 3611 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 3612 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 3614 |
TSQLRow *row2 = 0; |
TSQLRow *row2 = 0; |
| 3615 |
// |
// |
| 3616 |
UInt_t moved = 0; |
UInt_t moved = 0; |
|
// UInt_t timelim = 0; |
|
|
// TDatime *time = new TDatime(); |
|
| 3617 |
// |
// |
| 3618 |
stringstream oss; |
stringstream oss; |
| 3619 |
oss.str(""); |
oss.str(""); |
| 3620 |
// |
// |
| 3621 |
// |
if ( !strcmp(fcleanfile.Data(),"") ){ |
| 3622 |
// check if there are entries older than "olderthan" seconds from now |
// |
| 3623 |
// |
// check if there are entries older than "olderthan" seconds from now |
| 3624 |
oss.str(""); |
// |
| 3625 |
oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE" |
oss.str(""); |
| 3626 |
<< " INSERT_TIME <= '" << clean_time->AsSQLString() << "';"; |
oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE" |
| 3627 |
// |
<< " INSERT_TIME <= '" << clean_time->AsSQLString() << "';"; |
| 3628 |
if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str()); |
// |
| 3629 |
result = conn->Query(oss.str().c_str()); |
if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str()); |
| 3630 |
|
result = conn->Query(oss.str().c_str()); |
| 3631 |
|
// |
| 3632 |
|
} else { |
| 3633 |
|
oss.str(""); |
| 3634 |
|
oss << " SELECT ID FROM GL_ROOT WHERE NAME='" << fcleanfile.Data() << "';"; |
| 3635 |
|
if ( IsDebug() ) printf(" Getting ID_ROOT_L0 query %s \n",oss.str().c_str()); |
| 3636 |
|
result = conn->Query(oss.str().c_str()); |
| 3637 |
|
// |
| 3638 |
|
if ( result ){ |
| 3639 |
|
// |
| 3640 |
|
row = result->Next(); |
| 3641 |
|
// |
| 3642 |
|
if ( row ){ |
| 3643 |
|
oss.str(""); |
| 3644 |
|
oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE" |
| 3645 |
|
<< " ID_ROOT_L0=" << row->GetField(0) << ";"; |
| 3646 |
|
// |
| 3647 |
|
if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS for ROOT file query is \n %s \n",oss.str().c_str()); |
| 3648 |
|
result = conn->Query(oss.str().c_str()); |
| 3649 |
|
// |
| 3650 |
|
}; |
| 3651 |
|
} else { |
| 3652 |
|
return(2); |
| 3653 |
|
}; |
| 3654 |
|
}; |
| 3655 |
// |
// |
| 3656 |
if ( result ){ |
if ( result ){ |
| 3657 |
// |
// |
| 3728 |
if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0)); |
if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0)); |
| 3729 |
// |
// |
| 3730 |
// |
// |
| 3731 |
oss.str(""); |
glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS"); |
| 3732 |
oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0); |
// oss.str(""); |
| 3733 |
if ( IsDebug() ) printf(" Clean the GL_RUN_FRAGMENTS table: query is \n %s \n",oss.str().c_str()); |
// oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0); |
| 3734 |
result2 = conn->Query(oss.str().c_str()); |
// if ( IsDebug() ) printf(" Clean the GL_RUN_FRAGMENTS table: query is \n %s \n",oss.str().c_str()); |
| 3735 |
// |
// result2 = conn->Query(oss.str().c_str()); |
| 3736 |
if ( !result2 ) throw -4; |
// // |
| 3737 |
// |
// if ( !result2 ) throw -4; |
| 3738 |
|
// // |
| 3739 |
row = result->Next(); |
row = result->Next(); |
| 3740 |
}; |
}; |
| 3741 |
}; |
}; |
| 3747 |
* Check if runs are good, i.e. if the tracker calibration is correctly associated.. |
* Check if runs are good, i.e. if the tracker calibration is correctly associated.. |
| 3748 |
*/ |
*/ |
| 3749 |
Int_t PamelaDBOperations::ValidateRuns(){ |
Int_t PamelaDBOperations::ValidateRuns(){ |
| 3750 |
|
return(this->ValidateRuns("")); |
| 3751 |
|
}; |
| 3752 |
|
|
| 3753 |
|
/** |
| 3754 |
|
* Check if runs are good, i.e. if the tracker calibration is correctly associated.. |
| 3755 |
|
*/ |
| 3756 |
|
Int_t PamelaDBOperations::ValidateRuns(TString valfile){ |
| 3757 |
// |
// |
| 3758 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 3759 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
| 3760 |
// |
// |
| 3761 |
|
UInt_t calibtime = 50; |
| 3762 |
|
// |
| 3763 |
stringstream oss; |
stringstream oss; |
| 3764 |
oss.str(""); |
oss.str(""); |
| 3765 |
// |
// |
| 3768 |
// ======================================================= |
// ======================================================= |
| 3769 |
UInt_t t_stop = 0; |
UInt_t t_stop = 0; |
| 3770 |
UInt_t t_start = 0; |
UInt_t t_start = 0; |
| 3771 |
// -------------------------------------------------------------- |
if ( !strcmp(valfile.Data(),"") ) { |
| 3772 |
// 1) get the OBT of the last run inserted after clean-time limit |
// -------------------------------------------------------------- |
| 3773 |
// -------------------------------------------------------------- |
// 1) get the OBT of the last run inserted after clean-time limit |
| 3774 |
oss.str(""); |
// -------------------------------------------------------------- |
| 3775 |
oss << " SELECT * FROM GL_RUN WHERE INSERT_TIME <= '" << clean_time->AsSQLString() |
oss.str(""); |
| 3776 |
<< "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
oss << " SELECT * FROM GL_RUN WHERE INSERT_TIME <= '" << clean_time->AsSQLString() |
| 3777 |
if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str()); |
<< "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
| 3778 |
result = conn->Query(oss.str().c_str()); |
if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str()); |
| 3779 |
if ( !result ) throw -4; |
result = conn->Query(oss.str().c_str()); |
| 3780 |
if ( !result->GetRowCount() ) { |
if ( !result ) throw -4; |
| 3781 |
printf(" No runs to validate \n"); |
if ( !result->GetRowCount() ) { |
| 3782 |
return(1); |
printf(" No runs to validate \n"); |
| 3783 |
}else{ |
return(1); |
| 3784 |
row = result->Next(); |
}else{ |
| 3785 |
t_start = (UInt_t)atoll(row->GetField(4)); |
row = result->Next(); |
| 3786 |
}; |
t_start = (UInt_t)atoll(row->GetField(4)); |
| 3787 |
// -------------------------------------------------------------- |
}; |
| 3788 |
// 2) get the OBT of the last validated run |
// -------------------------------------------------------------- |
| 3789 |
// -------------------------------------------------------------- |
// 2) get the OBT of the last validated run |
| 3790 |
oss.str(""); |
// -------------------------------------------------------------- |
| 3791 |
oss << " SELECT * FROM GL_RUN WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start |
oss.str(""); |
| 3792 |
<<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
oss << " SELECT * FROM GL_RUN WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start |
| 3793 |
if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str()); |
<<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
| 3794 |
result = conn->Query(oss.str().c_str()); |
if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str()); |
| 3795 |
if ( !result ) throw -4; |
result = conn->Query(oss.str().c_str()); |
| 3796 |
if ( result->GetRowCount() ){ |
if ( !result ) throw -4; |
| 3797 |
row = result->Next(); |
if ( result->GetRowCount() ){ |
| 3798 |
t_stop = (UInt_t)atoll(row->GetField(4)); |
row = result->Next(); |
| 3799 |
|
t_stop = (UInt_t)atoll(row->GetField(4)); |
| 3800 |
|
}; |
| 3801 |
|
if ( IsDebug() ) printf("Validation interval: from time %u - to time %u \n\n",t_stop,t_start); |
| 3802 |
|
// -------------------------------------------------------------- |
| 3803 |
|
// now retrieves runs to be validated |
| 3804 |
|
// -------------------------------------------------------------- |
| 3805 |
|
oss.str(""); |
| 3806 |
|
oss << " SELECT * FROM GL_RUN WHERE RUNHEADER_TIME <=" << t_start; |
| 3807 |
|
oss << " AND RUNHEADER_TIME >="<< t_stop; |
| 3808 |
|
oss << " ORDER BY RUNHEADER_TIME DESC;"; |
| 3809 |
|
if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str()); |
| 3810 |
|
result = conn->Query(oss.str().c_str()); |
| 3811 |
|
} else { |
| 3812 |
|
// |
| 3813 |
|
stringstream myquery; |
| 3814 |
|
UInt_t myid = 0; |
| 3815 |
|
myquery.str(""); |
| 3816 |
|
myquery << " SELECT ID FROM GL_ROOT where NAME='"<<valfile.Data() <<"';"; |
| 3817 |
|
// |
| 3818 |
|
result = conn->Query(myquery.str().c_str()); |
| 3819 |
|
// |
| 3820 |
|
row = result->Next(); |
| 3821 |
|
if( !row ){ |
| 3822 |
|
if ( strcmp(valfile.Data(),GetRootName().Data()) ){ |
| 3823 |
|
if ( IsDebug() ) printf(" No file to be validated even if option \"-validate file\" was used!!\n"); |
| 3824 |
|
return(2); |
| 3825 |
|
}; |
| 3826 |
|
if ( IsDebug() ) printf(" No file to be validated (force mode)! \n"); |
| 3827 |
|
return(0); |
| 3828 |
|
}; |
| 3829 |
|
myid=(UInt_t)atoll(row->GetField(0)); |
| 3830 |
|
// |
| 3831 |
|
myquery.str(""); |
| 3832 |
|
myquery << " SELECT MAX(RUNTRAILER_TIME),MIN(RUNHEADER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< myid <<";"; |
| 3833 |
|
// |
| 3834 |
|
result = conn->Query(myquery.str().c_str()); |
| 3835 |
|
// |
| 3836 |
|
row = result->Next(); |
| 3837 |
|
if( !row->GetField(0) || !row->GetField(1)){ |
| 3838 |
|
// |
| 3839 |
|
if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n"); |
| 3840 |
|
// |
| 3841 |
|
return(0); |
| 3842 |
|
// |
| 3843 |
|
} else { |
| 3844 |
|
// |
| 3845 |
|
UInt_t runhtime = (UInt_t)atoll(row->GetField(0)); |
| 3846 |
|
UInt_t runttime = (UInt_t)atoll(row->GetField(1)); |
| 3847 |
|
UInt_t caltime = 0; |
| 3848 |
|
// |
| 3849 |
|
myquery.str(""); |
| 3850 |
|
myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime; |
| 3851 |
|
myquery << " order by FROM_TIME asc limit 1;"; |
| 3852 |
|
// |
| 3853 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 3854 |
|
// |
| 3855 |
|
// |
| 3856 |
|
result = conn->Query(myquery.str().c_str()); |
| 3857 |
|
// |
| 3858 |
|
row = result->Next(); |
| 3859 |
|
if( !row ){ |
| 3860 |
|
caltime = runhtime; |
| 3861 |
|
} else { |
| 3862 |
|
caltime = (UInt_t)atoll(row->GetField(0)); |
| 3863 |
|
}; |
| 3864 |
|
// |
| 3865 |
|
myquery.str(""); |
| 3866 |
|
myquery << " SELECT * from GL_RUN where RUNHEADER_TIME>="<< runttime <<" AND RUNHEADER_TIME<=" ; |
| 3867 |
|
myquery << caltime << " order by RUNHEADER_TIME DESC"; |
| 3868 |
|
// |
| 3869 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 3870 |
|
// |
| 3871 |
|
result = conn->Query(myquery.str().c_str()); |
| 3872 |
|
// |
| 3873 |
|
}; |
| 3874 |
}; |
}; |
| 3875 |
if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start); |
// |
|
// -------------------------------------------------------------- |
|
|
// now retrieves runs to be validated |
|
|
// -------------------------------------------------------------- |
|
|
oss.str(""); |
|
|
oss << " SELECT * FROM GL_RUN WHERE RUNHEADER_TIME <=" << t_start; |
|
|
oss << " AND RUNHEADER_TIME >="<< t_stop; |
|
|
oss << " ORDER BY RUNHEADER_TIME DESC;"; |
|
|
// if ( IsDebug() ) |
|
|
if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str()); |
|
|
result = conn->Query(oss.str().c_str()); |
|
| 3876 |
if ( !result ) throw -4; |
if ( !result ) throw -4; |
| 3877 |
if ( !result->GetRowCount() ) printf(" No runs to validate \n"); |
if ( !result->GetRowCount() && IsDebug() ) printf(" No runs to validate \n"); |
| 3878 |
// printf("------------------------------------------------------------------------------- \n"); |
// |
|
|
|
| 3879 |
Int_t nrow = 0; |
Int_t nrow = 0; |
| 3880 |
GL_RUN* this_run = new GL_RUN(); |
GL_RUN* this_run = new GL_RUN(); |
| 3881 |
GL_RUN* next_run = new GL_RUN(); |
GL_RUN* next_run = new GL_RUN(); |
| 3889 |
UInt_t t1=0,t2=0; |
UInt_t t1=0,t2=0; |
| 3890 |
// --------------------------------------------------------------------------------- |
// --------------------------------------------------------------------------------- |
| 3891 |
// - loop over runs, back in time, |
// - loop over runs, back in time, |
| 3892 |
// - select sequences of runs close in time (less than 60 s apart), |
// - select sequences of runs close in time (less than calibtime s apart), |
| 3893 |
// which could be preceeded by a calibration |
// which could be preceeded by a calibration |
| 3894 |
// - check if there might be a missing calibration |
// - check if there might be a missing calibration |
| 3895 |
// --------------------------------------------------------------------------------- |
// --------------------------------------------------------------------------------- |
| 3929 |
|
|
| 3930 |
if( this_run->ID == next_run->ID_RUN_FRAG ) interval = 0; //=> run fragments |
if( this_run->ID == next_run->ID_RUN_FRAG ) interval = 0; //=> run fragments |
| 3931 |
|
|
| 3932 |
if( interval >= 60 )CHECK = true; //more than 60 s => there might be a calibration |
if( interval >= calibtime )CHECK = true; //more than calibtime s => there might be a calibration |
| 3933 |
|
|
| 3934 |
if( !CHECK && this_run->VALIDATION ){ |
if( !CHECK && this_run->VALIDATION ){ |
| 3935 |
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true); |
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true); |
| 3969 |
}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); |
| 3970 |
}; |
}; |
| 3971 |
|
|
| 3972 |
if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID); |
if ( IsDebug() ) printf("%i Run %u \n",nrow,this_run->ID); |
| 3973 |
nrow++; |
nrow++; |
| 3974 |
|
|
| 3975 |
}; |
}; |
| 3998 |
// which should be equal to the time between ascending-nodes. |
// which should be equal to the time between ascending-nodes. |
| 3999 |
//============================================================== |
//============================================================== |
| 4000 |
if ( t2 - trkcalib->FROM_TIME > 5700) { |
if ( t2 - trkcalib->FROM_TIME > 5700) { |
| 4001 |
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); |
| 4002 |
//============================================================== |
//============================================================== |
| 4003 |
// there might be a missing calibration, due to: |
// there might be a missing calibration, due to: |
| 4004 |
// - MM full |
// - MM full |
| 4019 |
// it is enough to say that there are no missing calibrations |
// it is enough to say that there are no missing calibrations |
| 4020 |
//============================================================== |
//============================================================== |
| 4021 |
// the long time interval bewteen runs might be due to download |
// the long time interval bewteen runs might be due to download |
| 4022 |
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); |
| 4023 |
return(false); |
return(false); |
| 4024 |
|
|
| 4025 |
}; |
}; |
| 4127 |
|
|
| 4128 |
|
|
| 4129 |
// Insert tle in the table GL_TLE using the connection conn. |
// Insert tle in the table GL_TLE using the connection conn. |
| 4130 |
int PamelaDBOperations::insertTle(cTle *tle) |
Int_t PamelaDBOperations::insertTle(cTle *tle) |
| 4131 |
{ |
{ |
| 4132 |
stringstream oss; |
stringstream oss; |
| 4133 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
| 4210 |
|
|
| 4211 |
return date.str(); |
return date.str(); |
| 4212 |
} |
} |
| 4213 |
|
|
| 4214 |
|
/** |
| 4215 |
|
* Remove a file from the DB, delete on cascade all entries related to that file |
| 4216 |
|
* rearrange GL_RUN and GL_XXX_CALIB tables, turn off validation till the following |
| 4217 |
|
* calibration |
| 4218 |
|
**/ |
| 4219 |
|
Int_t PamelaDBOperations::removeFile(TString remfile){ |
| 4220 |
|
// |
| 4221 |
|
// Determine ID_ROOT_L0 and ID_RAW |
| 4222 |
|
// |
| 4223 |
|
TSQLResult *pResult; |
| 4224 |
|
TSQLRow *Row; |
| 4225 |
|
stringstream myquery; |
| 4226 |
|
// |
| 4227 |
|
myquery.str(""); |
| 4228 |
|
myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';"; |
| 4229 |
|
// |
| 4230 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4231 |
|
// |
| 4232 |
|
Row = pResult->Next(); |
| 4233 |
|
if( !Row ){ |
| 4234 |
|
if ( strcmp(remfile.Data(),GetRootName().Data()) ){ |
| 4235 |
|
if ( IsDebug() ) printf(" No file to be removed even if option \"-remove file\" was used!!\n"); |
| 4236 |
|
return(1); |
| 4237 |
|
}; |
| 4238 |
|
if ( IsDebug() ) printf(" No file to be removed (force mode)! \n"); |
| 4239 |
|
return(0); |
| 4240 |
|
}; |
| 4241 |
|
// |
| 4242 |
|
this->SetID_ROOT((UInt_t)atoll(Row->GetField(0))); |
| 4243 |
|
this->SetID_RAW((UInt_t)atoll(Row->GetField(1))); |
| 4244 |
|
// |
| 4245 |
|
this->ValidationOFF(); |
| 4246 |
|
// |
| 4247 |
|
this->RemoveCALIBS(); |
| 4248 |
|
// |
| 4249 |
|
this->RemoveRUNS(); |
| 4250 |
|
// |
| 4251 |
|
this->RemoveFILES(); |
| 4252 |
|
// |
| 4253 |
|
this->SetID_ROOT(0); |
| 4254 |
|
this->SetID_RAW(0); |
| 4255 |
|
// |
| 4256 |
|
return(0); |
| 4257 |
|
}; |
| 4258 |
|
|
| 4259 |
|
/** |
| 4260 |
|
* |
| 4261 |
|
* Set validation bit to zero for runs following the removing file till |
| 4262 |
|
* 1) a run with TRK_CALIB_USED=140 |
| 4263 |
|
* 2) a run with VALIDATION = 0 |
| 4264 |
|
* 3) the next calibration |
| 4265 |
|
* |
| 4266 |
|
**/ |
| 4267 |
|
void PamelaDBOperations::ValidationOFF(){ |
| 4268 |
|
TSQLResult *pResult; |
| 4269 |
|
TSQLRow *Row; |
| 4270 |
|
stringstream myquery; |
| 4271 |
|
Int_t unv = 0; |
| 4272 |
|
//select ID from GL_RUN where RUNHEADER_TIME>=1152671382 AND (VALIDATION=0 OR TRK_CALIB_USED=104) order by RUNHEADER_TIME asc limit 1; |
| 4273 |
|
myquery.str(""); |
| 4274 |
|
myquery << " SELECT MAX(RUNTRAILER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";"; |
| 4275 |
|
// |
| 4276 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4277 |
|
// |
| 4278 |
|
Row = pResult->Next(); |
| 4279 |
|
if( !Row->GetField(0) ){ |
| 4280 |
|
// |
| 4281 |
|
if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n"); |
| 4282 |
|
// |
| 4283 |
|
} else { |
| 4284 |
|
// |
| 4285 |
|
UInt_t runhtime = (UInt_t)atoll(Row->GetField(0)); |
| 4286 |
|
UInt_t caltime = 0; |
| 4287 |
|
// |
| 4288 |
|
myquery.str(""); |
| 4289 |
|
myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime; |
| 4290 |
|
myquery << " order by FROM_TIME asc limit 1;"; |
| 4291 |
|
// |
| 4292 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4293 |
|
// |
| 4294 |
|
// |
| 4295 |
|
delete pResult; |
| 4296 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4297 |
|
// |
| 4298 |
|
Row = pResult->Next(); |
| 4299 |
|
if( !Row ){ |
| 4300 |
|
caltime = runhtime; |
| 4301 |
|
} else { |
| 4302 |
|
caltime = (UInt_t)atoll(Row->GetField(0)); |
| 4303 |
|
}; |
| 4304 |
|
// |
| 4305 |
|
myquery.str(""); |
| 4306 |
|
myquery << " SELECT ID,RUNHEADER_TIME from GL_RUN where RUNHEADER_TIME>="<< runhtime <<" AND (VALIDATION=0 OR TRK_CALIB_USED=104 OR RUNHEADER_TIME>" ; |
| 4307 |
|
myquery << caltime << ") order by RUNHEADER_TIME asc LIMIT 1"; |
| 4308 |
|
// |
| 4309 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4310 |
|
// |
| 4311 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4312 |
|
// |
| 4313 |
|
Row = pResult->Next(); |
| 4314 |
|
if( !Row ){ |
| 4315 |
|
// |
| 4316 |
|
if ( IsDebug() ) printf(" NO RUN NEED TO BE UNVALIDATED \n"); |
| 4317 |
|
// |
| 4318 |
|
} else { |
| 4319 |
|
myquery.str(""); |
| 4320 |
|
myquery << " SELECT ID from GL_RUN where RUNHEADER_TIME<"<< Row->GetField(1) <<" AND "; |
| 4321 |
|
myquery << " RUNHEADER_TIME>=" <<runhtime; |
| 4322 |
|
myquery << " order by RUNHEADER_TIME asc;"; |
| 4323 |
|
// |
| 4324 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4325 |
|
// |
| 4326 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4327 |
|
// |
| 4328 |
|
Row = pResult->Next(); |
| 4329 |
|
while ( Row ){ |
| 4330 |
|
// |
| 4331 |
|
unv++; |
| 4332 |
|
this->assignVALIDATION((UInt_t)atoll(Row->GetField(0)), false); |
| 4333 |
|
Row = pResult->Next(); |
| 4334 |
|
// |
| 4335 |
|
}; |
| 4336 |
|
}; |
| 4337 |
|
}; |
| 4338 |
|
if ( IsDebug() ) printf(" %u runs have been unvalidated \n",unv); |
| 4339 |
|
}; |
| 4340 |
|
|
| 4341 |
|
/** |
| 4342 |
|
* |
| 4343 |
|
* Rearrange GL_RUN table and remove runs |
| 4344 |
|
* |
| 4345 |
|
**/ |
| 4346 |
|
void PamelaDBOperations::RemoveRUNS(){ |
| 4347 |
|
TSQLResult *pResult; |
| 4348 |
|
TSQLRow *Row; |
| 4349 |
|
stringstream myquery; |
| 4350 |
|
UInt_t drun = 0; |
| 4351 |
|
GL_RUN *delrun = new GL_RUN(); |
| 4352 |
|
// |
| 4353 |
|
myquery.str(""); |
| 4354 |
|
myquery << " SELECT ID FROM GL_RUN where ID_RUN_FRAG=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";"; |
| 4355 |
|
// |
| 4356 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4357 |
|
// |
| 4358 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4359 |
|
// |
| 4360 |
|
Row = pResult->Next(); |
| 4361 |
|
// |
| 4362 |
|
// |
| 4363 |
|
if ( !Row ){ |
| 4364 |
|
if ( IsDebug() ) printf(" No run with ID_RUN_FRAG=0 belonged to this file \n"); |
| 4365 |
|
} else { |
| 4366 |
|
if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n"); |
| 4367 |
|
while ( Row ){ |
| 4368 |
|
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
| 4369 |
|
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4370 |
|
drun++; |
| 4371 |
|
Row = pResult->Next(); |
| 4372 |
|
}; |
| 4373 |
|
}; |
| 4374 |
|
// |
| 4375 |
|
// |
| 4376 |
|
myquery.str(""); |
| 4377 |
|
myquery << " SELECT ID,ID_RUN_FRAG FROM GL_RUN where ID_RUN_FRAG!=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";"; |
| 4378 |
|
// |
| 4379 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4380 |
|
// |
| 4381 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4382 |
|
// |
| 4383 |
|
Row = pResult->Next(); |
| 4384 |
|
// |
| 4385 |
|
if ( !Row ){ |
| 4386 |
|
if ( IsDebug() ) printf(" No run with ID_RUN_FRAG!=0 belonged to this file \n"); |
| 4387 |
|
} else { |
| 4388 |
|
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n"); |
| 4389 |
|
while ( Row ){ |
| 4390 |
|
if ( IsDebug() ) printf(" restore run %u \n",(UInt_t)atoll(Row->GetField(1))); |
| 4391 |
|
delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS"); |
| 4392 |
|
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(1))); |
| 4393 |
|
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN"); |
| 4394 |
|
if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){ |
| 4395 |
|
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4396 |
|
delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN"); |
| 4397 |
|
}; |
| 4398 |
|
drun++; |
| 4399 |
|
Row = pResult->Next(); |
| 4400 |
|
}; |
| 4401 |
|
}; |
| 4402 |
|
// |
| 4403 |
|
if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN table \n",drun); |
| 4404 |
|
// |
| 4405 |
|
// |
| 4406 |
|
// |
| 4407 |
|
drun = 0; |
| 4408 |
|
// |
| 4409 |
|
myquery.str(""); |
| 4410 |
|
myquery << " SELECT ID_TRASH FROM GL_RUN_TRASH where BELONGED_TO='GL_RUN_FRAGMENTS' AND ID_ROOT_L0=" <<this->GetID_ROOT() <<";"; |
| 4411 |
|
// |
| 4412 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4413 |
|
// |
| 4414 |
|
Row = pResult->Next(); |
| 4415 |
|
// |
| 4416 |
|
if ( !Row ){ |
| 4417 |
|
if ( IsDebug() ) printf(" No run from GL_RUN_FRAGMENTS table in the trash table for this file \n"); |
| 4418 |
|
} else { |
| 4419 |
|
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n"); |
| 4420 |
|
while ( Row ){ |
| 4421 |
|
if ( IsDebug() ) printf(" del run idtrash %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4422 |
|
myquery.str(""); |
| 4423 |
|
myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";"; |
| 4424 |
|
conn->Query(myquery.str().c_str()); |
| 4425 |
|
drun++; |
| 4426 |
|
Row = pResult->Next(); |
| 4427 |
|
}; |
| 4428 |
|
}; |
| 4429 |
|
// |
| 4430 |
|
if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_TRASH table \n",drun); |
| 4431 |
|
// |
| 4432 |
|
// |
| 4433 |
|
// |
| 4434 |
|
drun = 0; |
| 4435 |
|
// |
| 4436 |
|
myquery.str(""); |
| 4437 |
|
myquery << " SELECT ID FROM GL_RUN_FRAGMENTS where ID_ROOT_L0=" <<this->GetID_ROOT() <<";"; |
| 4438 |
|
// |
| 4439 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4440 |
|
// |
| 4441 |
|
Row = pResult->Next(); |
| 4442 |
|
// |
| 4443 |
|
if ( !Row ){ |
| 4444 |
|
if ( IsDebug() ) printf(" No run in the GL_RUN_FRAGMENTS table for this file \n"); |
| 4445 |
|
} else { |
| 4446 |
|
if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n"); |
| 4447 |
|
while ( Row ){ |
| 4448 |
|
if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0))); |
| 4449 |
|
myquery.str(""); |
| 4450 |
|
myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";"; |
| 4451 |
|
conn->Query(myquery.str().c_str()); |
| 4452 |
|
drun++; |
| 4453 |
|
Row = pResult->Next(); |
| 4454 |
|
}; |
| 4455 |
|
}; |
| 4456 |
|
// |
| 4457 |
|
if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_FRAGMENTS table \n",drun); |
| 4458 |
|
// |
| 4459 |
|
// |
| 4460 |
|
// |
| 4461 |
|
delete delrun; |
| 4462 |
|
// |
| 4463 |
|
}; |
| 4464 |
|
|
| 4465 |
|
|
| 4466 |
|
/** |
| 4467 |
|
* |
| 4468 |
|
* Rearrange calibration tables |
| 4469 |
|
* |
| 4470 |
|
**/ |
| 4471 |
|
void PamelaDBOperations::RemoveFILES(){ |
| 4472 |
|
stringstream myquery; |
| 4473 |
|
// |
| 4474 |
|
myquery.str(""); |
| 4475 |
|
myquery << " DELETE FROM GL_RAW WHERE ID=" <<this->GetID_RAW() <<";"; |
| 4476 |
|
// |
| 4477 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
| 4478 |
|
// |
| 4479 |
|
conn->Query(myquery.str().c_str()); |
| 4480 |
|
// |
| 4481 |
|
}; |
| 4482 |
|
|
| 4483 |
|
/** |
| 4484 |
|
* |
| 4485 |
|
* Rearrange calibration tables |
| 4486 |
|
* |
| 4487 |
|
**/ |
| 4488 |
|
void PamelaDBOperations::RemoveCALIBS(){ |
| 4489 |
|
TSQLResult *pResult; |
| 4490 |
|
TSQLRow *Row; |
| 4491 |
|
stringstream myquery; |
| 4492 |
|
// |
| 4493 |
|
// |
| 4494 |
|
// Calorimeter |
| 4495 |
|
// |
| 4496 |
|
for (Int_t section = 0; section < 4; section++){ |
| 4497 |
|
myquery.str(""); |
| 4498 |
|
myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_CALO_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<" AND "; |
| 4499 |
|
myquery << " SECTION=" << section << ";"; |
| 4500 |
|
// |
| 4501 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4502 |
|
// |
| 4503 |
|
Row = pResult->Next(); |
| 4504 |
|
if( !Row->GetField(0) || !Row->GetField(1) ){ |
| 4505 |
|
// |
| 4506 |
|
if ( IsDebug() ) printf(" NO CALO CALIBRATION SECTION %i ASSOCIATED TO THIS FILE! \n",section); |
| 4507 |
|
// |
| 4508 |
|
} else { |
| 4509 |
|
// |
| 4510 |
|
myquery.str(""); |
| 4511 |
|
myquery << " UPDATE GL_CALO_CALIB SET TO_TIME=" << Row->GetField(1); |
| 4512 |
|
myquery << " WHERE TO_TIME="<< Row->GetField(0) << " AND "; |
| 4513 |
|
myquery << " SECTION=" << section << ";"; |
| 4514 |
|
// |
| 4515 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4516 |
|
// |
| 4517 |
|
if( !pResult ){ |
| 4518 |
|
// |
| 4519 |
|
if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n"); |
| 4520 |
|
// |
| 4521 |
|
throw -4; |
| 4522 |
|
// |
| 4523 |
|
}; |
| 4524 |
|
// |
| 4525 |
|
}; |
| 4526 |
|
}; |
| 4527 |
|
myquery.str(""); |
| 4528 |
|
myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT(); |
| 4529 |
|
// |
| 4530 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4531 |
|
// |
| 4532 |
|
if( !pResult ){ |
| 4533 |
|
// |
| 4534 |
|
if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n"); |
| 4535 |
|
// |
| 4536 |
|
throw -4; |
| 4537 |
|
// |
| 4538 |
|
}; |
| 4539 |
|
// |
| 4540 |
|
// Tracker |
| 4541 |
|
// |
| 4542 |
|
myquery.str(""); |
| 4543 |
|
myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_TRK_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";"; |
| 4544 |
|
// |
| 4545 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4546 |
|
// |
| 4547 |
|
Row = pResult->Next(); |
| 4548 |
|
if( !Row->GetField(0) || !Row->GetField(1) ){ |
| 4549 |
|
// |
| 4550 |
|
if ( IsDebug() ) printf(" NO TRK CALIBRATION ASSOCIATED TO THIS FILE! \n"); |
| 4551 |
|
// |
| 4552 |
|
} else { |
| 4553 |
|
// |
| 4554 |
|
myquery.str(""); |
| 4555 |
|
myquery << " UPDATE GL_TRK_CALIB SET TO_TIME=" << Row->GetField(1); |
| 4556 |
|
myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";"; |
| 4557 |
|
// |
| 4558 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4559 |
|
// |
| 4560 |
|
if( !pResult ){ |
| 4561 |
|
// |
| 4562 |
|
if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n"); |
| 4563 |
|
// |
| 4564 |
|
throw -4; |
| 4565 |
|
// |
| 4566 |
|
}; |
| 4567 |
|
// |
| 4568 |
|
myquery.str(""); |
| 4569 |
|
myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT(); |
| 4570 |
|
// |
| 4571 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4572 |
|
// |
| 4573 |
|
if( !pResult ){ |
| 4574 |
|
// |
| 4575 |
|
if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n"); |
| 4576 |
|
// |
| 4577 |
|
throw -4; |
| 4578 |
|
// |
| 4579 |
|
}; |
| 4580 |
|
}; |
| 4581 |
|
// |
| 4582 |
|
// |
| 4583 |
|
// S4 |
| 4584 |
|
// |
| 4585 |
|
myquery.str(""); |
| 4586 |
|
myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_S4_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";"; |
| 4587 |
|
// |
| 4588 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4589 |
|
// |
| 4590 |
|
Row = pResult->Next(); |
| 4591 |
|
if( !Row->GetField(0) || !Row->GetField(1) ){ |
| 4592 |
|
// |
| 4593 |
|
if ( IsDebug() ) printf(" NO S4 CALIBRATION ASSOCIATED TO THIS FILE! \n"); |
| 4594 |
|
// |
| 4595 |
|
} else { |
| 4596 |
|
// |
| 4597 |
|
myquery.str(""); |
| 4598 |
|
myquery << " UPDATE GL_S4_CALIB SET TO_TIME=" << Row->GetField(1); |
| 4599 |
|
myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";"; |
| 4600 |
|
// |
| 4601 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4602 |
|
// |
| 4603 |
|
if( !pResult ){ |
| 4604 |
|
// |
| 4605 |
|
if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n"); |
| 4606 |
|
// |
| 4607 |
|
throw -4; |
| 4608 |
|
// |
| 4609 |
|
}; |
| 4610 |
|
// |
| 4611 |
|
myquery.str(""); |
| 4612 |
|
myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT(); |
| 4613 |
|
// |
| 4614 |
|
pResult = conn->Query(myquery.str().c_str()); |
| 4615 |
|
// |
| 4616 |
|
if( !pResult ){ |
| 4617 |
|
// |
| 4618 |
|
if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n"); |
| 4619 |
|
// |
| 4620 |
|
throw -4; |
| 4621 |
|
// |
| 4622 |
|
}; |
| 4623 |
|
// |
| 4624 |
|
}; |
| 4625 |
|
}; |
| 4626 |
|
|
| 4627 |
|
/** |
| 4628 |
|
* |
| 4629 |
|
* Rearrange calibration tables |
| 4630 |
|
* |
| 4631 |
|
**/ |
| 4632 |
|
UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){ |
| 4633 |
|
|
| 4634 |
|
Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0}; |
| 4635 |
|
UInt_t timeaftercalib=120000; //2000; |
| 4636 |
|
// ---------- |
| 4637 |
|
// Check CRCs |
| 4638 |
|
// ---------- |
| 4639 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
| 4640 |
|
if( caltrk->crc_hcal[ipkt] )return 0; // :-( |
| 4641 |
|
for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-( |
| 4642 |
|
} |
| 4643 |
|
// ----------------------- |
| 4644 |
|
// Check missing packets: |
| 4645 |
|
// ----------------------- |
| 4646 |
|
// Readout order: |
| 4647 |
|
// ------------------ |
| 4648 |
|
// DSP packet board |
| 4649 |
|
// ------------------ |
| 4650 |
|
// 12 0 1 |
| 4651 |
|
// 10 1 1 |
| 4652 |
|
// 8 2 1 |
| 4653 |
|
// 4 3 1 |
| 4654 |
|
// 6 4 1 |
| 4655 |
|
// 2 5 1 |
| 4656 |
|
// ------------------ |
| 4657 |
|
// 11 0 2 |
| 4658 |
|
// 9 1 2 |
| 4659 |
|
// 7 2 2 |
| 4660 |
|
// 3 3 2 |
| 4661 |
|
// 5 4 2 |
| 4662 |
|
// 1 5 2 |
| 4663 |
|
// ------------------ |
| 4664 |
|
// ------------------------------------------------- |
| 4665 |
|
// Check if it is first or second calibration packet |
| 4666 |
|
// ------------------------------------------------- |
| 4667 |
|
UInt_t build=0; |
| 4668 |
|
TString classname = caltrk->GetName(); |
| 4669 |
|
UInt_t base=0; |
| 4670 |
|
UInt_t mask=0; |
| 4671 |
|
if(classname.Contains("CalibTrk1Event")){ |
| 4672 |
|
base=12; |
| 4673 |
|
mask=0x03F000; |
| 4674 |
|
} |
| 4675 |
|
if(classname.Contains("CalibTrk2Event")){ |
| 4676 |
|
base=18; |
| 4677 |
|
mask=0xFC0000; |
| 4678 |
|
} |
| 4679 |
|
// ------------------------------------------------- |
| 4680 |
|
// Count number of packets and set build variable |
| 4681 |
|
// ------------------------------------------------- |
| 4682 |
|
Int_t npkts=0; |
| 4683 |
|
for(Int_t ipkt=0; ipkt<6; ipkt++){ |
| 4684 |
|
if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){ |
| 4685 |
|
npkts++; |
| 4686 |
|
build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) ); |
| 4687 |
|
} |
| 4688 |
|
} |
| 4689 |
|
// if( npkts==6 )return 1; // :-) |
| 4690 |
|
|
| 4691 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl; |
| 4692 |
|
|
| 4693 |
|
// ----------------------------------------------- |
| 4694 |
|
// If missing packets: check the acq configuration |
| 4695 |
|
// (some DSPs might be excluded from acquisition) |
| 4696 |
|
// ----------------------------------------------- |
| 4697 |
|
|
| 4698 |
|
// ----------------------------------------------- |
| 4699 |
|
// retrieve the first run header after calib |
| 4700 |
|
// ----------------------------------------------- |
| 4701 |
|
PacketType *pctp; |
| 4702 |
|
EventCounter *cod; |
| 4703 |
|
cod = eh->GetCounter(); |
| 4704 |
|
Int_t irun = cod->Get(pctp->RunHeader); |
| 4705 |
|
TTree *rh=(TTree*)file->Get("RunHeader"); |
| 4706 |
|
if ( !rh || rh->IsZombie() ) throw -17; |
| 4707 |
|
if( rh->GetEntries() == irun ){ |
| 4708 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1) -- cannot validate :-( "<<endl; |
| 4709 |
|
return 0; // :-( |
| 4710 |
|
} |
| 4711 |
|
|
| 4712 |
|
RunHeaderEvent *run = 0; |
| 4713 |
|
EventHeader *hrun = 0; |
| 4714 |
|
rh->SetBranchAddress("RunHeader", &run); |
| 4715 |
|
rh->SetBranchAddress("Header", &hrun); |
| 4716 |
|
rh->GetEntry(irun); |
| 4717 |
|
// cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl; |
| 4718 |
|
|
| 4719 |
|
if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){ |
| 4720 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl; |
| 4721 |
|
return 0; // :-( |
| 4722 |
|
} |
| 4723 |
|
|
| 4724 |
|
if( !run->RM_ACQ_AFTER_CALIB ){ |
| 4725 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0 -- cannot validate :-( "<<endl; |
| 4726 |
|
return 0; // :-( |
| 4727 |
|
} |
| 4728 |
|
|
| 4729 |
|
UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime()); |
| 4730 |
|
if( dtime > timeaftercalib ){ |
| 4731 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl; |
| 4732 |
|
return 0; // :-( |
| 4733 |
|
} |
| 4734 |
|
|
| 4735 |
|
|
| 4736 |
|
|
| 4737 |
|
if( (run->ACQ_BUILD_INFO & mask) != build ){ |
| 4738 |
|
if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl; |
| 4739 |
|
return 0; // :-( |
| 4740 |
|
} |
| 4741 |
|
return 1; // :-) |
| 4742 |
|
|
| 4743 |
|
} |
| 4744 |
|
|
| 4745 |
|
/** |
| 4746 |
|
* |
| 4747 |
|
* Check the DB (only for overlapping runs at the moment) |
| 4748 |
|
* |
| 4749 |
|
**/ |
| 4750 |
|
UInt_t PamelaDBOperations::Check(){ |
| 4751 |
|
// |
| 4752 |
|
UInt_t test = 0; |
| 4753 |
|
// |
| 4754 |
|
UInt_t thisrht = 0; |
| 4755 |
|
UInt_t thisrtt = 0; |
| 4756 |
|
UInt_t thisid = 0; |
| 4757 |
|
UInt_t prevrht = 0; |
| 4758 |
|
UInt_t prevrtt = 0; |
| 4759 |
|
UInt_t previd = 0; |
| 4760 |
|
// |
| 4761 |
|
UInt_t prevl0id = 0; |
| 4762 |
|
UInt_t thisl0id = 0; |
| 4763 |
|
// |
| 4764 |
|
stringstream oss; |
| 4765 |
|
TSQLResult *result = 0; |
| 4766 |
|
TSQLRow *row = 0; |
| 4767 |
|
TSQLResult *result2 = 0; |
| 4768 |
|
TSQLRow *row2 = 0; |
| 4769 |
|
TSQLResult *result3 = 0; |
| 4770 |
|
TSQLRow *row3 = 0; |
| 4771 |
|
oss.str(""); |
| 4772 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME,NEVENTS FROM GL_RUN order by RUNHEADER_TIME asc;"; |
| 4773 |
|
// oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;"; |
| 4774 |
|
result = conn->Query(oss.str().c_str()); |
| 4775 |
|
// |
| 4776 |
|
if ( !result ) throw -4;; |
| 4777 |
|
// |
| 4778 |
|
row = result->Next(); |
| 4779 |
|
UInt_t nid = 0; |
| 4780 |
|
// |
| 4781 |
|
while ( row ){ |
| 4782 |
|
nid++; |
| 4783 |
|
if ( !(nid%1000) && nid ) printf(" %iK run scanned \n",nid/1000); |
| 4784 |
|
thisid = (UInt_t)atoll(row->GetField(0)); |
| 4785 |
|
thisl0id = (UInt_t)atoll(row->GetField(1)); |
| 4786 |
|
thisrht = (UInt_t)atoll(row->GetField(2)); |
| 4787 |
|
thisrtt = (UInt_t)atoll(row->GetField(3)); |
| 4788 |
|
// |
| 4789 |
|
if ( (UInt_t)atoll(row->GetField(4)) > 1 ){ |
| 4790 |
|
// |
| 4791 |
|
// |
| 4792 |
|
// |
| 4793 |
|
oss.str(""); |
| 4794 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!=" |
| 4795 |
|
<< thisid << " AND ( RUNHEADER_TIME=" |
| 4796 |
|
<< thisrht << " OR RUNTRAILER_TIME=" |
| 4797 |
|
<< thisrtt << " ) AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;"; |
| 4798 |
|
result3 = conn->Query(oss.str().c_str()); |
| 4799 |
|
if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str()); |
| 4800 |
|
if ( result3 ){ |
| 4801 |
|
// |
| 4802 |
|
oss.str(""); |
| 4803 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!=" |
| 4804 |
|
<< thisid << " AND RUNHEADER_TIME=" |
| 4805 |
|
<< thisrht << " AND RUNTRAILER_TIME!=" |
| 4806 |
|
<< thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;"; |
| 4807 |
|
result3 = conn->Query(oss.str().c_str()); |
| 4808 |
|
if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str()); |
| 4809 |
|
if ( result3 ){ |
| 4810 |
|
row3 = result3->Next(); |
| 4811 |
|
// |
| 4812 |
|
while ( row3 ){ |
| 4813 |
|
// |
| 4814 |
|
// 2 runs with same runheader |
| 4815 |
|
// |
| 4816 |
|
printf(" CHECK n.4 RUNs %u and %u HAVE SAME RUNHEADER \n",thisid,(UInt_t)atoll(row3->GetField(0))); |
| 4817 |
|
row3 = result3->Next(); |
| 4818 |
|
}; |
| 4819 |
|
}; |
| 4820 |
|
// |
| 4821 |
|
oss.str(""); |
| 4822 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!=" |
| 4823 |
|
<< thisid << " AND RUNHEADER_TIME!=" |
| 4824 |
|
<< thisrht << " AND RUNTRAILER_TIME=" |
| 4825 |
|
<< thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;"; |
| 4826 |
|
result3 = conn->Query(oss.str().c_str()); |
| 4827 |
|
if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str()); |
| 4828 |
|
if ( result3 ){ |
| 4829 |
|
row3 = result3->Next(); |
| 4830 |
|
// |
| 4831 |
|
while ( row3 ){ |
| 4832 |
|
// |
| 4833 |
|
// 2 runs with same runtrailer |
| 4834 |
|
// |
| 4835 |
|
printf(" CHECK n.5 RUNs %u and %u HAVE SAME RUNTRAILER \n",thisid,(UInt_t)atoll(row3->GetField(0))); |
| 4836 |
|
row3 = result3->Next(); |
| 4837 |
|
}; |
| 4838 |
|
}; |
| 4839 |
|
// |
| 4840 |
|
oss.str(""); |
| 4841 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!=" |
| 4842 |
|
<< thisid << " AND RUNHEADER_TIME=" |
| 4843 |
|
<< thisrht << " AND RUNTRAILER_TIME=" |
| 4844 |
|
<< thisrtt << " AND ID_RUN_FRAG!=" |
| 4845 |
|
<< thisid << " order by RUNHEADER_TIME asc;"; |
| 4846 |
|
result3 = conn->Query(oss.str().c_str()); |
| 4847 |
|
if ( result3 ){ |
| 4848 |
|
row3 = result3->Next(); |
| 4849 |
|
// |
| 4850 |
|
while ( row3 ){ |
| 4851 |
|
// |
| 4852 |
|
// duplicated run |
| 4853 |
|
// |
| 4854 |
|
printf(" CHECK n.7 RUNs %u and %u HAVE SAME RUNTRAILER AND RUNHEADER (ARE THE SAME?) \n",thisid,(UInt_t)atoll(row3->GetField(0))); |
| 4855 |
|
row3 = result3->Next(); |
| 4856 |
|
}; |
| 4857 |
|
}; |
| 4858 |
|
}; |
| 4859 |
|
// |
| 4860 |
|
oss.str(""); |
| 4861 |
|
oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!=" |
| 4862 |
|
<< thisid << " AND RUNHEADER_TIME>" |
| 4863 |
|
<< thisrht << " AND RUNTRAILER_TIME<" |
| 4864 |
|
<< thisrtt << " order by RUNHEADER_TIME asc;"; |
| 4865 |
|
result3 = conn->Query(oss.str().c_str()); |
| 4866 |
|
if ( result3 ){ |
| 4867 |
|
row3 = result3->Next(); |
| 4868 |
|
// |
| 4869 |
|
while ( row3 ){ |
| 4870 |
|
// |
| 4871 |
|
// run contained in the checked one |
| 4872 |
|
// |
| 4873 |
|
printf(" CHECK n.6 RUN %u CONTAINS RUN %u \n",thisid,(UInt_t)atoll(row3->GetField(0))); |
| 4874 |
|
row3 = result3->Next(); |
| 4875 |
|
}; |
| 4876 |
|
}; |
| 4877 |
|
// |
| 4878 |
|
}; |
| 4879 |
|
// |
| 4880 |
|
// if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){ |
| 4881 |
|
// if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){ |
| 4882 |
|
if ( (thisrht < prevrtt) && (thisrht != prevrht) ){ |
| 4883 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
| 4884 |
|
printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid); |
| 4885 |
|
TString prevf = ""; |
| 4886 |
|
TString thisf = ""; |
| 4887 |
|
oss.str(""); |
| 4888 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
| 4889 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4890 |
|
if ( !result2 ) throw -4;; |
| 4891 |
|
row2 = result2->Next(); |
| 4892 |
|
prevf = (TString)row2->GetField(0); |
| 4893 |
|
oss.str(""); |
| 4894 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
| 4895 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4896 |
|
if ( !result2 ) throw -4;; |
| 4897 |
|
row2 = result2->Next(); |
| 4898 |
|
thisf = (TString)row2->GetField(0); |
| 4899 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
| 4900 |
|
test = 1; |
| 4901 |
|
}; |
| 4902 |
|
// |
| 4903 |
|
if ( (thisrtt < prevrht) && (thisrht != prevrht) ){ |
| 4904 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
| 4905 |
|
printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid); |
| 4906 |
|
TString prevf = ""; |
| 4907 |
|
TString thisf = ""; |
| 4908 |
|
oss.str(""); |
| 4909 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
| 4910 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4911 |
|
if ( !result2 ) throw -4; |
| 4912 |
|
row2 = result2->Next(); |
| 4913 |
|
prevf = (TString)row2->GetField(0); |
| 4914 |
|
oss.str(""); |
| 4915 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
| 4916 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4917 |
|
if ( !result2 ) throw -4;; |
| 4918 |
|
row2 = result2->Next(); |
| 4919 |
|
thisf = (TString)row2->GetField(0); |
| 4920 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
| 4921 |
|
test = 1; |
| 4922 |
|
}; |
| 4923 |
|
// |
| 4924 |
|
if ( (thisrht > thisrtt) && (thisrht != prevrht) ){ |
| 4925 |
|
if ( IsDebug() ) printf(" IDprev %u ID %u prevrht %u prevrtt %u thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt); |
| 4926 |
|
printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid); |
| 4927 |
|
TString prevf = ""; |
| 4928 |
|
TString thisf = ""; |
| 4929 |
|
oss.str(""); |
| 4930 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";"; |
| 4931 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4932 |
|
if ( !result2 ) throw -4;; |
| 4933 |
|
row2 = result2->Next(); |
| 4934 |
|
prevf = (TString)row2->GetField(0); |
| 4935 |
|
oss.str(""); |
| 4936 |
|
oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";"; |
| 4937 |
|
result2 = conn->Query(oss.str().c_str()); |
| 4938 |
|
if ( !result2 ) throw -4;; |
| 4939 |
|
row2 = result2->Next(); |
| 4940 |
|
thisf = (TString)row2->GetField(0); |
| 4941 |
|
if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data()); |
| 4942 |
|
test = 1; |
| 4943 |
|
}; |
| 4944 |
|
|
| 4945 |
|
// |
| 4946 |
|
prevrht = thisrht; |
| 4947 |
|
prevrtt = thisrtt; |
| 4948 |
|
previd = thisid; |
| 4949 |
|
prevl0id = thisl0id; |
| 4950 |
|
row = result->Next(); |
| 4951 |
|
}; |
| 4952 |
|
// |
| 4953 |
|
return(test); |
| 4954 |
|
// |
| 4955 |
|
}; |