/[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.23 by mocchiut, Tue Jan 13 13:38:15 2009 UTC revision 1.29 by mocchiut, Wed Dec 23 07:19:32 2009 UTC
# Line 183  void PamelaDBOperations::OpenL0File(TStr Line 183  void PamelaDBOperations::OpenL0File(TStr
183            chboot = (UInt_t)atoll(row2->GetField(5));            chboot = (UInt_t)atoll(row2->GetField(5));
184            //            //
185            if ( !chboot ){            if ( !chboot ){
186                if ( debug ) printf(" mmm... no BOOT number determined by chewbacca, try to find it looking at DB \n");
187                //
188                // new feature: check in ROOT_TABLE (timewise) the previous and the next BOOT number != 0 , if it is the same we know this one!
189                // if it is different roll back to the old behaviour
190                //
191              qu.str("");              qu.str("");
192              //      qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<10000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<7200 group by BOOT_NUMBER order by BOOT_NUMBER asc;";              qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND REAL_TIME_INIT>"<< row2->GetField(10) <<" order by REAL_TIME_INIT asc limit 1;";
             qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<40000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<40000 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; // relax conditions... 090112 [8RED: error -29]  
193              if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());              if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
194                UInt_t chbootA = 0;
195                UInt_t chbootB = 1;
196              result3 = conn->Query(qu.str().c_str());              result3 = conn->Query(qu.str().c_str());
197              if ( result3 ){              if ( result3 ){
198                row3 = result3->Next();                row3 = result3->Next();
199                if ( row3 && result3->GetRowCount() == 1 ){                if ( row3 ){
200                  chboot = (UInt_t)atoll(row3->GetField(0));                  chbootA = (UInt_t)atoll(row3->GetField(0));
201                  if ( debug ) printf(" Found boot_number = %u \n",chboot);                  if ( debug ) printf(" Found boot_number A = %u \n",chbootA);
202                } else {                };
203                  if ( debug ) printf(" AGH CANNOT DETERMINE THE BOOT NUMBER... \n");              };
204                  throw -29;              delete result3;
205                qu.str("");
206                qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND REAL_TIME_INIT<"<< row2->GetField(10) <<" order by REAL_TIME_INIT desc limit 1;";
207                if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
208                result3 = conn->Query(qu.str().c_str());
209                if ( result3 ){
210                  row3 = result3->Next();
211                  if ( row3 ){
212                    chbootB = (UInt_t)atoll(row3->GetField(0));
213                    if ( debug ) printf(" Found boot_number B = %u \n",chbootB);
214                  };
215                };
216                if ( chbootA == chbootB ){
217                  chboot = chbootA;
218                  if ( debug ) printf(" Found boot_number! it is %u \n",chboot);
219                } else {        
220                  qu.str("");
221                  //            qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<10000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<7200 group by BOOT_NUMBER order by BOOT_NUMBER asc;";
222                  //            qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<40000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<40000 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; // relax conditions... 090112 [8RED: error -29]
223                  qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<50000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<50000 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; // even more relaxed conditions... 091214 [9REDtest3: error -29]
224                  if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str());
225                  result3 = conn->Query(qu.str().c_str());
226                  if ( result3 ){
227                    row3 = result3->Next();
228                    if ( row3 && result3->GetRowCount() == 1 ){
229                      chboot = (UInt_t)atoll(row3->GetField(0));
230                      if ( debug ) printf(" Found boot_number = %u \n",chboot);
231                    } else {
232                      if ( debug ) printf(" AGH CANNOT DETERMINE THE BOOT NUMBER... \n");
233                      throw -29;
234                    };
235                };                };
236              };              };
237            };            };
# Line 427  void PamelaDBOperations::SetOrbitNo(UInt Line 463  void PamelaDBOperations::SetOrbitNo(UInt
463        if ( IsDebug() ) printf(" Agh problems determining the orbit number! name = %s \n",name.Data());        if ( IsDebug() ) printf(" Agh problems determining the orbit number! name = %s \n",name.Data());
464        return;        return;
465      };      };
466      TString dwo = 0;      TString dwo = "";
467      for (Int_t i = 0; i<5; i++){      for (Int_t i = 0; i<5; i++){
468        dwo.Append(name[i],1);        dwo.Append(name[i],1);
469      };      };
# Line 857  Int_t PamelaDBOperations::SetUpperLimits Line 893  Int_t PamelaDBOperations::SetUpperLimits
893    //  if ( nevent < jump ) jump = int(nevent/10);    //  if ( nevent < jump ) jump = int(nevent/10);
894    //  if ( !jump ) jump = 1;    //  if ( !jump ) jump = 1;
895    //    //
896    if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) && nevent > deltapkt ){  
897      if ( ((PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt)) && nevent > deltapkt ){
898      //      //
899      if ( IsDebug() ) printf(" starting jump %i \n",jump);      if ( IsDebug() ) printf(" starting jump %i \n",jump);
900      //    if ( IsDebug() ) printf(" (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) %llu < %llu && %llu > %llu \n     OR \n     (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) %llu > %llu && %llu < %llu \n",PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst),PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst));      //    if ( IsDebug() ) printf(" (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) %llu < %llu && %llu > %llu \n     OR \n     (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) %llu > %llu && %llu < %llu \n",PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst),PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst));
# Line 4235  Int_t PamelaDBOperations::insertCALO_CAL Line 4272  Int_t PamelaDBOperations::insertCALO_CAL
4272    UInt_t totime = 0;    UInt_t totime = 0;
4273    UInt_t obt = 0;    UInt_t obt = 0;
4274    UInt_t pkt = 0;    UInt_t pkt = 0;
4275      Float_t totped = 0.;
4276    //    //
4277    tr = (TTree*)file->Get("CalibCalPed");    tr = (TTree*)file->Get("CalibCalPed");
4278    if ( !tr || tr->IsZombie() ) throw -21;    if ( !tr || tr->IsZombie() ) throw -21;
# Line 4252  Int_t PamelaDBOperations::insertCALO_CAL Line 4290  Int_t PamelaDBOperations::insertCALO_CAL
4290        if ( calibCalPed->cstwerr[section] ){        if ( calibCalPed->cstwerr[section] ){
4291          valid = 1;          valid = 1;
4292          if ( calibCalPed->cperror[section] ) valid = 0;          if ( calibCalPed->cperror[section] ) valid = 0;
4293            //
4294            // check pedestal values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.]
4295            //
4296            totped = 0.;
4297            Int_t ns = 0;
4298            if ( section == 2 ) ns = 3;
4299            if ( section == 3 ) ns = 1;
4300            if ( section == 1 ) ns = 2;
4301            for (UInt_t ss=0; ss<96; ss++){
4302              totped += fabs(calibCalPed->calped[ns][0][ss]);
4303            }
4304            if ( totped < 1. ){
4305              if ( IsDebug() ) printf(" Section %i totped %f -  No calibration data! Calorimeter power problems? \n",section,totped);
4306              valid = 0;
4307            };
4308            //
4309          ph = eh->GetPscuHeader();          ph = eh->GetPscuHeader();
4310          obt = ph->GetOrbitalTime();            obt = ph->GetOrbitalTime();  
4311          pkt = ph->GetCounter();            pkt = ph->GetCounter();  
# Line 4419  Int_t PamelaDBOperations::insertCALOPULS Line 4473  Int_t PamelaDBOperations::insertCALOPULS
4473    tr->SetBranchAddress("Header", &eh);    tr->SetBranchAddress("Header", &eh);
4474    nevents = tr->GetEntries();    nevents = tr->GetEntries();
4475    //    //
4476      Float_t totpul = 0.;
4477      //
4478    if ( nevents > 0 ){    if ( nevents > 0 ){
4479      //      //
4480      for (UInt_t i=0; i < nevents; i++){      for (UInt_t i=0; i < nevents; i++){
# Line 4428  Int_t PamelaDBOperations::insertCALOPULS Line 4484  Int_t PamelaDBOperations::insertCALOPULS
4484          if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){          if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){
4485            valid = 1;            valid = 1;
4486            if ( cp1->pperror[section] ) valid = 0;            if ( cp1->pperror[section] ) valid = 0;
4487              //
4488              // check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.]
4489              //
4490              totpul = 0.;
4491              Int_t ns = 0;
4492              if ( section == 2 ) ns = 3;
4493              if ( section == 3 ) ns = 1;
4494              if ( section == 1 ) ns = 2;
4495              for (UInt_t ss=0; ss<96; ss++){
4496                totpul += cp1->calpuls[ns][0][ss];
4497              }
4498              if ( totpul >= 3145632. ){
4499                if ( IsDebug() ) printf(" PULSE1 Section %i totpul %f -  No calibration data! Calorimeter power problems? \n",section,totpul);
4500                valid = 0;
4501              };
4502              //
4503            ph = eh->GetPscuHeader();            ph = eh->GetPscuHeader();
4504            obt = ph->GetOrbitalTime();              obt = ph->GetOrbitalTime();  
4505            pkt = ph->GetCounter();              pkt = ph->GetCounter();  
# Line 4585  Int_t PamelaDBOperations::insertCALOPULS Line 4657  Int_t PamelaDBOperations::insertCALOPULS
4657          if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){          if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){
4658            valid = 1;            valid = 1;
4659            if ( cp2->pperror[section] ) valid = 0;            if ( cp2->pperror[section] ) valid = 0;
4660              //
4661              // check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.]
4662              //
4663              totpul = 0.;
4664              Int_t ns = 0;
4665              if ( section == 2 ) ns = 3;
4666              if ( section == 3 ) ns = 1;
4667              if ( section == 1 ) ns = 2;
4668              for (UInt_t ss=0; ss<96; ss++){
4669                totpul += cp2->calpuls[ns][0][ss];
4670              }
4671              if ( totpul >= 3145632. ){
4672                if ( IsDebug() ) printf(" PULSE2 Section %i totpul %f -  No calibration data! Calorimeter power problems? \n",section,totpul);
4673                valid = 0;
4674              };
4675              //
4676            ph = eh->GetPscuHeader();            ph = eh->GetPscuHeader();
4677            obt = ph->GetOrbitalTime();              obt = ph->GetOrbitalTime();  
4678            pkt = ph->GetCounter();              pkt = ph->GetCounter();  
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.23