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

Diff of /YodaProfiler/src/PamelaDBOperations.cpp

Parent Directory | Revision Log | View Patch Patch

-revision 1.24 by mocchiut,
Mon Dec  4 11:32:07 2006 UTC
+revision 1.56 by mocchiut,
Mon Mar 10 10:58:01 2008 UTC
 Line 24
  #include <RunHeaderEvent.h>
  #include <RunTrailerEvent.h>
  #include <CalibCalPedEvent.h>
+ #include <CalibCalPulse1Event.h>
+ #include <CalibCalPulse2Event.h>
  #include <CalibS4Event.h>
  #include <CalibTrk1Event.h>
  #include <CalibTrk2Event.h>
-Line 56 
 string getTleDatetime(cTle*);
+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){
+ 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){
    //
+   STATIC=false;
+   if ( staticp ) STATIC=true;
    //
    SetConnection(host,user,password);
    //
-Line 66 
 PamelaDBOperations::PamelaDBOperations(T
+Line 70 
 PamelaDBOperations::PamelaDBOperations(T
    glrun = new GL_RUN();
    //
    if ( !boot ) SetNOBOOT(false);
-   SetBOOTnumber(boot);
+   SetTsync(tsync,gpamela);
-   SetTsync(tsync);
+   SetBOOTnumber(boot,gpamela);
    SetObt0(obt0);
    //
    SetTLEPath(tlefilename);
-Line 79 
 PamelaDBOperations::PamelaDBOperations(T
+Line 83 
 PamelaDBOperations::PamelaDBOperations(T
    INSERT_ROOT = !filerootname.IsNull();
    if( INSERT_ROOT ){
      this->SetRootName(filerootname);
-     this->SetOrbitNo();
+     this->SetOrbitNo(dwinput);
      file = TFile::Open(this->GetRootName().Data());
    } else {
      this->SetRootName("");
-Line 165 
 void PamelaDBOperations::SetAutoBoot(Boo
+Line 169 
 void PamelaDBOperations::SetAutoBoot(Boo
  };
  /**
+  * Set the pedantic flag
+  *
+  */
+ void PamelaDBOperations::SetPedantic(Bool_t dbg){
+   PEDANTIC = dbg;
+ };
+ /**
   * Set the nofrag flag
   *
   */
-Line 177 
 void PamelaDBOperations::SetNoFrag(Bool_
+Line 189 
 void PamelaDBOperations::SetNoFrag(Bool_
   * @param boot    BOOT number of the RAW file
   */
  void PamelaDBOperations::SetBOOTnumber(UInt_t boot){
+   this->SetBOOTnumber(boot,false);
+ };
+ /**
+  * Store the BOOT number of the RAW file.
+  * @param boot    BOOT number of the RAW file
+  */
+ void PamelaDBOperations::SetBOOTnumber(UInt_t boot, Bool_t gpamela){
    BOOTNO=boot;
+   if ( gpamela ){
+     stringstream oss;
+     TSQLResult *result = 0;
+     TSQLRow    *row    = 0;
+     if ( !boot ){
+       //
+       BOOTNO = 1;
+       //
+       // look in the DB for the last timesync and the last run
+       //
+       oss.str("");
+       oss << "SELECT BOOT_NUMBER FROM GL_RUN order by RUNHEADER_TIME desc limit 1;";
+       result = conn->Query(oss.str().c_str());
+       if ( result ){
+         row = result->Next();
+         if ( row ){
+           BOOTNO = (UInt_t)atoll(row->GetField(0)) + 1;
+         };
+       };
+     };
+   };
  };
  /**
-Line 185 
 void PamelaDBOperations::SetBOOTnumber(U
+Line 226 
 void PamelaDBOperations::SetBOOTnumber(U
   * @param boot    time sync
   */
  void PamelaDBOperations::SetTsync(UInt_t ts){
+   this->SetTsync(ts,false);
+ };
+ /**
+  * Store the time sync of the RAW file.
+  * @param boot    time sync
+  */
+ void PamelaDBOperations::SetTsync(UInt_t ts, Bool_t gpamela){
+   //
+   // if not gpamela or given tsync file set ts
+   //
    tsync=ts;
+   if ( gpamela ){
+     stringstream oss;
+     TSQLResult *result = 0;
+     TSQLRow    *row    = 0;
+     TSQLResult *result2 = 0;
+     TSQLRow    *row2    = 0;
+     if ( !ts ){
+       //
+       tsync = 1;
+       //
+       // look in the DB for the last timesync and the last run
+       //
+       oss.str("");
+       oss << "SELECT TIMESYNC FROM GL_TIMESYNC order by TIMESYNC desc limit 1;";
+       result = conn->Query(oss.str().c_str());
+       if ( result ){
+         row = result->Next();
+         if ( row ){
+           tsync = (UInt_t)atoll(row->GetField(0)) + 1;
+           oss.str("");
+           oss << "SELECT (RUNTRAILER_TIME-RUNHEADER_TIME) FROM GL_RUN order by RUNHEADER_TIME desc limit 1;";
+           result2 = conn->Query(oss.str().c_str());
+           if ( result2 ){
+             row2 = result2->Next();
+             if ( row2 ){
+               tsync += (UInt_t)atoll(row2->GetField(0));
+             };
+           }
+         };
+       };
+     };
+   };
  };
  /**
-Line 215 
 void PamelaDBOperations::SetRootName(TSt
+Line 299 
 void PamelaDBOperations::SetRootName(TSt
  /**
   * Store the downlink orbit number from filename.
   */
- void PamelaDBOperations::SetOrbitNo(){
+ void PamelaDBOperations::SetOrbitNo(UInt_t dwinput){
    dworbit = 0;
+   //
+   if ( dwinput ){
+     dworbit = dwinput;
+     if ( IsDebug() ) printf(" Downlink orbit given by hand: %i  \n",dworbit);
+     return;
+   };
+   //
    TString name = this->GetRootFile();
    Int_t nlength = name.Length();
    if ( nlength < 5 ) return;
-Line 254 
 void PamelaDBOperations::SetTLEPath(TStr
+Line 345 
 void PamelaDBOperations::SetTLEPath(TStr
    tlefilename = str;
  };
+ TString PamelaDBOperations::GetRawPath(){
+   if ( STATIC ){
+     return((TString)gSystem->DirName(filerawname.Data())+'/');
+   } else {
+     return((TString)gSystem->ExpandPathName("$PAM_RAW")+'/');
+   };
+ };
+ TString PamelaDBOperations::GetRootPath(){
+   if ( STATIC ){
+     return((TString)gSystem->DirName(filerootname.Data())+'/');
+   } else {
+     return((TString)gSystem->ExpandPathName("$PAM_L0")+'/');
+   };
+ };
  /**
   * Store the olderthan variable
   * @param olderthan
-Line 270 
 Bool_t PamelaDBOperations::SetID_RAW(){
+Line 377 
 Bool_t PamelaDBOperations::SetID_RAW(){
    TSQLResult *result = 0;
    TSQLRow    *row    = 0;
    oss.str("");
-   oss << "SELECT ID FROM GL_RAW WHERE "
+   if ( STATIC ){
-       << " PATH = '" << this->GetRawPath().Data() << "' AND "
+     oss << "SELECT ID FROM GL_RAW WHERE "
-       << " NAME = '" << this->GetRawFile().Data() << "' ";
+         << " PATH = '" << this->GetRawPath().Data() << "' AND "
+         << " NAME = '" << this->GetRawFile().Data() << "' ";
+   } else {
+     oss << "SELECT ID FROM GL_RAW WHERE "
+         << " PATH = '$PAM_RAW' AND "
+         << " NAME = '" << this->GetRawFile().Data() << "' ";
+   }
    result = conn->Query(oss.str().c_str());
    if ( result == NULL ) throw -4;
    row = result->Next();
    if ( !row ) return(false);
-   delete result;
    id = (UInt_t)atoll(row->GetField(0));
+   delete result;
    return(true);
  }
-Line 315 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 428 
 Int_t PamelaDBOperations::SetUpperLimits
    Long64_t deltapkt = 5000LL;
    Long64_t deltaobt = 50000LL;
    //
-   //   pcksList packetsNames;
+   pcksList packetsNames;
-   //   pcksList::iterator Iter;
+   pcksList::iterator Iter;
-   //   getPacketsNames(packetsNames);
+   getPacketsNames(packetsNames);
    //
    pktfirst = 0;
    obtfirst = 0;
+   ppktfirst = 0;
+   pobtfirst = 0;
    //
    TTree *T = 0;
    T = (TTree*)file->Get("Physics");
-Line 334 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 449 
 Int_t PamelaDBOperations::SetUpperLimits
    ph = eh->GetPscuHeader();
    pktfirst = ph->GetCounter();
    obtfirst = ph->GetOrbitalTime();
+   ppktfirst = pktfirst;
+   pobtfirst = obtfirst;
    //
-   //   code = eh->GetCounter();
+   code = eh->GetCounter();
-   //   UInt_t en = 0;
+   UInt_t en = 0;
-   //   for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){
+   for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){
-   //     en = code->Get(GetPacketType(*Iter));
+     en = code->Get(GetPacketType(*Iter));
-   //     if ( en ) printf(" Packet type is %s, entries: %i \n",*Iter,en);
+     if ( !strcmp("CalibCalPed",*Iter) || !strcmp("CalibTrk1",*Iter) || !strcmp("CalibTrk2",*Iter) || !strcmp("CalibS4",*Iter) ){
-   //};
+       if ( IsDebug() ) printf(" We have a calibration before the first physic packet: packet type is %s, entries: %i \n",*Iter,en);
+       //
+       TTree *TC = 0;
+       TC = (TTree*)file->Get("CalibCalPed");
+       if ( !TC || TC->IsZombie() ) throw -16;
+       EventHeader *ehc = 0;
+       PscuHeader *phc = 0;
+       TC->SetBranchAddress("Header", &ehc);
+       TC->GetEntry(0);
+       phc = ehc->GetPscuHeader();
+       pktfirst = phc->GetCounter();
+       obtfirst = phc->GetOrbitalTime();
+       //
+     };
+   };
    //
    T->GetEntry(nevent-1);
    ph = eh->GetPscuHeader();
-Line 351 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 482 
 Int_t PamelaDBOperations::SetUpperLimits
    upperobt = OBT(obtlast);
    upperentry = nevent-1;
    //
-   if ( IsDebug() ) printf(" First entries are: OBT %i pkt_num %i \n",obtfirst,pktfirst);
+   if ( IsDebug() ) printf(" First entries are: OBT %u pkt_num %u \n",obtfirst,pktfirst);
    //
    if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld entry %i\n",upperobt,upperpkt,upperentry);
    //
-   if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(1);
+   if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(32);
    //
-   if ( !nevent ) return(2);
+   if ( !nevent ) return(64);
    //
-   if ( nevent < 2 ) return(4);
+   if ( nevent < 2 ) return(128);
+   if ( nevent < jump ) jump = 1;
+   //  if ( nevent < jump ) jump = int(nevent/10);
+   //  if ( !jump ) jump = 1;
    //
-   if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) ){
+   if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) && nevent > deltapkt ){
      //
+     if ( IsDebug() ) printf(" starting jump %i \n",jump);
+     if ( PEDANTIC ) throw -66;
      // go back
      zomp = nevent - 2;
      //
-Line 394 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 530 
 Int_t PamelaDBOperations::SetUpperLimits
          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){
            zomp = i + jump + 1;
            if ( zomp > nevent-2 ) zomp = nevent - 2;
-           if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
+           if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %u upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
            break;
          };
          //
-Line 438 
 Int_t PamelaDBOperations::SetUpperLimits
+Line 574 
 Int_t PamelaDBOperations::SetUpperLimits
    Long64_t pkth = 0LL;
    Long64_t obth = 0LL;
    //
-   T->GetEntry(upperentry);
+   if ( rhev || rtev ){
-   code = eh->GetCounter();
-   Int_t lasttrail = code->Get(pctp->RunTrailer);
+     T->GetEntry(upperentry);
-   Int_t lasthead = code->Get(pctp->RunHeader);
+     code = eh->GetCounter();
-   if ( lasttrail < rtev ){
+     Int_t lasttrail = code->Get(pctp->RunTrailer);
-     rt->GetEntry(lasttrail);
+     Int_t lasthead = code->Get(pctp->RunHeader);
-     pht = eht->GetPscuHeader();
+     if ( lasttrail < rtev ){
-     pktt = PKT(pht->GetCounter());
+       rt->GetEntry(lasttrail);
-     obtt = OBT(pht->GetOrbitalTime());
+       pht = eht->GetPscuHeader();
-   };
+       pktt = PKT(pht->GetCounter());
-   //
+       obtt = OBT(pht->GetOrbitalTime());
-   if ( lasthead < rhev ){
+     };
-     rh->GetEntry(lasthead);
+     //
-     phh = ehh->GetPscuHeader();
+     if ( lasthead < rhev ){
-     pkth = PKT(phh->GetCounter());
+       rh->GetEntry(lasthead);
-     obth = OBT(phh->GetOrbitalTime());
+       phh = ehh->GetPscuHeader();
-   };
+       pkth = PKT(phh->GetCounter());
-   //
+       obth = OBT(phh->GetOrbitalTime());
-   if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
+     };
-   if ( pkth > upperpkt && obth > upperobt ){
+     //
-     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
+     if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
-     upperpkt = pkth;
+     if ( pkth > upperpkt && obth > upperobt ){
-     upperobt = obth;
+       if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
-     rhev = lasthead+1;
+       upperpkt = pkth;
-   } else {
+       upperobt = obth;
-     rhev = lasthead;
+       rhev = lasthead+1;
-   };
+     } else {
-   if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
+       rhev = lasthead;
-   //
+     };
-   if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+     if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
-   if ( pktt > upperpkt && obtt > upperobt ){
+     //
-     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
+     if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
-     upperpkt = pktt;
+     if ( pktt > upperpkt && obtt > upperobt ){
-     upperobt = obtt;
+       if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
-     rtev = lasttrail+1;
+       upperpkt = pktt;
-   } else {
+       upperobt = obtt;
-     rtev = lasttrail;
+       rtev = lasttrail+1;
-   };
+     } else {
-   if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+       rtev = lasttrail;
-   //  goto kikko;
+     };
-   //
+     if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
-   //
+     //  goto kikko;
-   // Check if runtrailer/runheader are within lower limits
+     //
-   //
+     //
-   //
+     // Check if runtrailer/runheader are within lower limits
-   pkth = 0LL;
+     //
-   obth = 0LL;
+     //
-   spkth = 0LL;
+     pkth = 0LL;
-   sobth = 0LL;
+     obth = 0LL;
-   for (Int_t k=0; k<rhev; k++){
+     spkth = 0LL;
-     if ( k > 0 ){
+     sobth = 0LL;
-       spkth = pkth;
+     for (Int_t k=0; k<rhev; k++){
-       sobth = obth;
+       if ( k > 0 ){
-     };
+         spkth = pkth;
-     rh->GetEntry(k);
+         sobth = obth;
-     phh = ehh->GetPscuHeader();
+       };
-     pkth = PKT(phh->GetCounter());
+       rh->GetEntry(k);
-     obth = OBT(phh->GetOrbitalTime());
+       phh = ehh->GetPscuHeader();
-     //
+       pkth = PKT(phh->GetCounter());
- //    if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth);
+       obth = OBT(phh->GetOrbitalTime());
-     //
+       //
-     if ( pkth < spkth && obth < sobth ){
+       //    if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth);
-       if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);
+       //
-       //
+       if ( pkth < spkth && obth < sobth ){
-       rhev = k-1;
+         if ( PEDANTIC ) throw -66;
-       rh->GetEntry(rhev);
+         if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);
-       pkth = spkth;
+         //
-       obth = sobth;
+         rhev = k-1;
-       //
+         rh->GetEntry(rhev);
-       UInt_t evbefh = 0;
+         pkth = spkth;
-       code = ehh->GetCounter();
+         obth = sobth;
-       evbefh = code->Get(pctp->Physics);
+         //
-       if ( evbefh >= 0 ){
+         UInt_t evbefh = 0;
-         T->GetEntry(evbefh);
+         code = ehh->GetCounter();
-         ph = eh->GetPscuHeader();
+         evbefh = code->Get(pctp->Physics);
-         t_pktlast = PKT(ph->GetCounter());
+         if ( evbefh >= 0 ){
-         t_obtlast = OBT(ph->GetOrbitalTime());
+           T->GetEntry(evbefh);
-         if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump
-           upperpkt = pkth;
-           upperobt = obth;
-           upperentry = evbefh-1;
-         } else {
-           while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){
-             evbefh++;
-             T->GetEntry(evbefh);
-             ph = eh->GetPscuHeader();
-             t_pktlast = PKT(ph->GetCounter());
-             t_obtlast = OBT(ph->GetOrbitalTime());
-           };
-           T->GetEntry(evbefh-1);
            ph = eh->GetPscuHeader();
-           upperpkt = PKT(ph->GetCounter());
+           t_pktlast = PKT(ph->GetCounter());
-           upperobt = OBT(ph->GetOrbitalTime());
+           t_obtlast = OBT(ph->GetOrbitalTime());
-           upperentry = evbefh-1;
+           if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump
-         };
+             upperpkt = pkth;
+             upperobt = obth;
+             upperentry = evbefh-1;
+           } else {
+             while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){
+               evbefh++;
+               T->GetEntry(evbefh);
+               ph = eh->GetPscuHeader();
+               t_pktlast = PKT(ph->GetCounter());
+               t_obtlast = OBT(ph->GetOrbitalTime());
+             };
+             T->GetEntry(evbefh-1);
+             ph = eh->GetPscuHeader();
+             upperpkt = PKT(ph->GetCounter());
+             upperobt = OBT(ph->GetOrbitalTime());
+             upperentry = evbefh-1;
+           };
+         };
+         if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
+         goto kikko0;
        };
-       if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
-       goto kikko0;
-     };
-   };
-  kikko0:
-   //
-   //
-   //
-   pktt = 0LL;
-   obtt = 0LL;
-   spktt = 0LL;
-   sobtt = 0LL;
-   for (Int_t k=0; k<rtev; k++){
-     if ( k > 0 ){
-       spktt = pktt;
-       sobtt = obtt;
      };
-     rt->GetEntry(k);
+   kikko0:
-     pht = eht->GetPscuHeader();
-     pktt = PKT(pht->GetCounter());
-     obtt = OBT(pht->GetOrbitalTime());
      //
- //    if ( IsDebug() ) printf(" k %i rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt);
      //
-     if ( pktt < spktt && obtt < sobtt ){
+     //
-       if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);
+     pktt = 0LL;
+     obtt = 0LL;
+     spktt = 0LL;
+     sobtt = 0LL;
+     for (Int_t k=0; k<rtev; k++){
+       if ( k > 0 ){
+         spktt = pktt;
+         sobtt = obtt;
+       };
+       rt->GetEntry(k);
+       pht = eht->GetPscuHeader();
+       pktt = PKT(pht->GetCounter());
+       obtt = OBT(pht->GetOrbitalTime());
        //
-       rtev = k-1;
+       //    if ( IsDebug() ) printf(" k %i rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt);
-       rt->GetEntry(rtev);
-       pktt = spktt;
-       obtt = sobtt;
-       if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
        //
-       UInt_t evbeft = 0;
+       if ( pktt < spktt && obtt < sobtt ){
-       code = eht->GetCounter();
+         if ( PEDANTIC ) throw -66;
-       evbeft = code->Get(pctp->Physics);
+         if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);
-       if ( evbeft >= 0 ){
+         //
-         T->GetEntry(evbeft);
+         rtev = k-1;
-         ph = eh->GetPscuHeader();
+         rt->GetEntry(rtev);
-         t_pktlast = PKT(ph->GetCounter());
+         pktt = spktt;
-         t_obtlast = OBT(ph->GetOrbitalTime());
+         obtt = sobtt;
-         if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump
+         if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
-           upperpkt = pktt;
+         //
-           upperobt = obtt;
+         UInt_t evbeft = 0;
-           upperentry = evbeft-1;
+         code = eht->GetCounter();
-         } else {
+         evbeft = code->Get(pctp->Physics);
-           while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){
+         if ( evbeft >= 0 ){
-             evbeft++;
+           T->GetEntry(evbeft);
-             T->GetEntry(evbeft);
+           ph = eh->GetPscuHeader();
+           t_pktlast = PKT(ph->GetCounter());
+           t_obtlast = OBT(ph->GetOrbitalTime());
+           if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump
+             upperpkt = pktt;
+             upperobt = obtt;
+             upperentry = evbeft-1;
+           } else {
+             while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){
+               evbeft++;
+               T->GetEntry(evbeft);
+               ph = eh->GetPscuHeader();
+               t_pktlast = PKT(ph->GetCounter());
+               t_obtlast = OBT(ph->GetOrbitalTime());
+             };
+             T->GetEntry(evbeft-1);
              ph = eh->GetPscuHeader();
-             t_pktlast = PKT(ph->GetCounter());
+             upperpkt = PKT(ph->GetCounter());
-             t_obtlast = OBT(ph->GetOrbitalTime());
+             upperobt = OBT(ph->GetOrbitalTime());
+             upperentry = evbeft-1;
            };
-           T->GetEntry(evbeft-1);
-           ph = eh->GetPscuHeader();
-           upperpkt = PKT(ph->GetCounter());
-           upperobt = OBT(ph->GetOrbitalTime());
-           upperentry = evbeft-1;
          };
+         if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+         goto kikko;
+         //      break;
+         //
        };
-       if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+       //
-       goto kikko;
-       //      break;
-       //
      };
      //
-   };
+   kikko:
-   //
+     //
-  kikko:
+     T->GetEntry(upperentry);
-   //
+     code = eh->GetCounter();
-   T->GetEntry(upperentry);
+     lasttrail = code->Get(pctp->RunTrailer);
-   code = eh->GetCounter();
+     lasthead = code->Get(pctp->RunHeader);
-   lasttrail = code->Get(pctp->RunTrailer);
+     if ( lasttrail < rtev ){
-   lasthead = code->Get(pctp->RunHeader);
+       rt->GetEntry(lasttrail);
-   if ( lasttrail < rtev ){
+       pht = eht->GetPscuHeader();
-     rt->GetEntry(lasttrail);
+       pktt = PKT(pht->GetCounter());
-     pht = eht->GetPscuHeader();
+       obtt = OBT(pht->GetOrbitalTime());
-     pktt = PKT(pht->GetCounter());
+     };
-     obtt = OBT(pht->GetOrbitalTime());
+     //
-   };
+     if ( lasthead < rhev ){
-   //
+       rh->GetEntry(lasthead);
-   if ( lasthead < rhev ){
+       phh = ehh->GetPscuHeader();
-     rh->GetEntry(lasthead);
+       pkth = PKT(phh->GetCounter());
-     phh = ehh->GetPscuHeader();
+       obth = OBT(phh->GetOrbitalTime());
-     pkth = PKT(phh->GetCounter());
+     };
-     obth = OBT(phh->GetOrbitalTime());
+     //
-   };
+     if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
-   //
+     if ( pkth > upperpkt && obth > upperobt ){
-   if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
+       if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
-   if ( pkth > upperpkt && obth > upperobt ){
+       upperpkt = pkth;
-     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
+       upperobt = obth;
-     upperpkt = pkth;
+       rhev = lasthead+1;
-     upperobt = obth;
+     } else {
-     rhev = lasthead+1;
+       rhev = lasthead;
-   } else {
+     };
-     rhev = lasthead;
+     if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
-   };
+     //
-   if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
+     if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
-   //
+     if ( pktt > upperpkt && obtt > upperobt ){
-   if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+       if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
-   if ( pktt > upperpkt && obtt > upperobt ){
+       upperpkt = pktt;
-     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
+       upperobt = obtt;
-     upperpkt = pktt;
+       rtev = lasttrail+1;
-     upperobt = obtt;
+     } else {
-     rtev = lasttrail+1;
+       rtev = lasttrail;
-   } else {
+     };
-     rtev = lasttrail;
+     if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
+     //
    };
-   if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
    //
    if ( IsDebug() ) printf(" Upper limits are: OBT %lld pkt_num %lld upper entry %i \n",upperobt,upperpkt,upperentry);
    //
-Line 778 
 void PamelaDBOperations::CheckConnection
+Line 920 
 void PamelaDBOperations::CheckConnection
    if( !conn ) throw -1;
    bool connect = conn->IsConnected();
    if( !connect ) throw -1;
+   //
+   if ( IsDebug() ) printf("\n DB INFORMATIONS:\n SQL: %s Version: %s Host %s Port %i \n\n",conn->GetDBMS(),conn->ServerInfo(),conn->GetHost(),conn->GetPort());
+   //
    if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27;
    //
    // set DB timezone to UTC
-Line 789 
 void PamelaDBOperations::CheckConnection
+Line 934 
 void PamelaDBOperations::CheckConnection
    TSQLResult *result = 0;
    result = conn->Query(oss.str().c_str());
    if ( !result ) throw -10;
+   oss.str("");
+   oss << "SET wait_timeout=173000;";
+   conn->Query(oss.str().c_str());
    //
  };
-Line 797 
 void PamelaDBOperations::CheckConnection
+Line 945 
 void PamelaDBOperations::CheckConnection
   */
  Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){
    //
-   if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u  (UInt_t)(16777214/2)) is %u  \n",pkt_num,pktfirst,(UInt_t)(16777214/2));
+   if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u  (UInt_t)(16777214/2)) is %u  \n",pkt_num,ppktfirst,(UInt_t)(16777214/2));
    //
-   if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2)  ){
+   if ( pkt_num < (ppktfirst/2) && ppktfirst > (16777214/2)  ){
      if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL);
      return((Long64_t)pkt_num+16777215LL);
    };
    //
-   if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){
+   if ( pkt_num > ((Long64_t)ppktfirst*2) && pkt_num > (16777214/2) ){
      if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL);
      return((Long64_t)pkt_num-16777215LL);
    };
-Line 819 
 Long64_t PamelaDBOperations::PKT(UInt_t
+Line 967 
 Long64_t PamelaDBOperations::PKT(UInt_t
   */
  Long64_t PamelaDBOperations::OBT(UInt_t obt){
    //
-   if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((Long64_t)(obt+numeric_limits<UInt_t>::max()));
+   if ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u  (numeric_limits<UInt_t>::max()/2) is %u  \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2));
+   //
+   if ( obt < (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());
+   };
    //
-   if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
+   if ( obt > ((Long64_t)pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
+     if ( IsDebug() ) printf(" pobtfirst*2 %lld \n",((Long64_t)pobtfirst*2));
+     if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
      return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
    };
    //
+   if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)obt);
    return((Long64_t)obt);
  };
-Line 940 
 Int_t PamelaDBOperations::insertPamelaRa
+Line 1096 
 Int_t PamelaDBOperations::insertPamelaRa
    if ( idr ) return(1);
    //
    oss.str("");
-   oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
+   if ( STATIC ){
-       << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";
+     oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
+         << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";
+   } else {
+     oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')";
+   };
    if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
    //
    idr =  this->SetID_RAW();
-Line 960 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1120 
 Int_t PamelaDBOperations::insertPamelaGL
    TSQLResult *result = 0;
    TSQLRow *row = 0;
    UInt_t t0 = 0;
+   Int_t signal = 0;
+   UInt_t idresof = 0;
+   //
+   signal = this->SetUpperLimits();
    //
    stringstream   oss;
    //
    if ( this->GetID_RAW() == 0 )  throw -11;
    //
    oss.str("");
-   oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='"
+   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 SPECIAL_FILE='"
        << this->GetRawFile().Data() << "';";
    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
    result = conn->Query(oss.str().c_str());
-Line 975 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1139 
 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) FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "
+     oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE),ID FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "
          << dworbit << " order by FROM_ORBIT desc limit 1;";
      if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
      result = conn->Query(oss.str().c_str());
-Line 984 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1148 
 Int_t PamelaDBOperations::insertPamelaGL
      if ( !row ) throw -10;
    };
    //
+   idresof = (UInt_t)atoll(row->GetField(6));
+   //
    TTimeStamp tu = TTimeStamp((UInt_t)atoi(row->GetField(0)),(UInt_t)atoi(row->GetField(1)),(UInt_t)atoi(row->GetField(2)),(UInt_t)atoi(row->GetField(3)),(UInt_t)atoi(row->GetField(4)),(UInt_t)atoi(row->GetField(5)),0,true,0);
    t0 = (UInt_t)tu.GetSec();
-   if ( IsDebug() ) printf(" t0 is %u ti is %s %s %s %s %s %s\n",t0,row->GetField(0),row->GetField(1),row->GetField(2),row->GetField(3),row->GetField(4),row->GetField(5));
+   if ( IsDebug() ) printf(" t0 is %u ti is %s %s %s %s %s %s %s\n",t0,row->GetField(0),row->GetField(1),row->GetField(2),row->GetField(3),row->GetField(4),row->GetField(5),row->GetField(6));
    //
    /*
     * Verify that the TIMESYNC have been not already processed
-Line 1008 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1174 
 Int_t PamelaDBOperations::insertPamelaGL
      tsync = (UInt_t)atoll(row->GetField(2));
      obt0 = (UInt_t)atoll(row->GetField(1));
      //
+     if ( PEDANTIC ) throw -67;
      return(1);
    };
    //
    TTree *T = 0;
-   Int_t signal = 0;
    //
    UInt_t nevent = 0;
    UInt_t recEntries = 0;
-Line 1178 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1344 
 Int_t PamelaDBOperations::insertPamelaGL
    //
    if ( !existsts && obt0 ){ // insert timesync by hand
      //
+     if ( PEDANTIC ) throw -68;
      if ( IsDebug() ) printf(" No incl mcmd \n");
      signal = 30;
      //
-Line 1193 
 Int_t PamelaDBOperations::insertPamelaGL
+Line 1360 
 Int_t PamelaDBOperations::insertPamelaGL
    if ( !existsts ) throw -3;
    //
    oss.str("");
-   oss << "INSERT INTO GL_TIMESYNC (ID_RAW,TYPE,OBT0,TIMESYNC) VALUES ('"
+   oss << "INSERT INTO GL_TIMESYNC (ID_RAW,TYPE,OBT0,TIMESYNC,ID_RESURS_OFFSET) VALUES ('"
        << this->GetID_RAW() << "','"//224'"
        << dec << (UInt_t)TYPE << "','"
        << dec << (UInt_t)OBT << "','"
-       << dec << (UInt_t)TSYNC << "');";
+       << dec << (UInt_t)TSYNC << "','"
+       << dec << (UInt_t)idresof << "');";
    conn->Query(oss.str().c_str());
    if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str());
+   if ( conn->GetErrorCode() ){
+     printf(" OK, you got an error because the database structure you are using is not up to date\n Using backward compability code, hence you can continue safetly \n");
+     oss.str("");
+     oss << "INSERT INTO GL_TIMESYNC (ID_RAW,TYPE,OBT0,TIMESYNC) VALUES ('"
+         << this->GetID_RAW() << "','"//224'"
+         << dec << (UInt_t)TYPE << "','"
+         << dec << (UInt_t)OBT << "','"
+         << dec << (UInt_t)TSYNC << "');";
+     conn->Query(oss.str().c_str());
+     if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str());
+   };
    //
    if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
+   //
    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;
    //
    tsync = TSYNC;
-Line 1222 
 Int_t PamelaDBOperations::insertPamelaRo
+Line 1402 
 Int_t PamelaDBOperations::insertPamelaRo
    UInt_t idtimesync = 0;
    //
    oss.str("");
-   oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
+   if ( STATIC ){
-       << " LEFT JOIN GL_ROOT "
+     oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
-       << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
+         << " LEFT JOIN GL_ROOT "
-       << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "
+         << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
-       << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
+         << " 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;
-Line 1253 
 Int_t PamelaDBOperations::insertPamelaRo
+Line 1441 
 Int_t PamelaDBOperations::insertPamelaRo
    idtimesync = (UInt_t)atoll(row->GetField(0));
    //
    oss.str("");
-   oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
+   if ( STATIC ){
-       << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";
+     oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
+         << 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 (conn->Query(oss.str().c_str()) == 0) throw -4;
    //
-Line 1281 
 Int_t PamelaDBOperations::assignBOOT_NUM
+Line 1474 
 Int_t PamelaDBOperations::assignBOOT_NUM
    TSQLResult *result = 0;
    TSQLRow    *row    = 0;
    oss.str("");
-   oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
+   if ( STATIC ){
-       << " PATH = '" << this->GetRawPath().Data() << "' AND "
+     oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
-       << " NAME = '" << this->GetRawFile().Data() << "' ";
+         << " PATH = '" << this->GetRawPath().Data() << "' AND "
+         << " NAME = '" << this->GetRawFile().Data() << "' ";
+   } else {
+     oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
+         << " PATH = '$PAM_RAW' AND "
+         << " NAME = '" << this->GetRawFile().Data() << "' ";
+   };
    result = conn->Query(oss.str().c_str());
    //
    if ( !result ) throw -4;;
-Line 1373 
 Int_t PamelaDBOperations::assignBOOT_NUM
+Line 1572 
 Int_t PamelaDBOperations::assignBOOT_NUM
    //
    Int_t sgn = 0;
    //
-   if ( !found ){
+   if ( !found && !BOOTNO ){
      throw -29;
    } else {
      if ( afound ){
-Line 1402 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1601 
 Int_t PamelaDBOperations::insertPamelaRU
    stringstream oss;
    oss.str("");
    //
-   signal = this->SetUpperLimits();
+   //  signal = this->SetUpperLimits();
    //
    // loop on runheader and runtrailer events
    //
-Line 1420 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1619 
 Int_t PamelaDBOperations::insertPamelaRU
    rt->SetBranchAddress("RunTrailer", &runt);
    rt->SetBranchAddress("Header", &eht);
    //
+   TTree *T = (TTree*)file->Get("Physics");
+   if ( !T || T->IsZombie() ) throw -16;
+   EventHeader *eh = 0;
+   T->SetBranchAddress("Header", &eh);
+   //
+   if ( !(rh->GetEntries()) && !(rt->GetEntries()) && !(T->GetEntries()) ) return(16);
+   //
    UInt_t obtt = 0;
    UInt_t obth = 0;
    UInt_t pktt = 0;
-Line 1453 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1659 
 Int_t PamelaDBOperations::insertPamelaRU
        //
        if ( !ptt &&  !(ptht+1) ){
          //
-         if ( IsDebug() ) printf(" Piece of run at the beginning of the file %i %i %i \n",ptht,pth,ptt);
+         if ( IsDebug() ) printf(" Piece of run at the beginning of the file %u %u %u \n",ptht,pth,ptt);
          //
          this->HandleRunFragments(true,false,0,(evbeft-1));
          //
          //
        } else if ( pth == ptht ){
          //
-         if ( IsDebug() ) printf(" Missing header %i %i %i\n",ptht,pth,ptt);
+         if ( IsDebug() ) printf(" Missing header %u %u %u\n",ptht,pth,ptt);
          //
          if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
          rt->GetEntry(ptt-1);
-Line 1469 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1675 
 Int_t PamelaDBOperations::insertPamelaRU
          rt->GetEntry(ptt);
          pht = eht->GetPscuHeader();
          //
-         if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %i %i %i \n",ptht,pth,ptt);
+         if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %u %u %u \n",ptht,pth,ptt);
-         if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %i %i %i \n",pkth,obth,obtt);
+         if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %u %u %u \n",pkth,obth,obtt);
          //
          this->HandleMissingHoT(true,false,evbefh,evbeft-1);
          //
        } else {
          //
-         if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %i %i %i\n",ptht,pth,ptt);
+         if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %u %u %u\n",ptht,pth,ptt);
          //
          rh->GetEntry(ptht);
          phh = ehh->GetPscuHeader();
-Line 1484 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1690 
 Int_t PamelaDBOperations::insertPamelaRU
          obth = phh->GetOrbitalTime();
          cod = ehh->GetCounter();
          evbefh = cod->Get(pctp->Physics);
-         if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);
+         if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
          //
          // handle this run
          //
-Line 1494 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1700 
 Int_t PamelaDBOperations::insertPamelaRU
          //
          if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){
            //
-           if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER \n");
+           if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER evbefh = %u \n",evbefh);
-           //
-           this->HandleRunFragments(true,true,0,(evbefh-1));
            //
+           if ( evbefh == 0 ) {
+             //
+             signal = 8;
+             if ( IsDebug() ) printf(" Not supported yet: run with no events, no runtrailer, no runheader \n");
+             //
+           } else {
+             //
+             this->HandleRunFragments(true,true,0,(evbefh-1));
+             //
+           };
          };
          //
          //
          if ( (ptht - pth) > 1 ){
            //
            if ( IsDebug() ) printf(" Missing runtrailers! \n");
-           if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %i %i %i \n",ptht,pth,ptt);
+           if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %u %u %u \n",ptht,pth,ptt);
            // is not the consecutive header
            while ( pth != ptht ){
              //
-Line 1525 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1739 
 Int_t PamelaDBOperations::insertPamelaRU
              obth = phh->GetOrbitalTime();
              evbefh = cod->Get(pctp->Physics);
              //
-             if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %i %i %i \n",ptht,pth,ptt);
+             if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %u %u %u \n",ptht,pth,ptt);
-             if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);
+             if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
              //
              this->HandleMissingHoT(false,true,evbefh,evbeft-1);
              //
-Line 1535 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1749 
 Int_t PamelaDBOperations::insertPamelaRU
          } else if ( !(ptht - pth) ){
            //
            if ( IsDebug() ) printf(" Missing runheader! \n");
-           if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %i %i %i \n",ptht,pth,ptt);
+           if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %u %u %u \n",ptht,pth,ptt);
            if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n");
            //
          } else {
-Line 1556 
 Int_t PamelaDBOperations::insertPamelaRU
+Line 1770 
 Int_t PamelaDBOperations::insertPamelaRU
            obth = phh->GetOrbitalTime();
            cod = ehh->GetCounter();
            evbefh = cod->Get(pctp->Physics);
-           if ( IsDebug() ) printf(" Piece of run at the end of file %i %i %i \n",pkth,obth,obtt);
+           if ( IsDebug() ) printf(" Piece of run at the end of file %u %u %u \n",pkth,obth,obtt);
-           if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i %i %i \n",ptht,pth,ptt);
+           if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u %u %u \n",ptht,pth,ptt);
-           if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i \n",rhev);
+           if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u \n",rhev);
            //
            this->HandleRunFragments(false,true,evbefh,upperentry);
          } else {
-Line 1606 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1820 
 Bool_t PamelaDBOperations::IsRunAlreadyI
    //                                  ( runhead_time <= _our_runhead_time &&  runtrail_time >= _our_runtrail_time) &&
    //                                  ( runhead_obt <= _our_runheadobt || runhead_pkt <= _our_runheadpkt ) &&
    //                                  ( runtrail_obt <= _our_runtrailobt || runtrail_pkt <= _our_runtrailpkt ) )
+   //                                  ||
+   //                                  ( runhead_time = _our_runhead_time &&  runtrail_time = _our_runtrail_time && nevents > 100 )
+   //                                  )
+   //                                  ||
+   //                                  ( runhead_time = _our_runhead_time &&  runtrail_time > _our_runtrail_time && nevents > 100 )
+   //                                  )
+   //                                  ||
+   //                                  ( runhead_time < _our_runhead_time &&  runtrail_time = _our_runtrail_time && nevents > 100 )
    //                                  )
    //
    oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE "
-Line 1621 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1843 
 Bool_t PamelaDBOperations::IsRunAlreadyI
        << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "
        << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("
        << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "
-       << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));";
+       << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") OR "
+       << " (RUNHEADER_TIME=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "                //  these two lines in a certain way disable the patch below...
+       << " RUNTRAILER_TIME=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100) OR"    //
+       << " (RUNHEADER_TIME=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "                //
+       << " RUNTRAILER_TIME>" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100) OR"    //
+       << " (RUNHEADER_TIME<" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "                //
+       << " RUNTRAILER_TIME=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100)"    //
+       << " ));";
    //
    if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());
    result = conn->Query(oss.str().c_str());
-Line 1643 
 Bool_t PamelaDBOperations::IsRunAlreadyI
+Line 1872 
 Bool_t PamelaDBOperations::IsRunAlreadyI
      //
      // the run has already been inserted
      //
-     //    return(true);  //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here
+     if ( signal && IsDebug() ) printf(" The run has already been inserted \n");
+     return(true);  //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here
      //
      // PATCH!
      // we keep the processing run if (in order of growing importance) 1) we have the runtrailer while the old run doesn't have it 2) we have the runheader
-Line 1821 
 void PamelaDBOperations::HandleRunFragme
+Line 2051 
 void PamelaDBOperations::HandleRunFragme
      //
    } else {
      if ( IsDebug() ) printf(" The run is already present in the fragment table \n");
+     if ( PEDANTIC ) throw -69;
      return;
    };
    //
-Line 1958 
 void PamelaDBOperations::HandleRunFragme
+Line 2189 
 void PamelaDBOperations::HandleRunFragme
        ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT());
        UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT());
        ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT());
-       if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
+       if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
        TTree *T= 0;
        T = (TTree*)file->Get("Physics");
        if ( !T || T->IsZombie() ) throw -16;
-Line 1971 
 void PamelaDBOperations::HandleRunFragme
+Line 2202 
 void PamelaDBOperations::HandleRunFragme
          bpkt = PKT(ph->GetCounter());
          bobt = OBT(ph->GetOrbitalTime());
          firstev++;
+         if ( PEDANTIC ) throw -71;
        };
-       if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
+       if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
        //
        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
-Line 2155 
 void PamelaDBOperations::HandleRunFragme
+Line 2387 
 void PamelaDBOperations::HandleRunFragme
        ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT());
        UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT());
        ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT());
-       if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
+       if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
        TTree *T= 0;
        T = (TTree*)file->Get("Physics");
        if ( !T || T->IsZombie() ) throw -16;
-Line 2168 
 void PamelaDBOperations::HandleRunFragme
+Line 2400 
 void PamelaDBOperations::HandleRunFragme
          apkt = PKT(ph->GetCounter());
          aobt = OBT(ph->GetOrbitalTime());
          lastev--;
+         if ( PEDANTIC ) throw -72;
        };
-       if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
+       if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
        //
        glrun->SetEV_TO(lastev);
        glrun->SetNEVENTS(lastev-firstev+1);
-Line 2254 
 void PamelaDBOperations::HandleRunFragme
+Line 2487 
 void PamelaDBOperations::HandleRunFragme
      //
      if ( row ){
        if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n");
+       if ( PEDANTIC ) throw -70;
      } else {
        if ( NoFrag() ){
          glrun->SetID_RUN_FRAG(glrun->GetID());
-Line 2776 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 3010 
 Int_t PamelaDBOperations::insertCALO_CAL
          fromtime = this->GetAbsTime(ph->GetOrbitalTime());
          if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
            //
-           if ( IsDebug() ) printf(" Calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);
+           if ( IsDebug() ) printf(" Calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
            //
            // check if the calibration has already been inserted
            //
-Line 2797 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 3031 
 Int_t PamelaDBOperations::insertCALO_CAL
            if ( row ){
              //
              if ( IsDebug() ) printf(" Calo calibration already inserted in the DB\n");
+             if ( PEDANTIC ) throw -73;
              //
            } else {
              //
-Line 2821 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 3056 
 Int_t PamelaDBOperations::insertCALO_CAL
                //
                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
                if ( fromtime < 1150871000 ){ //1150866904
-                 if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
+                 if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
                };
                //
-Line 2883 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 3118 
 Int_t PamelaDBOperations::insertCALO_CAL
            //
          } else {
            //
-           if ( IsDebug() ) printf(" Repetead calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);
+           if ( IsDebug() ) printf(" Repeated calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
+           if ( PEDANTIC ) throw -74;
            //
          };
          //
-Line 2894 
 Int_t PamelaDBOperations::insertCALO_CAL
+Line 3130 
 Int_t PamelaDBOperations::insertCALO_CAL
    return(0);
  };
+ /**
+  * Scan calorimeter calibrations packets, fill the GL_CALO_CALIB table
+  */
+ Int_t PamelaDBOperations::insertCALOPULSE_CALIB(){
+   //
+   TSQLResult *result = 0;
+   TSQLRow    *row    = 0;
+   //
+   stringstream oss;
+   oss.str("");
+   //
+   oss << " DESCRIBE GL_CALOPULSE_CALIB;";
+   if ( IsDebug() ) printf(" Check if the GL_CALOPULSE_CALIB table exists: query is \n %s \n",oss.str().c_str());
+   result = conn->Query(oss.str().c_str());
+   //
+   if ( conn->GetErrorCode() ){
+     if ( IsDebug() ) printf(" The GL_CALOPULSE_CALIB table does not exists! \n");
+     throw -30;
+   };
+   //
+   // CaloPulse1
+   //
+   CalibCalPulse1Event *cp1 = 0;
+   TTree *tr = 0;
+   EventHeader *eh = 0;
+   PscuHeader *ph = 0;
+   //
+   UInt_t nevents = 0;
+   UInt_t fromtime = 0;
+   UInt_t totime = 0;
+   UInt_t obt = 0;
+   UInt_t pkt = 0;
+   //
+   tr = (TTree*)file->Get("CalibCalPulse1");
+   if ( !tr || tr->IsZombie() ) throw -31;
+   //
+   tr->SetBranchAddress("CalibCalPulse1", &cp1);
+   tr->SetBranchAddress("Header", &eh);
+   nevents = tr->GetEntries();
+   //
+   if ( nevents > 0 ){
+     //
+     for (UInt_t i=0; i < nevents; i++){
+       tr->GetEntry(i);
+       for (UInt_t section = 0; section < 4; section++){
+         //
+         if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){
+           valid = 1;
+           if ( cp1->pperror[section] ) valid = 0;
+           ph = eh->GetPscuHeader();
+           obt = ph->GetOrbitalTime();
+           pkt = ph->GetCounter();
+           fromtime = this->GetAbsTime(ph->GetOrbitalTime());
+           if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
+             //
+             if ( IsDebug() ) printf(" Calo pulse1 calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
+             //
+             // check if the calibration has already been inserted
+             //
+             oss.str("");
+             oss << " SELECT ID FROM GL_CALOPULSE_CALIB WHERE "
+                 << " SECTION = "<< section << " AND "
+                 << " PULSE_AMPLITUDE = 0 AND "
+                 << " BOOT_NUMBER = "<< this->GetBOOTnumber() << " AND "
+                 << " OBT = "<< obt << " AND "
+                 << " PKT = "<< pkt << ";";
+             //
+             if ( IsDebug() ) printf(" Check if the calo pulse1 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(" Calo pulse1 calibration already inserted in the DB\n");
+               if ( PEDANTIC ) throw -75;
+               //
+             } else {
+               //
+               // we have to insert a new calibration, check where to place it
+               //
+               oss.str("");
+               oss << " SELECT ID,TO_TIME FROM GL_CALOPULSE_CALIB WHERE "
+                   << " SECTION = "<< section << " AND "
+                   << " PULSE_AMPLITUDE = 0 AND "
+                   << " SECTION = "<< section << " AND "
+                   << " FROM_TIME < "<< fromtime << " AND "
+                   << " TO_TIME > "<< fromtime << ";";
+               //
+               if ( IsDebug() ) printf(" Check where to place the pulse1 calo 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 ){
+                 //
+                 // no calibrations in the db contain our calibration
+                 //
+                 if ( IsDebug() ) printf(" Pulse1 calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
+                 if ( fromtime < 1150871000 ){ //1150866904
+                   if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
+                   fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
+                 };
+                 //
+                 oss.str("");
+                 oss << " SELECT FROM_TIME FROM GL_CALOPULSE_CALIB WHERE "
+                     << " PULSE_AMPLITUDE = 0 AND "
+                     << " SECTION = "<< section << " AND "
+                     << " FROM_TIME > "<< fromtime << " ORDER BY FROM_TIME ASC LIMIT 1;";
+                 //
+                 if ( IsDebug() ) printf(" Check the upper limit for pulse1 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_CALOPULSE_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 pulse1 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_CALOPULSE_CALIB (ID,ID_ROOT_L0,EV_ROOT,FROM_TIME,TO_TIME,SECTION,PULSED_STRIP,PULSE_AMPLITUDE,OBT,PKT,BOOT_NUMBER,VALIDATION) "
+                   << " VALUES (NULL,' "
+                   << idroot << "','"
+                   << i << "','"
+                   << fromtime << "','"
+                   << totime << "','"
+                   << section << "',NULL,'0','"
+                   << obt << "','"
+                   << pkt << "','"
+                   << this->GetBOOTnumber() << "','"
+                   << valid << "');";
+               //
+               if ( IsDebug() ) printf(" Insert the new pulse1 calibration: query is \n %s \n",oss.str().c_str());
+               //
+               result = conn->Query(oss.str().c_str());
+               //
+               if ( !result ) throw -4;
+               //
+             };
+             //
+           } else {
+             //
+             if ( IsDebug() ) printf(" Repeated pulse1 calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
+             if ( PEDANTIC ) throw -76;
+             //
+           };
+           //
+         };
+       };
+     };
+   };
+   //
+   // CaloPulse2
+   //
+   tr->Reset();
+   CalibCalPulse2Event *cp2 = 0;
+   tr = 0;
+   //
+   nevents = 0;
+   fromtime = 0;
+   totime = 0;
+   obt = 0;
+   pkt = 0;
+   //
+   tr = (TTree*)file->Get("CalibCalPulse2");
+   if ( !tr || tr->IsZombie() ) throw -32;
+   //
+   tr->SetBranchAddress("CalibCalPulse2", &cp2);
+   tr->SetBranchAddress("Header", &eh);
+   nevents = tr->GetEntries();
+   //
+   if ( nevents > 0 ){
+     //
+     for (UInt_t i=0; i < nevents; i++){
+       tr->GetEntry(i);
+       for (UInt_t section = 0; section < 4; section++){
+         //
+         if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){
+           valid = 1;
+           if ( cp2->pperror[section] ) valid = 0;
+           ph = eh->GetPscuHeader();
+           obt = ph->GetOrbitalTime();
+           pkt = ph->GetCounter();
+           fromtime = this->GetAbsTime(ph->GetOrbitalTime());
+           if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
+             //
+             if ( IsDebug() ) printf(" Calo pulse2 calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
+             //
+             // check if the calibration has already been inserted
+             //
+             oss.str("");
+             oss << " SELECT ID FROM GL_CALOPULSE_CALIB WHERE "
+                 << " SECTION = "<< section << " AND "
+                 << " PULSE_AMPLITUDE != 0 AND "
+                 << " BOOT_NUMBER = "<< this->GetBOOTnumber() << " AND "
+                 << " OBT = "<< obt << " AND "
+                 << " PKT = "<< pkt << ";";
+             //
+             if ( IsDebug() ) printf(" Check if the calo pulse2 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(" Calo pulse2 calibration already inserted in the DB\n");
+               if ( PEDANTIC ) throw -77;
+               //
+             } else {
+               //
+               // we have to insert a new calibration
+               //
+               //
+               // Determine the amplitude of the pulse
+               //
+               UInt_t pampli = 1;
+               UInt_t pstrip = 0;
+               UInt_t se = 0;
+               if ( section == 1 ) se = 2;
+               if ( section == 2 ) se = 3;
+               if ( section == 3 ) se = 1;
+               for (Int_t ii=0;ii<16;ii++){
+                 if ( cp2->calpuls[se][0][ii] > 10000. ){
+                   pampli = 2;
+                   pstrip = ii;
+                 };
+               };
+               if ( pampli == 1 ){
+                 Bool_t found = false;
+                 Float_t delta=0.;
+                 UInt_t cstr = 0;
+                 while ( !found && cstr < 16 ){
+                   for (Int_t ii=0;ii<16;ii++){
+                     delta = cp2->calpuls[se][0][ii] - cp2->calpuls[se][0][cstr];
+                     if ( IsDebug() ) printf(" cstr is %u ii is %i delta is %f \n",cstr,ii,delta);
+                     if ( delta > 500. ){
+                       pampli = 1;
+                       pstrip = ii;
+                       found = true;
+                       if ( IsDebug() ) printf(" FOUND cstr is %u ii is %i delta is %f \n",cstr,ii,delta);
+                     };
+                   };
+                   cstr++;
+                 };
+                 if ( !found ) pstrip = 100;
+               };
+               if ( IsDebug() ) printf(" The amplitude of the pulser is %u (where 1 = low pulse, 2 = high pulse), pulsed strip is %u \n",pampli,pstrip);
+               //
+               // we have to insert a new calibration, check where to place it
+               //
+               oss.str("");
+               oss << " SELECT ID,TO_TIME FROM GL_CALOPULSE_CALIB WHERE "
+                   << " SECTION = "<< section << " AND "
+                   << " PULSE_AMPLITUDE = " << pampli << " AND "
+                   << " SECTION = "<< section << " AND "
+                   << " FROM_TIME < "<< fromtime << " AND "
+                   << " TO_TIME > "<< fromtime << ";";
+               //
+               if ( IsDebug() ) printf(" Check where to place the pulse2 calo 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 ){
+                 //
+                 // no calibrations in the db contain our calibration
+                 //
+                 if ( IsDebug() ) printf(" Pulse2 calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
+                 if ( fromtime < 1150871000 ){ //1150866904
+                   if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
+                   fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
+                 };
+                 //
+                 oss.str("");
+                 oss << " SELECT FROM_TIME FROM GL_CALOPULSE_CALIB WHERE "
+                     << " PULSE_AMPLITUDE = " << pampli << " AND "
+                     << " SECTION = "<< section << " AND "
+                     << " FROM_TIME > "<< fromtime << " ORDER BY FROM_TIME ASC LIMIT 1;";
+                 //
+                 if ( IsDebug() ) printf(" Check the upper limit for pulse2 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_CALOPULSE_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 pulse2 calibration: query is \n %s \n",oss.str().c_str());
+                 result = conn->Query(oss.str().c_str());
+                 //
+                 if ( !result ) throw -4;
+                 //
+               };
+               //
+               // Fill the DB
+               //
+               oss.str("");
+               //              oss << " INSERT INTO GL_CALOPULSE_CALIB (ID,ID_ROOT_L0,EV_ROOT,FROM_TIME,TO_TIME,SECTION,OBT,PKT,BOOT_NUMBER,VALIDATION) "
+               oss << " INSERT INTO GL_CALOPULSE_CALIB (ID,ID_ROOT_L0,EV_ROOT,FROM_TIME,TO_TIME,SECTION,PULSED_STRIP,PULSE_AMPLITUDE,OBT,PKT,BOOT_NUMBER,VALIDATION) "
+                   << " VALUES (NULL,' "
+                   << idroot << "','"
+                   << i << "','"
+                   << fromtime << "','"
+                   << totime << "','"
+                   << section << "','"
+                   << pstrip << "','"
+                   << pampli << "','"
+                   << obt << "','"
+                   << pkt << "','"
+                   << this->GetBOOTnumber() << "','"
+                   << valid << "');";
+               //
+               if ( IsDebug() ) printf(" Insert the new pulse2 calibration: query is \n %s \n",oss.str().c_str());
+               //
+               result = conn->Query(oss.str().c_str());
+               //
+               if ( !result ) throw -4;
+               //
+             };
+             //
+           } else {
+             //
+             if ( IsDebug() ) printf(" Repeated pulse2 calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
+             if ( PEDANTIC ) throw -78;
+             //
+           };
+           //
+         };
+       };
+     };
+   };
+   //
+   return(0);
+ };
  /**
   * Fill the GL_TRK_CALIB table
   */
-Line 2943 
 void PamelaDBOperations::HandleTRK_CALIB
+Line 3559 
 void PamelaDBOperations::HandleTRK_CALIB
    if ( row ){
      //
      if ( IsDebug() ) printf(" Trk calibration already inserted in the DB\n");
+     if ( PEDANTIC ) throw -80;
      //
    } else {
      //
-Line 3101 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3718 
 Int_t PamelaDBOperations::insertTRK_CALI
      //
      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
        //
-       if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);
+       if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
        //
        valid = ValidateTrkCalib( caltrk1, eh1 );
        if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
-Line 3144 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3761 
 Int_t PamelaDBOperations::insertTRK_CALI
          //
        };
        //
-       if ( IsDebug() ) printf(" Found trk calibration2 at obt %i pkt %i t2 is %i \n",obt2,pkt2,t2);
+       if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2);
        //
        // The calibration is good
        //
        if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){
          //
-         if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2);
+         if ( IsDebug() ) printf(" The trk calibration2 at obt %u pkt %u t2 is %u is good \n",obt2,pkt2,t2);
          //
          UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
          if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
-Line 3164 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3781 
 Int_t PamelaDBOperations::insertTRK_CALI
          //
          if ( t2 != pret2+1 ){
            //
-           if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i pret2 is %i \n",obt2,pkt2,t2,pret2);
+           if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u pret2 is %u \n",obt2,pkt2,t2,pret2);
            //
            while ( t2 > pret2+1 ){
              //
-Line 3193 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3810 
 Int_t PamelaDBOperations::insertTRK_CALI
          //
          // Check for missing calibtrk2
          //
-         if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);
+         if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
          t2 = pret2;
          //
          // handle missing calib2
-Line 3207 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3824 
 Int_t PamelaDBOperations::insertTRK_CALI
        //
      } else {
        //
-       if ( IsDebug() ) printf(" Repetead trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);
+       if ( IsDebug() ) printf(" Repetead trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
+       if ( PEDANTIC ) throw -79;
        //
      };
      //
-Line 3232 
 Int_t PamelaDBOperations::insertTRK_CALI
+Line 3850 
 Int_t PamelaDBOperations::insertTRK_CALI
      valid = 0;
      if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
        //
-       if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);
+       if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
        //
        this->HandleTRK_CALIB(false,true);
        //
-Line 3286 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 3904 
 Int_t PamelaDBOperations::insertS4_CALIB
      fromtime = this->GetAbsTime(ph->GetOrbitalTime());
      if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
        //
-       if ( IsDebug() ) printf(" S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);
+       if ( IsDebug() ) printf(" S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
        //
        // check if the calibration has already been inserted
        //
-Line 3306 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 3924 
 Int_t PamelaDBOperations::insertS4_CALIB
        if ( row ){
          //
          if ( IsDebug() ) printf(" S4 calibration already inserted in the DB\n");
+         if ( PEDANTIC ) throw -81;
          //
        } else {
          //
-Line 3329 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 3948 
 Int_t PamelaDBOperations::insertS4_CALIB
            //
            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");
            if ( fromtime < 1150871000 ){
-             if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
+             if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
            };
            //
-Line 3388 
 Int_t PamelaDBOperations::insertS4_CALIB
+Line 4007 
 Int_t PamelaDBOperations::insertS4_CALIB
        //
      } else {
        //
-       if ( IsDebug() ) printf(" Repetead S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);
+       if ( IsDebug() ) printf(" Repeated S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
+       if ( PEDANTIC ) throw -82;
        //
      };
      //
-Line 3525 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
+Line 4145 
 Int_t PamelaDBOperations::CleanGL_RUN_FR
          //
        } else {
          if ( IsDebug() ) printf(" The already exist in the GL_RUN table! \n");
+         if ( PEDANTIC ) throw -83;
        };
        if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0));
        //
-Line 3599 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 4220 
 Int_t PamelaDBOperations::ValidateRuns(T
        row = result->Next();
        t_stop = (UInt_t)atoll(row->GetField(4));
      };
-     if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);
+     if ( IsDebug() ) printf("Validation interval: from time %u - to time %u \n\n",t_stop,t_start);
      // --------------------------------------------------------------
      // now retrieves runs to be validated
      // --------------------------------------------------------------
-Line 3770 
 Int_t PamelaDBOperations::ValidateRuns(T
+Line 4391 
 Int_t PamelaDBOperations::ValidateRuns(T
                    }else printf("ValidateRuns ***WARNING*** : run sequence exceed assumed size (%i) \n",nseq_max);
            };
-           if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID);
+           if ( IsDebug() ) printf("%i Run %u \n",nrow,this_run->ID);
            nrow++;
    };
-Line 3799 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 4420 
 Bool_t PamelaDBOperations::MissingTRK_CA
          // which should be equal to the time between ascending-nodes.
          //==============================================================
          if ( t2 - trkcalib->FROM_TIME > 5700) {
-                 if ( IsDebug() )printf("Long time between calib and run start %i :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);
+                 if ( IsDebug() )printf("Long time between calib and run start %u :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);
          //==============================================================
          // there might be a missing calibration, due to:
          // - MM full
-Line 3820 
 Bool_t PamelaDBOperations::MissingTRK_CA
+Line 4441 
 Bool_t PamelaDBOperations::MissingTRK_CA
          // it is enough to say that there are no missing calibrations
          //==============================================================
          // the long time interval bewteen runs might be due to download
-         if ( IsDebug() )printf("Short time between calib and run start %i :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
+         if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
          return(false);
  };
-Line 4136 
 void PamelaDBOperations::ValidationOFF()
+Line 4757 
 void PamelaDBOperations::ValidationOFF()
        };
      };
    };
-   if ( IsDebug() ) printf(" %i runs have been unvalidated \n",unv);
+   if ( IsDebug() ) printf(" %u runs have been unvalidated \n",unv);
  };
  /**
-Line 4167 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4788 
 void PamelaDBOperations::RemoveRUNS(){
      if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");
      while ( Row ){
        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
-       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+       if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
        drun++;
        Row = pResult->Next();
      };
-Line 4188 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4809 
 void PamelaDBOperations::RemoveRUNS(){
    } else {
      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");
      while ( Row ){
-       if ( IsDebug() ) printf(" restore run %i \n",(UInt_t)atoll(Row->GetField(1)));
+       if ( IsDebug() ) printf(" restore run %u \n",(UInt_t)atoll(Row->GetField(1)));
        delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");
-       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(1)));
+       if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(1)));
        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");
        if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){
-         if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+         if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
        };
        drun++;
-Line 4219 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4840 
 void PamelaDBOperations::RemoveRUNS(){
    } else {
      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");
      while ( Row ){
-       if ( IsDebug() ) printf(" del run idtrash %i \n",(UInt_t)atoll(Row->GetField(0)));
+       if ( IsDebug() ) printf(" del run idtrash %u \n",(UInt_t)atoll(Row->GetField(0)));
        myquery.str("");
        myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";
        conn->Query(myquery.str().c_str());
-Line 4228 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4849 
 void PamelaDBOperations::RemoveRUNS(){
      };
    };
    //
-   if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_TRASH table \n",drun);
+   if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_TRASH table \n",drun);
    //
    //
    //
-Line 4246 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4867 
 void PamelaDBOperations::RemoveRUNS(){
    } else {
      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");
      while ( Row ){
-       if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
+       if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
        myquery.str("");
        myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";
        conn->Query(myquery.str().c_str());
-Line 4255 
 void PamelaDBOperations::RemoveRUNS(){
+Line 4876 
 void PamelaDBOperations::RemoveRUNS(){
      };
    };
    //
-   if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_FRAGMENTS table \n",drun);
+   if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_FRAGMENTS table \n",drun);
    //
    //
    //
-Line 4325 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 4946 
 void PamelaDBOperations::RemoveCALIBS(){
        //
      };
    };
+   Bool_t OLDDB = false;
+   for (Int_t section = 0; section < 4; section++){
+     myquery.str("");
+     myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_CALOPULSE_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<" AND ";
+     myquery << " SECTION=" << section << ";";
+     //
+     pResult = conn->Query(myquery.str().c_str());
+     //
+     if ( conn->GetErrorCode() ){
+       printf(" Section %i : warning, old databse structure no GL_CALOPULSE_CALIB table!\n",section);
+       OLDDB=true;
+     } else {
+       Row = pResult->Next();
+       if( !Row->GetField(0) || !Row->GetField(1) ){
+         //
+         if ( IsDebug() ) printf(" NO PULSE CALO CALIBRATION SECTION %i ASSOCIATED TO THIS FILE! \n",section);
+         //
+       } else {
+         //
+         myquery.str("");
+         myquery << " UPDATE GL_CALOPULSE_CALIB SET TO_TIME=" << Row->GetField(1);
+         myquery << " WHERE TO_TIME="<< Row->GetField(0) << " AND ";
+         myquery << " SECTION=" << section << ";";
+         //
+         pResult = conn->Query(myquery.str().c_str());
+         //
+         if( !pResult ){
+           //
+           if ( IsDebug() ) printf(" ERROR DELETING CALO PULSE CALIBRATIONS \n");
+           //
+           throw -4;
+           //
+         };
+         //
+       };
+     };
+   };
    myquery.str("");
-   myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+   myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT() << ";";
    //
    pResult = conn->Query(myquery.str().c_str());
    //
-Line 4338 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 4996 
 void PamelaDBOperations::RemoveCALIBS(){
      //
    };
    //
+   myquery.str("");
+   myquery << " DELETE FROM GL_CALOPULSE_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT() << ";";
+   //
+   pResult = conn->Query(myquery.str().c_str());
+   if ( IsDebug() ) printf(" Delete from GL_CALOPULSE_CALIB query is %s \n",myquery.str().c_str());
+   if ( !OLDDB ){
+     //
+     if( !pResult ){
+       //
+       if ( IsDebug() ) printf(" ERROR DELETING PULSE CALO CALIBRATIONS \n");
+       //
+       throw -4;
+       //
+     };
+   };
+   //
    // Tracker
    //
    myquery.str("");
-Line 4367 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 5041 
 void PamelaDBOperations::RemoveCALIBS(){
      };
      //
      myquery.str("");
-     myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+     myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT() << ";";
      //
      pResult = conn->Query(myquery.str().c_str());
      //
-Line 4410 
 void PamelaDBOperations::RemoveCALIBS(){
+Line 5084 
 void PamelaDBOperations::RemoveCALIBS(){
      };
      //
      myquery.str("");
-     myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
+     myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT() << ";";
      //
      pResult = conn->Query(myquery.str().c_str());
      //
-Line 4549 
 UInt_t PamelaDBOperations::ValidateTrkCa
+Line 5223 
 UInt_t PamelaDBOperations::ValidateTrkCa
   *
   **/
  UInt_t PamelaDBOperations::Check(){
+         return(this->Check(0,0));
+ }
+ UInt_t PamelaDBOperations::Check(UInt_t from, UInt_t to){
+   //
+         if ( IsDebug() ) printf(" from %u to %u \n",from,to);
    //
    UInt_t test = 0;
    //
-Line 4567 
 UInt_t PamelaDBOperations::Check(){
+Line 5247 
 UInt_t PamelaDBOperations::Check(){
    TSQLRow    *row    = 0;
    TSQLResult *result2 = 0;
    TSQLRow    *row2    = 0;
+   TSQLResult *result3 = 0;
+   TSQLRow    *row3    = 0;
    oss.str("");
-   oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN order by RUNHEADER_TIME asc;";
+   oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME,NEVENTS FROM GL_RUN order by RUNHEADER_TIME asc;";
    //  oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;";
    result = conn->Query(oss.str().c_str());
    //
    if ( !result ) throw -4;;
    //
    row = result->Next();
+   UInt_t nid = 0;
    //
    while ( row ){
+     nid++;
+     if ( !(nid%1000) && nid ) printf(" %iK run scanned \n",nid/1000);
      thisid = (UInt_t)atoll(row->GetField(0));
      thisl0id = (UInt_t)atoll(row->GetField(1));
      thisrht = (UInt_t)atoll(row->GetField(2));
      thisrtt = (UInt_t)atoll(row->GetField(3));
      //
+     if ( from > 0 && nid <= from ) goto ss;
+     if ( to > 0 && nid >= to ) goto ss;
+ //
+     if ( (UInt_t)atoll(row->GetField(4)) > 1 ){
+       //
+       //
+       //
+       oss.str("");
+       oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
+           << thisid << " AND ( RUNHEADER_TIME="
+           << thisrht << " OR RUNTRAILER_TIME="
+           << thisrtt << " ) AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;";
+       result3 = conn->Query(oss.str().c_str());
+       if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str());
+       if ( result3 ){
+         //
+         oss.str("");
+         oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
+             << thisid << " AND RUNHEADER_TIME="
+             << thisrht << " AND RUNTRAILER_TIME!="
+             << thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;";
+         result3 = conn->Query(oss.str().c_str());
+         if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str());
+         if ( result3 ){
+           row3 = result3->Next();
+           //
+           while ( row3 ){
+             //
+             // 2 runs with same runheader
+             //
+             printf(" CHECK n.4 RUNs %u and %u HAVE SAME RUNHEADER \n",thisid,(UInt_t)atoll(row3->GetField(0)));
+             row3 = result3->Next();
+           };
+ //      delete result3;
+         };
+         //
+         oss.str("");
+         oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
+             << thisid << " AND RUNHEADER_TIME!="
+             << thisrht << " AND RUNTRAILER_TIME="
+             << thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;";
+         result3 = conn->Query(oss.str().c_str());
+         if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str());
+         if ( result3 ){
+           row3 = result3->Next();
+           //
+           while ( row3 ){
+             //
+             // 2 runs with same runtrailer
+             //
+             printf(" CHECK n.5 RUNs %u and %u HAVE SAME RUNTRAILER \n",thisid,(UInt_t)atoll(row3->GetField(0)));
+             row3 = result3->Next();
+           };
+ //      delete result3;
+         };
+         //
+         oss.str("");
+         oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
+             << thisid << " AND RUNHEADER_TIME="
+             << thisrht << " AND RUNTRAILER_TIME="
+             << thisrtt << " AND ID_RUN_FRAG!="
+             << thisid << " order by RUNHEADER_TIME asc;";
+         result3 = conn->Query(oss.str().c_str());
+         if ( result3 ){
+           row3 = result3->Next();
+           //
+           while ( row3 ){
+             //
+             // duplicated run
+             //
+             printf(" CHECK n.7 RUNs %u and %u HAVE SAME RUNTRAILER AND RUNHEADER (ARE THE SAME?) \n",thisid,(UInt_t)atoll(row3->GetField(0)));
+             row3 = result3->Next();
+           };
+  //     delete result3;
+         };
+       };
+       //
+       oss.str("");
+       oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
+           << thisid << " AND RUNHEADER_TIME>"
+           << thisrht << " AND RUNTRAILER_TIME<"
+           << thisrtt << " order by RUNHEADER_TIME asc;";
+       result3 = conn->Query(oss.str().c_str());
+       if ( result3 ){
+         row3 = result3->Next();
+         //
+         while ( row3 ){
+           //
+           // run contained in the checked one
+           //
+           printf(" CHECK n.6 RUN %u CONTAINS RUN %u  \n",thisid,(UInt_t)atoll(row3->GetField(0)));
+           row3 = result3->Next();
+         };
+ //      delete result3;
+       };
+       //
+     };
+     //
      //    if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){
      //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){
      if ( (thisrht < prevrtt) && (thisrht != prevrht) ){
-Line 4603 
 UInt_t PamelaDBOperations::Check(){
+Line 5388 
 UInt_t PamelaDBOperations::Check(){
        thisf = (TString)row2->GetField(0);
        if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
        test = 1;
+ //      delete result2;
      };
      //
      if ( (thisrtt < prevrht) && (thisrht != prevrht) ){
-Line 4613 
 UInt_t PamelaDBOperations::Check(){
+Line 5399 
 UInt_t PamelaDBOperations::Check(){
        oss.str("");
        oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
        result2 = conn->Query(oss.str().c_str());
-       if ( !result2 ) throw -4;;
+       if ( !result2 ) throw -4;
        row2 = result2->Next();
        prevf = (TString)row2->GetField(0);
        oss.str("");
-Line 4624 
 UInt_t PamelaDBOperations::Check(){
+Line 5410 
 UInt_t PamelaDBOperations::Check(){
        thisf = (TString)row2->GetField(0);
        if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
        test = 1;
+ //      delete result2;
      };
      //
      if ( (thisrht > thisrtt) && (thisrht != prevrht) ){
-Line 4645 
 UInt_t PamelaDBOperations::Check(){
+Line 5432 
 UInt_t PamelaDBOperations::Check(){
        thisf = (TString)row2->GetField(0);
        if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
        test = 1;
+ //      delete result2;
      };
+ ss:
      //
      prevrht = thisrht;
      prevrtt = thisrtt;
      previd = thisid;
      prevl0id = thisl0id;
      row = result->Next();
+ //    if ( result2 ) delete result2;
+ //    if ( result3 ) delete result3;
    };
    //
    return(test);

 Legend:



Removed from v.1.24
 


changed lines


 
Added in v.1.56
 Legend:



Removed from v.1.24
 


changed lines


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

	ViewVC Help
Powered by ViewVC 1.1.23