UgliDbiScintMdl.cxx

Go to the documentation of this file.

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

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

#include "UgliGeometry/UgliDbiTableDescr.h"

#include <cassert>

ClassImp(UgliDbiScintMdl)

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

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

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

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

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

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


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

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

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

UInt_t UgliDbiScintMdl::HashToIndex(const PlexScintMdlId scintmdlid)
{
   // make a natural index

   UInt_t hvalue; 
   int ierr = 0;

   // VERY WASTEFUL
   // assume 8 modules/plane independent of detector

   int module = scintmdlid.GetModule();
   int plane  = scintmdlid.GetPlane();

   hvalue = (plane << 3) + module;
   if ( module > 7 ) ierr += BIT(0);
   // breaks under compressed far geometry
   // if ( hvalue > 3896-1 ) ierr += BIT(1); 
   const int mxplns = 768; // 323*2;  // 243*2;
   if ( hvalue > (mxplns<<3)-1 ) ierr += BIT(1);

   if (ierr) {
      MSG("Ugli",Msg::kWarning) 
         << " module " << scintmdlid.AsString() 
         << " breaks hashing assumptions " 
         << "    err 0x" << hex << ierr << dec 
         << " : scintmdlid " << scintmdlid
         << " = hash " << hvalue
         << endl;
   }

   return hvalue;
   
}

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

void UgliDbiScintMdl::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 module = -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 == "MDLINPLANE"    ) rs >> module;
         // "MODULE" is obsolete
         else if ( colName == "MODULE"        ) rs >> module;
         else if ( colName == "WIDTH"         ) rs >> fWidth;
         else if ( colName == "TPOSRELPLN"    ) rs >> fTPosRelPln;
         else if ( colName == "LPOSRELPLN"    ) rs >> fLPosRelPln;
         else if ( colName == "ZROTRELPLNRAD" ) rs >> fZRotRelPlnRad;
         else if ( colName == "CLEARLENEAST"  ) rs >> fClearLenEast;
         else if ( colName == "CLEARLENWEST"  ) rs >> fClearLenWest;
         else if ( colName == "WLSLENEAST"    ) rs >> fWlsLenEast;
         else if ( colName == "WLSLENWEST"    ) rs >> fWlsLenWest;
         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;
      rs >> hashv
         >> plane
         >> module
         >> fWidth
         >> fTPosRelPln
         >> fLPosRelPln
         >> fZRotRelPlnRad
         >> fClearLenEast
         >> fClearLenWest
         >> fWlsLenEast
         >> fWlsLenWest;
          
   }

   fModuleId = PlexScintMdlId(det,plane,module,view,coverage);

}

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

void UgliDbiScintMdl::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)fModuleId.GetDetector();
      else if ( colName == "PLANE"         ) ors << fModuleId.GetPlane();
      else if ( colName == "PLNVIEW"       ) 
                        ors << (int)fModuleId.GetPlaneView();
      // "VIEW" is obsolete
      else if ( colName == "VIEW"          ) 
                        ors << (int)fModuleId.GetPlaneView();
      else if ( colName == "COVERAGE"      ) 
                        ors << (int)fModuleId.GetPlaneCoverage();
      else if ( colName == "MDLINPLANE"    ) ors << fModuleId.GetModule();
      // "MODULE" is obsolete
      else if ( colName == "MODULE"        ) ors << fModuleId.GetModule();
      else if ( colName == "WIDTH"         ) ors << fWidth;
      else if ( colName == "TPOSRELPLN"    ) ors << fTPosRelPln;
      else if ( colName == "LPOSRELPLN"    ) ors << fLPosRelPln;
      else if ( colName == "ZROTRELPLNRAD" ) ors << fZRotRelPlnRad;
      else if ( colName == "CLEARLENEAST"  ) ors << fClearLenEast;
      else if ( colName == "CLEARLENWEST"  ) ors << fClearLenWest;
      else if ( colName == "WLSLENEAST"    ) ors << fWlsLenEast;
      else if ( colName == "WLSLENWEST"    ) ors << fWlsLenWest;

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

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

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

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

std::ostream& UgliDbiScintMdl::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   << GetModule() << sep  // within aggregate(=plane) module => row_counter
       << -1          << sep  // Hashval is obsolete
       << GetPlane()            << sep
       << GetModule()           << sep
       << fWidth                << sep
       << fTPosRelPln           << sep
       << fLPosRelPln           << sep
       << fZRotRelPlnRad        << sep
       << fClearLenEast         << sep
       << fClearLenWest         << sep
       << fWlsLenEast           << sep
       << fWlsLenWest           << endl;

  return os;
}

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

const char* UgliDbiScintMdl::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, \
       MDLINPLANE    smallint, \
       WIDTH         float,    \
       TPOSRELPLN    float,    \
       LPOSRELPLN    float,    \
       ZROTRELPLNRAD float,    \
       CLEARLENEAST  float,    \
       CLEARLENWEST  float,    \
       WLSLENEAST    float,    \
       WLSLENWEST    float,    \
       primary key (SEQNO,ROW_COUNTER) \
       )";

   return const_tabledescr;
}


/*    Template for New Member Function

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

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

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

//  None.


}

*/

---