UgliDbiStrip.cxx

Go to the documentation of this file.

// $Id: UgliDbiStrip.cxx,v 1.21 2007/03/09 05:59:47 rhatcher Exp $
//
// UgliDbiStrip
//
// Package: UgliGeometry
//
// R. Hatcher  2001-08-08
//
// Concept:
//   
// Purpose:
//

#include "UgliGeometry/UgliDbiStrip.h"
#include "MessageService/MsgService.h"
#include "DatabaseInterface/DbiOutRowStream.h"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiValidityRec.h"

#include "UgliGeometry/UgliDbiTableDescr.h"

#include <cassert>

ClassImp(UgliDbiStrip)

//   Definition of static data members
//   *********************************

CVSID("$Id: UgliDbiStrip.cxx,v 1.21 2007/03/09 05:59:47 rhatcher Exp $ CVSID_DBIRESULTPTR ");
Bool_t UgliDbiStrip::fgDefensiveUnpkg = kFALSE;

//  Instantiate associated Result Pointer class.
//  *******************************************

#include "DatabaseInterface/DbiResultPtr.tpl"
template class  DbiResultPtr<UgliDbiStrip>;

#include "DatabaseInterface/DbiWriter.tpl"
template class  DbiWriter<UgliDbiStrip>;

// Definition of member functions (alphabetical order)
// ***************************************************

//.....................................................................

std::ostream& operator<<(std::ostream& os, const UgliDbiStrip& uds)
{ return uds.FormatToOStream(os); }

//.....................................................................

UInt_t UgliDbiStrip::HashToIndex(const PlexStripEndId seid)
{
   // make a natural index

   UInt_t hvalue; 
   int ierr = 0;

   // VERY WASTEFUL
   // assume 191 strips/plane
   // all planes occupied (ie. no sparsification of near spectrometer)

   int strip = seid.GetStrip();

   hvalue = seid.GetPlane()*192 + strip;
   if ( strip > 191 ) ierr += BIT(0);

   if (ierr) MSG("Ugli",Msg::kWarning) 
      << seid.AsString("c")
      << " breaks hashing assumptions " 
      << "    err 0x" << hex << ierr << dec 
      << endl;

   return hvalue;
}

//.....................................................................

void UgliDbiStrip::Fill(DbiResultSet& rs,
                                      const DbiValidityRec* vrec) {
//
//
//  Purpose:  Fill object from Result Set
//
//  Arguments: 
//    rs           in    Result Set used to fill object
//    vrec         in    Associated validity record (or 0 if filling
//                                                    DbiValidityRec)
//
//  Return:    
//
//  Contact:   R. Hatcher
//
//  Specification:-
//  =============
//
//  o Fill object from current row of Result Set.

//  Program Notes:-
//  =============

//  None.

   Detector::Detector_t det =
      (Detector::Detector_t) 
      ((vrec) ? vrec->GetVldRange().GetDetectorMask() : 0);
   
   Int_t plane  = -1;
   Int_t strip  = -1;
   PlaneView::PlaneView_t view             = PlaneView::kUnknown;
   PlaneCoverage::PlaneCoverage_t coverage = PlaneCoverage::kUnknown;
   Int_t tmp;

   if (fgDefensiveUnpkg) {

      // unpack allowing for various column names and ordering

      Int_t numCol = rs.NumCols();
      //  The first column (SeqNo) has already been processed.
      for (Int_t curCol = rs.HasRowCounter() ? 3 : 2; curCol <= numCol; ++curCol) {
         string colName = rs.CurColName();
         if      ( colName == "HASHVAL"       ) rs.IncrementCurCol();
         // "HASH" is obsolete
         else if ( colName == "HASH"          ) rs.IncrementCurCol();
         else if ( colName == "DETECTOR"      ) {
            rs >> tmp; det = (Detector::Detector_t)tmp;
         }
         else if ( colName == "PLANE"         ) rs >> plane;
         else if ( colName == "PLNVIEW"       ) { 
            rs >> tmp; view = (PlaneView::PlaneView_t)tmp;
         }
         // "VIEW" is obsolete
         else if ( colName == "VIEW"          ) { 
            rs >> tmp; view = (PlaneView::PlaneView_t)tmp;
         }
         else if ( colName == "COVERAGE"      ) {
            rs >> tmp; coverage = (PlaneCoverage::PlaneCoverage_t)tmp;
         }
         else if ( colName == "STRIP"         ) rs >> strip;
         else if ( colName == "TPOSRELMDL"    ) rs >> fTPosRelMdl;
         else if ( colName == "LPOSRELMDL"    ) rs >> fLPosRelMdl;
         else if ( colName == "ZROTRELMDLRAD" ) rs >> fZRotRelMdlRad;
         else if ( colName == "TCURVFORM"     ) rs >> fTCurvForm;
         else if ( colName == "TCURVPARM0"    ) rs >> fTCurvParm[0];
         else if ( colName == "TCURVPARM1"    ) rs >> fTCurvParm[1];
         else if ( colName == "TCURVPARM2"    ) rs >> fTCurvParm[2];
         else if ( colName == "TCURVPARM3"    ) rs >> fTCurvParm[3];
         else if ( colName == "ZCURVFORM"     ) rs >> fZCurvForm;
         else if ( colName == "ZCURVPARM0"    ) rs >> fZCurvParm[0];
         else if ( colName == "ZCURVPARM1"    ) rs >> fZCurvParm[1];
         else if ( colName == "ZCURVPARM2"    ) rs >> fZCurvParm[2];
         else if ( colName == "ZCURVPARM3"    ) rs >> fZCurvParm[3];
         else {
            static int nmsg = 25;
            if (nmsg-- > 0) MSG("Ugli",Msg::kWarning) 
               << "Ignoring column " << rs.CurColNum() 
               << " (of " << rs.NumCols() << ")"
               << " \""<< colName << "\"; not part of "
               << ClassName() << endl;
            rs.IncrementCurCol();
         }
      }
   } // end defensive unpacking
   else {
      Int_t hashv, tmpview, tmpcoverage;
      rs >> hashv
         >> plane
         >> strip
         >> tmpview
         >> tmpcoverage
         >> fTPosRelMdl
         >> fLPosRelMdl
         >> fZRotRelMdlRad
         >> fTCurvForm
         >> fZCurvForm
         >> fTCurvParm[0]
         >> fTCurvParm[1]
         >> fTCurvParm[2]
         >> fTCurvParm[3]
         >> fZCurvParm[0]
         >> fZCurvParm[1]
         >> fZCurvParm[2]
         >> fZCurvParm[3];

      view = (PlaneView::PlaneView_t)tmpview;
      coverage = (PlaneCoverage::PlaneCoverage_t)tmpcoverage;

   }

   fStripId = 
      PlexStripEndId(det,plane,strip,StripEnd::kWhole,
                     StripEnd::kWhole,view,coverage);
}

//.....................................................................

void UgliDbiStrip::Store(DbiOutRowStream& ors,
                         const DbiValidityRec* /* vrec */) const {
//
//
//  Purpose:  Stream object to output row stream
//
//  Arguments: 
//    ors          in     Output row stream.
//    vrec         in     Associated validity record (or 0 if filling
//                                                    DbiValidityRec)
//
//  Return:    
//
//  Contact:   R. Hatcher
//
//  Specification:-
//  =============
//
//  o  Stream object to output row stream.

//  Program Notes:-
//  =============

//  None.

   const char* tabledescr = GetTableDescr();
   vector< pair<string,string> > components = 
      UgliDbiTableDescr::ParseTableDescr(tabledescr);
   int n = components.size();

   for (int i=0; i<n; i++) {
      pair<string,string> colPair = components[i];
      string colName = colPair.first;
      if      ( colName == "HASHVAL"       ) ors << (int)GetIndex(0);
      // "HASH" is obsolete
      else if ( colName == "HASH"          ) ors << (int)GetIndex(0);
      else if ( colName == "DETECTOR"      ) 
                        ors << (int)fStripId.GetDetector();
      else if ( colName == "PLANE"         ) ors << fStripId.GetPlane();
      else if ( colName == "PLNVIEW"       ) 
                        ors << (int)fStripId.GetPlaneView();
      // "VIEW" is obsolete
      else if ( colName == "VIEW"          ) 
                        ors << (int)fStripId.GetPlaneView();
      else if ( colName == "COVERAGE"      ) 
                        ors << (int)fStripId.GetPlaneCoverage();
      else if ( colName == "STRIP"         ) ors << fStripId.GetStrip();
      else if ( colName == "TPOSRELMDL"    ) ors << fTPosRelMdl;
      else if ( colName == "LPOSRELMDL"    ) ors << fLPosRelMdl;
      else if ( colName == "ZROTRELMDLRAD" ) ors << fZRotRelMdlRad;
      else if ( colName == "TCURVFORM"     ) ors << fTCurvForm;
      else if ( colName == "TCURVPARM0"    ) ors << fTCurvParm[0];
      else if ( colName == "TCURVPARM1"    ) ors << fTCurvParm[1];
      else if ( colName == "TCURVPARM2"    ) ors << fTCurvParm[2];
      else if ( colName == "TCURVPARM3"    ) ors << fTCurvParm[3];
      else if ( colName == "ZCURVFORM"     ) ors << fZCurvForm;
      else if ( colName == "ZCURVPARM0"    ) ors << fZCurvParm[0];
      else if ( colName == "ZCURVPARM1"    ) ors << fZCurvParm[1];
      else if ( colName == "ZCURVPARM2"    ) ors << fZCurvParm[2];
      else if ( colName == "ZCURVPARM3"    ) ors << fZCurvParm[3];

      else if ( colName == "ROW_COUNTER"   ) {;}
      else if ( colName == "SEQNO"         ) {;}
      else {
         MSG("Ugli",Msg::kFatal)
            << "Column name '" << colName << "' not correctly handled" << endl;
         assert(0);
      }
   }
}

//.....................................................................

void UgliDbiStrip::Print(Option_t* option) const
{
    FormatToOStream(cout,option);
}

//.....................................................................

std::ostream& UgliDbiStrip::FormatToOStream(std::ostream& os,
                                            Option_t* option,
                                            const DbiValidityRec* dvr) const
{
//
//
//  Purpose:  Print the current row
//
//  Return:   nothing
//
//  Contact:   R. Hatcher
//
//  Specification:-
//  =============
//

//  Program Notes:-
//  =============

//  None.

  string opt = option;

  if ( opt.find("H") != string::npos ) {
     string tabledescr = GetTableDescr();
     tabledescr = UgliDbiTableDescr::TextTableDescrLine(tabledescr);
     os << tabledescr << endl;
  }
  else if ( opt.find("h") != string::npos ) {
    string tabledescr = GetTableDescr();
    vector< pair<string,string> > components = 
        UgliDbiTableDescr::ParseTableDescr(tabledescr);
    for (unsigned int ic = 0; ic < components.size(); ++ic) 
        os << components[ic].first << '\t';
    os << endl;
  }

  char sep = '\t';
  if ( opt.find(",") != string::npos ) sep = ',';

  if ( ! dvr ) {
      os << "SeqNo"         << sep;  // no way of knowing the SEQNO?
  }
  else {
      os << dvr->GetSeqNo() << sep;
  }

  os   << GetStrip() << sep  // within aggregate(=plane) strip => row_counter
       << -1         << sep  // Hashval is obsolete
       << fStripId.GetPlane()              << sep
       << fStripId.GetStrip()              << sep
       << (int)fStripId.GetPlaneView()     << sep
       << (int)fStripId.GetPlaneCoverage() << sep
       << fTPosRelMdl           << sep
       << fLPosRelMdl           << sep
       << fZRotRelMdlRad        << sep
       << fTCurvForm            << sep
       << fZCurvForm            << sep
       << fTCurvParm[0]         << sep
       << fTCurvParm[1]         << sep
       << fTCurvParm[2]         << sep
       << fTCurvParm[3]         << sep
       << fZCurvParm[0]         << sep
       << fZCurvParm[1]         << sep
       << fZCurvParm[2]         << sep
       << fZCurvParm[3]         << endl;

  return os;
}

//.....................................................................

const char* UgliDbiStrip::GetTableDescr()
{
//
//
//  Purpose:  Return a string describing rows in the table
//            Used in creating temporary tables
//
//  Return:   const char* to parenthesized comma separated list
//            of column name and type pairs
//
//  Contact:   R. Hatcher
//
//  Specification:-
//  =============
//

//  Program Notes:-
//  =============

//  None.

   const char* const_tabledescr = "(\
       SEQNO         int,      \
       ROW_COUNTER   int,      \
       HASHVAL       int,      \
       PLANE         smallint, \
       STRIP         smallint, \
       PLNVIEW       tinyint,  \
       COVERAGE      tinyint,  \
       TPOSRELMDL    float,    \
       LPOSRELMDL    float,    \
       ZROTRELMDLRAD float,    \
       TCURVFORM     tinyint,  \
       ZCURVFORM     tinyint,  \
       TCURVPARM0    float,    \
       TCURVPARM1    float,    \
       TCURVPARM2    float,    \
       TCURVPARM3    float,    \
       ZCURVPARM0    float,    \
       ZCURVPARM1    float,    \
       ZCURVPARM2    float,    \
       ZCURVPARM3    float,    \
       primary key (SEQNO,ROW_COUNTER) \
       )";

   return const_tabledescr;
}


/*    Template for New Member Function

//.....................................................................

UgliDbiStrip:: {
//
//
//  Purpose:  
//
//  Arguments: 
//    xxxxxxxxx    in    yyyyyy
//
//  Return:    
//
//  Contact:   R. Hatcher
//
//  Specification:-
//  =============
//
//  o 

//  Program Notes:-
//  =============

//  None.


}

*/

---