/[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.5 by pam-fi, Fri Sep 8 14:42:48 2006 UTC revision 1.13 by mocchiut, Fri Oct 20 11:11:12 2006 UTC
# 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 <cTle.h>
35    #include <cEci.h>
36    #include <cJulian.h>
37    
38  #include <PamelaDBOperations.h>  #include <PamelaDBOperations.h>
39  //  //
40  using namespace std;  using namespace std;
41  using namespace pamela;  using namespace pamela;
42    
43    // Some function to work with cTle stuff.
44    bool compTLE(cTle* tle1, cTle *tle2);
45    float getTleJulian(cTle *);
46    string getTleDatetime(cTle*);
47    
48  /**  /**
49   * Constructor.   * Constructor.
50   * @param host         hostname for the SQL connection.   * @param host         hostname for the SQL connection.
# Line 47  using namespace pamela; Line 56  using namespace pamela;
56   * @param obt0         file obt0.   * @param obt0         file obt0.
57   * @param tsync        file timesync.   * @param tsync        file timesync.
58   * @param debug        debug flag.   * @param debug        debug flag.
59     * @param tlefilename  ascii file with TLE 3 line elements.
60   */   */
61  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){
62    //    //
63    //    //
64    SetConnection(host,user,password);    SetConnection(host,user,password);
# 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    //    //
   SetRootName(filerootname);  
   SetRawName(filerawname);  
77    //    //
78    this->OpenFile();    INSERT_RAW =!filerawname.IsNull();
79      if(INSERT_RAW)SetRawName(filerawname);
80      //
81      INSERT_ROOT = !filerootname.IsNull();
82      if( INSERT_ROOT ){
83        this->SetRootName(filerootname);
84        this->SetOrbitNo();
85        file = TFile::Open(this->GetRootName().Data());
86      };
87    //    //
88    this->SetID_RAW(0);    this->SetID_RAW(0);
89    this->SetID_ROOT(0);    this->SetID_ROOT(0);
90    
91      VALIDATE = false;
92      
93    //    //
94  };  };
95    
# Line 78  PamelaDBOperations::PamelaDBOperations(T Line 98  PamelaDBOperations::PamelaDBOperations(T
98   */   */
99  void PamelaDBOperations::Close(){  void PamelaDBOperations::Close(){
100    if( conn && conn->IsConnected() ) conn->Close();    if( conn && conn->IsConnected() ) conn->Close();
101      delete clean_time;
102    delete glrun;    delete glrun;
103    delete this;    delete this;
104  };  };
# Line 86  void PamelaDBOperations::Close(){ Line 107  void PamelaDBOperations::Close(){
107  // SETTERS  // SETTERS
108  //  //
109    
110    //
111    // must be out of the constructor in order to FORCE the validation of the latest runs in case you run the validation together with the latest file
112    //
113    void PamelaDBOperations::CheckValidate(Long64_t olderthan){
114      clean_time = new TDatime();
115      if(olderthan >= 0){
116        VALIDATE = true;
117        UInt_t timelim = 0;
118        timelim =  (UInt_t)clean_time->Convert() - olderthan;
119        clean_time->Set(timelim,false);
120      };
121    };
122    
123  /**  /**
124   * Open the DB connection   * Open the DB connection
125   * @param host         hostname for the SQL connection.   * @param host         hostname for the SQL connection.
# Line 93  void PamelaDBOperations::Close(){ Line 127  void PamelaDBOperations::Close(){
127   * @param password     password for the SQL connection.   * @param password     password for the SQL connection.
128   */   */
129  void PamelaDBOperations::SetConnection(TString host, TString user, TString password){  void PamelaDBOperations::SetConnection(TString host, TString user, TString password){
130      if ( IsDebug() ) printf(" Connecting using host = %s user = %s password = %s \n",host.Data(),user.Data(),password.Data());
131    conn = TSQLServer::Connect(host.Data(),user.Data(),password.Data());    conn = TSQLServer::Connect(host.Data(),user.Data(),password.Data());
132  };  };
133    
# Line 161  void PamelaDBOperations::SetRootName(TSt Line 196  void PamelaDBOperations::SetRootName(TSt
196  };  };
197    
198  /**  /**
199     * Store the downlink orbit number from filename.
200     */
201    void PamelaDBOperations::SetOrbitNo(){
202      dworbit = 0;
203      TString name = this->GetRootFile();
204      Int_t nlength = name.Length();
205      if ( nlength < 5 ) return;
206      TString dwo = 0;
207      for (Int_t i = 0; i<5; i++){
208        dwo.Append(name[i],1);
209      };
210      if ( dwo.IsDigit() ){
211        dworbit = (UInt_t)dwo.Atoi();
212      } else {
213        dwo="";
214        for (Int_t i = 8; i<13; i++){
215          dwo.Append(name[i],1);
216        };    
217        if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
218      };
219      if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());
220      return;
221    };
222    
223    
224    
225    /**
226   * Store the NOBOOT flag.   * Store the NOBOOT flag.
227   * @param noboot    true/false.   * @param noboot    true/false.
228   */   */
# Line 169  void PamelaDBOperations::SetNOBOOT(Bool_ Line 231  void PamelaDBOperations::SetNOBOOT(Bool_
231  };  };
232    
233  /**  /**
234     * Store path to the TLE file.
235     */
236    void PamelaDBOperations::SetTLEPath(TString str){
237      tlefilename = str;
238    };
239    
240    /**
241   * Store the olderthan variable   * Store the olderthan variable
242   * @param olderthan   * @param olderthan
243   */   */
244  void PamelaDBOperations::SetOlderThan(Long64_t oldthan){  // void PamelaDBOperations::SetOlderThan(Long64_t oldthan){
245    olderthan = oldthan;  //   olderthan = oldthan;
246  };  // };
247    
248  /**  /**
249   * Retrieve the ID_RAW, if exists, returns NULL if does not exist.   * Retrieve the ID_RAW, if exists, returns NULL if does not exist.
# Line 225  Int_t PamelaDBOperations::SetUpperLimits Line 294  Int_t PamelaDBOperations::SetUpperLimits
294    UInt_t jump = 50000; // was 5000    UInt_t jump = 50000; // was 5000
295    EventCounter *code=0;    EventCounter *code=0;
296    //    //
297      UInt_t deltapkt = 5000;
298      ULong64_t deltaobt = 50000;
299      //
300    //   pcksList packetsNames;    //   pcksList packetsNames;
301    //   pcksList::iterator Iter;    //   pcksList::iterator Iter;
302    //   getPacketsNames(packetsNames);    //   getPacketsNames(packetsNames);
# Line 271  Int_t PamelaDBOperations::SetUpperLimits Line 343  Int_t PamelaDBOperations::SetUpperLimits
343    //    //
344    if ( nevent < 2 ) return(4);    if ( nevent < 2 ) return(4);
345    //    //
346    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) ){
347        //
348      // go back      // go back
349      zomp = nevent - 2;      zomp = nevent - 2;
350      //      //
# Line 320  Int_t PamelaDBOperations::SetUpperLimits Line 393  Int_t PamelaDBOperations::SetUpperLimits
393    // check if last runtrailer is within limits, if not extend limits (one should check for all packets but we need only runtrailer)    // check if last runtrailer is within limits, if not extend limits (one should check for all packets but we need only runtrailer)
394    //    //
395    PacketType *pctp=0;    PacketType *pctp=0;
   T->GetEntry(upperentry);  
   code = eh->GetCounter();  
   Int_t lasttrail = code->Get(pctp->RunTrailer);  
   Int_t lasthead = code->Get(pctp->RunHeader);  
396    TTree *rh=(TTree*)file->Get("RunHeader");    TTree *rh=(TTree*)file->Get("RunHeader");
397    if ( !rh || rh->IsZombie() ) throw -17;    if ( !rh || rh->IsZombie() ) throw -17;
398    TTree *rt=(TTree*)file->Get("RunTrailer");    TTree *rt=(TTree*)file->Get("RunTrailer");
# Line 346  Int_t PamelaDBOperations::SetUpperLimits Line 415  Int_t PamelaDBOperations::SetUpperLimits
415    UInt_t pkth = 0;    UInt_t pkth = 0;
416    ULong64_t obth = 0;    ULong64_t obth = 0;
417    //    //
418      T->GetEntry(upperentry);
419      code = eh->GetCounter();
420      Int_t lasttrail = code->Get(pctp->RunTrailer);
421      Int_t lasthead = code->Get(pctp->RunHeader);
422    if ( lasttrail < rtev ){    if ( lasttrail < rtev ){
423      rt->GetEntry(lasttrail);      rt->GetEntry(lasttrail);
424      pht = eht->GetPscuHeader();      pht = eht->GetPscuHeader();
# Line 356  Int_t PamelaDBOperations::SetUpperLimits Line 429  Int_t PamelaDBOperations::SetUpperLimits
429    if ( lasthead < rhev ){    if ( lasthead < rhev ){
430      rh->GetEntry(lasthead);      rh->GetEntry(lasthead);
431      phh = ehh->GetPscuHeader();      phh = ehh->GetPscuHeader();
432      pkth = PKT(pht->GetCounter());      pkth = PKT(phh->GetCounter());
433      obth = OBT(pht->GetOrbitalTime());      obth = OBT(phh->GetOrbitalTime());
434    };    };
435    //    //
436    if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);    if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
# Line 401  Int_t PamelaDBOperations::SetUpperLimits Line 474  Int_t PamelaDBOperations::SetUpperLimits
474      pkth = PKT(phh->GetCounter());      pkth = PKT(phh->GetCounter());
475      obth = OBT(phh->GetOrbitalTime());      obth = OBT(phh->GetOrbitalTime());
476      //      //
477    //    if ( IsDebug() ) printf(" k %i rhev before %i ph %u upperp %u oh %u uppero %u \n",k,rhev,pkth,spkth,obth,sobth);
478        //
479      if ( pkth < spkth && obth < sobth ){      if ( pkth < spkth && obth < sobth ){
       if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);  
480        if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);        if ( IsDebug() ) printf(" RH PROBLEMS determining the event repetition at the end of the file lasthead %i  \n",rhev);
481        //        //
482        rhev--;        rhev = k-1;
483        rh->GetEntry(rhev);        rh->GetEntry(rhev);
484        pkth = spkth;        pkth = spkth;
485        obth = sobth;        obth = sobth;
       if ( IsDebug() ) printf(" lasthead %i pt %i p1 %i ot %u o1 %u \n",rhev,pkth,spkth,obth,sobth);  
486        //        //
487        UInt_t evbefh = 0;        UInt_t evbefh = 0;
488        code = ehh->GetCounter();        code = ehh->GetCounter();
# Line 419  Int_t PamelaDBOperations::SetUpperLimits Line 492  Int_t PamelaDBOperations::SetUpperLimits
492          ph = eh->GetPscuHeader();          ph = eh->GetPscuHeader();
493          t_pktlast = PKT(ph->GetCounter());          t_pktlast = PKT(ph->GetCounter());
494          t_obtlast = OBT(ph->GetOrbitalTime());          t_obtlast = OBT(ph->GetOrbitalTime());
495          if ( t_pktlast < spkth && t_obtlast < sobth ){ // jump          if ( t_pktlast <= spkth && t_obtlast <= sobth ){ // jump
496            upperpkt = pkth;            upperpkt = pkth;
497            upperobt = obth;            upperobt = obth;
498            upperentry = evbefh-1;            upperentry = evbefh-1;
# Line 439  Int_t PamelaDBOperations::SetUpperLimits Line 512  Int_t PamelaDBOperations::SetUpperLimits
512          };              };    
513        };        };
514        if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);        if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
515          goto kikko0;
516      };      };
517    };    };
518     kikko0:
519    //    //
520    //    //
521    //    //
# Line 458  Int_t PamelaDBOperations::SetUpperLimits Line 533  Int_t PamelaDBOperations::SetUpperLimits
533      pktt = PKT(pht->GetCounter());      pktt = PKT(pht->GetCounter());
534      obtt = OBT(pht->GetOrbitalTime());      obtt = OBT(pht->GetOrbitalTime());
535      //      //
536    //    if ( IsDebug() ) printf(" k %i rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",k,rtev,pktt,spktt,obtt,sobtt);
537        //
538      if ( pktt < spktt && obtt < sobtt ){      if ( pktt < spktt && obtt < sobtt ){
       if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);  
539        if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);        if ( IsDebug() ) printf(" RT PROBLEMS determining the event repetition at the end of the file lasttrail %i \n",rtev);
540        //        //
541        rtev--;        rtev = k-1;
542        rt->GetEntry(rtev);        rt->GetEntry(rtev);
543        pktt = spktt;        pktt = spktt;
544        obtt = sobtt;        obtt = sobtt;
# Line 476  Int_t PamelaDBOperations::SetUpperLimits Line 552  Int_t PamelaDBOperations::SetUpperLimits
552          ph = eh->GetPscuHeader();          ph = eh->GetPscuHeader();
553          t_pktlast = PKT(ph->GetCounter());          t_pktlast = PKT(ph->GetCounter());
554          t_obtlast = OBT(ph->GetOrbitalTime());          t_obtlast = OBT(ph->GetOrbitalTime());
555          if ( t_pktlast < spktt && t_obtlast < sobtt ){ // jump          if ( t_pktlast <= spktt && t_obtlast <= sobtt ){ // jump
556            upperpkt = pktt;            upperpkt = pktt;
557            upperobt = obtt;            upperobt = obtt;
558            upperentry = evbeft-1;            upperentry = evbeft-1;
# Line 496  Int_t PamelaDBOperations::SetUpperLimits Line 572  Int_t PamelaDBOperations::SetUpperLimits
572          };          };
573        };        };
574        if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);        if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
575        break;        goto kikko;
576          //      break;
577        //        //
578      };      };
579      //        //  
580    };    };
581    //    //
582    // kikko:   kikko:
583      //
584      T->GetEntry(upperentry);
585      code = eh->GetCounter();
586      lasttrail = code->Get(pctp->RunTrailer);
587      lasthead = code->Get(pctp->RunHeader);
588      if ( lasttrail < rtev ){
589        rt->GetEntry(lasttrail);
590        pht = eht->GetPscuHeader();
591        pktt = PKT(pht->GetCounter());
592        obtt = OBT(pht->GetOrbitalTime());
593      };
594      //
595      if ( lasthead < rhev ){
596        rh->GetEntry(lasthead);
597        phh = ehh->GetPscuHeader();
598        pkth = PKT(phh->GetCounter());
599        obth = OBT(phh->GetOrbitalTime());
600      };
601      //
602      if ( IsDebug() ) printf(" rhev before %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
603      if ( pkth > upperpkt && obth > upperobt ){
604        if ( IsDebug() ) printf(" Upper limits extended to include last header: ph %i upperp %i oh %llu uppero %llu \n",pkth,upperpkt,obth,upperobt);
605        upperpkt = pkth;
606        upperobt = obth;
607        rhev = lasthead+1;
608      } else {
609        rhev = lasthead;
610      };
611      if ( IsDebug() ) printf(" rhev after %i ph %i upperp %i oh %llu uppero %llu \n",rhev,pkth,upperpkt,obth,upperobt);
612      //
613      if ( IsDebug() ) printf(" rtev beforev %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
614      if ( pktt > upperpkt && obtt > upperobt ){
615        if ( IsDebug() ) printf(" Upper limits extended to include last trailer: pt %i upperp %i ot %llu uppero %llu \n",pktt,upperpkt,obtt,upperobt);
616        upperpkt = pktt;
617        upperobt = obtt;
618        rtev = lasttrail+1;
619      } else {
620        rtev = lasttrail;
621      };
622      if ( IsDebug() ) printf(" rtev after %i  pt %i upperp %i ot %llu uppero %llu \n",rtev,pktt,upperpkt,obtt,upperobt);
623    //    //
624    if ( IsDebug() ) printf(" Upper limits are: OBT %llu pkt_num %i upper entry %i \n",upperobt,upperpkt,upperentry);    if ( IsDebug() ) printf(" Upper limits are: OBT %llu pkt_num %i upper entry %i \n",upperobt,upperpkt,upperentry);
625    //    //
626    return(0);    return(0);
627  }  }
628    
629    /**
630     *
631     * Trick to have unique RUN ID numbers even when runs are deleted and mysql deamon restarted.
632     * Entries in the _RUNID_GEN table are never deleted.
633     *
634     **/
635    UInt_t PamelaDBOperations::AssignRunID(){
636      //
637      TSQLResult *result = 0;
638      TSQLRow *row = 0;
639      UInt_t runid = 0;
640      //
641      stringstream   oss;
642      //  
643      oss.str("");
644      oss << "INSERT INTO _RUNID_GEN VALUES (NULL);";
645      result = conn->Query(oss.str().c_str());
646      if ( !result ) throw -10;
647      oss.str("");
648      oss << "SELECT ID FROM _RUNID_GEN ORDER BY ID DESC LIMIT 1;";
649      result = conn->Query(oss.str().c_str());
650      if ( !result ) throw -10;
651      //
652      row = result->Next();
653      //
654      if ( !row ) throw -28;
655      //
656      runid = (UInt_t)atoll(row->GetField(0));
657      //
658      return(runid);
659    };
660    
661  //  //
662  // GETTERS  // GETTERS
663  //  //
# Line 584  const PacketType* PamelaDBOperations::Ge Line 733  const PacketType* PamelaDBOperations::Ge
733  // PRIVATE FUNCTIONS  // PRIVATE FUNCTIONS
734  //  //
735    
736  /**  // /**
737   * Open the ROOT filename for reading  // * Open the ROOT filename for reading
738   */  // */
739  void PamelaDBOperations::OpenFile(){  // void PamelaDBOperations::OpenFile(){
740    file = TFile::Open(this->GetRootName().Data());  //   file = TFile::Open(this->GetRootName().Data());
741    //   //
742    
743    void PamelaDBOperations::CheckFile(){  
744      if ( !file ) throw -12;
745  };  };
746    
747    
748  /**  /**
749   * Check if LEVEL0 file and DB connection have really be opened   * Check if LEVEL0 file and DB connection have really be opened
750   */   */
751  void PamelaDBOperations::CheckFile(){    void PamelaDBOperations::CheckConnection(){  
   //  
   if ( !file ) throw -12;  
752    //    //
753    // check connection    // check connection
754    //    //
755    if( !conn ) throw -1;        if( !conn ) throw -1;
756    bool connect = conn->IsConnected();    bool connect = conn->IsConnected();
757    if( !connect ) throw -1;        if( !connect ) throw -1;
758      if ( !dworbit ) throw -27;
759  };  };
760    
761  /**  /**
# Line 783  Int_t PamelaDBOperations::insertPamelaGL Line 935  Int_t PamelaDBOperations::insertPamelaGL
935    if ( this->GetID_RAW() == 0 )  throw -11;    if ( this->GetID_RAW() == 0 )  throw -11;
936    //    //
937    oss.str("");    oss.str("");
938    oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ID_RAW<= "    oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='"
939        << id << " AND TO_ID_RAW >= "        << this->GetRawFile().Data() << "';";
       << id << ";";  
940    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
941    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
942    if ( !result ) throw -10;    if ( !result ) throw -10;
943    row = result->Next();    row = result->Next();
944    if ( !row ) throw -10;    //
945      if ( !row ){
946        oss.str("");
947        oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "
948            << dworbit << " order by FROM_ORBIT desc limit 1;";
949        if ( IsDebug() ) printf(" %s \n",oss.str().c_str());
950        result = conn->Query(oss.str().c_str());
951        if ( !result ) throw -10;
952        row = result->Next();
953        if ( !row ) throw -10;
954      };
955    //    //
956    t0 = (UInt_t)TDatime(row->GetField(0)).Convert();    t0 = (UInt_t)TDatime(row->GetField(0)).Convert();
957    /*    /*
# Line 864  Int_t PamelaDBOperations::insertPamelaGL Line 1025  Int_t PamelaDBOperations::insertPamelaGL
1025          //          //
1026          TYPE = 55;//224;          TYPE = 55;//224;
1027          //          //
1028            if ( IsDebug() ) printf("mcmd tsync %i tsync %u obt %u \n",i,TSYNC,OBT);
1029            //
1030          if ( TSYNC && OBT ){          if ( TSYNC && OBT ){
1031            existsts = true;            existsts = true;
1032            goto eout;            goto eout;
# Line 886  Int_t PamelaDBOperations::insertPamelaGL Line 1049  Int_t PamelaDBOperations::insertPamelaGL
1049      //      //
1050      rt->SetBranchAddress("RunTrailer", &runt);      rt->SetBranchAddress("RunTrailer", &runt);
1051      //      //
1052      if ( rhev > 0 ){      Int_t nrhev = rh->GetEntries();
1053        rh->GetEntry(0);      Int_t nrtev = rt->GetEntries();
1054        //      if ( IsDebug() ) printf(" ou nevent %i rhev %i rtev %i \n",nevent,nrhev,nrtev);
1055        TSYNC = runh->LAST_TIME_SYNC_INFO;      //
1056        OBT = runh->OBT_TIME_SYNC * 1000;      if ( nrhev > 0 ){
1057        //        for (Int_t i=0; i<nrhev; i++){
1058        TYPE = 20;          //
1059        //          rh->GetEntry(i);
1060        if ( TSYNC && OBT ){          //
1061          existsts = true;          TSYNC = runh->LAST_TIME_SYNC_INFO;
1062          goto eout;          OBT = runh->OBT_TIME_SYNC * 1000;
1063            //
1064            TYPE = 20;
1065            //
1066            if ( IsDebug() ) printf("runheader %i tsync %u obt %u \n",i,TSYNC,OBT);
1067            //
1068            if ( TSYNC && OBT ){
1069              existsts = true;
1070              goto eout;
1071            };
1072        };        };
1073        //        //
1074      };      };
1075      if ( rtev > 0 ){      if ( nrtev > 0 ){
1076        //        //
1077        if ( IsDebug() ) printf(" No runheader \n");        if ( IsDebug() ) printf(" No runheader \n");
1078        signal = 6;        signal = 6;
1079        //        //
1080        rt->GetEntry(0);        for (Int_t i=0; i<nrtev; i++){
1081        //          //
1082        TSYNC = runt->LAST_TYME_SYNC_INFO;          rt->GetEntry(i);
1083        OBT = runt->OBT_TYME_SYNC * 1000;          //
1084        //          TSYNC = runt->LAST_TYME_SYNC_INFO;
1085        TYPE = 21;          OBT = runt->OBT_TYME_SYNC * 1000;
1086        //          //
1087        if ( TSYNC && OBT ){          TYPE = 21;
1088          existsts = true;          //
1089          goto eout;          if ( IsDebug() ) printf("runtrailer %i tsync %u obt %u \n",i,TSYNC,OBT);
1090            //
1091            if ( TSYNC && OBT ){
1092              existsts = true;
1093              goto eout;
1094            };
1095        };        };
1096        //        //
1097      } else {      } else {
# Line 1454  void PamelaDBOperations::HandleRun(){   Line 1631  void PamelaDBOperations::HandleRun(){  
1631      //      //
1632      this->FillClass();      this->FillClass();
1633      //      //
1634      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);      if ( !IsRunAlreadyInserted() ){
1635          glrun->SetID(this->AssignRunID());
1636          glrun->SetID_RUN_FRAG(0);
1637          glrun->Fill_GL_RUN(conn);
1638        };
1639    } else {    } else {
1640      //      //
1641      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 1626  void PamelaDBOperations::HandleRunFragme Line 1807  void PamelaDBOperations::HandleRunFragme
1807        };        };
1808        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
1809        //        //
       glrun1->SetID(0);  
1810        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
1811        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
1812        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());        glrun1->SetRUNTRAILER_TIME(glrun->GetRUNTRAILER_TIME());
# Line 1654  void PamelaDBOperations::HandleRunFragme Line 1834  void PamelaDBOperations::HandleRunFragme
1834        //        //
1835        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
1836          //          //
1837            glrun->SetID(this->AssignRunID());
1838            glrun->SetID_RUN_FRAG(glrun1->GetID());
1839          glrun->Fill_GL_RUN(conn);          glrun->Fill_GL_RUN(conn);
1840          //          //
1841          // get id number          // get id number
1842          //          //
1843          oss.str("");  //      oss.str("");
1844          oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "  //      oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "
1845              << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "  //          << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "
1846              << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";  //          << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";
1847          //  //      //
1848          if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());  //      if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());
1849          //  //      //
1850          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
1851          if ( !result ) throw -4;  //      if ( !result ) throw -4;
1852          row = result->Next();  //      row = result->Next();
1853          //  //      //
1854          UInt_t idrun0 = 0;  //      UInt_t idrun0 = 0;
1855          if ( !row ){  //      if ( !row ){
1856            throw -10;  //        throw -10;
1857          } else {  //      } else {
1858            idrun0 = (UInt_t)atoll(row->GetField(0));  //        idrun0 = (UInt_t)atoll(row->GetField(0));
1859          };  //      };
         //  
         glrun1->SetID_RUN_FRAG(idrun0);  
1860          //          //
1861            //      glrun1->SetID_RUN_FRAG(idrun0);
1862            glrun1->SetID_RUN_FRAG(glrun->GetID());
1863          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
1864          //          //
1865          oss.str("");  //      oss.str("");
1866          oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";  //      oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";
1867          //  //      //
1868          if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());  //      if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());
1869          //  //      //
1870          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
1871          if ( !result ) throw -4;  //      if ( !result ) throw -4;
1872          row = result->Next();  //      row = result->Next();
1873          //  //      //
1874          UInt_t idrun1 = 0;  //      UInt_t idrun1 = 0;
1875          if ( !row ){  //      if ( !row ){
1876            throw -10;  //        throw -10;
1877          } else {  //      } else {
1878            idrun1 = (UInt_t)atoll(row->GetField(0));  //        idrun1 = (UInt_t)atoll(row->GetField(0));
1879          };  //      };
1880          //  //      //
1881          oss.str("");  //      oss.str("");
1882          oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";  //      oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";
1883          //  //      //
1884          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
1885          if ( !result ) throw -4;  //      if ( !result ) throw -4;
1886          //          //
1887        };        };
1888        //        //
# Line 1850  void PamelaDBOperations::HandleRunFragme Line 2032  void PamelaDBOperations::HandleRunFragme
2032        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());        glrun->SetRUNTRAILER_OBT(glrun1->GetRUNTRAILER_OBT());
2033        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());        glrun->SetRUNTRAILER_PKT(glrun1->GetRUNTRAILER_PKT());
2034        //        //
       glrun1->SetID(0);  
2035        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());        glrun1->SetRUNHEADER_TIME(glrun->GetRUNHEADER_TIME());
2036        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());        glrun1->SetRUNHEADER_OBT(glrun->GetRUNHEADER_OBT());
2037        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());        glrun1->SetRUNHEADER_PKT(glrun->GetRUNHEADER_PKT());
# Line 1870  void PamelaDBOperations::HandleRunFragme Line 2051  void PamelaDBOperations::HandleRunFragme
2051        //        //
2052        if ( !IsRunAlreadyInserted() ){        if ( !IsRunAlreadyInserted() ){
2053          //          //
2054            glrun->SetID(this->AssignRunID());
2055            //
2056            glrun->SetID_RUN_FRAG(glrun1->GetID());
2057          glrun->Fill_GL_RUN(conn);          glrun->Fill_GL_RUN(conn);
2058          //          //
2059          // get id number          // get id number
2060          //          //
2061          oss.str("");  //      oss.str("");
2062          oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "  //      oss << " SELECT ID FROM GL_RUN WHERE BOOT_NUMBER="<<this->GetBOOTnumber()<<" AND "
2063              << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "  //          << " RUNHEADER_OBT="<<glrun->GetRUNHEADER_OBT()<<" AND "
2064              << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";  //          << " RUNTRAILER_OBT="<<glrun->GetRUNTRAILER_OBT()<<";";
2065          //  //      //
2066          if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());  //      if ( IsDebug() ) printf(" search for idrun0 , query is : \n%s\n",oss.str().c_str());
2067          //  //      //
2068          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
2069          if ( !result ) throw -4;  //      if ( !result ) throw -4;
2070          row = result->Next();  //      row = result->Next();
2071          //  //      //
2072          UInt_t idrun0 = 0;  //      UInt_t idrun0 = 0;
2073          if ( !row ){  //      if ( !row ){
2074            throw -10;  //        throw -10;
2075          } else {  //      } else {
2076            idrun0 = (UInt_t)atoll(row->GetField(0));  //        idrun0 = (UInt_t)atoll(row->GetField(0));
2077          };  //      };
2078          //  //      //
2079          glrun1->SetID_RUN_FRAG(idrun0);  //      glrun1->SetID_RUN_FRAG(idrun0);
2080          //  //
2081            glrun1->SetID_RUN_FRAG(glrun->GetID());
2082          glrun1->Fill_GL_RUN(conn);          glrun1->Fill_GL_RUN(conn);
2083          //          //
2084          oss.str("");  //      oss.str("");
2085          oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";  //      oss << " SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG="<<idrun0<<" ;";
2086          //  //      //
2087          if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());  //      if ( IsDebug() ) printf(" search for idrun1 , query is : \n%s\n",oss.str().c_str());
2088          //  //      //
2089          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
2090          if ( !result ) throw -4;  //      if ( !result ) throw -4;
2091          row = result->Next();  //      row = result->Next();
2092          //  //      //
2093          UInt_t idrun1 = 0;  //      UInt_t idrun1 = 0;
2094          if ( !row ){  //      if ( !row ){
2095            throw -10;  //        throw -10;
2096          } else {  //      } else {
2097            idrun1 = (UInt_t)atoll(row->GetField(0));  //        idrun1 = (UInt_t)atoll(row->GetField(0));
2098          };  //      };
2099          //  //      //
2100          oss.str("");  //      oss.str("");
2101          oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";  //      oss << " UPDATE GL_RUN SET ID_RUN_FRAG="<<idrun1<<" WHERE ID="<<idrun0<<" ;";
2102          //  //      //
2103          result =  conn->Query(oss.str().c_str());  //      result =  conn->Query(oss.str().c_str());
2104          if ( !result ) throw -4;  //      if ( !result ) throw -4;
2105          //          //
2106        };        };
2107        //        //
# Line 2000  void PamelaDBOperations::HandleRunFragme Line 2185  void PamelaDBOperations::HandleRunFragme
2185          if ( IsDebug() ) printf(" The run is new \n");          if ( IsDebug() ) printf(" The run is new \n");
2186          if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n");          if ( IsDebug() ) printf(" -> fill the GL_RUNFRAGMENTS table \n");
2187          //          //
2188            glrun->SetID(this->AssignRunID());
2189            glrun->SetID_RUN_FRAG(0);
2190          glrun->Fill_GL_RUN_FRAGMENTS(conn);          glrun->Fill_GL_RUN_FRAGMENTS(conn);
2191          //          //
2192        } else {        } else {
# Line 2032  void PamelaDBOperations::HandleMissingHo Line 2219  void PamelaDBOperations::HandleMissingHo
2219      //      //
2220      this->FillClass(mishead,mistrail,firstev,lastev);      this->FillClass(mishead,mistrail,firstev,lastev);
2221      //      //
2222      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2223          glrun->SetID(this->AssignRunID());
2224          glrun->SetID_RUN_FRAG(0);
2225          glrun->Fill_GL_RUN(conn);    
2226        };
2227      //      //
2228    };    };
2229    //    //
# Line 2199  Bool_t PamelaDBOperations::IsRunConsiste Line 2390  Bool_t PamelaDBOperations::IsRunConsiste
2390            //                  //      
2391            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2392            //            //
2393            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2394                glrun->SetID(this->AssignRunID());
2395                glrun->SetID_RUN_FRAG(0);
2396                glrun->Fill_GL_RUN(conn);      
2397              };
2398            //            //
2399            firstevno = lastentry + 1;            firstevno = lastentry + 1;
2400            //            //
# Line 2275  void PamelaDBOperations::HandleSuspiciou Line 2470  void PamelaDBOperations::HandleSuspiciou
2470      if ( IsDebug() ) printf(" Checking but no events in the run! \n");      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
2471      //      //
2472      this->FillClass();      this->FillClass();
2473      if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);          if ( !IsRunAlreadyInserted() ){
2474          glrun->SetID(this->AssignRunID());
2475          glrun->SetID_RUN_FRAG(0);
2476          glrun->Fill_GL_RUN(conn);    
2477        };
2478      //      //
2479    } else {    } else {
2480      //      //
# Line 2303  void PamelaDBOperations::HandleSuspiciou Line 2502  void PamelaDBOperations::HandleSuspiciou
2502        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");
2503        //        //
2504        this->FillClass();        this->FillClass();
2505        if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                if ( !IsRunAlreadyInserted() ){
2506            glrun->SetID(this->AssignRunID());
2507            glrun->SetID_RUN_FRAG(0);
2508            glrun->Fill_GL_RUN(conn);        
2509          };
2510        //        //
2511      } else {      } else {
2512        //        //
# Line 2401  void PamelaDBOperations::HandleSuspiciou Line 2604  void PamelaDBOperations::HandleSuspiciou
2604            //                  //      
2605            this->SetCommonGLRUN(firstTime,lastTime);            this->SetCommonGLRUN(firstTime,lastTime);
2606            //            //
2607            if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                  if ( !IsRunAlreadyInserted() ){
2608                glrun->SetID(this->AssignRunID());
2609                glrun->SetID_RUN_FRAG(0);
2610                glrun->Fill_GL_RUN(conn);      
2611              };
2612            //            //
2613            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...
2614              //              //
# Line 2436  void PamelaDBOperations::HandleSuspiciou Line 2643  void PamelaDBOperations::HandleSuspiciou
2643              //                  //    
2644              this->SetCommonGLRUN(firstTime,lastTime);              this->SetCommonGLRUN(firstTime,lastTime);
2645              //              //
2646              if ( !IsRunAlreadyInserted() ) glrun->Fill_GL_RUN(conn);                    if ( !IsRunAlreadyInserted() ){
2647                  glrun->SetID(this->AssignRunID());
2648                  glrun->SetID_RUN_FRAG(0);
2649                  glrun->Fill_GL_RUN(conn);      
2650                };
2651            };            };
2652            //            //
2653            firstevno = lastentry + 1;            firstevno = lastentry + 1;
# Line 3124  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3335  Int_t PamelaDBOperations::CleanGL_RUN_FR
3335    TSQLRow    *row2   = 0;    TSQLRow    *row2   = 0;
3336    //    //
3337    UInt_t moved = 0;    UInt_t moved = 0;
3338    UInt_t timelim = 0;  //  UInt_t timelim = 0;
3339    TDatime *time = new TDatime();  //  TDatime *time = new TDatime();
3340    //    //
3341    stringstream oss;    stringstream oss;
3342    oss.str("");    oss.str("");
3343    //    //
3344    //    //
   //  
   if ( olderthan < 0 ){  
     if ( IsDebug() ) printf(" Do not clean GL_RUN_FRAGMENTS table \n");  
     return(1);  
   };  
   //  
   // timelim = now - olderthan  
   //  
   time->Set();  
   timelim =  (UInt_t)time->Convert() - olderthan;  
   time->Set(timelim,false);  
   //  
3345    // check if there are entries older than "olderthan" seconds from now    // check if there are entries older than "olderthan" seconds from now
3346    //    //
3347    oss.str("");    oss.str("");
3348    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"    oss << " SELECT ID FROM GL_RUN_FRAGMENTS WHERE"
3349        << " INSERT_TIME <= '" << time->AsSQLString() << "';";        << " INSERT_TIME <= '" << clean_time->AsSQLString() << "';";
3350    //    //
3351    if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",time->AsSQLString(),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());
3352    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
3353    //    //
3354    if ( result ){    if ( result ){
3355      //      //
3356      row = result->Next();      row = result->Next();
3357      //      //
3358      if ( row ){      while ( row ){
3359        //        //
3360        oss.str("");        oss.str("");
3361        oss << " ID= "<< row->GetField(0);        oss << " ID= "<< row->GetField(0);
# Line 3191  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3390  Int_t PamelaDBOperations::CleanGL_RUN_FR
3390          if ( IsDebug() ) printf(" The run is new \n");          if ( IsDebug() ) printf(" The run is new \n");
3391          if ( IsDebug() ) printf(" -> fill the DB \n");                if ( IsDebug() ) printf(" -> fill the DB \n");      
3392          //          //
3393          glrun->SetID(0);          //      glrun->SetID(this->AssignRunID()); we use the old run number!
3394            glrun->SetID_RUN_FRAG(glrun->GetID());
3395          glrun->Fill_GL_RUN(conn);            glrun->Fill_GL_RUN(conn);  
3396          //          //
3397          oss.str("");  //      oss.str("");
3398          oss << " SELECT ID FROM GL_RUN WHERE "  //      oss << " SELECT ID FROM GL_RUN WHERE "
3399              << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "  //          << " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND "
3400              << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "  //          << " RUNHEADER_PKT=" << (UInt_t)glrun->GetRUNHEADER_PKT() << " AND "
3401              << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "  //          << " RUNTRAILER_PKT=" << (UInt_t)glrun->GetRUNTRAILER_PKT() << " AND "
3402              << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "  //          << " RUNHEADER_OBT=" << (UInt_t)glrun->GetRUNHEADER_OBT() << " AND "
3403              << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";  //          << " RUNTRAILER_OBT=" << (UInt_t)glrun->GetRUNTRAILER_OBT() << "; ";
3404          //  //      //
3405          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());
3406          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3407          //  //      //
3408          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3409          //  //      //
3410          row2 = result2->Next();  //      row2 = result2->Next();
3411          //  //      //
3412          if ( !row2 ) throw -25;  //      if ( !row2 ) throw -25;
3413          //  //      //
3414          oss.str("");  //      oss.str("");
3415          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);
3416          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());
3417          result2 = conn->Query(oss.str().c_str());  //      result2 = conn->Query(oss.str().c_str());
3418          //  //      //
3419          if ( !result2 ) throw -4;  //      if ( !result2 ) throw -4;
3420          //          //
3421          moved++;          moved++;
3422          //          //
# Line 3233  Int_t PamelaDBOperations::CleanGL_RUN_FR Line 3433  Int_t PamelaDBOperations::CleanGL_RUN_FR
3433        //        //
3434        if ( !result2 ) throw -4;        if ( !result2 ) throw -4;
3435        //        //
3436          row = result->Next();
3437      };      };
3438    };    };
3439    if ( IsDebug() ) printf(" Moved %u runs\n",moved);    if ( IsDebug() ) printf(" Moved %u runs\n",moved);
# Line 3247  Int_t PamelaDBOperations::ValidateRuns() Line 3448  Int_t PamelaDBOperations::ValidateRuns()
3448    TSQLResult *result = 0;    TSQLResult *result = 0;
3449    TSQLRow    *row    = 0;    TSQLRow    *row    = 0;
3450    //    //
   UInt_t timelim = 0;  
   TDatime *time = new TDatime();  
   //  
3451    stringstream oss;    stringstream oss;
3452    oss.str("");    oss.str("");
3453    //    //
   //  
   //  
   if ( olderthan < 0 ){  
     if ( IsDebug() ) printf(" Do not validate runs \n");  
     return(1);  
   };  
   //  
   // timelim = now - olderthan  
   //  
   time->Set();  
   timelim =  (UInt_t)time->Convert() - olderthan;  
   time->Set(timelim,false);  
   
   printf("Validate runs older than %s \n",time->AsSQLString());  
   
3454    // =======================================================    // =======================================================
3455    // validate runs by checking missing calibrations    // validate runs by checking missing calibrations
3456    // =======================================================    // =======================================================
# Line 3277  Int_t PamelaDBOperations::ValidateRuns() Line 3460  Int_t PamelaDBOperations::ValidateRuns()
3460    // 1) get the OBT of the last run inserted after clean-time limit    // 1) get the OBT of the last run inserted after clean-time limit
3461    // --------------------------------------------------------------    // --------------------------------------------------------------
3462    oss.str("");    oss.str("");
3463    oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << time->AsSQLString()    oss << " SELECT * FROM GL_RUN  WHERE INSERT_TIME <= '" << clean_time->AsSQLString()
3464        << "' ORDER BY INSERT_TIME DESC LIMIT 1;";            << "' ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3465    printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" Get start validation-time: query is \n %s \n",oss.str().c_str());
3466    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
3467    if ( !result ) throw -4;    if ( !result ) throw -4;
3468    if ( !result->GetRowCount() ) {    if ( !result->GetRowCount() ) {
# Line 3288  Int_t PamelaDBOperations::ValidateRuns() Line 3471  Int_t PamelaDBOperations::ValidateRuns()
3471    }else{    }else{
3472          row = result->Next();          row = result->Next();
3473          t_start = (UInt_t)atoll(row->GetField(4));          t_start = (UInt_t)atoll(row->GetField(4));
         printf("t_start %i\n",t_start);  
3474    };      };  
3475    // --------------------------------------------------------------    // --------------------------------------------------------------
3476    // 2) get the OBT of the last validated run    // 2) get the OBT of the last validated run
# Line 3296  Int_t PamelaDBOperations::ValidateRuns() Line 3478  Int_t PamelaDBOperations::ValidateRuns()
3478    oss.str("");    oss.str("");
3479    oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start    oss << " SELECT * FROM GL_RUN  WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start
3480        <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";        <<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;";
3481    printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());    if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str());
3482    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
3483    if ( !result ) throw -4;    if ( !result ) throw -4;
3484    if ( result->GetRowCount() ){    if ( result->GetRowCount() ){
3485            row = result->Next();            row = result->Next();
3486            t_stop = (UInt_t)atoll(row->GetField(4));            t_stop = (UInt_t)atoll(row->GetField(4));
3487    };    };
3488    printf("t_stop %i\n",t_stop);    if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);
3489    // --------------------------------------------------------------    // --------------------------------------------------------------
3490    // now retrieves runs to be validated    // now retrieves runs to be validated
3491    // --------------------------------------------------------------    // --------------------------------------------------------------
# Line 3312  Int_t PamelaDBOperations::ValidateRuns() Line 3494  Int_t PamelaDBOperations::ValidateRuns()
3494    oss << " AND RUNHEADER_TIME >="<< t_stop;    oss << " AND RUNHEADER_TIME >="<< t_stop;
3495    oss << " ORDER BY RUNHEADER_TIME DESC;";    oss << " ORDER BY RUNHEADER_TIME DESC;";
3496  //  if ( IsDebug() )  //  if ( IsDebug() )
3497    printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());    if ( IsDebug() )printf(" Check runs for validation: query is \n %s \n",oss.str().c_str());
3498    result = conn->Query(oss.str().c_str());    result = conn->Query(oss.str().c_str());
3499    if ( !result ) throw -4;    if ( !result ) throw -4;
3500    if ( !result->GetRowCount() ) printf(" No runs to validate \n");    if ( !result->GetRowCount() ) printf(" No runs to validate \n");
3501    printf("------------------------------------------------------------------------------- \n");  //  printf("------------------------------------------------------------------------------- \n");
3502        
3503    Int_t nrow = 0;    Int_t nrow = 0;
3504    GL_RUN* this_run = new GL_RUN();    GL_RUN* this_run = new GL_RUN();
3505    GL_RUN* next_run = new GL_RUN();    GL_RUN* next_run = new GL_RUN();
3506    Int_t   nseq_max = 100;    Int_t   nseq_max = 1000;
3507  //  UInt_t* sequence = new UInt_t[100];  //  UInt_t* sequence = new UInt_t[100];
3508    vector<UInt_t> sequence(nseq_max);    vector<UInt_t> sequence(nseq_max);
3509    Int_t   nseq = 0;    Int_t   nseq = 0;
# Line 3344  Int_t PamelaDBOperations::ValidateRuns() Line 3526  Int_t PamelaDBOperations::ValidateRuns()
3526            //get run info            //get run info
3527            //------------            //------------
3528            this_run->Set_GL_RUN(row);            this_run->Set_GL_RUN(row);
3529                                
           printf(" RUN ID %i --- TRK_CALIB_USED %i --- RM_ACQ_AFTER_CALIB %i --- TIME %i %i \n",this_run->ID,this_run->TRK_CALIB_USED,this_run->RM_ACQ_AFTER_CALIB,this_run->RUNHEADER_TIME, this_run->RUNTRAILER_TIME);  
             
3530            Bool_t this_BAD = false;            Bool_t this_BAD = false;
3531            if(this_run->GetTRK_CALIB_USED() == 1 || this_run->GetTRK_CALIB_USED() == 2) this_ONLINE = true;            if(this_run->GetTRK_CALIB_USED() == 1 || this_run->GetTRK_CALIB_USED() == 2) this_ONLINE = true;
3532            else if (this_run->GetTRK_CALIB_USED() == 104)                          this_ONLINE = false;            else if (this_run->GetTRK_CALIB_USED() == 104)                          this_ONLINE = false;
3533            else{            else{
3534                    printf("Missing or corrupted header!! \n");  //                printf("Missing or corrupted header!! \n");
3535                    this_ONLINE = false;                    this_ONLINE = false;
3536                    this_BAD = true;                    this_BAD = true;
3537            };            };
# Line 3375  Int_t PamelaDBOperations::ValidateRuns() Line 3555  Int_t PamelaDBOperations::ValidateRuns()
3555                                                        
3556                            if( interval >= 60 )CHECK = true;                     //more than 60 s => there might be a calibration                            if( interval >= 60 )CHECK = true;                     //more than 60 s => there might be a calibration
3557                                                        
 /*                        if( CHECK && !next_run->GetRM_ACQ_AFTER_CALIB() )  
                                   printf(" ValidateRuns ***WARNING*** : DT = %i but RM_ACQ_AFTER_CALIB = %i \n",(t2-t1),next_run->GetRM_ACQ_AFTER_CALIB());*/  
                                                     
                             
3558                            if( !CHECK && this_run->VALIDATION ){                            if( !CHECK && this_run->VALIDATION ){
3559                                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);                                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true);
3560                                    nseq=0;                                    nseq=0;
# Line 3387  Int_t PamelaDBOperations::ValidateRuns() Line 3563  Int_t PamelaDBOperations::ValidateRuns()
3563                    }else if( !this_ONLINE && next_ONLINE) {              // this: DEFAULT + next:ON-LINE                    }else if( !this_ONLINE && next_ONLINE) {              // this: DEFAULT + next:ON-LINE
3564                                                        
3565                            CHECK = true;                            CHECK = true;
 //                        if( interval < 60 ) printf(" ValidateRuns ***WARNING*** : kkkkkk DT = %i \n",interval);  
3566    
3567                    }else if( !next_ONLINE ){                                             // this:ANY + next:DEFAULT                    }else if( !next_ONLINE ){                                             // this:ANY + next:DEFAULT
3568                                                        
# Line 3401  Int_t PamelaDBOperations::ValidateRuns() Line 3576  Int_t PamelaDBOperations::ValidateRuns()
3576            //----------------------------            //----------------------------
3577            if( CHECK ){            if( CHECK ){
3578                    // check if calibration exists                    // check if calibration exists
3579                    printf("DT %i ===> CHECK Missing calibration\n",interval);                    if ( IsDebug() )printf("DT %i ===> CHECK Missing calibration\n",interval);
3580                    Bool_t MISSING = MissingTRK_CALIB(t1,t2);                    Bool_t MISSING = MissingTRK_CALIB(t1,t2);
3581                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],!MISSING);                    for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],!MISSING);
3582                    nseq=0;                    nseq=0;
3583            } else printf("DT %i\n",interval);            };
3584            //--------------            //--------------
3585            //store run info            //store run info
3586            //--------------            //--------------
# Line 3418  Int_t PamelaDBOperations::ValidateRuns() Line 3593  Int_t PamelaDBOperations::ValidateRuns()
3593                    }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);
3594            };            };
3595                        
3596              if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID);
3597            nrow++;            nrow++;
3598                        
3599    };    };
# Line 3446  Bool_t PamelaDBOperations::MissingTRK_CA Line 3622  Bool_t PamelaDBOperations::MissingTRK_CA
3622          // which should be equal to the time between ascending-nodes.          // which should be equal to the time between ascending-nodes.
3623          //==============================================================          //==============================================================
3624          if ( t2 - trkcalib->FROM_TIME > 5700) {          if ( t2 - trkcalib->FROM_TIME > 5700) {
3625                  printf("Long time to previous calib %i :-( \n",t2 - trkcalib->FROM_TIME);                  if ( IsDebug() )printf("Long time between calib and run start %i :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);
3626          //==============================================================          //==============================================================
3627          // there might be a missing calibration, due to:          // there might be a missing calibration, due to:
3628          // - MM full          // - MM full
# Line 3467  Bool_t PamelaDBOperations::MissingTRK_CA Line 3643  Bool_t PamelaDBOperations::MissingTRK_CA
3643          // it is enough to say that there are no missing calibrations          // it is enough to say that there are no missing calibrations
3644          //==============================================================          //==============================================================
3645          // the long time interval bewteen runs might be due to download          // the long time interval bewteen runs might be due to download
3646          printf("Short time to previous calib %i :-) \n",t2 - trkcalib->FROM_TIME);          if ( IsDebug() )printf("Short time between calib and run start %i :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
3647          return(false);          return(false);
3648                    
3649  };  };
# Line 3480  Int_t PamelaDBOperations::assignVALIDATI Line 3656  Int_t PamelaDBOperations::assignVALIDATI
3656          TSQLResult *result = 0;          TSQLResult *result = 0;
3657          stringstream oss;          stringstream oss;
3658          oss.str("");          oss.str("");
3659          oss << " UPDATE GL_RUN SET VALIDATION="<< (UInt_t)validation <<" WHERE ID= " << idrun << "';";          oss << " UPDATE GL_RUN SET VALIDATION="<< (UInt_t)validation <<" WHERE ID= " << idrun << ";";
3660          //          //
3661  //      if ( IsDebug() )  //      if ( IsDebug() )
3662          printf(" Set VALIDATION = %i for run %i \n",validation,idrun);  //      printf(" Set VALIDATION = %i for run %i \n",validation,idrun);
3663  //      printf(" Set VALIDATION = %i for run %i : query is \n %s \n",validation,idrun,oss.str().c_str());          if ( IsDebug() )printf(" Query: %s \n",oss.str().c_str());
3664  //      result = conn->Query(oss.str().c_str());          result = conn->Query(oss.str().c_str());
3665  //      if ( !result ) throw -4;          if ( !result ) throw -4;
3666          return(0);          return(0);
3667  }  }
3668    
3669    
3670    
3671    // Insert TLEs from file tlefilename in the table GL_TLE in the db
3672    // opened by conn, sorting them by date from older to newer, if each
3673    // TLE has not been alread inserted.
3674    Int_t PamelaDBOperations::populateTLE()//(TSQLServer *conn, char *tleFile)
3675    {
3676      fstream tlefile(tlefilename, ios::in);
3677    
3678      if ( !tlefile ) throw -7;
3679    
3680      vector<cTle*> ctles;
3681      vector<cTle*>::iterator iter;
3682      int present = 0;
3683    
3684      // Get three lines from tlefile, create a cTle object and put it
3685      // into ctles
3686      while(1) {
3687        cTle *tlef;
3688        string str1, str2, str3;
3689    
3690        getline(tlefile, str1);
3691        if(tlefile.eof()) break;
3692    
3693        getline(tlefile, str2);
3694        if(tlefile.eof()) break;
3695    
3696        getline(tlefile, str3);
3697        if(tlefile.eof()) break;
3698    
3699        // We now have three good lines for a cTle.
3700        tlef = new cTle(str1, str2, str3);
3701        ctles.push_back(tlef);
3702      }
3703    
3704      tlefile.close();
3705    
3706      // Sort by date
3707      sort(ctles.begin(), ctles.end(), compTLE);
3708    
3709      // Now we insert each TLE into the db
3710      for(iter = ctles.begin(); iter != ctles.end(); iter++) {
3711        cTle *tle = *iter;
3712    
3713        // Do nothing if it's already present in the db.  Just increase
3714        // the counter present.
3715        if (! isTlePresent(tle))
3716          {
3717            int status = insertTle(tle);
3718    
3719            // Insert query failed.  Return 1.
3720            if(status == EXIT_FAILURE) {
3721              
3722              if( IsDebug() ) {
3723                cerr << "Error: inserting TLE:" << endl
3724                     << tle->getName() << endl
3725                     << tle->getLine1() << endl
3726                     << tle->getLine2() << endl;
3727              }
3728    
3729              throw -4;
3730              return 1;
3731            }
3732    
3733          }
3734        else
3735          present++;
3736    
3737      }
3738    
3739      int inserted = ctles.size() - present;  // Number of inserted TLE.
3740      if ( IsDebug() )
3741        cout << "\nProcessed TLEs ranging from " << getTleDatetime(ctles[0]) << " to " << getTleDatetime(ctles[ctles.size()-1]) << "." << endl
3742             << inserted << " newly inserted TLEs out of " << ctles.size() << " processed." << endl;
3743    
3744      ctles.clear();
3745    
3746    
3747      // Return 2 if no new TLE has been inserted.  0 otherwise.
3748      if(! inserted ) return 2;
3749      return 0;
3750    }
3751    
3752    
3753    // Insert tle in the table GL_TLE using the connection conn.
3754    int PamelaDBOperations::insertTle(cTle *tle)
3755    {
3756      stringstream oss;
3757      TSQLResult *result = 0;
3758    
3759      oss.str("");
3760      oss << " INSERT INTO GL_TLE (TLE1, TLE2, TLE3, FROM_TIME)"
3761          << " VALUES ( '"
3762          << tle->getName() << "', '"
3763          << tle->getLine1() << "', '"
3764          << tle->getLine2() << "', '"
3765          << getTleDatetime(tle) << "')";
3766    
3767      //  cout << oss.str().c_str() << endl;
3768      result = conn->Query(oss.str().c_str());
3769      if (result == NULL)
3770        return EXIT_FAILURE;
3771    
3772      return EXIT_SUCCESS;
3773    }
3774    
3775    
3776    // Return whether tle is already in the db connected by conn.
3777    bool PamelaDBOperations::isTlePresent(cTle *tle)
3778    {
3779      stringstream oss;
3780      TSQLResult *result = 0;
3781    
3782      oss.str("");
3783      oss << "SELECT * FROM GL_TLE WHERE FROM_TIME = '"
3784          << getTleDatetime(tle) << "'";
3785    
3786      result = conn->Query(oss.str().c_str());
3787      if (result == NULL) throw -4;
3788    
3789      if (result->GetRowCount())
3790        return true;
3791      else
3792        return false;
3793    }
3794    
3795    
3796    // Return whether the first TLE is dated early than the second
3797    bool compTLE (cTle *tle1, cTle *tle2)
3798    {
3799      return getTleJulian(tle1) < getTleJulian(tle2);
3800    }
3801    
3802    
3803    // Return the date of the tle using the format (year-2000)*1e3 +
3804    // julian day.  e.g. 6365 is the 31th Dec 2006.
3805    // It does *not* return a cJulian date.
3806    float getTleJulian(cTle *tle) {
3807      return tle->getField(cTle::FLD_EPOCHYEAR)*1e3 + tle->getField(cTle::FLD_EPOCHDAY);
3808    }
3809    
3810    
3811    // Return a string like YYYY-MM-DD hh:mm:ss, usable for mysql datetime
3812    // format.
3813    string getTleDatetime(cTle *tle)
3814    {
3815      int year, mon, day, hh, mm, ss;
3816      double dom; // day of month (is double!)
3817      stringstream date; // date in datetime format
3818    
3819      // create a cJulian from the date in tle
3820      cJulian jdate = cJulian( 2000 + (int) tle->getField(cTle::FLD_EPOCHYEAR), tle->getField(cTle::FLD_EPOCHDAY));
3821    
3822      // get year, month, day of month
3823      jdate.getComponent(&year, &mon, &dom);
3824    
3825      // build a datetime YYYY-MM-DD hh:mm:ss
3826      date.str("");
3827      day = (int) floor(dom);
3828      hh = (int) floor( (dom - day) * 24);
3829      mm = (int) floor( ((dom - day) * 24 - hh) * 60);
3830      ss = (int) floor( ((((dom - day) * 24 - hh) * 60 - mm) * 60));
3831      //  ms = (int) floor( (((((dom - day) * 24 - hh) * 60 - mm) * 60) - ss) * 1000);
3832    
3833      date << year << "-" << mon << "-" << day << " " << hh << ":" << mm << ":" << ss;
3834    
3835      return date.str();
3836    }
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23