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

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

Parent Directory | Revision Log | View Patch Patch

-revision 1.2 by mocchiut,
Thu Sep 25 12:13:54 2008 UTC
+revision 1.20 by mocchiut,
Tue Dec 23 20:43:10 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 99 
 PamelaDBOperations::PamelaDBOperations(T
+Line 101 
 PamelaDBOperations::PamelaDBOperations(T
    if (INSERT_RAW) SetRawName(filerawname);
    //
    INSERT_ROOT = !filerootname.IsNull();
+   if ( INSERT_ROOT ) this->SetRootName(filerootname);
    this->SetOrbitNo(dwinput);
    //
    this->SetID_RAW(0);
-Line 132 
 void PamelaDBOperations::OpenL0File(TStr
+Line 135 
 void PamelaDBOperations::OpenL0File(TStr
        stringstream qu;
        TSQLResult *result = 0;
        TSQLResult *result2 = 0;
+       TSQLResult *result3 = 0;
        TSQLRow    *row    = 0;
        TSQLRow    *row2    = 0;
+       TSQLRow    *row3    = 0;
        TString chpath;
        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 154 
 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 157 
 void PamelaDBOperations::OpenL0File(TStr
+Line 163 
 void PamelaDBOperations::OpenL0File(TStr
        };
        delete result;
        qu.str("");
-       qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL from ROOT_TABLE where ID_N=" << ridn << ";";
+       qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL,REAL_TIME_INIT from ROOT_TABLE where ID_N=" << ridn << ";";
        if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
        result2 = conn->Query(qu.str().c_str());
        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));
            chresursts = (UInt_t)atoll(row2->GetField(4));
            chboot = (UInt_t)atoll(row2->GetField(5));
            //
+           if ( !chboot ){
+             qu.str("");
+             qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<10000 AND (REAL_TIME_INIT-"<< row2->GetField(10) <<")<7200 group by BOOT_NUMBER order by BOOT_NUMBER asc;";
+             if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
+             result3 = conn->Query(qu.str().c_str());
+             if ( result3 ){
+               row3 = result3->Next();
+               if ( row3 && result3->GetRowCount() == 1 ){
+                 chboot = (UInt_t)atoll(row3->GetField(0));
+                 if ( debug ) printf(" Found boot_number = %u \n",chboot);
+               } else {
+                 if ( debug ) printf(" AGH CANNOT DETERMINE THE BOOT NUMBER... \n");
+                 throw -29;
+               };
+             };
+           };
+           //
            chpktinit = (UInt_t)atoll(row2->GetField(6));
            chpktfinal = (UInt_t)atoll(row2->GetField(7));
            chobtinit = (UInt_t)atoll(row2->GetField(8));
-Line 182 
 void PamelaDBOperations::OpenL0File(TStr
+Line 209 
 void PamelaDBOperations::OpenL0File(TStr
          throw -85;
        };
        filerootname = chpath + chfile;// + ".root";
-       if ( debug ) printf(" chewbacca: filename is %s \n",filerootname.Data());
      };
      this->SetRootName(filerootname);
-     file = TFile::Open(this->GetRootName().Data());
-   } else {
-     this->SetRootName("");
    };
+   if ( debug ) printf(" Filename is %s \n",filerootname.Data());
+   file = TFile::Open(this->GetRootName().Data());
+   //  } else {
+   //    this->SetRootName("");
+   //    this->SetRootName(filerootname);
+   //  };
  }
  //
-Line 388 
 void PamelaDBOperations::SetOrbitNo(UInt
+Line 417 
 void PamelaDBOperations::SetOrbitNo(UInt
      return;
    };
    //
-   TString name = this->GetRootFile();
+   if ( !chewbacca ){
-   Int_t nlength = name.Length();
+     TString name = this->GetRootFile();
-   if ( nlength < 5 ) return;
+     Int_t nlength = name.Length();
-   TString dwo = 0;
+     if ( nlength < 5 ){
-   for (Int_t i = 0; i<5; i++){
+       if ( IsDebug() ) printf(" Agh problems determining the orbit number! name = %s \n",name.Data());
-     dwo.Append(name[i],1);
+       return;
-   };
+     };
-   if ( dwo.IsDigit() ){
+     TString dwo = 0;
-     dworbit = (UInt_t)dwo.Atoi();
+     for (Int_t i = 0; i<5; i++){
-   } else {
-     dwo="";
-     for (Int_t i = 8; i<13; i++){
        dwo.Append(name[i],1);
      };
-     if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
+     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());
    };
-   if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());
    return;
  };
-Line 435 
 TString PamelaDBOperations::GetRawPath()
+Line 469 
 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 649 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 687 
 Int_t PamelaDBOperations::SetUpperLimits
      UInt_t tjump = 50000;
      //UInt_t tjump = 100;
      while ( tjump > 0 ){
-       pktlast = numeric_limits<UInt_t>::max();
+       //      pktlast = numeric_limits<UInt_t>::max();
-       while ( pktlast > chpktmax && (Int_t)(nevent-1-it) >= 0 ){
+       pktlast = chpktmax + 1;
+       while ( PKT(pktlast) > PKT(chpktmax) && (Int_t)(nevent-1-it) >= 0 ){
          if ( (Int_t)(nevent-1-it) >= 0 ){
            T->GetEntry(nevent-1-it);
            ph = eh->GetPscuHeader();
-Line 658 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 697 
 Int_t PamelaDBOperations::SetUpperLimits
          } else {
            pktlast = chpktmax + 1;
          };
-         if ( (!(it%1000) || abs((int)pktlast - (int)chpktmax)<1000 ) && debug ) printf(" look for up %i %i %i nevent %u (nevent-1-it) %i \n",it,pktlast,chpktmax,nevent,(Int_t)(nevent-1-it));
+         if ( (!(it%1000) || abs((int)(PKT(pktlast) - PKT(chpktmax)))<1000 ) && debug ) printf(" look for up %i %i %i nevent %u (nevent-1-it) %i \n",it,pktlast,chpktmax,nevent,(Int_t)(nevent-1-it));
          it += tjump;
        };
        if ( tjump > 1 ) it -= 2*tjump;
-Line 682 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 721 
 Int_t PamelaDBOperations::SetUpperLimits
      tjump = 50000;
      //tjump = 100;
      while ( tjump > 0 ){
-       pktlast = 0;
+       //      pktlast = 0;
-       while ( pktlast < chpktmin && it < (Int_t)nevent ){
+       pktlast = chpktmin - 1;
+       if ( debug ) printf("LLlook for down %i %i %llu %llu \n",it,pktlast,PKT(pktlast),PKT(chpktmin));
+       while ( PKT(pktlast) < PKT(chpktmin) && it < (Int_t)nevent ){
          if ( it < (Int_t)nevent ){
            T->GetEntry(it);
            ph = eh->GetPscuHeader();
-Line 1125 
 UInt_t PamelaDBOperations::AssignRunID()
+Line 1166 
 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 1322 
 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 1330 
 Long64_t PamelaDBOperations::OBT(UInt_t
+Line 1379 
 Long64_t PamelaDBOperations::OBT(UInt_t
    //
    //  if ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u  (numeric_limits<UInt_t>::max()/2) is %u  \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2));
    //
-   if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){
+   if ( obt < ((Long64_t)pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){
      //    if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max());
      return((Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max());
    };
-Line 1377 
 void PamelaDBOperations::FillClass(Bool_
+Line 1426 
 void PamelaDBOperations::FillClass(Bool_
    UInt_t lastPkt = 0;
    UInt_t rhtime = 0;
    UInt_t rttime = 0;
+   //
+   if ( IsDebug() ) printf(" A firstev %i lastev %i nevents %i \n",firstev,lastev,lastev-firstev+1);
+   //
    if ( !mishead ){
      codh = ehh->GetCounter();
-     if ( lastev+1 == firstev || (lastev == firstev && lastev != -1) ){
+     if ( (lastev+1 == firstev && lastev != -1) || (lastev == firstev && lastev != -1) ){
+       if ( IsDebug() ) printf(" B firstev %i lastev %i nevents %i \n",firstev,lastev,lastev-firstev+1);
        firstev = 1;
        lastev = 0;
      } else {
-Line 1392 
 void PamelaDBOperations::FillClass(Bool_
+Line 1445 
 void PamelaDBOperations::FillClass(Bool_
    };
    if ( !mistrail ){
      codt = eht->GetCounter();
-     if ( lastev+1 == firstev || (lastev == firstev && lastev != -1)){
+     if ( (lastev+1 == firstev && lastev != -1) || (lastev == firstev && lastev != -1)){
+       if ( IsDebug() ) printf(" C firstev %i lastev %i nevents %i \n",firstev,lastev,lastev-firstev+1);
        lastev = 0;
        firstev = lastev+1;
      } else {
-Line 1409 
 void PamelaDBOperations::FillClass(Bool_
+Line 1463 
 void PamelaDBOperations::FillClass(Bool_
    if ( mishead ) {
      glrun->Set_GL_RUNH0();
      //
-     if ( lastev+1 == firstev || (lastev == firstev && lastev != -1) ){
+     if ( (lastev+1 == firstev && lastev != -1) || (lastev == firstev && lastev != -1) ){
        firstObt = lastObt;
        firstPkt = lastPkt;
        rhtime = rttime;
-Line 1428 
 void PamelaDBOperations::FillClass(Bool_
+Line 1482 
 void PamelaDBOperations::FillClass(Bool_
    if ( mistrail ){
      glrun->Set_GL_RUNT0();
      //
-     if ( lastev+1 == firstev || (lastev == firstev && lastev != -1) ){
+     if ( (lastev+1 == firstev && lastev != -1) || (lastev == firstev && lastev != -1) ){
        lastObt = firstObt;
        lastPkt = firstPkt;
        rttime = rhtime;
-Line 1463 
 void PamelaDBOperations::FillClass(Bool_
+Line 1517 
 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 1514 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1604 
 Int_t PamelaDBOperations::insertPamelaGL
      result = conn->Query(oss.str().c_str());
      if ( !result ) throw -10;
      row = result->Next();
+     if ( !row ) throw -92;
      idresof = (UInt_t)atoll(row->GetField(0));
      existsts = true;
      goto eout;
-Line 1535 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1626 
 Int_t PamelaDBOperations::insertPamelaGL
      if ( !row ){
        oss.str("");
        oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE),ID FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "
-           << dworbit << " order by FROM_ORBIT desc limit 1;";
+           << dworbit << " AND SPECIAL_FILE='' 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;
-Line 1793 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1884 
 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 2138 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2267 
 Int_t PamelaDBOperations::insertPamelaRU
          ptht = cod->Get(pctp->RunHeader) - 1;
          //      evbeft = cod->Get(pctp->Physics);
          tcod = (UInt_t)cod->Get(pctp->Physics);
-         if ( !tcod ) tcod = 1;
+         if ( !tcod ) tcod = 1;    // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
          evbeft = TMath::Min(upperentry,(tcod-1));
          if ( debug ) printf(" Loop in runtrailers, evbeft is %u upperentry %u cod->getetc %u \n",evbeft,upperentry,cod->Get(pctp->Physics));
          //
-Line 2177 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2306 
 Int_t PamelaDBOperations::insertPamelaRU
            obth = phh->GetOrbitalTime();
            cod = ehh->GetCounter();
            tcod = (UInt_t)cod->Get(pctp->Physics);
-           if ( !tcod ) tcod = 1;
+           if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
            evbefh = TMath::Max(chminentry,(tcod-1));
            //
            if ( PKT(pkth) >= PKT(pktfirst) && PKT(pkth) <= upperpkt ){
-Line 2216 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2345 
 Int_t PamelaDBOperations::insertPamelaRU
                //
                this->HandleRunFragments(true,true,chminentry,(evbefh));
                //
+               // what about pth=ptht-1
              };
            };
            //
-Line 2238 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2368 
 Int_t PamelaDBOperations::insertPamelaRU
                cod = ehh->GetCounter();
                //            evbeft = cod->Get(pctp->Physics);
                tcod = (UInt_t)cod->Get(pctp->Physics);
-               if ( !tcod ) tcod = 1;
+               if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
                evbeft = TMath::Min(upperentry,(tcod-1));
                rh->GetEntry(pth);
                phh = ehh->GetPscuHeader();
-Line 2278 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2408 
 Int_t PamelaDBOperations::insertPamelaRU
          //      if ( ptt+1 == rtev){
          if ( conptt+1 == nrtev ){
            //    if ( conptt+1 == (nrtev+nrtbef )){
+           if ( IsDebug() ) printf(" >>>>>>>>>>> %i %u %i %u \n",ptht,rhev,nrtev,conptt);
            ptht++;
            if ( ptht < rhev ){
              rh->GetEntry(ptht);
+             //pth++;
+             //if ( pth < rhev ){
+             //rh->GetEntry(pth);
              phh = ehh->GetPscuHeader();
              pkth = phh->GetCounter();
              obth = phh->GetOrbitalTime();
-Line 2304 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 2438 
 Int_t PamelaDBOperations::insertPamelaRU
              if ( (UInt_t)evbeft < upperentry-1 && upperentry>0 ){
                if ( IsDebug() ) printf(" Piece of run at the end of the file with NO RUNHEADER! evbeft %u upperentry-1 %u \n",(UInt_t)evbeft,upperentry-1);
                //
-               if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
+               //              if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
-               rt->GetEntry(ptt-1);
+               if ( (rtev-1) < 0 || ptt < 0 ) throw -15; // should never arrive here!
+               //              rt->GetEntry(ptt-1);
+               rt->GetEntry(rtev-1);
                cod = eht->GetCounter();
-               tcod = (UInt_t)cod->Get(pctp->Physics);
+               tcod = (UInt_t)cod->Get(pctp->Physics)+1;
                evbefh = TMath::Max(chminentry,tcod);
                //            evbefh = cod->Get(pctp->Physics);
                rt->GetEntry(ptt);
-Line 2538 
 void PamelaDBOperations::HandleRunFragme
+Line 2674 
 void PamelaDBOperations::HandleRunFragme
      //
      // we have now the good first piece of a run, fill the glrun object
      //
-     if ( rhfirstev != firstev && !mishead ) mishead = true;
+     if ( firstev != (lastev+1) ){ // could be a problem sometimes (?)
-     if ( rtlastev != lastev && !mistrail ) mistrail = true;
+       if ( rhfirstev != firstev && !mishead ) mishead = true;
+       if ( rtlastev != lastev && !mistrail ) mistrail = true;
+     };
      //
+     if ( IsDebug() ) printf(" bhere firstev is %i lastev is %i \n",firstev,lastev);
      this->FillClass(mishead,mistrail,firstev,lastev);
+     if ( IsDebug() ) printf(" chere firstev is %i lastev is %i \n",firstev,lastev);
      //
      if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n");
      if ( IsDebug() ) printf(" C THIS RUN: RUNHEADER_OBT %u RUNTRAILER_OBT %u RUNHEADER_PKT %u RUNTRAILER_PKT %u \n", glrun->GetRUNHEADER_OBT(),glrun->GetRUNTRAILER_OBT(),glrun->GetRUNHEADER_PKT(),glrun->GetRUNTRAILER_PKT());
-Line 2636 
 void PamelaDBOperations::HandleRunFragme
+Line 2776 
 void PamelaDBOperations::HandleRunFragme
          //
          found = false; // default value
          //
-         if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");
+         if ( IsDebug() ) printf(" A Found a possible candidate, checking if it is the good one... \n");
          //
          // if we have both runheader and runtrailer we can check with pkt_counter:
          //
-Line 2652 
 void PamelaDBOperations::HandleRunFragme
+Line 2792 
 void PamelaDBOperations::HandleRunFragme
              if ( IsDebug() ) printf(" FOUND!!! check %llu pktt %llu \n",chkpkt,pktt);
              //
              found = true;
+             if ( IsDebug() ) printf(" where firstev is %i lastev is %i \n",firstev,lastev);
              //
            } else {
              //
-Line 2717 
 void PamelaDBOperations::HandleRunFragme
+Line 2858 
 void PamelaDBOperations::HandleRunFragme
          oss << " ID="<<row->GetField(0)<<";";
          //
          glrun1->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); // here we have runheader infos
+         if ( IsDebug() ) printf(" there firstev is %i lastev is %i \n",firstev,lastev);
          //
          // merge infos
          //
-Line 2747 
 void PamelaDBOperations::HandleRunFragme
+Line 2889 
 void PamelaDBOperations::HandleRunFragme
          glrun1->SetRUNTRAILER_OBT(glrun->GetRUNTRAILER_OBT());
          glrun1->SetRUNTRAILER_PKT(glrun->GetRUNTRAILER_PKT());
          //
+         if ( IsDebug() ) printf(" here firstev is %i lastev is %i \n",firstev,lastev);
+         //
          glrun->SetEV_FROM(firstev);
          glrun->SetNEVENTS(lastev-firstev+1);
          //
-Line 2839 
 void PamelaDBOperations::HandleRunFragme
+Line 2983 
 void PamelaDBOperations::HandleRunFragme
          //
          found = false; // default value
          //
-         if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");
+         if ( IsDebug() ) printf(" B Found a possible candidate, checking if it is the good one... \n");
          //
          // if we have both runheader and runtrailer we can check with pkt_counter:
          //
-Line 3068 
 void PamelaDBOperations::HandleMissingHo
+Line 3212 
 void PamelaDBOperations::HandleMissingHo
        glrun->SetID(this->AssignRunID());
        glrun->SetID_RUN_FRAG(0);
        glrun->Fill_GL_RUN(conn);    // it'ok we arrive here only inside a file hence in the middle of the runs...
-     };
+       };
      //
    };
    //
-Line 3112 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3256 
 Bool_t PamelaDBOperations::IsRunConsiste
    if ( firstev == lastev+1 || lastev == firstev ) { // no events inside the run!
      //if ( firstev <= lastev+1 ) { // no events inside the run!
      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
+     firstev = lastev+1; // this is necessary for files with no Physics entries, should have no influence on other files
      // return true is correct
      return(true);
      //
-Line 3141 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3286 
 Bool_t PamelaDBOperations::IsRunConsiste
        //
      } else {
        //
-       if ( IsDebug() ) printf(" There are no-phyics packets inside the run!\n");
+       if ( IsDebug() ) printf(" There are no-physics packets inside the run!\n");
        //
        // HERE WE MUST HANDLE THAT RUNS AND GO BACK
        //
-Line 3151 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3296 
 Bool_t PamelaDBOperations::IsRunConsiste
        UInt_t check = 0;
        UInt_t lastevtemp = lastev;
        UInt_t firstevno = firstev;
+       UInt_t rhchko=0;
+       UInt_t rhchk=0;
        //
        for (UInt_t i=firstev; i<=lastev; i++){
          //
-Line 3159 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3306 
 Bool_t PamelaDBOperations::IsRunConsiste
          //
          check = 0;
          //
+         // if we have a runheader set lastev then exit
+         //
          for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){
            if ( strcmp(*Iter,"Physics") ) check += code->Get(GetPacketType(*Iter));
          };
+         // check here if we have a runheader
+         rhchko = rhchk;
+         rhchk = code->Get(GetPacketType("RunHeader"));
          //
          if ( checkfirst < check || i == lastev ){
            //
-Line 3245 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3398 
 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 3371 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3512 
 Bool_t PamelaDBOperations::IsRunConsiste
                  //
                  found = false; // default value
                  //
-                 if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");
+                 if ( IsDebug() ) printf(" C Found a possible candidate, checking if it is the good one... \n");
                  //
                  // if we have both runheader and runtrailer we can check with pkt_counter:
                  //
-Line 3576 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3717 
 Bool_t PamelaDBOperations::IsRunConsiste
              }; // EEE
            } else {
              if ( !IsRunAlreadyInserted() ){
                glrun->SetID(this->AssignRunID());
-Line 3621 
 Bool_t PamelaDBOperations::IsRunConsiste
+Line 3731 
 Bool_t PamelaDBOperations::IsRunConsiste
            //
          };
          //
+         if ( i > firstev ){
+           if ( rhchko != rhchk ){
+             if ( IsDebug() ) printf("oh oh... we have a runheader! stop here and handle later the remaining piece\n");
+             lastev = i;
+             return(false);
+           };
+         };
+         //
          if ( check == checklast && i != lastev ){
            lastevtemp = i - 1;
            i = lastev - 1;
-Line 4441 
 Int_t PamelaDBOperations::insertCALOPULS
+Line 4559 
 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");
+     glcal->ID         = 0;
-     if ( PEDANTIC ) throw -80;
+     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;
      //
-   } else {
+     HandleTRK_CALIB(glcal);
+     //
+     delete glcal;
+ }
+ /**
+  * 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;
+     //
+     TSQLResult *result = 0;
+     TSQLRow    *row    = 0;
+     //
+     stringstream oss;
+     oss.str("");
      //
-     // we have to insert a new calibration, check where to place it
+     //
+     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,TO_TIME FROM GL_TRK_CALIB WHERE "
+     oss << " SELECT ID FROM GL_TRK_CALIB WHERE "
-         << " FROM_TIME < "<< fromtime << " AND "
+         << " BOOT_NUMBER = "<< boot_number; //
-         << " TO_TIME > "<< fromtime << ";";
+     oss << " AND FROM_TIME="<<fromtime; /// NEWNEWNEW -- VA BENE ?!?!?!?!
+     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 where to place the 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();
      //
-     if ( !row ){
+     if ( row ){
-       //
+         //
-       // no calibrations in the db contain our calibration
+         if ( IsDebug() ) printf(" Trk calibration already inserted in the DB\n");
-       //
+         if ( PEDANTIC ) throw -80;
-       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
+         // we have to insert a new calibration, check where to place it
-       //
+         //
-       totime = (UInt_t)atoll(row->GetField(1));
+         oss.str("");
-       //
+         oss << " SELECT ID,TO_TIME FROM GL_TRK_CALIB WHERE "
-       oss.str("");
+             << " FROM_TIME < "<< fromtime << " AND "
-       oss << " UPDATE GL_TRK_CALIB SET "
+             << " TO_TIME > "<< fromtime << ";";
-           << " 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(" Check where to place the trk calibration: query is \n %s \n",oss.str().c_str());
-       //
+         result = conn->Query(oss.str().c_str());
-       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;
-       //
+         //
-       if ( !result ) throw -4;
+         row = result->Next();
-       //
+         //
+         if ( !row ){
+             //
+             // 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 < 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("");
+         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 << "NULL,";
+         } else {
+             oss << "'"
+                 << t1 << "',";
+         };
+         //
+         if ( !pk2 ){
+             oss << "NULL,'";
+         } else {
+             oss << "'"
+                 << t2 << "','";
+         };
+         //
+         oss << fromtime << "','"
+             << totime << "','"
+             << obt1 << "','"
+             << pkt1 << "','"
+             << obt2 << "','"
+             << pkt2 << "','"
+             << boot_number << "','"
+             << valid << "');";
+         //
+         if ( IsDebug() ) printf(" Insert 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("");
-     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) "
+     oss << " SELECT ID FROM GL_TRK_CALIB ORDER BY ID DESC LIMIT 1 ;";
-         << " VALUES (NULL,' "
+     if ( IsDebug() ) cout << oss.str().c_str() << endl;
-         << idroot << "',";
-     //
-     if ( !pk1 ){
-       oss << "NULL,";
-     } else {
-       oss << "'"
-           << t1 << "',";
-     };
-     //
-     if ( !pk2 ){
-       oss << "NULL,'";
-     } else {
-       oss << "'"
-           << t2 << "','";
-     };
-     //
-     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());
-     //
      result = conn->Query(oss.str().c_str());
+     if ( !result ) throw -4;;
+     row = result->Next();
+     if(row)glcal->ID = (UInt_t)atoll(row->GetField(0));
      //
-     if ( !result ) throw -4;
-     //
-   };
-   //
  };
  /**
   * 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;
+     CalibTrk1Event *caltrk1 = 0;
-     tr1->GetEntry(t1);
+     CalibTrk2Event *caltrk2 = 0;
+     TTree *tr1 = 0;
+     TTree *tr2 = 0;
+     EventHeader *eh1 = 0;
+     PscuHeader *ph1 = 0;
+     EventHeader *eh2 = 0;
+     PscuHeader *ph2 = 0;
      //
-     ph1 = eh1->GetPscuHeader();
+     PacketType *pctp=0;
-     obt1 = ph1->GetOrbitalTime();
+     EventCounter *codt2=0;
-     pkt1 = ph1->GetCounter();
-     fromtime = this->GetAbsTime(ph1->GetOrbitalTime());
      //
-     //     valid = 1;
+     Int_t nevents1 = 0;
-     //     //
+     Int_t nevents2 = 0;
-     //     if ( caltrk1->unpackError != 0 && caltrk1->good0 == 0 ) valid = 0;// CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT1
      //
+     fromtime = 0;
      //
-     if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->PKT(pkt1) <= upperpkt && this->OBT(obt1) >= this->OBT(obtfirst) && this->OBT(obt1) <= upperobt ){
+     obt1 = 0;
-       //    if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
+     pkt1 = 0;
-       //
+     obt2 = 0;
-       if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
+     pkt2 = 0;
-       //
+     //
-       valid = ValidateTrkCalib( caltrk1, eh1 );
+     tr1 = (TTree*)file->Get("CalibTrk1");
-       if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
+     if ( !tr1 || tr1->IsZombie() ) throw -22;
-       //
+     tr2 = (TTree*)file->Get("CalibTrk2");
-       // Do we have the second calibration packet?
+     if ( !tr2 || tr2->IsZombie() ) throw -23;
-       //
+     //
-       while ( t2t1cal < t1+1 ){ // get the calibration packet2 that follows the packet1
+     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;
+     //
+     bool MISSING_pkt1 = true;
+     bool MISSING_pkt2 = true;
+     int  ncalib = 0;
+     bool try_to_recover = false;
+     //
+     for (t1=0; t1 < nevents1; t1++){//loop over packet1
          //
-         t2++;
+         pret2 = t2;
+         tr1->GetEntry(t1);
          //
-         pret2 = t2 - 1; // EMILIANO
+         ph1 = eh1->GetPscuHeader();
+         obt1 = ph1->GetOrbitalTime();
+         pkt1 = ph1->GetCounter();
+         fromtime = GetAbsTime(ph1->GetOrbitalTime());
          //
-         if ( t2 < nevents2 ){
+         // chek if the packet number and obt are consistent with the other packets ???
-           tr2->GetEntry(t2);
+         //
-           codt2 = eh2->GetCounter();
+         if ( PKT(pkt1) >= PKT(pktfirst) && PKT(pkt1) <= upperpkt && OBT(obt1) >= OBT(obtfirst) && OBT(obt1) <= upperobt ){
-           t2t1cal = codt2->Get(pctp->CalibTrk1);
+             //    if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
-           //
+             //
-           ph2 = eh2->GetPscuHeader();
+             if ( IsDebug() ) printf("\n Trk calibration1 %u at time %u obt %u pkt %u \n",t1,fromtime,obt1,pkt1);
-           obt2 = ph2->GetOrbitalTime();
+             //
-           pkt2 = ph2->GetCounter();
+             valid = ValidateTrkCalib( caltrk1, eh1 );
-           //
+             if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
-           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
+             // 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 {
            //
-           // running out of vector without finding the corresponding calibration, sig
+             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;
-           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
+     }; //end loop on pkt1
-         //      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;
-         //      };
+     //
+     // we have one more calib pkt2 !
+     //
+     t2++;
+     while ( t2 < nevents2 ){
          //
-       };
+         // handle missing calib1
-       //
-       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);
-         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 );
-         if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
-         valid = valid & valid2;
-         //
-         // Handle good calib
          //
-         this->HandleTRK_CALIB(true,true);
+         if ( IsDebug() ) printf(" t2 is %u nevents2 is %u \n",t2,nevents2);
+         obt1 = 0;
+         pkt1 = 0;
          //
-         // Check for missing calibtrk1
+         tr2->GetEntry(t2);
+         ph2 = eh2->GetPscuHeader();
+         obt2 = ph2->GetOrbitalTime();
+         pkt2 = ph2->GetCounter();
          //
-         if ( t2 != pret2+1 ){
+         fromtime = this->GetAbsTime(ph2->GetOrbitalTime());
-           //
+         valid = 0;
-           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);
+         //  if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->PKT(pkt1) <= upperpkt && this->OBT(obt1) >= this->OBT(obtfirst) && this->OBT(obt1) <= upperobt ){
-           //
+         // EMILIANO
-           while ( t2 > pret2+1 ){
+         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) ){
-             // 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());
+             if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
              //
-             valid = 0;
+             MISSING_pkt1 = true;
-             this->HandleTRK_CALIB(false,true);
+             MISSING_pkt2 = false;
+             this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2);
+             ncalib++;
+             if( MISSING_pkt1||MISSING_pkt2||!valid )try_to_recover=true;
              //
-           };
-           //
          };
          //
-       } else if ( this->PKT(pkt2) > this->PKT(pkt1)+1 ){
+         t2++;
-         //
-         // 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
+     };
+ // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ // -----------------------------------------------------------------
+ // 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;
          //
-         obt2 = 0;
+         stringstream oss;
-         pkt2 = 0;
+         oss.str("");
-         valid = 0;
-         this->HandleTRK_CALIB(true,false);
          //
-       };
-       //
+         ////////////////////////////////////////////////////////////////////////
-     } else {
+         // retrieve the name of the current file:
-       //
+         ////////////////////////////////////////////////////////////////////////
-       if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u OUTSIDE the considered time interval \n",fromtime,obt1,pkt1);
+         oss.str("");
-       //      if ( PEDANTIC ) throw -79;
+         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();
-   // we have one more calib pkt2 !
+         TString thisfilename = (TString)row->GetField(0);
-   //
+         if ( IsDebug() ) cout << "Current file ==> "<<thisfilename<<endl;
-   t2++;
-   while ( t2 < nevents2 ){
+         ////////////////////////////////////////////////////////////////////////
-     //
+         // read all the calibrations inserted
-     // handle missing calib1
+         ////////////////////////////////////////////////////////////////////////
-     //
+         oss.str("");
-     if ( IsDebug() ) printf(" t2 is %u nevents2 is %u \n",t2,nevents2);
+         oss << " SELECT ";
-     obt1 = 0;
+         oss << " ID,FROM_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,TO_TIME";
-     pkt1 = 0;
+         oss << " FROM GL_TRK_CALIB ";
-     //
+         oss << " ORDER BY ID DESC LIMIT "<<ncalib<<"; ";
-     tr2->GetEntry(t2);
+         if ( IsDebug() ) cout << oss.str().c_str() << endl;
-     ph2 = eh2->GetPscuHeader();
-     obt2 = ph2->GetOrbitalTime();
+         result = conn->Query(oss.str().c_str());
-     pkt2 = ph2->GetCounter();
+         if ( !result ) throw -4;;
-     //
+         if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl;
-     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 ){
+         // loop over calibrations ...
-     // EMILIANO
+         // -----------------------------------
-     if ( this->PKT(pkt2) >= this->PKT(pktfirst) && this->PKT(pkt2 <= upperpkt) && this->OBT(obt2) >= this->OBT(obtfirst) && this->OBT(obt2) <= upperobt ){
+         UInt_t nn=0;
-       //    if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
+         do {
-       //
+             row = result->Next();
-       if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
+             if(!row)break;
-       //
-       this->HandleTRK_CALIB(false,true);
+             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));
-     t2++;
+             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 5449 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 6120 
 Int_t PamelaDBOperations::ValidateRuns(T
      // 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 << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
+     oss << " SELECT * FROM GL_RUN  WHERE VALIDATION>0 AND RUNHEADER_TIME<="<< t_start
          <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
      if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
      result = conn->Query(oss.str().c_str());
-Line 5592 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 6264 
 Int_t PamelaDBOperations::ValidateRuns(T
          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);
+           //      for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);
+           for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],this_run->VALIDATION);
            nseq=0;
          }
-Line 5613 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 6286 
 Int_t PamelaDBOperations::ValidateRuns(T
      if( CHECK ){
        // check if calibration exists
        if ( IsDebug() )printf("DT %i ===> CHECK Missing calibration\n",interval);
-       Bool_t MISSING = MissingTRK_CALIB(t1,t2);
+       //      Bool_t MISSING = MissingTRK_CALIB(t1,t2);
-       for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],!MISSING);
+       UInt_t MISSING = MissingTRK_CALIB(t1,t2);
+       UInt_t val = 0;
+       if ( MISSING == 1 ) val = 0;
+       if ( MISSING == 0 ) val = 1;
+       if ( MISSING == 2 ) val = 2;
+       for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],val);
        nseq=0;
      };
      //--------------
-Line 5644 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 6322 
 Int_t PamelaDBOperations::ValidateRuns(T
   * @param t2 To absolute time
   * @return true if there might be a missing calibration
   */
- Bool_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){
+ //Bool_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){
+ UInt_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){
    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->Query_GL_TRK_CALIB(t2, conn) )return(true);        //>>> missing
+   if ( trkcalib->Query_GL_TRK_CALIB(t2, conn) )return(1);       //>>> missing
-   if ( trkcalib->TO_TIME  < t2 ) return(true);                                  //>>> missing
+   //  if ( trkcalib->TO_TIME  < t2 ) return(true);                                      //>>> missing
+   if ( trkcalib->TO_TIME  < t2 ) return(1);                                     //>>> missing
    //==============================================================
    // Check is done first on the basis of time between calibration,
-Line 5668 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 6349 
 Bool_t PamelaDBOperations::MissingTRK_CA
      //==============================================================
      Bool_t DOWNLOAD = false;
      // check if the calib was skipped becouse of download .... DA FARE!!
-     if(DOWNLOAD)return(false);
+     //    if(DOWNLOAD)return(false);
+     if(DOWNLOAD)return(0);
-     return(true);                                       //>>> missing
+     //    return(true);                                         //>>> missing
+     return(1);                                  //>>> missing
    };
-Line 5680 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 6363 
 Bool_t PamelaDBOperations::MissingTRK_CA
    //==============================================================
    // the long time interval bewteen runs might be due to download
    if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
-   return(false);
+   if ( trkcalib->VALIDATION ) return(0);
+   if ( IsDebug() )printf("Calibration is not validated... :-/ ==> OK but with VALIDATION=2! \n");
+   return(2);
  };
  /**
-Line 5689 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 6374 
 Bool_t PamelaDBOperations::MissingTRK_CA
   * @param validation true/false
   */
  Int_t PamelaDBOperations::assignVALIDATION(UInt_t idrun, Bool_t validation){
+   return(this->assignVALIDATION(idrun,(UInt_t)validation));
+ }
+ Int_t PamelaDBOperations::assignVALIDATION(UInt_t idrun, UInt_t validation){
    TSQLResult *result = 0;
    stringstream oss;
    oss.str("");
-Line 5885 
 Int_t PamelaDBOperations::removeFile(TSt
+Line 6574 
 Int_t PamelaDBOperations::removeFile(TSt
    stringstream myquery;
    //
    myquery.str("");
-   myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';";
+   myquery << " SELECT ID,ID_RAW,ID_TIMESYNC FROM GL_ROOT where NAME='"<<remfile.Data() <<"';";
    //
    pResult = conn->Query(myquery.str().c_str());
    //
-Line 5901 
 Int_t PamelaDBOperations::removeFile(TSt
+Line 6590 
 Int_t PamelaDBOperations::removeFile(TSt
    //
    this->SetID_ROOT((UInt_t)atoll(Row->GetField(0)));
    this->SetID_RAW((UInt_t)atoll(Row->GetField(1)));
+   UInt_t idtsy=(UInt_t)atoll(Row->GetField(2));
    //
    this->ValidationOFF();
    //
-Line 5908 
 Int_t PamelaDBOperations::removeFile(TSt
+Line 6598 
 Int_t PamelaDBOperations::removeFile(TSt
    //
    this->RemoveRUNS();
    //
-   this->RemoveFILES();
+   this->RemoveFILES(idtsy);
    //
    this->SetID_ROOT(0);
    this->SetID_RAW(0);
-Line 6128 
 void PamelaDBOperations::RemoveRUNS(){
+Line 6818 
 void PamelaDBOperations::RemoveRUNS(){
   * Rearrange calibration tables
   *
   **/
- void PamelaDBOperations::RemoveFILES(){
+ void PamelaDBOperations::RemoveFILES(UInt_t idtsy){
    stringstream myquery;
    //
    myquery.str("");
-Line 6138 
 void PamelaDBOperations::RemoveFILES(){
+Line 6828 
 void PamelaDBOperations::RemoveFILES(){
    //
    conn->Query(myquery.str().c_str());
    //
+   myquery.str("");
+   myquery << " DELETE FROM GL_ROOT WHERE ID=" <<this->GetID_ROOT() <<";";
+   //
+   if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+   //
+   conn->Query(myquery.str().c_str());  //
+   //
+   if ( !chewbacca ){
+     myquery.str("");
+     myquery << " DELETE FROM GL_TIMESYNC WHERE ID=" << idtsy <<";";
+     //
+     if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
+     //
+     conn->Query(myquery.str().c_str());
+   };
+   //
  };
  /**
-Line 6342 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 7048 
 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.2
 


changed lines


 
Added in v.1.20
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.20
-Removed from v.1.2
+Added in v.1.20

	ViewVC Help
Powered by ViewVC 1.1.23