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

Diff of /YodaProfiler/src/PamelaDBOperations.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.8 by pam-fi, Mon Sep 11 16:38:18 2006 UTC revision 1.29 by mocchiut, Tue Apr 17 12:07:01 2007 UTC
# Line 14  Line 14 
14  #include <TSQLRow.h>  #include <TSQLRow.h>
15  #include <TTree.h>  #include <TTree.h>
16  #include <TGraph.h>  #include <TGraph.h>
17  #include <TDatime.h>  #include <TTimeStamp.h>
18  #include <TF1.h>  #include <TF1.h>
19  //  //
20  #include <EventHeader.h>  #include <EventHeader.h>
# Line 31  Line 31 
31  #include <varDump/VarDumpRecord.h>  #include <varDump/VarDumpRecord.h>
32  #include <physics/S4/S4Event.h>  #include <physics/S4/S4Event.h>
33  //  //
34    #include <sgp4.h>
35    
36  #include <PamelaDBOperations.h>  #include <PamelaDBOperations.h>
37  //  //
38  using namespace std;  using namespace std;
39  using namespace pamela;  using namespace pamela;
40    
41    // Some function to work with cTle stuff.
42    bool compTLE(cTle* tle1, cTle *tle2);
43    float getTleJulian(cTle *);
44    string getTleDatetime(cTle*);
45    
46  /**  /**
47   * Constructor.   * Constructor.
48   * @param host         hostname for the SQL connection.   * @param host         hostname for the SQL connection.
# Line 47  using namespace pamela; Line 54  using namespace pamela;
54   * @param obt0         file obt0.   * @param obt0         file obt0.
55   * @param tsync        file timesync.   * @param tsync        file timesync.
56   * @param debug        debug flag.   * @param debug        debug flag.
57     * @param tlefilename  ascii file with TLE 3 line elements.
58   */   */
59  PamelaDBOperations::PamelaDBOperations(TString host, TString user, TString password, TString filerawname, TString filerootname, UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t debug, Long64_t olderthan){  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){
60    //    //
61      STATIC=false;
62      if ( staticp ) STATIC=true;
63    //    //
64    SetConnection(host,user,password);    SetConnection(host,user,password);
65    //    //
# Line 62  PamelaDBOperations::PamelaDBOperations(T Line 72  PamelaDBOperations::PamelaDBOperations(T
72    SetTsync(tsync);    SetTsync(tsync);
73    SetObt0(obt0);    SetObt0(obt0);
74    //    //
75      SetTLEPath(tlefilename);
76      //
77    //    //
78    INSERT_RAW =!filerawname.IsNull();    INSERT_RAW =!filerawname.IsNull();
79    if(INSERT_RAW)SetRawName(filerawname);    if(INSERT_RAW)SetRawName(filerawname);
80    //    //
81    INSERT_ROOT = !filerootname.IsNull();    INSERT_ROOT = !filerootname.IsNull();
82    if(INSERT_ROOT){    if( INSERT_ROOT ){
83            file = TFile::Open(this->GetRootName().Data());      this->SetRootName(filerootname);
84            if ( !file ) INSERT_ROOT = false;      this->SetOrbitNo(dwinput);
85        file = TFile::Open(this->GetRootName().Data());
86      } else {
87        this->SetRootName("");
88    };    };
 //  this->OpenFile();  
89    //    //
90    this->SetID_RAW(0);    this->SetID_RAW(0);
91    this->SetID_ROOT(0);    this->SetID_ROOT(0);
92    
93    VALIDATE = false;    VALIDATE = false;
   clean_time = new TDatime();  
   if(olderthan >= 0){  
         VALIDATE = true;  
         UInt_t timelim = 0;  
         // timelim = now - olderthan  
         timelim =  (UInt_t)clean_time->Convert() - olderthan;  
         clean_time->Set(timelim,false);  
   };  
94        
95      //    //
96  };  };
97    
98  /**  /**
# Line 103  void PamelaDBOperations::Close(){ Line 109  void PamelaDBOperations::Close(){
109  // SETTERS  // SETTERS
110  //  //
111    
112    //
113    // must be out of the constructor in order to FORCE the validation of the latest runs in case you run the validation together with the latest file
114    //
115    void PamelaDBOperations::CheckValidate(Long64_t olderthan){
116      clean_time = new TDatime();
117      //
118      if(olderthan >= 0){
119        VALIDATE = true;
120        UInt_t timelim = 0;
121        timelim =  (UInt_t)clean_time->Convert(true) - olderthan;
122        clean_time->Set(timelim,false);
123      };
124    };
125    
126  /**  /**
127   * Open the DB connection   * Open the DB connection
128   * @param host         hostname for the SQL connection.   * @param host         hostname for the SQL connection.
# Line 110  void PamelaDBOperations::Close(){ Line 130  void PamelaDBOperations::Close(){
130   * @param password     password for the SQL connection.   * @param password     password for the SQL connection.
131   */   */
132  void PamelaDBOperations::SetConnection(TString host, TString user, TString password){  void PamelaDBOperations::SetConnection(TString host, TString user, TString password){
133      if ( IsDebug() ) printf(" Connecting using host = %s user = %s password = %s \n",host.Data(),user.Data(),password.Data());
134    conn = TSQLServer::Connect(host.Data(),user.Data(),password.Data());    conn = TSQLServer::Connect(host.Data(),user.Data(),password.Data());
135  };  };
136    
# Line 138  void PamelaDBOperations::SetDebugFlag(Bo Line 159  void PamelaDBOperations::SetDebugFlag(Bo
159  };  };
160    
161  /**  /**
162     * Set the autoboot flag
163     *
164     */
165    void PamelaDBOperations::SetAutoBoot(Bool_t dbg){
166      AUTOBOOT = dbg;
167    };
168    
169    /**
170     * Set the nofrag flag
171     *
172     */
173    void PamelaDBOperations::SetNoFrag(Bool_t nf){
174      NOFRAG = nf;
175    };
176    
177    /**
178   * Store the BOOT number of the RAW file.   * Store the BOOT number of the RAW file.
179   * @param boot    BOOT number of the RAW file   * @param boot    BOOT number of the RAW file
180   */   */
# Line 178  void PamelaDBOperations::SetRootName(TSt Line 215  void PamelaDBOperations::SetRootName(TSt
215  };  };
216    
217  /**  /**
218     * Store the downlink orbit number from filename.
219     */
220    void PamelaDBOperations::SetOrbitNo(UInt_t dwinput){
221      dworbit = 0;
222      //
223      if ( dwinput ){
224        dworbit = dwinput;
225        if ( IsDebug() ) printf(" Downlink orbit given by hand: %i  \n",dworbit);
226        return;
227      };
228      //
229      TString name = this->GetRootFile();
230      Int_t nlength = name.Length();
231      if ( nlength < 5 ) return;
232      TString dwo = 0;
233      for (Int_t i = 0; i<5; i++){
234        dwo.Append(name[i],1);
235      };
236      if ( dwo.IsDigit() ){
237        dworbit = (UInt_t)dwo.Atoi();
238      } else {
239        dwo="";
240        for (Int_t i = 8; i<13; i++){
241          dwo.Append(name[i],1);
242        };    
243        if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
244      };
245      if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());
246      return;
247    };
248    
249    
250    
251    /**
252   * Store the NOBOOT flag.   * Store the NOBOOT flag.
253   * @param noboot    true/false.   * @param noboot    true/false.
254   */   */
# Line 186  void PamelaDBOperations::SetNOBOOT(Bool_ Line 257  void PamelaDBOperations::SetNOBOOT(Bool_
257  };  };
258    
259  /**  /**
260     * Store path to the TLE file.
261     */
262    void PamelaDBOperations::SetTLEPath(TString str){
263      tlefilename = str;
264    };
265    
266    TString PamelaDBOperations::GetRawPath(){
267      if ( STATIC ){
268        return((TString)gSystem->DirName(filerawname.Data())+'/');
269      } else {
270        return((TString)gSystem->ExpandPathName("$PAM_RAW")+'/');
271      };
272    };
273    
274    TString PamelaDBOperations::GetRootPath(){
275      if ( STATIC ){
276        return((TString)gSystem->DirName(filerootname.Data())+'/');
277      } else {
278        return((TString)gSystem->ExpandPathName("$PAM_L0")+'/');
279      };
280    };
281    
282    /**
283   * Store the olderthan variable   * Store the olderthan variable
284   * @param olderthan   * @param olderthan
285   */   */
# Line 201  Bool_t PamelaDBOperations::SetID_RAW(){ Line 295  Bool_t PamelaDBOperations::SetID_RAW(){
295    TSQLResult *result = 0;    TSQLResult *result = 0;
296    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
297    oss.str("");    oss.str("");
298    oss << "SELECT ID FROM GL_RAW WHERE "    if ( STATIC ){
299        << " PATH = '" << this->GetRawPath().Data() << "' AND "      oss << "SELECT ID FROM GL_RAW WHERE "
300        << " NAME = '" << this->GetRawFile().Data() << "' ";          << " PATH = '" << this->GetRawPath().Data() << "' AND "
301            << " NAME = '" << this->GetRawFile().Data() << "' ";  
302      } else {
303        oss << "SELECT ID FROM GL_RAW WHERE "
304            << " PATH = '$PAM_RAW' AND "
305            << " NAME = '" << this->GetRawFile().Data() << "' ";  
306      }
307    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
308    if ( result == NULL ) throw -4;    if ( result == NULL ) throw -4;
309    row = result->Next();    row = result->Next();
# Line 234  Int_t PamelaDBOperations::SetUpperLimits Line 334  Int_t PamelaDBOperations::SetUpperLimits
334    UInt_t nevent = 0;    UInt_t nevent = 0;
335    UInt_t pktlast = 0;    UInt_t pktlast = 0;
336    UInt_t obtlast = 0;    UInt_t obtlast = 0;
337    UInt_t t_pktlast = 0;    Long64_t t_pktlast = 0LL;
338    UInt_t t_obtlast = 0;    //  UInt_t t_obtlast = 0;
339    UInt_t upperpkt2 = 0;    Long64_t t_obtlast = 0LL;
340    ULong64_t upperobt2 = 0;    Long64_t upperpkt2 = 0LL;
341      Long64_t upperobt2 = 0LL;
342    UInt_t zomp = 0;    UInt_t zomp = 0;
343    UInt_t jump = 50000; // was 5000    UInt_t jump = 50000; // was 5000
344    EventCounter *code=0;    EventCounter *code=0;
345    //    //
346    //   pcksList packetsNames;    Long64_t deltapkt = 5000LL;
347    //   pcksList::iterator Iter;    Long64_t deltaobt = 50000LL;
348    //   getPacketsNames(packetsNames);    //
349      pcksList packetsNames;
350      pcksList::iterator Iter;
351      getPacketsNames(packetsNames);
352    //    //
353    pktfirst = 0;    pktfirst = 0;
354    obtfirst = 0;    obtfirst = 0;
355      ppktfirst = 0;
356      pobtfirst = 0;
357    //    //
358    TTree *T = 0;    TTree *T = 0;
359    T = (TTree*)file->Get("Physics");    T = (TTree*)file->Get("Physics");
# Line 261  Int_t PamelaDBOperations::SetUpperLimits Line 367  Int_t PamelaDBOperations::SetUpperLimits
367    ph = eh->GetPscuHeader();    ph = eh->GetPscuHeader();
368    pktfirst = ph->GetCounter();    pktfirst = ph->GetCounter();
369    obtfirst = ph->GetOrbitalTime();      obtfirst = ph->GetOrbitalTime();  
370      ppktfirst = pktfirst;
371      pobtfirst = obtfirst;
372    //    //
373    //   code = eh->GetCounter();    code = eh->GetCounter();
374    //   UInt_t en = 0;    UInt_t en = 0;
375    //   for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){    for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){
376    //     en = code->Get(GetPacketType(*Iter));      en = code->Get(GetPacketType(*Iter));
377    //     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) ){
378    //};          if ( IsDebug() ) printf(" We have a calibration before the first physic packet: packet type is %s, entries: %i \n",*Iter,en);
379          //
380          TTree *TC = 0;
381          TC = (TTree*)file->Get("CalibCalPed");
382          if ( !TC || TC->IsZombie() ) throw -16;
383          EventHeader *ehc = 0;
384          PscuHeader *phc = 0;
385          TC->SetBranchAddress("Header", &ehc);
386          TC->GetEntry(0);
387          phc = ehc->GetPscuHeader();
388          pktfirst = phc->GetCounter();
389          obtfirst = phc->GetOrbitalTime();  
390          //
391        };
392      };  
393    //    //
394    T->GetEntry(nevent-1);    T->GetEntry(nevent-1);
395    ph = eh->GetPscuHeader();    ph = eh->GetPscuHeader();
# Line 278  Int_t PamelaDBOperations::SetUpperLimits Line 400  Int_t PamelaDBOperations::SetUpperLimits
400    upperobt = OBT(obtlast);    upperobt = OBT(obtlast);
401    upperentry = nevent-1;    upperentry = nevent-1;
402    //    //
403    if ( IsDebug() ) printf(" First entries are: OBT %llu pkt_num %i \n",obtfirst,pktfirst);    if ( IsDebug() ) printf(" First entries are: OBT %i pkt_num %i \n",obtfirst,pktfirst);
404    //    //
405    if ( IsDebug() ) printf(" Last entries are: OBT %llu pkt_num %i entry %i\n",upperobt,upperpkt,upperentry);    if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld entry %i\n",upperobt,upperpkt,upperentry);
406    //    //
407    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(1);
408    //    //
409    if ( !nevent ) return(2);    if ( !nevent ) return(2);
410    //    //
411    if ( nevent < 2 ) return(4);    if ( nevent < 2 ) return(4);
412      if ( nevent < jump ) jump = 1;
413      //  if ( nevent < jump ) jump = int(nevent/10);
414      //  if ( !jump ) jump = 1;
415    //    //
416    if ( PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst) ){    if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) && nevent > deltapkt ){
417        //
418        if ( IsDebug() ) printf(" starting jump %i \n",jump);
419      // go back      // go back
420      zomp = nevent - 2;      zomp = nevent - 2;
421      //      //
# Line 310  Int_t PamelaDBOperations::SetUpperLimits Line 437  Int_t PamelaDBOperations::SetUpperLimits
437          upperpkt2 = PKT(ph->GetCounter());          upperpkt2 = PKT(ph->GetCounter());
438          upperobt2 = OBT(ph->GetOrbitalTime());          upperobt2 = OBT(ph->GetOrbitalTime());
439          //          //
440          if ( (t_pktlast < upperpkt && t_obtlast > upperobt) || (t_pktlast < upperpkt2 && t_obtlast > upperobt2) ) throw -13;              if ( (t_pktlast < upperpkt && t_obtlast > upperobt) || (t_pktlast < upperpkt2 && t_obtlast > upperobt2) ){
441              if ( IsDebug() ) printf(" .-. upperpkt2 %lld upperobt2 %lld \n",upperpkt2,upperobt2);  
442              if ( IsDebug() ) printf(" .-. upperpkt %lld t_pktlast %lld upperobt %lld t_obtlast %lld \n",upperpkt,t_pktlast,upperobt,t_obtlast);    
443              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);    
444              throw -13;    
445            };
446          //          //
447          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){
448            zomp = i + jump + 1;            zomp = i + jump + 1;
449            if ( zomp > nevent-2 ) zomp = nevent - 2;            if ( zomp > nevent-2 ) zomp = nevent - 2;
450            if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %i pktlast %i upperobt %llu obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);            if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
451            break;            break;
452          };          };
453          //          //
# Line 350  Int_t PamelaDBOperations::SetUpperLimits Line 482  Int_t PamelaDBOperations::SetUpperLimits
482    //    //
483    rhev = rh->GetEntries();    rhev = rh->GetEntries();
484    rtev = rt->GetEntries();    rtev = rt->GetEntries();
485    UInt_t sobtt = 0;    Long64_t sobtt = 0LL;
486    UInt_t sobth = 0;    Long64_t sobth = 0LL;
487    UInt_t spktt = 0;    Long64_t spktt = 0LL;
488    UInt_t spkth = 0;    Long64_t spkth = 0LL;
489    UInt_t pktt = 0;    Long64_t pktt = 0LL;
490    ULong64_t obtt = 0;    Long64_t obtt = 0LL;
491    UInt_t pkth = 0;    Long64_t pkth = 0LL;
492    ULong64_t obth = 0;    Long64_t obth = 0LL;
   //  
   T->GetEntry(upperentry);  
   code = eh->GetCounter();  
   Int_t lasttrail = code->Get(pctp->RunTrailer);  
   Int_t lasthead = code->Get(pctp->RunHeader);  
   if ( lasttrail < rtev ){  
     rt->GetEntry(lasttrail);  
     pht = eht->GetPscuHeader();  
     pktt = PKT(pht->GetCounter());  
     obtt = OBT(pht->GetOrbitalTime());  
   };  
   //  
   if ( lasthead < rhev ){  
     rh->GetEntry(lasthead);  
     phh = ehh->GetPscuHeader();  
     pkth = PKT(phh->GetCounter());  
     obth = OBT(phh->GetOrbitalTime());  
   };  
   //  
   if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);  
   if ( pkth > upperpkt && obth > upperobt ){  
     if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt);  
     upperpkt = pkth;  
     upperobt = obth;  
     rhev = lasthead+1;  
   } else {  
     rhev = lasthead;  
   };  
   if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);  
   //  
   if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);  
   if ( pktt > upperpkt && obtt > upperobt ){  
     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt);  
     upperpkt = pktt;  
     upperobt = obtt;  
     rtev = lasttrail+1;  
   } else {  
     rtev = lasttrail;  
   };  
   if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);  
   //  goto kikko;  
   //  
   //  
   // Check if runtrailer/runheader are within lower limits  
   //  
493    //    //
494    pkth = 0;    if ( rhev || rtev ){
495    obth = 0;  
496    spkth = 0;      T->GetEntry(upperentry);
497    sobth = 0;      code = eh->GetCounter();
498    for (Int_t k=0; k<rhev; k++){      Int_t lasttrail = code->Get(pctp->RunTrailer);
499      if ( k > 0 ){      Int_t lasthead = code->Get(pctp->RunHeader);
500        spkth = pkth;      if ( lasttrail < rtev ){
501        sobth = obth;        rt->GetEntry(lasttrail);
502      };        pht = eht->GetPscuHeader();
503      rh->GetEntry(k);        pktt = PKT(pht->GetCounter());
504      phh = ehh->GetPscuHeader();        obtt = OBT(pht->GetOrbitalTime());
505      pkth = PKT(phh->GetCounter());      };
506      obth = OBT(phh->GetOrbitalTime());      //
507      //      if ( lasthead < rhev ){
508  //    if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth);        rh->GetEntry(lasthead);
509      //        phh = ehh->GetPscuHeader();
510      if ( pkth < spkth && obth < sobth ){        pkth = PKT(phh->GetCounter());
511        if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);        obth = OBT(phh->GetOrbitalTime());
512        //      };
513        rhev = k-1;      //
514        rh->GetEntry(rhev);      if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
515        pkth = spkth;      if ( pkth > upperpkt && obth > upperobt ){
516        obth = sobth;        if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
517        //        upperpkt = pkth;
518        UInt_t evbefh = 0;        upperobt = obth;
519        code = ehh->GetCounter();        rhev = lasthead+1;
520        evbefh = code->Get(pctp->Physics);          } else {
521        if ( evbefh >= 0 ){        rhev = lasthead;
522          T->GetEntry(evbefh);      };
523          ph = eh->GetPscuHeader();      if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
524          t_pktlast = PKT(ph->GetCounter());      //
525          t_obtlast = OBT(ph->GetOrbitalTime());      if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
526          if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump      if ( pktt > upperpkt && obtt > upperobt ){
527            upperpkt = pkth;        if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
528            upperobt = obth;        upperpkt = pktt;
529            upperentry = evbefh-1;        upperobt = obtt;
530          } else {        rtev = lasttrail+1;
531            while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){      } else {
532              evbefh++;          rtev = lasttrail;
533              T->GetEntry(evbefh);      };
534              ph = eh->GetPscuHeader();      if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
535              t_pktlast = PKT(ph->GetCounter());      //  goto kikko;
536              t_obtlast = OBT(ph->GetOrbitalTime());      //
537            };      //
538            T->GetEntry(evbefh-1);      // Check if runtrailer/runheader are within lower limits
539        //
540        //
541        pkth = 0LL;
542        obth = 0LL;
543        spkth = 0LL;
544        sobth = 0LL;
545        for (Int_t k=0; k<rhev; k++){
546          if ( k > 0 ){
547            spkth = pkth;
548            sobth = obth;
549          };
550          rh->GetEntry(k);
551          phh = ehh->GetPscuHeader();
552          pkth = PKT(phh->GetCounter());
553          obth = OBT(phh->GetOrbitalTime());
554          //
555          //    if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth);
556          //
557          if ( pkth < spkth && obth < sobth ){
558            if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);
559            //
560            rhev = k-1;
561            rh->GetEntry(rhev);
562            pkth = spkth;
563            obth = sobth;
564            //
565            UInt_t evbefh = 0;
566            code = ehh->GetCounter();
567            evbefh = code->Get(pctp->Physics);    
568            if ( evbefh >= 0 ){
569              T->GetEntry(evbefh);
570            ph = eh->GetPscuHeader();            ph = eh->GetPscuHeader();
571            upperpkt = PKT(ph->GetCounter());            t_pktlast = PKT(ph->GetCounter());
572            upperobt = OBT(ph->GetOrbitalTime());            t_obtlast = OBT(ph->GetOrbitalTime());
573            upperentry = evbefh-1;            if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump
574          };                  upperpkt = pkth;
575                upperobt = obth;
576                upperentry = evbefh-1;
577              } else {
578                while ( t_pktlast > spkth && t_obtlast > sobth && evbefh < nevent ){
579                  evbefh++;
580                  T->GetEntry(evbefh);
581                  ph = eh->GetPscuHeader();
582                  t_pktlast = PKT(ph->GetCounter());
583                  t_obtlast = OBT(ph->GetOrbitalTime());
584                };
585                T->GetEntry(evbefh-1);
586                ph = eh->GetPscuHeader();
587                upperpkt = PKT(ph->GetCounter());
588                upperobt = OBT(ph->GetOrbitalTime());
589                upperentry = evbefh-1;
590              };    
591            };
592            if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
593            goto kikko0;
594        };        };
       if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);  
       goto kikko0;  
     };  
   };  
  kikko0:  
   //  
   //  
   //  
   pktt = 0;  
   obtt = 0;  
   spktt = 0;  
   sobtt = 0;  
   for (Int_t k=0; k<rtev; k++){  
     if ( k > 0 ){  
       spktt = pktt;  
       sobtt = obtt;  
595      };      };
596      rt->GetEntry(k);    kikko0:
     pht = eht->GetPscuHeader();  
     pktt = PKT(pht->GetCounter());  
     obtt = OBT(pht->GetOrbitalTime());  
597      //      //
 //    if ( IsDebug() ) printf(" k %i rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt);  
598      //      //
599      if ( pktt < spktt && obtt < sobtt ){      //
600        if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);      pktt = 0LL;
601        obtt = 0LL;
602        spktt = 0LL;
603        sobtt = 0LL;
604        for (Int_t k=0; k<rtev; k++){
605          if ( k > 0 ){
606            spktt = pktt;
607            sobtt = obtt;
608          };
609          rt->GetEntry(k);
610          pht = eht->GetPscuHeader();
611          pktt = PKT(pht->GetCounter());
612          obtt = OBT(pht->GetOrbitalTime());
613        //        //
614        rtev = k-1;        //    if ( IsDebug() ) printf(" k %i rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt);
       rt->GetEntry(rtev);  
       pktt = spktt;  
       obtt = sobtt;  
       if ( IsDebug() ) printf(" lasttrail %i pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);  
615        //        //
616        UInt_t evbeft = 0;        if ( pktt < spktt && obtt < sobtt ){
617        code = eht->GetCounter();          if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);
618        evbeft = code->Get(pctp->Physics);              //
619        if ( evbeft >= 0 ){          rtev = k-1;
620          T->GetEntry(evbeft);          rt->GetEntry(rtev);
621          ph = eh->GetPscuHeader();          pktt = spktt;
622          t_pktlast = PKT(ph->GetCounter());          obtt = sobtt;
623          t_obtlast = OBT(ph->GetOrbitalTime());          if ( IsDebug() ) printf(" lasttrail %i pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
624          if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump          //
625            upperpkt = pktt;          UInt_t evbeft = 0;
626            upperobt = obtt;          code = eht->GetCounter();
627            upperentry = evbeft-1;          evbeft = code->Get(pctp->Physics);    
628          } else {          if ( evbeft >= 0 ){
629            while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){            T->GetEntry(evbeft);
630              evbeft++;              ph = eh->GetPscuHeader();
631              T->GetEntry(evbeft);            t_pktlast = PKT(ph->GetCounter());
632              t_obtlast = OBT(ph->GetOrbitalTime());
633              if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump
634                upperpkt = pktt;
635                upperobt = obtt;
636                upperentry = evbeft-1;
637              } else {
638                while ( t_pktlast > spktt && t_obtlast > sobtt && evbeft < nevent ){
639                  evbeft++;
640                  T->GetEntry(evbeft);
641                  ph = eh->GetPscuHeader();
642                  t_pktlast = PKT(ph->GetCounter());
643                  t_obtlast = OBT(ph->GetOrbitalTime());
644                };
645                T->GetEntry(evbeft-1);
646              ph = eh->GetPscuHeader();              ph = eh->GetPscuHeader();
647              t_pktlast = PKT(ph->GetCounter());              upperpkt = PKT(ph->GetCounter());
648              t_obtlast = OBT(ph->GetOrbitalTime());              upperobt = OBT(ph->GetOrbitalTime());
649                upperentry = evbeft-1;
650            };            };
           T->GetEntry(evbeft-1);  
           ph = eh->GetPscuHeader();  
           upperpkt = PKT(ph->GetCounter());  
           upperobt = OBT(ph->GetOrbitalTime());  
           upperentry = evbeft-1;  
651          };          };
652            if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
653            goto kikko;
654            //      break;
655            //
656        };        };
657        if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);        //  
       goto kikko;  
       //      break;  
       //  
658      };      };
659      //        //
660      kikko:
661        //
662        T->GetEntry(upperentry);
663        code = eh->GetCounter();
664        lasttrail = code->Get(pctp->RunTrailer);
665        lasthead = code->Get(pctp->RunHeader);
666        if ( lasttrail < rtev ){
667          rt->GetEntry(lasttrail);
668          pht = eht->GetPscuHeader();
669          pktt = PKT(pht->GetCounter());
670          obtt = OBT(pht->GetOrbitalTime());
671        };
672        //
673        if ( lasthead < rhev ){
674          rh->GetEntry(lasthead);
675          phh = ehh->GetPscuHeader();
676          pkth = PKT(phh->GetCounter());
677          obth = OBT(phh->GetOrbitalTime());
678        };
679        //
680        if ( IsDebug() ) printf(" rhev before %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
681        if ( pkth > upperpkt && obth > upperobt ){
682          if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %lld upperp %lld oh %lld uppero %lld \n",pkth,upperpkt,obth,upperobt);
683          upperpkt = pkth;
684          upperobt = obth;
685          rhev = lasthead+1;
686        } else {
687          rhev = lasthead;
688        };
689        if ( IsDebug() ) printf(" rhev after %i ph %lld upperp %lld oh %lld uppero %lld \n",rhev,pkth,upperpkt,obth,upperobt);
690        //
691        if ( IsDebug() ) printf(" rtev beforev %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
692        if ( pktt > upperpkt && obtt > upperobt ){
693          if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %lld upperp %lld ot %lld uppero %lld \n",pktt,upperpkt,obtt,upperobt);
694          upperpkt = pktt;
695          upperobt = obtt;
696          rtev = lasttrail+1;
697        } else {
698          rtev = lasttrail;
699        };
700        if ( IsDebug() ) printf(" rtev after %i  pt %lld upperp %lld ot %lld uppero %lld \n",rtev,pktt,upperpkt,obtt,upperobt);
701        //
702    };    };
703    //    //
704   kikko:    if ( IsDebug() ) printf(" Upper limits are: OBT %lld pkt_num %lld upper entry %i \n",upperobt,upperpkt,upperentry);
705    //    //
706    T->GetEntry(upperentry);    return(0);
707    code = eh->GetCounter();  }
708    lasttrail = code->Get(pctp->RunTrailer);  
709    lasthead = code->Get(pctp->RunHeader);  /**
710    if ( lasttrail < rtev ){   *
711      rt->GetEntry(lasttrail);   * Trick to have unique RUN ID numbers even when runs are deleted and mysql deamon restarted.
712      pht = eht->GetPscuHeader();   * Entries in the _RUNID_GEN table are never deleted.
713      pktt = PKT(pht->GetCounter());   *
714      obtt = OBT(pht->GetOrbitalTime());   **/
715    };  UInt_t PamelaDBOperations::AssignRunID(){
716    //    //
717    if ( lasthead < rhev ){    TSQLResult *result = 0;
718      rh->GetEntry(lasthead);    TSQLRow *row = 0;
719      phh = ehh->GetPscuHeader();    UInt_t runid = 0;
     pkth = PKT(phh->GetCounter());  
     obth = OBT(phh->GetOrbitalTime());  
   };  
720    //    //
721    if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);    stringstream   oss;
722    if ( pkth > upperpkt && obth > upperobt ){    //  
723      if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt);    oss.str("");
724      upperpkt = pkth;    oss << "INSERT INTO _RUNID_GEN VALUES (NULL);";
725      upperobt = obth;    result = conn->Query(oss.str().c_str());
726      rhev = lasthead+1;    if ( !result ) throw -10;
727    } else {    oss.str("");
728      rhev = lasthead;    oss << "SELECT ID FROM _RUNID_GEN ORDER BY ID DESC LIMIT 1;";
729    };    result = conn->Query(oss.str().c_str());
730    if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);    if ( !result ) throw -10;
731    //    //
732    if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);    row = result->Next();
   if ( pktt > upperpkt && obtt > upperobt ){  
     if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt);  
     upperpkt = pktt;  
     upperobt = obtt;  
     rtev = lasttrail+1;  
   } else {  
     rtev = lasttrail;  
   };  
   if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);  
733    //    //
734    if ( IsDebug() ) printf(" Upper limits are: OBT %llu pkt_num %i upper entry %i \n",upperobt,upperpkt,upperentry);    if ( !row ) throw -28;
735    //    //
736    return(0);    runid = (UInt_t)atoll(row->GetField(0));
737  }    //
738      return(runid);
739    };
740    
741  //  //
742  // GETTERS  // GETTERS
# Line 651  const PacketType* PamelaDBOperations::Ge Line 819  const PacketType* PamelaDBOperations::Ge
819  // void PamelaDBOperations::OpenFile(){  // void PamelaDBOperations::OpenFile(){
820  //   file = TFile::Open(this->GetRootName().Data());  //   file = TFile::Open(this->GetRootName().Data());
821  //   //  //   //
822  //   if ( !file ) throw -12;  
823  // };  void PamelaDBOperations::CheckFile(){  
824      if ( !file ) throw -12;
825    };
826    
827    
828  /**  /**
# Line 665  void PamelaDBOperations::CheckConnection Line 835  void PamelaDBOperations::CheckConnection
835    if( !conn ) throw -1;    if( !conn ) throw -1;
836    bool connect = conn->IsConnected();    bool connect = conn->IsConnected();
837    if( !connect ) throw -1;    if( !connect ) throw -1;
838      //
839      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());
840      //
841      if ( !dworbit && strcmp(this->GetRootName().Data(),"") ) throw -27;
842      //
843      // set DB timezone to UTC
844      //
845      stringstream oss;
846      //  
847      oss.str("");
848      oss << "SET time_zone='+0:00';";
849      TSQLResult *result = 0;
850      result = conn->Query(oss.str().c_str());
851      if ( !result ) throw -10;
852      oss.str("");
853      oss << "SET wait_timeout=173000;";
854      conn->Query(oss.str().c_str());
855      //
856  };  };
857    
858  /**  /**
859   * Return the correct packet number if we went back to zero   * Return the correct packet number if we went back to zero
860   */   */
861  UInt_t PamelaDBOperations::PKT(UInt_t pkt_num){    Long64_t PamelaDBOperations::PKT(UInt_t pkt_num){  
862    //    //
863    //  if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u (pktfirst - (UInt_t)(16777212/2)) is %u  \n",pkt_num,pktfirst,(pktfirst - (UInt_t)(16777212/2)));    if ( IsDebug() ) printf(" pkt conversion: pkt_num is %u pktfirst is %u  (UInt_t)(16777214/2)) is %u  \n",pkt_num,ppktfirst,(UInt_t)(16777214/2));
864    //    //
865    if ( pkt_num < (pktfirst/2) && pktfirst > (16777214/2)  ) return((pkt_num+16777215));    if ( pkt_num < (ppktfirst/2) && ppktfirst > (16777214/2)  ){
866        if ( IsDebug() ) printf(" rise up pktnum %lld \n",(Long64_t)pkt_num+16777215LL);
867        return((Long64_t)pkt_num+16777215LL);
868      };
869    //    //
870    if ( pkt_num > pktfirst*2 && pkt_num > (16777214/2) ){    if ( pkt_num > ppktfirst*2 && pkt_num > (16777214/2) ){
871      if ( (pkt_num-16777215) < 0 ){      if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL);
872        return((16777215-pkt_num));      return((Long64_t)pkt_num-16777215LL);
     } else {  
       return((pkt_num-16777215));  
     };  
873    };    };
874    //    //
875    return(pkt_num);    if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)pkt_num);
876      return((Long64_t)pkt_num);
877    //    //
878  };  };
879    
880  /**  /**
881   * Return the correct On Board Time if we went back to zero   * Return the correct On Board Time if we went back to zero
882   */   */
883  ULong64_t PamelaDBOperations::OBT(UInt_t obt){    Long64_t PamelaDBOperations::OBT(UInt_t obt){  
884    //    //
885    if ( obt < (obtfirst/2) && obtfirst > (numeric_limits<UInt_t>::max()/2) ) return((ULong64_t)(obt+numeric_limits<UInt_t>::max()));    if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ) return((Long64_t)(obt+numeric_limits<UInt_t>::max()));
886    //    //
887    if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){    if ( obt > (pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
888      if ( (obt-numeric_limits<UInt_t>::max()) < 0  ){      return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
       return((ULong64_t)(numeric_limits<UInt_t>::max()-obt));  
     } else {  
       return((ULong64_t)(obt-numeric_limits<UInt_t>::max()));  
     };  
889    };    };
890    //    //
891    return((ULong64_t)obt);    return((Long64_t)obt);
892  };  };
893    
894  /**  /**
# Line 818  Int_t PamelaDBOperations::insertPamelaRa Line 1003  Int_t PamelaDBOperations::insertPamelaRa
1003    if ( idr ) return(1);    if ( idr ) return(1);
1004    //    //
1005    oss.str("");    oss.str("");
1006    oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"    if ( STATIC ){
1007        << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";      oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
1008            << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";
1009      } else {
1010        oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')";
1011      };
1012    if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;    if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
1013    //    //
1014    idr =  this->SetID_RAW();    idr =  this->SetID_RAW();
# Line 844  Int_t PamelaDBOperations::insertPamelaGL Line 1033  Int_t PamelaDBOperations::insertPamelaGL
1033    if ( this->GetID_RAW() == 0 )  throw -11;    if ( this->GetID_RAW() == 0 )  throw -11;
1034    //    //
1035    oss.str("");    oss.str("");
1036    oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ID_RAW<= "    oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='"
1037        << id << " AND TO_ID_RAW >= "        << this->GetRawFile().Data() << "';";
       << id << ";";  
1038    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
1039    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1040    if ( !result ) throw -10;    if ( !result ) throw -10;
1041    row = result->Next();    row = result->Next();
   if ( !row ) throw -10;  
1042    //    //
1043    t0 = (UInt_t)TDatime(row->GetField(0)).Convert();    if ( !row ){
1044        oss.str("");
1045        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< "
1046            << dworbit << " order by FROM_ORBIT desc limit 1;";
1047        if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
1048        result = conn->Query(oss.str().c_str());
1049        if ( !result ) throw -10;
1050        row = result->Next();
1051        if ( !row ) throw -10;
1052      };
1053      //
1054      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);
1055      t0 = (UInt_t)tu.GetSec();
1056      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));
1057      //
1058    /*    /*
1059     * Verify that the TIMESYNC have been not already processed     * Verify that the TIMESYNC have been not already processed
1060     */     */
# Line 868  Int_t PamelaDBOperations::insertPamelaGL Line 1069  Int_t PamelaDBOperations::insertPamelaGL
1069    if (result == NULL) throw -10;    if (result == NULL) throw -10;
1070    row = result->Next();    row = result->Next();
1071    if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){    if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){
1072        if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
1073      toffset = (UInt_t)atoll(row->GetField(2)) - (UInt_t)(this->OBT((UInt_t)atoll(row->GetField(1)))/1000) + t0;      toffset = (UInt_t)atoll(row->GetField(2)) - (UInt_t)(this->OBT((UInt_t)atoll(row->GetField(1)))/1000) + t0;
1074        //
1075        tsync = (UInt_t)atoll(row->GetField(2));
1076        obt0 = (UInt_t)atoll(row->GetField(1));
1077        //
1078      return(1);      return(1);
1079    };    };
1080    //    //
# Line 925  Int_t PamelaDBOperations::insertPamelaGL Line 1131  Int_t PamelaDBOperations::insertPamelaGL
1131          //          //
1132          TYPE = 55;//224;          TYPE = 55;//224;
1133          //          //
1134            if ( IsDebug() ) printf("mcmd tsync %i tsync %u obt %u \n",i,TSYNC,OBT);
1135            //
1136          if ( TSYNC && OBT ){          if ( TSYNC && OBT ){
1137            existsts = true;            existsts = true;
1138            goto eout;            goto eout;
# Line 947  Int_t PamelaDBOperations::insertPamelaGL Line 1155  Int_t PamelaDBOperations::insertPamelaGL
1155      //      //
1156      rt->SetBranchAddress("RunTrailer", &runt);      rt->SetBranchAddress("RunTrailer", &runt);
1157      //      //
1158      if ( rhev > 0 ){      Int_t nrhev = rh->GetEntries();
1159        rh->GetEntry(0);      Int_t nrtev = rt->GetEntries();
1160        //      if ( IsDebug() ) printf(" ou nevent %i rhev %i rtev %i \n",nevent,nrhev,nrtev);
1161        TSYNC = runh->LAST_TIME_SYNC_INFO;      //
1162        OBT = runh->OBT_TIME_SYNC * 1000;      if ( nrhev > 0 ){
1163        //        for (Int_t i=0; i<nrhev; i++){
1164        TYPE = 20;          //
1165        //          rh->GetEntry(i);
1166        if ( TSYNC && OBT ){          //
1167          existsts = true;          TSYNC = runh->LAST_TIME_SYNC_INFO;
1168          goto eout;          OBT = runh->OBT_TIME_SYNC * 1000;
1169            //
1170            TYPE = 20;
1171            //
1172            if ( IsDebug() ) printf("runheader %i tsync %u obt %u \n",i,TSYNC,OBT);
1173            //
1174            if ( TSYNC && OBT ){
1175              existsts = true;
1176              goto eout;
1177            };
1178        };        };
1179        //        //
1180      };      };
1181      if ( rtev > 0 ){      if ( nrtev > 0 ){
1182        //        //
1183        if ( IsDebug() ) printf(" No runheader \n");        if ( IsDebug() ) printf(" No runheader \n");
1184        signal = 6;        signal = 6;
1185        //        //
1186        rt->GetEntry(0);        for (Int_t i=0; i<nrtev; i++){
1187        //          //
1188        TSYNC = runt->LAST_TYME_SYNC_INFO;          rt->GetEntry(i);
1189        OBT = runt->OBT_TYME_SYNC * 1000;          //
1190        //          TSYNC = runt->LAST_TYME_SYNC_INFO;
1191        TYPE = 21;          OBT = runt->OBT_TYME_SYNC * 1000;
1192        //          //
1193        if ( TSYNC && OBT ){          TYPE = 21;
1194          existsts = true;          //
1195          goto eout;          if ( IsDebug() ) printf("runtrailer %i tsync %u obt %u \n",i,TSYNC,OBT);
1196            //
1197            if ( TSYNC && OBT ){
1198              existsts = true;
1199              goto eout;
1200            };
1201        };        };
1202        //        //
1203      } else {      } else {
# Line 1046  Int_t PamelaDBOperations::insertPamelaGL Line 1268  Int_t PamelaDBOperations::insertPamelaGL
1268    conn->Query(oss.str().c_str());    conn->Query(oss.str().c_str());
1269    if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to fill it:\n %s \n",oss.str().c_str());
1270    //    //
1271      if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
1272    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;
1273    //    //
1274      tsync = TSYNC;
1275      obt0 = OBT;
1276      //
1277    delete result;    delete result;
1278    return(signal);    return(signal);
1279  }  }
# Line 1063  Int_t PamelaDBOperations::insertPamelaRo Line 1289  Int_t PamelaDBOperations::insertPamelaRo
1289    UInt_t idtimesync = 0;    UInt_t idtimesync = 0;
1290    //    //
1291    oss.str("");    oss.str("");
1292    oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "    if ( STATIC ){
1293        << " LEFT JOIN GL_ROOT "      oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
1294        << " ON GL_RAW.ID = GL_ROOT.ID_RAW "          << " LEFT JOIN GL_ROOT "
1295        << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "          << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
1296        << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";          << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "
1297            << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
1298      } else {
1299        oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
1300            << " LEFT JOIN GL_ROOT "
1301            << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
1302            << " WHERE GL_RAW.PATH = '$PAM_RAW' AND "
1303            << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
1304      };
1305    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1306    //    //
1307    if ( !result ) throw -12;    if ( !result ) throw -12;
# Line 1094  Int_t PamelaDBOperations::insertPamelaRo Line 1328  Int_t PamelaDBOperations::insertPamelaRo
1328    idtimesync = (UInt_t)atoll(row->GetField(0));    idtimesync = (UInt_t)atoll(row->GetField(0));
1329    //    //
1330    oss.str("");    oss.str("");
1331    oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"    if ( STATIC ){
1332        << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";      oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
1333            << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";
1334      } else {
1335        oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
1336            << this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')";
1337      };
1338    //    //
1339    if (conn->Query(oss.str().c_str()) == 0) throw -4;    if (conn->Query(oss.str().c_str()) == 0) throw -4;
1340    //    //
# Line 1122  Int_t PamelaDBOperations::assignBOOT_NUM Line 1361  Int_t PamelaDBOperations::assignBOOT_NUM
1361    TSQLResult *result = 0;    TSQLResult *result = 0;
1362    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
1363    oss.str("");    oss.str("");
1364    oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "    if ( STATIC ){
1365        << " PATH = '" << this->GetRawPath().Data() << "' AND "      oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
1366        << " NAME = '" << this->GetRawFile().Data() << "' ";          << " PATH = '" << this->GetRawPath().Data() << "' AND "
1367            << " NAME = '" << this->GetRawFile().Data() << "' ";
1368      } else {
1369        oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
1370            << " PATH = '$PAM_RAW' AND "
1371            << " NAME = '" << this->GetRawFile().Data() << "' ";
1372      };
1373    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1374    //    //
1375    if ( !result ) throw -4;;    if ( !result ) throw -4;;
# Line 1147  Int_t PamelaDBOperations::assignBOOT_NUM Line 1392  Int_t PamelaDBOperations::assignBOOT_NUM
1392    VarDumpEvent  *vde = 0;    VarDumpEvent  *vde = 0;
1393    VarDumpRecord *vdr = 0;    VarDumpRecord *vdr = 0;
1394    //    //
1395      Bool_t found = false;
1396    trDumpEv->SetBranchAddress("VarDump", &vde);    trDumpEv->SetBranchAddress("VarDump", &vde);
1397    if ( trDumpEv->GetEntries() > 0 ){    if ( trDumpEv->GetEntries() > 0 ){
1398      Bool_t found = false;      found = false;
1399      for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){      for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){
1400        trDumpEv->GetEntry(i);        trDumpEv->GetEntry(i);
1401        vde->Records->GetEntries();        //      vde->Records->GetEntries();
1402        if ( vde->Records->GetEntries()>0 ){        if ( vde->Records->GetEntries()>5 ){
1403          found = true;          found = true;
1404          goto fill;          goto fill;
1405        };        };
# Line 1166  Int_t PamelaDBOperations::assignBOOT_NUM Line 1412  Int_t PamelaDBOperations::assignBOOT_NUM
1412        this->SetBOOTnumber((Int_t)vdr->VAR_VALUE);        this->SetBOOTnumber((Int_t)vdr->VAR_VALUE);
1413        //        //
1414      } else {      } else {
1415        if ( !this->GetBOOTnumber() ) return(4);        if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4);
1416        };
1417      } else {
1418        if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2);
1419      };
1420      //
1421      UInt_t bn = 0;
1422      Bool_t afound = false;
1423      if ( !found && this->AutoBoot()){
1424        afound = true;
1425        //
1426        // Search for other files with similar timesync
1427        //
1428        if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0);
1429        UInt_t upperts = tsync-(obt0/1000)+5;
1430        UInt_t lowerts = tsync-(obt0/1000)-5;
1431        oss.str("");
1432        oss << "select GL_RAW.BOOT_NUMBER from GL_TIMESYNC LEFT JOIN GL_RAW ON GL_RAW.ID = GL_TIMESYNC.ID_RAW where TIMESYNC-(OBT0/1000)<"
1433            << upperts
1434            << " AND TIMESYNC-(OBT0/1000)>"
1435            << lowerts
1436            << " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;";
1437        result = conn->Query(oss.str().c_str());
1438        if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str());
1439        //
1440        if ( !result ) throw -4;;
1441        found = true;
1442        if ( result->GetRowCount()<3 ){
1443          if ( IsDebug() ) printf(" AGH! no results!\n");
1444          found = false;
1445        } else {
1446          row = result->Next();      
1447          bn = (UInt_t)atoll(row->GetField(0));
1448          for ( Int_t r=1; r<result->GetRowCount() ;r++){
1449            if ( !row ) throw -4;
1450            if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0));
1451            if ( bn != (UInt_t)atoll(row->GetField(0)) ){
1452              if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0)));
1453              found = false;
1454            };
1455            row = result->Next();
1456          };
1457      };      };
1458      };
1459      //
1460      Int_t sgn = 0;
1461      //
1462      if ( !found && !BOOTNO ){
1463        throw -29;
1464    } else {    } else {
1465      if ( !this->GetBOOTnumber() ) return(2);      if ( afound ){
1466          this->SetBOOTnumber(bn);
1467          sgn = 8;
1468        };
1469    };    };
1470    //    //
1471    oss.str("");    oss.str("");
# Line 1179  Int_t PamelaDBOperations::assignBOOT_NUM Line 1475  Int_t PamelaDBOperations::assignBOOT_NUM
1475    conn->Query(oss.str().c_str());          conn->Query(oss.str().c_str());      
1476    //    //
1477    delete result;    delete result;
1478    return(0);    return(sgn);
1479  };  };
1480    
1481  /**  /**
# Line 1442  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1738  Bool_t PamelaDBOperations::IsRunAlreadyI
1738      if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){      if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){
1739        //        //
1740        if ( IsDebug() ) printf(" The new run has more events than the old one \n");        if ( IsDebug() ) printf(" The new run has more events than the old one \n");
1741        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1742        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1743        if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());        //       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
1744        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1745    //       conn->Query(oss.str().c_str());
1746        if ( signal ) signal = false;        if ( signal ) signal = false;
1747        goto gonext;              goto gonext;      
1748        //        //
# Line 1459  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1756  Bool_t PamelaDBOperations::IsRunAlreadyI
1756        //        //
1757        if ( IsDebug() ) printf(" The new run has the same number of events and the runheader the old one miss the runheader \n");        if ( IsDebug() ) printf(" The new run has the same number of events and the runheader the old one miss the runheader \n");
1758        //        //
1759        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1760        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1761        if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());        //       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
1762        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1763    //       conn->Query(oss.str().c_str());
1764        //        //
1765        if ( signal ) signal = false;        if ( signal ) signal = false;
1766        goto gonext;        goto gonext;
# Line 1476  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1774  Bool_t PamelaDBOperations::IsRunAlreadyI
1774        //        //
1775        if ( IsDebug() ) printf(" The new run has the same number of events, the runheader and the runtrailer the old one miss the runtrailer \n");        if ( IsDebug() ) printf(" The new run has the same number of events, the runheader and the runtrailer the old one miss the runtrailer \n");
1776        //        //
1777        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1778        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1779        if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());        //       oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";
1780        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1781    //       conn->Query(oss.str().c_str());
1782        if ( signal ) signal = false;        if ( signal ) signal = false;
1783        //        //
1784      };      };
# Line 1494  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1793  Bool_t PamelaDBOperations::IsRunAlreadyI
1793    delete result;    delete result;
1794    //    //
1795    if ( signal && IsDebug() ) printf(" The run has already been inserted \n");        if ( signal && IsDebug() ) printf(" The run has already been inserted \n");    
1796      if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n");
1797    return(signal);    return(signal);
1798  };  };
1799    
# Line 1515  void PamelaDBOperations::HandleRun(){   Line 1815  void PamelaDBOperations::HandleRun(){  
1815      //      //
1816      this->FillClass();      this->FillClass();
1817      //      //
1818      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);      if ( !IsRunAlreadyInserted() ){
1819          glrun->SetID(this->AssignRunID());
1820          glrun->SetID_RUN_FRAG(0);
1821          glrun->Fill_GL_RUN(conn);
1822        };
1823    } else {    } else {
1824      //      //
1825      if ( IsDebug() ) printf(" oh no! the distance between runheader and runtrailer seems wrong: check %llu pktt %llu \n",chkpkt,pktt);      if ( IsDebug() ) printf(" oh no! the distance between runheader and runtrailer seems wrong: check %llu pktt %llu \n",chkpkt,pktt);
# Line 1538  void PamelaDBOperations::HandleRunFragme Line 1842  void PamelaDBOperations::HandleRunFragme
1842    UInt_t rhfirstev = firstev;    UInt_t rhfirstev = firstev;
1843    UInt_t rtlastev = lastev;    UInt_t rtlastev = lastev;
1844    Bool_t found = false;    Bool_t found = false;
1845      Bool_t foundinrun = false;
1846    //    //
1847    TSQLResult *result = 0;    TSQLResult *result = 0;
1848    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
# Line 1564  void PamelaDBOperations::HandleRunFragme Line 1869  void PamelaDBOperations::HandleRunFragme
1869    //    //
1870    if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n");    if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n");
1871    //    //
1872    // can we find the other piece of the run in the GL_RUN_FRAGMENTS table?    // First of all insert the run in the fragment table...
1873    //        //
1874    if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is    oss.str("");
1875                                             // missing it no way we can found a piece in the frag table    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE "
1876      //        << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND ("
1877      oss.str("");        << " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND "
1878          << " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND ("
1879          << " RUNHEADER_OBT>=" << glrun->GetRUNHEADER_OBT() << " OR "
1880          << " RUNHEADER_PKT>=" << glrun->GetRUNHEADER_PKT() << ") AND ("
1881          << " RUNTRAILER_OBT<=" << glrun->GetRUNTRAILER_OBT() << " OR "
1882          << " RUNTRAILER_PKT<=" << glrun->GetRUNTRAILER_PKT() << ") ) OR "
1883          << " (RUNHEADER_TIME<=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1884          << " RUNTRAILER_TIME>=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND ("
1885          << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "
1886          << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("
1887          << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "
1888          << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));";
1889      //
1890      if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());
1891      result = conn->Query(oss.str().c_str());
1892      //
1893      if ( !result ) throw -4;
1894      //
1895      row = result->Next();
1896      //
1897      if ( !row ){
1898        //
1899        // no, insert this run in the GL_RUN_FRAGMENTS table (check if exist before!)
1900        //
1901        if ( IsDebug() ) printf(" The run is new \n");
1902        if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n");
1903        //
1904        glrun->SetID(this->AssignRunID());
1905        glrun->SetID_RUN_FRAG(0);
1906        glrun->Fill_GL_RUN_FRAGMENTS(conn);
1907        //
1908      } else {
1909        if ( IsDebug() ) printf(" The run is already present in the fragment table \n");
1910        return;
1911      };
1912      //
1913      //
1914      // can we find the other piece of the run in the GL_RUN_FRAGMENTS table?
1915      //    
1916      if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is
1917                                               // missing it no way we can found a piece in the frag table
1918        //
1919        oss.str("");
1920      oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE "      oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE "
1921          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
1922          << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME()          << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1923            << " ID != " << glrun->ID
1924          << " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!          << " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
1925      //      //
1926      if ( IsDebug() ) printf(" look for runheader in the fragments table: query is \n %s \n",oss.str().c_str());      if ( IsDebug() ) printf(" look for runheader in the fragments table: query is \n %s \n",oss.str().c_str());
# Line 1582  void PamelaDBOperations::HandleRunFragme Line 1930  void PamelaDBOperations::HandleRunFragme
1930      //      //
1931      row = result->Next();      row = result->Next();
1932      //      //
1933        if ( !row && NoFrag() ){
1934          //
1935          oss.str("");
1936          oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN WHERE "
1937              << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
1938              << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1939              << " ID != " << glrun->ID
1940              << " AND ID=ID_RUN_FRAG ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
1941          //
1942          if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
1943          result = conn->Query(oss.str().c_str());
1944          //
1945          if ( !result ) throw -4;
1946          //
1947          foundinrun = true;
1948          //
1949          row = result->Next();
1950          //
1951        };
1952        //
1953      if ( !row ){      if ( !row ){
1954        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
1955        found = false;        found = false;
# Line 1656  void PamelaDBOperations::HandleRunFragme Line 2024  void PamelaDBOperations::HandleRunFragme
2024        //        //
2025        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");
2026        //        //
2027          if ( foundinrun ){
2028            glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
2029            glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
2030          };
2031          //
2032        GL_RUN *glrun1 = new GL_RUN();        GL_RUN *glrun1 = new GL_RUN();
2033        //        //
2034        UInt_t idfrag = (UInt_t)atoll(row->GetField(0));        //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
2035        //        //
2036        oss.str("");        oss.str("");
2037        oss << " ID="<<row->GetField(0)<<";";        oss << " ID="<<row->GetField(0)<<";";
# Line 1687  void PamelaDBOperations::HandleRunFragme Line 2060  void PamelaDBOperations::HandleRunFragme
2060        };        };
2061        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: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
2062        //        //
       glrun1->SetID(0);  
2063        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
2064        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
2065        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());
# Line 1696  void PamelaDBOperations::HandleRunFragme Line 2068  void PamelaDBOperations::HandleRunFragme
2068        //        //
2069        glrun->SetEV_FROM(firstev);        glrun->SetEV_FROM(firstev);
2070        glrun->SetNEVENTS(lastev-firstev+1);        glrun->SetNEVENTS(lastev-firstev+1);
2071          //
2072        glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME());        glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME());
2073        glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT());        glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT());
2074        glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT());        glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT());
# Line 1715  void PamelaDBOperations::HandleRunFragme Line 2088  void PamelaDBOperations::HandleRunFragme
2088        //        //
2089        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
2090          //          //
2091            //      glrun->SetID(this->AssignRunID());
2092            glrun->SetID_RUN_FRAG(glrun1->GetID());
2093          glrun->Fill_GL_RUN(conn);          glrun->Fill_GL_RUN(conn);
2094          //          //
2095          // get id number          // set id number
         //  
         oss.str("");  
         oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "  
             << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "  
             << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";  
         //  
         if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         row = result->Next();  
         //  
         UInt_t idrun0 = 0;  
         if ( !row ){  
           throw -10;  
         } else {  
           idrun0 = (UInt_t)atoll(row->GetField(0));  
         };  
         //  
         glrun1->SetID_RUN_FRAG(idrun0);  
2096          //          //
2097            glrun1->SetID_RUN_FRAG(glrun->GetID());
2098          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
2099          //          //
         oss.str("");  
         oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";  
         //  
         if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         row = result->Next();  
         //  
         UInt_t idrun1 = 0;  
         if ( !row ){  
           throw -10;  
         } else {  
           idrun1 = (UInt_t)atoll(row->GetField(0));  
         };  
         //  
         oss.str("");  
         oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         //  
2100        };        };
       //  
       delete glrun1;  
       //  
2101        // delete old entry in fragment table        // delete old entry in fragment table
2102        //        //
2103        oss.str("");        glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2104          glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2105        //        //
2106        oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";";        delete glrun1;
2107        //        //
       if ( IsDebug() ) printf(" Delete from frag table the old run :\n query is \n %s \n",oss.str().c_str());  
       result =  conn->Query(oss.str().c_str());  
       if ( !result ) throw -4;  
2108        //        //
2109        return;        return;
2110        //        //
# Line 1784  void PamelaDBOperations::HandleRunFragme Line 2113  void PamelaDBOperations::HandleRunFragme
2113    };    };
2114    //    //
2115    if ( mistrail && rtlastev == lastev ) { // look for runtrailer (only when at the end of the file, if at the beginning and the runh is    if ( mistrail && rtlastev == lastev ) { // look for runtrailer (only when at the end of the file, if at the beginning and the runh is
2116                                             // missing it no way we can found a piece in the frag table      // missing it no way we can found a piece in the frag table
2117      //      //
2118      oss.str("");      oss.str("");
2119      oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN_FRAGMENTS WHERE "      oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN_FRAGMENTS WHERE "
2120          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
2121          << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME()          << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
2122            << " ID != " << glrun->ID
2123          << " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";          << " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
2124      //      //
2125      if ( IsDebug() ) printf(" look for runtrailer in the fragments table: query is \n %s \n",oss.str().c_str());      if ( IsDebug() ) printf(" look for runtrailer in the fragments table: query is \n %s \n",oss.str().c_str());
# Line 1799  void PamelaDBOperations::HandleRunFragme Line 2129  void PamelaDBOperations::HandleRunFragme
2129      //      //
2130      row = result->Next();      row = result->Next();
2131      //      //
2132        if ( !row && NoFrag() ){
2133          //
2134          oss.str("");
2135          oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN WHERE "
2136              << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
2137              << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
2138              << " ID != " << glrun->ID
2139              << " AND ID=ID_RUN_FRAG ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
2140          //
2141          if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
2142          result = conn->Query(oss.str().c_str());
2143          //
2144          if ( !result ) throw -4;
2145          //
2146          foundinrun = true;
2147          row = result->Next();
2148          //
2149        };
2150        //
2151      if ( !row ){      if ( !row ){
2152        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
2153        found = false;        found = false;
# Line 1872  void PamelaDBOperations::HandleRunFragme Line 2221  void PamelaDBOperations::HandleRunFragme
2221        //        //
2222        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");        if ( IsDebug() ) printf(" now you can handle the piece of the run \n ");
2223        //        //
2224          if ( foundinrun ){
2225            glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
2226            glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
2227          };
2228          //
2229        GL_RUN *glrun1 = new GL_RUN();        GL_RUN *glrun1 = new GL_RUN();
2230        //        //
2231        UInt_t idfrag = (UInt_t)atoll(row->GetField(0));        //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
2232        //        //
2233        oss.str("");        oss.str("");
2234        oss << " ID="<<row->GetField(0)<<";";        oss << " ID="<<row->GetField(0)<<";";
# Line 1911  void PamelaDBOperations::HandleRunFragme Line 2265  void PamelaDBOperations::HandleRunFragme
2265        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());
2266        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());
2267        //        //
       glrun1->SetID(0);  
2268        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());
2269        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());
2270        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());
# Line 1931  void PamelaDBOperations::HandleRunFragme Line 2284  void PamelaDBOperations::HandleRunFragme
2284        //        //
2285        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
2286          //          //
2287          glrun->Fill_GL_RUN(conn);          //      glrun->SetID(this->AssignRunID());
         //  
         // get id number  
         //  
         oss.str("");  
         oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "  
             << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "  
             << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";  
         //  
         if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());  
2288          //          //
2289          result =  conn->Query(oss.str().c_str());          glrun->SetID_RUN_FRAG(glrun1->GetID());
2290          if ( !result ) throw -4;          glrun->Fill_GL_RUN(conn);
         row = result->Next();  
         //  
         UInt_t idrun0 = 0;  
         if ( !row ){  
           throw -10;  
         } else {  
           idrun0 = (UInt_t)atoll(row->GetField(0));  
         };  
2291          //          //
2292          glrun1->SetID_RUN_FRAG(idrun0);          // set id number
2293          //          //
2294            glrun1->SetID_RUN_FRAG(glrun->GetID());
2295          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
2296          //          //
         oss.str("");  
         oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";  
         //  
         if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         row = result->Next();  
         //  
         UInt_t idrun1 = 0;  
         if ( !row ){  
           throw -10;  
         } else {  
           idrun1 = (UInt_t)atoll(row->GetField(0));  
         };  
         //  
         oss.str("");  
         oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         //  
2297        };        };
2298        //        //
2299        delete glrun1;        // delete old entries in fragment table
       //  
       // delete old entry in fragment table  
       //  
       oss.str("");  
       //  
       oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";";  
2300        //        //
2301        if ( IsDebug() ) printf(" Delete from frag table the old run :\n query is \n %s \n",oss.str().c_str());        glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2302        result =  conn->Query(oss.str().c_str());        glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
       if ( !result ) throw -4;  
2303        //        //
2304          delete glrun1;
2305        //        //
2306        return;        return;
2307        //        //
# Line 2029  void PamelaDBOperations::HandleRunFragme Line 2338  void PamelaDBOperations::HandleRunFragme
2338      //      //
2339      row = result->Next();      row = result->Next();
2340      //      //
2341      if ( !row ){      if ( row ){
       //  
       oss.str("");  
       oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE "  
           << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND ("  
           << " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND "  
           << " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND ("  
           << " RUNHEADER_OBT>=" << glrun->GetRUNHEADER_OBT() << " OR "  
           << " RUNHEADER_PKT>=" << glrun->GetRUNHEADER_PKT() << ") AND ("  
           << " RUNTRAILER_OBT<=" << glrun->GetRUNTRAILER_OBT() << " OR "  
           << " RUNTRAILER_PKT<=" << glrun->GetRUNTRAILER_PKT() << ") ) OR "  
           << " (RUNHEADER_TIME<=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "  
           << " RUNTRAILER_TIME>=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND ("  
           << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "  
           << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("  
           << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "  
           << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));";  
       //  
       if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());  
       result = conn->Query(oss.str().c_str());  
       //  
       if ( !result ) throw -4;  
       //  
       row = result->Next();  
       //  
       if ( !row ){  
         //  
         // no, insert this run in the GL_RUN_FRAGMENTS table (check if exist before!)  
         //  
         if ( IsDebug() ) printf(" The run is new \n");  
         if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n");  
         //  
         glrun->Fill_GL_RUN_FRAGMENTS(conn);  
         //  
       } else {  
         if ( IsDebug() ) printf(" The run is already present in the fragment table \n");  
       };  
     } else {  
2342        if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n");        if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n");
2343        } else {
2344          if ( NoFrag() ){
2345            glrun->SetID_RUN_FRAG(glrun->GetID());
2346            glrun->Fill_GL_RUN(conn);      
2347            glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2348          };    
2349      };      };
2350    };    };
2351    //    //
# Line 2093  void PamelaDBOperations::HandleMissingHo Line 2371  void PamelaDBOperations::HandleMissingHo
2371      //      //
2372      this->FillClass(mishead,mistrail,firstev,lastev);      this->FillClass(mishead,mistrail,firstev,lastev);
2373      //      //
2374      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2375          glrun->SetID(this->AssignRunID());
2376          glrun->SetID_RUN_FRAG(0);
2377          glrun->Fill_GL_RUN(conn);    
2378        };
2379      //      //
2380    };    };
2381    //    //
# Line 2260  Bool_t PamelaDBOperations::IsRunConsiste Line 2542  Bool_t PamelaDBOperations::IsRunConsiste
2542            //                  //      
2543            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2544            //            //
2545            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2546                glrun->SetID(this->AssignRunID());
2547                glrun->SetID_RUN_FRAG(0);
2548                glrun->Fill_GL_RUN(conn);      
2549              };
2550            //            //
2551            firstevno = lastentry + 1;            firstevno = lastentry + 1;
2552            //            //
# Line 2336  void PamelaDBOperations::HandleSuspiciou Line 2622  void PamelaDBOperations::HandleSuspiciou
2622      if ( IsDebug() ) printf(" Checking but no events in the run! \n");      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
2623      //      //
2624      this->FillClass();      this->FillClass();
2625      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2626          glrun->SetID(this->AssignRunID());
2627          glrun->SetID_RUN_FRAG(0);
2628          glrun->Fill_GL_RUN(conn);    
2629        };
2630      //      //
2631    } else {    } else {
2632      //      //
# Line 2364  void PamelaDBOperations::HandleSuspiciou Line 2654  void PamelaDBOperations::HandleSuspiciou
2654        if ( IsDebug() ) printf(" No packets but physics inside the run, I will consider it as good\n");        if ( IsDebug() ) printf(" No packets but physics inside the run, I will consider it as good\n");
2655        //        //
2656        this->FillClass();        this->FillClass();
2657        if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                if ( !IsRunAlreadyInserted() ){
2658            glrun->SetID(this->AssignRunID());
2659            glrun->SetID_RUN_FRAG(0);
2660            glrun->Fill_GL_RUN(conn);        
2661          };
2662        //        //
2663      } else {      } else {
2664        //        //
# Line 2462  void PamelaDBOperations::HandleSuspiciou Line 2756  void PamelaDBOperations::HandleSuspiciou
2756            //                  //      
2757            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2758            //            //
2759            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2760                glrun->SetID(this->AssignRunID());
2761                glrun->SetID_RUN_FRAG(0);
2762                glrun->Fill_GL_RUN(conn);      
2763              };
2764            //            //
2765            if ( i == lastev && checkfirst < check ){ // if the last event gives a wrong check...            if ( i == lastev && checkfirst < check ){ // if the last event gives a wrong check...
2766              //              //
# Line 2497  void PamelaDBOperations::HandleSuspiciou Line 2795  void PamelaDBOperations::HandleSuspiciou
2795              //                  //    
2796              this->SetCommonGLRUN(firstTime,lastTime);              this->SetCommonGLRUN(firstTime,lastTime);
2797              //              //
2798              if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                    if ( !IsRunAlreadyInserted() ){
2799                  glrun->SetID(this->AssignRunID());
2800                  glrun->SetID_RUN_FRAG(0);
2801                  glrun->Fill_GL_RUN(conn);      
2802                };
2803            };            };
2804            //            //
2805            firstevno = lastentry + 1;            firstevno = lastentry + 1;
# Line 2604  Int_t PamelaDBOperations::insertCALO_CAL Line 2906  Int_t PamelaDBOperations::insertCALO_CAL
2906                // no calibrations in the db contain our calibration                // no calibrations in the db contain our calibration
2907                //                //
2908                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
2909                if ( fromtime < 1150863400 ){                if ( fromtime < 1150871000 ){ //1150866904
2910                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
2911                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
2912                };                };
# Line 2667  Int_t PamelaDBOperations::insertCALO_CAL Line 2969  Int_t PamelaDBOperations::insertCALO_CAL
2969            //            //
2970          } else {          } else {
2971            //            //
2972            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 %i obt %i pkt %i \n",section,fromtime,obt,pkt);
2973            //            //
2974          };                };      
2975          //                //      
# Line 2749  void PamelaDBOperations::HandleTRK_CALIB Line 3051  void PamelaDBOperations::HandleTRK_CALIB
3051        // no calibrations in the db contain our calibration        // no calibrations in the db contain our calibration
3052        //        //
3053        if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n");        if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n");
3054        if ( fromtime < 1150863400 ) fromtime = 0; // the first flight calibration was taken at about 1150863300 s, this line allows to analyze first runs in raw mode        if ( fromtime < 1150871000 ) fromtime = 0; // the first flight calibration was taken at about 1150863300 s, this line allows to analyze first runs in raw mode
3055        //        //
3056        oss.str("");        oss.str("");
3057        oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "        oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "
# Line 2878  Int_t PamelaDBOperations::insertTRK_CALI Line 3180  Int_t PamelaDBOperations::insertTRK_CALI
3180      pkt1 = ph1->GetCounter();        pkt1 = ph1->GetCounter();  
3181      fromtime = this->GetAbsTime(ph1->GetOrbitalTime());        fromtime = this->GetAbsTime(ph1->GetOrbitalTime());  
3182      //      //
3183      valid = 1;  //     valid = 1;
3184      //  //     //
3185      if ( caltrk1->unpackError != 0 && caltrk1->good0 == 0 ) valid = 0;// CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT1  //     if ( caltrk1->unpackError != 0 && caltrk1->good0 == 0 ) valid = 0;// CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT1
3186      //      //
3187      //      //
3188      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
3189        //        //
3190        if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);        if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);
3191          //      
3192          valid = ValidateTrkCalib( caltrk1, eh1 );
3193          if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
3194        //        //
3195        // Do we have the second calibration packet?        // Do we have the second calibration packet?
3196        //        //
# Line 2902  Int_t PamelaDBOperations::insertTRK_CALI Line 3207  Int_t PamelaDBOperations::insertTRK_CALI
3207            obt2 = ph2->GetOrbitalTime();              obt2 = ph2->GetOrbitalTime();  
3208            pkt2 = ph2->GetCounter();              pkt2 = ph2->GetCounter();  
3209            //            //
3210            if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2  //        if ( caltrk2->unpackError != 0 || caltrk2->good0 == 0 ) valid = 0; // CONDITIONS ON THE GOODNESS OF THE CALIBRATION PKT2
3211            //                  //      
3212          } else {          } else {
3213            //            //
# Line 2933  Int_t PamelaDBOperations::insertTRK_CALI Line 3238  Int_t PamelaDBOperations::insertTRK_CALI
3238          //          //
3239          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 %i pkt %i t2 is %i is good \n",obt2,pkt2,t2);
3240          //          //
3241            UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
3242            if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
3243            valid = valid & valid2;
3244            //
3245          // Handle good calib          // Handle good calib
3246          //          //
3247          this->HandleTRK_CALIB(true,true);          this->HandleTRK_CALIB(true,true);
# Line 3105  Int_t PamelaDBOperations::insertS4_CALIB Line 3414  Int_t PamelaDBOperations::insertS4_CALIB
3414            // no calibrations in the db contain our calibration            // no calibrations in the db contain our calibration
3415            //            //
3416            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");
3417            if ( fromtime < 1150863400 ){            if ( fromtime < 1150871000 ){  
3418              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);
3419              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
3420            };            };
# Line 3178  Int_t PamelaDBOperations::insertS4_CALIB Line 3487  Int_t PamelaDBOperations::insertS4_CALIB
3487   * Scan the fragment table and move old fragments to the GL_RUN table   * Scan the fragment table and move old fragments to the GL_RUN table
3488   */   */
3489  Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){  Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){
3490      return(this->CleanGL_RUN_FRAGMENTS(""));
3491    };
3492    
3493    /**
3494     * Scan the fragment table and move old fragments to the GL_RUN table
3495     */
3496    Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(TString fcleanfile){
3497    //    //
3498    TSQLResult *result = 0;    TSQLResult *result = 0;
3499    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
# Line 3185  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3501  Int_t PamelaDBOperations::CleanGL_RUN_FR
3501    TSQLRow    *row2   = 0;    TSQLRow    *row2   = 0;
3502    //    //
3503    UInt_t moved = 0;    UInt_t moved = 0;
 //  UInt_t timelim = 0;  
 //  TDatime *time = new TDatime();  
3504    //    //
3505    stringstream oss;    stringstream oss;
3506    oss.str("");    oss.str("");
3507    //    //
3508    //    if ( !strcmp(fcleanfile.Data(),"") ){
3509    // check if there are entries older than "olderthan" seconds from now      //
3510    //      // check if there are entries older than "olderthan" seconds from now
3511    oss.str("");      //
3512    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"      oss.str("");
3513        << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";      oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
3514    //          << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";
3515    if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str());      //
3516    result = conn->Query(oss.str().c_str());      if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str());
3517        result = conn->Query(oss.str().c_str());
3518        //
3519      } else {
3520        oss.str("");
3521        oss << " SELECT ID FROM GL_ROOT WHERE NAME='" << fcleanfile.Data() << "';";
3522        if ( IsDebug() ) printf(" Getting ID_ROOT_L0 query %s \n",oss.str().c_str());
3523        result = conn->Query(oss.str().c_str());
3524        //    
3525        if ( result ){
3526          //
3527          row = result->Next();
3528          //
3529          if ( row ){
3530            oss.str("");
3531            oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
3532                << " ID_ROOT_L0=" << row->GetField(0) << ";";
3533            //
3534            if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS for ROOT file query is \n %s \n",oss.str().c_str());
3535            result = conn->Query(oss.str().c_str());
3536            //    
3537          };
3538        } else {
3539          return(2);
3540        };
3541      };
3542    //    //
3543    if ( result ){    if ( result ){
3544      //      //
3545      row = result->Next();      row = result->Next();
3546      //      //
3547      if ( row ){      while ( row ){
3548        //        //
3549        oss.str("");        oss.str("");
3550        oss << " ID= "<< row->GetField(0);        oss << " ID= "<< row->GetField(0);
# Line 3240  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3579  Int_t PamelaDBOperations::CleanGL_RUN_FR
3579          if ( IsDebug() ) printf(" The run is new \n");          if ( IsDebug() ) printf(" The run is new \n");
3580          if ( IsDebug() ) printf(" -> fill the DB \n");                if ( IsDebug() ) printf(" -> fill the DB \n");      
3581          //          //
3582          glrun->SetID(0);          //      glrun->SetID(this->AssignRunID()); we use the old run number!
3583            glrun->SetID_RUN_FRAG(glrun->GetID());
3584          glrun->Fill_GL_RUN(conn);            glrun->Fill_GL_RUN(conn);  
3585          //          //
3586          oss.str("");  //      oss.str("");
3587          oss << " SELECT ID FROM GL_RUN WHERE "  //      oss << " SELECT ID FROM GL_RUN WHERE "
3588              << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "  //          << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "
3589              << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "  //          << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "
3590              << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "  //          << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "
3591              << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "  //          << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "
3592              << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";  //          << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";
3593          //  //      //
3594          if ( IsDebug() ) printf(" Look for the ID of the inserted run: query is \n %s \n",oss.str().c_str());  //      if ( IsDebug() ) printf(" Look for the ID of the inserted run: query is \n %s \n",oss.str().c_str());
3595          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3596          //  //      //
3597          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3598          //  //      //
3599          row2 = result2->Next();  //      row2 = result2->Next();
3600          //  //      //
3601          if ( !row2 ) throw -25;  //      if ( !row2 ) throw -25;
3602          //  //      //
3603          oss.str("");  //      oss.str("");
3604          oss << " UPDATE GL_RUN SET ID_RUN_FRAG = " << row2->GetField(0) << " WHERE ID = " << row2->GetField(0);  //      oss << " UPDATE GL_RUN SET ID_RUN_FRAG = " << row2->GetField(0) << " WHERE ID = " << row2->GetField(0);
3605          if ( IsDebug() ) printf(" Update the ID_RUN_FRAG of the inserted run: query is \n %s \n",oss.str().c_str());  //      if ( IsDebug() ) printf(" Update the ID_RUN_FRAG of the inserted run: query is \n %s \n",oss.str().c_str());
3606          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3607          //  //      //
3608          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3609          //          //
3610          moved++;          moved++;
3611          //          //
# Line 3275  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3615  Int_t PamelaDBOperations::CleanGL_RUN_FR
3615        if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0));              if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0));      
3616        //        //
3617        //        //
3618        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
3619        oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0);  //      oss.str("");
3620        if ( IsDebug() ) printf(" Clean the GL_RUN_FRAGMENTS table: query is \n %s \n",oss.str().c_str());  //       oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0);
3621        result2 = conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" Clean the GL_RUN_FRAGMENTS table: query is \n %s \n",oss.str().c_str());
3622        //  //       result2 = conn->Query(oss.str().c_str());
3623        if ( !result2 ) throw -4;  //      //
3624        //  //      if ( !result2 ) throw -4;
3625    //      //
3626          row = result->Next();
3627      };      };
3628    };    };
3629    if ( IsDebug() ) printf(" Moved %u runs\n",moved);    if ( IsDebug() ) printf(" Moved %u runs\n",moved);
# Line 3292  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3634  Int_t PamelaDBOperations::CleanGL_RUN_FR
3634   * Check if runs are good, i.e. if the tracker calibration is correctly associated..   * Check if runs are good, i.e. if the tracker calibration is correctly associated..
3635   */   */
3636  Int_t PamelaDBOperations::ValidateRuns(){  Int_t PamelaDBOperations::ValidateRuns(){
3637      return(this->ValidateRuns(""));
3638    };
3639    
3640    /**
3641     * Check if runs are good, i.e. if the tracker calibration is correctly associated..
3642     */
3643    Int_t PamelaDBOperations::ValidateRuns(TString valfile){
3644    //    //
3645    TSQLResult *result = 0;    TSQLResult *result = 0;
3646    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
3647    //    //
3648      UInt_t calibtime = 50;
3649      //
3650    stringstream oss;    stringstream oss;
3651    oss.str("");    oss.str("");
3652    //    //
# Line 3304  Int_t PamelaDBOperations::ValidateRuns() Line 3655  Int_t PamelaDBOperations::ValidateRuns()
3655    // =======================================================    // =======================================================
3656    UInt_t t_stop  = 0;    UInt_t t_stop  = 0;
3657    UInt_t t_start = 0;    UInt_t t_start = 0;
3658    // --------------------------------------------------------------    if ( !strcmp(valfile.Data(),"") ) {
3659    // 1) get the OBT of the last run inserted after clean-time limit      // --------------------------------------------------------------
3660    // --------------------------------------------------------------      // 1) get the OBT of the last run inserted after clean-time limit
3661    oss.str("");      // --------------------------------------------------------------
3662    oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()      oss.str("");
3663            << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";      oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()
3664    if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());          << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3665    result = conn->Query(oss.str().c_str());      if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());
3666    if ( !result ) throw -4;      result = conn->Query(oss.str().c_str());
3667    if ( !result->GetRowCount() ) {      if ( !result ) throw -4;
3668            printf(" No runs to validate \n");      if ( !result->GetRowCount() ) {
3669            return(1);        printf(" No runs to validate \n");
3670    }else{        return(1);
3671          row = result->Next();      }else{
3672          t_start = (UInt_t)atoll(row->GetField(4));        row = result->Next();
3673    };          t_start = (UInt_t)atoll(row->GetField(4));
3674    // --------------------------------------------------------------      };  
3675    // 2) get the OBT of the last validated run      // --------------------------------------------------------------
3676    // --------------------------------------------------------------      // 2) get the OBT of the last validated run
3677    oss.str("");      // --------------------------------------------------------------
3678    oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start      oss.str("");
3679        <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";      oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
3680    if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());          <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3681    result = conn->Query(oss.str().c_str());      if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
3682    if ( !result ) throw -4;      result = conn->Query(oss.str().c_str());
3683    if ( result->GetRowCount() ){      if ( !result ) throw -4;
3684            row = result->Next();      if ( result->GetRowCount() ){
3685            t_stop = (UInt_t)atoll(row->GetField(4));        row = result->Next();
3686    };        t_stop = (UInt_t)atoll(row->GetField(4));
3687    if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);      };
3688    // --------------------------------------------------------------      if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);
3689    // now retrieves runs to be validated      // --------------------------------------------------------------
3690    // --------------------------------------------------------------      // now retrieves runs to be validated
3691    oss.str("");      // --------------------------------------------------------------
3692    oss << " SELECT * FROM GL_RUN  WHERE  RUNHEADER_TIME <=" << t_start;      oss.str("");
3693    oss << " AND RUNHEADER_TIME >="<< t_stop;      oss << " SELECT * FROM GL_RUN  WHERE  RUNHEADER_TIME <=" << t_start;
3694    oss << " ORDER BY RUNHEADER_TIME DESC;";      oss << " AND RUNHEADER_TIME >="<< t_stop;
3695  //  if ( IsDebug() )      oss << " ORDER BY RUNHEADER_TIME DESC;";
3696    if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());      if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());
3697    result = conn->Query(oss.str().c_str());      result = conn->Query(oss.str().c_str());
3698    if ( !result ) throw -4;    } else {
3699    if ( !result->GetRowCount() ) printf(" No runs to validate \n");      //
3700  //  printf("------------------------------------------------------------------------------- \n");      stringstream myquery;
3701          UInt_t myid = 0;
3702    Int_t nrow = 0;      myquery.str("");
3703    GL_RUN* this_run = new GL_RUN();      myquery << " SELECT ID FROM GL_ROOT where NAME='"<<valfile.Data() <<"';";
3704    GL_RUN* next_run = new GL_RUN();      //
3705    Int_t   nseq_max = 1000;      result = conn->Query(myquery.str().c_str());
3706  //  UInt_t* sequence = new UInt_t[100];      //
3707    vector<UInt_t> sequence(nseq_max);      row = result->Next();      
3708    Int_t   nseq = 0;      if( !row ){
3709    Bool_t CHECK = false;        if ( strcmp(valfile.Data(),GetRootName().Data()) ){
3710    Bool_t this_ONLINE = false;          if ( IsDebug() ) printf(" No file to be validated even if option \"-validate file\" was used!!\n");
3711    Bool_t next_ONLINE = false;          return(2);
3712    UInt_t t1=0,t2=0;        };
3713    // ---------------------------------------------------------------------------------        if ( IsDebug() ) printf(" No file to be validated (force mode)! \n");
3714    // - loop over runs, back in time,        return(0);
3715    // - select sequences of runs close in time (less than 60 s apart),      };
3716    //   which could be preceeded by a calibration      myid=(UInt_t)atoll(row->GetField(0));
3717    // - check if there might be a missing calibration      //
3718    // ---------------------------------------------------------------------------------      myquery.str("");
3719    while(1){      myquery << " SELECT MAX(RUNTRAILER_TIME),MIN(RUNHEADER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< myid <<";";
3720        //
3721        result = conn->Query(myquery.str().c_str());
3722        //
3723        row = result->Next();      
3724        if( !row->GetField(0) || !row->GetField(1)){
3725          //
3726          if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
3727          //
3728          return(0);
3729          //
3730        } else {  
3731          //
3732          UInt_t runhtime = (UInt_t)atoll(row->GetField(0));
3733          UInt_t runttime = (UInt_t)atoll(row->GetField(1));
3734          UInt_t caltime = 0;
3735          //
3736          myquery.str("");
3737          myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
3738          myquery << " order by FROM_TIME asc limit 1;";
3739          //
3740          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
3741          //
3742          //
3743          result = conn->Query(myquery.str().c_str());
3744          //
3745          row = result->Next();      
3746          if( !row ){
3747            caltime = runhtime;
3748          } else {
3749            caltime = (UInt_t)atoll(row->GetField(0));
3750          };
3751          //
3752          myquery.str("");
3753          myquery << " SELECT * from GL_RUN where RUNHEADER_TIME>="<< runttime <<" AND RUNHEADER_TIME<=" ;
3754          myquery << caltime << " order by RUNHEADER_TIME DESC";
3755          //
3756          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
3757          //
3758          result = conn->Query(myquery.str().c_str());
3759          //
3760        };
3761      };
3762      //
3763      if ( !result ) throw -4;
3764      if ( !result->GetRowCount() && IsDebug() ) printf(" No runs to validate \n");
3765      //
3766      Int_t nrow = 0;
3767      GL_RUN* this_run = new GL_RUN();
3768      GL_RUN* next_run = new GL_RUN();
3769      Int_t   nseq_max = 1000;
3770    //  UInt_t* sequence = new UInt_t[100];
3771      vector<UInt_t> sequence(nseq_max);
3772      Int_t   nseq = 0;
3773      Bool_t CHECK = false;
3774      Bool_t this_ONLINE = false;
3775      Bool_t next_ONLINE = false;
3776      UInt_t t1=0,t2=0;
3777      // ---------------------------------------------------------------------------------
3778      // - loop over runs, back in time,
3779      // - select sequences of runs close in time (less than calibtime s apart),
3780      //   which could be preceeded by a calibration
3781      // - check if there might be a missing calibration
3782      // ---------------------------------------------------------------------------------
3783      while(1){
3784                        
3785            row = result->Next();            row = result->Next();
3786            if( row == NULL ) break;            if( row == NULL ) break;
# Line 3401  Int_t PamelaDBOperations::ValidateRuns() Line 3816  Int_t PamelaDBOperations::ValidateRuns()
3816                                                        
3817                            if( this_run->ID == next_run->ID_RUN_FRAG ) interval = 0;     //=> run fragments                            if( this_run->ID == next_run->ID_RUN_FRAG ) interval = 0;     //=> run fragments
3818                                                        
3819                            if( interval >= 60 )CHECK = true;                     //more than 60 s => there might be a calibration                            if( interval >= calibtime )CHECK = true;                      //more than calibtime s => there might be a calibration
3820                                                        
3821                            if( !CHECK && this_run->VALIDATION ){                            if( !CHECK && this_run->VALIDATION ){
3822                                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);                                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);
# Line 3516  Int_t PamelaDBOperations::assignVALIDATI Line 3931  Int_t PamelaDBOperations::assignVALIDATI
3931    
3932    
3933    
3934    // Insert TLEs from file tlefilename in the table GL_TLE in the db
3935    // opened by conn, sorting them by date from older to newer, if each
3936    // TLE has not been alread inserted.
3937    Int_t PamelaDBOperations::populateTLE()//(TSQLServer *conn, char *tleFile)
3938    {
3939      fstream tlefile(tlefilename, ios::in);
3940    
3941      if ( !tlefile ) throw -7;
3942    
3943      vector<cTle*> ctles;
3944      vector<cTle*>::iterator iter;
3945      int present = 0;
3946    
3947      // Get three lines from tlefile, create a cTle object and put it
3948      // into ctles
3949      while(1) {
3950        cTle *tlef;
3951        string str1, str2, str3;
3952    
3953        getline(tlefile, str1);
3954        if(tlefile.eof()) break;
3955    
3956        getline(tlefile, str2);
3957        if(tlefile.eof()) break;
3958    
3959        getline(tlefile, str3);
3960        if(tlefile.eof()) break;
3961    
3962        // We now have three good lines for a cTle.
3963        tlef = new cTle(str1, str2, str3);
3964        ctles.push_back(tlef);
3965      }
3966    
3967      tlefile.close();
3968    
3969      // Sort by date
3970      sort(ctles.begin(), ctles.end(), compTLE);
3971    
3972      // Now we insert each TLE into the db
3973      for(iter = ctles.begin(); iter != ctles.end(); iter++) {
3974        cTle *tle = *iter;
3975    
3976        // Do nothing if it's already present in the db.  Just increase
3977        // the counter present.
3978        if (! isTlePresent(tle))
3979          {
3980            int status = insertTle(tle);
3981    
3982            // Insert query failed.  Return 1.
3983            if(status == EXIT_FAILURE) {
3984              
3985              if( IsDebug() ) {
3986                cerr << "Error: inserting TLE:" << endl
3987                     << tle->getName() << endl
3988                     << tle->getLine1() << endl
3989                     << tle->getLine2() << endl;
3990              }
3991    
3992              throw -4;
3993              return 1;
3994            }
3995    
3996          }
3997        else
3998          present++;
3999    
4000      }
4001    
4002      int inserted = ctles.size() - present;  // Number of inserted TLE.
4003      if ( IsDebug() )
4004        cout << "\nProcessed TLEs ranging from " << getTleDatetime(ctles[0]) << " to " << getTleDatetime(ctles[ctles.size()-1]) << "." << endl
4005             << inserted << " newly inserted TLEs out of " << ctles.size() << " processed." << endl;
4006    
4007      ctles.clear();
4008    
4009    
4010      // Return 2 if no new TLE has been inserted.  0 otherwise.
4011      if(! inserted ) return 2;
4012      return 0;
4013    }
4014    
4015    
4016    // Insert tle in the table GL_TLE using the connection conn.
4017    Int_t PamelaDBOperations::insertTle(cTle *tle)
4018    {
4019      stringstream oss;
4020      TSQLResult *result = 0;
4021    
4022      oss.str("");
4023      oss << " INSERT INTO GL_TLE (TLE1, TLE2, TLE3, FROM_TIME)"
4024          << " VALUES ( '"
4025          << tle->getName() << "', '"
4026          << tle->getLine1() << "', '"
4027          << tle->getLine2() << "', '"
4028          << getTleDatetime(tle) << "')";
4029    
4030      //  cout << oss.str().c_str() << endl;
4031      result = conn->Query(oss.str().c_str());
4032      if (result == NULL)
4033        return EXIT_FAILURE;
4034    
4035      return EXIT_SUCCESS;
4036    }
4037    
4038    
4039    // Return whether tle is already in the db connected by conn.
4040    bool PamelaDBOperations::isTlePresent(cTle *tle)
4041    {
4042      stringstream oss;
4043      TSQLResult *result = 0;
4044    
4045      oss.str("");
4046      oss << "SELECT * FROM GL_TLE WHERE FROM_TIME = '"
4047          << getTleDatetime(tle) << "'";
4048    
4049      result = conn->Query(oss.str().c_str());
4050      if (result == NULL) throw -4;
4051    
4052      if (result->GetRowCount())
4053        return true;
4054      else
4055        return false;
4056    }
4057    
4058    
4059    // Return whether the first TLE is dated early than the second
4060    bool compTLE (cTle *tle1, cTle *tle2)
4061    {
4062      return getTleJulian(tle1) < getTleJulian(tle2);
4063    }
4064    
4065    
4066    // Return the date of the tle using the format (year-2000)*1e3 +
4067    // julian day.  e.g. 6365 is the 31th Dec 2006.
4068    // It does *not* return a cJulian date.
4069    float getTleJulian(cTle *tle) {
4070      return tle->getField(cTle::FLD_EPOCHYEAR)*1e3 + tle->getField(cTle::FLD_EPOCHDAY);
4071    }
4072    
4073    
4074    // Return a string like YYYY-MM-DD hh:mm:ss, usable for mysql datetime
4075    // format.
4076    string getTleDatetime(cTle *tle)
4077    {
4078      int year, mon, day, hh, mm, ss;
4079      double dom; // day of month (is double!)
4080      stringstream date; // date in datetime format
4081    
4082      // create a cJulian from the date in tle
4083      cJulian jdate = cJulian( 2000 + (int) tle->getField(cTle::FLD_EPOCHYEAR), tle->getField(cTle::FLD_EPOCHDAY));
4084    
4085      // get year, month, day of month
4086      jdate.getComponent(&year, &mon, &dom);
4087    
4088      // build a datetime YYYY-MM-DD hh:mm:ss
4089      date.str("");
4090      day = (int) floor(dom);
4091      hh = (int) floor( (dom - day) * 24);
4092      mm = (int) floor( ((dom - day) * 24 - hh) * 60);
4093      ss = (int) floor( ((((dom - day) * 24 - hh) * 60 - mm) * 60));
4094      //  ms = (int) floor( (((((dom - day) * 24 - hh) * 60 - mm) * 60) - ss) * 1000);
4095    
4096      date << year << "-" << mon << "-" << day << " " << hh << ":" << mm << ":" << ss;
4097    
4098      return date.str();
4099    }
4100    
4101    /**
4102     * Remove a file from the DB, delete on cascade all entries related to that file
4103     * rearrange GL_RUN and GL_XXX_CALIB tables, turn off validation till the following
4104     * calibration
4105     **/
4106    Int_t PamelaDBOperations::removeFile(TString remfile){
4107      //
4108      // Determine ID_ROOT_L0 and ID_RAW
4109      //
4110      TSQLResult *pResult;
4111      TSQLRow *Row;
4112      stringstream myquery;
4113      //  
4114      myquery.str("");
4115      myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';";
4116      //
4117      pResult = conn->Query(myquery.str().c_str());
4118      //
4119      Row = pResult->Next();      
4120      if( !Row ){
4121        if ( strcmp(remfile.Data(),GetRootName().Data()) ){
4122          if ( IsDebug() ) printf(" No file to be removed even if option \"-remove file\" was used!!\n");
4123          return(1);
4124        };
4125        if ( IsDebug() ) printf(" No file to be removed (force mode)! \n");
4126        return(0);
4127      };
4128      //
4129      this->SetID_ROOT((UInt_t)atoll(Row->GetField(0)));
4130      this->SetID_RAW((UInt_t)atoll(Row->GetField(1)));
4131      //
4132      this->ValidationOFF();
4133      //
4134      this->RemoveCALIBS();
4135      //
4136      this->RemoveRUNS();
4137      //
4138      this->RemoveFILES();
4139      //
4140      this->SetID_ROOT(0);
4141      this->SetID_RAW(0);
4142      //
4143      return(0);
4144    };
4145    
4146    /**
4147     *
4148     * Set validation bit to zero for runs following the removing file till
4149     * 1) a run with TRK_CALIB_USED=140
4150     * 2) a run with VALIDATION = 0
4151     * 3) the next calibration
4152     *
4153     **/
4154    void PamelaDBOperations::ValidationOFF(){
4155      TSQLResult *pResult;
4156      TSQLRow *Row;
4157      stringstream myquery;
4158      Int_t unv = 0;
4159      //select ID from GL_RUN where RUNHEADER_TIME>=1152671382 AND (VALIDATION=0 OR TRK_CALIB_USED=104) order by RUNHEADER_TIME asc limit 1;
4160      myquery.str("");
4161      myquery << " SELECT MAX(RUNTRAILER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4162      //
4163      pResult = conn->Query(myquery.str().c_str());
4164      //
4165      Row = pResult->Next();      
4166      if( !Row->GetField(0) ){
4167        //
4168        if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
4169        //
4170      } else {  
4171        //
4172        UInt_t runhtime = (UInt_t)atoll(Row->GetField(0));
4173        UInt_t caltime = 0;
4174        //
4175        myquery.str("");
4176        myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
4177        myquery << " order by FROM_TIME asc limit 1;";
4178        //
4179        if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4180        //
4181        //
4182        delete pResult;
4183        pResult = conn->Query(myquery.str().c_str());
4184        //
4185        Row = pResult->Next();      
4186        if( !Row ){
4187          caltime = runhtime;
4188        } else {
4189          caltime = (UInt_t)atoll(Row->GetField(0));
4190        };
4191        //
4192        myquery.str("");
4193        myquery << " SELECT ID,RUNHEADER_TIME from GL_RUN where RUNHEADER_TIME>="<< runhtime <<" AND (VALIDATION=0 OR TRK_CALIB_USED=104 OR RUNHEADER_TIME>" ;
4194        myquery << caltime << ") order by RUNHEADER_TIME asc LIMIT 1";
4195        //
4196        if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4197        //
4198        pResult = conn->Query(myquery.str().c_str());
4199        //
4200        Row = pResult->Next();      
4201        if( !Row ){
4202          //
4203          if ( IsDebug() ) printf(" NO RUN NEED TO BE UNVALIDATED \n");
4204          //
4205        } else {
4206          myquery.str("");
4207          myquery << " SELECT ID from GL_RUN where RUNHEADER_TIME<"<< Row->GetField(1) <<" AND ";
4208          myquery << " RUNHEADER_TIME>=" <<runhtime;
4209          myquery << " order by RUNHEADER_TIME asc;";
4210          //
4211          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4212          //
4213          pResult = conn->Query(myquery.str().c_str());
4214          //
4215          Row = pResult->Next();  
4216          while ( Row ){
4217            //      
4218            unv++;
4219            this->assignVALIDATION((UInt_t)atoll(Row->GetField(0)), false);
4220            Row = pResult->Next();  
4221            //
4222          };
4223        };
4224      };
4225      if ( IsDebug() ) printf(" %i runs have been unvalidated \n",unv);
4226    };
4227    
4228    /**
4229     *
4230     * Rearrange GL_RUN table and remove runs
4231     *
4232     **/
4233    void PamelaDBOperations::RemoveRUNS(){
4234      TSQLResult *pResult;
4235      TSQLRow *Row;
4236      stringstream myquery;
4237      UInt_t drun = 0;
4238      GL_RUN *delrun = new GL_RUN();
4239      //
4240      myquery.str("");
4241      myquery << " SELECT ID FROM GL_RUN where ID_RUN_FRAG=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4242      //
4243      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4244      //
4245      pResult = conn->Query(myquery.str().c_str());
4246      //
4247      Row = pResult->Next();    
4248      //
4249      //
4250      if ( !Row ){
4251        if ( IsDebug() ) printf(" No run with ID_RUN_FRAG=0 belonged to this file \n");
4252      } else {
4253        if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");    
4254        while ( Row ){
4255          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4256          if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
4257          drun++;
4258          Row = pResult->Next();    
4259        };
4260      };
4261      //
4262      //
4263      myquery.str("");
4264      myquery << " SELECT ID,ID_RUN_FRAG FROM GL_RUN where ID_RUN_FRAG!=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4265      //
4266      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4267      //
4268      pResult = conn->Query(myquery.str().c_str());
4269      //
4270      Row = pResult->Next();    
4271      //
4272      if ( !Row ){
4273        if ( IsDebug() ) printf(" No run with ID_RUN_FRAG!=0 belonged to this file \n");
4274      } else {
4275        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");    
4276        while ( Row ){
4277          if ( IsDebug() ) printf(" restore run %i \n",(UInt_t)atoll(Row->GetField(1)));
4278          delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");
4279          if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(1)));
4280          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");
4281          if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){
4282            if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
4283            delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4284          };
4285          drun++;
4286          Row = pResult->Next();    
4287        };
4288      };
4289      //
4290      if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN table \n",drun);    
4291      //
4292      //
4293      //
4294      drun = 0;
4295      //
4296      myquery.str("");
4297      myquery << " SELECT ID_TRASH FROM GL_RUN_TRASH where BELONGED_TO='GL_RUN_FRAGMENTS' AND ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4298      //
4299      pResult = conn->Query(myquery.str().c_str());
4300      //
4301      Row = pResult->Next();    
4302      //
4303      if ( !Row ){
4304        if ( IsDebug() ) printf(" No run from GL_RUN_FRAGMENTS table in the trash table for this file \n");
4305      } else {
4306        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");    
4307        while ( Row ){
4308          if ( IsDebug() ) printf(" del run idtrash %i \n",(UInt_t)atoll(Row->GetField(0)));
4309          myquery.str("");
4310          myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";      
4311          conn->Query(myquery.str().c_str());
4312          drun++;
4313          Row = pResult->Next();    
4314        };
4315      };
4316      //  
4317      if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_TRASH table \n",drun);    
4318      //
4319      //
4320      //
4321      drun = 0;
4322      //
4323      myquery.str("");
4324      myquery << " SELECT ID FROM GL_RUN_FRAGMENTS where ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4325      //
4326      pResult = conn->Query(myquery.str().c_str());
4327      //
4328      Row = pResult->Next();    
4329      //
4330      if ( !Row ){
4331        if ( IsDebug() ) printf(" No run in the GL_RUN_FRAGMENTS table for this file \n");
4332      } else {
4333        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");    
4334        while ( Row ){
4335          if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));
4336          myquery.str("");
4337          myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";
4338          conn->Query(myquery.str().c_str());
4339          drun++;
4340          Row = pResult->Next();    
4341        };
4342      };
4343      //  
4344      if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_FRAGMENTS table \n",drun);    
4345      //
4346      //
4347      //
4348      delete delrun;
4349      //
4350    };
4351    
4352    
4353    /**
4354     *
4355     * Rearrange calibration tables
4356     *
4357     **/
4358    void PamelaDBOperations::RemoveFILES(){
4359      stringstream myquery;
4360      //
4361      myquery.str("");
4362      myquery << " DELETE FROM GL_RAW WHERE ID=" <<this->GetID_RAW() <<";";
4363      //
4364      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4365      //
4366      conn->Query(myquery.str().c_str());
4367      //
4368    };
4369    
4370    /**
4371     *
4372     * Rearrange calibration tables
4373     *
4374     **/
4375    void PamelaDBOperations::RemoveCALIBS(){
4376      TSQLResult *pResult;
4377      TSQLRow *Row;
4378      stringstream myquery;
4379      //  
4380      //
4381      // Calorimeter
4382      //
4383      for (Int_t section = 0; section < 4; section++){
4384        myquery.str("");
4385        myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_CALO_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<" AND ";
4386        myquery << " SECTION=" << section << ";";
4387        //
4388        pResult = conn->Query(myquery.str().c_str());
4389        //
4390        Row = pResult->Next();      
4391        if( !Row->GetField(0) || !Row->GetField(1) ){
4392          //
4393          if ( IsDebug() ) printf(" NO CALO CALIBRATION SECTION %i ASSOCIATED TO THIS FILE! \n",section);
4394          //
4395        } else {
4396          //
4397          myquery.str("");
4398          myquery << " UPDATE GL_CALO_CALIB SET TO_TIME=" << Row->GetField(1);
4399          myquery << " WHERE TO_TIME="<< Row->GetField(0) << " AND ";
4400          myquery << " SECTION=" << section << ";";
4401          //
4402          pResult = conn->Query(myquery.str().c_str());
4403          //
4404          if( !pResult ){
4405            //
4406            if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
4407            //
4408            throw -4;
4409            //
4410          };
4411          //
4412        };
4413      };
4414      myquery.str("");
4415      myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4416      //
4417      pResult = conn->Query(myquery.str().c_str());
4418      //
4419      if( !pResult ){
4420        //
4421        if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
4422        //
4423        throw -4;
4424        //
4425      };
4426      //
4427      // Tracker
4428      //
4429      myquery.str("");
4430      myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_TRK_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4431      //
4432      pResult = conn->Query(myquery.str().c_str());
4433      //
4434      Row = pResult->Next();      
4435      if( !Row->GetField(0) || !Row->GetField(1) ){
4436        //
4437        if ( IsDebug() ) printf(" NO TRK CALIBRATION ASSOCIATED TO THIS FILE! \n");
4438        //
4439      } else {
4440        //
4441        myquery.str("");
4442        myquery << " UPDATE GL_TRK_CALIB SET TO_TIME=" << Row->GetField(1);
4443        myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
4444        //
4445        pResult = conn->Query(myquery.str().c_str());
4446        //
4447        if( !pResult ){
4448          //
4449          if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
4450          //
4451          throw -4;
4452          //
4453        };
4454        //
4455        myquery.str("");
4456        myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4457        //
4458        pResult = conn->Query(myquery.str().c_str());
4459        //
4460        if( !pResult ){
4461          //
4462          if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
4463          //
4464          throw -4;
4465          //
4466        };
4467      };
4468      //
4469      //
4470      // S4
4471      //
4472      myquery.str("");
4473      myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_S4_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4474      //
4475      pResult = conn->Query(myquery.str().c_str());
4476      //
4477      Row = pResult->Next();      
4478      if( !Row->GetField(0) || !Row->GetField(1) ){
4479        //
4480        if ( IsDebug() ) printf(" NO S4 CALIBRATION ASSOCIATED TO THIS FILE! \n");
4481        //
4482      } else {
4483        //
4484        myquery.str("");
4485        myquery << " UPDATE GL_S4_CALIB SET TO_TIME=" << Row->GetField(1);
4486        myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
4487        //
4488        pResult = conn->Query(myquery.str().c_str());
4489        //
4490        if( !pResult ){
4491          //
4492          if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
4493          //
4494          throw -4;
4495          //
4496        };
4497        //
4498        myquery.str("");
4499        myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4500        //
4501        pResult = conn->Query(myquery.str().c_str());
4502        //
4503        if( !pResult ){
4504          //
4505          if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
4506          //
4507          throw -4;
4508          //
4509        };
4510        //
4511      };
4512    };
4513    
4514    /**
4515     *
4516     * Rearrange calibration tables
4517     *
4518     **/
4519    UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){
4520    
4521        Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0};
4522        UInt_t timeaftercalib=120000; //2000;
4523    //  ----------
4524    //  Check CRCs
4525    //  ----------
4526        for(Int_t ipkt=0; ipkt<6; ipkt++){
4527            if( caltrk->crc_hcal[ipkt] )return 0; // :-(
4528            for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-(
4529        }
4530    //  -----------------------
4531    //  Check missing packets:
4532    //  -----------------------
4533    //    Readout order:
4534    //    ------------------
4535    //    DSP   packet board
4536    //    ------------------
4537    //    12    0      1
4538    //    10    1      1
4539    //     8    2      1
4540    //     4    3      1
4541    //     6    4      1
4542    //     2    5      1
4543    //    ------------------
4544    //    11    0      2
4545    //     9    1      2
4546    //     7    2      2
4547    //     3    3      2
4548    //     5    4      2
4549    //     1    5      2
4550    //    ------------------
4551    //  -------------------------------------------------
4552    //  Check if it is first or second calibration packet
4553    //  -------------------------------------------------
4554        UInt_t build=0;
4555        TString classname = caltrk->GetName();
4556        UInt_t base=0;
4557        UInt_t mask=0;
4558        if(classname.Contains("CalibTrk1Event")){
4559            base=12;
4560            mask=0x03F000;
4561        }
4562        if(classname.Contains("CalibTrk2Event")){
4563            base=18;
4564            mask=0xFC0000;
4565        }
4566    //  -------------------------------------------------
4567    //  Count number of packets and set build variable
4568    //  -------------------------------------------------
4569        Int_t  npkts=0;
4570        for(Int_t ipkt=0; ipkt<6; ipkt++){
4571            if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){
4572                npkts++;
4573                build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) );
4574            }
4575        }
4576    //    if( npkts==6 )return 1; // :-)
4577    
4578    //    cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl;
4579    
4580    //  -----------------------------------------------
4581    //  If missing packets: check the acq configuration
4582    //  (some DSPs might be excluded from acquisition)
4583    //  -----------------------------------------------
4584    
4585    //  -----------------------------------------------
4586    //  retrieve the first run header after calib
4587    //  -----------------------------------------------
4588         PacketType *pctp;
4589         EventCounter *cod;
4590         cod = eh->GetCounter();
4591         Int_t irun = cod->Get(pctp->RunHeader);
4592         TTree *rh=(TTree*)file->Get("RunHeader");
4593         if ( !rh || rh->IsZombie() ) throw -17;
4594         if( rh->GetEntries() == irun ){
4595              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1)  -- cannot validate :-( "<<endl;
4596             return 0; // :-(
4597         }
4598    
4599         RunHeaderEvent *run  = 0;
4600         EventHeader    *hrun = 0;
4601         rh->SetBranchAddress("RunHeader", &run);
4602         rh->SetBranchAddress("Header", &hrun);
4603         rh->GetEntry(irun);
4604    //     cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl;
4605    
4606         if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){
4607              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl;
4608             return 0; // :-(
4609         }
4610        
4611         if( !run->RM_ACQ_AFTER_CALIB ){
4612              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0    -- cannot validate :-( "<<endl;
4613             return 0; // :-(
4614         }
4615    
4616         UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime());
4617         if( dtime > timeaftercalib ){
4618             if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl;
4619             return 0; // :-(
4620         }
4621        
4622    
4623    
4624         if( (run->ACQ_BUILD_INFO & mask) != build ){
4625             if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl;
4626             return 0; // :-(
4627         }
4628         return 1; // :-)
4629    
4630    }
4631    
4632    /**
4633     *
4634     * Check the DB (only for overlapping runs at the moment)
4635     *
4636     **/
4637    UInt_t PamelaDBOperations::Check(){
4638      //
4639      UInt_t test = 0;
4640      //
4641      UInt_t thisrht = 0;
4642      UInt_t thisrtt = 0;
4643      UInt_t thisid = 0;
4644      UInt_t prevrht = 0;
4645      UInt_t prevrtt = 0;
4646      UInt_t previd = 0;
4647      //
4648      UInt_t prevl0id = 0;
4649      UInt_t thisl0id = 0;
4650      //
4651      stringstream oss;
4652      TSQLResult *result = 0;
4653      TSQLRow    *row    = 0;
4654      TSQLResult *result2 = 0;
4655      TSQLRow    *row2    = 0;
4656      oss.str("");
4657      oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN order by RUNHEADER_TIME asc;";
4658      //  oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;";
4659      result = conn->Query(oss.str().c_str());
4660      //
4661      if ( !result ) throw -4;;
4662      //
4663      row = result->Next();
4664      //
4665      while ( row ){
4666        thisid = (UInt_t)atoll(row->GetField(0));
4667        thisl0id = (UInt_t)atoll(row->GetField(1));
4668        thisrht = (UInt_t)atoll(row->GetField(2));
4669        thisrtt = (UInt_t)atoll(row->GetField(3));
4670        //
4671        //    if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){
4672        //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){
4673        if ( (thisrht < prevrtt) && (thisrht != prevrht) ){
4674          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4675          printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid);
4676          TString prevf = "";
4677          TString thisf = "";
4678          oss.str("");
4679          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4680          result2 = conn->Query(oss.str().c_str());
4681          if ( !result2 ) throw -4;;
4682          row2 = result2->Next();
4683          prevf = (TString)row2->GetField(0);
4684          oss.str("");
4685          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4686          result2 = conn->Query(oss.str().c_str());
4687          if ( !result2 ) throw -4;;
4688          row2 = result2->Next();
4689          thisf = (TString)row2->GetField(0);
4690          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4691          test = 1;
4692        };
4693        //
4694        if ( (thisrtt < prevrht) && (thisrht != prevrht) ){
4695          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4696          printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid);
4697          TString prevf = "";
4698          TString thisf = "";
4699          oss.str("");
4700          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4701          result2 = conn->Query(oss.str().c_str());
4702          if ( !result2 ) throw -4;;
4703          row2 = result2->Next();
4704          prevf = (TString)row2->GetField(0);
4705          oss.str("");
4706          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4707          result2 = conn->Query(oss.str().c_str());
4708          if ( !result2 ) throw -4;;
4709          row2 = result2->Next();
4710          thisf = (TString)row2->GetField(0);
4711          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4712          test = 1;
4713        };
4714        //
4715        if ( (thisrht > thisrtt) && (thisrht != prevrht) ){
4716          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4717          printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid);
4718          TString prevf = "";
4719          TString thisf = "";
4720          oss.str("");
4721          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4722          result2 = conn->Query(oss.str().c_str());
4723          if ( !result2 ) throw -4;;
4724          row2 = result2->Next();
4725          prevf = (TString)row2->GetField(0);
4726          oss.str("");
4727          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4728          result2 = conn->Query(oss.str().c_str());
4729          if ( !result2 ) throw -4;;
4730          row2 = result2->Next();
4731          thisf = (TString)row2->GetField(0);
4732          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4733          test = 1;
4734        };
4735    
4736        //
4737        prevrht = thisrht;
4738        prevrtt = thisrtt;
4739        previd = thisid;
4740        prevl0id = thisl0id;
4741        row = result->Next();
4742      };
4743      //
4744      return(test);
4745      //
4746    };

Legend:
Removed from v.1.8  
changed lines
  Added in v.1.29

  ViewVC Help
Powered by ViewVC 1.1.23