/[PAMELA software]/YodaProfiler/src/PamelaDBOperations.cpp

Diff of /YodaProfiler/src/PamelaDBOperations.cpp

Parent Directory | Revision Log | View Patch Patch

-revision 1.12 by mocchiut,
Tue Oct 17 15:20:04 2006 UTC
+revision 1.24 by mocchiut,
Mon Dec  4 11:32:07 2006 UTC
 Line 14
  #include <TSQLRow.h>
  #include <TTree.h>
  #include <TGraph.h>
- #include <TDatime.h>
+ #include <TTimeStamp.h>
  #include <TF1.h>
  //
  #include <EventHeader.h>
 Line 31
  #include <varDump/VarDumpRecord.h>
  #include <physics/S4/S4Event.h>
  //
- #include <cTle.h>
+ #include <sgp4.h>
- #include <cEci.h>
- #include <cJulian.h>
  #include <PamelaDBOperations.h>
  //
-Line 81 
 PamelaDBOperations::PamelaDBOperations(T
+Line 79 
 PamelaDBOperations::PamelaDBOperations(T
    INSERT_ROOT = !filerootname.IsNull();
    if( INSERT_ROOT ){
      this->SetRootName(filerootname);
+     this->SetOrbitNo();
      file = TFile::Open(this->GetRootName().Data());
+   } else {
+     this->SetRootName("");
    };
    //
    this->SetID_RAW(0);
-Line 111 
 void PamelaDBOperations::Close(){
+Line 112 
 void PamelaDBOperations::Close(){
  //
  void PamelaDBOperations::CheckValidate(Long64_t olderthan){
    clean_time = new TDatime();
+   //
    if(olderthan >= 0){
      VALIDATE = true;
      UInt_t timelim = 0;
-     timelim =  (UInt_t)clean_time->Convert() - olderthan;
+     timelim =  (UInt_t)clean_time->Convert(true) - olderthan;
      clean_time->Set(timelim,false);
    };
  };
-Line 126 
 void PamelaDBOperations::CheckValidate(L
+Line 128 
 void PamelaDBOperations::CheckValidate(L
   * @param password     password for the SQL connection.
   */
  void PamelaDBOperations::SetConnection(TString host, TString user, TString password){
+   if ( IsDebug() ) printf(" Connecting using host = %s user = %s password = %s \n",host.Data(),user.Data(),password.Data());
    conn = TSQLServer::Connect(host.Data(),user.Data(),password.Data());
  };
-Line 154 
 void PamelaDBOperations::SetDebugFlag(Bo
+Line 157 
 void PamelaDBOperations::SetDebugFlag(Bo
  };
  /**
+  * Set the autoboot flag
+  *
+  */
+ void PamelaDBOperations::SetAutoBoot(Bool_t dbg){
+   AUTOBOOT = dbg;
+ };
+ /**
+  * Set the nofrag flag
+  *
+  */
+ void PamelaDBOperations::SetNoFrag(Bool_t nf){
+   NOFRAG = nf;
+ };
+ /**
   * Store the BOOT number of the RAW file.
   * @param boot    BOOT number of the RAW file
   */
-Line 194 
 void PamelaDBOperations::SetRootName(TSt
+Line 213 
 void PamelaDBOperations::SetRootName(TSt
  };
  /**
+  * Store the downlink orbit number from filename.
+  */
+ void PamelaDBOperations::SetOrbitNo(){
+   dworbit = 0;
+   TString name = this->GetRootFile();
+   Int_t nlength = name.Length();
+   if ( nlength < 5 ) return;
+   TString dwo = 0;
+   for (Int_t i = 0; i<5; i++){
+     dwo.Append(name[i],1);
+   };
+   if ( dwo.IsDigit() ){
+     dworbit = (UInt_t)dwo.Atoi();
+   } else {
+     dwo="";
+     for (Int_t i = 8; i<13; i++){
+       dwo.Append(name[i],1);
+     };
+     if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
+   };
+   if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());
+   return;
+ };
+ /**
   * Store the NOBOOT flag.
   * @param noboot    true/false.
   */
-Line 257 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 303 
 Int_t PamelaDBOperations::SetUpperLimits
    UInt_t nevent = 0;
    UInt_t pktlast = 0;
    UInt_t obtlast = 0;
-   UInt_t t_pktlast = 0;
+   Long64_t t_pktlast = 0LL;
-   UInt_t t_obtlast = 0;
+   //  UInt_t t_obtlast = 0;
-   UInt_t upperpkt2 = 0;
+   Long64_t t_obtlast = 0LL;
-   ULong64_t upperobt2 = 0;
+   Long64_t upperpkt2 = 0LL;
+   Long64_t upperobt2 = 0LL;
    UInt_t zomp = 0;
    UInt_t jump = 50000; // was 5000
    EventCounter *code=0;
    //
+   Long64_t deltapkt = 5000LL;
+   Long64_t deltaobt = 50000LL;
+   //
    //   pcksList packetsNames;
    //   pcksList::iterator Iter;
    //   getPacketsNames(packetsNames);
-Line 301 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 351 
 Int_t PamelaDBOperations::SetUpperLimits
    upperobt = OBT(obtlast);
    upperentry = nevent-1;
    //
-   if ( IsDebug() ) printf(" First entries are: OBT %llu pkt_num %i \n",obtfirst,pktfirst);
+   if ( IsDebug() ) printf(" First entries are: OBT %i pkt_num %i \n",obtfirst,pktfirst);
    //
-   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);
    //
    if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(1);
    //
-Line 311 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 361 
 Int_t PamelaDBOperations::SetUpperLimits
    //
    if ( nevent < 2 ) return(4);
    //
-   if ( PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst) ){
+   if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) ){
+     //
      // go back
      zomp = nevent - 2;
      //
-Line 333 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 384 
 Int_t PamelaDBOperations::SetUpperLimits
          upperpkt2 = PKT(ph->GetCounter());
          upperobt2 = OBT(ph->GetOrbitalTime());
          //
-         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) ){
+           if ( IsDebug() ) printf(" .-. upperpkt2 %lld upperobt2 %lld \n",upperpkt2,upperobt2);
+           if ( IsDebug() ) printf(" .-. upperpkt %lld t_pktlast %lld upperobt %lld t_obtlast %lld \n",upperpkt,t_pktlast,upperobt,t_obtlast);
+           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);
+           throw -13;
+         };
          //
          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){
            zomp = i + jump + 1;
            if ( zomp > nevent-2 ) zomp = nevent - 2;
-           if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %i pktlast %i upperobt %llu obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
+           if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
            break;
          };
          //
-Line 373 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 429 
 Int_t PamelaDBOperations::SetUpperLimits
    //
    rhev = rh->GetEntries();
    rtev = rt->GetEntries();
-   UInt_t sobtt = 0;
+   Long64_t sobtt = 0LL;
-   UInt_t sobth = 0;
+   Long64_t sobth = 0LL;
-   UInt_t spktt = 0;
+   Long64_t spktt = 0LL;
-   UInt_t spkth = 0;
+   Long64_t spkth = 0LL;
-   UInt_t pktt = 0;
+   Long64_t pktt = 0LL;
-   ULong64_t obtt = 0;
+   Long64_t obtt = 0LL;
-   UInt_t pkth = 0;
+   Long64_t pkth = 0LL;
-   ULong64_t obth = 0;
+   Long64_t obth = 0LL;
    //
    T->GetEntry(upperentry);
    code = eh->GetCounter();
-Line 400 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 456 
 Int_t PamelaDBOperations::SetUpperLimits
      obth = OBT(phh->GetOrbitalTime());
    };
    //
-   if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
+   if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
    if ( pkth > upperpkt && obth > upperobt ){
-     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt);
+     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
      upperpkt = pkth;
      upperobt = obth;
      rhev = lasthead+1;
    } else {
      rhev = lasthead;
    };
-   if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
+   if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
    //
-   if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
+   if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
    if ( pktt > upperpkt && obtt > upperobt ){
-     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt);
+     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
      upperpkt = pktt;
      upperobt = obtt;
      rtev = lasttrail+1;
    } else {
      rtev = lasttrail;
    };
-   if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
+   if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
    //  goto kikko;
    //
    //
    // Check if runtrailer/runheader are within lower limits
    //
    //
-   pkth = 0;
+   pkth = 0LL;
-   obth = 0;
+   obth = 0LL;
-   spkth = 0;
+   spkth = 0LL;
-   sobth = 0;
+   sobth = 0LL;
    for (Int_t k=0; k<rhev; k++){
      if ( k > 0 ){
        spkth = pkth;
-Line 478 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 534 
 Int_t PamelaDBOperations::SetUpperLimits
            upperentry = evbefh-1;
          };
        };
-       if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
+       if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
        goto kikko0;
      };
    };
-Line 486 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 542 
 Int_t PamelaDBOperations::SetUpperLimits
    //
    //
    //
-   pktt = 0;
+   pktt = 0LL;
-   obtt = 0;
+   obtt = 0LL;
-   spktt = 0;
+   spktt = 0LL;
-   sobtt = 0;
+   sobtt = 0LL;
    for (Int_t k=0; k<rtev; k++){
      if ( k > 0 ){
        spktt = pktt;
-Line 509 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 565 
 Int_t PamelaDBOperations::SetUpperLimits
        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);
+       if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
        //
        UInt_t evbeft = 0;
        code = eht->GetCounter();
-Line 538 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 594 
 Int_t PamelaDBOperations::SetUpperLimits
            upperentry = evbeft-1;
          };
        };
-       if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
+       if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
        goto kikko;
        //      break;
        //
-Line 566 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 622 
 Int_t PamelaDBOperations::SetUpperLimits
      obth = OBT(phh->GetOrbitalTime());
    };
    //
-   if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
+   if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
    if ( pkth > upperpkt && obth > upperobt ){
-     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt);
+     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
      upperpkt = pkth;
      upperobt = obth;
      rhev = lasthead+1;
    } else {
      rhev = lasthead;
    };
-   if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
+   if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
    //
-   if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
+   if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
    if ( pktt > upperpkt && obtt > upperobt ){
-     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt);
+     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
      upperpkt = pktt;
      upperobt = obtt;
      rtev = lasttrail+1;
    } else {
      rtev = lasttrail;
    };
-   if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
+   if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
    //
-   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);
    //
    return(0);
  }
+ /**
+  *
+  * Trick to have unique RUN ID numbers even when runs are deleted and mysql deamon restarted.
+  * Entries in the _RUNID_GEN table are never deleted.
+  *
+  **/
+ UInt_t PamelaDBOperations::AssignRunID(){
+   //
+   TSQLResult *result = 0;
+   TSQLRow *row = 0;
+   UInt_t runid = 0;
+   //
+   stringstream   oss;
+   //
+   oss.str("");
+   oss << "INSERT INTO _RUNID_GEN VALUES (NULL);";
+   result = conn->Query(oss.str().c_str());
+   if ( !result ) throw -10;
+   oss.str("");
+   oss << "SELECT ID FROM _RUNID_GEN ORDER BY ID DESC LIMIT 1;";
+   result = conn->Query(oss.str().c_str());
+   if ( !result ) throw -10;
+   //
+   row = result->Next();
+   //
+   if ( !row ) throw -28;
+   //
+   runid = (UInt_t)atoll(row->GetField(0));
+   //
+   return(runid);
+ };
  //
  // GETTERS
  //
-Line 690 
 void PamelaDBOperations::CheckConnection
+Line 778 
 void PamelaDBOperations::CheckConnection
    if( !conn ) throw -1;
    bool connect = conn->IsConnected();
    if( !connect ) throw -1;
+   if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27;
+   //
+   // set DB timezone to UTC
+   //
+   stringstream oss;
+   //
+   oss.str("");
+   oss << "SET time_zone='+0:00';";
+   TSQLResult *result = 0;
+   result = conn->Query(oss.str().c_str());
+   if ( !result ) throw -10;
+   //
  };
  /**
   * Return the correct packet number if we went back to zero
   */
- UInt_t PamelaDBOperations::PKT(UInt_t pkt_num){
+ Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){
    //
-   //  if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u (pktfirst - (UInt_t)(16777212/2)) is %u  \n",pkt_num,pktfirst,(pktfirst - (UInt_t)(16777212/2)));
+   if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u  (UInt_t)(16777214/2)) is %u  \n",pkt_num,pktfirst,(UInt_t)(16777214/2));
    //
-   if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2)  ) return((pkt_num+16777215));
+   if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2)  ){
+     if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL);
+     return((Long64_t)pkt_num+16777215LL);
+   };
    //
    if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){
-     if ( (pkt_num-16777215) < 0 ){
+     if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL);
-       return((16777215-pkt_num));
+     return((Long64_t)pkt_num-16777215LL);
-     } else {
-       return((pkt_num-16777215));
-     };
    };
    //
-   return(pkt_num);
+   if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)pkt_num);
+   return((Long64_t)pkt_num);
    //
  };
  /**
   * Return the correct On Board Time if we went back to zero
   */
- ULong64_t PamelaDBOperations::OBT(UInt_t obt){
+ Long64_t PamelaDBOperations::OBT(UInt_t obt){
    //
-   if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((ULong64_t)(obt+numeric_limits<UInt_t>::max()));
+   if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((Long64_t)(obt+numeric_limits<UInt_t>::max()));
    //
    if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
-     if ( (obt-numeric_limits<UInt_t>::max()) < 0  ){
+     return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
-       return((ULong64_t)(numeric_limits<UInt_t>::max()-obt));
-     } else {
-       return((ULong64_t)(obt-numeric_limits<UInt_t>::max()));
-     };
    };
    //
-   return((ULong64_t)obt);
+   return((Long64_t)obt);
  };
  /**
-Line 869 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 966 
 Int_t PamelaDBOperations::insertPamelaGL
    if ( this->GetID_RAW() == 0 )  throw -11;
    //
    oss.str("");
-   oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ID_RAW<= "
+   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='"
-       << id << " AND TO_ID_RAW >= "
+       << this->GetRawFile().Data() << "';";
-       << id << ";";
    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
    result = conn->Query(oss.str().c_str());
    if ( !result ) throw -10;
    row = result->Next();
-   if ( !row ) throw -10;
    //
-   t0 = (UInt_t)TDatime(row->GetField(0)).Convert();
+   if ( !row ){
+     oss.str("");
+     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< "
+         << dworbit << " order by FROM_ORBIT desc limit 1;";
+     if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
+     result = conn->Query(oss.str().c_str());
+     if ( !result ) throw -10;
+     row = result->Next();
+     if ( !row ) throw -10;
+   };
+   //
+   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);
+   t0 = (UInt_t)tu.GetSec();
+   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));
+   //
    /*
     * Verify that the TIMESYNC have been not already processed
     */
-Line 893 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1002 
 Int_t PamelaDBOperations::insertPamelaGL
    if (result == NULL) throw -10;
    row = result->Next();
    if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){
+     if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
      toffset = (UInt_t)atoll(row->GetField(2)) - (UInt_t)(this->OBT((UInt_t)atoll(row->GetField(1)))/1000) + t0;
+     //
+     tsync = (UInt_t)atoll(row->GetField(2));
+     obt0 = (UInt_t)atoll(row->GetField(1));
+     //
      return(1);
    };
    //
-Line 1087 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1201 
 Int_t PamelaDBOperations::insertPamelaGL
    conn->Query(oss.str().c_str());
    if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str());
    //
+   if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;
    //
+   tsync = TSYNC;
+   obt0 = OBT;
+   //
    delete result;
    return(signal);
  }
-Line 1188 
 Int_t PamelaDBOperations::assignBOOT_NUM
+Line 1306 
 Int_t PamelaDBOperations::assignBOOT_NUM
    VarDumpEvent  *vde = 0;
    VarDumpRecord *vdr = 0;
    //
+   Bool_t found = false;
    trDumpEv->SetBranchAddress("VarDump", &vde);
    if ( trDumpEv->GetEntries() > 0 ){
-     Bool_t found = false;
+     found = false;
      for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){
        trDumpEv->GetEntry(i);
-       vde->Records->GetEntries();
+       //      vde->Records->GetEntries();
-       if ( vde->Records->GetEntries()>0 ){
+       if ( vde->Records->GetEntries()>5 ){
          found = true;
          goto fill;
        };
-Line 1207 
 Int_t PamelaDBOperations::assignBOOT_NUM
+Line 1326 
 Int_t PamelaDBOperations::assignBOOT_NUM
        this->SetBOOTnumber((Int_t)vdr->VAR_VALUE);
        //
      } else {
-       if ( !this->GetBOOTnumber() ) return(4);
+       if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4);
+     };
+   } else {
+     if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2);
+   };
+   //
+   UInt_t bn = 0;
+   Bool_t afound = false;
+   if ( !found && this->AutoBoot()){
+     afound = true;
+     //
+     // Search for other files with similar timesync
+     //
+     if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0);
+     UInt_t upperts = tsync-(obt0/1000)+5;
+     UInt_t lowerts = tsync-(obt0/1000)-5;
+     oss.str("");
+     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)<"
+         << upperts
+         << " AND TIMESYNC-(OBT0/1000)>"
+         << lowerts
+         << " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;";
+     result = conn->Query(oss.str().c_str());
+     if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str());
+     //
+     if ( !result ) throw -4;;
+     found = true;
+     if ( result->GetRowCount()<3 ){
+       if ( IsDebug() ) printf(" AGH! no results!\n");
+       found = false;
+     } else {
+       row = result->Next();
+       bn = (UInt_t)atoll(row->GetField(0));
+       for ( Int_t r=1; r<result->GetRowCount() ;r++){
+         if ( !row ) throw -4;
+         if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0));
+         if ( bn != (UInt_t)atoll(row->GetField(0)) ){
+           if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0)));
+           found = false;
+         };
+         row = result->Next();
+       };
      };
+   };
+   //
+   Int_t sgn = 0;
+   //
+   if ( !found ){
+     throw -29;
    } else {
-     if ( !this->GetBOOTnumber() ) return(2);
+     if ( afound ){
+       this->SetBOOTnumber(bn);
+       sgn = 8;
+     };
    };
    //
    oss.str("");
-Line 1220 
 Int_t PamelaDBOperations::assignBOOT_NUM
+Line 1389 
 Int_t PamelaDBOperations::assignBOOT_NUM
    conn->Query(oss.str().c_str());
    //
    delete result;
-   return(0);
+   return(sgn);
  };
  /**
-Line 1483 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1652 
 Bool_t PamelaDBOperations::IsRunAlreadyI
      if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){
        //
        if ( IsDebug() ) printf(" The new run has more events than the old one \n");
-       oss.str("");
+       glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
-       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
+ //       oss.str("");
-       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) <<";";
-       conn->Query(oss.str().c_str());
+ //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());
+ //       conn->Query(oss.str().c_str());
        if ( signal ) signal = false;
        goto gonext;
        //
-Line 1500 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1670 
 Bool_t PamelaDBOperations::IsRunAlreadyI
        //
        if ( IsDebug() ) printf(" The new run has the same number of events and the runheader the old one miss the runheader \n");
        //
-       oss.str("");
+       glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
-       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
+ //       oss.str("");
-       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) <<";";
-       conn->Query(oss.str().c_str());
+ //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());
+ //       conn->Query(oss.str().c_str());
        //
        if ( signal ) signal = false;
        goto gonext;
-Line 1517 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1688 
 Bool_t PamelaDBOperations::IsRunAlreadyI
        //
        if ( IsDebug() ) printf(" The new run has the same number of events, the runheader and the runtrailer the old one miss the runtrailer \n");
        //
-       oss.str("");
+       glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
-       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
+ //       oss.str("");
-       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) <<";";
-       conn->Query(oss.str().c_str());
+ //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());
+ //       conn->Query(oss.str().c_str());
        if ( signal ) signal = false;
        //
      };
-Line 1535 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1707 
 Bool_t PamelaDBOperations::IsRunAlreadyI
    delete result;
    //
    if ( signal && IsDebug() ) printf(" The run has already been inserted \n");
+   if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n");
    return(signal);
  };
-Line 1556 
 void PamelaDBOperations::HandleRun(){
+Line 1729 
 void PamelaDBOperations::HandleRun(){
      //
      this->FillClass();
      //
-     if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+     if ( !IsRunAlreadyInserted() ){
+       glrun->SetID(this->AssignRunID());
+       glrun->SetID_RUN_FRAG(0);
+       glrun->Fill_GL_RUN(conn);
+     };
    } else {
      //
      if ( IsDebug() ) printf(" oh no! the distance between runheader and runtrailer seems wrong: check %llu pktt %llu \n",chkpkt,pktt);
-Line 1579 
 void PamelaDBOperations::HandleRunFragme
+Line 1756 
 void PamelaDBOperations::HandleRunFragme
    UInt_t rhfirstev = firstev;
    UInt_t rtlastev = lastev;
    Bool_t found = false;
+   Bool_t foundinrun = false;
    //
    TSQLResult *result = 0;
    TSQLRow    *row    = 0;
-Line 1605 
 void PamelaDBOperations::HandleRunFragme
+Line 1783 
 void PamelaDBOperations::HandleRunFragme
    //
    if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n");
    //
+   // First of all insert the run in the fragment table...
+   //
+   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");
+     return;
+   };
+   //
+   //
    // can we find the other piece of the run in the GL_RUN_FRAGMENTS table?
    //
    if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is
-Line 1613 
 void PamelaDBOperations::HandleRunFragme
+Line 1833 
 void PamelaDBOperations::HandleRunFragme
      oss.str("");
      oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE "
          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
-         << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME()
+         << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
+         << " ID != " << glrun->ID
          << " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
      //
      if ( IsDebug() ) printf(" look for runheader in the fragments table: query is \n %s \n",oss.str().c_str());
-Line 1623 
 void PamelaDBOperations::HandleRunFragme
+Line 1844 
 void PamelaDBOperations::HandleRunFragme
      //
      row = result->Next();
      //
+     if ( !row && NoFrag() ){
+       //
+       oss.str("");
+       oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN WHERE "
+           << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
+           << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
+           << " ID != " << glrun->ID
+           << " AND ID=ID_RUN_FRAG ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
+       //
+       if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
+       result = conn->Query(oss.str().c_str());
+       //
+       if ( !result ) throw -4;
+       //
+       foundinrun = true;
+       //
+       row = result->Next();
+       //
+     };
+     //
      if ( !row ){
        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
        found = false;
-Line 1697 
 void PamelaDBOperations::HandleRunFragme
+Line 1938 
 void PamelaDBOperations::HandleRunFragme
        //
        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");
        //
+       if ( foundinrun ){
+         glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
+         glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
+       };
+       //
        GL_RUN *glrun1 = new GL_RUN();
        //
-       UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
+       //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
        //
        oss.str("");
        oss << " ID="<<row->GetField(0)<<";";
-Line 1728 
 void PamelaDBOperations::HandleRunFragme
+Line 1974 
 void PamelaDBOperations::HandleRunFragme
        };
        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
        //
-       glrun1->SetID(0);
        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());
-Line 1737 
 void PamelaDBOperations::HandleRunFragme
+Line 1982 
 void PamelaDBOperations::HandleRunFragme
        //
        glrun->SetEV_FROM(firstev);
        glrun->SetNEVENTS(lastev-firstev+1);
+       //
        glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME());
        glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT());
        glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT());
-Line 1756 
 void PamelaDBOperations::HandleRunFragme
+Line 2002 
 void PamelaDBOperations::HandleRunFragme
        //
        if ( !IsRunAlreadyInserted() ){
          //
+         //      glrun->SetID(this->AssignRunID());
+         glrun->SetID_RUN_FRAG(glrun1->GetID());
          glrun->Fill_GL_RUN(conn);
          //
-         // 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));
-         };
-         //
-         glrun1->SetID_RUN_FRAG(idrun0);
          //
+         glrun1->SetID_RUN_FRAG(glrun->GetID());
          glrun1->Fill_GL_RUN(conn);
          //
-         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;
-         //
        };
-       //
-       delete glrun1;
-       //
        // delete old entry in fragment table
        //
-       oss.str("");
+       glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
+       glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
        //
-       oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";";
+       delete glrun1;
        //
-       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;
        //
        return;
        //
-Line 1825 
 void PamelaDBOperations::HandleRunFragme
+Line 2027 
 void PamelaDBOperations::HandleRunFragme
    };
    //
    if ( mistrail && rtlastev == lastev ) { // look for runtrailer (only when at the end of the file, if at the beginning and the runh is
-                                            // 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
      //
      oss.str("");
      oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN_FRAGMENTS WHERE "
          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
-         << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME()
+         << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
+         << " ID != " << glrun->ID
          << " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
      //
      if ( IsDebug() ) printf(" look for runtrailer in the fragments table: query is \n %s \n",oss.str().c_str());
-Line 1840 
 void PamelaDBOperations::HandleRunFragme
+Line 2043 
 void PamelaDBOperations::HandleRunFragme
      //
      row = result->Next();
      //
+     if ( !row && NoFrag() ){
+       //
+       oss.str("");
+       oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN WHERE "
+           << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
+           << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
+           << " ID != " << glrun->ID
+           << " AND ID=ID_RUN_FRAG ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
+       //
+       if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
+       result = conn->Query(oss.str().c_str());
+       //
+       if ( !result ) throw -4;
+       //
+       foundinrun = true;
+       row = result->Next();
+       //
+     };
+     //
      if ( !row ){
        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
        found = false;
-Line 1913 
 void PamelaDBOperations::HandleRunFragme
+Line 2135 
 void PamelaDBOperations::HandleRunFragme
        //
        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");
        //
+       if ( foundinrun ){
+         glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
+         glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
+       };
+       //
        GL_RUN *glrun1 = new GL_RUN();
        //
-       UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
+       //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
        //
        oss.str("");
        oss << " ID="<<row->GetField(0)<<";";
-Line 1952 
 void PamelaDBOperations::HandleRunFragme
+Line 2179 
 void PamelaDBOperations::HandleRunFragme
        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());
        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());
        //
-       glrun1->SetID(0);
        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());
        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());
        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());
-Line 1972 
 void PamelaDBOperations::HandleRunFragme
+Line 2198 
 void PamelaDBOperations::HandleRunFragme
        //
        if ( !IsRunAlreadyInserted() ){
          //
-         glrun->Fill_GL_RUN(conn);
+         //      glrun->SetID(this->AssignRunID());
-         //
-         // get 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;
+         glrun->SetID_RUN_FRAG(glrun1->GetID());
-         if ( !row ){
+         glrun->Fill_GL_RUN(conn);
-           throw -10;
-         } else {
-           idrun0 = (UInt_t)atoll(row->GetField(0));
-         };
          //
-         glrun1->SetID_RUN_FRAG(idrun0);
+         // set id number
          //
+         glrun1->SetID_RUN_FRAG(glrun->GetID());
          glrun1->Fill_GL_RUN(conn);
          //
-         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;
-         //
        };
        //
-       delete glrun1;
+       // delete old entries in fragment table
-       //
-       // delete old entry in fragment table
-       //
-       oss.str("");
-       //
-       oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";";
        //
-       if ( IsDebug() ) printf(" Delete from frag table the old run :\n query is \n %s \n",oss.str().c_str());
+       glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
-       result =  conn->Query(oss.str().c_str());
+       glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
-       if ( !result ) throw -4;
        //
+       delete glrun1;
        //
        return;
        //
-Line 2070 
 void PamelaDBOperations::HandleRunFragme
+Line 2252 
 void PamelaDBOperations::HandleRunFragme
      //
      row = result->Next();
      //
-     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->Fill_GL_RUN_FRAGMENTS(conn);
-         //
-       } else {
-         if ( IsDebug() ) printf(" The run is already present in the fragment table \n");
-       };
-     } else {
        if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n");
+     } else {
+       if ( NoFrag() ){
+         glrun->SetID_RUN_FRAG(glrun->GetID());
+         glrun->Fill_GL_RUN(conn);
+         glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
+       };
      };
    };
    //
-Line 2134 
 void PamelaDBOperations::HandleMissingHo
+Line 2285 
 void PamelaDBOperations::HandleMissingHo
      //
      this->FillClass(mishead,mistrail,firstev,lastev);
      //
-     if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+     if ( !IsRunAlreadyInserted() ){
+       glrun->SetID(this->AssignRunID());
+       glrun->SetID_RUN_FRAG(0);
+       glrun->Fill_GL_RUN(conn);
+     };
      //
    };
    //
-Line 2301 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 2456 
 Bool_t PamelaDBOperations::IsRunConsiste
            //
            this->SetCommonGLRUN(firstTime,lastTime);
            //
-           if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+           if ( !IsRunAlreadyInserted() ){
+             glrun->SetID(this->AssignRunID());
+             glrun->SetID_RUN_FRAG(0);
+             glrun->Fill_GL_RUN(conn);
+           };
            //
            firstevno = lastentry + 1;
            //
-Line 2377 
 void PamelaDBOperations::HandleSuspiciou
+Line 2536 
 void PamelaDBOperations::HandleSuspiciou
      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
      //
      this->FillClass();
-     if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+     if ( !IsRunAlreadyInserted() ){
+       glrun->SetID(this->AssignRunID());
+       glrun->SetID_RUN_FRAG(0);
+       glrun->Fill_GL_RUN(conn);
+     };
      //
    } else {
      //
-Line 2405 
 void PamelaDBOperations::HandleSuspiciou
+Line 2568 
 void PamelaDBOperations::HandleSuspiciou
        if ( IsDebug() ) printf(" No packets but physics inside the run, I will consider it as good\n");
        //
        this->FillClass();
-       if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+       if ( !IsRunAlreadyInserted() ){
+         glrun->SetID(this->AssignRunID());
+         glrun->SetID_RUN_FRAG(0);
+         glrun->Fill_GL_RUN(conn);
+       };
        //
      } else {
        //
-Line 2503 
 void PamelaDBOperations::HandleSuspiciou
+Line 2670 
 void PamelaDBOperations::HandleSuspiciou
            //
            this->SetCommonGLRUN(firstTime,lastTime);
            //
-           if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+           if ( !IsRunAlreadyInserted() ){
+             glrun->SetID(this->AssignRunID());
+             glrun->SetID_RUN_FRAG(0);
+             glrun->Fill_GL_RUN(conn);
+           };
            //
            if ( i == lastev && checkfirst < check ){ // if the last event gives a wrong check...
              //
-Line 2538 
 void PamelaDBOperations::HandleSuspiciou
+Line 2709 
 void PamelaDBOperations::HandleSuspiciou
              //
              this->SetCommonGLRUN(firstTime,lastTime);
              //
-             if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);
+             if ( !IsRunAlreadyInserted() ){
+               glrun->SetID(this->AssignRunID());
+               glrun->SetID_RUN_FRAG(0);
+               glrun->Fill_GL_RUN(conn);
+             };
            };
            //
            firstevno = lastentry + 1;
-Line 2645 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 2820 
 Int_t PamelaDBOperations::insertCALO_CAL
                // no calibrations in the db contain our calibration
                //
                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
-               if ( fromtime < 1150863400 ){
+               if ( fromtime < 1150871000 ){ //1150866904
                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
                };
-Line 2790 
 void PamelaDBOperations::HandleTRK_CALIB
+Line 2965 
 void PamelaDBOperations::HandleTRK_CALIB
        // no calibrations in the db contain our calibration
        //
        if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n");
-       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
        //
        oss.str("");
        oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "
-Line 2919 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3094 
 Int_t PamelaDBOperations::insertTRK_CALI
      pkt1 = ph1->GetCounter();
      fromtime = this->GetAbsTime(ph1->GetOrbitalTime());
      //
-     valid = 1;
+ //     valid = 1;
-     //
+ //     //
-     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
      //
      //
      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
        //
        if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);
+       //
+       valid = ValidateTrkCalib( caltrk1, eh1 );
+       if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
        //
        // Do we have the second calibration packet?
        //
-Line 2943 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3121 
 Int_t PamelaDBOperations::insertTRK_CALI
            obt2 = ph2->GetOrbitalTime();
            pkt2 = ph2->GetCounter();
            //
-           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
            //
          } else {
            //
-Line 2974 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3152 
 Int_t PamelaDBOperations::insertTRK_CALI
          //
          if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2);
          //
+         UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
+         if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
+         valid = valid & valid2;
+         //
          // Handle good calib
          //
          this->HandleTRK_CALIB(true,true);
-Line 3146 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 3328 
 Int_t PamelaDBOperations::insertS4_CALIB
            // no calibrations in the db contain our calibration
            //
            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");
-           if ( fromtime < 1150863400 ){
+           if ( fromtime < 1150871000 ){
              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
            };
-Line 3219 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 3401 
 Int_t PamelaDBOperations::insertS4_CALIB
   * Scan the fragment table and move old fragments to the GL_RUN table
   */
  Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){
+   return(this->CleanGL_RUN_FRAGMENTS(""));
+ };
+ /**
+  * Scan the fragment table and move old fragments to the GL_RUN table
+  */
+ Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(TString fcleanfile){
    //
    TSQLResult *result = 0;
    TSQLRow    *row    = 0;
-Line 3226 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
+Line 3415 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
    TSQLRow    *row2   = 0;
    //
    UInt_t moved = 0;
- //  UInt_t timelim = 0;
- //  TDatime *time = new TDatime();
    //
    stringstream oss;
    oss.str("");
    //
-   //
+   if ( !strcmp(fcleanfile.Data(),"") ){
-   // check if there are entries older than "olderthan" seconds from now
+     //
-   //
+     // check if there are entries older than "olderthan" seconds from now
-   oss.str("");
+     //
-   oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
+     oss.str("");
-       << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";
+     oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
-   //
+         << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";
-   if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str());
+     //
-   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());
+     result = conn->Query(oss.str().c_str());
+     //
+   } else {
+     oss.str("");
+     oss << " SELECT ID FROM GL_ROOT WHERE NAME='" << fcleanfile.Data() << "';";
+     if ( IsDebug() ) printf(" Getting ID_ROOT_L0 query %s \n",oss.str().c_str());
+     result = conn->Query(oss.str().c_str());
+     //
+     if ( result ){
+       //
+       row = result->Next();
+       //
+       if ( row ){
+         oss.str("");
+         oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
+             << " ID_ROOT_L0=" << row->GetField(0) << ";";
+         //
+         if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS for ROOT file query is \n %s \n",oss.str().c_str());
+         result = conn->Query(oss.str().c_str());
+         //
+       };
+     } else {
+       return(2);
+     };
+   };
    //
    if ( result ){
      //
-Line 3281 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
+Line 3493 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
          if ( IsDebug() ) printf(" The run is new \n");
          if ( IsDebug() ) printf(" -> fill the DB \n");
          //
-         glrun->SetID(0);
+         //      glrun->SetID(this->AssignRunID()); we use the old run number!
+         glrun->SetID_RUN_FRAG(glrun->GetID());
          glrun->Fill_GL_RUN(conn);
          //
-         oss.str("");
+ //      oss.str("");
-         oss << " SELECT ID FROM GL_RUN WHERE "
+ //      oss << " SELECT ID FROM GL_RUN WHERE "
-             << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "
+ //          << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "
-             << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "
+ //          << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "
-             << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "
+ //          << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "
-             << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "
+ //          << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "
-             << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";
+ //          << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";
-         //
+ //      //
-         if ( IsDebug() ) printf(" Look for the ID of the inserted run: query is \n %s \n",oss.str().c_str());
+ //      if ( IsDebug() ) printf(" Look for the ID of the inserted run: query is \n %s \n",oss.str().c_str());
-         result2 = conn->Query(oss.str().c_str());
+ //      result2 = conn->Query(oss.str().c_str());
-         //
+ //      //
-         if ( !result2 ) throw -4;
+ //      if ( !result2 ) throw -4;
-         //
+ //      //
-         row2 = result2->Next();
+ //      row2 = result2->Next();
-         //
+ //      //
-         if ( !row2 ) throw -25;
+ //      if ( !row2 ) throw -25;
-         //
+ //      //
-         oss.str("");
+ //      oss.str("");
-         oss << " UPDATE GL_RUN SET ID_RUN_FRAG = " << row2->GetField(0) << " WHERE ID = " << row2->GetField(0);
+ //      oss << " UPDATE GL_RUN SET ID_RUN_FRAG = " << row2->GetField(0) << " WHERE ID = " << row2->GetField(0);
-         if ( IsDebug() ) printf(" Update the ID_RUN_FRAG of the inserted run: query is \n %s \n",oss.str().c_str());
+ //      if ( IsDebug() ) printf(" Update the ID_RUN_FRAG of the inserted run: query is \n %s \n",oss.str().c_str());
-         result2 = conn->Query(oss.str().c_str());
+ //      result2 = conn->Query(oss.str().c_str());
-         //
+ //      //
-         if ( !result2 ) throw -4;
+ //      if ( !result2 ) throw -4;
          //
          moved++;
          //
-Line 3316 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
+Line 3529 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
        if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0));
        //
        //
-       oss.str("");
+       glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
-       oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0);
+ //      oss.str("");
-       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);
-       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());
-       //
+ //       result2 = conn->Query(oss.str().c_str());
-       if ( !result2 ) throw -4;
+ //      //
-       //
+ //      if ( !result2 ) throw -4;
+ //      //
        row = result->Next();
      };
    };
-Line 3334 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
+Line 3548 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
   * Check if runs are good, i.e. if the tracker calibration is correctly associated..
   */
  Int_t PamelaDBOperations::ValidateRuns(){
+   return(this->ValidateRuns(""));
+ };
+ /**
+  * Check if runs are good, i.e. if the tracker calibration is correctly associated..
+  */
+ Int_t PamelaDBOperations::ValidateRuns(TString valfile){
    //
    TSQLResult *result = 0;
    TSQLRow    *row    = 0;
    //
+   UInt_t calibtime = 50;
+   //
    stringstream oss;
    oss.str("");
    //
-Line 3346 
 Int_t PamelaDBOperations::ValidateRuns()
+Line 3569 
 Int_t PamelaDBOperations::ValidateRuns()
    // =======================================================
    UInt_t t_stop  = 0;
    UInt_t t_start = 0;
-   // --------------------------------------------------------------
+   if ( !strcmp(valfile.Data(),"") ) {
-   // 1) get the OBT of the last run inserted after clean-time limit
+     // --------------------------------------------------------------
-   // --------------------------------------------------------------
+     // 1) get the OBT of the last run inserted after clean-time limit
-   oss.str("");
+     // --------------------------------------------------------------
-   oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()
+     oss.str("");
-           << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
+     oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()
-   if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());
+         << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
-   result = conn->Query(oss.str().c_str());
+     if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());
-   if ( !result ) throw -4;
+     result = conn->Query(oss.str().c_str());
-   if ( !result->GetRowCount() ) {
+     if ( !result ) throw -4;
-           printf(" No runs to validate \n");
+     if ( !result->GetRowCount() ) {
-           return(1);
+       printf(" No runs to validate \n");
-   }else{
+       return(1);
-         row = result->Next();
+     }else{
-         t_start = (UInt_t)atoll(row->GetField(4));
+       row = result->Next();
-   };
+       t_start = (UInt_t)atoll(row->GetField(4));
-   // --------------------------------------------------------------
+     };
-   // 2) get the OBT of the last validated run
+     // --------------------------------------------------------------
-   // --------------------------------------------------------------
+     // 2) get the OBT of the last validated run
-   oss.str("");
+     // --------------------------------------------------------------
-   oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
+     oss.str("");
-       <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
+     oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
-   if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
+         <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
-   result = conn->Query(oss.str().c_str());
+     if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
-   if ( !result ) throw -4;
+     result = conn->Query(oss.str().c_str());
-   if ( result->GetRowCount() ){
+     if ( !result ) throw -4;
-           row = result->Next();
+     if ( result->GetRowCount() ){
-           t_stop = (UInt_t)atoll(row->GetField(4));
+       row = result->Next();
+       t_stop = (UInt_t)atoll(row->GetField(4));
+     };
+     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() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());
+     result = conn->Query(oss.str().c_str());
+   } else {
+     //
+     stringstream myquery;
+     UInt_t myid = 0;
+     myquery.str("");
+     myquery << " SELECT ID FROM GL_ROOT where NAME='"<<valfile.Data() <<"';";
+     //
+     result = conn->Query(myquery.str().c_str());
+     //
+     row = result->Next();
+     if( !row ){
+       if ( strcmp(valfile.Data(),GetRootName().Data()) ){
+         if ( IsDebug() ) printf(" No file to be validated even if option \"-validate file\" was used!!\n");
+         return(2);
+       };
+       if ( IsDebug() ) printf(" No file to be validated (force mode)! \n");
+       return(0);
+     };
+     myid=(UInt_t)atoll(row->GetField(0));
+     //
+     myquery.str("");
+     myquery << " SELECT MAX(RUNTRAILER_TIME),MIN(RUNHEADER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< myid <<";";
+     //
+     result = conn->Query(myquery.str().c_str());
+     //
+     row = result->Next();
+     if( !row->GetField(0) || !row->GetField(1)){
+       //
+       if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
+       //
+       return(0);
+       //
+     } else {
+       //
+       UInt_t runhtime = (UInt_t)atoll(row->GetField(0));
+       UInt_t runttime = (UInt_t)atoll(row->GetField(1));
+       UInt_t caltime = 0;
+       //
+       myquery.str("");
+       myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
+       myquery << " order by FROM_TIME asc limit 1;";
+       //
+       if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+       //
+       //
+       result = conn->Query(myquery.str().c_str());
+       //
+       row = result->Next();
+       if( !row ){
+         caltime = runhtime;
+       } else {
+         caltime = (UInt_t)atoll(row->GetField(0));
+       };
+       //
+       myquery.str("");
+       myquery << " SELECT * from GL_RUN where RUNHEADER_TIME>="<< runttime <<" AND RUNHEADER_TIME<=" ;
+       myquery << caltime << " order by RUNHEADER_TIME DESC";
+       //
+       if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+       //
+       result = conn->Query(myquery.str().c_str());
+       //
+     };
    };
-   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());
    if ( !result ) throw -4;
-   if ( !result->GetRowCount() ) printf(" No runs to validate \n");
+   if ( !result->GetRowCount() && IsDebug() ) printf(" No runs to validate \n");
- //  printf("------------------------------------------------------------------------------- \n");
+   //
    Int_t nrow = 0;
    GL_RUN* this_run = new GL_RUN();
    GL_RUN* next_run = new GL_RUN();
-Line 3403 
 Int_t PamelaDBOperations::ValidateRuns()
+Line 3690 
 Int_t PamelaDBOperations::ValidateRuns()
    UInt_t t1=0,t2=0;
    // ---------------------------------------------------------------------------------
    // - loop over runs, back in time,
-   // - select sequences of runs close in time (less than 60 s apart),
+   // - select sequences of runs close in time (less than calibtime s apart),
    //   which could be preceeded by a calibration
    // - check if there might be a missing calibration
    // ---------------------------------------------------------------------------------
-Line 3443 
 Int_t PamelaDBOperations::ValidateRuns()
+Line 3730 
 Int_t PamelaDBOperations::ValidateRuns()
                            if( this_run->ID == next_run->ID_RUN_FRAG ) interval = 0;     //=> run fragments
-                           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
                            if( !CHECK && this_run->VALIDATION ){
                                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);
-Line 3641 
 Int_t PamelaDBOperations::populateTLE()/
+Line 3928 
 Int_t PamelaDBOperations::populateTLE()/
  // Insert tle in the table GL_TLE using the connection conn.
- int PamelaDBOperations::insertTle(cTle *tle)
+ Int_t PamelaDBOperations::insertTle(cTle *tle)
  {
    stringstream oss;
    TSQLResult *result = 0;
-Line 3724 
 string getTleDatetime(cTle *tle)
+Line 4011 
 string getTleDatetime(cTle *tle)
    return date.str();
  }
+ /**
+  * Remove a file from the DB, delete on cascade all entries related to that file
+  * rearrange GL_RUN and GL_XXX_CALIB tables, turn off validation till the following
+  * calibration
+  **/
+ Int_t PamelaDBOperations::removeFile(TString remfile){
+   //
+   // Determine ID_ROOT_L0 and ID_RAW
+   //
+   TSQLResult *pResult;
+   TSQLRow *Row;
+   stringstream myquery;
+   //
+   myquery.str("");
+   myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   if( !Row ){
+     if ( strcmp(remfile.Data(),GetRootName().Data()) ){
+       if ( IsDebug() ) printf(" No file to be removed even if option \"-remove file\" was used!!\n");
+       return(1);
+     };
+     if ( IsDebug() ) printf(" No file to be removed (force mode)! \n");
+     return(0);
+   };
+   //
+   this->SetID_ROOT((UInt_t)atoll(Row->GetField(0)));
+   this->SetID_RAW((UInt_t)atoll(Row->GetField(1)));
+   //
+   this->ValidationOFF();
+   //
+   this->RemoveCALIBS();
+   //
+   this->RemoveRUNS();
+   //
+   this->RemoveFILES();
+   //
+   this->SetID_ROOT(0);
+   this->SetID_RAW(0);
+   //
+   return(0);
+ };
+ /**
+  *
+  * Set validation bit to zero for runs following the removing file till
+  * 1) a run with TRK_CALIB_USED=140
+  * 2) a run with VALIDATION = 0
+  * 3) the next calibration
+  *
+  **/
+ void PamelaDBOperations::ValidationOFF(){
+   TSQLResult *pResult;
+   TSQLRow *Row;
+   stringstream myquery;
+   Int_t unv = 0;
+   //select ID from GL_RUN where RUNHEADER_TIME>=1152671382 AND (VALIDATION=0 OR TRK_CALIB_USED=104) order by RUNHEADER_TIME asc limit 1;
+   myquery.str("");
+   myquery << " SELECT MAX(RUNTRAILER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   if( !Row->GetField(0) ){
+     //
+     if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
+     //
+   } else {
+     //
+     UInt_t runhtime = (UInt_t)atoll(Row->GetField(0));
+     UInt_t caltime = 0;
+     //
+     myquery.str("");
+     myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
+     myquery << " order by FROM_TIME asc limit 1;";
+     //
+     if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+     //
+     //
+     delete pResult;
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     Row = pResult->Next();
+     if( !Row ){
+       caltime = runhtime;
+     } else {
+       caltime = (UInt_t)atoll(Row->GetField(0));
+     };
+     //
+     myquery.str("");
+     myquery << " SELECT ID,RUNHEADER_TIME from GL_RUN where RUNHEADER_TIME>="<< runhtime <<" AND (VALIDATION=0 OR TRK_CALIB_USED=104 OR RUNHEADER_TIME>" ;
+     myquery << caltime << ") order by RUNHEADER_TIME asc LIMIT 1";
+     //
+     if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     Row = pResult->Next();
+     if( !Row ){
+       //
+       if ( IsDebug() ) printf(" NO RUN NEED TO BE UNVALIDATED \n");
+       //
+     } else {
+       myquery.str("");
+       myquery << " SELECT ID from GL_RUN where RUNHEADER_TIME<"<< Row->GetField(1) <<" AND ";
+       myquery << " RUNHEADER_TIME>=" <<runhtime;
+       myquery << " order by RUNHEADER_TIME asc;";
+       //
+       if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+       //
+       pResult = conn->Query(myquery.str().c_str());
+       //
+       Row = pResult->Next();
+       while ( Row ){
+         //
+         unv++;
+         this->assignVALIDATION((UInt_t)atoll(Row->GetField(0)), false);
+         Row = pResult->Next();
+         //
+       };
+     };
+   };
+   if ( IsDebug() ) printf(" %i runs have been unvalidated \n",unv);
+ };
+ /**
+  *
+  * Rearrange GL_RUN table and remove runs
+  *
+  **/
+ void PamelaDBOperations::RemoveRUNS(){
+   TSQLResult *pResult;
+   TSQLRow *Row;
+   stringstream myquery;
+   UInt_t drun = 0;
+   GL_RUN *delrun = new GL_RUN();
+   //
+   myquery.str("");
+   myquery << " SELECT ID FROM GL_RUN where ID_RUN_FRAG=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
+   //
+   if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   //
+   //
+   if ( !Row ){
+     if ( IsDebug() ) printf(" No run with ID_RUN_FRAG=0 belonged to this file \n");
+   } else {
+     if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");
+     while ( Row ){
+       delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
+       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+       drun++;
+       Row = pResult->Next();
+     };
+   };
+   //
+   //
+   myquery.str("");
+   myquery << " SELECT ID,ID_RUN_FRAG FROM GL_RUN where ID_RUN_FRAG!=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
+   //
+   if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   //
+   if ( !Row ){
+     if ( IsDebug() ) printf(" No run with ID_RUN_FRAG!=0 belonged to this file \n");
+   } else {
+     if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");
+     while ( Row ){
+       if ( IsDebug() ) printf(" restore run %i \n",(UInt_t)atoll(Row->GetField(1)));
+       delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");
+       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(1)));
+       delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");
+       if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){
+         if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+         delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
+       };
+       drun++;
+       Row = pResult->Next();
+     };
+   };
+   //
+   if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN table \n",drun);
+   //
+   //
+   //
+   drun = 0;
+   //
+   myquery.str("");
+   myquery << " SELECT ID_TRASH FROM GL_RUN_TRASH where BELONGED_TO='GL_RUN_FRAGMENTS' AND ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   //
+   if ( !Row ){
+     if ( IsDebug() ) printf(" No run from GL_RUN_FRAGMENTS table in the trash table for this file \n");
+   } else {
+     if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");
+     while ( Row ){
+       if ( IsDebug() ) printf(" del run idtrash %i \n",(UInt_t)atoll(Row->GetField(0)));
+       myquery.str("");
+       myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";
+       conn->Query(myquery.str().c_str());
+       drun++;
+       Row = pResult->Next();
+     };
+   };
+   //
+   if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_TRASH table \n",drun);
+   //
+   //
+   //
+   drun = 0;
+   //
+   myquery.str("");
+   myquery << " SELECT ID FROM GL_RUN_FRAGMENTS where ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   //
+   if ( !Row ){
+     if ( IsDebug() ) printf(" No run in the GL_RUN_FRAGMENTS table for this file \n");
+   } else {
+     if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");
+     while ( Row ){
+       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+       myquery.str("");
+       myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";
+       conn->Query(myquery.str().c_str());
+       drun++;
+       Row = pResult->Next();
+     };
+   };
+   //
+   if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_FRAGMENTS table \n",drun);
+   //
+   //
+   //
+   delete delrun;
+   //
+ };
+ /**
+  *
+  * Rearrange calibration tables
+  *
+  **/
+ void PamelaDBOperations::RemoveFILES(){
+   stringstream myquery;
+   //
+   myquery.str("");
+   myquery << " DELETE FROM GL_RAW WHERE ID=" <<this->GetID_RAW() <<";";
+   //
+   if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+   //
+   conn->Query(myquery.str().c_str());
+   //
+ };
+ /**
+  *
+  * Rearrange calibration tables
+  *
+  **/
+ void PamelaDBOperations::RemoveCALIBS(){
+   TSQLResult *pResult;
+   TSQLRow *Row;
+   stringstream myquery;
+   //
+   //
+   // Calorimeter
+   //
+   for (Int_t section = 0; section < 4; section++){
+     myquery.str("");
+     myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_CALO_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<" AND ";
+     myquery << " SECTION=" << section << ";";
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     Row = pResult->Next();
+     if( !Row->GetField(0) || !Row->GetField(1) ){
+       //
+       if ( IsDebug() ) printf(" NO CALO CALIBRATION SECTION %i ASSOCIATED TO THIS FILE! \n",section);
+       //
+     } else {
+       //
+       myquery.str("");
+       myquery << " UPDATE GL_CALO_CALIB SET TO_TIME=" << Row->GetField(1);
+       myquery << " WHERE TO_TIME="<< Row->GetField(0) << " AND ";
+       myquery << " SECTION=" << section << ";";
+       //
+       pResult = conn->Query(myquery.str().c_str());
+       //
+       if( !pResult ){
+         //
+         if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
+         //
+         throw -4;
+         //
+       };
+       //
+     };
+   };
+   myquery.str("");
+   myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   if( !pResult ){
+     //
+     if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
+     //
+     throw -4;
+     //
+   };
+   //
+   // Tracker
+   //
+   myquery.str("");
+   myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_TRK_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   if( !Row->GetField(0) || !Row->GetField(1) ){
+     //
+     if ( IsDebug() ) printf(" NO TRK CALIBRATION ASSOCIATED TO THIS FILE! \n");
+     //
+   } else {
+     //
+     myquery.str("");
+     myquery << " UPDATE GL_TRK_CALIB SET TO_TIME=" << Row->GetField(1);
+     myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     if( !pResult ){
+       //
+       if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
+       //
+       throw -4;
+       //
+     };
+     //
+     myquery.str("");
+     myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     if( !pResult ){
+       //
+       if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
+       //
+       throw -4;
+       //
+     };
+   };
+   //
+   //
+   // S4
+   //
+   myquery.str("");
+   myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_S4_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   //
+   Row = pResult->Next();
+   if( !Row->GetField(0) || !Row->GetField(1) ){
+     //
+     if ( IsDebug() ) printf(" NO S4 CALIBRATION ASSOCIATED TO THIS FILE! \n");
+     //
+   } else {
+     //
+     myquery.str("");
+     myquery << " UPDATE GL_S4_CALIB SET TO_TIME=" << Row->GetField(1);
+     myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     if( !pResult ){
+       //
+       if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
+       //
+       throw -4;
+       //
+     };
+     //
+     myquery.str("");
+     myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     if( !pResult ){
+       //
+       if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
+       //
+       throw -4;
+       //
+     };
+     //
+   };
+ };
+ /**
+  *
+  * Rearrange calibration tables
+  *
+  **/
+ UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){
+     Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0};
+     UInt_t timeaftercalib=120000; //2000;
+ //  ----------
+ //  Check CRCs
+ //  ----------
+     for(Int_t ipkt=0; ipkt<6; ipkt++){
+         if( caltrk->crc_hcal[ipkt] )return 0; // :-(
+         for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-(
+     }
+ //  -----------------------
+ //  Check missing packets:
+ //  -----------------------
+ //    Readout order:
+ //    ------------------
+ //    DSP   packet board
+ //    ------------------
+ //    12    0      1
+ //    10    1      1
+ //     8    2      1
+ //     4    3      1
+ //     6    4      1
+ //     2    5      1
+ //    ------------------
+ //    11    0      2
+ //     9    1      2
+ //     7    2      2
+ //     3    3      2
+ //     5    4      2
+ //     1    5      2
+ //    ------------------
+ //  -------------------------------------------------
+ //  Check if it is first or second calibration packet
+ //  -------------------------------------------------
+     UInt_t build=0;
+     TString classname = caltrk->GetName();
+     UInt_t base=0;
+     UInt_t mask=0;
+     if(classname.Contains("CalibTrk1Event")){
+         base=12;
+         mask=0x03F000;
+     }
+     if(classname.Contains("CalibTrk2Event")){
+         base=18;
+         mask=0xFC0000;
+     }
+ //  -------------------------------------------------
+ //  Count number of packets and set build variable
+ //  -------------------------------------------------
+     Int_t  npkts=0;
+     for(Int_t ipkt=0; ipkt<6; ipkt++){
+         if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){
+             npkts++;
+             build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) );
+         }
+     }
+ //    if( npkts==6 )return 1; // :-)
+ //    cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl;
+ //  -----------------------------------------------
+ //  If missing packets: check the acq configuration
+ //  (some DSPs might be excluded from acquisition)
+ //  -----------------------------------------------
+ //  -----------------------------------------------
+ //  retrieve the first run header after calib
+ //  -----------------------------------------------
+      PacketType *pctp;
+      EventCounter *cod;
+      cod = eh->GetCounter();
+      Int_t irun = cod->Get(pctp->RunHeader);
+      TTree *rh=(TTree*)file->Get("RunHeader");
+      if ( !rh || rh->IsZombie() ) throw -17;
+      if( rh->GetEntries() == irun ){
+           if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1)  -- cannot validate :-( "<<endl;
+          return 0; // :-(
+      }
+      RunHeaderEvent *run  = 0;
+      EventHeader    *hrun = 0;
+      rh->SetBranchAddress("RunHeader", &run);
+      rh->SetBranchAddress("Header", &hrun);
+      rh->GetEntry(irun);
+ //     cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl;
+      if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){
+           if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl;
+          return 0; // :-(
+      }
+      if( !run->RM_ACQ_AFTER_CALIB ){
+           if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0    -- cannot validate :-( "<<endl;
+          return 0; // :-(
+      }
+      UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime());
+      if( dtime > timeaftercalib ){
+          if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl;
+          return 0; // :-(
+      }
+      if( (run->ACQ_BUILD_INFO & mask) != build ){
+          if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl;
+          return 0; // :-(
+      }
+      return 1; // :-)
+ }
+ /**
+  *
+  * Check the DB (only for overlapping runs at the moment)
+  *
+  **/
+ UInt_t PamelaDBOperations::Check(){
+   //
+   UInt_t test = 0;
+   //
+   UInt_t thisrht = 0;
+   UInt_t thisrtt = 0;
+   UInt_t thisid = 0;
+   UInt_t prevrht = 0;
+   UInt_t prevrtt = 0;
+   UInt_t previd = 0;
+   //
+   UInt_t prevl0id = 0;
+   UInt_t thisl0id = 0;
+   //
+   stringstream oss;
+   TSQLResult *result = 0;
+   TSQLRow    *row    = 0;
+   TSQLResult *result2 = 0;
+   TSQLRow    *row2    = 0;
+   oss.str("");
+   oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN order by RUNHEADER_TIME asc;";
+   //  oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;";
+   result = conn->Query(oss.str().c_str());
+   //
+   if ( !result ) throw -4;;
+   //
+   row = result->Next();
+   //
+   while ( row ){
+     thisid = (UInt_t)atoll(row->GetField(0));
+     thisl0id = (UInt_t)atoll(row->GetField(1));
+     thisrht = (UInt_t)atoll(row->GetField(2));
+     thisrtt = (UInt_t)atoll(row->GetField(3));
+     //
+     //    if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){
+     //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){
+     if ( (thisrht < prevrtt) && (thisrht != prevrht) ){
+       if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
+       printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid);
+       TString prevf = "";
+       TString thisf = "";
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       prevf = (TString)row2->GetField(0);
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       thisf = (TString)row2->GetField(0);
+       if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
+       test = 1;
+     };
+     //
+     if ( (thisrtt < prevrht) && (thisrht != prevrht) ){
+       if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
+       printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid);
+       TString prevf = "";
+       TString thisf = "";
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       prevf = (TString)row2->GetField(0);
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       thisf = (TString)row2->GetField(0);
+       if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
+       test = 1;
+     };
+     //
+     if ( (thisrht > thisrtt) && (thisrht != prevrht) ){
+       if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
+       printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid);
+       TString prevf = "";
+       TString thisf = "";
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       prevf = (TString)row2->GetField(0);
+       oss.str("");
+       oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
+       result2 = conn->Query(oss.str().c_str());
+       if ( !result2 ) throw -4;;
+       row2 = result2->Next();
+       thisf = (TString)row2->GetField(0);
+       if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
+       test = 1;
+     };
+     //
+     prevrht = thisrht;
+     prevrtt = thisrtt;
+     previd = thisid;
+     prevl0id = thisl0id;
+     row = result->Next();
+   };
+   //
+   return(test);
+   //
+ };

 Legend:



Removed from v.1.12
 


changed lines


 
Added in v.1.24
 Legend:



Removed from v.1.12
 


changed lines


 
Added in v.1.24
-Removed from v.1.12
+Added in v.1.24

	ViewVC Help
Powered by ViewVC 1.1.23