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

Diff of /chewbacca/YodaProfiler/src/PamelaDBOperations.cpp

Parent Directory | Revision Log | View Patch Patch

-revision 1.3 by mocchiut,
Thu Sep 25 15:30:10 2008 UTC
+revision 1.7 by mocchiut,
Tue Nov  4 13:19:36 2008 UTC
 Line 58 
 string getTleDatetime(cTle*);
   * @param debug        debug flag.
   * @param tlefilename  ascii file with TLE 3 line elements.
   */
- 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){
+ 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, Bool_t keepenv){
    //
    chewbacca = false;
    chminentry = 0;
 Line 76 
 PamelaDBOperations::PamelaDBOperations(T
      //
    };
    //
+   KEEPENV = false;
+   if ( keepenv ) KEEPENV=true;
    STATIC=false;
    if ( staticp ) STATIC=true;
    //
-Line 138 
 void PamelaDBOperations::OpenL0File(TStr
+Line 140 
 void PamelaDBOperations::OpenL0File(TStr
        TString chfile;
        UInt_t ridn = 0;
        qu.str("");
-       qu << "SELECT ROOT_ID_N,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL from ROOT_TABLE_MERGING where ID_N=" << chID << ";";
+       qu << "SELECT ROOT_ID_N,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL,INSERTED_BY from ROOT_TABLE_MERGING where ID_N=" << chID << ";";
        if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
        result = conn->Query(qu.str().c_str());
        if ( result ){
-Line 149 
 void PamelaDBOperations::OpenL0File(TStr
+Line 151 
 void PamelaDBOperations::OpenL0File(TStr
            chpktmax = (UInt_t)atoll(row->GetField(2));
            chobtmin = (UInt_t)atoll(row->GetField(3));
            chobtmax = (UInt_t)atoll(row->GetField(4));
+           chiby = (TString)(row->GetField(5));
          } else {
            throw -84;
          };
-Line 163 
 void PamelaDBOperations::OpenL0File(TStr
+Line 166 
 void PamelaDBOperations::OpenL0File(TStr
        if ( result2 ){
          row2 = result2->Next();
          if ( row2 ){
-           chpath = (TString)gSystem->ExpandPathName(row2->GetField(0))+'/';
+           if ( KEEPENV ){
+             chpath = (TString)(row2->GetField(0))+'/';
+           } else {
+             chpath = (TString)gSystem->ExpandPathName(row2->GetField(0))+'/';
+           };
            chfile = (TString)(row2->GetField(1));
            chobtts = (UInt_t)atoll(row2->GetField(2));
            chlastts = (UInt_t)atoll(row2->GetField(3));
-Line 435 
 TString PamelaDBOperations::GetRawPath()
+Line 442 
 TString PamelaDBOperations::GetRawPath()
  TString PamelaDBOperations::GetRootPath(){
    if ( STATIC ){
-     return((TString)gSystem->DirName(filerootname.Data())+'/');
+     return((TString)gSystem->DirName(this->GetRootName().Data())+'/');
    } else {
-     return((TString)gSystem->ExpandPathName("$PAM_L0")+'/');
+     if ( KEEPENV ){
+       return((TString)gSystem->ExpandPathName(gSystem->DirName(filerootname.Data()))+'/');
+     } else {
+       return((TString)gSystem->ExpandPathName("$PAM_L0")+'/');
+     };
    };
  };
-Line 1125 
 UInt_t PamelaDBOperations::AssignRunID()
+Line 1136 
 UInt_t PamelaDBOperations::AssignRunID()
    stringstream   oss;
    //
    oss.str("");
-   oss << "INSERT INTO _RUNID_GEN VALUES (NULL);";
+   if ( chewbacca ){// if chewbacca and tag=none then use chewbacca tag (chiby = chewbacca inserted by), if chewbacca and tag!=none then use tag, if not chewbacca use tag.
+     if ( !strcmp(tag.Data(),"NONE") ){
+       oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< chiby.Data() <<"');";
+     } else {
+       oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< tag.Data() <<"');";
+     };
+   } else {
+     oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< tag.Data() <<"');";
+   };
    result = conn->Query(oss.str().c_str());
    if ( !result ) throw -10;
    oss.str("");
-Line 1273 
 void PamelaDBOperations::LockTables(){
+Line 1292 
 void PamelaDBOperations::LockTables(){
    stringstream oss;
    //
    oss.str("");
-   oss << "lock table GL_RUN write, GL_ROOT write, GL_RAW write, GL_TIMESYNC write, GL_RESURS_OFFSET write, GL_PARAM write, GL_TLE write, GL_RUN_FRAGMENTS write, GL_RUN_TRASH write, GL_CALO_CALIB write, GL_CALOPULSE_CALIB write, GL_TRK_CALIB write, GL_S4_CALIB write, ROOT_TABLE_MERGING write, ROOT_TABLE write, _RUNID_GEN write;";
+   oss << "lock table GL_RUN write, GL_ROOT write, GL_RAW write, GL_TIMESYNC write, GL_RESURS_OFFSET write, GL_PARAM write, GL_TLE write, GL_RUN_FRAGMENTS write, GL_RUN_TRASH write, GL_CALO_CALIB write, GL_CALOPULSE_CALIB write, GL_TRK_CALIB write, GL_S4_CALIB write, ROOT_TABLE_MERGING write, ROOT_TABLE_BAD write, ROOT_TABLE write, _RUNID_GEN write;";
    TSQLResult *result = 0;
    result = conn->Query(oss.str().c_str());
    if ( !result ) throw -10;
-Line 1468 
 void PamelaDBOperations::FillClass(Bool_
+Line 1487 
 void PamelaDBOperations::FillClass(Bool_
   * Insert a new row into GL_RAW table.
   */
  Int_t PamelaDBOperations::insertPamelaRawFile(){
+     //
+     Bool_t idr =  this->SetID_RAW();
+     if ( idr ) return(1);
+     GL_RAW glraw = GL_RAW();
+     glraw.PATH = GetRawPath();
+     glraw.NAME = GetRawFile();
+     //    glraw.BOOT_NUMBER = 0;//???
+     glraw.BOOT_NUMBER = this->GetBOOTnumber();
+     if( insertPamelaRawFile(&glraw) )return(1);
+     //
+     idr =  this->SetID_RAW();
+     if ( !idr ) throw -11;
+     return(0);
+ }
+ /**
+  * Insert a new row into GL_RAW table.
+  */
+ Int_t PamelaDBOperations::insertPamelaRawFile(GL_RAW *glraw){
    //
-   stringstream oss;
+     if(!glraw)return(1);//?? ok I think
-   //
+     //
-   Bool_t idr =  this->SetID_RAW();
+     stringstream oss;
-   if ( idr ) return(1);
+     //
-   //
+     oss.str("");
-   oss.str("");
+     if ( STATIC ){
-   if ( STATIC ){
+         oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
-     oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
+             << glraw->PATH << "', '" << glraw->NAME << "')";
-         << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";
+     } else {
-   } else {
+         oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << glraw->NAME << "')";
-     oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')";
+     };
-   };
+     if ( debug ) cout <<oss.str().c_str() <<endl;
-   if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
+     if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
-   //
+     //
-   idr =  this->SetID_RAW();
+     oss.str("");
-   if ( !idr ) throw -11;
+     oss << "SELECT ID FROM GL_RAW WHERE NAME=\""<<glraw->NAME<<"\";";
-   //
+     if ( debug ) cout << oss.str().c_str()<<endl;
-   return(0);
+     if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
+     //
+     TSQLResult *result = 0;
+     TSQLRow    *row    = 0;
+     result = conn->Query(oss.str().c_str());
+     if ( result == NULL ) throw -4;
+     row = result->Next();
+     if ( !row ) return(1);
+     glraw->ID = (UInt_t)atoll(row->GetField(0));
+     if ( debug ) printf(" The ID of the RAW file is %u \n",glraw->ID);
+     delete result;
+     delete row;
+     //
+     return(0);
  }
-Line 1798 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1853 
 Int_t PamelaDBOperations::insertPamelaGL
   * The raw file indicates in the parameters should be already been stored in the database.
   */
  Int_t PamelaDBOperations::insertPamelaRootFile(){
-   stringstream oss;
-   TSQLResult *result = 0;
+     stringstream oss;
-   TSQLRow    *row    = 0;
+     TSQLResult *result = 0;
-   UInt_t idtimesync = 0;
+     TSQLRow    *row    = 0;
-   //
-   //
-   if ( chewbacca ){
-     oss.str("");
-     oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;";
-     if ( debug ) printf(" query is %s \n",oss.str().c_str());
-     result = conn->Query(oss.str().c_str());
-     //
-     if ( !result ) throw -3;
      //
-     row = result->Next();
+     // ----------------------
+     // determine the timesync
+     // ----------------------
+     UInt_t idtimesync = 0;
      //
-     if ( !row ) throw -3;
+     if ( chewbacca ){
-     idtimesync = (UInt_t)atoll(row->GetField(0));
+         oss.str("");
-   } else {
+         oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;";
-     oss.str("");
+         if ( debug ) printf(" %s \n",oss.str().c_str());
-     if ( STATIC ){
+         result = conn->Query(oss.str().c_str());
-       oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
+         //
-           << " LEFT JOIN GL_ROOT "
+         if ( !result ) throw -3;
-           << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
+         //
-           << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "
+         row = result->Next();
-           << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
+         //
+         if ( !row ) throw -3;
+         idtimesync = (UInt_t)atoll(row->GetField(0));
      } else {
-       oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
+         oss.str("");
-           << " LEFT JOIN GL_ROOT "
+         if ( STATIC ){
-           << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
+             oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
-           << " WHERE GL_RAW.PATH = '$PAM_RAW' AND "
+                 << " LEFT JOIN GL_ROOT "
-           << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
+                 << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
+                 << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "
+                 << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
+         } else {
+             oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
+                 << " LEFT JOIN GL_ROOT "
+                 << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
+                 << " WHERE GL_RAW.PATH = '$PAM_RAW' AND "
+                 << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
+         };
+         result = conn->Query(oss.str().c_str());
+         //
+         if ( !result ) throw -12;
+         //
+         row = result->Next();
+         //
+         if ( !row ) throw -10;
+         if ( row != NULL && (UInt_t)atoll(row->GetField(0))>0 ){
+             idroot = (UInt_t)atoll(row->GetField(2));
+             delete row;
+             delete result;
+             return(1);
+         };
+         //
+         // determine which timesync has to be used
+         //
+         oss.str("");
+         oss << "SELECT GL_TIMESYNC.ID FROM GL_TIMESYNC LEFT JOIN GL_RAW ON GL_RAW.ID = GL_TIMESYNC.ID_RAW ORDER BY GL_TIMESYNC.ID DESC LIMIT 1;";
+         result = conn->Query(oss.str().c_str());
+         //
+         if ( !result ) throw -3;
+         //
+         row = result->Next();
+         //
+         if ( !row ) throw -3;
+         idtimesync = (UInt_t)atoll(row->GetField(0));
      };
-     result = conn->Query(oss.str().c_str());
-     //
+     delete row;
-     if ( !result ) throw -12;
+     delete result;
-     //
-     row = result->Next();
+     // ----------------------
-     //
+     // insert root file
-     if ( !row ) throw -10;
+     // ----------------------
-     if ( row != NULL && (UInt_t)atoll(row->GetField(0))>0 ){
-       idroot = (UInt_t)atoll(row->GetField(2));
+     GL_ROOT glroot = GL_ROOT();
-       return(1);
+     glroot.ID_RAW = GetID_RAW();
+     glroot.ID_TIMESYNC = idtimesync;
+     if ( STATIC ){
+       glroot.PATH = GetRootPath();
+     } else {
+       if ( KEEPENV ){
+         glroot.PATH = gSystem->DirName(filerootname.Data());
+       } else {
+         glroot.PATH = "$PAM_L0";
+       };
      };
-     //
+     glroot.NAME = GetRootFile();
-     // determine which timesync has to be used
-     //
+     if ( insertPamelaRootFile(&glroot) )return 1;
-     oss.str("");
-     oss << "SELECT GL_TIMESYNC.ID FROM GL_TIMESYNC LEFT JOIN GL_RAW ON GL_RAW.ID = GL_TIMESYNC.ID_RAW ORDER BY GL_TIMESYNC.ID DESC LIMIT 1;";
+     SetID_ROOT(glroot.ID);
-     result = conn->Query(oss.str().c_str());
-     //
-     if ( !result ) throw -3;
+     return (0);
-     //
+ }
-     row = result->Next();
+ /**
-     //
+  * Insert all the new rows into GL_ROOT.
-     if ( !row ) throw -3;
+  * The raw file indicates in the parameters should be already been stored in the database.
-     idtimesync = (UInt_t)atoll(row->GetField(0));
+  */
-   };
+ Int_t PamelaDBOperations::insertPamelaRootFile(GL_ROOT *glroot){
+   stringstream oss;
+   TSQLResult *result = 0;
+   TSQLRow    *row    = 0;
+   //
    //
    oss.str("");
-   if ( STATIC ){
+   oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
-     oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
+       << glroot->ID_RAW << "', '" << glroot->ID_TIMESYNC << "', '" << glroot->PATH << "', '" << glroot->NAME << "')";
-         << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";
-   } else {
-     oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
-         << this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')";
-   };
    //
-   if ( debug ) printf(" query is %s \n",oss.str().c_str());
+   if ( debug ) printf("%s \n",oss.str().c_str());
    if (conn->Query(oss.str().c_str()) == 0) throw -4;
    //
    delete result;
    //
    oss.str("");
    //  oss << "SELECT ID FROM GL_ROOT WHERE ID_RAW=" << this->GetID_RAW() << ";";
-   if ( STATIC ){
+   oss << "SELECT ID FROM GL_ROOT WHERE PATH='" << glroot->PATH << "' and NAME='"<< glroot->NAME <<"';";
-     oss << "SELECT ID FROM GL_ROOT WHERE PATH='" << this->GetRootPath().Data() << "' and NAME='"<< this->GetRootFile().Data() <<"';";
-   } else {
-     oss << "SELECT ID FROM GL_ROOT WHERE PATH='$PAM_L0' and NAME='"<< this->GetRootFile().Data() <<"';";
-   };
    //
-   if ( debug ) printf(" query is %s \n",oss.str().c_str());
+   if ( debug ) printf("%s \n",oss.str().c_str());
    result = conn->Query(oss.str().c_str());
    if ( !result ) throw -12;
    row = result->Next();
    if ( !row ) throw -3;
-   this->SetID_ROOT((UInt_t)atoll(row->GetField(0)));
-   if ( debug ) printf(" The ID of the ROOT file is %u \n",this->GetID_ROOT());
+   glroot->ID = (UInt_t)atoll(row->GetField(0));
+   if ( debug ) printf(" The ID of the ROOT file is %u \n",glroot->ID);
    //
    delete result;
    //
-Line 3250 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3343 
 Bool_t PamelaDBOperations::IsRunConsiste
              //      this->HandleRunFragments(true,mistrail,firstentry,lastentry); // cannot call it here since it enters a loop which will destroy the already stored variables if we arrive here from HandleRunFragments
              //
              mishead = true;
-Line 3581 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3662 
 Bool_t PamelaDBOperations::IsRunConsiste
              }; // EEE
            } else {
              if ( !IsRunAlreadyInserted() ){
                glrun->SetID(this->AssignRunID());
-Line 4446 
 Int_t PamelaDBOperations::insertCALOPULS
+Line 4496 
 Int_t PamelaDBOperations::insertCALOPULS
   * Fill the GL_TRK_CALIB table
   */
  void PamelaDBOperations::HandleTRK_CALIB(Bool_t pk1, Bool_t pk2){
-   //
-   TSQLResult *result = 0;
+     GL_TRK_CALIB *glcal = new GL_TRK_CALIB();
-   TSQLRow    *row    = 0;
-   //
-   stringstream oss;
-   oss.str("");
-   //
-   UInt_t totime = 0;
-   //
-   if ( !pk1 && !pk2 ){
-     if ( IsDebug() ) printf(" Cannot handle trk calibration with both packet missing!\n");
-     return;
-   };
-   //
-   // check if the calibration has already been inserted
-   //
-   oss.str("");
-   oss << " SELECT ID FROM GL_TRK_CALIB WHERE "
-       << " BOOT_NUMBER = "<< this->GetBOOTnumber(); //
-   oss << " AND ( ( ";
-   if ( pk1 ){
-     oss << " OBT1 = "<< obt1 << " AND "
-         << " PKT1 = "<< pkt1
-         << " ) OR ( ";
-   } else {
-     oss << " PKT1 = "<< pkt2-1
-         << " ) OR ( ";
-   };
-   if ( pk2 ){
-     oss << " OBT2 = "<< obt2 << " AND "
-         << " PKT2 = "<< pkt2;
-   } else {
-     oss << " PKT2 = "<< pkt1+1;
-   };
-   oss << " ) );";
-   //
-   if ( IsDebug() ) printf(" Check if the trk calibration has already 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 ){
-     //
-     if ( IsDebug() ) printf(" Trk calibration already inserted in the DB\n");
-     if ( PEDANTIC ) throw -80;
      //
-   } else {
+     glcal->ID         = 0;
+     glcal->ID_ROOT_L0 = GetID_ROOT();
+     glcal->EV_ROOT_CALIBTRK1 = t1;
+     glcal->EV_ROOT_CALIBTRK2 = t2;
+     glcal->FROM_TIME = fromtime;
+     glcal->TO_TIME   = 0;
+     glcal->OBT1      = obt1;
+     glcal->OBT2      = obt2;
+     glcal->PKT1      = pkt1;
+     glcal->PKT2      = pkt2;
+     glcal->BOOT_NUMBER = GetBOOTnumber();
+     glcal->VALIDATION = valid;
      //
-     // we have to insert a new calibration, check where to place it
+     HandleTRK_CALIB(glcal);
      //
-     oss.str("");
+     delete glcal;
-     oss << " SELECT ID,TO_TIME FROM GL_TRK_CALIB WHERE "
+ }
-         << " FROM_TIME < "<< fromtime << " AND "
+ /**
-         << " TO_TIME > "<< fromtime << ";";
+  * Fill the GL_TRK_CALIB table
+  */
+ void PamelaDBOperations::HandleTRK_CALIB(GL_TRK_CALIB *glcal){
+     Bool_t pk1 = (glcal->OBT1>0&&glcal->PKT1>0);
+     Bool_t pk2 = (glcal->OBT2>0&&glcal->PKT2>0);
+     UInt_t boot_number = glcal->BOOT_NUMBER;
+     UInt_t obt1 = glcal->OBT1;
+     UInt_t obt2 = glcal->OBT2;
+     UInt_t pkt1 = glcal->PKT1;
+     UInt_t pkt2 = glcal->PKT2;
+     UInt_t fromtime = glcal->FROM_TIME;
+     UInt_t totime = 0;
+     UInt_t idroot = glcal->ID_ROOT_L0;
+     UInt_t t1 = glcal->EV_ROOT_CALIBTRK1;
+     UInt_t t2 = glcal->EV_ROOT_CALIBTRK2;
+     UInt_t valid = glcal->VALIDATION;
      //
-     if ( IsDebug() ) printf(" Check where to place the trk calibration: query is \n %s \n",oss.str().c_str());
+     TSQLResult *result = 0;
-     result = conn->Query(oss.str().c_str());
+     TSQLRow    *row    = 0;
      //
-     if ( !result ) throw -4;
+     stringstream oss;
+     oss.str("");
      //
-     row = result->Next();
      //
-     if ( !row ){
+     if ( !pk1 && !pk2 ){
-       //
+         if ( IsDebug() ) printf(" Cannot handle trk calibration with both packet missing!\n");
-       // no calibrations in the db contain our calibration
+         return;
-       //
-       if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n");
-       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 "
-           << " FROM_TIME > "<< fromtime << " ORDER BY FROM_TIME ASC LIMIT 1;";
-       //
-       if ( IsDebug() ) printf(" Check the upper limit for calibration: 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 ){
-         totime = numeric_limits<UInt_t>::max();
-       } else {
-         totime = (UInt_t)atoll(row->GetField(0));
-       };
-       //
-     } else {
-       //
-       // determine upper and lower limits and make space for the new calibration
-       //
-       totime = (UInt_t)atoll(row->GetField(1));
-       //
-       oss.str("");
-       oss << " UPDATE GL_TRK_CALIB SET "
-           << " TO_TIME = "<< fromtime << " WHERE "  // NOTICE: to_time is equal to from_time of the calibration before, so the interval is: [from_time,to_time[
-           << " ID = "<< row->GetField(0) << ";";
-       //
-       if ( IsDebug() ) printf(" Make space for the new trk calibration: query is \n %s \n",oss.str().c_str());
-       result = conn->Query(oss.str().c_str());
-       //
-       if ( !result ) throw -4;
-       //
      };
      //
-     oss.str("");
+     // check if the calibration has already been inserted
-     oss << " INSERT INTO GL_TRK_CALIB (ID,ID_ROOT_L0,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,FROM_TIME,TO_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION) "
-         << " VALUES (NULL,' "
-         << idroot << "',";
      //
-     if ( !pk1 ){
+     oss.str("");
-       oss << "NULL,";
+     oss << " SELECT ID FROM GL_TRK_CALIB WHERE "
+         << " BOOT_NUMBER = "<< boot_number; //
+     oss << " AND FROM_TIME="<<fromtime; /// NEWNEWNEW -- VA BENE ?!?!?!?!
+     oss << " AND ( ( ";
+     if ( pk1 ){
+         oss << " OBT1 = "<< obt1 << " AND "
+             << " PKT1 = "<< pkt1
+             << " ) OR ( ";
      } else {
-       oss << "'"
+         oss << " PKT1 = "<< pkt2-1
-           << t1 << "',";
+             << " ) OR ( ";
      };
-     //
+     if ( pk2 ){
-     if ( !pk2 ){
+         oss << " OBT2 = "<< obt2 << " AND "
-       oss << "NULL,'";
+             << " PKT2 = "<< pkt2;
      } else {
-       oss << "'"
+         oss << " PKT2 = "<< pkt1+1;
-           << t2 << "','";
+     };
-     };
+     oss << " ) );";
-     //
-     oss << fromtime << "','"
-         << totime << "','"
-         << obt1 << "','"
-         << pkt1 << "','"
-         << obt2 << "','"
-         << pkt2 << "','"
-         << this->GetBOOTnumber() << "','"
-         << valid << "');";
-     //
-     if ( IsDebug() ) printf(" Insert the new trk calibration: query is \n %s \n",oss.str().c_str());
      //
+     if ( IsDebug() ) printf(" Check if the trk calibration has already 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();
-   //
- };
- /**
-  * Scan tracker calibrations packets, fill the GL_TRK_CALIB table
-  */
- Int_t PamelaDBOperations::insertTRK_CALIB(){
-   //
-   CalibTrk1Event *caltrk1 = 0;
-   CalibTrk2Event *caltrk2 = 0;
-   TTree *tr1 = 0;
-   TTree *tr2 = 0;
-   EventHeader *eh1 = 0;
-   PscuHeader *ph1 = 0;
-   EventHeader *eh2 = 0;
-   PscuHeader *ph2 = 0;
-   //
-   PacketType *pctp=0;
-   EventCounter *codt2=0;
-   //
-   Int_t nevents1 = 0;
-   Int_t nevents2 = 0;
-   //
-   fromtime = 0;
-   //
-   obt1 = 0;
-   pkt1 = 0;
-   obt2 = 0;
-   pkt2 = 0;
-   //
-   tr1 = (TTree*)file->Get("CalibTrk1");
-   if ( !tr1 || tr1->IsZombie() ) throw -22;
-   tr2 = (TTree*)file->Get("CalibTrk2");
-   if ( !tr2 || tr2->IsZombie() ) throw -23;
-   //
-   tr1->SetBranchAddress("CalibTrk1", &caltrk1);
-   tr1->SetBranchAddress("Header", &eh1);
-   nevents1 = tr1->GetEntries();
-   tr2->SetBranchAddress("CalibTrk2", &caltrk2);
-   tr2->SetBranchAddress("Header", &eh2);
-   nevents2 = tr2->GetEntries();
-   //
-   if ( !nevents1 && !nevents2 ) return(1);
-   //
-   t2 = -1;
-   Int_t pret2 = 0;
-   Int_t t2t1cal = 0;
-   //
-   for (t1=0; t1 < nevents1; t1++){
-     //
-     pret2 = t2;
-     tr1->GetEntry(t1);
-     //
-     ph1 = eh1->GetPscuHeader();
-     obt1 = ph1->GetOrbitalTime();
-     pkt1 = ph1->GetCounter();
-     fromtime = this->GetAbsTime(ph1->GetOrbitalTime());
-     //
-     //     valid = 1;
-     //     //
-     //     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->PKT(pkt1) <= upperpkt && this->OBT(obt1) >= this->OBT(obtfirst) && this->OBT(obt1) <= upperobt ){
+     if ( row ){
-       //    if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
-       //
-       if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
-       //
-       valid = ValidateTrkCalib( caltrk1, eh1 );
-       if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
-       //
-       // Do we have the second calibration packet?
-       //
-       while ( t2t1cal < t1+1 ){ // get the calibration packet2 that follows the packet1
-         //
-         t2++;
          //
-         pret2 = t2 - 1; // EMILIANO
+         if ( IsDebug() ) printf(" Trk calibration already inserted in the DB\n");
+         if ( PEDANTIC ) throw -80;
          //
-         if ( t2 < nevents2 ){
+     } else {
-           tr2->GetEntry(t2);
-           codt2 = eh2->GetCounter();
-           t2t1cal = codt2->Get(pctp->CalibTrk1);
-           //
-           ph2 = eh2->GetPscuHeader();
-           obt2 = ph2->GetOrbitalTime();
-           pkt2 = ph2->GetCounter();
-           //
-           if ( IsDebug() ) printf(" This is a trk calibration2 at obt %u pkt %u t2 is %u , t2t1cal is %u \n",obt2,pkt2,t2,t2t1cal);
-           //      if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2
-           //
-         } else {
-           //
-           // running out of vector without finding the corresponding calibration, sig
-           //
-           if ( IsDebug() ) printf(" t2 >= nevents2 \n");
-           pret2 = t2;
-           obt2 = 0;
-           //      pkt2 = pkt1+2;
-           pkt2 = 0;
-           t2t1cal = t1+1;
-         };
-         //      if ( (this->PKT(pkt2) < this->PKT(pktfirst) || this->PKT(pkt2) > upperpkt) && (this->OBT(obt2) < this->OBT(obtfirst) || this->OBT(obt2) > upperobt) ){
-         // EMILIANO
-         //      if ( (this->PKT(pkt2) < this->PKT(pktfirst) || this->PKT(pkt2) > upperpkt) || (this->OBT(obt2) < this->OBT(obtfirst) || this->OBT(obt2) > upperobt) ){
-         //        //    if ( this->PKT(pkt2) < this->PKT(pktfirst) && this->OBT(obt2) < this->OBT(obtfirst) ){
-         //        if ( IsDebug() ) printf(" running out of vector without finding the corresponding calibration, sig \n");
-         //        //
-         //        // running out of vector without finding the corresponding calibration, sig
-         //        //
-         //        pret2 = t2;
-         //        obt2 = 0;
-         //        //      pkt2 = pkt1+2;
-         //        pkt2 = 0;
-         //        t2t1cal = t1+1;
-         //      };
-         //
-       };
-       //
-       if ( IsDebug() ) printf(" Check if trk calibration2 exists \n");
-       //
-       // EMILIANO
-       if ( (this->PKT(pkt2) < this->PKT(pktfirst) || this->PKT(pkt2) > upperpkt) || (this->OBT(obt2) < this->OBT(obtfirst) || this->OBT(obt2) > upperobt) ){
-         //      if ( this->PKT(pkt2) < this->PKT(pktfirst) && this->OBT(obt2) < this->OBT(obtfirst) ){
-         if ( IsDebug() ) printf(" The calibration found is outside the interval, sig \n");
-         //
-         // running out of vector without finding the corresponding calibration, sig
-         //
-         pret2 = t2;
-         obt2 = 0;
-         pkt2 = 0;
-       };
-       //
-       // The calibration is good
-       //
-       if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){
          //
-         if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2);
+         // we have to insert a new calibration, check where to place it
-         if ( IsDebug() ) printf(" The trk calibration2 at obt %u pkt %u t2 is %u is good \n",obt2,pkt2,t2);
          //
-         UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
+         oss.str("");
-         if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
+         oss << " SELECT ID,TO_TIME FROM GL_TRK_CALIB WHERE "
-         valid = valid & valid2;
+             << " FROM_TIME < "<< fromtime << " AND "
+             << " TO_TIME > "<< fromtime << ";";
          //
-         // Handle good calib
+         if ( IsDebug() ) printf(" Check where to place the trk calibration: query is \n %s \n",oss.str().c_str());
+         result = conn->Query(oss.str().c_str());
          //
-         this->HandleTRK_CALIB(true,true);
+         if ( !result ) throw -4;
          //
-         // Check for missing calibtrk1
+         row = result->Next();
          //
-         if ( t2 != pret2+1 ){
+         if ( !row ){
-           //
-           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);
-           //
-           while ( t2 > pret2+1 ){
              //
-             // handle missing calib1
+             // no calibrations in the db contain our calibration
              //
-             pret2++;
+             if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n");
+             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
              //
-             obt1 = 0;
+             oss.str("");
-             pkt1 = 0;
+             oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "
+                 << " FROM_TIME > "<< fromtime << " ORDER BY FROM_TIME ASC LIMIT 1;";
              //
-             tr2->GetEntry(pret2);
+             if ( IsDebug() ) printf(" Check the upper limit for calibration: query is \n %s \n",oss.str().c_str());
-             ph2 = eh2->GetPscuHeader();
+             result = conn->Query(oss.str().c_str());
-             obt2 = ph2->GetOrbitalTime();
-             pkt2 = ph2->GetCounter();
              //
-             fromtime = this->GetAbsTime(ph2->GetOrbitalTime());
+             if ( !result ) throw -4;
              //
-             valid = 0;
+             row = result->Next();
-             this->HandleTRK_CALIB(false,true);
+             if ( !row ){
+                 totime = numeric_limits<UInt_t>::max();
+             } else {
+                 totime = (UInt_t)atoll(row->GetField(0));
+             };
+             //
+         } else {
+             //
+             // determine upper and lower limits and make space for the new calibration
+             //
+             totime = (UInt_t)atoll(row->GetField(1));
+             //
+             oss.str("");
+             oss << " UPDATE GL_TRK_CALIB SET "
+                 << " TO_TIME = "<< fromtime << " WHERE "  // NOTICE: to_time is equal to from_time of the calibration before, so the interval is: [from_time,to_time[
+                 << " ID = "<< row->GetField(0) << ";";
+             //
+             if ( IsDebug() ) printf(" Make space for the new trk calibration: query is \n %s \n",oss.str().c_str());
+             result = conn->Query(oss.str().c_str());
+             //
+             if ( !result ) throw -4;
              //
-           };
-           //
          };
          //
-       } else if ( this->PKT(pkt2) > this->PKT(pkt1)+1 ){
+         oss.str("");
+         oss << " INSERT INTO GL_TRK_CALIB (ID,ID_ROOT_L0,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,FROM_TIME,TO_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION) "
+             << " VALUES (NULL,' "
+             << idroot << "',";
          //
-         // Check for missing calibtrk2
+         if ( !pk1 ){
+             oss << "NULL,";
+         } else {
+             oss << "'"
+                 << t1 << "',";
+         };
          //
-         if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
+         if ( !pk2 ){
-         t2 = pret2;
+             oss << "NULL,'";
+         } else {
+             oss << "'"
+                 << t2 << "','";
+         };
          //
-         // handle missing calib2
+         oss << fromtime << "','"
+             << totime << "','"
+             << obt1 << "','"
+             << pkt1 << "','"
+             << obt2 << "','"
+             << pkt2 << "','"
+             << boot_number << "','"
+             << valid << "');";
          //
-         obt2 = 0;
+         if ( IsDebug() ) printf(" Insert the new trk calibration: query is \n %s \n",oss.str().c_str());
-         pkt2 = 0;
+         //
-         valid = 0;
+         result = conn->Query(oss.str().c_str());
-         this->HandleTRK_CALIB(true,false);
+         //
+         if ( !result ) throw -4;
          //
-       };
-       //
-     } else {
-       //
-       if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u OUTSIDE the considered time interval \n",fromtime,obt1,pkt1);
-       //      if ( PEDANTIC ) throw -79;
-       //
      };
-     //
-   };
+     oss.str("");
-   //
+     oss << " SELECT ID FROM GL_TRK_CALIB ORDER BY ID DESC LIMIT 1 ;";
-   // we have one more calib pkt2 !
+     if ( IsDebug() ) cout << oss.str().c_str() << endl;
-   //
+     result = conn->Query(oss.str().c_str());
-   t2++;
+     if ( !result ) throw -4;;
-   while ( t2 < nevents2 ){
+     row = result->Next();
+     if(row)glcal->ID = (UInt_t)atoll(row->GetField(0));
      //
-     // handle missing calib1
+ };
+ /**
+  * Scan tracker calibrations packets, fill the GL_TRK_CALIB table
+  */
+ Int_t PamelaDBOperations::insertTRK_CALIB(){
+     //
+     CalibTrk1Event *caltrk1 = 0;
+     CalibTrk2Event *caltrk2 = 0;
+     TTree *tr1 = 0;
+     TTree *tr2 = 0;
+     EventHeader *eh1 = 0;
+     PscuHeader *ph1 = 0;
+     EventHeader *eh2 = 0;
+     PscuHeader *ph2 = 0;
+     //
+     PacketType *pctp=0;
+     EventCounter *codt2=0;
+     //
+     Int_t nevents1 = 0;
+     Int_t nevents2 = 0;
+     //
+     fromtime = 0;
      //
-     if ( IsDebug() ) printf(" t2 is %u nevents2 is %u \n",t2,nevents2);
      obt1 = 0;
      pkt1 = 0;
+     obt2 = 0;
+     pkt2 = 0;
      //
-     tr2->GetEntry(t2);
+     tr1 = (TTree*)file->Get("CalibTrk1");
-     ph2 = eh2->GetPscuHeader();
+     if ( !tr1 || tr1->IsZombie() ) throw -22;
-     obt2 = ph2->GetOrbitalTime();
+     tr2 = (TTree*)file->Get("CalibTrk2");
-     pkt2 = ph2->GetCounter();
+     if ( !tr2 || tr2->IsZombie() ) throw -23;
      //
-     fromtime = this->GetAbsTime(ph2->GetOrbitalTime());
+     tr1->SetBranchAddress("CalibTrk1", &caltrk1);
-     valid = 0;
+     tr1->SetBranchAddress("Header", &eh1);
-     //  if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->PKT(pkt1) <= upperpkt && this->OBT(obt1) >= this->OBT(obtfirst) && this->OBT(obt1) <= upperobt ){
+     nevents1 = tr1->GetEntries();
-     // EMILIANO
+     tr2->SetBranchAddress("CalibTrk2", &caltrk2);
-     if ( this->PKT(pkt2) >= this->PKT(pktfirst) && this->PKT(pkt2 <= upperpkt) && this->OBT(obt2) >= this->OBT(obtfirst) && this->OBT(obt2) <= upperobt ){
+     tr2->SetBranchAddress("Header", &eh2);
-       //    if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
+     nevents2 = tr2->GetEntries();
-       //
+     //
-       if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
+     if ( !nevents1 && !nevents2 ) return(1);
-       //
+     //
-       this->HandleTRK_CALIB(false,true);
+     t2 = -1;
-       //
+     Int_t pret2 = 0;
-     };
+     Int_t t2t1cal = 0;
+     //
+     bool MISSING_pkt1 = true;
+     bool MISSING_pkt2 = true;
+     int  ncalib = 0;
+     bool try_to_recover = false;
      //
-     t2++;
+     for (t1=0; t1 < nevents1; t1++){//loop over packet1
+         //
+         pret2 = t2;
+         tr1->GetEntry(t1);
+         //
+         ph1 = eh1->GetPscuHeader();
+         obt1 = ph1->GetOrbitalTime();
+         pkt1 = ph1->GetCounter();
+         fromtime = GetAbsTime(ph1->GetOrbitalTime());
+         //
+         // chek if the packet number and obt are consistent with the other packets ???
+         //
+         if ( PKT(pkt1) >= PKT(pktfirst) && PKT(pkt1) <= upperpkt && OBT(obt1) >= OBT(obtfirst) && OBT(obt1) <= upperobt ){
+             //    if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
+             //
+             if ( IsDebug() ) printf("\n Trk calibration1 %u at time %u obt %u pkt %u \n",t1,fromtime,obt1,pkt1);
+             //
+             valid = ValidateTrkCalib( caltrk1, eh1 );
+             if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
+             //
+             // Do we have the second calibration packet?
+             //
+             if ( IsDebug() ) cout << " Loop over calibration2 to search associated calibration: "<<endl;
+             while ( t2t1cal < t1+1 ){ // get the calibration packet2 that follows the packet1
+                 //
+                 t2++;
+                 //
+                 pret2 = t2 - 1; // EMILIANO
+                 //
+                 if ( t2 < nevents2 ){
+                     tr2->GetEntry(t2);
+                     codt2 = eh2->GetCounter();
+                     t2t1cal = codt2->Get(pctp->CalibTrk1);
+                     //
+                     ph2 = eh2->GetPscuHeader();
+                     obt2 = ph2->GetOrbitalTime();
+                     pkt2 = ph2->GetCounter();
+                     //
+                     if ( IsDebug() ) printf(" >> trk calibration2 at obt %u pkt %u t2 is %u , t2t1cal is %u \n",obt2,pkt2,t2,t2t1cal);
+                     //    if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2
+                     //
+                 } else {
+                     //
+                     // running out of vector without finding the corresponding calibration, sig
+                     //
+                     if ( IsDebug() ) printf(" t2 >= nevents2 \n");
+                     pret2 = t2;
+                     obt2 = 0;
+                     //    pkt2 = pkt1+2;
+                     pkt2 = 0;
+                     t2t1cal = t1+1;
+                 };
+                 //      if ( (this->PKT(pkt2) < this->PKT(pktfirst) || this->PKT(pkt2) > upperpkt) && (this->OBT(obt2) < this->OBT(obtfirst) || this->OBT(obt2) > upperobt) ){
+                 // EMILIANO
+                 //      if ( (this->PKT(pkt2) < this->PKT(pktfirst) || this->PKT(pkt2) > upperpkt) || (this->OBT(obt2) < this->OBT(obtfirst) || this->OBT(obt2) > upperobt) ){
+                 //        //    if ( this->PKT(pkt2) < this->PKT(pktfirst) && this->OBT(obt2) < this->OBT(obtfirst) ){
+                 //        if ( IsDebug() ) printf(" running out of vector without finding the corresponding calibration, sig \n");
+                 //        //
+                 //        // running out of vector without finding the corresponding calibration, sig
+                 //        //
+                 //        pret2 = t2;
+                 //        obt2 = 0;
+                 //        //      pkt2 = pkt1+2;
+                 //        pkt2 = 0;
+                 //        t2t1cal = t1+1;
+                 //      };
+                 //
+             };
+             //
+             if ( IsDebug() ) printf(" Check if trk calibration2 is the right one \n");
+             //
+             // EMILIANO
+             if ( ( PKT(pkt2) < PKT(pktfirst) || PKT(pkt2) > upperpkt) || (OBT(obt2) < OBT(obtfirst) || OBT(obt2) > upperobt) ){
+                 //      if ( this->PKT(pkt2) < this->PKT(pktfirst) && this->OBT(obt2) < this->OBT(obtfirst) ){
+                 if ( IsDebug() ) printf(" *WARNING*  The calibration found is outside the interval, sig \n");
+                 //
+                 // running out of vector without finding the corresponding calibration, sig
+                 //
+                 pret2 = t2;
+                 obt2 = 0;
+                 pkt2 = 0;
+             };
+             if ( PKT(pkt2) == PKT(pkt1)+1 ){
+                 if ( IsDebug() ) cout << " ...OK"<<endl;
+                 // =======================
+                 // The calibration is good
+                 // =======================
+                 //
+ //      if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2);
+ //      if ( IsDebug() ) printf(" Trk calibration2 at obt %u pkt %u t2 is %u is good \n",obt2,pkt2,t2);
+ //        if ( IsDebug() ) printf("\n Trk calibration2 at time %u obt %u pkt %u \n",fromtime,obt2,pkt2);
+                 if ( IsDebug() ) printf(" Trk calibration2 %u at time %u obt %u pkt %u \n",t2,fromtime,obt2,pkt2);
+                 //
+                 UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
+                 if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
+ //      valid = valid & valid2;
+                 valid = valid & valid2; //QUESTO VA CAMBIATO
+                 //
+                 // Handle good calib
+                 //
+                 MISSING_pkt1 = false;
+                 MISSING_pkt2 = false;
+ //              this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+                 //
+                 // Check for missing calibtrk1
+                 //
+                 if ( t2 != pret2+1 ){
+                     //
+                     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);
+                     //
+                     while ( t2 > pret2+1 ){
+                         //
+                         // handle missing calib1
+                         //
+                         pret2++;
+                         //
+                         obt1 = 0;
+                         pkt1 = 0;
+                         //
+                         tr2->GetEntry(pret2);
+                         ph2 = eh2->GetPscuHeader();
+                         obt2 = ph2->GetOrbitalTime();
+                         pkt2 = ph2->GetCounter();
+                         //
+                         fromtime = this->GetAbsTime(ph2->GetOrbitalTime());
+                         //
+                         valid = 0;
+                         MISSING_pkt1 = true;
+                         MISSING_pkt2 = false;
+ //                      this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+                         //
+                     };
+                     //
+                 };
+                 //
+             } else if ( this->PKT(pkt2) > this->PKT(pkt1)+1 ){
+                 //
+                 // Check for missing calibtrk2
+                 //
+                 if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
+                 t2 = pret2;
+                 //
+                 // handle missing calib2
+                 //
+                 obt2 = 0;
+                 pkt2 = 0;
+                 valid = 0;
+                 MISSING_pkt1 = false;
+                 MISSING_pkt2 = true;
+ //              this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+                 //
+             };
+             //
+             if( !(MISSING_pkt1&MISSING_pkt2) ){
+               this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+               ncalib++;
+               if( MISSING_pkt1||MISSING_pkt2||!valid )try_to_recover=true;
+             }
+         } else {
+           //
+             if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u OUTSIDE the considered time interval \n",fromtime,obt1,pkt1);
+             //      if ( PEDANTIC ) throw -79;
+             //
+         };
+         //
+     }; //end loop on pkt1
      //
-   };
+     // we have one more calib pkt2 !
+     //
+     t2++;
+     while ( t2 < nevents2 ){
+         //
+         // handle missing calib1
+         //
+         if ( IsDebug() ) printf(" t2 is %u nevents2 is %u \n",t2,nevents2);
+         obt1 = 0;
+         pkt1 = 0;
+         //
+         tr2->GetEntry(t2);
+         ph2 = eh2->GetPscuHeader();
+         obt2 = ph2->GetOrbitalTime();
+         pkt2 = ph2->GetCounter();
+         //
+         fromtime = this->GetAbsTime(ph2->GetOrbitalTime());
+         valid = 0;
+         //  if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->PKT(pkt1) <= upperpkt && this->OBT(obt1) >= this->OBT(obtfirst) && this->OBT(obt1) <= upperobt ){
+         // EMILIANO
+         if ( this->PKT(pkt2) >= this->PKT(pktfirst) && this->PKT(pkt2 <= upperpkt) && this->OBT(obt2) >= this->OBT(obtfirst) && this->OBT(obt2) <= upperobt ){
+             //    if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
+             //
+             if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
+             //
+             MISSING_pkt1 = true;
+             MISSING_pkt2 = false;
+             this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+             ncalib++;
+             if( MISSING_pkt1||MISSING_pkt2||!valid )try_to_recover=true;
+             //
+         };
+         //
+         t2++;
+         //
+     };
+ // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ // -----------------------------------------------------------------
+ // in case of corruption, check if the calibration can be recovered
+ // from another chewbacca file
+ // -----------------------------------------------------------------
+ // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ //    cout <<" TRY TO RECOVER ?? "<<try_to_recover<<endl;
+     if(chewbacca&&try_to_recover){
+         if ( IsDebug() ) cout << endl << ">>>> TRY TO RECOVER TRACKER CALIBRATIONS <<<<"<<endl;
+         TSQLResult *result = 0;
+         TSQLRow    *row    = 0;
+         //
+         stringstream oss;
+         oss.str("");
+         //
+         ////////////////////////////////////////////////////////////////////////
+         // retrieve the name of the current file:
+         ////////////////////////////////////////////////////////////////////////
+         oss.str("");
+         oss << "SELECT NAME FROM GL_ROOT where ID=" << GetID_ROOT() <<";";
+         if ( IsDebug() ) cout << oss.str().c_str() << endl;
+         result = conn->Query(oss.str().c_str());
+         if ( !result ) throw -4;;
+         row = result->Next();
+         TString thisfilename = (TString)row->GetField(0);
+         if ( IsDebug() ) cout << "Current file ==> "<<thisfilename<<endl;
+         ////////////////////////////////////////////////////////////////////////
+         // read all the calibrations inserted
+         ////////////////////////////////////////////////////////////////////////
+         oss.str("");
+         oss << " SELECT ";
+         oss << " ID,FROM_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,TO_TIME";
+         oss << " FROM GL_TRK_CALIB ";
+         oss << " ORDER BY ID DESC LIMIT "<<ncalib<<"; ";
+         if ( IsDebug() ) cout << oss.str().c_str() << endl;
+         result = conn->Query(oss.str().c_str());
+         if ( !result ) throw -4;;
+         if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl;
+         // -----------------------------------
+         // loop over calibrations ...
+         // -----------------------------------
+         UInt_t nn=0;
+         do {
+             row = result->Next();
+             if(!row)break;
+             UInt_t id       = (UInt_t)atoll(row->GetField(0));
+             UInt_t fromtime = (UInt_t)atoll(row->GetField(1));
+             UInt_t obt1     = (UInt_t)atoll(row->GetField(2));
+             UInt_t pkt1     = (UInt_t)atoll(row->GetField(3));
+             UInt_t obt2     = (UInt_t)atoll(row->GetField(4));
+             UInt_t pkt2     = (UInt_t)atoll(row->GetField(5));
+             UInt_t boot     = (UInt_t)atoll(row->GetField(6));
+             UInt_t valid    = (UInt_t)atoll(row->GetField(7));
+             bool MISSING_pkt1 = (row->GetFieldLength(8)==0);
+             bool MISSING_pkt2 = (row->GetFieldLength(9)==0);
+             UInt_t totime   = (UInt_t)atoll(row->GetField(10));
+             // -------------------------------------
+             // ...check if the entry is corrupted...
+             // -------------------------------------
+             cout <<"*** "<< MISSING_pkt1 << MISSING_pkt2 << valid <<endl;
+             bool CORRUPTED = (MISSING_pkt1||MISSING_pkt2||!valid);
+             if ( IsDebug() ) cout << "("<<nn<<")  ID = "<<id<<" from GL_TRK_CALIB  ==>  corrupted ? "<<CORRUPTED<<endl;
+ //          if( !CORRUPTED  )continue; // nothing to do
+             /////////////////////////////////////////////////////////
+             // if it is corrupted, ...look for ather chewbacca files
+             // containing the same calibrations ...
+             /////////////////////////////////////////////////////////
+             bool this_MISSING_pkt1 = false;
+             bool this_MISSING_pkt2 = false;
+             int  this_t1=0;
+             int  this_t2=0;;
+             UInt_t this_valid = 0;
+             TString path       = "";
+             TString name       = "";
+             TString raw        = "";
+             UInt_t obt0        = 0;
+             UInt_t timesync    = 0;
+             UInt_t boot_number = 0;
+             bool   FOUND       = false;
+             if ( IsDebug() ) cout << "------------------------------------------------------------" <<endl;
+ //          for(int itable=0; itable<2; itable++){
+             for(int itable=0; itable<1; itable++){
+                 // ------------------------------------------------------
+                 // loop over both ROOT_TABLE and ROOT_TABLE_BAD
+                 // ------------------------------------------------------
+                 TString table = "ROOT_TABLE";
+                 if(itable==1)table = "ROOT_TABLE_BAD";
+                 oss.str("");
+                 oss << " SELECT ";
+                 oss << " FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,BOOT_NUMBER,INPUT_NAME ";
+                 oss << " FROM "<<table;
+                 oss << " WHERE 1 " << endl;
+                 oss << " AND FILE_NAME != \""<< thisfilename<<"\"";
+                 if( !MISSING_pkt1 ){
+                     oss << " AND ";
+                     oss << " PKT_NUMBER_INIT < "<<pkt1;
+                     oss << " AND ";
+                     oss << " PKT_NUMBER_FINAL > "<<pkt1;
+                     oss << " AND ";
+                     oss << " PKT_OBT_INIT < "<<obt1;
+                     oss << " AND ";
+                     oss << " PKT_OBT_FINAL > "<<obt1;
+                 }else{
+                     if(pkt2>1) pkt1 = pkt2-1;//serve dopo
+                 }
+                 if( !MISSING_pkt2 ){
+                     oss << " AND ";
+                     oss << " PKT_NUMBER_INIT < "<<pkt2;
+                     oss << " AND ";
+                     oss << " PKT_NUMBER_FINAL > "<<pkt2;
+                     oss << " AND ";
+                     oss << " PKT_OBT_INIT < "<<obt2;
+                     oss << " AND ";
+                     oss << " PKT_OBT_FINAL > "<<obt2;
+                 }else{
+                     if(pkt1>0) pkt2 = pkt1+1;//serve dopo
+                 }
+                 if( boot> 0 ){
+                     oss << " AND ";
+                     oss << " BOOT_NUMBER = "<<boot;
+                 }else{
+                 }
+                 oss << " ORDER BY BAD_PKT_CALREAD ASC; ";
+                 TSQLResult *result2 = 0;
+                 TSQLRow    *row2    = 0;
+                 if ( IsDebug() ) cout << oss.str().c_str() << endl;
+                 result2 = conn->Query(oss.str().c_str());
+                 if ( !result2 ) throw -4;;
+                 if ( IsDebug() ) cout <<"Rows: "<<result2->GetRowCount()<<endl;
+                 // ------------------------------------------------------
+                 // loop over files containing repetition (if any)
+                 // ------------------------------------------------------
+                 do {
+                     row2 = result2->Next();
+                     if(!row2)break;
+                     // ------------------------------------------------------
+                     // ... a repetition is found ...
+                     // ------------------------------------------------------
+                     path       = (TString)row2->GetField(0);
+                     name       = (TString)row2->GetField(1);
+                     raw        = (TString)row2->GetField(5);
+                     obt0        = (UInt_t)atoll(row2->GetField(2));
+                     timesync    = (UInt_t)atoll(row2->GetField(3));
+                     boot_number = (UInt_t)atoll(row2->GetField(4));
+                     if ( IsDebug() ) cout << "- - - - - - - - - - -" <<endl;
+ //                  cout << path <<endl;
+ //                  cout << "File    : " <<name <<endl;
+ //                  cout << obt0 <<endl;
+ //                  cout << timesync <<endl;
+ //                  cout << "boot n. : "<<boot_number <<endl;
+ //                  cout << raw <<endl;
+                     // ------------------------------------------------------
+                     // ... retrieve the calibration packets.
+                     // ------------------------------------------------------
+                     if ( IsDebug() ) printf(" file is %s/%s \n",((TString)gSystem->ExpandPathName(path.Data())).Data(),name.Data());
+                     TFile *file = new TFile(((TString)gSystem->ExpandPathName(path.Data()))+"/"+name); // EM, path could be symbolic and we must expand it
+                     if(!file)throw -100;
+                     if(file->IsZombie())throw -100;
+                     //
+                     tr1 = (TTree*)file->Get("CalibTrk1");
+                     if ( !tr1 || tr1->IsZombie() ) throw -22;
+                     tr2 = (TTree*)file->Get("CalibTrk2");
+                     if ( !tr2 || tr2->IsZombie() ) throw -23;
+                     //
+                     tr1->SetBranchAddress("CalibTrk1", &caltrk1);
+                     tr1->SetBranchAddress("Header", &eh1);
+                     nevents1 = tr1->GetEntries();
+                     tr2->SetBranchAddress("CalibTrk2", &caltrk2);
+                     tr2->SetBranchAddress("Header", &eh2);
+                     nevents2 = tr2->GetEntries();
+                     for(this_t1=0; this_t1<nevents1; this_t1++){
+                         tr1->GetEntry(this_t1);
+                         if(
+                             (UInt_t)eh1->GetPscuHeader()->GetCounter() == pkt1 &&
+                             true) break;
+                         this_MISSING_pkt1 = true;
+                     }
+                     for(this_t2=0; this_t2<nevents2; this_t2++){
+                         tr2->GetEntry(this_t2);
+                         if(
+                             (UInt_t)eh2->GetPscuHeader()->GetCounter() == pkt2 &&
+                             true) break;
+                         this_MISSING_pkt2 = true;
+                     }
+                     this_valid =
+                         ValidateTrkCalib( caltrk1, eh1, file )
+                         *
+                         ValidateTrkCalib( caltrk2, eh2, file );
+                     // ---------------------------------------------------------------------
+                     // accept the calibration if it is better than the previous:
+                     //
+                     // - if the new calibration is perfect (both valid packets)
+                     // - if the new calibration has both the packets and the previous does not
+                     // ---------------------------------------------------------------------
+                     if(
+                         ( !this_MISSING_pkt1&&!this_MISSING_pkt2&&this_valid )||
+                         ( (MISSING_pkt1||MISSING_pkt2) && (!this_MISSING_pkt1&&!this_MISSING_pkt2) )||
+                         false)FOUND=true;
+                     if(file)file->Close();
+                     if(FOUND)break;
+                 }while(1);//endl loop over root table entries
+                 if(FOUND)break;
+             }//end loop over tables
+             if(FOUND){
+             if ( IsDebug() ) cout << " >>> REPETITION FOUND :-) <<<" <<endl;
+                 ////////////////////////////////////////////
+                 // insert a new entry in GL_TRK_CALIB and
+                 // modify the time-tag of the previous one
+                 ////////////////////////////////////////////
+                 // ---------------------------------------------------------------------
+                 // step 1: insert a new raw file in GL_RAW
+                 // ---------------------------------------------------------------------
+                 //
+                 // check if the raw file already exist
+                 //
+                 UInt_t id_raw = 0;   // EM GL_RAW is there only for backward compatibility so we do not need to fill it when in "chewbacca" mode
+ //              oss.str("");
+ //              oss << "SELECT ID FROM GL_RAW where NAME=\"" << gSystem->BaseName(raw.Data()) <<"\";";
+ //              if ( IsDebug() ) cout << oss.str().c_str() << endl;
+ //              result = conn->Query(oss.str().c_str());
+ //              if ( !result ) throw -4;;
+ //              if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl;
+ //              if( result->GetRowCount() == 0){
+ //                  if ( IsDebug() ) cout << " << Insert new RAW file >> "<<endl;
+ //                  // - - - - - - - - - - -
+ //                  // insert new raw file
+ //                  // - - - - - - - - - - -
+ //                  GL_RAW glraw = GL_RAW();
+ //                  glraw.PATH        = gSystem->DirName(raw.Data());
+ //                  glraw.NAME        = gSystem->BaseName(raw.Data());
+ //                  glraw.BOOT_NUMBER = boot_number;
+ //                  //
+ //                  insertPamelaRawFile( &glraw );
+ //                  //
+ //                  id_raw = glraw.ID;
+ //              }else{
+ //                  row = result->Next();
+ //                  id_raw = (UInt_t)atoll(row->GetField(0));
+ //              }
+ //              if ( IsDebug() ) cout << "ID_RAW = "<<id_raw<<endl;
+                 // ---------------------------------------------------------------------
+                 // step 1(bis): retrieve the timesync id associated to the file
+                 // (NB, uso lo stesso associato al file iniziale)
+                 // ---------------------------------------------------------------------
+                 UInt_t idtimesync = 0;
+                 oss.str("");
+                 oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;";
+                 if ( debug ) printf(" %s \n",oss.str().c_str());
+                 result = conn->Query(oss.str().c_str());
+                 if ( !result ) throw -3;
+                 row = result->Next();
+                 if ( !row ) throw -3;
+                 idtimesync = (UInt_t)atoll(row->GetField(0));
+                 if ( IsDebug() ) cout << "ID_TIMESYNC = "<<idtimesync<<endl;
+                 // ---------------------------------------------------------------------
+                 // step 2: insert a new root file in GL_ROOT
+                 // ---------------------------------------------------------------------
+                 //
+                 // check if the root file already exist
+                 //
+                 UInt_t id_root = 0;
+                 oss.str("");
+                 oss << "SELECT ID FROM GL_ROOT where NAME=\"" << gSystem->BaseName(name.Data()) <<"\";";
+                 if ( IsDebug() ) cout << oss.str().c_str() << endl;
+                 result = conn->Query(oss.str().c_str());
+                 if ( !result ) throw -4;;
+                 if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl;
+                 if( result->GetRowCount() == 0){
+                     if ( IsDebug() ) cout << " << Insert new ROOT file >> "<<endl;
+                     // - - - - - - - - - - -
+                     // insert new root file
+                     // - - - - - - - - - - -
+                     GL_ROOT glroot = GL_ROOT();
+                     glroot.ID_RAW      = id_raw;
+                     glroot.ID_TIMESYNC = idtimesync;
+                     //
+                     // EM STATIC = the full expanded path must be put in the DB, KEEPENV = the path given as input (or found in ROOT_TABLE) must be used,
+                     //             NOT STATIC NOT KEEPENV = $PAM_L0 must be used in the DB
+                     //
+                     if ( STATIC ){
+                       glroot.PATH        = (TString)gSystem->ExpandPathName(path);
+                     } else {
+                       if ( KEEPENV ){
+                         glroot.PATH      = path;
+                       } else {
+                         glroot.PATH      = "$PAM_L0";
+                       };
+                     };
+                     //              glroot.PATH        = path;
+                     glroot.NAME        = name;
+                     //
+                     insertPamelaRootFile( &glroot );
+                     //
+                     id_root = glroot.ID;
+                 }else{
+                     row = result->Next();
+                     if(row)id_root = (UInt_t)atoll(row->GetField(0));
+                 }
+                 if ( IsDebug() ) cout << "ID_ROOT = "<<id_root<<endl;
+                 // ---------------------------------------------------------------------
+                 // step 3: modify time-tag of corrupted GL_TRK_CALIB entry
+                 // ---------------------------------------------------------------------
+                 if ( IsDebug() ) cout << " << Modify time-tag of calibration ID="<<id<<" >> "<<endl;
+                 oss.str("");
+                 oss << " UPDATE GL_TRK_CALIB SET "
+                     << " TO_TIME=0 , FROM_TIME=0  WHERE "
+                     << " ID = "<< id  << ";";
+                 if ( IsDebug() ) cout << oss.str().c_str() << endl;
+                 result = conn->Query(oss.str().c_str());
+                 if ( !result ) throw -4;;
+                 // ---------------------------------------------------------------------
+                 // step 4: insert the new calibration:
+                 // ---------------------------------------------------------------------
+                 if ( IsDebug() ) cout << " << Insert new TRK calibration >> "<<endl;
+                 //
+                 GL_TRK_CALIB glcal = GL_TRK_CALIB();
+                 //
+                 glcal.ID_ROOT_L0        = id_root;
+                 glcal.EV_ROOT_CALIBTRK1 = this_t1;
+                 glcal.EV_ROOT_CALIBTRK2 = this_t2;
+                 glcal.FROM_TIME         = fromtime;
+                 glcal.TO_TIME           = totime;
+                 glcal.OBT1              = obt1;
+                 glcal.OBT2              = obt2;
+                 glcal.PKT1              = pkt1;
+                 glcal.PKT2              = pkt1;
+                 glcal.BOOT_NUMBER       = GetBOOTnumber();
+                 glcal.VALIDATION        = this_valid;
+                 //
+                 HandleTRK_CALIB(&glcal);
+                 if ( IsDebug() ) cout << "ID = "<<glcal.ID<<endl;
+                 //
+             }
+             if ( IsDebug() ) cout << "------------------------------------------------------------" <<endl;
+         }while(1);//end loop over calibrations
+         if( result )delete result;
+         if( row    )delete row;
+     }
+ //      // ------------------------------
+ //      // try to recover the calibration
+ //      // ------------------------------
+ //      cout << "TRY TO RECOVER TRACKER CALIBRATION"<<endl;
+ //      //
+ //      ULong64_t time  = 0; //absolute time
+ //      string    path[100]; //mettere un limite massimo
+ //      int       nrows = 0;
+ //      UInt_t    pkt   = 0;
+ //      UInt_t    obt   = 0;
+ //      char     *type  = "";
+ //      EventHeader *eh   = new EventHeader();
+ //      CalibTrk1Event *c = new CalibTrk1Event();
+ //       //
+ //       if(which_is_not_valid==1 || which_is_not_valid==3){
+ //        //
+ //        cout << "PKT1 --> missing or corrupted "<<endl;
+ //        type = "CalibTrk1";
+ //        pkt  = pkt1;
+ //        obt  = obt1;
+ //        time = this->GetAbsTime(obt1);
+ //        if( pkt1 == 0 ){//missing
+ //            time = this->GetAbsTime(obt2);
+ //            pkt  = pkt2-1;
+ //        }
+ //        //
+ //       }else if (which_is_not_valid==2 || which_is_not_valid==3){
+ //        //
+ //        cout << "PKT2--> missing or corrupted "<<endl;
+ //        type = "CalibTrk2 ";
+ //        pkt  = pkt2;
+ //        obt  = obt2;
+ //        time = this->GetAbsTime(obt2);
+ //        if( pkt2 == 0 ){//missing
+ //            time = this->GetAbsTime(obt1);
+ //            pkt  = pkt1+1;
+ //        }
+ //        //
+ //       }else{
+ //        cout << "this should not happen!!! "<<endl;
+ //        trow -666;
+ //       }
+ //       nrows = Query_ROOT_TABLE(time,conn,path);// get the list of file which might contain the packet
+ //       for(int r=0; r<nrows; r++){ //loop over rows
+ //        if(path)cout << r << " >>>> "<<(path+r)->c_str() << endl;
+ //        /// verifica che il file non sia quello gia` aperto
+ //       }
+ //       ////////////////////////////////////////////////////////////////////////
+ //       TSQLResult *result = 0;
+ //       TSQLRow    *row    = 0;
+ //       //
+ //       stringstream oss;
+ //       oss.str("");
+ //       // ----------------------------------------
+ //       // read the id of last calibration inserted
+ //       // ----------------------------------------
+ //       oss.str("");
+ //       oss << " SELECT ";
+ //       oss << " (ID,ID_ROOT_L0,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,FROM_TIME,TO_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION) ";
+ //       oss << " ORDER BY ID DESC LIMIT 1; ";
+ //       result = conn->Query(oss.str().c_str());
+ //       row = result->Next();
+ //       if( !row )throw -666;
+ //       if( result )delete result;
+ //       if( row    )delete row;
+ //       UInt_t id = (UInt_t)atoll(row->GetField(0));
+ //       // -------------------------------------
+ //       // ...and modify it with new parameters
+ //       // -------------------------------------
+ //   }
    //
    return(0);
  };
-Line 5653 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 6256 
 Bool_t PamelaDBOperations::MissingTRK_CA
    GL_TRK_CALIB* trkcalib = new GL_TRK_CALIB();
-   // get the closest VALIDATED calibration before the run start (t2)
+   // get the closest calibration before the run start (t2)
    if ( trkcalib->Query_GL_TRK_CALIB(t2, conn) )return(true);    //>>> missing
    if ( trkcalib->TO_TIME  < t2 ) return(true);                                  //>>> missing
-Line 6347 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 6950 
 void PamelaDBOperations::RemoveCALIBS(){
   * Rearrange calibration tables
   *
   **/
- UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){
+ UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh , TFile *file){
-   Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0};
+     if(!caltrk) return 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 ( IsDebug() ) cout << "ValidateTrkCalib:"<<endl;
-   //  -----------------------------------------------
+     UInt_t validate = 1;
-   //  If missing packets: check the acq configuration
+     Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0};
-   //  (some DSPs might be excluded from acquisition)
+     UInt_t timeaftercalib=120000; //2000;
-   //  -----------------------------------------------
+     TString classname = caltrk->GetName();
-   //  -----------------------------------------------
+ //  ----------------------------------
-   //  retrieve the first run header after calib
+ //  Check CRCs and failed calibrations
-   //  -----------------------------------------------
+ //  ----------------------------------
-   PacketType *pctp;
+     for(Int_t ipkt=0; ipkt<6; ipkt++){
-   EventCounter *cod;
+         if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){
-   cod = eh->GetCounter();
+             if( caltrk->crc_hcal[ipkt] ){
-   Int_t irun = cod->Get(pctp->RunHeader);
+ //              if(IsDebug())cout<<"(CRC Header)";
-   TTree *rh=(TTree*)file->Get("RunHeader");
+                 validate = 0;
-   if ( !rh || rh->IsZombie() ) throw -17;
+                 if(IsDebug())cout <<endl<<" *** CRC *** (header DSPn "<<caltrk->DSPnumber[ipkt]<<")";
-   if( rh->GetEntries() == irun ){
-     if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1)  -- cannot validate :-( "<<endl;
+             }
-     return 0; // :-(
+             for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] ){
-   }
+ //              if(IsDebug())cout<<"(CRC Pkt-"<<ilad<<")";
+                 if(IsDebug())cout <<endl<<" *** CRC *** (data DSPn "<<caltrk->DSPnumber[ipkt]<<" ladder "<<ilad<<")";
+                 validate = 0;
+             }
+             if( !(caltrk->ncalib_event[ipkt]==0 && caltrk->cal_flag[ipkt]==0) ){
+                 if(IsDebug())cout <<endl<<" *** FAILURE ***     (data DSPn "<<caltrk->DSPnumber[ipkt]<<")";
+                 validate = 0;
+             }
+         }else{
+ //          validate=0;
+             if(IsDebug())cout <<endl<<" *** DSPn *** ("<< caltrk->DSPnumber[ipkt] <<" @pkt "<<ipkt<<")";
+         }
+     }
-   RunHeaderEvent *run  = 0;
+ //  -----------------------
-   EventHeader    *hrun = 0;
+ //  Check missing packets:
-   rh->SetBranchAddress("RunHeader", &run);
+ //  -----------------------
-   rh->SetBranchAddress("Header", &hrun);
+ //    Readout order:
-   rh->GetEntry(irun);
+ //    ------------------
-   //     cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl;
+ //    DSP   packet board
+ //    ------------------
-   if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){
+ //    12    0      1
-     if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl;
+ //    10    1      1
-     return 0; // :-(
+ //     8    2      1
-   }
+ //     4    3      1
+ //     6    4      1
-   if( !run->RM_ACQ_AFTER_CALIB ){
+ //     2    5      1
-     if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0    -- cannot validate :-( "<<endl;
+ //    ------------------
-     return 0; // :-(
+ //    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;
+     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 valid packets and set build variable
+ //  ----------------------------------------------------
+     if(IsDebug())cout <<endl<< " DSP: ";
+     Int_t  npkts=0;
+     for(Int_t ipkt=0; ipkt<6; ipkt++){
+         if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){
+             if(IsDebug())cout <<" "<<caltrk->DSPnumber[ipkt];
+             npkts++;
+             build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) );
+ //          cout << caltrk->DSPnumber[ipkt]
+         };
+     }
+     if(IsDebug())cout << " ==> "<< hex << build << dec;
+ //  ----------------------------------------------------
+ //  If the number of valid packets is 6, ok exit...
+ //  ----------------------------------------------------
+     if( npkts==6 ){
+         return validate; // exit
+     }
+ ////////////////////////////////////////////////////////
+ //  ...otherwise there might be some missing packets
+ //
+ //  In this case check the acq configuration
+ //  (some DSPs might be excluded from acquisition)
+ ////////////////////////////////////////////////////////
+     if(!eh || !file || (file&&file->IsZombie()) ){
+         if ( IsDebug() )cout << " *** MISSING VIEW *** eh="<<eh<<" file="<<file<<" cannot validate"<<endl;
+         return (0);
+     }
+ //  -----------------------------------------------
+ //  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 << " *** MISSING-PKT *** 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);
+     if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){
+         if ( IsDebug() ) cout << " *** MISSING-PKT *** no runs after calib (2) -- cannot validate :-( "<<endl;
+         return 0; // :-(
+     }
-   UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime());
+     UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime());
-   if( dtime > timeaftercalib ){
+     if( dtime > timeaftercalib ){
-     if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl;
+         if ( IsDebug() ) cout << " *** MISSING-PKT *** run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl;
-     return 0; // :-(
+         return 0; // :-(
-   }
+     }
+     if ( IsDebug() ) cout <<endl<< " ACQ_BUILD_INFO ==> "<<hex<<(run->ACQ_BUILD_INFO & mask)<<dec;
+     if( (run->ACQ_BUILD_INFO & mask) != build ){
+         validate=0; // :-(
+         cout <<endl<< " *** MISSING-PKT *** packet mismatch: ACQ_BUILD_INFO="<<hex<<(run->ACQ_BUILD_INFO&mask)<<" != "<<build<<dec;
+     };
+     return validate;
-   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; // :-)
  }

 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.7
 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.7
-Removed from v.1.3
+Added in v.1.7

	ViewVC Help
Powered by ViewVC 1.1.23