/[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.11 by mocchiut, Wed Sep 13 08:34:01 2006 UTC revision 1.46 by mocchiut, Thu Aug 30 15:04:44 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){  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);
# Line 69  PamelaDBOperations::PamelaDBOperations(T Line 81  PamelaDBOperations::PamelaDBOperations(T
81    INSERT_ROOT = !filerootname.IsNull();    INSERT_ROOT = !filerootname.IsNull();
82    if( INSERT_ROOT ){    if( INSERT_ROOT ){
83      this->SetRootName(filerootname);      this->SetRootName(filerootname);
84        this->SetOrbitNo(dwinput);
85      file = TFile::Open(this->GetRootName().Data());      file = TFile::Open(this->GetRootName().Data());
86      } else {
87        this->SetRootName("");
88    };    };
89    //    //
90    this->SetID_RAW(0);    this->SetID_RAW(0);
# Line 99  void PamelaDBOperations::Close(){ Line 114  void PamelaDBOperations::Close(){
114  //  //
115  void PamelaDBOperations::CheckValidate(Long64_t olderthan){  void PamelaDBOperations::CheckValidate(Long64_t olderthan){
116    clean_time = new TDatime();    clean_time = new TDatime();
117      //
118    if(olderthan >= 0){    if(olderthan >= 0){
119      VALIDATE = true;      VALIDATE = true;
120      UInt_t timelim = 0;      UInt_t timelim = 0;
121      timelim =  (UInt_t)clean_time->Convert() - olderthan;      timelim =  (UInt_t)clean_time->Convert(true) - olderthan;
122      clean_time->Set(timelim,false);      clean_time->Set(timelim,false);
123    };    };
124  };  };
# Line 114  void PamelaDBOperations::CheckValidate(L Line 130  void PamelaDBOperations::CheckValidate(L
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 142  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 182  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 190  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 205  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 238  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 265  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 282  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 %u pkt_num %u \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(32);
408    //    //
409    if ( !nevent ) return(2);    if ( !nevent ) return(64);
410    //    //
411    if ( nevent < 2 ) return(4);    if ( nevent < 2 ) return(128);
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 314  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 %u upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
451            break;            break;
452          };          };
453          //          //
# Line 354  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);        //  
658        goto kikko;      };
659        //      break;      //
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 671  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 > ((Long64_t)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 ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u  (numeric_limits<UInt_t>::max()/2) is %u  \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2));
886    //    //
887    if ( obt > (obtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){    if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){
888      if ( (obt-numeric_limits<UInt_t>::max()) < 0  ){      if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)(obt+numeric_limits<UInt_t>::max()));
889        return((ULong64_t)(numeric_limits<UInt_t>::max()-obt));      return((Long64_t)(obt+numeric_limits<UInt_t>::max()));
890      } else {    };
891        return((ULong64_t)(obt-numeric_limits<UInt_t>::max()));    //
892      };    if ( obt > ((Long64_t)pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
893        if ( IsDebug() ) printf(" pobtfirst*2 %lld \n",((Long64_t)pobtfirst*2));
894        if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
895        return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
896    };    };
897    //    //
898    return((ULong64_t)obt);    if ( IsDebug() ) printf(" as it is %lld \n",(Long64_t)obt);
899      return((Long64_t)obt);
900  };  };
901    
902  /**  /**
# Line 824  Int_t PamelaDBOperations::insertPamelaRa Line 1011  Int_t PamelaDBOperations::insertPamelaRa
1011    if ( idr ) return(1);    if ( idr ) return(1);
1012    //    //
1013    oss.str("");    oss.str("");
1014    oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"    if ( STATIC ){
1015        << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";      oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('"
1016            << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')";
1017      } else {
1018        oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')";
1019      };
1020    if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;    if ( conn->Query(oss.str().c_str()) == 0 ) throw -4;
1021    //    //
1022    idr =  this->SetID_RAW();    idr =  this->SetID_RAW();
# Line 844  Int_t PamelaDBOperations::insertPamelaGL Line 1035  Int_t PamelaDBOperations::insertPamelaGL
1035    TSQLResult *result = 0;    TSQLResult *result = 0;
1036    TSQLRow *row = 0;    TSQLRow *row = 0;
1037    UInt_t t0 = 0;    UInt_t t0 = 0;
1038      Int_t signal = 0;
1039      //
1040      signal = this->SetUpperLimits();
1041    //    //
1042    stringstream   oss;    stringstream   oss;
1043    //      //  
1044    if ( this->GetID_RAW() == 0 )  throw -11;    if ( this->GetID_RAW() == 0 )  throw -11;
1045    //    //
1046    oss.str("");    oss.str("");
1047    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='"
1048        << id << " AND TO_ID_RAW >= "        << this->GetRawFile().Data() << "';";
       << id << ";";  
1049    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
1050    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1051    if ( !result ) throw -10;    if ( !result ) throw -10;
1052    row = result->Next();    row = result->Next();
   if ( !row ) throw -10;  
1053    //    //
1054    t0 = (UInt_t)TDatime(row->GetField(0)).Convert();    if ( !row ){
1055        oss.str("");
1056        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< "
1057            << dworbit << " order by FROM_ORBIT desc limit 1;";
1058        if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
1059        result = conn->Query(oss.str().c_str());
1060        if ( !result ) throw -10;
1061        row = result->Next();
1062        if ( !row ) throw -10;
1063      };
1064      //
1065      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);
1066      t0 = (UInt_t)tu.GetSec();
1067      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));
1068      //
1069    /*    /*
1070     * Verify that the TIMESYNC have been not already processed     * Verify that the TIMESYNC have been not already processed
1071     */     */
# Line 874  Int_t PamelaDBOperations::insertPamelaGL Line 1080  Int_t PamelaDBOperations::insertPamelaGL
1080    if (result == NULL) throw -10;    if (result == NULL) throw -10;
1081    row = result->Next();    row = result->Next();
1082    if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){    if ((row != NULL) && ((UInt_t)atoll(row->GetField(0)) > 0)){
1083        if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
1084      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;
1085        //
1086        tsync = (UInt_t)atoll(row->GetField(2));
1087        obt0 = (UInt_t)atoll(row->GetField(1));
1088        //
1089      return(1);      return(1);
1090    };    };
1091    //    //
1092    TTree *T = 0;    TTree *T = 0;
   Int_t signal = 0;  
1093    //    //
1094    UInt_t nevent = 0;    UInt_t nevent = 0;
1095    UInt_t recEntries = 0;    UInt_t recEntries = 0;
# Line 1068  Int_t PamelaDBOperations::insertPamelaGL Line 1278  Int_t PamelaDBOperations::insertPamelaGL
1278    conn->Query(oss.str().c_str());    conn->Query(oss.str().c_str());
1279    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());
1280    //    //
1281      if ( IsDebug() ) printf(" found a timesync t0 is %u \n",t0);
1282    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;    toffset = (UInt_t)TSYNC - (UInt_t)(this->OBT(OBT)/1000) + t0;
1283    //    //
1284      tsync = TSYNC;
1285      obt0 = OBT;
1286      //
1287    delete result;    delete result;
1288    return(signal);    return(signal);
1289  }  }
# Line 1085  Int_t PamelaDBOperations::insertPamelaRo Line 1299  Int_t PamelaDBOperations::insertPamelaRo
1299    UInt_t idtimesync = 0;    UInt_t idtimesync = 0;
1300    //    //
1301    oss.str("");    oss.str("");
1302    oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "    if ( STATIC ){
1303        << " LEFT JOIN GL_ROOT "      oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
1304        << " ON GL_RAW.ID = GL_ROOT.ID_RAW "          << " LEFT JOIN GL_ROOT "
1305        << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "          << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
1306        << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";          << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND "
1307            << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
1308      } else {
1309        oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW "
1310            << " LEFT JOIN GL_ROOT "
1311            << " ON GL_RAW.ID = GL_ROOT.ID_RAW "
1312            << " WHERE GL_RAW.PATH = '$PAM_RAW' AND "
1313            << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID ";
1314      };
1315    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1316    //    //
1317    if ( !result ) throw -12;    if ( !result ) throw -12;
# Line 1116  Int_t PamelaDBOperations::insertPamelaRo Line 1338  Int_t PamelaDBOperations::insertPamelaRo
1338    idtimesync = (UInt_t)atoll(row->GetField(0));    idtimesync = (UInt_t)atoll(row->GetField(0));
1339    //    //
1340    oss.str("");    oss.str("");
1341    oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"    if ( STATIC ){
1342        << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";      oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
1343            << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')";
1344      } else {
1345        oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('"
1346            << this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')";
1347      };
1348    //    //
1349    if (conn->Query(oss.str().c_str()) == 0) throw -4;    if (conn->Query(oss.str().c_str()) == 0) throw -4;
1350    //    //
# Line 1144  Int_t PamelaDBOperations::assignBOOT_NUM Line 1371  Int_t PamelaDBOperations::assignBOOT_NUM
1371    TSQLResult *result = 0;    TSQLResult *result = 0;
1372    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
1373    oss.str("");    oss.str("");
1374    oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "    if ( STATIC ){
1375        << " PATH = '" << this->GetRawPath().Data() << "' AND "      oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
1376        << " NAME = '" << this->GetRawFile().Data() << "' ";          << " PATH = '" << this->GetRawPath().Data() << "' AND "
1377            << " NAME = '" << this->GetRawFile().Data() << "' ";
1378      } else {
1379        oss << "SELECT ID, BOOT_NUMBER FROM GL_RAW WHERE "
1380            << " PATH = '$PAM_RAW' AND "
1381            << " NAME = '" << this->GetRawFile().Data() << "' ";
1382      };
1383    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
1384    //    //
1385    if ( !result ) throw -4;;    if ( !result ) throw -4;;
# Line 1169  Int_t PamelaDBOperations::assignBOOT_NUM Line 1402  Int_t PamelaDBOperations::assignBOOT_NUM
1402    VarDumpEvent  *vde = 0;    VarDumpEvent  *vde = 0;
1403    VarDumpRecord *vdr = 0;    VarDumpRecord *vdr = 0;
1404    //    //
1405      Bool_t found = false;
1406    trDumpEv->SetBranchAddress("VarDump", &vde);    trDumpEv->SetBranchAddress("VarDump", &vde);
1407    if ( trDumpEv->GetEntries() > 0 ){    if ( trDumpEv->GetEntries() > 0 ){
1408      Bool_t found = false;      found = false;
1409      for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){      for ( Int_t i = 0; i < trDumpEv->GetEntries(); i++){
1410        trDumpEv->GetEntry(i);        trDumpEv->GetEntry(i);
1411        vde->Records->GetEntries();        //      vde->Records->GetEntries();
1412        if ( vde->Records->GetEntries()>0 ){        if ( vde->Records->GetEntries()>5 ){
1413          found = true;          found = true;
1414          goto fill;          goto fill;
1415        };        };
# Line 1188  Int_t PamelaDBOperations::assignBOOT_NUM Line 1422  Int_t PamelaDBOperations::assignBOOT_NUM
1422        this->SetBOOTnumber((Int_t)vdr->VAR_VALUE);        this->SetBOOTnumber((Int_t)vdr->VAR_VALUE);
1423        //        //
1424      } else {      } else {
1425        if ( !this->GetBOOTnumber() ) return(4);        if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(4);
1426        };
1427      } else {
1428        if ( !this->GetBOOTnumber() && !this->AutoBoot()) return(2);
1429      };
1430      //
1431      UInt_t bn = 0;
1432      Bool_t afound = false;
1433      if ( !found && this->AutoBoot()){
1434        afound = true;
1435        //
1436        // Search for other files with similar timesync
1437        //
1438        if ( IsDebug() ) printf(" tsync %u obt0 %u \n",tsync,obt0);
1439        UInt_t upperts = tsync-(obt0/1000)+5;
1440        UInt_t lowerts = tsync-(obt0/1000)-5;
1441        oss.str("");
1442        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)<"
1443            << upperts
1444            << " AND TIMESYNC-(OBT0/1000)>"
1445            << lowerts
1446            << " AND GL_RAW.BOOT_NUMBER>0 GROUP BY GL_TIMESYNC.OBT0;";
1447        result = conn->Query(oss.str().c_str());
1448        if ( IsDebug() ) printf(" Query the GL_TIMESYNC table to find boot number:\n %s \n",oss.str().c_str());
1449        //
1450        if ( !result ) throw -4;;
1451        found = true;
1452        if ( result->GetRowCount()<3 ){
1453          if ( IsDebug() ) printf(" AGH! no results!\n");
1454          found = false;
1455        } else {
1456          row = result->Next();      
1457          bn = (UInt_t)atoll(row->GetField(0));
1458          for ( Int_t r=1; r<result->GetRowCount() ;r++){
1459            if ( !row ) throw -4;
1460            if ( IsDebug() ) printf(" BOOT number is %s \n",row->GetField(0));
1461            if ( bn != (UInt_t)atoll(row->GetField(0)) ){
1462              if ( IsDebug() ) printf(" AGH! bn = %u here instead %u \n",bn,(UInt_t)atoll(row->GetField(0)));
1463              found = false;
1464            };
1465            row = result->Next();
1466          };
1467      };      };
1468      };
1469      //
1470      Int_t sgn = 0;
1471      //
1472      if ( !found && !BOOTNO ){
1473        throw -29;
1474    } else {    } else {
1475      if ( !this->GetBOOTnumber() ) return(2);      if ( afound ){
1476          this->SetBOOTnumber(bn);
1477          sgn = 8;
1478        };
1479    };    };
1480    //    //
1481    oss.str("");    oss.str("");
# Line 1201  Int_t PamelaDBOperations::assignBOOT_NUM Line 1485  Int_t PamelaDBOperations::assignBOOT_NUM
1485    conn->Query(oss.str().c_str());          conn->Query(oss.str().c_str());      
1486    //    //
1487    delete result;    delete result;
1488    return(0);    return(sgn);
1489  };  };
1490    
1491  /**  /**
# Line 1214  Int_t PamelaDBOperations::insertPamelaRU Line 1498  Int_t PamelaDBOperations::insertPamelaRU
1498    stringstream oss;    stringstream oss;
1499    oss.str("");    oss.str("");
1500    //    //
1501    signal = this->SetUpperLimits();    //  signal = this->SetUpperLimits();
1502    //    //
1503    // loop on runheader and runtrailer events    // loop on runheader and runtrailer events
1504    //    //
# Line 1232  Int_t PamelaDBOperations::insertPamelaRU Line 1516  Int_t PamelaDBOperations::insertPamelaRU
1516    rt->SetBranchAddress("RunTrailer", &runt);    rt->SetBranchAddress("RunTrailer", &runt);
1517    rt->SetBranchAddress("Header", &eht);    rt->SetBranchAddress("Header", &eht);
1518    //    //
1519      TTree *T = (TTree*)file->Get("Physics");
1520      if ( !T || T->IsZombie() ) throw -16;
1521      EventHeader *eh = 0;
1522      T->SetBranchAddress("Header", &eh);
1523      //
1524      if ( !(rh->GetEntries()) && !(rt->GetEntries()) && !(T->GetEntries()) ) return(16);
1525      //
1526    UInt_t obtt = 0;    UInt_t obtt = 0;
1527    UInt_t obth = 0;    UInt_t obth = 0;
1528    UInt_t pktt = 0;    UInt_t pktt = 0;
# Line 1265  Int_t PamelaDBOperations::insertPamelaRU Line 1556  Int_t PamelaDBOperations::insertPamelaRU
1556        //        //
1557        if ( !ptt &&  !(ptht+1) ){        if ( !ptt &&  !(ptht+1) ){
1558          //          //
1559          if ( IsDebug() ) printf(" Piece of run at the beginning of the file %i %i %i \n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Piece of run at the beginning of the file %u %u %u \n",ptht,pth,ptt);
1560          //          //
1561          this->HandleRunFragments(true,false,0,(evbeft-1));          this->HandleRunFragments(true,false,0,(evbeft-1));
1562          //          //
1563          //          //
1564        } else if ( pth == ptht ){        } else if ( pth == ptht ){
1565          //          //
1566          if ( IsDebug() ) printf(" Missing header %i %i %i\n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Missing header %u %u %u\n",ptht,pth,ptt);
1567          //          //
1568          if ( (ptt-1) < 0 ) throw -15; // should never arrive here!          if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
1569          rt->GetEntry(ptt-1);          rt->GetEntry(ptt-1);
# Line 1281  Int_t PamelaDBOperations::insertPamelaRU Line 1572  Int_t PamelaDBOperations::insertPamelaRU
1572          rt->GetEntry(ptt);          rt->GetEntry(ptt);
1573          pht = eht->GetPscuHeader();          pht = eht->GetPscuHeader();
1574          //          //
1575          if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %i %i %i \n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %u %u %u \n",ptht,pth,ptt);
1576          if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %i %i %i \n",pkth,obth,obtt);          if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %u %u %u \n",pkth,obth,obtt);
1577          //          //
1578          this->HandleMissingHoT(true,false,evbefh,evbeft-1);          this->HandleMissingHoT(true,false,evbefh,evbeft-1);
1579          //          //
1580        } else {        } else {
1581          //          //
1582          if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %i %i %i\n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %u %u %u\n",ptht,pth,ptt);
1583          //          //
1584          rh->GetEntry(ptht);          rh->GetEntry(ptht);
1585          phh = ehh->GetPscuHeader();          phh = ehh->GetPscuHeader();
# Line 1296  Int_t PamelaDBOperations::insertPamelaRU Line 1587  Int_t PamelaDBOperations::insertPamelaRU
1587          obth = phh->GetOrbitalTime();          obth = phh->GetOrbitalTime();
1588          cod = ehh->GetCounter();          cod = ehh->GetCounter();
1589          evbefh = cod->Get(pctp->Physics);          evbefh = cod->Get(pctp->Physics);
1590          if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);          if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
1591          //          //
1592          // handle this run          // handle this run
1593          //          //
# Line 1306  Int_t PamelaDBOperations::insertPamelaRU Line 1597  Int_t PamelaDBOperations::insertPamelaRU
1597          //          //
1598          if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){          if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){
1599            //            //
1600            if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER \n");            if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER evbefh = %u \n",evbefh);        
           //  
           this->HandleRunFragments(true,true,0,(evbefh-1));  
1601            //            //
1602              if ( evbefh == 0 ) {
1603                //
1604                signal = 8;
1605                if ( IsDebug() ) printf(" Not supported yet: run with no events, no runtrailer, no runheader \n");
1606                //
1607              } else {
1608                //
1609                this->HandleRunFragments(true,true,0,(evbefh-1));
1610                //
1611              };
1612          };          };
1613          //          //
1614          //          //
1615          if ( (ptht - pth) > 1 ){          if ( (ptht - pth) > 1 ){
1616            //            //
1617            if ( IsDebug() ) printf(" Missing runtrailers! \n");            if ( IsDebug() ) printf(" Missing runtrailers! \n");
1618            if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %u %u %u \n",ptht,pth,ptt);
1619            // is not the consecutive header            // is not the consecutive header
1620            while ( pth != ptht ){                      while ( pth != ptht ){          
1621              //              //
# Line 1337  Int_t PamelaDBOperations::insertPamelaRU Line 1636  Int_t PamelaDBOperations::insertPamelaRU
1636              obth = phh->GetOrbitalTime();              obth = phh->GetOrbitalTime();
1637              evbefh = cod->Get(pctp->Physics);              evbefh = cod->Get(pctp->Physics);
1638              //              //
1639              if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %i %i %i \n",ptht,pth,ptt);              if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %u %u %u \n",ptht,pth,ptt);
1640              if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);              if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
1641              //              //
1642              this->HandleMissingHoT(false,true,evbefh,evbeft-1);              this->HandleMissingHoT(false,true,evbefh,evbeft-1);
1643              //              //
# Line 1347  Int_t PamelaDBOperations::insertPamelaRU Line 1646  Int_t PamelaDBOperations::insertPamelaRU
1646          } else if ( !(ptht - pth) ){          } else if ( !(ptht - pth) ){
1647            //            //
1648            if ( IsDebug() ) printf(" Missing runheader! \n");            if ( IsDebug() ) printf(" Missing runheader! \n");
1649            if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %u %u %u \n",ptht,pth,ptt);
1650            if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n");            if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n");
1651            //            //
1652          } else {          } else {
# Line 1368  Int_t PamelaDBOperations::insertPamelaRU Line 1667  Int_t PamelaDBOperations::insertPamelaRU
1667            obth = phh->GetOrbitalTime();            obth = phh->GetOrbitalTime();
1668            cod = ehh->GetCounter();            cod = ehh->GetCounter();
1669            evbefh = cod->Get(pctp->Physics);            evbefh = cod->Get(pctp->Physics);
1670            if ( IsDebug() ) printf(" Piece of run at the end of file %i %i %i \n",pkth,obth,obtt);            if ( IsDebug() ) printf(" Piece of run at the end of file %u %u %u \n",pkth,obth,obtt);
1671            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u %u %u \n",ptht,pth,ptt);
1672            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i \n",rhev);            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u \n",rhev);
1673            //            //
1674            this->HandleRunFragments(false,true,evbefh,upperentry);            this->HandleRunFragments(false,true,evbefh,upperentry);
1675          } else {          } else {
# Line 1418  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1717  Bool_t PamelaDBOperations::IsRunAlreadyI
1717    //                                  ( runhead_time <= _our_runhead_time &&  runtrail_time >= _our_runtrail_time) &&    //                                  ( runhead_time <= _our_runhead_time &&  runtrail_time >= _our_runtrail_time) &&
1718    //                                  ( runhead_obt <= _our_runheadobt || runhead_pkt <= _our_runheadpkt ) &&    //                                  ( runhead_obt <= _our_runheadobt || runhead_pkt <= _our_runheadpkt ) &&
1719    //                                  ( runtrail_obt <= _our_runtrailobt || runtrail_pkt <= _our_runtrailpkt ) )    //                                  ( runtrail_obt <= _our_runtrailobt || runtrail_pkt <= _our_runtrailpkt ) )
1720      //                                  ||
1721      //                                  ( runhead_time = _our_runhead_time &&  runtrail_time = _our_runtrail_time && nevents > 100 )
1722    //                                  )    //                                  )
1723    //    //
1724    oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE "    oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE "
# Line 1433  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1734  Bool_t PamelaDBOperations::IsRunAlreadyI
1734        << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "        << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "
1735        << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("        << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("
1736        << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "        << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "
1737        << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));";        << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") OR "
1738          << " (RUNHEADER_TIME=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "                //  these two lines in a certain way disable the patch below...
1739          << " RUNTRAILER_TIME=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND NEVENTS>100)"    //
1740          << " ));";
1741    //    //
1742    if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());
1743    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
# Line 1455  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1759  Bool_t PamelaDBOperations::IsRunAlreadyI
1759      //      //
1760      // the run has already been inserted      // the run has already been inserted
1761      //      //
1762      //    return(true);  //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here      if ( signal && IsDebug() ) printf(" The run has already been inserted \n");    
1763        return(true);  //<<<<<<<<<<<<<<<<<<<<<<<< patch follows, uncomment here
1764      //      //
1765      // PATCH!      // PATCH!
1766      // we keep the processing run if (in order of growing importance) 1) we have the runtrailer while the old run doesn't have it 2) we have the runheader      // we keep the processing run if (in order of growing importance) 1) we have the runtrailer while the old run doesn't have it 2) we have the runheader
# Line 1464  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1769  Bool_t PamelaDBOperations::IsRunAlreadyI
1769      if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){      if ( glrun->GetNEVENTS() > (UInt_t)atoll(row->GetField(1)) ){
1770        //        //
1771        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");
1772        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1773        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1774        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) <<";";
1775        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1776    //       conn->Query(oss.str().c_str());
1777        if ( signal ) signal = false;        if ( signal ) signal = false;
1778        goto gonext;              goto gonext;      
1779        //        //
# Line 1481  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1787  Bool_t PamelaDBOperations::IsRunAlreadyI
1787        //        //
1788        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");
1789        //        //
1790        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1791        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1792        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) <<";";
1793        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1794    //       conn->Query(oss.str().c_str());
1795        //        //
1796        if ( signal ) signal = false;        if ( signal ) signal = false;
1797        goto gonext;        goto gonext;
# Line 1498  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1805  Bool_t PamelaDBOperations::IsRunAlreadyI
1805        //        //
1806        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");
1807        //        //
1808        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
1809        oss << "DELETE FROM GL_RUN WHERE ID=" << row->GetField(0) <<";";  //       oss.str("");
1810        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) <<";";
1811        conn->Query(oss.str().c_str());  //       if ( IsDebug() ) printf(" delete the run entry: query is \n %s \n",oss.str().c_str());      
1812    //       conn->Query(oss.str().c_str());
1813        if ( signal ) signal = false;        if ( signal ) signal = false;
1814        //        //
1815      };      };
# Line 1516  Bool_t PamelaDBOperations::IsRunAlreadyI Line 1824  Bool_t PamelaDBOperations::IsRunAlreadyI
1824    delete result;    delete result;
1825    //    //
1826    if ( signal && IsDebug() ) printf(" The run has already been inserted \n");        if ( signal && IsDebug() ) printf(" The run has already been inserted \n");    
1827      if ( !signal && IsDebug() ) printf(" The run existed and was deleted, fill the DB \n");
1828    return(signal);    return(signal);
1829  };  };
1830    
# Line 1537  void PamelaDBOperations::HandleRun(){   Line 1846  void PamelaDBOperations::HandleRun(){  
1846      //      //
1847      this->FillClass();      this->FillClass();
1848      //      //
1849      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);      if ( !IsRunAlreadyInserted() ){
1850          glrun->SetID(this->AssignRunID());
1851          glrun->SetID_RUN_FRAG(0);
1852          glrun->Fill_GL_RUN(conn);
1853        };
1854    } else {    } else {
1855      //      //
1856      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 1560  void PamelaDBOperations::HandleRunFragme Line 1873  void PamelaDBOperations::HandleRunFragme
1873    UInt_t rhfirstev = firstev;    UInt_t rhfirstev = firstev;
1874    UInt_t rtlastev = lastev;    UInt_t rtlastev = lastev;
1875    Bool_t found = false;    Bool_t found = false;
1876      Bool_t foundinrun = false;
1877    //    //
1878    TSQLResult *result = 0;    TSQLResult *result = 0;
1879    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
# Line 1586  void PamelaDBOperations::HandleRunFragme Line 1900  void PamelaDBOperations::HandleRunFragme
1900    //    //
1901    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");
1902    //    //
1903    // 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...
1904    //        //
1905    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("");
1906                                             // missing it no way we can found a piece in the frag table    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE "
1907          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND ("
1908          << " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND "
1909          << " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND ("
1910          << " RUNHEADER_OBT>=" << glrun->GetRUNHEADER_OBT() << " OR "
1911          << " RUNHEADER_PKT>=" << glrun->GetRUNHEADER_PKT() << ") AND ("
1912          << " RUNTRAILER_OBT<=" << glrun->GetRUNTRAILER_OBT() << " OR "
1913          << " RUNTRAILER_PKT<=" << glrun->GetRUNTRAILER_PKT() << ") ) OR "
1914          << " (RUNHEADER_TIME<=" << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1915          << " RUNTRAILER_TIME>=" << (UInt_t)glrun->GetRUNTRAILER_TIME() <<" AND ("
1916          << " RUNHEADER_OBT<=" << glrun->GetRUNHEADER_OBT() << " OR "
1917          << " RUNHEADER_PKT<=" << glrun->GetRUNHEADER_PKT() << ") AND ("
1918          << " RUNTRAILER_OBT>=" << glrun->GetRUNTRAILER_OBT() << " OR "
1919          << " RUNTRAILER_PKT>=" << glrun->GetRUNTRAILER_PKT() << ") ));";
1920      //
1921      if ( IsDebug() ) printf(" check if run has been inserted: query is \n %s \n",oss.str().c_str());
1922      result = conn->Query(oss.str().c_str());
1923      //
1924      if ( !result ) throw -4;
1925      //
1926      row = result->Next();
1927      //
1928      if ( !row ){
1929      //      //
1930      oss.str("");      // no, insert this run in the GL_RUN_FRAGMENTS table (check if exist before!)
1931      oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE "      //
1932          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "      if ( IsDebug() ) printf(" The run is new \n");
1933          << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME()      if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n");
1934          << " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!      //
1935        glrun->SetID(this->AssignRunID());
1936        glrun->SetID_RUN_FRAG(0);
1937        glrun->Fill_GL_RUN_FRAGMENTS(conn);
1938        //
1939      } else {
1940        if ( IsDebug() ) printf(" The run is already present in the fragment table \n");
1941        return;
1942      };
1943      //
1944      //
1945      // can we find the other piece of the run in the GL_RUN_FRAGMENTS table?
1946      //    
1947      if ( mishead && rhfirstev == firstev ) { // look for runheader (only when at the beginning of the file, if at the end and the runh is
1948                                               // missing it no way we can found a piece in the frag table
1949        //
1950        oss.str("");
1951        oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN_FRAGMENTS WHERE "
1952            << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
1953            << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1954            << " ID != " << glrun->ID
1955            << " ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
1956      //      //
1957      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());
1958      result = conn->Query(oss.str().c_str());      result = conn->Query(oss.str().c_str());
# Line 1604  void PamelaDBOperations::HandleRunFragme Line 1961  void PamelaDBOperations::HandleRunFragme
1961      //      //
1962      row = result->Next();      row = result->Next();
1963      //      //
1964        if ( !row && NoFrag() ){
1965          //
1966          oss.str("");
1967          oss << " SELECT ID,TRK_CALIB_USED,RUNTRAILER_TIME,RUNTRAILER_OBT,RUNHEADER_PKT,RUNTRAILER_PKT FROM GL_RUN WHERE "
1968              << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
1969              << " RUNHEADER_TIME <= " << (UInt_t)glrun->GetRUNHEADER_TIME() << " AND "
1970              << " ID != " << glrun->ID
1971              << " AND ID=ID_RUN_FRAG ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; // DESC NOT ASC!!
1972          //
1973          if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
1974          result = conn->Query(oss.str().c_str());
1975          //
1976          if ( !result ) throw -4;
1977          //
1978          foundinrun = true;
1979          //
1980          row = result->Next();
1981          //
1982        };
1983        //
1984      if ( !row ){      if ( !row ){
1985        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
1986        found = false;        found = false;
# Line 1678  void PamelaDBOperations::HandleRunFragme Line 2055  void PamelaDBOperations::HandleRunFragme
2055        //        //
2056        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 ");
2057        //        //
2058          if ( foundinrun ){
2059            glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
2060            glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
2061          };
2062          //
2063        GL_RUN *glrun1 = new GL_RUN();        GL_RUN *glrun1 = new GL_RUN();
2064        //        //
2065        UInt_t idfrag = (UInt_t)atoll(row->GetField(0));        //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
2066        //        //
2067        oss.str("");        oss.str("");
2068        oss << " ID="<<row->GetField(0)<<";";        oss << " ID="<<row->GetField(0)<<";";
# Line 1693  void PamelaDBOperations::HandleRunFragme Line 2075  void PamelaDBOperations::HandleRunFragme
2075        ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT());        ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT());
2076        UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT());        UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT());
2077        ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT());        ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT());
2078        if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);        if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
2079        TTree *T= 0;        TTree *T= 0;
2080        T = (TTree*)file->Get("Physics");        T = (TTree*)file->Get("Physics");
2081        if ( !T || T->IsZombie() ) throw -16;        if ( !T || T->IsZombie() ) throw -16;
# Line 1707  void PamelaDBOperations::HandleRunFragme Line 2089  void PamelaDBOperations::HandleRunFragme
2089          bobt = OBT(ph->GetOrbitalTime());                bobt = OBT(ph->GetOrbitalTime());      
2090          firstev++;          firstev++;
2091        };        };
2092        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);        if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
2093        //        //
       glrun1->SetID(0);  
2094        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
2095        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
2096        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());
# Line 1718  void PamelaDBOperations::HandleRunFragme Line 2099  void PamelaDBOperations::HandleRunFragme
2099        //        //
2100        glrun->SetEV_FROM(firstev);        glrun->SetEV_FROM(firstev);
2101        glrun->SetNEVENTS(lastev-firstev+1);        glrun->SetNEVENTS(lastev-firstev+1);
2102          //
2103        glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME());        glrun->SetRUNHEADER_TIME(glrun1->GetRUNHEADER_TIME());
2104        glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT());        glrun->SetRUNHEADER_OBT(glrun1->GetRUNHEADER_OBT());
2105        glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT());        glrun->SetRUNHEADER_PKT(glrun1->GetRUNHEADER_PKT());
# Line 1737  void PamelaDBOperations::HandleRunFragme Line 2119  void PamelaDBOperations::HandleRunFragme
2119        //        //
2120        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
2121          //          //
2122            //      glrun->SetID(this->AssignRunID());
2123            glrun->SetID_RUN_FRAG(glrun1->GetID());
2124          glrun->Fill_GL_RUN(conn);          glrun->Fill_GL_RUN(conn);
2125          //          //
2126          // 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);  
2127          //          //
2128            glrun1->SetID_RUN_FRAG(glrun->GetID());
2129          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
2130          //          //
         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;  
         //  
2131        };        };
       //  
       delete glrun1;  
       //  
2132        // delete old entry in fragment table        // delete old entry in fragment table
2133        //        //
2134        oss.str("");        glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2135          glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2136        //        //
2137        oss << " DELETE FROM GL_RUN_FRAGMENTS WHERE ID = " << idfrag << ";";        delete glrun1;
2138        //        //
       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;  
2139        //        //
2140        return;        return;
2141        //        //
# Line 1806  void PamelaDBOperations::HandleRunFragme Line 2144  void PamelaDBOperations::HandleRunFragme
2144    };    };
2145    //    //
2146    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
2147                                             // 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
2148      //      //
2149      oss.str("");      oss.str("");
2150      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 "
2151          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "          << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
2152          << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME()          << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
2153            << " ID != " << glrun->ID
2154          << " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";          << " ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
2155      //      //
2156      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 1821  void PamelaDBOperations::HandleRunFragme Line 2160  void PamelaDBOperations::HandleRunFragme
2160      //      //
2161      row = result->Next();      row = result->Next();
2162      //      //
2163        if ( !row && NoFrag() ){
2164          //
2165          oss.str("");
2166          oss << " SELECT ID,PKT_COUNTER,RUNHEADER_TIME,RUNHEADER_OBT,RUNTRAILER_PKT,RUNHEADER_PKT FROM GL_RUN WHERE "
2167              << " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND "
2168              << " RUNTRAILER_TIME >= " << (UInt_t)glrun->GetRUNTRAILER_TIME() << " AND "
2169              << " ID != " << glrun->ID
2170              << " AND ID=ID_RUN_FRAG ORDER BY RUNTRAILER_TIME ASC LIMIT 1;";
2171          //
2172          if ( IsDebug() ) printf(" look for runheader in the GL_RUN table: query is \n %s \n",oss.str().c_str());
2173          result = conn->Query(oss.str().c_str());
2174          //
2175          if ( !result ) throw -4;
2176          //
2177          foundinrun = true;
2178          row = result->Next();
2179          //
2180        };
2181        //
2182      if ( !row ){      if ( !row ){
2183        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");        if ( IsDebug() ) printf(" the corresponding piece has NOT been found \n");
2184        found = false;        found = false;
# Line 1894  void PamelaDBOperations::HandleRunFragme Line 2252  void PamelaDBOperations::HandleRunFragme
2252        //        //
2253        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 ");
2254        //        //
2255          if ( foundinrun ){
2256            glrun->RestoreRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
2257            glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN");
2258          };
2259          //
2260        GL_RUN *glrun1 = new GL_RUN();        GL_RUN *glrun1 = new GL_RUN();
2261        //        //
2262        UInt_t idfrag = (UInt_t)atoll(row->GetField(0));        //      UInt_t idfrag = (UInt_t)atoll(row->GetField(0));
2263        //        //
2264        oss.str("");        oss.str("");
2265        oss << " ID="<<row->GetField(0)<<";";        oss << " ID="<<row->GetField(0)<<";";
# Line 1909  void PamelaDBOperations::HandleRunFragme Line 2272  void PamelaDBOperations::HandleRunFragme
2272        ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT());        ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT());
2273        UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT());        UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT());
2274        ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT());        ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT());
2275        if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);        if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
2276        TTree *T= 0;        TTree *T= 0;
2277        T = (TTree*)file->Get("Physics");        T = (TTree*)file->Get("Physics");
2278        if ( !T || T->IsZombie() ) throw -16;        if ( !T || T->IsZombie() ) throw -16;
# Line 1923  void PamelaDBOperations::HandleRunFragme Line 2286  void PamelaDBOperations::HandleRunFragme
2286          aobt = OBT(ph->GetOrbitalTime());                aobt = OBT(ph->GetOrbitalTime());      
2287          lastev--;          lastev--;
2288        };        };
2289        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);        if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
2290        //        //
2291        glrun->SetEV_TO(lastev);        glrun->SetEV_TO(lastev);
2292        glrun->SetNEVENTS(lastev-firstev+1);        glrun->SetNEVENTS(lastev-firstev+1);
# Line 1933  void PamelaDBOperations::HandleRunFragme Line 2296  void PamelaDBOperations::HandleRunFragme
2296        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());
2297        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());
2298        //        //
       glrun1->SetID(0);  
2299        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());
2300        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());
2301        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());
# Line 1953  void PamelaDBOperations::HandleRunFragme Line 2315  void PamelaDBOperations::HandleRunFragme
2315        //        //
2316        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
2317          //          //
2318          glrun->Fill_GL_RUN(conn);          //      glrun->SetID(this->AssignRunID());
         //  
         // get id number  
         //  
         oss.str("");  
         oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "  
             << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "  
             << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";  
         //  
         if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());  
         //  
         result =  conn->Query(oss.str().c_str());  
         if ( !result ) throw -4;  
         row = result->Next();  
2319          //          //
2320          UInt_t idrun0 = 0;          glrun->SetID_RUN_FRAG(glrun1->GetID());
2321          if ( !row ){          glrun->Fill_GL_RUN(conn);
           throw -10;  
         } else {  
           idrun0 = (UInt_t)atoll(row->GetField(0));  
         };  
2322          //          //
2323          glrun1->SetID_RUN_FRAG(idrun0);          // set id number
2324          //          //
2325            glrun1->SetID_RUN_FRAG(glrun->GetID());
2326          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
2327          //          //
         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;  
         //  
2328        };        };
2329        //        //
2330        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 << ";";  
2331        //        //
2332        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");
2333        result =  conn->Query(oss.str().c_str());        glrun1->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
       if ( !result ) throw -4;  
2334        //        //
2335          delete glrun1;
2336        //        //
2337        return;        return;
2338        //        //
# Line 2051  void PamelaDBOperations::HandleRunFragme Line 2369  void PamelaDBOperations::HandleRunFragme
2369      //      //
2370      row = result->Next();      row = result->Next();
2371      //      //
2372      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 {  
2373        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");
2374        } else {
2375          if ( NoFrag() ){
2376            glrun->SetID_RUN_FRAG(glrun->GetID());
2377            glrun->Fill_GL_RUN(conn);      
2378            glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS");
2379          };    
2380      };      };
2381    };    };
2382    //    //
# Line 2115  void PamelaDBOperations::HandleMissingHo Line 2402  void PamelaDBOperations::HandleMissingHo
2402      //      //
2403      this->FillClass(mishead,mistrail,firstev,lastev);      this->FillClass(mishead,mistrail,firstev,lastev);
2404      //      //
2405      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2406          glrun->SetID(this->AssignRunID());
2407          glrun->SetID_RUN_FRAG(0);
2408          glrun->Fill_GL_RUN(conn);    
2409        };
2410      //      //
2411    };    };
2412    //    //
# Line 2282  Bool_t PamelaDBOperations::IsRunConsiste Line 2573  Bool_t PamelaDBOperations::IsRunConsiste
2573            //                  //      
2574            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2575            //            //
2576            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2577                glrun->SetID(this->AssignRunID());
2578                glrun->SetID_RUN_FRAG(0);
2579                glrun->Fill_GL_RUN(conn);      
2580              };
2581            //            //
2582            firstevno = lastentry + 1;            firstevno = lastentry + 1;
2583            //            //
# Line 2358  void PamelaDBOperations::HandleSuspiciou Line 2653  void PamelaDBOperations::HandleSuspiciou
2653      if ( IsDebug() ) printf(" Checking but no events in the run! \n");      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
2654      //      //
2655      this->FillClass();      this->FillClass();
2656      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2657          glrun->SetID(this->AssignRunID());
2658          glrun->SetID_RUN_FRAG(0);
2659          glrun->Fill_GL_RUN(conn);    
2660        };
2661      //      //
2662    } else {    } else {
2663      //      //
# Line 2386  void PamelaDBOperations::HandleSuspiciou Line 2685  void PamelaDBOperations::HandleSuspiciou
2685        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");
2686        //        //
2687        this->FillClass();        this->FillClass();
2688        if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                if ( !IsRunAlreadyInserted() ){
2689            glrun->SetID(this->AssignRunID());
2690            glrun->SetID_RUN_FRAG(0);
2691            glrun->Fill_GL_RUN(conn);        
2692          };
2693        //        //
2694      } else {      } else {
2695        //        //
# Line 2484  void PamelaDBOperations::HandleSuspiciou Line 2787  void PamelaDBOperations::HandleSuspiciou
2787            //                  //      
2788            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2789            //            //
2790            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2791                glrun->SetID(this->AssignRunID());
2792                glrun->SetID_RUN_FRAG(0);
2793                glrun->Fill_GL_RUN(conn);      
2794              };
2795            //            //
2796            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...
2797              //              //
# Line 2519  void PamelaDBOperations::HandleSuspiciou Line 2826  void PamelaDBOperations::HandleSuspiciou
2826              //                  //    
2827              this->SetCommonGLRUN(firstTime,lastTime);              this->SetCommonGLRUN(firstTime,lastTime);
2828              //              //
2829              if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                    if ( !IsRunAlreadyInserted() ){
2830                  glrun->SetID(this->AssignRunID());
2831                  glrun->SetID_RUN_FRAG(0);
2832                  glrun->Fill_GL_RUN(conn);      
2833                };
2834            };            };
2835            //            //
2836            firstevno = lastentry + 1;            firstevno = lastentry + 1;
# Line 2582  Int_t PamelaDBOperations::insertCALO_CAL Line 2893  Int_t PamelaDBOperations::insertCALO_CAL
2893          fromtime = this->GetAbsTime(ph->GetOrbitalTime());            fromtime = this->GetAbsTime(ph->GetOrbitalTime());  
2894          if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){          if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
2895            //            //
2896            if ( IsDebug() ) printf(" Calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);            if ( IsDebug() ) printf(" Calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
2897            //            //
2898            // check if the calibration has already been inserted            // check if the calibration has already been inserted
2899            //            //
# Line 2626  Int_t PamelaDBOperations::insertCALO_CAL Line 2937  Int_t PamelaDBOperations::insertCALO_CAL
2937                // no calibrations in the db contain our calibration                // no calibrations in the db contain our calibration
2938                //                //
2939                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);
2940                if ( fromtime < 1150863400 ){                if ( fromtime < 1150871000 ){ //1150866904
2941                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
2942                  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
2943                };                };
2944                //                //
# Line 2689  Int_t PamelaDBOperations::insertCALO_CAL Line 3000  Int_t PamelaDBOperations::insertCALO_CAL
3000            //            //
3001          } else {          } else {
3002            //            //
3003            if ( IsDebug() ) printf(" Repetead calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);            if ( IsDebug() ) printf(" Repeated calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
3004            //            //
3005          };                };      
3006          //                //      
# Line 2771  void PamelaDBOperations::HandleTRK_CALIB Line 3082  void PamelaDBOperations::HandleTRK_CALIB
3082        // no calibrations in the db contain our calibration        // no calibrations in the db contain our calibration
3083        //        //
3084        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");
3085        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
3086        //        //
3087        oss.str("");        oss.str("");
3088        oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "        oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE "
# Line 2900  Int_t PamelaDBOperations::insertTRK_CALI Line 3211  Int_t PamelaDBOperations::insertTRK_CALI
3211      pkt1 = ph1->GetCounter();        pkt1 = ph1->GetCounter();  
3212      fromtime = this->GetAbsTime(ph1->GetOrbitalTime());        fromtime = this->GetAbsTime(ph1->GetOrbitalTime());  
3213      //      //
3214      valid = 1;  //     valid = 1;
3215      //  //     //
3216      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
3217      //      //
3218      //      //
3219      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) ){
3220        //        //
3221        if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);        if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
3222          //      
3223          valid = ValidateTrkCalib( caltrk1, eh1 );
3224          if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
3225        //        //
3226        // Do we have the second calibration packet?        // Do we have the second calibration packet?
3227        //        //
# Line 2924  Int_t PamelaDBOperations::insertTRK_CALI Line 3238  Int_t PamelaDBOperations::insertTRK_CALI
3238            obt2 = ph2->GetOrbitalTime();              obt2 = ph2->GetOrbitalTime();  
3239            pkt2 = ph2->GetCounter();              pkt2 = ph2->GetCounter();  
3240            //            //
3241            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
3242            //                  //      
3243          } else {          } else {
3244            //            //
# Line 2947  Int_t PamelaDBOperations::insertTRK_CALI Line 3261  Int_t PamelaDBOperations::insertTRK_CALI
3261          //          //
3262        };        };
3263        //        //
3264        if ( IsDebug() ) printf(" Found trk calibration2 at obt %i pkt %i t2 is %i \n",obt2,pkt2,t2);        if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2);
3265        //        //
3266        // The calibration is good        // The calibration is good
3267        //        //
3268        if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){        if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){
3269          //          //
3270          if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2);          if ( IsDebug() ) printf(" The trk calibration2 at obt %u pkt %u t2 is %u is good \n",obt2,pkt2,t2);
3271            //
3272            UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
3273            if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
3274            valid = valid & valid2;
3275          //          //
3276          // Handle good calib          // Handle good calib
3277          //          //
# Line 2963  Int_t PamelaDBOperations::insertTRK_CALI Line 3281  Int_t PamelaDBOperations::insertTRK_CALI
3281          //          //
3282          if ( t2 != pret2+1 ){          if ( t2 != pret2+1 ){
3283            //            //
3284            if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i pret2 is %i \n",obt2,pkt2,t2,pret2);            if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u pret2 is %u \n",obt2,pkt2,t2,pret2);
3285            //            //
3286            while ( t2 > pret2+1 ){            while ( t2 > pret2+1 ){
3287              //              //
# Line 2992  Int_t PamelaDBOperations::insertTRK_CALI Line 3310  Int_t PamelaDBOperations::insertTRK_CALI
3310          //          //
3311          // Check for missing calibtrk2          // Check for missing calibtrk2
3312          //          //
3313          if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);          if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
3314          t2 = pret2;          t2 = pret2;
3315          //          //
3316          // handle missing calib2          // handle missing calib2
# Line 3006  Int_t PamelaDBOperations::insertTRK_CALI Line 3324  Int_t PamelaDBOperations::insertTRK_CALI
3324        //        //
3325      } else {      } else {
3326        //        //
3327        if ( IsDebug() ) printf(" Repetead trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);        if ( IsDebug() ) printf(" Repetead trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
3328        //        //
3329      };      };
3330      //          //    
# Line 3031  Int_t PamelaDBOperations::insertTRK_CALI Line 3349  Int_t PamelaDBOperations::insertTRK_CALI
3349      valid = 0;      valid = 0;
3350      if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){      if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
3351        //        //
3352        if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);        if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
3353        //        //
3354        this->HandleTRK_CALIB(false,true);        this->HandleTRK_CALIB(false,true);
3355        //        //
# Line 3085  Int_t PamelaDBOperations::insertS4_CALIB Line 3403  Int_t PamelaDBOperations::insertS4_CALIB
3403      fromtime = this->GetAbsTime(ph->GetOrbitalTime());        fromtime = this->GetAbsTime(ph->GetOrbitalTime());  
3404      if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){      if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
3405        //        //
3406        if ( IsDebug() ) printf(" S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);        if ( IsDebug() ) printf(" S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
3407        //        //
3408        // check if the calibration has already been inserted        // check if the calibration has already been inserted
3409        //        //
# Line 3127  Int_t PamelaDBOperations::insertS4_CALIB Line 3445  Int_t PamelaDBOperations::insertS4_CALIB
3445            // no calibrations in the db contain our calibration            // no calibrations in the db contain our calibration
3446            //            //
3447            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");
3448            if ( fromtime < 1150863400 ){            if ( fromtime < 1150871000 ){  
3449              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
3450              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
3451            };            };
3452            //            //
# Line 3187  Int_t PamelaDBOperations::insertS4_CALIB Line 3505  Int_t PamelaDBOperations::insertS4_CALIB
3505        //        //
3506      } else {      } else {
3507        //        //
3508        if ( IsDebug() ) printf(" Repetead S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);        if ( IsDebug() ) printf(" Repeated S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
3509        //        //
3510      };            };      
3511      //            //      
# Line 3200  Int_t PamelaDBOperations::insertS4_CALIB Line 3518  Int_t PamelaDBOperations::insertS4_CALIB
3518   * 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
3519   */   */
3520  Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){  Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(){
3521      return(this->CleanGL_RUN_FRAGMENTS(""));
3522    };
3523    
3524    /**
3525     * Scan the fragment table and move old fragments to the GL_RUN table
3526     */
3527    Int_t PamelaDBOperations::CleanGL_RUN_FRAGMENTS(TString fcleanfile){
3528    //    //
3529    TSQLResult *result = 0;    TSQLResult *result = 0;
3530    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
# Line 3207  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3532  Int_t PamelaDBOperations::CleanGL_RUN_FR
3532    TSQLRow    *row2   = 0;    TSQLRow    *row2   = 0;
3533    //    //
3534    UInt_t moved = 0;    UInt_t moved = 0;
 //  UInt_t timelim = 0;  
 //  TDatime *time = new TDatime();  
3535    //    //
3536    stringstream oss;    stringstream oss;
3537    oss.str("");    oss.str("");
3538    //    //
3539    //    if ( !strcmp(fcleanfile.Data(),"") ){
3540    // check if there are entries older than "olderthan" seconds from now      //
3541    //      // check if there are entries older than "olderthan" seconds from now
3542    oss.str("");      //
3543    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"      oss.str("");
3544        << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";      oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
3545    //          << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";
3546    if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str());      //
3547    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());
3548        result = conn->Query(oss.str().c_str());
3549        //
3550      } else {
3551        oss.str("");
3552        oss << " SELECT ID FROM GL_ROOT WHERE NAME='" << fcleanfile.Data() << "';";
3553        if ( IsDebug() ) printf(" Getting ID_ROOT_L0 query %s \n",oss.str().c_str());
3554        result = conn->Query(oss.str().c_str());
3555        //    
3556        if ( result ){
3557          //
3558          row = result->Next();
3559          //
3560          if ( row ){
3561            oss.str("");
3562            oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
3563                << " ID_ROOT_L0=" << row->GetField(0) << ";";
3564            //
3565            if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS for ROOT file query is \n %s \n",oss.str().c_str());
3566            result = conn->Query(oss.str().c_str());
3567            //    
3568          };
3569        } else {
3570          return(2);
3571        };
3572      };
3573    //    //
3574    if ( result ){    if ( result ){
3575      //      //
3576      row = result->Next();      row = result->Next();
3577      //      //
3578      if ( row ){      while ( row ){
3579        //        //
3580        oss.str("");        oss.str("");
3581        oss << " ID= "<< row->GetField(0);        oss << " ID= "<< row->GetField(0);
# Line 3262  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3610  Int_t PamelaDBOperations::CleanGL_RUN_FR
3610          if ( IsDebug() ) printf(" The run is new \n");          if ( IsDebug() ) printf(" The run is new \n");
3611          if ( IsDebug() ) printf(" -> fill the DB \n");                if ( IsDebug() ) printf(" -> fill the DB \n");      
3612          //          //
3613          glrun->SetID(0);          //      glrun->SetID(this->AssignRunID()); we use the old run number!
3614            glrun->SetID_RUN_FRAG(glrun->GetID());
3615          glrun->Fill_GL_RUN(conn);            glrun->Fill_GL_RUN(conn);  
3616          //          //
3617          oss.str("");  //      oss.str("");
3618          oss << " SELECT ID FROM GL_RUN WHERE "  //      oss << " SELECT ID FROM GL_RUN WHERE "
3619              << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "  //          << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "
3620              << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "  //          << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "
3621              << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "  //          << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "
3622              << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "  //          << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "
3623              << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";  //          << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";
3624          //  //      //
3625          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());
3626          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3627          //  //      //
3628          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3629          //  //      //
3630          row2 = result2->Next();  //      row2 = result2->Next();
3631          //  //      //
3632          if ( !row2 ) throw -25;  //      if ( !row2 ) throw -25;
3633          //  //      //
3634          oss.str("");  //      oss.str("");
3635          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);
3636          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());
3637          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3638          //  //      //
3639          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3640          //          //
3641          moved++;          moved++;
3642          //          //
# Line 3297  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3646  Int_t PamelaDBOperations::CleanGL_RUN_FR
3646        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));      
3647        //        //
3648        //        //
3649        oss.str("");        glrun->DeleteRun(conn,(UInt_t)atoll(row->GetField(0)),"GL_RUN_FRAGMENTS");
3650        oss << " DELETE from GL_RUN_FRAGMENTS where ID = " << row->GetField(0);  //      oss.str("");
3651        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);
3652        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());
3653        //  //       result2 = conn->Query(oss.str().c_str());
3654        if ( !result2 ) throw -4;  //      //
3655        //  //      if ( !result2 ) throw -4;
3656    //      //
3657          row = result->Next();
3658      };      };
3659    };    };
3660    if ( IsDebug() ) printf(" Moved %u runs\n",moved);    if ( IsDebug() ) printf(" Moved %u runs\n",moved);
# Line 3314  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3665  Int_t PamelaDBOperations::CleanGL_RUN_FR
3665   * 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..
3666   */   */
3667  Int_t PamelaDBOperations::ValidateRuns(){  Int_t PamelaDBOperations::ValidateRuns(){
3668      return(this->ValidateRuns(""));
3669    };
3670    
3671    /**
3672     * Check if runs are good, i.e. if the tracker calibration is correctly associated..
3673     */
3674    Int_t PamelaDBOperations::ValidateRuns(TString valfile){
3675    //    //
3676    TSQLResult *result = 0;    TSQLResult *result = 0;
3677    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
3678    //    //
3679      UInt_t calibtime = 50;
3680      //
3681    stringstream oss;    stringstream oss;
3682    oss.str("");    oss.str("");
3683    //    //
# Line 3326  Int_t PamelaDBOperations::ValidateRuns() Line 3686  Int_t PamelaDBOperations::ValidateRuns()
3686    // =======================================================    // =======================================================
3687    UInt_t t_stop  = 0;    UInt_t t_stop  = 0;
3688    UInt_t t_start = 0;    UInt_t t_start = 0;
3689    // --------------------------------------------------------------    if ( !strcmp(valfile.Data(),"") ) {
3690    // 1) get the OBT of the last run inserted after clean-time limit      // --------------------------------------------------------------
3691    // --------------------------------------------------------------      // 1) get the OBT of the last run inserted after clean-time limit
3692    oss.str("");      // --------------------------------------------------------------
3693    oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()      oss.str("");
3694            << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";      oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()
3695    if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());          << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3696    result = conn->Query(oss.str().c_str());      if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());
3697    if ( !result ) throw -4;      result = conn->Query(oss.str().c_str());
3698    if ( !result->GetRowCount() ) {      if ( !result ) throw -4;
3699            printf(" No runs to validate \n");      if ( !result->GetRowCount() ) {
3700            return(1);        printf(" No runs to validate \n");
3701    }else{        return(1);
3702          row = result->Next();      }else{
3703          t_start = (UInt_t)atoll(row->GetField(4));        row = result->Next();
3704    };          t_start = (UInt_t)atoll(row->GetField(4));
3705    // --------------------------------------------------------------      };  
3706    // 2) get the OBT of the last validated run      // --------------------------------------------------------------
3707    // --------------------------------------------------------------      // 2) get the OBT of the last validated run
3708    oss.str("");      // --------------------------------------------------------------
3709    oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start      oss.str("");
3710        <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";      oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
3711    if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());          <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3712    result = conn->Query(oss.str().c_str());      if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
3713    if ( !result ) throw -4;      result = conn->Query(oss.str().c_str());
3714    if ( result->GetRowCount() ){      if ( !result ) throw -4;
3715            row = result->Next();      if ( result->GetRowCount() ){
3716            t_stop = (UInt_t)atoll(row->GetField(4));        row = result->Next();
3717          t_stop = (UInt_t)atoll(row->GetField(4));
3718        };
3719        if ( IsDebug() ) printf("Validation interval: from time %u - to time %u \n\n",t_stop,t_start);
3720        // --------------------------------------------------------------
3721        // now retrieves runs to be validated
3722        // --------------------------------------------------------------
3723        oss.str("");
3724        oss << " SELECT * FROM GL_RUN  WHERE  RUNHEADER_TIME <=" << t_start;
3725        oss << " AND RUNHEADER_TIME >="<< t_stop;
3726        oss << " ORDER BY RUNHEADER_TIME DESC;";
3727        if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());
3728        result = conn->Query(oss.str().c_str());
3729      } else {
3730        //
3731        stringstream myquery;
3732        UInt_t myid = 0;
3733        myquery.str("");
3734        myquery << " SELECT ID FROM GL_ROOT where NAME='"<<valfile.Data() <<"';";
3735        //
3736        result = conn->Query(myquery.str().c_str());
3737        //
3738        row = result->Next();      
3739        if( !row ){
3740          if ( strcmp(valfile.Data(),GetRootName().Data()) ){
3741            if ( IsDebug() ) printf(" No file to be validated even if option \"-validate file\" was used!!\n");
3742            return(2);
3743          };
3744          if ( IsDebug() ) printf(" No file to be validated (force mode)! \n");
3745          return(0);
3746        };
3747        myid=(UInt_t)atoll(row->GetField(0));
3748        //
3749        myquery.str("");
3750        myquery << " SELECT MAX(RUNTRAILER_TIME),MIN(RUNHEADER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< myid <<";";
3751        //
3752        result = conn->Query(myquery.str().c_str());
3753        //
3754        row = result->Next();      
3755        if( !row->GetField(0) || !row->GetField(1)){
3756          //
3757          if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
3758          //
3759          return(0);
3760          //
3761        } else {  
3762          //
3763          UInt_t runhtime = (UInt_t)atoll(row->GetField(0));
3764          UInt_t runttime = (UInt_t)atoll(row->GetField(1));
3765          UInt_t caltime = 0;
3766          //
3767          myquery.str("");
3768          myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
3769          myquery << " order by FROM_TIME asc limit 1;";
3770          //
3771          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
3772          //
3773          //
3774          result = conn->Query(myquery.str().c_str());
3775          //
3776          row = result->Next();      
3777          if( !row ){
3778            caltime = runhtime;
3779          } else {
3780            caltime = (UInt_t)atoll(row->GetField(0));
3781          };
3782          //
3783          myquery.str("");
3784          myquery << " SELECT * from GL_RUN where RUNHEADER_TIME>="<< runttime <<" AND RUNHEADER_TIME<=" ;
3785          myquery << caltime << " order by RUNHEADER_TIME DESC";
3786          //
3787          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
3788          //
3789          result = conn->Query(myquery.str().c_str());
3790          //
3791        };
3792    };    };
3793    if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);    //
   // --------------------------------------------------------------  
   // now retrieves runs to be validated  
   // --------------------------------------------------------------  
   oss.str("");  
   oss << " SELECT * FROM GL_RUN  WHERE  RUNHEADER_TIME <=" << t_start;  
   oss << " AND RUNHEADER_TIME >="<< t_stop;  
   oss << " ORDER BY RUNHEADER_TIME DESC;";  
 //  if ( IsDebug() )  
   if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());  
   result = conn->Query(oss.str().c_str());  
3794    if ( !result ) throw -4;    if ( !result ) throw -4;
3795    if ( !result->GetRowCount() ) printf(" No runs to validate \n");    if ( !result->GetRowCount() && IsDebug() ) printf(" No runs to validate \n");
3796  //  printf("------------------------------------------------------------------------------- \n");    //
     
3797    Int_t nrow = 0;    Int_t nrow = 0;
3798    GL_RUN* this_run = new GL_RUN();    GL_RUN* this_run = new GL_RUN();
3799    GL_RUN* next_run = new GL_RUN();    GL_RUN* next_run = new GL_RUN();
# Line 3383  Int_t PamelaDBOperations::ValidateRuns() Line 3807  Int_t PamelaDBOperations::ValidateRuns()
3807    UInt_t t1=0,t2=0;    UInt_t t1=0,t2=0;
3808    // ---------------------------------------------------------------------------------    // ---------------------------------------------------------------------------------
3809    // - loop over runs, back in time,    // - loop over runs, back in time,
3810    // - select sequences of runs close in time (less than 60 s apart),    // - select sequences of runs close in time (less than calibtime s apart),
3811    //   which could be preceeded by a calibration    //   which could be preceeded by a calibration
3812    // - check if there might be a missing calibration    // - check if there might be a missing calibration
3813    // ---------------------------------------------------------------------------------    // ---------------------------------------------------------------------------------
# Line 3423  Int_t PamelaDBOperations::ValidateRuns() Line 3847  Int_t PamelaDBOperations::ValidateRuns()
3847                                                        
3848                            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
3849                                                        
3850                            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
3851                                                        
3852                            if( !CHECK && this_run->VALIDATION ){                            if( !CHECK && this_run->VALIDATION ){
3853                                    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 3463  Int_t PamelaDBOperations::ValidateRuns() Line 3887  Int_t PamelaDBOperations::ValidateRuns()
3887                    }else printf("ValidateRuns ***WARNING*** : run sequence exceed assumed size (%i) \n",nseq_max);                    }else printf("ValidateRuns ***WARNING*** : run sequence exceed assumed size (%i) \n",nseq_max);
3888            };            };
3889                        
3890            if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID);            if ( IsDebug() ) printf("%i Run %u \n",nrow,this_run->ID);
3891            nrow++;            nrow++;
3892                        
3893    };    };
# Line 3492  Bool_t PamelaDBOperations::MissingTRK_CA Line 3916  Bool_t PamelaDBOperations::MissingTRK_CA
3916          // which should be equal to the time between ascending-nodes.          // which should be equal to the time between ascending-nodes.
3917          //==============================================================          //==============================================================
3918          if ( t2 - trkcalib->FROM_TIME > 5700) {          if ( t2 - trkcalib->FROM_TIME > 5700) {
3919                  if ( IsDebug() )printf("Long time between calib and run start %i :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);                  if ( IsDebug() )printf("Long time between calib and run start %u :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);
3920          //==============================================================          //==============================================================
3921          // there might be a missing calibration, due to:          // there might be a missing calibration, due to:
3922          // - MM full          // - MM full
# Line 3513  Bool_t PamelaDBOperations::MissingTRK_CA Line 3937  Bool_t PamelaDBOperations::MissingTRK_CA
3937          // it is enough to say that there are no missing calibrations          // it is enough to say that there are no missing calibrations
3938          //==============================================================          //==============================================================
3939          // the long time interval bewteen runs might be due to download          // the long time interval bewteen runs might be due to download
3940          if ( IsDebug() )printf("Short time between calib and run start %i :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);          if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
3941          return(false);          return(false);
3942                    
3943  };  };
# Line 3538  Int_t PamelaDBOperations::assignVALIDATI Line 3962  Int_t PamelaDBOperations::assignVALIDATI
3962    
3963    
3964    
3965    // Insert TLEs from file tlefilename in the table GL_TLE in the db
3966    // opened by conn, sorting them by date from older to newer, if each
3967    // TLE has not been alread inserted.
3968    Int_t PamelaDBOperations::populateTLE()//(TSQLServer *conn, char *tleFile)
3969    {
3970      fstream tlefile(tlefilename, ios::in);
3971    
3972      if ( !tlefile ) throw -7;
3973    
3974      vector<cTle*> ctles;
3975      vector<cTle*>::iterator iter;
3976      int present = 0;
3977    
3978      // Get three lines from tlefile, create a cTle object and put it
3979      // into ctles
3980      while(1) {
3981        cTle *tlef;
3982        string str1, str2, str3;
3983    
3984        getline(tlefile, str1);
3985        if(tlefile.eof()) break;
3986    
3987        getline(tlefile, str2);
3988        if(tlefile.eof()) break;
3989    
3990        getline(tlefile, str3);
3991        if(tlefile.eof()) break;
3992    
3993        // We now have three good lines for a cTle.
3994        tlef = new cTle(str1, str2, str3);
3995        ctles.push_back(tlef);
3996      }
3997    
3998      tlefile.close();
3999    
4000      // Sort by date
4001      sort(ctles.begin(), ctles.end(), compTLE);
4002    
4003      // Now we insert each TLE into the db
4004      for(iter = ctles.begin(); iter != ctles.end(); iter++) {
4005        cTle *tle = *iter;
4006    
4007        // Do nothing if it's already present in the db.  Just increase
4008        // the counter present.
4009        if (! isTlePresent(tle))
4010          {
4011            int status = insertTle(tle);
4012    
4013            // Insert query failed.  Return 1.
4014            if(status == EXIT_FAILURE) {
4015              
4016              if( IsDebug() ) {
4017                cerr << "Error: inserting TLE:" << endl
4018                     << tle->getName() << endl
4019                     << tle->getLine1() << endl
4020                     << tle->getLine2() << endl;
4021              }
4022    
4023              throw -4;
4024              return 1;
4025            }
4026    
4027          }
4028        else
4029          present++;
4030    
4031      }
4032    
4033      int inserted = ctles.size() - present;  // Number of inserted TLE.
4034      if ( IsDebug() )
4035        cout << "\nProcessed TLEs ranging from " << getTleDatetime(ctles[0]) << " to " << getTleDatetime(ctles[ctles.size()-1]) << "." << endl
4036             << inserted << " newly inserted TLEs out of " << ctles.size() << " processed." << endl;
4037    
4038      ctles.clear();
4039    
4040    
4041      // Return 2 if no new TLE has been inserted.  0 otherwise.
4042      if(! inserted ) return 2;
4043      return 0;
4044    }
4045    
4046    
4047    // Insert tle in the table GL_TLE using the connection conn.
4048    Int_t PamelaDBOperations::insertTle(cTle *tle)
4049    {
4050      stringstream oss;
4051      TSQLResult *result = 0;
4052    
4053      oss.str("");
4054      oss << " INSERT INTO GL_TLE (TLE1, TLE2, TLE3, FROM_TIME)"
4055          << " VALUES ( '"
4056          << tle->getName() << "', '"
4057          << tle->getLine1() << "', '"
4058          << tle->getLine2() << "', '"
4059          << getTleDatetime(tle) << "')";
4060    
4061      //  cout << oss.str().c_str() << endl;
4062      result = conn->Query(oss.str().c_str());
4063      if (result == NULL)
4064        return EXIT_FAILURE;
4065    
4066      return EXIT_SUCCESS;
4067    }
4068    
4069    
4070    // Return whether tle is already in the db connected by conn.
4071    bool PamelaDBOperations::isTlePresent(cTle *tle)
4072    {
4073      stringstream oss;
4074      TSQLResult *result = 0;
4075    
4076      oss.str("");
4077      oss << "SELECT * FROM GL_TLE WHERE FROM_TIME = '"
4078          << getTleDatetime(tle) << "'";
4079    
4080      result = conn->Query(oss.str().c_str());
4081      if (result == NULL) throw -4;
4082    
4083      if (result->GetRowCount())
4084        return true;
4085      else
4086        return false;
4087    }
4088    
4089    
4090    // Return whether the first TLE is dated early than the second
4091    bool compTLE (cTle *tle1, cTle *tle2)
4092    {
4093      return getTleJulian(tle1) < getTleJulian(tle2);
4094    }
4095    
4096    
4097    // Return the date of the tle using the format (year-2000)*1e3 +
4098    // julian day.  e.g. 6365 is the 31th Dec 2006.
4099    // It does *not* return a cJulian date.
4100    float getTleJulian(cTle *tle) {
4101      return tle->getField(cTle::FLD_EPOCHYEAR)*1e3 + tle->getField(cTle::FLD_EPOCHDAY);
4102    }
4103    
4104    
4105    // Return a string like YYYY-MM-DD hh:mm:ss, usable for mysql datetime
4106    // format.
4107    string getTleDatetime(cTle *tle)
4108    {
4109      int year, mon, day, hh, mm, ss;
4110      double dom; // day of month (is double!)
4111      stringstream date; // date in datetime format
4112    
4113      // create a cJulian from the date in tle
4114      cJulian jdate = cJulian( 2000 + (int) tle->getField(cTle::FLD_EPOCHYEAR), tle->getField(cTle::FLD_EPOCHDAY));
4115    
4116      // get year, month, day of month
4117      jdate.getComponent(&year, &mon, &dom);
4118    
4119      // build a datetime YYYY-MM-DD hh:mm:ss
4120      date.str("");
4121      day = (int) floor(dom);
4122      hh = (int) floor( (dom - day) * 24);
4123      mm = (int) floor( ((dom - day) * 24 - hh) * 60);
4124      ss = (int) floor( ((((dom - day) * 24 - hh) * 60 - mm) * 60));
4125      //  ms = (int) floor( (((((dom - day) * 24 - hh) * 60 - mm) * 60) - ss) * 1000);
4126    
4127      date << year << "-" << mon << "-" << day << " " << hh << ":" << mm << ":" << ss;
4128    
4129      return date.str();
4130    }
4131    
4132    /**
4133     * Remove a file from the DB, delete on cascade all entries related to that file
4134     * rearrange GL_RUN and GL_XXX_CALIB tables, turn off validation till the following
4135     * calibration
4136     **/
4137    Int_t PamelaDBOperations::removeFile(TString remfile){
4138      //
4139      // Determine ID_ROOT_L0 and ID_RAW
4140      //
4141      TSQLResult *pResult;
4142      TSQLRow *Row;
4143      stringstream myquery;
4144      //  
4145      myquery.str("");
4146      myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';";
4147      //
4148      pResult = conn->Query(myquery.str().c_str());
4149      //
4150      Row = pResult->Next();      
4151      if( !Row ){
4152        if ( strcmp(remfile.Data(),GetRootName().Data()) ){
4153          if ( IsDebug() ) printf(" No file to be removed even if option \"-remove file\" was used!!\n");
4154          return(1);
4155        };
4156        if ( IsDebug() ) printf(" No file to be removed (force mode)! \n");
4157        return(0);
4158      };
4159      //
4160      this->SetID_ROOT((UInt_t)atoll(Row->GetField(0)));
4161      this->SetID_RAW((UInt_t)atoll(Row->GetField(1)));
4162      //
4163      this->ValidationOFF();
4164      //
4165      this->RemoveCALIBS();
4166      //
4167      this->RemoveRUNS();
4168      //
4169      this->RemoveFILES();
4170      //
4171      this->SetID_ROOT(0);
4172      this->SetID_RAW(0);
4173      //
4174      return(0);
4175    };
4176    
4177    /**
4178     *
4179     * Set validation bit to zero for runs following the removing file till
4180     * 1) a run with TRK_CALIB_USED=140
4181     * 2) a run with VALIDATION = 0
4182     * 3) the next calibration
4183     *
4184     **/
4185    void PamelaDBOperations::ValidationOFF(){
4186      TSQLResult *pResult;
4187      TSQLRow *Row;
4188      stringstream myquery;
4189      Int_t unv = 0;
4190      //select ID from GL_RUN where RUNHEADER_TIME>=1152671382 AND (VALIDATION=0 OR TRK_CALIB_USED=104) order by RUNHEADER_TIME asc limit 1;
4191      myquery.str("");
4192      myquery << " SELECT MAX(RUNTRAILER_TIME) FROM GL_RUN WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4193      //
4194      pResult = conn->Query(myquery.str().c_str());
4195      //
4196      Row = pResult->Next();      
4197      if( !Row->GetField(0) ){
4198        //
4199        if ( IsDebug() ) printf(" NO RUN ASSOCIATED TO THIS FILE! \n");
4200        //
4201      } else {  
4202        //
4203        UInt_t runhtime = (UInt_t)atoll(Row->GetField(0));
4204        UInt_t caltime = 0;
4205        //
4206        myquery.str("");
4207        myquery << " SELECT FROM_TIME FROM GL_TRK_CALIB where FROM_TIME>" <<runhtime;
4208        myquery << " order by FROM_TIME asc limit 1;";
4209        //
4210        if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4211        //
4212        //
4213        delete pResult;
4214        pResult = conn->Query(myquery.str().c_str());
4215        //
4216        Row = pResult->Next();      
4217        if( !Row ){
4218          caltime = runhtime;
4219        } else {
4220          caltime = (UInt_t)atoll(Row->GetField(0));
4221        };
4222        //
4223        myquery.str("");
4224        myquery << " SELECT ID,RUNHEADER_TIME from GL_RUN where RUNHEADER_TIME>="<< runhtime <<" AND (VALIDATION=0 OR TRK_CALIB_USED=104 OR RUNHEADER_TIME>" ;
4225        myquery << caltime << ") order by RUNHEADER_TIME asc LIMIT 1";
4226        //
4227        if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4228        //
4229        pResult = conn->Query(myquery.str().c_str());
4230        //
4231        Row = pResult->Next();      
4232        if( !Row ){
4233          //
4234          if ( IsDebug() ) printf(" NO RUN NEED TO BE UNVALIDATED \n");
4235          //
4236        } else {
4237          myquery.str("");
4238          myquery << " SELECT ID from GL_RUN where RUNHEADER_TIME<"<< Row->GetField(1) <<" AND ";
4239          myquery << " RUNHEADER_TIME>=" <<runhtime;
4240          myquery << " order by RUNHEADER_TIME asc;";
4241          //
4242          if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4243          //
4244          pResult = conn->Query(myquery.str().c_str());
4245          //
4246          Row = pResult->Next();  
4247          while ( Row ){
4248            //      
4249            unv++;
4250            this->assignVALIDATION((UInt_t)atoll(Row->GetField(0)), false);
4251            Row = pResult->Next();  
4252            //
4253          };
4254        };
4255      };
4256      if ( IsDebug() ) printf(" %u runs have been unvalidated \n",unv);
4257    };
4258    
4259    /**
4260     *
4261     * Rearrange GL_RUN table and remove runs
4262     *
4263     **/
4264    void PamelaDBOperations::RemoveRUNS(){
4265      TSQLResult *pResult;
4266      TSQLRow *Row;
4267      stringstream myquery;
4268      UInt_t drun = 0;
4269      GL_RUN *delrun = new GL_RUN();
4270      //
4271      myquery.str("");
4272      myquery << " SELECT ID FROM GL_RUN where ID_RUN_FRAG=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4273      //
4274      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4275      //
4276      pResult = conn->Query(myquery.str().c_str());
4277      //
4278      Row = pResult->Next();    
4279      //
4280      //
4281      if ( !Row ){
4282        if ( IsDebug() ) printf(" No run with ID_RUN_FRAG=0 belonged to this file \n");
4283      } else {
4284        if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");    
4285        while ( Row ){
4286          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4287          if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4288          drun++;
4289          Row = pResult->Next();    
4290        };
4291      };
4292      //
4293      //
4294      myquery.str("");
4295      myquery << " SELECT ID,ID_RUN_FRAG FROM GL_RUN where ID_RUN_FRAG!=0 and ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4296      //
4297      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4298      //
4299      pResult = conn->Query(myquery.str().c_str());
4300      //
4301      Row = pResult->Next();    
4302      //
4303      if ( !Row ){
4304        if ( IsDebug() ) printf(" No run with ID_RUN_FRAG!=0 belonged to this file \n");
4305      } else {
4306        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");    
4307        while ( Row ){
4308          if ( IsDebug() ) printf(" restore run %u \n",(UInt_t)atoll(Row->GetField(1)));
4309          delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");
4310          if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(1)));
4311          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");
4312          if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){
4313            if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4314            delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4315          };
4316          drun++;
4317          Row = pResult->Next();    
4318        };
4319      };
4320      //
4321      if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN table \n",drun);    
4322      //
4323      //
4324      //
4325      drun = 0;
4326      //
4327      myquery.str("");
4328      myquery << " SELECT ID_TRASH FROM GL_RUN_TRASH where BELONGED_TO='GL_RUN_FRAGMENTS' AND ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4329      //
4330      pResult = conn->Query(myquery.str().c_str());
4331      //
4332      Row = pResult->Next();    
4333      //
4334      if ( !Row ){
4335        if ( IsDebug() ) printf(" No run from GL_RUN_FRAGMENTS table in the trash table for this file \n");
4336      } else {
4337        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");    
4338        while ( Row ){
4339          if ( IsDebug() ) printf(" del run idtrash %u \n",(UInt_t)atoll(Row->GetField(0)));
4340          myquery.str("");
4341          myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";      
4342          conn->Query(myquery.str().c_str());
4343          drun++;
4344          Row = pResult->Next();    
4345        };
4346      };
4347      //  
4348      if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_TRASH table \n",drun);    
4349      //
4350      //
4351      //
4352      drun = 0;
4353      //
4354      myquery.str("");
4355      myquery << " SELECT ID FROM GL_RUN_FRAGMENTS where ID_ROOT_L0=" <<this->GetID_ROOT() <<";";
4356      //
4357      pResult = conn->Query(myquery.str().c_str());
4358      //
4359      Row = pResult->Next();    
4360      //
4361      if ( !Row ){
4362        if ( IsDebug() ) printf(" No run in the GL_RUN_FRAGMENTS table for this file \n");
4363      } else {
4364        if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");    
4365        while ( Row ){
4366          if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4367          myquery.str("");
4368          myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";
4369          conn->Query(myquery.str().c_str());
4370          drun++;
4371          Row = pResult->Next();    
4372        };
4373      };
4374      //  
4375      if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_FRAGMENTS table \n",drun);    
4376      //
4377      //
4378      //
4379      delete delrun;
4380      //
4381    };
4382    
4383    
4384    /**
4385     *
4386     * Rearrange calibration tables
4387     *
4388     **/
4389    void PamelaDBOperations::RemoveFILES(){
4390      stringstream myquery;
4391      //
4392      myquery.str("");
4393      myquery << " DELETE FROM GL_RAW WHERE ID=" <<this->GetID_RAW() <<";";
4394      //
4395      if ( IsDebug() ) printf("  query is \n %s \n",myquery.str().c_str());
4396      //
4397      conn->Query(myquery.str().c_str());
4398      //
4399    };
4400    
4401    /**
4402     *
4403     * Rearrange calibration tables
4404     *
4405     **/
4406    void PamelaDBOperations::RemoveCALIBS(){
4407      TSQLResult *pResult;
4408      TSQLRow *Row;
4409      stringstream myquery;
4410      //  
4411      //
4412      // Calorimeter
4413      //
4414      for (Int_t section = 0; section < 4; section++){
4415        myquery.str("");
4416        myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_CALO_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<" AND ";
4417        myquery << " SECTION=" << section << ";";
4418        //
4419        pResult = conn->Query(myquery.str().c_str());
4420        //
4421        Row = pResult->Next();      
4422        if( !Row->GetField(0) || !Row->GetField(1) ){
4423          //
4424          if ( IsDebug() ) printf(" NO CALO CALIBRATION SECTION %i ASSOCIATED TO THIS FILE! \n",section);
4425          //
4426        } else {
4427          //
4428          myquery.str("");
4429          myquery << " UPDATE GL_CALO_CALIB SET TO_TIME=" << Row->GetField(1);
4430          myquery << " WHERE TO_TIME="<< Row->GetField(0) << " AND ";
4431          myquery << " SECTION=" << section << ";";
4432          //
4433          pResult = conn->Query(myquery.str().c_str());
4434          //
4435          if( !pResult ){
4436            //
4437            if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
4438            //
4439            throw -4;
4440            //
4441          };
4442          //
4443        };
4444      };
4445      myquery.str("");
4446      myquery << " DELETE FROM GL_CALO_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4447      //
4448      pResult = conn->Query(myquery.str().c_str());
4449      //
4450      if( !pResult ){
4451        //
4452        if ( IsDebug() ) printf(" ERROR DELETING CALO CALIBRATIONS \n");
4453        //
4454        throw -4;
4455        //
4456      };
4457      //
4458      // Tracker
4459      //
4460      myquery.str("");
4461      myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_TRK_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4462      //
4463      pResult = conn->Query(myquery.str().c_str());
4464      //
4465      Row = pResult->Next();      
4466      if( !Row->GetField(0) || !Row->GetField(1) ){
4467        //
4468        if ( IsDebug() ) printf(" NO TRK CALIBRATION ASSOCIATED TO THIS FILE! \n");
4469        //
4470      } else {
4471        //
4472        myquery.str("");
4473        myquery << " UPDATE GL_TRK_CALIB SET TO_TIME=" << Row->GetField(1);
4474        myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
4475        //
4476        pResult = conn->Query(myquery.str().c_str());
4477        //
4478        if( !pResult ){
4479          //
4480          if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
4481          //
4482          throw -4;
4483          //
4484        };
4485        //
4486        myquery.str("");
4487        myquery << " DELETE FROM GL_TRK_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4488        //
4489        pResult = conn->Query(myquery.str().c_str());
4490        //
4491        if( !pResult ){
4492          //
4493          if ( IsDebug() ) printf(" ERROR DELETING TRK CALIBRATIONS \n");
4494          //
4495          throw -4;
4496          //
4497        };
4498      };
4499      //
4500      //
4501      // S4
4502      //
4503      myquery.str("");
4504      myquery << " SELECT MIN(FROM_TIME),MAX(TO_TIME) FROM GL_S4_CALIB WHERE ID_ROOT_L0="<< this->GetID_ROOT() <<";";
4505      //
4506      pResult = conn->Query(myquery.str().c_str());
4507      //
4508      Row = pResult->Next();      
4509      if( !Row->GetField(0) || !Row->GetField(1) ){
4510        //
4511        if ( IsDebug() ) printf(" NO S4 CALIBRATION ASSOCIATED TO THIS FILE! \n");
4512        //
4513      } else {
4514        //
4515        myquery.str("");
4516        myquery << " UPDATE GL_S4_CALIB SET TO_TIME=" << Row->GetField(1);
4517        myquery << " WHERE TO_TIME="<< Row->GetField(0) << ";";
4518        //
4519        pResult = conn->Query(myquery.str().c_str());
4520        //
4521        if( !pResult ){
4522          //
4523          if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
4524          //
4525          throw -4;
4526          //
4527        };
4528        //
4529        myquery.str("");
4530        myquery << " DELETE FROM GL_S4_CALIB WHERE ID_ROOT_L0=" << this->GetID_ROOT();
4531        //
4532        pResult = conn->Query(myquery.str().c_str());
4533        //
4534        if( !pResult ){
4535          //
4536          if ( IsDebug() ) printf(" ERROR DELETING S4 CALIBRATIONS \n");
4537          //
4538          throw -4;
4539          //
4540        };
4541        //
4542      };
4543    };
4544    
4545    /**
4546     *
4547     * Rearrange calibration tables
4548     *
4549     **/
4550    UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){
4551    
4552        Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0};
4553        UInt_t timeaftercalib=120000; //2000;
4554    //  ----------
4555    //  Check CRCs
4556    //  ----------
4557        for(Int_t ipkt=0; ipkt<6; ipkt++){
4558            if( caltrk->crc_hcal[ipkt] )return 0; // :-(
4559            for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-(
4560        }
4561    //  -----------------------
4562    //  Check missing packets:
4563    //  -----------------------
4564    //    Readout order:
4565    //    ------------------
4566    //    DSP   packet board
4567    //    ------------------
4568    //    12    0      1
4569    //    10    1      1
4570    //     8    2      1
4571    //     4    3      1
4572    //     6    4      1
4573    //     2    5      1
4574    //    ------------------
4575    //    11    0      2
4576    //     9    1      2
4577    //     7    2      2
4578    //     3    3      2
4579    //     5    4      2
4580    //     1    5      2
4581    //    ------------------
4582    //  -------------------------------------------------
4583    //  Check if it is first or second calibration packet
4584    //  -------------------------------------------------
4585        UInt_t build=0;
4586        TString classname = caltrk->GetName();
4587        UInt_t base=0;
4588        UInt_t mask=0;
4589        if(classname.Contains("CalibTrk1Event")){
4590            base=12;
4591            mask=0x03F000;
4592        }
4593        if(classname.Contains("CalibTrk2Event")){
4594            base=18;
4595            mask=0xFC0000;
4596        }
4597    //  -------------------------------------------------
4598    //  Count number of packets and set build variable
4599    //  -------------------------------------------------
4600        Int_t  npkts=0;
4601        for(Int_t ipkt=0; ipkt<6; ipkt++){
4602            if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){
4603                npkts++;
4604                build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) );
4605            }
4606        }
4607    //    if( npkts==6 )return 1; // :-)
4608    
4609    //    cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl;
4610    
4611    //  -----------------------------------------------
4612    //  If missing packets: check the acq configuration
4613    //  (some DSPs might be excluded from acquisition)
4614    //  -----------------------------------------------
4615    
4616    //  -----------------------------------------------
4617    //  retrieve the first run header after calib
4618    //  -----------------------------------------------
4619         PacketType *pctp;
4620         EventCounter *cod;
4621         cod = eh->GetCounter();
4622         Int_t irun = cod->Get(pctp->RunHeader);
4623         TTree *rh=(TTree*)file->Get("RunHeader");
4624         if ( !rh || rh->IsZombie() ) throw -17;
4625         if( rh->GetEntries() == irun ){
4626              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1)  -- cannot validate :-( "<<endl;
4627             return 0; // :-(
4628         }
4629    
4630         RunHeaderEvent *run  = 0;
4631         EventHeader    *hrun = 0;
4632         rh->SetBranchAddress("RunHeader", &run);
4633         rh->SetBranchAddress("Header", &hrun);
4634         rh->GetEntry(irun);
4635    //     cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl;
4636    
4637         if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){
4638              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl;
4639             return 0; // :-(
4640         }
4641        
4642         if( !run->RM_ACQ_AFTER_CALIB ){
4643              if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0    -- cannot validate :-( "<<endl;
4644             return 0; // :-(
4645         }
4646    
4647         UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime());
4648         if( dtime > timeaftercalib ){
4649             if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl;
4650             return 0; // :-(
4651         }
4652        
4653    
4654    
4655         if( (run->ACQ_BUILD_INFO & mask) != build ){
4656             if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl;
4657             return 0; // :-(
4658         }
4659         return 1; // :-)
4660    
4661    }
4662    
4663    /**
4664     *
4665     * Check the DB (only for overlapping runs at the moment)
4666     *
4667     **/
4668    UInt_t PamelaDBOperations::Check(){
4669      //
4670      UInt_t test = 0;
4671      //
4672      UInt_t thisrht = 0;
4673      UInt_t thisrtt = 0;
4674      UInt_t thisid = 0;
4675      UInt_t prevrht = 0;
4676      UInt_t prevrtt = 0;
4677      UInt_t previd = 0;
4678      //
4679      UInt_t prevl0id = 0;
4680      UInt_t thisl0id = 0;
4681      //
4682      stringstream oss;
4683      TSQLResult *result = 0;
4684      TSQLRow    *row    = 0;
4685      TSQLResult *result2 = 0;
4686      TSQLRow    *row2    = 0;
4687      TSQLResult *result3 = 0;
4688      TSQLRow    *row3    = 0;
4689      oss.str("");
4690      oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME,NEVENTS FROM GL_RUN order by RUNHEADER_TIME asc;";
4691      //  oss << "SELECT ID,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN where ID>10170 and ID<10190 order by RUNHEADER_TIME asc;";
4692      result = conn->Query(oss.str().c_str());
4693      //
4694      if ( !result ) throw -4;;
4695      //
4696      row = result->Next();
4697      //
4698      while ( row ){
4699        thisid = (UInt_t)atoll(row->GetField(0));
4700        thisl0id = (UInt_t)atoll(row->GetField(1));
4701        thisrht = (UInt_t)atoll(row->GetField(2));
4702        thisrtt = (UInt_t)atoll(row->GetField(3));
4703        //
4704        if ( (UInt_t)atoll(row->GetField(4)) > 1 ){
4705          //
4706          //
4707          //
4708          oss.str("");
4709          oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
4710              << thisid << " AND RUNHEADER_TIME="
4711              << thisrht << " AND RUNTRAILER_TIME!="
4712              << thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;";
4713          result3 = conn->Query(oss.str().c_str());
4714          if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str());
4715          if ( result3 ){
4716            row3 = result3->Next();
4717            //
4718            while ( row3 ){
4719              //
4720              // 2 runs with same runheader
4721              //
4722              printf(" CHECK n.4 RUNs %u and %u HAVE SAME RUNHEADER \n",thisid,(UInt_t)atoll(row3->GetField(0)));    
4723              row3 = result3->Next();
4724            };
4725          };
4726          //
4727          oss.str("");
4728          oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
4729              << thisid << " AND RUNHEADER_TIME!="
4730              << thisrht << " AND RUNTRAILER_TIME="
4731              << thisrtt << " AND NEVENTS!=0 AND NEVENTS!=1 order by RUNHEADER_TIME asc;";
4732          result3 = conn->Query(oss.str().c_str());
4733          if ( IsDebug() ) printf(" query is %s \n",oss.str().c_str());
4734          if ( result3 ){
4735            row3 = result3->Next();
4736            //
4737            while ( row3 ){
4738              //
4739              // 2 runs with same runtrailer
4740              //
4741              printf(" CHECK n.5 RUNs %u and %u HAVE SAME RUNTRAILER \n",thisid,(UInt_t)atoll(row3->GetField(0)));
4742              row3 = result3->Next();
4743            };
4744          };
4745          //
4746          oss.str("");
4747          oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
4748              << thisid << " AND RUNHEADER_TIME>"
4749              << thisrht << " AND RUNTRAILER_TIME<"
4750              << thisrtt << " order by RUNHEADER_TIME asc;";
4751          result3 = conn->Query(oss.str().c_str());
4752          if ( result3 ){
4753            row3 = result3->Next();
4754            //
4755            while ( row3 ){
4756              //
4757              // run contained in the checked one
4758              //
4759              printf(" CHECK n.6 RUN %u CONTAINS RUN %u  \n",thisid,(UInt_t)atoll(row3->GetField(0)));
4760              row3 = result3->Next();
4761            };
4762          };
4763          //
4764          oss.str("");
4765          oss << "SELECT ID,ID_ROOT_L0,RUNHEADER_TIME,RUNTRAILER_TIME FROM GL_RUN WHERE ID!="
4766              << thisid << " AND RUNHEADER_TIME="
4767              << thisrht << " AND RUNTRAILER_TIME="
4768              << thisrtt << " AND ID_RUN_FRAG!="
4769              << thisid << " order by RUNHEADER_TIME asc;";
4770          result3 = conn->Query(oss.str().c_str());
4771          if ( result3 ){
4772            row3 = result3->Next();
4773            //
4774            while ( row3 ){
4775              //
4776              // duplicated run
4777              //
4778              printf(" CHECK n.7 RUNs %u and %u HAVE SAME RUNTRAILER AND RUNHEADER (ARE THE SAME?) \n",thisid,(UInt_t)atoll(row3->GetField(0)));
4779              row3 = result3->Next();
4780            };
4781          };
4782        };
4783        //
4784        //    if ( thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt && !(!prevrht && !prevrtt &&!previd) ){
4785        //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){
4786        if ( (thisrht < prevrtt) && (thisrht != prevrht) ){
4787          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4788          printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid);
4789          TString prevf = "";
4790          TString thisf = "";
4791          oss.str("");
4792          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4793          result2 = conn->Query(oss.str().c_str());
4794          if ( !result2 ) throw -4;;
4795          row2 = result2->Next();
4796          prevf = (TString)row2->GetField(0);
4797          oss.str("");
4798          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4799          result2 = conn->Query(oss.str().c_str());
4800          if ( !result2 ) throw -4;;
4801          row2 = result2->Next();
4802          thisf = (TString)row2->GetField(0);
4803          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4804          test = 1;
4805        };
4806        //
4807        if ( (thisrtt < prevrht) && (thisrht != prevrht) ){
4808          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4809          printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid);
4810          TString prevf = "";
4811          TString thisf = "";
4812          oss.str("");
4813          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4814          result2 = conn->Query(oss.str().c_str());
4815          if ( !result2 ) throw -4;
4816          row2 = result2->Next();
4817          prevf = (TString)row2->GetField(0);
4818          oss.str("");
4819          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4820          result2 = conn->Query(oss.str().c_str());
4821          if ( !result2 ) throw -4;;
4822          row2 = result2->Next();
4823          thisf = (TString)row2->GetField(0);
4824          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4825          test = 1;
4826        };
4827        //
4828        if ( (thisrht > thisrtt) && (thisrht != prevrht) ){
4829          if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4830          printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid);
4831          TString prevf = "";
4832          TString thisf = "";
4833          oss.str("");
4834          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)prevl0id <<";";
4835          result2 = conn->Query(oss.str().c_str());
4836          if ( !result2 ) throw -4;;
4837          row2 = result2->Next();
4838          prevf = (TString)row2->GetField(0);
4839          oss.str("");
4840          oss << "SELECT NAME FROM GL_ROOT where ID=" << (UInt_t)thisl0id <<";";
4841          result2 = conn->Query(oss.str().c_str());
4842          if ( !result2 ) throw -4;;
4843          row2 = result2->Next();
4844          thisf = (TString)row2->GetField(0);
4845          if ( IsDebug() ) printf(" ==> files %s and %s \n",prevf.Data(),thisf.Data());
4846          test = 1;
4847        };
4848    
4849        //
4850        prevrht = thisrht;
4851        prevrtt = thisrtt;
4852        previd = thisid;
4853        prevl0id = thisl0id;
4854        row = result->Next();
4855      };
4856      //
4857      return(test);
4858      //
4859    };
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.46
  ViewVC Help
Powered by ViewVC 1.1.23