UgliDbiStructHash.cxx

Go to the documentation of this file.

// $Id: UgliDbiStructHash.cxx,v 1.6 2005/08/26 19:00:02 rhatcher Exp $
//
// UgliDbiStructHash
//
// Package: UgliGeometry
//
// R. Hatcher  2002-08-02
//
// Concept:
//   
// Purpose:
//

#include "UgliGeometry/UgliDbiStructHash.h"
#include "MessageService/MsgService.h"

ClassImp(UgliDbiStructHash)

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

CVSID("$Id: UgliDbiStructHash.cxx,v 1.6 2005/08/26 19:00:02 rhatcher Exp $ CVSID_DBIRESULTPTR ");

UInt_t arbitraryPln(Detector::Detector_t det,
                    PlaneView::PlaneView_t view,
                    PlaneCoverage::PlaneCoverage_t cover) {
  // return an arbitrary plane ... execpt in case of veto stuff
  if (Detector::kFar != det) return 0;
  // more work if far
  if (view < PlaneView::kVSTopFlat && cover < PlaneCoverage::kVScN) return 0;
  // should check that both view,cover are veto related
  return 512;
}


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

bool operator==(const UgliDbiStructHash::viewcover_t& vc1, 
                const UgliDbiStructHash::viewcover_t& vc2) {
  if ( vc1.view  != vc2.view  ) return false;
  if ( vc1.cover != vc2.cover ) return false;
  return true;
}

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

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

UInt_t UgliDbiStructHash::HashAsPlane() const
{ UShort_t maxpart = 0; return HashToIndex(maxpart); }

UInt_t UgliDbiStructHash::HashAsScintMdl() const
{ 
  UInt_t apln = arbitraryPln(fDetector,fPlaneView,fPlaneCoverage);
  PlexPlaneId plnid(fDetector,apln,false,fPlaneView,fPlaneCoverage);
  UShort_t maxpart = plnid.GetNumScintMdls();
  return HashToIndex(maxpart); 
}

UInt_t UgliDbiStructHash::HashAsStrip() const
{ 
  UInt_t apln = arbitraryPln(fDetector,fPlaneView,fPlaneCoverage);
  PlexPlaneId plnid(fDetector,apln,false,fPlaneView,fPlaneCoverage);
  UShort_t maxpart = plnid.GetNumStrips();
  return HashToIndex(maxpart); 
}

UInt_t UgliDbiStructHash::HashToIndex(UShort_t maxpart) const
{
   // make a natural index where the "part" takes on values [0..maxpart]

   // Different hashes based on detector

  static viewcover_t gNearCombos[] = {
    {PlaneView::kUnknown,PlaneCoverage::kNoActive},
    {PlaneView::kUnknown,PlaneCoverage::kUnknown},
    {PlaneView::kU,PlaneCoverage::kNearPartial},
    {PlaneView::kU,PlaneCoverage::kNearFull},
    {PlaneView::kV,PlaneCoverage::kNearPartial},
    {PlaneView::kV,PlaneCoverage::kNearFull}
  };
  static unsigned int gNumNearCombos =
    sizeof(gNearCombos)/sizeof(viewcover_t);

  static viewcover_t gFarCombos[] = {
    {PlaneView::kUnknown,PlaneCoverage::kNoActive},
    {PlaneView::kUnknown,PlaneCoverage::kUnknown},
    {PlaneView::kU,PlaneCoverage::kComplete},
    {PlaneView::kV,PlaneCoverage::kComplete},

    // veto shield top (C/c-N, E/e-S)

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVScN},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVScN},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVScN},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSCN},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSCN},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSCN},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSeS},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSeS},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSeS},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSES},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSES},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSES},

    // veto shield walls (C/c-N, E/e-S)

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVScN},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVScN},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVScN},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSCN},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSCN},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSCN},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSeS},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSeS},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSeS},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSES},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSES},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSES},

    // veto shield top (C/c-S, E/e-N)

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVScS},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVScS},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVScS},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSCS},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSCS},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSCS},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSeN},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSeN},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSeN},

    {PlaneView::kVSTopFlat,       PlaneCoverage::kVSEN},
    {PlaneView::kVSTopEastSlant,  PlaneCoverage::kVSEN},
    {PlaneView::kVSTopWestSlant,  PlaneCoverage::kVSEN},

    // veto shield walls (C/c-S, E/e-N)

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVScS},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVScS},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVScS},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSCS},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSCS},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSCS},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSeN},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSeN},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSeN},

    {PlaneView::kVSWallOnEdge,    PlaneCoverage::kVSEN},
    {PlaneView::kVSWallEastSlant, PlaneCoverage::kVSEN},
    {PlaneView::kVSWallWestSlant, PlaneCoverage::kVSEN},


  };
  static unsigned int gNumFarCombos =
    sizeof(gFarCombos)/sizeof(viewcover_t);


  static viewcover_t gCalDetCombos[] = {
    {PlaneView::kUnknown,PlaneCoverage::kNoActive},
    {PlaneView::kUnknown,PlaneCoverage::kUnknown},
    {PlaneView::kU,PlaneCoverage::kComplete},
    {PlaneView::kV,PlaneCoverage::kComplete},
    {PlaneView::kA,PlaneCoverage::kComplete},
    {PlaneView::kB,PlaneCoverage::kComplete}
  };
  static unsigned int gNumCalDetCombos =
    sizeof(gCalDetCombos)/sizeof(viewcover_t);


  static viewcover_t gAllCombos[] = {
    {PlaneView::kUnknown,PlaneCoverage::kNoActive},
    {PlaneView::kUnknown,PlaneCoverage::kUnknown},
    {PlaneView::kU,PlaneCoverage::kComplete},
    {PlaneView::kV,PlaneCoverage::kComplete},
    {PlaneView::kA,PlaneCoverage::kComplete},
    {PlaneView::kB,PlaneCoverage::kComplete}
  };
  static unsigned int gNumAllCombos =
    sizeof(gAllCombos)/sizeof(viewcover_t);

   switch (fDetector) {
   case (Detector::kNear):   
     return HashToIndex(gNearCombos,
                        gNumNearCombos,
                        maxpart);
   case (Detector::kFar):    
     return HashToIndex(gFarCombos,
                        gNumFarCombos,
                        maxpart);
   case (Detector::kCalDet): 
     return HashToIndex(gCalDetCombos,
                        gNumCalDetCombos,
                        maxpart);
   default:                 
     return HashToIndex(gAllCombos,
                        gNumAllCombos,
                        maxpart);
   }
}

UInt_t UgliDbiStructHash::HashToIndex(viewcover_t* vclist, UShort_t nlist, 
                                      UShort_t /* maxpart */ ) const
{

  viewcover_t thisvc;
  thisvc.view  = fPlaneView;
  thisvc.cover = fPlaneCoverage;

  int indx_fnd = -1;
  for (unsigned int indx = 0; indx<nlist; ++indx) {
    if (thisvc == vclist[indx]) 
      { indx_fnd = indx; break; }
  }
  if (indx_fnd == -1) {
     MSG("Ugli",Msg::kWarning)
       << "Hash det " << Detector::AsString(fDetector) 
       << " view " << PlaneView::AsString(fPlaneView)
       << " cover " << PlaneCoverage::AsString(fPlaneCoverage)
       << " breaks hashing scheme" << endl;
     indx_fnd = 0;
  }

  UInt_t hvalue = fPart*nlist + indx_fnd;
  return hvalue;

}

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

---