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

Diff of /YodaProfiler/src/GLTables.cpp

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

revision 1.15 by mocchiut, Tue Oct 31 15:36:04 2006 UTC revision 1.18 by mocchiut, Thu Nov 9 16:47:47 2006 UTC
# Line 10  Line 10 
10  //  //
11  #include <TFile.h>  #include <TFile.h>
12  #include <TTree.h>  #include <TTree.h>
13    #include <TTimeStamp.h>
14  #include <EventHeader.h>  #include <EventHeader.h>
15  #include <PscuHeader.h>  #include <PscuHeader.h>
16  //  //
# Line 153  GL_TIMESYNC::GL_TIMESYNC(){ Line 154  GL_TIMESYNC::GL_TIMESYNC(){
154    obtfirst = 0;    obtfirst = 0;
155    pktfirst = 0;    pktfirst = 0;
156    toffset  = 0;    toffset  = 0;
157    ID       = 0;    //  ID       = 0;
158    ID_RAW   = 0;    //  ID_RAW   = 0;
159    OBT0     = 0;    OBT0     = 0;
160    TIMESYNC = 0;    TIMESYNC = 0;
161    TYPE     = 0;    TYPE     = 0;
# Line 861  Int_t GL_RUN::Query_GL_RUN(UInt_t run, T Line 862  Int_t GL_RUN::Query_GL_RUN(UInt_t run, T
862        if( Row == NULL ) break;        if( Row == NULL ) break;
863  //        Set_GL_RUN(Row);  //        Set_GL_RUN(Row);
864        for( t = 0; t < pResult->GetFieldCount(); t++){        for( t = 0; t < pResult->GetFieldCount(); t++){
865                    if (t== 0) ID                = (UInt_t)atoll(Row->GetField(t));          if (t== 0) ID                = (UInt_t)atoll(Row->GetField(t));
866                    if (t== 1) ID_RUN_FRAG       = (UInt_t)atoll(Row->GetField(t));          if (t== 1) ID_RUN_FRAG       = (UInt_t)atoll(Row->GetField(t));
867                    if (t== 2) ID_ROOT_L0        = (UInt_t)atoll(Row->GetField(t));          if (t== 2) ID_ROOT_L0        = (UInt_t)atoll(Row->GetField(t));
868                    if (t== 3) ID_ROOT_L2        = (UInt_t)atoll(Row->GetField(t));          if (t== 3) ID_ROOT_L2        = (UInt_t)atoll(Row->GetField(t));
869                    if (t== 4) RUNHEADER_TIME    = (UInt_t)atoll(Row->GetField(t));          if (t== 4) RUNHEADER_TIME    = (UInt_t)atoll(Row->GetField(t));
870                    if (t== 5) RUNTRAILER_TIME   = (UInt_t)atoll(Row->GetField(t));          if (t== 5) RUNTRAILER_TIME   = (UInt_t)atoll(Row->GetField(t));
871                    if (t== 6) RUNHEADER_OBT     = (UInt_t)atoll(Row->GetField(t));          if (t== 6) RUNHEADER_OBT     = (UInt_t)atoll(Row->GetField(t));
872                    if (t== 7) RUNTRAILER_OBT    = (UInt_t)atoll(Row->GetField(t));          if (t== 7) RUNTRAILER_OBT    = (UInt_t)atoll(Row->GetField(t));
873                    if (t== 8) RUNHEADER_PKT     = (UInt_t)atoll(Row->GetField(t));          if (t== 8) RUNHEADER_PKT     = (UInt_t)atoll(Row->GetField(t));
874                    if (t== 9) RUNTRAILER_PKT    = (UInt_t)atoll(Row->GetField(t));          if (t== 9) RUNTRAILER_PKT    = (UInt_t)atoll(Row->GetField(t));
875                    if (t==10) EV_FROM           = (UInt_t)atoll(Row->GetField(t));          if (t==10) EV_FROM           = (UInt_t)atoll(Row->GetField(t));
876                    if (t==11) EV_TO             = (UInt_t)atoll(Row->GetField(t));          if (t==11) EV_TO             = (UInt_t)atoll(Row->GetField(t));
877                    if (t==12) NEVENTS           = (UInt_t)atoll(Row->GetField(t));          if (t==12) NEVENTS           = (UInt_t)atoll(Row->GetField(t));
878                    if (t==13) LAST_TIMESYNC     = (UInt_t)atoll(Row->GetField(t));          if (t==13) LAST_TIMESYNC     = (UInt_t)atoll(Row->GetField(t));
879                    if (t==14) OBT_TIMESYNC      = (UInt_t)atoll(Row->GetField(t));          if (t==14) OBT_TIMESYNC      = (UInt_t)atoll(Row->GetField(t));
880                    if (t==15) COMPILATIONTIMESTAMP = (UInt_t)atoll(Row->GetField(t));          if (t==15) COMPILATIONTIMESTAMP = (UInt_t)atoll(Row->GetField(t));
881                    if (t==16) FAV_WRK_SCHEDULE  = (UInt_t)atoll(Row->GetField(t));          if (t==16) FAV_WRK_SCHEDULE  = (UInt_t)atoll(Row->GetField(t));
882                    if (t==17) EFF_WRK_SCHEDULE  = (UInt_t)atoll(Row->GetField(t));          if (t==17) EFF_WRK_SCHEDULE  = (UInt_t)atoll(Row->GetField(t));
883                    if (t==18) PRH_VAR_TRG_MODE_A= (UInt_t)atoll(Row->GetField(t));          if (t==18) PRH_VAR_TRG_MODE_A= (UInt_t)atoll(Row->GetField(t));
884                    if (t==19) PRH_VAR_TRG_MODE_B= (UInt_t)atoll(Row->GetField(t));          if (t==19) PRH_VAR_TRG_MODE_B= (UInt_t)atoll(Row->GetField(t));
885                    if (t==20) ACQ_BUILD_INFO    = (UInt_t)atoll(Row->GetField(t));          if (t==20) ACQ_BUILD_INFO    = (UInt_t)atoll(Row->GetField(t));
886                    if (t==21) ACQ_VAR_INFO      = (UInt_t)atoll(Row->GetField(t));          if (t==21) ACQ_VAR_INFO      = (UInt_t)atoll(Row->GetField(t));
887                    if (t==22) RM_ACQ_AFTER_CALIB= (UInt_t)atoll(Row->GetField(t));          if (t==22) RM_ACQ_AFTER_CALIB= (UInt_t)atoll(Row->GetField(t));
888                    if (t==23) RM_ACQ_SETTING_MODE = (UInt_t)atoll(Row->GetField(t));          if (t==23) RM_ACQ_SETTING_MODE = (UInt_t)atoll(Row->GetField(t));
889                    if (t==24) PKT_COUNTER       = (UInt_t)atoll(Row->GetField(t));          if (t==24) PKT_COUNTER       = (UInt_t)atoll(Row->GetField(t));
890                    if (t==25) PKT_READY_COUNTER = (UInt_t)atoll(Row->GetField(t));          if (t==25) PKT_READY_COUNTER = (UInt_t)atoll(Row->GetField(t));
891                    if (t==26) TRK_CALIB_USED    = (UInt_t)atoll(Row->GetField(t));          if (t==26) TRK_CALIB_USED    = (UInt_t)atoll(Row->GetField(t));
892                    if (t==27) CAL_DSP_MASK      = (UInt_t)atoll(Row->GetField(t));          if (t==27) CAL_DSP_MASK      = (UInt_t)atoll(Row->GetField(t));
893                    if (t==28) BOOT_NUMBER       = (UInt_t)atoll(Row->GetField(t));          if (t==28) BOOT_NUMBER       = (UInt_t)atoll(Row->GetField(t));
894                    if (t==29) VALIDATION        = (UInt_t)atoll(Row->GetField(t));          if (t==29) VALIDATION        = (UInt_t)atoll(Row->GetField(t));
895        };        };
896    };    };
897    //  delete pResult;    //  delete pResult;
# Line 1303  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr Line 1304  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr
1304    //    //
1305    // look for Resurs offset    // look for Resurs offset
1306    //    //
1307    UInt_t t0 = 0;    T0 = 0;
1308    //    //
1309    //    //
1310    stringstream oss;    stringstream oss;
1311    oss.str("");    oss.str("");
1312    oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='"    oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE SPECIAL_FILE='"
1313        << rname.str().c_str() << "';";        << rname.str().c_str() << "';";
   //  if ( IsDebug() ) printf(" %s \n",oss.str().c_str());  
1314    pResult = dbc->Query(oss.str().c_str());    pResult = dbc->Query(oss.str().c_str());
1315    Row = pResult->Next();    Row = pResult->Next();
1316    //    //
# Line 1331  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr Line 1331  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr
1331      };          };    
1332      if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();      if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi();
1333    };    };
   // if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data());    
1334    //    //
1335    if ( !Row ){    if ( !Row ){
1336      oss.str("");      oss.str("");
1337      oss << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "      oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE) FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< "
1338          << dworbit << " order by FROM_ORBIT desc limit 1;";          << dworbit << " order by FROM_ORBIT desc limit 1;";
     //    if ( IsDebug() ) printf(" %s \n",oss.str().c_str());  
1339      pResult = dbc->Query(oss.str().c_str());      pResult = dbc->Query(oss.str().c_str());
1340      Row = pResult->Next();      Row = pResult->Next();
1341      if ( !Row ){      if ( !Row ){
# Line 1346  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr Line 1344  GL_TIMESYNC::GL_TIMESYNC(UInt_t id, TStr
1344      };      };
1345    };    };
1346    //    //
1347    t0 = (UInt_t)TDatime(Row->GetField(0)).Convert();    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);
1348    //    T0 = (UInt_t)tu.GetSec();
   //  myquery.str("");  
   //   myquery << "SELECT OFFSET_DATE FROM GL_RESURS_OFFSET WHERE FROM_ID_RAW<= "  
   //       << idraw << " AND TO_ID_RAW >= "  
   //       << idraw << ";";  
   //   pResult = dbc->Query(myquery.str().c_str());  
   //   if ( pResult ){  
   //     Row = pResult->Next();  
   //     if ( Row ){  
   //       //  
   //       t0 = (UInt_t)TDatime(Row->GetField(0)).Convert();  
   //       //  
   //     };  
   //   };  
   //  
1349    //    //
1350    // look for the correct timesync entry    // look for the correct timesync entry
1351    //    //
1352    myquery.str("");    myquery.str("");
1353    myquery << " SELECT OBT0,TIMESYNC FROM GL_TIMESYNC "    myquery << " SELECT OBT0,TIMESYNC,TYPE FROM GL_TIMESYNC "
1354        << " WHERE ID_RAW = " << idraw        << " WHERE ID_RAW = " << idraw
1355        << ";";        << ";";
1356    pResult = dbc->Query(myquery.str().c_str());    pResult = dbc->Query(myquery.str().c_str());
1357    if ( pResult ){    if ( pResult ){
1358      Row = pResult->Next();      Row = pResult->Next();
1359      if ( (Row != NULL) && ((UInt_t)atoll(Row->GetField(0)) > 0 ) ){      if ( (Row != NULL) && ((UInt_t)atoll(Row->GetField(0)) > 0 ) ){
1360        toffset = (UInt_t)atoll(Row->GetField(1)) - (UInt_t)(this->DBobt((UInt_t)atoll(Row->GetField(0)))/1000) + t0;        OBT0 = (UInt_t)atoll(Row->GetField(0));
1361          TIMESYNC = (UInt_t)atoll(Row->GetField(1));
1362          TYPE = (UInt_t)atoll(Row->GetField(2));
1363          toffset = (UInt_t)atoll(Row->GetField(1)) - (UInt_t)(this->DBobt((UInt_t)atoll(Row->GetField(0)))/1000) + T0;
1364      };      };
1365    };    };
1366    //    //
# Line 1392  UInt_t GL_TIMESYNC::DBabsTime(UInt_t OBT Line 1379  UInt_t GL_TIMESYNC::DBabsTime(UInt_t OBT
1379    //    //
1380  };  };
1381    
1382    /**
1383     *
1384     * Returns the Resurs time given the OBT needed to process inclination and orbital infos
1385     *
1386     */
1387    UInt_t GL_TIMESYNC::ResursTime(UInt_t OBT){  
1388      //
1389      return(((UInt_t)((this->DBobt(OBT)-this->DBobt(OBT0))/1000)+TIMESYNC));
1390      //
1391    };
1392    
1393    
1394  ULong64_t GL_TIMESYNC::DBobt(UInt_t obt){    ULong64_t GL_TIMESYNC::DBobt(UInt_t obt){  
1395    //    //
# Line 1435  UInt_t GL_TIMESYNC::DBpkt(UInt_t pkt_num Line 1433  UInt_t GL_TIMESYNC::DBpkt(UInt_t pkt_num
1433   */   */
1434  TString GL_TIMESYNC::ConvertTime(TString &tzone, UInt_t dbt){  TString GL_TIMESYNC::ConvertTime(TString &tzone, UInt_t dbt){
1435    //    //
1436    TDatime *time = new TDatime();    Int_t offset = 0;
1437    TString rtime;    TString rtime;
1438    Bool_t found = false;    Bool_t found = false;
1439    //    //
   time->Set(dbt,false); // UTC (Coordinated Universal Time)    
   //  
1440    if ( !strcmp(tzone.Data(),"MSK") || !strcmp(tzone.Data(),"MWT") ){    if ( !strcmp(tzone.Data(),"MSK") || !strcmp(tzone.Data(),"MWT") ){
1441      //      //
1442      UInt_t timeMSK = time->Convert() + 60*60*3; // UTC (Coordinated Universal Time) + 3hs = Moscow Winter Time      offset = 60*60*3; // UTC (Coordinated Universal Time) + 3hs = Moscow Winter Time
     time->Set(timeMSK,false);  
1443      found = true;      found = true;
1444      //      //
1445    };    };
1446    //    //
1447    if ( !strcmp(tzone.Data(),"CET") ){    if ( !strcmp(tzone.Data(),"CET") ){
1448      //      //
1449      UInt_t timeCET = time->Convert() + 60*60*1; // CET (Central European Time) = UTC + 1 hs      offset = 60*60*1; // CET (Central European Time) = UTC + 1 hs
     time->Set(timeCET,false);  
1450      found = true;      found = true;
1451      //      //
1452    };    };
1453    //    //
1454    if ( !strcmp(tzone.Data(),"CEST") ){    if ( !strcmp(tzone.Data(),"CEST") ){
1455      //      //
1456      UInt_t timeCEST = time->Convert() + 60*60*2; // CEST (Central European Summer Time) = UTC + 2 h      offset = 60*60*2; // CEST (Central European Summer Time) = UTC + 2 h
     time->Set(timeCEST,false);  
1457      found = true;      found = true;
1458      //      //
1459    };    };
1460    //    //
1461    if ( !strcmp(tzone.Data(),"MSD") || !strcmp(tzone.Data(),"MST")){    if ( !strcmp(tzone.Data(),"MSD") || !strcmp(tzone.Data(),"MST")){
1462      //      //
1463      UInt_t timeMSD = time->Convert() + 60*60*4; // MSD (Moscow Summer Time) = UTC + 4 h      offset = 60*60*4; // MSD (Moscow Summer Time) = UTC + 4 h
     time->Set(timeMSD,false);  
1464      found = true;      found = true;
1465      //      //
1466    };    };
# Line 1478  TString GL_TIMESYNC::ConvertTime(TString Line 1470  TString GL_TIMESYNC::ConvertTime(TString
1470      tzone = "UTC";      tzone = "UTC";
1471    };    };
1472    //    //
1473    rtime = time->AsSQLString();    dbt += offset;
1474      //
1475      TTimeStamp *time = new TTimeStamp((time_t)dbt,0);
1476      //
1477      rtime = time->AsString("s");
1478      //
1479      delete time;
1480    //    //
1481    return(rtime);    return(rtime);
1482  }  }
1483    
1484  /*  /*
1485   *   *
1486   * Convert the time in the DB from UInt_t to a string   * Convert the time from TZONE to UTC
1487   *   *
1488   * @param dbt time in the DB   * @param dbt time in the DB
1489   * @param tzone Time Zone, can be UTC,GMT,CET,CEST,MSD default is MSK   * @param tzone Time Zone, can be UTC,GMT,CET,CEST,MSD default is MSK
# Line 1493  TString GL_TIMESYNC::ConvertTime(TString Line 1491  TString GL_TIMESYNC::ConvertTime(TString
1491   */   */
1492  TString GL_TIMESYNC::UnConvertTime(TString &tzone, UInt_t dbt){  TString GL_TIMESYNC::UnConvertTime(TString &tzone, UInt_t dbt){
1493    //    //
1494    TDatime *time = new TDatime();    Int_t offset = 0;
1495    TString rtime;    TString rtime;
1496    //    //
1497    Bool_t found = false;    Bool_t found = false;
1498    //    //
   time->Set(dbt,false); // tzone  
   //  
1499    if ( !strcmp(tzone.Data(),"MSK") || !strcmp(tzone.Data(),"MWT") ){    if ( !strcmp(tzone.Data(),"MSK") || !strcmp(tzone.Data(),"MWT") ){
1500      //      //
1501      UInt_t timeUTC = time->Convert() - 60*60*3; // UTC (Coordinated Universal Time) = Moscow Winter Time - 3hs      offset = -60*60*3; // UTC (Coordinated Universal Time) = Moscow Winter Time - 3hs
     time->Set(timeUTC,false);  
1502      found = true;      found = true;
1503      //      //
1504    };    };
1505    //    //
1506    if ( !strcmp(tzone.Data(),"CET") ){    if ( !strcmp(tzone.Data(),"CET") ){
1507      //      //
1508      UInt_t timeCET = time->Convert() - 60*60*1; // CET (Central European Time) - 1 hs = UTC        offset = -60*60*1; // CET (Central European Time) - 1 hs = UTC  
     time->Set(timeCET,false);  
1509      found = true;      found = true;
1510      //      //
1511    };    };
1512    //    //
1513    if ( !strcmp(tzone.Data(),"CEST") ){    if ( !strcmp(tzone.Data(),"CEST") ){
1514      //      //
1515      UInt_t timeCEST = time->Convert() - 60*60*2; // CEST (Central European Summer Time) - 2 h = UTC      offset = -60*60*2; // CEST (Central European Summer Time) - 2 h = UTC
     time->Set(timeCEST,false);  
1516      found = true;      found = true;
1517      //      //
1518    };    };
1519    //    //
1520    if ( !strcmp(tzone.Data(),"MSD") || !strcmp(tzone.Data(),"MST") ){    if ( !strcmp(tzone.Data(),"MSD") || !strcmp(tzone.Data(),"MST") ){
1521      //      //
1522      UInt_t timeMSD = time->Convert() - 60*60*4; // MSD (Moscow Summer Time) - 4 h = UTC      offset = -60*60*4; // MSD (Moscow Summer Time) - 4 h = UTC
     time->Set(timeMSD,false);  
1523      found = true;      found = true;
1524      //      //
1525    };    };
1526    //    //
1527    if ( !found && strcmp(tzone.Data(),"UTC") && strcmp(tzone.Data(),"GMT") && strcmp(tzone.Data(),"") ){    if ( !found && strcmp(tzone.Data(),"UTC") && strcmp(tzone.Data(),"GMT") && strcmp(tzone.Data(),"") ){
1528      //      //
1529        offset = 0;
1530      printf("\n Unknown time zone %s using UTC \n",tzone.Data());      printf("\n Unknown time zone %s using UTC \n",tzone.Data());
1531      tzone = "UTC";      tzone = "UTC";
1532    };    };
1533    //    //
1534    rtime = time->AsSQLString();    dbt += offset;
1535      TTimeStamp *time = new TTimeStamp((time_t)dbt,0);
1536      //
1537      rtime = time->AsString("s");
1538    //    //
1539    return(rtime);    return(rtime);
1540  }  }
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23