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

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

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

revision 1.12 by mocchiut, Fri Nov 28 10:57:44 2008 UTC revision 1.18 by mocchiut, Mon Dec 22 10:38:15 2008 UTC
# Line 101  PamelaDBOperations::PamelaDBOperations(T Line 101  PamelaDBOperations::PamelaDBOperations(T
101    if (INSERT_RAW) SetRawName(filerawname);    if (INSERT_RAW) SetRawName(filerawname);
102    //    //
103    INSERT_ROOT = !filerootname.IsNull();    INSERT_ROOT = !filerootname.IsNull();
104      if ( INSERT_ROOT ) this->SetRootName(filerootname);
105    this->SetOrbitNo(dwinput);    this->SetOrbitNo(dwinput);
106    //    //
107    this->SetID_RAW(0);    this->SetID_RAW(0);
# Line 134  void PamelaDBOperations::OpenL0File(TStr Line 135  void PamelaDBOperations::OpenL0File(TStr
135        stringstream qu;        stringstream qu;
136        TSQLResult *result = 0;        TSQLResult *result = 0;
137        TSQLResult *result2 = 0;        TSQLResult *result2 = 0;
138          TSQLResult *result3 = 0;
139        TSQLRow    *row    = 0;            TSQLRow    *row    = 0;    
140        TSQLRow    *row2    = 0;            TSQLRow    *row2    = 0;    
141          TSQLRow    *row3    = 0;    
142        TString chpath;        TString chpath;
143        TString chfile;        TString chfile;
144        UInt_t ridn = 0;        UInt_t ridn = 0;
# Line 160  void PamelaDBOperations::OpenL0File(TStr Line 163  void PamelaDBOperations::OpenL0File(TStr
163        };        };
164        delete result;        delete result;
165        qu.str("");        qu.str("");
166        qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL from ROOT_TABLE where ID_N=" << ridn << ";";        qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL,REAL_TIME_INIT from ROOT_TABLE where ID_N=" << ridn << ";";
167        if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());        if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
168        result2 = conn->Query(qu.str().c_str());        result2 = conn->Query(qu.str().c_str());
169        if ( result2 ){        if ( result2 ){
# Line 177  void PamelaDBOperations::OpenL0File(TStr Line 180  void PamelaDBOperations::OpenL0File(TStr
180            chresursts = (UInt_t)atoll(row2->GetField(4));            chresursts = (UInt_t)atoll(row2->GetField(4));
181            chboot = (UInt_t)atoll(row2->GetField(5));            chboot = (UInt_t)atoll(row2->GetField(5));
182            //            //
183              if ( !chboot ){
184                qu.str("");
185                qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<10000 AND (REAL_TIME_INIT-"<< row2->GetField(10) <<")<7200 group by BOOT_NUMBER order by BOOT_NUMBER asc limit 1;";
186                if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
187                result3 = conn->Query(qu.str().c_str());
188                if ( result3 ){
189                  row3 = result3->Next();
190                  if ( row3 ){
191                    chboot = (UInt_t)atoll(row3->GetField(0));
192                    if ( debug ) printf(" Found boot_number = %u \n",chboot);
193                  } else {
194                    if ( debug ) printf(" AGH CANNOT DETERMINE THE BOOT NUMBER... \n");
195                  };
196                };
197              };
198              //
199            chpktinit = (UInt_t)atoll(row2->GetField(6));            chpktinit = (UInt_t)atoll(row2->GetField(6));
200            chpktfinal = (UInt_t)atoll(row2->GetField(7));            chpktfinal = (UInt_t)atoll(row2->GetField(7));
201            chobtinit = (UInt_t)atoll(row2->GetField(8));            chobtinit = (UInt_t)atoll(row2->GetField(8));
# Line 189  void PamelaDBOperations::OpenL0File(TStr Line 208  void PamelaDBOperations::OpenL0File(TStr
208          throw -85;          throw -85;
209        };        };
210        filerootname = chpath + chfile;// + ".root";        filerootname = chpath + chfile;// + ".root";
       if ( debug ) printf(" chewbacca: filename is %s \n",filerootname.Data());  
211      };      };
212      this->SetRootName(filerootname);      this->SetRootName(filerootname);
     file = TFile::Open(this->GetRootName().Data());  
   } else {  
     this->SetRootName("");  
213    };    };
214      if ( debug ) printf(" Filename is %s \n",filerootname.Data());
215      file = TFile::Open(this->GetRootName().Data());
216      //  } else {
217      //    this->SetRootName("");
218      //    this->SetRootName(filerootname);
219      //  };
220  }  }
221    
222  //  //
# Line 395  void PamelaDBOperations::SetOrbitNo(UInt Line 416  void PamelaDBOperations::SetOrbitNo(UInt
416      return;      return;
417    };    };
418    //    //
419    TString name = this->GetRootFile();    if ( !chewbacca ){
420    Int_t nlength = name.Length();      TString name = this->GetRootFile();
421    if ( nlength < 5 ) return;      Int_t nlength = name.Length();
422    TString dwo = 0;      if ( nlength < 5 ){
423    for (Int_t i = 0; i<5; i++){        if ( IsDebug() ) printf(" Agh problems determining the orbit number! name = %s \n",name.Data());
424      dwo.Append(name[i],1);        return;
425    };      };
426    if ( dwo.IsDigit() ){      TString dwo = 0;
427      dworbit = (UInt_t)dwo.Atoi();      for (Int_t i = 0; i<5; i++){
   } else {  
     dwo="";  
     for (Int_t i = 8; i<13; i++){  
428        dwo.Append(name[i],1);        dwo.Append(name[i],1);
429      };          };
430      if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();      if ( dwo.IsDigit() ){
431          dworbit = (UInt_t)dwo.Atoi();
432        } else {
433          dwo="";
434          for (Int_t i = 8; i<13; i++){
435            dwo.Append(name[i],1);
436          };    
437          if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
438        };
439        if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());
440    };    };
   if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());  
441    return;    return;
442  };  };
443    
# Line 1577  Int_t PamelaDBOperations::insertPamelaGL Line 1603  Int_t PamelaDBOperations::insertPamelaGL
1603      result = conn->Query(oss.str().c_str());      result = conn->Query(oss.str().c_str());
1604      if ( !result ) throw -10;      if ( !result ) throw -10;
1605      row = result->Next();      row = result->Next();
1606        if ( !row ) throw -92;
1607      idresof = (UInt_t)atoll(row->GetField(0));      idresof = (UInt_t)atoll(row->GetField(0));
1608      existsts = true;      existsts = true;
1609      goto eout;      goto eout;
# Line 2239  Int_t PamelaDBOperations::insertPamelaRU Line 2266  Int_t PamelaDBOperations::insertPamelaRU
2266          ptht = cod->Get(pctp->RunHeader) - 1;          ptht = cod->Get(pctp->RunHeader) - 1;
2267          //      evbeft = cod->Get(pctp->Physics);          //      evbeft = cod->Get(pctp->Physics);
2268          tcod = (UInt_t)cod->Get(pctp->Physics);          tcod = (UInt_t)cod->Get(pctp->Physics);
2269          if ( !tcod ) tcod = 1;          if ( !tcod ) tcod = 1;    // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
2270          evbeft = TMath::Min(upperentry,(tcod-1));          evbeft = TMath::Min(upperentry,(tcod-1));
2271          if ( debug ) printf(" Loop in runtrailers, evbeft is %u upperentry %u cod->getetc %u \n",evbeft,upperentry,cod->Get(pctp->Physics));          if ( debug ) printf(" Loop in runtrailers, evbeft is %u upperentry %u cod->getetc %u \n",evbeft,upperentry,cod->Get(pctp->Physics));
2272          //          //
# Line 2278  Int_t PamelaDBOperations::insertPamelaRU Line 2305  Int_t PamelaDBOperations::insertPamelaRU
2305            obth = phh->GetOrbitalTime();            obth = phh->GetOrbitalTime();
2306            cod = ehh->GetCounter();            cod = ehh->GetCounter();
2307            tcod = (UInt_t)cod->Get(pctp->Physics);            tcod = (UInt_t)cod->Get(pctp->Physics);
2308            if ( !tcod ) tcod = 1;            if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
2309            evbefh = TMath::Max(chminentry,(tcod-1));            evbefh = TMath::Max(chminentry,(tcod-1));
2310            //            //
2311            if ( PKT(pkth) >= PKT(pktfirst) && PKT(pkth) <= upperpkt ){            if ( PKT(pkth) >= PKT(pktfirst) && PKT(pkth) <= upperpkt ){
# Line 2340  Int_t PamelaDBOperations::insertPamelaRU Line 2367  Int_t PamelaDBOperations::insertPamelaRU
2367                cod = ehh->GetCounter();                cod = ehh->GetCounter();
2368                //            evbeft = cod->Get(pctp->Physics);                //            evbeft = cod->Get(pctp->Physics);
2369                tcod = (UInt_t)cod->Get(pctp->Physics);                tcod = (UInt_t)cod->Get(pctp->Physics);
2370                if ( !tcod ) tcod = 1;                if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT
2371                evbeft = TMath::Min(upperentry,(tcod-1));                evbeft = TMath::Min(upperentry,(tcod-1));
2372                rh->GetEntry(pth);                rh->GetEntry(pth);
2373                phh = ehh->GetPscuHeader();                phh = ehh->GetPscuHeader();
# Line 2382  Int_t PamelaDBOperations::insertPamelaRU Line 2409  Int_t PamelaDBOperations::insertPamelaRU
2409            //    if ( conptt+1 == (nrtev+nrtbef )){              //    if ( conptt+1 == (nrtev+nrtbef )){  
2410            if ( IsDebug() ) printf(" >>>>>>>>>>> %i %u %i %u \n",ptht,rhev,nrtev,conptt);            if ( IsDebug() ) printf(" >>>>>>>>>>> %i %u %i %u \n",ptht,rhev,nrtev,conptt);
2411            ptht++;            ptht++;
2412            if ( ptht < rhev ){            if ( ptht < rhev ){
2413              rh->GetEntry(ptht);              rh->GetEntry(ptht);
2414              //pth++;              //pth++;
2415              //if ( pth < rhev ){              //if ( pth < rhev ){
# Line 2410  Int_t PamelaDBOperations::insertPamelaRU Line 2437  Int_t PamelaDBOperations::insertPamelaRU
2437              if ( (UInt_t)evbeft < upperentry-1 && upperentry>0 ){              if ( (UInt_t)evbeft < upperentry-1 && upperentry>0 ){
2438                if ( IsDebug() ) printf(" Piece of run at the end of the file with NO RUNHEADER! evbeft %u upperentry-1 %u \n",(UInt_t)evbeft,upperentry-1);                if ( IsDebug() ) printf(" Piece of run at the end of the file with NO RUNHEADER! evbeft %u upperentry-1 %u \n",(UInt_t)evbeft,upperentry-1);
2439                //                //
2440                if ( (ptt-1) < 0 ) throw -15; // should never arrive here!                //              if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
2441                  if ( (rtev-1) < 0 || ptt < 0 ) throw -15; // should never arrive here!
2442                //              rt->GetEntry(ptt-1);                //              rt->GetEntry(ptt-1);
2443                rt->GetEntry(rtev-1);                rt->GetEntry(rtev-1);
2444                cod = eht->GetCounter();                cod = eht->GetCounter();
# Line 2645  void PamelaDBOperations::HandleRunFragme Line 2673  void PamelaDBOperations::HandleRunFragme
2673      //      //
2674      // we have now the good first piece of a run, fill the glrun object      // we have now the good first piece of a run, fill the glrun object
2675      //      //
2676      if ( rhfirstev != firstev && !mishead ) mishead = true;      if ( firstev != (lastev+1) ){ // could be a problem sometimes (?)
2677      if ( rtlastev != lastev && !mistrail ) mistrail = true;        if ( rhfirstev != firstev && !mishead ) mishead = true;
2678          if ( rtlastev != lastev && !mistrail ) mistrail = true;
2679        };
2680      //      //
2681        if ( IsDebug() ) printf(" bhere firstev is %i lastev is %i \n",firstev,lastev);
2682      this->FillClass(mishead,mistrail,firstev,lastev);      this->FillClass(mishead,mistrail,firstev,lastev);
2683        if ( IsDebug() ) printf(" chere firstev is %i lastev is %i \n",firstev,lastev);
2684      //      //
2685      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");
2686      if ( IsDebug() ) printf(" C THIS RUN: RUNHEADER_OBT %u RUNTRAILER_OBT %u RUNHEADER_PKT %u RUNTRAILER_PKT %u \n", glrun->GetRUNHEADER_OBT(),glrun->GetRUNTRAILER_OBT(),glrun->GetRUNHEADER_PKT(),glrun->GetRUNTRAILER_PKT());      if ( IsDebug() ) printf(" C THIS RUN: RUNHEADER_OBT %u RUNTRAILER_OBT %u RUNHEADER_PKT %u RUNTRAILER_PKT %u \n", glrun->GetRUNHEADER_OBT(),glrun->GetRUNTRAILER_OBT(),glrun->GetRUNHEADER_PKT(),glrun->GetRUNTRAILER_PKT());
# Line 2743  void PamelaDBOperations::HandleRunFragme Line 2775  void PamelaDBOperations::HandleRunFragme
2775          //          //
2776          found = false; // default value          found = false; // default value
2777          //          //
2778          if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");          if ( IsDebug() ) printf(" A Found a possible candidate, checking if it is the good one... \n");
2779          //          //
2780          // if we have both runheader and runtrailer we can check with pkt_counter:          // if we have both runheader and runtrailer we can check with pkt_counter:
2781          //          //
# Line 2759  void PamelaDBOperations::HandleRunFragme Line 2791  void PamelaDBOperations::HandleRunFragme
2791              if ( IsDebug() ) printf(" FOUND!!! check %llu pktt %llu \n",chkpkt,pktt);              if ( IsDebug() ) printf(" FOUND!!! check %llu pktt %llu \n",chkpkt,pktt);
2792              //              //
2793              found = true;              found = true;
2794                if ( IsDebug() ) printf(" where firstev is %i lastev is %i \n",firstev,lastev);
2795              //              //
2796            } else {            } else {
2797              //              //
# Line 2824  void PamelaDBOperations::HandleRunFragme Line 2857  void PamelaDBOperations::HandleRunFragme
2857          oss << " ID="<<row->GetField(0)<<";";          oss << " ID="<<row->GetField(0)<<";";
2858          //          //
2859          glrun1->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); // here we have runheader infos          glrun1->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); // here we have runheader infos
2860            if ( IsDebug() ) printf(" there firstev is %i lastev is %i \n",firstev,lastev);
2861          //          //
2862          // merge infos          // merge infos
2863          //          //
# Line 2854  void PamelaDBOperations::HandleRunFragme Line 2888  void PamelaDBOperations::HandleRunFragme
2888          glrun1->SetRUNTRAILER_OBT(glrun->GetRUNTRAILER_OBT());          glrun1->SetRUNTRAILER_OBT(glrun->GetRUNTRAILER_OBT());
2889          glrun1->SetRUNTRAILER_PKT(glrun->GetRUNTRAILER_PKT());          glrun1->SetRUNTRAILER_PKT(glrun->GetRUNTRAILER_PKT());
2890          //          //
2891            if ( IsDebug() ) printf(" here firstev is %i lastev is %i \n",firstev,lastev);
2892            //
2893          glrun->SetEV_FROM(firstev);          glrun->SetEV_FROM(firstev);
2894          glrun->SetNEVENTS(lastev-firstev+1);          glrun->SetNEVENTS(lastev-firstev+1);
2895          //          //
# Line 2946  void PamelaDBOperations::HandleRunFragme Line 2982  void PamelaDBOperations::HandleRunFragme
2982          //          //
2983          found = false; // default value          found = false; // default value
2984          //          //
2985          if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");          if ( IsDebug() ) printf(" B Found a possible candidate, checking if it is the good one... \n");
2986          //          //
2987          // if we have both runheader and runtrailer we can check with pkt_counter:          // if we have both runheader and runtrailer we can check with pkt_counter:
2988          //          //
# Line 3175  void PamelaDBOperations::HandleMissingHo Line 3211  void PamelaDBOperations::HandleMissingHo
3211        glrun->SetID(this->AssignRunID());        glrun->SetID(this->AssignRunID());
3212        glrun->SetID_RUN_FRAG(0);        glrun->SetID_RUN_FRAG(0);
3213        glrun->Fill_GL_RUN(conn);    // it'ok we arrive here only inside a file hence in the middle of the runs...        glrun->Fill_GL_RUN(conn);    // it'ok we arrive here only inside a file hence in the middle of the runs...
3214      };        };
3215      //      //
3216    };    };
3217    //    //
# Line 3219  Bool_t PamelaDBOperations::IsRunConsiste Line 3255  Bool_t PamelaDBOperations::IsRunConsiste
3255    if ( firstev == lastev+1 || lastev == firstev ) { // no events inside the run!    if ( firstev == lastev+1 || lastev == firstev ) { // no events inside the run!
3256      //if ( firstev <= lastev+1 ) { // no events inside the run!      //if ( firstev <= lastev+1 ) { // no events inside the run!
3257      if ( IsDebug() ) printf(" Checking but no events in the run! \n");      if ( IsDebug() ) printf(" Checking but no events in the run! \n");
3258        firstev = lastev+1; // this is necessary for files with no Physics entries, should have no influence on other files
3259      // return true is correct      // return true is correct
3260      return(true);      return(true);
3261      //      //
# Line 3474  Bool_t PamelaDBOperations::IsRunConsiste Line 3511  Bool_t PamelaDBOperations::IsRunConsiste
3511                  //                  //
3512                  found = false; // default value                  found = false; // default value
3513                  //                  //
3514                  if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n");                  if ( IsDebug() ) printf(" C Found a possible candidate, checking if it is the good one... \n");
3515                  //                  //
3516                  // if we have both runheader and runtrailer we can check with pkt_counter:                  // if we have both runheader and runtrailer we can check with pkt_counter:
3517                  //                  //
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23