BfldHandlerVoronoi.cxx

Go to the documentation of this file.

// $Id: BfldHandlerVoronoi.cxx,v 1.10 2005/02/03 22:35:49 rhatcher Exp $
// 
// BfldHandlerVoronoi 
// 
// BfldHandlerVoronoi searches and interpolates BfldMeshVoronoi grids
//
// Author:  R. Hatcher 2000.06.20
//

#include "BField/BfldHandlerVoronoi.h"

#include "Conventions/Munits.h"
#include "BField/BfldMapVoronoi.h"
#include "BField/BfldMeshVoronoi.h"

#include "BField/BFLCache.h"

#include "MessageService/MsgService.h"
#include "MessageService/MsgFormat.h"
CVSID("$Id: BfldHandlerVoronoi.cxx,v 1.10 2005/02/03 22:35:49 rhatcher Exp $");

ClassImp(BfldHandlerVoronoi)

//_____________________________________________________________________________
BfldHandlerVoronoi::BfldHandlerVoronoi()
   : BfldHandler(), fCostasHandler(0)
{
   MSG("Bfld",Msg::kWarning) 
      << "BfldHandlerVoronoi::ctor" << endl;
}

//_____________________________________________________________________________
BfldHandlerVoronoi::~BfldHandlerVoronoi()
{

   if (fCostasHandler) delete fCostasHandler;
}

//_____________________________________________________________________________
void BfldHandlerVoronoi::SetInterpMethod(const BfldInterpMethod::InterpMethod_t method) 
{
   // set the interpolation method for this handler
   // the Voronoi specialization translates "kDefault" to "kNatural"
   // 
   // and should also disallow illegal choices

   static Int_t nwarn = 10;

   BfldInterpMethod::InterpMethod_t chosen = method;

   switch (chosen) {
   case BfldInterpMethod::kBilinear:  // not yet!
     break;
   case BfldInterpMethod::kPlanar:
   case BfldInterpMethod::kPlanarVec:
      // can't do these methods
      if (nwarn > 0) {
         nwarn--;
         MSG("Bfld",Msg::kWarning)
            << "BfldHandlerRect2d::SetInterpMethod '"
            << BfldInterpMethod::AsString(chosen) 
            << "' not supported"<< endl;
      }
      // fall through to our choice of default
   case BfldInterpMethod::kDefault:
      // pick a default method
      chosen = BfldInterpMethod::kNatural;
      break;
   default:
      // no need to adjust anything for most methods
      break;
   }
   
   BfldHandler::SetInterpMethod(chosen); // do the normal thing

   SetCostasInterpMethod();  // make sure Costas handler is notified
   
}

//_____________________________________________________________________________
void BfldHandlerVoronoi::SetMap(BfldMap* bmap)
{
  BfldHandler::SetMap(bmap); // do the normal thing
  SetCostasHandler(bmap->GetGrid());  // make sure Costas handler is attached
}

//_____________________________________________________________________________
void BfldHandlerVoronoi::SetMesh(BfldMesh* mesh)
{
  BfldHandler::SetMesh(mesh); // do the normal thing
  SetCostasHandler(mesh->GetGrid());  // make sure Costas handler is attached
}

//_____________________________________________________________________________
void BfldHandlerVoronoi::SetCostasHandler(BfldGrid::Grid_t grid)
{

  // this isn't really the right way to do this ...but a hack
  if ( ! fCostasHandler ) {

     MSG("Bfld",Msg::kDebug) 
        << "BfldHandlerVoronoi::SetCostasHandler create BFLHander" << endl;

    switch (grid) {
    case BfldGrid::kNearCoarseV: fCostasHandler = new BFLHandler(0);
      break;
    case BfldGrid::kFarCoarseV:  fCostasHandler = new BFLHandler(1);
      break;
    default:
      MSG("Bfld",Msg::kFatal) 
        << "BfldHandlerVoronoi::SetCostasHandler couldn't attach BFLHandler for " 
        << BfldGrid::AsString(grid) << " grid " << endl;
      break;
    }
    
    SetCostasInterpMethod();

    //fCostasHandler->SetLastPolygAsSeed();

    //  BFLHandler creates a BFLInterpolation which creates a BFLCache
    //  no need to explicitly do so -- it is a memory leak
    //BFLCache *mycached = new BFLCache;
    //fCostasHandler->SetCache(mycached);
    
  }

}

//_____________________________________________________________________________
void BfldHandlerVoronoi::SetCostasInterpMethod()
{
   // inform fCostasHandler about the interpolation approach requested

//   if (fCostasHandler) fCostasHandler->SetInterpolant((int)fInterpMethod);
   if (fCostasHandler) fCostasHandler->SetInterpolant(fInterpMethod);
}

//_____________________________________________________________________________
TVector3 BfldHandlerVoronoi::GetBFieldMeshCoord(const TVector3 &xyz)
{
   // search for the right location
   // and then interpolate

   // Costas expects position in mm ...
   // Costas currently returns in Tesla ... 

   // convert position to "mm"
   TVector3 xyz_mm = xyz * (1./Munits::mm);
   // return field in "Tesla"
   return fCostasHandler->GetBField(xyz_mm,kTRUE) * Munits::tesla;
}

//_____________________________________________________________________________

---