/[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.32 by mocchiut, Mon Apr 30 08:19:36 2007 UTC revision 1.41 by mocchiut, Mon Apr 30 10:26:37 2007 UTC
# Line 400  Int_t PamelaDBOperations::SetUpperLimits Line 400  Int_t PamelaDBOperations::SetUpperLimits
400    upperobt = OBT(obtlast);    upperobt = OBT(obtlast);
401    upperentry = nevent-1;    upperentry = nevent-1;
402    //    //
403    if ( IsDebug() ) printf(" First entries are: OBT %i pkt_num %i \n",obtfirst,pktfirst);    if ( IsDebug() ) printf(" First entries are: OBT %u pkt_num %u \n",obtfirst,pktfirst);
404    //    //
405    if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld entry %i\n",upperobt,upperpkt,upperentry);    if ( IsDebug() ) printf(" Last entries are: OBT %lld pkt_num %lld entry %i\n",upperobt,upperpkt,upperentry);
406    //    //
407    if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(1);    if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) || (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) ) return(32);
408    //    //
409    if ( !nevent ) return(2);    if ( !nevent ) return(64);
410    //    //
411    if ( nevent < 2 ) return(4);    if ( nevent < 2 ) return(128);
412    if ( nevent < jump ) jump = 1;    if ( nevent < jump ) jump = 1;
413    //  if ( nevent < jump ) jump = int(nevent/10);    //  if ( nevent < jump ) jump = int(nevent/10);
414    //  if ( !jump ) jump = 1;    //  if ( !jump ) jump = 1;
# Line 447  Int_t PamelaDBOperations::SetUpperLimits Line 447  Int_t PamelaDBOperations::SetUpperLimits
447          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){          if ( t_pktlast < upperpkt && t_obtlast < upperobt && t_pktlast < upperpkt2 && t_obtlast < upperobt2 ){
448            zomp = i + jump + 1;            zomp = i + jump + 1;
449            if ( zomp > nevent-2 ) zomp = nevent - 2;            if ( zomp > nevent-2 ) zomp = nevent - 2;
450            if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %i upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);            if ( IsDebug() ) printf(" .-. jump %i zomp %i upperpkt %lld pktlast %u upperobt %lld obtlast %u last entry is %i \n",jump,zomp,upperpkt,pktlast,upperobt,obtlast,i);
451            break;            break;
452          };          };
453          //          //
# Line 867  Long64_t PamelaDBOperations::PKT(UInt_t Line 867  Long64_t PamelaDBOperations::PKT(UInt_t
867      return((Long64_t)pkt_num+16777215LL);      return((Long64_t)pkt_num+16777215LL);
868    };    };
869    //    //
870    if ( pkt_num > (Long64_t)ppktfirst*2 && pkt_num > (16777214/2) ){    if ( pkt_num > ((Long64_t)ppktfirst*2) && pkt_num > (16777214/2) ){
871      if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL);      if ( IsDebug() ) printf(" rise down pktnum %lld \n",(Long64_t)pkt_num-16777215LL);
872      return((Long64_t)pkt_num-16777215LL);      return((Long64_t)pkt_num-16777215LL);
873    };    };
# Line 882  Long64_t PamelaDBOperations::PKT(UInt_t Line 882  Long64_t PamelaDBOperations::PKT(UInt_t
882   */   */
883  Long64_t PamelaDBOperations::OBT(UInt_t obt){    Long64_t PamelaDBOperations::OBT(UInt_t obt){  
884    //    //
885    if ( IsDebug() ) printf(" obt conversion: ob is %u obtfirst is %u  (numeric_limits<UInt_t>::max()/2) is %u  \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2));    if ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u  (numeric_limits<UInt_t>::max()/2) is %u  \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2));
886    //    //
887    if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){    if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){
888      if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)(obt+numeric_limits<UInt_t>::max()));      if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)(obt+numeric_limits<UInt_t>::max()));
889      return((Long64_t)(obt+numeric_limits<UInt_t>::max()));      return((Long64_t)(obt+numeric_limits<UInt_t>::max()));
890    };    };
891    //    //
892    if ( obt > (Long64_t)(pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){    if ( obt > ((Long64_t)pobtfirst*2) && obt > (numeric_limits<UInt_t>::max()/2) ){
893        if ( IsDebug() ) printf(" pobtfirst*2 %lld \n",((Long64_t)pobtfirst*2));
894      if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());      if ( IsDebug() ) printf(" rise down pktnum %lld \n", (Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
895      return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());      return((Long64_t)obt-(Long64_t)numeric_limits<UInt_t>::max());
896    };    };
# Line 1034  Int_t PamelaDBOperations::insertPamelaGL Line 1035  Int_t PamelaDBOperations::insertPamelaGL
1035    TSQLResult *result = 0;    TSQLResult *result = 0;
1036    TSQLRow *row = 0;    TSQLRow *row = 0;
1037    UInt_t t0 = 0;    UInt_t t0 = 0;
1038      Int_t signal = 0;
1039      //
1040      signal = this->SetUpperLimits();
1041    //    //
1042    stringstream   oss;    stringstream   oss;
1043    //      //  
# Line 1086  Int_t PamelaDBOperations::insertPamelaGL Line 1090  Int_t PamelaDBOperations::insertPamelaGL
1090    };    };
1091    //    //
1092    TTree *T = 0;    TTree *T = 0;
   Int_t signal = 0;  
1093    //    //
1094    UInt_t nevent = 0;    UInt_t nevent = 0;
1095    UInt_t recEntries = 0;    UInt_t recEntries = 0;
# Line 1495  Int_t PamelaDBOperations::insertPamelaRU Line 1498  Int_t PamelaDBOperations::insertPamelaRU
1498    stringstream oss;    stringstream oss;
1499    oss.str("");    oss.str("");
1500    //    //
1501    signal = this->SetUpperLimits();    //  signal = this->SetUpperLimits();
1502    //    //
1503    // loop on runheader and runtrailer events    // loop on runheader and runtrailer events
1504    //    //
# Line 1546  Int_t PamelaDBOperations::insertPamelaRU Line 1549  Int_t PamelaDBOperations::insertPamelaRU
1549        //        //
1550        if ( !ptt &&  !(ptht+1) ){        if ( !ptt &&  !(ptht+1) ){
1551          //          //
1552          if ( IsDebug() ) printf(" Piece of run at the beginning of the file %i %i %i \n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Piece of run at the beginning of the file %u %u %u \n",ptht,pth,ptt);
1553          //          //
1554          this->HandleRunFragments(true,false,0,(evbeft-1));          this->HandleRunFragments(true,false,0,(evbeft-1));
1555          //          //
1556          //          //
1557        } else if ( pth == ptht ){        } else if ( pth == ptht ){
1558          //          //
1559          if ( IsDebug() ) printf(" Missing header %i %i %i\n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Missing header %u %u %u\n",ptht,pth,ptt);
1560          //          //
1561          if ( (ptt-1) < 0 ) throw -15; // should never arrive here!          if ( (ptt-1) < 0 ) throw -15; // should never arrive here!
1562          rt->GetEntry(ptt-1);          rt->GetEntry(ptt-1);
# Line 1562  Int_t PamelaDBOperations::insertPamelaRU Line 1565  Int_t PamelaDBOperations::insertPamelaRU
1565          rt->GetEntry(ptt);          rt->GetEntry(ptt);
1566          pht = eht->GetPscuHeader();          pht = eht->GetPscuHeader();
1567          //          //
1568          if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %i %i %i \n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Try to find the beginning of a run which has only the runtrailer   %u %u %u \n",ptht,pth,ptt);
1569          if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %i %i %i \n",pkth,obth,obtt);          if ( IsDebug() ) printf(" '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''  %u %u %u \n",pkth,obth,obtt);
1570          //          //
1571          this->HandleMissingHoT(true,false,evbefh,evbeft-1);          this->HandleMissingHoT(true,false,evbefh,evbeft-1);
1572          //          //
1573        } else {        } else {
1574          //          //
1575          if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %i %i %i\n",ptht,pth,ptt);          if ( IsDebug() ) printf(" Could be a good run, we have a runheader followed by a runtrailer %u %u %u\n",ptht,pth,ptt);
1576          //          //
1577          rh->GetEntry(ptht);          rh->GetEntry(ptht);
1578          phh = ehh->GetPscuHeader();          phh = ehh->GetPscuHeader();
# Line 1577  Int_t PamelaDBOperations::insertPamelaRU Line 1580  Int_t PamelaDBOperations::insertPamelaRU
1580          obth = phh->GetOrbitalTime();          obth = phh->GetOrbitalTime();
1581          cod = ehh->GetCounter();          cod = ehh->GetCounter();
1582          evbefh = cod->Get(pctp->Physics);          evbefh = cod->Get(pctp->Physics);
1583          if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);          if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
1584          //          //
1585          // handle this run          // handle this run
1586          //          //
# Line 1587  Int_t PamelaDBOperations::insertPamelaRU Line 1590  Int_t PamelaDBOperations::insertPamelaRU
1590          //          //
1591          if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){          if ( PKT(pkth)>PKT(pktfirst) && OBT(obth)>OBT(obtfirst) && !ptt ){
1592            //            //
1593            if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER \n");            if ( IsDebug() ) printf(" Piece of run at the beginning of the file WITH NO RUNTRAILER evbefh = %u \n",evbefh);        
           //  
           this->HandleRunFragments(true,true,0,(evbefh-1));  
1594            //            //
1595              if ( evbefh == 0 ) {
1596                //
1597                signal = 8;
1598                if ( IsDebug() ) printf(" Not supported yet: run with no events, no runtrailer, no runheader \n");
1599                //
1600              } else {
1601                //
1602                this->HandleRunFragments(true,true,0,(evbefh-1));
1603                //
1604              };
1605          };          };
1606          //          //
1607          //          //
1608          if ( (ptht - pth) > 1 ){          if ( (ptht - pth) > 1 ){
1609            //            //
1610            if ( IsDebug() ) printf(" Missing runtrailers! \n");            if ( IsDebug() ) printf(" Missing runtrailers! \n");
1611            if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" Attention there is a jump in the runheader counter %u %u %u \n",ptht,pth,ptt);
1612            // is not the consecutive header            // is not the consecutive header
1613            while ( pth != ptht ){                      while ( pth != ptht ){          
1614              //              //
# Line 1618  Int_t PamelaDBOperations::insertPamelaRU Line 1629  Int_t PamelaDBOperations::insertPamelaRU
1629              obth = phh->GetOrbitalTime();              obth = phh->GetOrbitalTime();
1630              evbefh = cod->Get(pctp->Physics);              evbefh = cod->Get(pctp->Physics);
1631              //              //
1632              if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %i %i %i \n",ptht,pth,ptt);              if ( IsDebug() ) printf(" Try to find the end of a run which has only the runheader   %u %u %u \n",ptht,pth,ptt);
1633              if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %i %i %i \n",pkth,obth,obtt);              if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' %u %u %u \n",pkth,obth,obtt);
1634              //              //
1635              this->HandleMissingHoT(false,true,evbefh,evbeft-1);              this->HandleMissingHoT(false,true,evbefh,evbeft-1);
1636              //              //
# Line 1628  Int_t PamelaDBOperations::insertPamelaRU Line 1639  Int_t PamelaDBOperations::insertPamelaRU
1639          } else if ( !(ptht - pth) ){          } else if ( !(ptht - pth) ){
1640            //            //
1641            if ( IsDebug() ) printf(" Missing runheader! \n");            if ( IsDebug() ) printf(" Missing runheader! \n");
1642            if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" Attention! the runheader counter did not changed %u %u %u \n",ptht,pth,ptt);
1643            if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n");            if ( IsDebug() ) printf(" The run should have already been handled by HandleRun() \n");
1644            //            //
1645          } else {          } else {
# Line 1649  Int_t PamelaDBOperations::insertPamelaRU Line 1660  Int_t PamelaDBOperations::insertPamelaRU
1660            obth = phh->GetOrbitalTime();            obth = phh->GetOrbitalTime();
1661            cod = ehh->GetCounter();            cod = ehh->GetCounter();
1662            evbefh = cod->Get(pctp->Physics);            evbefh = cod->Get(pctp->Physics);
1663            if ( IsDebug() ) printf(" Piece of run at the end of file %i %i %i \n",pkth,obth,obtt);            if ( IsDebug() ) printf(" Piece of run at the end of file %u %u %u \n",pkth,obth,obtt);
1664            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i %i %i \n",ptht,pth,ptt);            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u %u %u \n",ptht,pth,ptt);
1665            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %i \n",rhev);            if ( IsDebug() ) printf(" ''''''''''''''''''''''''''''''' %u \n",rhev);
1666            //            //
1667            this->HandleRunFragments(false,true,evbefh,upperentry);            this->HandleRunFragments(false,true,evbefh,upperentry);
1668          } else {          } else {
# Line 2051  void PamelaDBOperations::HandleRunFragme Line 2062  void PamelaDBOperations::HandleRunFragme
2062        ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT());        ULong64_t aobt = OBT(glrun1->GetRUNTRAILER_OBT());
2063        UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT());        UInt_t bpkt = PKT(glrun->GetRUNHEADER_PKT());
2064        ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT());        ULong64_t bobt = OBT(glrun->GetRUNHEADER_OBT());
2065        if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);        if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
2066        TTree *T= 0;        TTree *T= 0;
2067        T = (TTree*)file->Get("Physics");        T = (TTree*)file->Get("Physics");
2068        if ( !T || T->IsZombie() ) throw -16;        if ( !T || T->IsZombie() ) throw -16;
# Line 2065  void PamelaDBOperations::HandleRunFragme Line 2076  void PamelaDBOperations::HandleRunFragme
2076          bobt = OBT(ph->GetOrbitalTime());                bobt = OBT(ph->GetOrbitalTime());      
2077          firstev++;          firstev++;
2078        };        };
2079        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);        if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu firstev is %i\n",apkt,bpkt,aobt,bobt,firstev);
2080        //        //
2081        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());        glrun1->SetPKT_COUNTER(glrun->GetPKT_COUNTER());
2082        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());        glrun1->SetPKT_READY_COUNTER(glrun->GetPKT_READY_COUNTER());
# Line 2248  void PamelaDBOperations::HandleRunFragme Line 2259  void PamelaDBOperations::HandleRunFragme
2259        ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT());        ULong64_t aobt = OBT(glrun->GetRUNTRAILER_OBT());
2260        UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT());        UInt_t bpkt = PKT(glrun1->GetRUNHEADER_PKT());
2261        ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT());        ULong64_t bobt = OBT(glrun1->GetRUNHEADER_OBT());
2262        if ( IsDebug() ) printf(" Check overlapping events: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);        if ( IsDebug() ) printf(" Check overlapping events: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
2263        TTree *T= 0;        TTree *T= 0;
2264        T = (TTree*)file->Get("Physics");        T = (TTree*)file->Get("Physics");
2265        if ( !T || T->IsZombie() ) throw -16;        if ( !T || T->IsZombie() ) throw -16;
# Line 2262  void PamelaDBOperations::HandleRunFragme Line 2273  void PamelaDBOperations::HandleRunFragme
2273          aobt = OBT(ph->GetOrbitalTime());                aobt = OBT(ph->GetOrbitalTime());      
2274          lastev--;          lastev--;
2275        };        };
2276        if ( IsDebug() ) printf(" Check overlapping events done: %i %i %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);        if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev);
2277        //        //
2278        glrun->SetEV_TO(lastev);        glrun->SetEV_TO(lastev);
2279        glrun->SetNEVENTS(lastev-firstev+1);        glrun->SetNEVENTS(lastev-firstev+1);
# Line 2869  Int_t PamelaDBOperations::insertCALO_CAL Line 2880  Int_t PamelaDBOperations::insertCALO_CAL
2880          fromtime = this->GetAbsTime(ph->GetOrbitalTime());            fromtime = this->GetAbsTime(ph->GetOrbitalTime());  
2881          if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){          if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
2882            //            //
2883            if ( IsDebug() ) printf(" Calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);            if ( IsDebug() ) printf(" Calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
2884            //            //
2885            // check if the calibration has already been inserted            // check if the calibration has already been inserted
2886            //            //
# Line 2914  Int_t PamelaDBOperations::insertCALO_CAL Line 2925  Int_t PamelaDBOperations::insertCALO_CAL
2925                //                //
2926                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);                if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB for section %i \n",section);
2927                if ( fromtime < 1150871000 ){ //1150866904                if ( fromtime < 1150871000 ){ //1150866904
2928                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);                  if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
2929                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode                  fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
2930                };                };
2931                //                //
# Line 2976  Int_t PamelaDBOperations::insertCALO_CAL Line 2987  Int_t PamelaDBOperations::insertCALO_CAL
2987            //            //
2988          } else {          } else {
2989            //            //
2990            if ( IsDebug() ) printf(" Repeated calo calibration for section %i at time %i obt %i pkt %i \n",section,fromtime,obt,pkt);            if ( IsDebug() ) printf(" Repeated calo calibration for section %i at time %u obt %u pkt %u \n",section,fromtime,obt,pkt);
2991            //            //
2992          };                };      
2993          //                //      
# Line 3194  Int_t PamelaDBOperations::insertTRK_CALI Line 3205  Int_t PamelaDBOperations::insertTRK_CALI
3205      //      //
3206      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){      if ( this->PKT(pkt1) >= this->PKT(pktfirst) && this->OBT(obt1) >= this->OBT(obtfirst) ){
3207        //        //
3208        if ( IsDebug() ) printf(" Trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);        if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
3209        //              //      
3210        valid = ValidateTrkCalib( caltrk1, eh1 );        valid = ValidateTrkCalib( caltrk1, eh1 );
3211        if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;        if ( IsDebug() ) cout << " pkt1 validation --> "<<valid<<endl;
# Line 3237  Int_t PamelaDBOperations::insertTRK_CALI Line 3248  Int_t PamelaDBOperations::insertTRK_CALI
3248          //          //
3249        };        };
3250        //        //
3251        if ( IsDebug() ) printf(" Found trk calibration2 at obt %i pkt %i t2 is %i \n",obt2,pkt2,t2);        if ( IsDebug() ) printf(" Found trk calibration2 at obt %u pkt %u t2 is %u \n",obt2,pkt2,t2);
3252        //        //
3253        // The calibration is good        // The calibration is good
3254        //        //
3255        if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){        if ( this->PKT(pkt2) == this->PKT(pkt1)+1 ){
3256          //          //
3257          if ( IsDebug() ) printf(" The trk calibration2 at obt %i pkt %i t2 is %i is good \n",obt2,pkt2,t2);          if ( IsDebug() ) printf(" The trk calibration2 at obt %u pkt %u t2 is %u is good \n",obt2,pkt2,t2);
3258          //          //
3259          UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );          UInt_t valid2 = ValidateTrkCalib( caltrk2, eh2 );
3260          if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;          if ( IsDebug() ) cout << " pkt2 validation --> "<<valid2<<endl;
# Line 3257  Int_t PamelaDBOperations::insertTRK_CALI Line 3268  Int_t PamelaDBOperations::insertTRK_CALI
3268          //          //
3269          if ( t2 != pret2+1 ){          if ( t2 != pret2+1 ){
3270            //            //
3271            if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i pret2 is %i \n",obt2,pkt2,t2,pret2);            if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u pret2 is %u \n",obt2,pkt2,t2,pret2);
3272            //            //
3273            while ( t2 > pret2+1 ){            while ( t2 > pret2+1 ){
3274              //              //
# Line 3286  Int_t PamelaDBOperations::insertTRK_CALI Line 3297  Int_t PamelaDBOperations::insertTRK_CALI
3297          //          //
3298          // Check for missing calibtrk2          // Check for missing calibtrk2
3299          //          //
3300          if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);          if ( IsDebug() ) printf(" Missing the trk calibration2! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
3301          t2 = pret2;          t2 = pret2;
3302          //          //
3303          // handle missing calib2          // handle missing calib2
# Line 3300  Int_t PamelaDBOperations::insertTRK_CALI Line 3311  Int_t PamelaDBOperations::insertTRK_CALI
3311        //        //
3312      } else {      } else {
3313        //        //
3314        if ( IsDebug() ) printf(" Repetead trk calibration1 at time %i obt %i pkt %i \n",fromtime,obt1,pkt1);        if ( IsDebug() ) printf(" Repetead trk calibration1 at time %u obt %u pkt %u \n",fromtime,obt1,pkt1);
3315        //        //
3316      };      };
3317      //          //    
# Line 3325  Int_t PamelaDBOperations::insertTRK_CALI Line 3336  Int_t PamelaDBOperations::insertTRK_CALI
3336      valid = 0;      valid = 0;
3337      if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){      if ( this->PKT(pkt2) > this->PKT(pktfirst) || this->OBT(obt2) > this->OBT(obtfirst) ){
3338        //        //
3339        if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %i pkt %i t2 is %i\n",obt2,pkt2,t2);        if ( IsDebug() ) printf(" Missing the trk calibration1! Next one at obt %u pkt %u t2 is %u\n",obt2,pkt2,t2);
3340        //        //
3341        this->HandleTRK_CALIB(false,true);        this->HandleTRK_CALIB(false,true);
3342        //        //
# Line 3379  Int_t PamelaDBOperations::insertS4_CALIB Line 3390  Int_t PamelaDBOperations::insertS4_CALIB
3390      fromtime = this->GetAbsTime(ph->GetOrbitalTime());        fromtime = this->GetAbsTime(ph->GetOrbitalTime());  
3391      if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){      if ( this->PKT(pkt) >= this->PKT(pktfirst) && this->OBT(obt) >= this->OBT(obtfirst) ){
3392        //        //
3393        if ( IsDebug() ) printf(" S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);        if ( IsDebug() ) printf(" S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
3394        //        //
3395        // check if the calibration has already been inserted        // check if the calibration has already been inserted
3396        //        //
# Line 3422  Int_t PamelaDBOperations::insertS4_CALIB Line 3433  Int_t PamelaDBOperations::insertS4_CALIB
3433            //            //
3434            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");            if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB \n");
3435            if ( fromtime < 1150871000 ){              if ( fromtime < 1150871000 ){  
3436              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %i \n",fromtime);              if ( IsDebug() ) printf(" First PAMELA flight calibration at time %u \n",fromtime);
3437              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode              fromtime = 0;// the first flight calibration was taken at about  1156429100 s, this line allow to analyze first runs in raw mode
3438            };            };
3439            //            //
# Line 3481  Int_t PamelaDBOperations::insertS4_CALIB Line 3492  Int_t PamelaDBOperations::insertS4_CALIB
3492        //        //
3493      } else {      } else {
3494        //        //
3495        if ( IsDebug() ) printf(" Repetead S4 calibration at time %i obt %i pkt %i \n",fromtime,obt,pkt);        if ( IsDebug() ) printf(" Repeated S4 calibration at time %u obt %u pkt %u \n",fromtime,obt,pkt);
3496        //        //
3497      };            };      
3498      //            //      
# Line 3692  Int_t PamelaDBOperations::ValidateRuns(T Line 3703  Int_t PamelaDBOperations::ValidateRuns(T
3703        row = result->Next();        row = result->Next();
3704        t_stop = (UInt_t)atoll(row->GetField(4));        t_stop = (UInt_t)atoll(row->GetField(4));
3705      };      };
3706      if ( IsDebug() ) printf("Validation interval: from time %i - to time %i \n\n",t_stop,t_start);      if ( IsDebug() ) printf("Validation interval: from time %u - to time %u \n\n",t_stop,t_start);
3707      // --------------------------------------------------------------      // --------------------------------------------------------------
3708      // now retrieves runs to be validated      // now retrieves runs to be validated
3709      // --------------------------------------------------------------      // --------------------------------------------------------------
# Line 3863  Int_t PamelaDBOperations::ValidateRuns(T Line 3874  Int_t PamelaDBOperations::ValidateRuns(T
3874                    }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);
3875            };            };
3876                        
3877            if ( IsDebug() ) printf("%i Run %i \n",nrow,this_run->ID);            if ( IsDebug() ) printf("%i Run %u \n",nrow,this_run->ID);
3878            nrow++;            nrow++;
3879                        
3880    };    };
# Line 3892  Bool_t PamelaDBOperations::MissingTRK_CA Line 3903  Bool_t PamelaDBOperations::MissingTRK_CA
3903          // which should be equal to the time between ascending-nodes.          // which should be equal to the time between ascending-nodes.
3904          //==============================================================          //==============================================================
3905          if ( t2 - trkcalib->FROM_TIME > 5700) {          if ( t2 - trkcalib->FROM_TIME > 5700) {
3906                  if ( IsDebug() )printf("Long time between calib and run start %i :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);                  if ( IsDebug() )printf("Long time between calib and run start %u :-( ==> there might be a missing calib \n",t2 - trkcalib->FROM_TIME);
3907          //==============================================================          //==============================================================
3908          // there might be a missing calibration, due to:          // there might be a missing calibration, due to:
3909          // - MM full          // - MM full
# Line 3913  Bool_t PamelaDBOperations::MissingTRK_CA Line 3924  Bool_t PamelaDBOperations::MissingTRK_CA
3924          // it is enough to say that there are no missing calibrations          // it is enough to say that there are no missing calibrations
3925          //==============================================================          //==============================================================
3926          // the long time interval bewteen runs might be due to download          // the long time interval bewteen runs might be due to download
3927          if ( IsDebug() )printf("Short time between calib and run start %i :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);          if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME);
3928          return(false);          return(false);
3929                    
3930  };  };
# Line 4229  void PamelaDBOperations::ValidationOFF() Line 4240  void PamelaDBOperations::ValidationOFF()
4240        };        };
4241      };      };
4242    };    };
4243    if ( IsDebug() ) printf(" %i runs have been unvalidated \n",unv);    if ( IsDebug() ) printf(" %u runs have been unvalidated \n",unv);
4244  };  };
4245    
4246  /**  /**
# Line 4260  void PamelaDBOperations::RemoveRUNS(){ Line 4271  void PamelaDBOperations::RemoveRUNS(){
4271      if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");          if ( IsDebug() ) printf(" Deleting run from GL_RUN table \n");    
4272      while ( Row ){      while ( Row ){
4273        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4274        if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));        if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4275        drun++;        drun++;
4276        Row = pResult->Next();            Row = pResult->Next();    
4277      };      };
# Line 4281  void PamelaDBOperations::RemoveRUNS(){ Line 4292  void PamelaDBOperations::RemoveRUNS(){
4292    } else {    } else {
4293      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");          if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN table \n");    
4294      while ( Row ){      while ( Row ){
4295        if ( IsDebug() ) printf(" restore run %i \n",(UInt_t)atoll(Row->GetField(1)));        if ( IsDebug() ) printf(" restore run %u \n",(UInt_t)atoll(Row->GetField(1)));
4296        delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");        delrun->RestoreRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN_FRAGMENTS");
4297        if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(1)));        if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(1)));
4298        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");        delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(1)),"GL_RUN");
4299        if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){        if ( (UInt_t)atoll(Row->GetField(1)) != (UInt_t)atoll(Row->GetField(0)) ){
4300          if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));          if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4301          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");          delrun->DeleteRun(conn,(UInt_t)atoll(Row->GetField(0)),"GL_RUN");
4302        };        };
4303        drun++;        drun++;
# Line 4312  void PamelaDBOperations::RemoveRUNS(){ Line 4323  void PamelaDBOperations::RemoveRUNS(){
4323    } else {    } else {
4324      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");          if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_TRASH table \n");    
4325      while ( Row ){      while ( Row ){
4326        if ( IsDebug() ) printf(" del run idtrash %i \n",(UInt_t)atoll(Row->GetField(0)));        if ( IsDebug() ) printf(" del run idtrash %u \n",(UInt_t)atoll(Row->GetField(0)));
4327        myquery.str("");        myquery.str("");
4328        myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";              myquery << " DELETE FROM GL_RUN_TRASH where ID_TRASH=" << Row->GetField(0) <<";";      
4329        conn->Query(myquery.str().c_str());        conn->Query(myquery.str().c_str());
# Line 4321  void PamelaDBOperations::RemoveRUNS(){ Line 4332  void PamelaDBOperations::RemoveRUNS(){
4332      };      };
4333    };    };
4334    //      //  
4335    if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_TRASH table \n",drun);        if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_TRASH table \n",drun);    
4336    //    //
4337    //    //
4338    //    //
# Line 4339  void PamelaDBOperations::RemoveRUNS(){ Line 4350  void PamelaDBOperations::RemoveRUNS(){
4350    } else {    } else {
4351      if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");          if ( IsDebug() ) printf(" Deleting run fragments from GL_RUN_FRAGMENTS table \n");    
4352      while ( Row ){      while ( Row ){
4353        if ( IsDebug() ) printf(" del run %i \n",(UInt_t)atoll(Row->GetField(0)));        if ( IsDebug() ) printf(" del run %u \n",(UInt_t)atoll(Row->GetField(0)));
4354        myquery.str("");        myquery.str("");
4355        myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";        myquery << " DELETE FROM GL_RUN_FRAGMENTS where ID=" << Row->GetField(0) <<";";
4356        conn->Query(myquery.str().c_str());        conn->Query(myquery.str().c_str());
# Line 4348  void PamelaDBOperations::RemoveRUNS(){ Line 4359  void PamelaDBOperations::RemoveRUNS(){
4359      };      };
4360    };    };
4361    //      //  
4362    if ( IsDebug() ) printf(" Deleted %i run(s) from GL_RUN_FRAGMENTS table \n",drun);        if ( IsDebug() ) printf(" Deleted %u run(s) from GL_RUN_FRAGMENTS table \n",drun);    
4363    //    //
4364    //    //
4365    //    //
# Line 4679  UInt_t PamelaDBOperations::Check(){ Line 4690  UInt_t PamelaDBOperations::Check(){
4690      //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){      //    if ( (thisrht < prevrtt || thisrtt < prevrht || thisrht > thisrtt) && (thisrht != prevrht) ){
4691      if ( (thisrht < prevrtt) && (thisrht != prevrht) ){      if ( (thisrht < prevrtt) && (thisrht != prevrht) ){
4692        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4693        printf(" CHECK n.1 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid);        printf(" CHECK n.1 TIME SCREW of %u s AROUND RUNs %u and %u \n",(thisrht-prevrtt),previd,thisid);
4694        TString prevf = "";        TString prevf = "";
4695        TString thisf = "";        TString thisf = "";
4696        oss.str("");        oss.str("");
# Line 4700  UInt_t PamelaDBOperations::Check(){ Line 4711  UInt_t PamelaDBOperations::Check(){
4711      //      //
4712      if ( (thisrtt < prevrht) && (thisrht != prevrht) ){      if ( (thisrtt < prevrht) && (thisrht != prevrht) ){
4713        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4714        printf(" CHECK n.2 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid);        printf(" CHECK n.2 TIME SCREW of %u s AROUND RUNs %u and %u \n",(thisrtt-prevrht),previd,thisid);
4715        TString prevf = "";        TString prevf = "";
4716        TString thisf = "";        TString thisf = "";
4717        oss.str("");        oss.str("");
# Line 4721  UInt_t PamelaDBOperations::Check(){ Line 4732  UInt_t PamelaDBOperations::Check(){
4732      //      //
4733      if ( (thisrht > thisrtt) && (thisrht != prevrht) ){      if ( (thisrht > thisrtt) && (thisrht != prevrht) ){
4734        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);        if ( IsDebug() ) printf(" IDprev %u ID %u    prevrht %u prevrtt %u  thisrht %u thisrtt %u \n",previd,thisid,prevrht,prevrtt,thisrht,thisrtt);
4735        printf(" CHECK n.3 TIME SCREW of %i s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid);        printf(" CHECK n.3 TIME SCREW of %u s AROUND RUNs %u and %u \n",(thisrht-thisrtt),previd,thisid);
4736        TString prevf = "";        TString prevf = "";
4737        TString thisf = "";        TString thisf = "";
4738        oss.str("");        oss.str("");
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.23