BFLInterpolation.h

Go to the documentation of this file.

//                                                                   //
//authors: Costas Andreopoulos, Elias Athanasopoulos, George Tzanakos//
//                Physics Department, University of Athens           //

#ifndef BFLINTERPOLATION_H
#define BFLINTERPOLATION_H

#include <iostream>

#include "TObject.h"
#include "TVector3.h"
#include "TObjArray.h"

//#include "BField/BFLWingedEdge.h"
//#include "BField/BFLVorOperator.h"
//#include "BField/BFLCache.h"
#include "BField/BfldInterpMethod.h"
#include "BField/BfldGrid.h"

class BFLWingedEdge;
class BFLVorOperator;
class BFLCache;
class BFLNode;
class BFLNode2ACell;
class TIntList;

const Int_t kFail = 1;
const Int_t kMaxVtx = 8;

//const Int_t kNaturalNeighbors = 1;

class BFLInterpolation: public TObject
{
public: 

  //BFLInterpolation();
  BFLInterpolation(TObjArray * bfield, BFLWingedEdge * voronoi);
  ~BFLInterpolation();

  TVector3 GetB(const TVector3& PosVector,
                const Bool_t updateSeedPolygon = kFALSE);

  virtual BFLCache* GetCache(void) { return fCache; }
  virtual void SetInterpolant(BfldInterpMethod::InterpMethod_t intp) {
      fHowToInterpolate = intp; 
  }
  virtual void SetGrid(BfldGrid::Grid_t grid) { fGrid = grid;  }
  virtual void SetVoronoi(BFLWingedEdge * vor) { fVor = vor; }
  virtual void SetBField(TObjArray * bfield) { fBField = bfield; }
  virtual void SetCache(BFLCache * cache) { fCache = cache; }
  virtual void SetSearchSeed(Int_t seed) { fSearchSeed = seed; }
  virtual Int_t GetLastPolygon(void) { return fCurrentPolygID; }
    
private: 

  Int_t FormVoronoiCell(BFLNode * point);
  Double_t CellArea(TObjArray * vertices);
  TVector3 NNInterpolation(BFLNode * point);
  TVector3 CNInterpolation(BFLNode * point);
  TVector3 BilinearInterpolation(BFLNode * point);
  TVector3 PlanarInterpolation(BFLNode *point); //Added by A.G.

  void ResetForNextQuery(void);
  Bool_t IsInsideANSYSCell(BFLNode2ACell * TableEntry, BFLNode * point);
  void IntBI1stOrder(void);

  BfldInterpMethod::InterpMethod_t   fHowToInterpolate;
  BfldGrid::Grid_t                   fGrid;
  Int_t                              fCurrentPolygID;
  Int_t                              fSearchSeed;
  Int_t                              fNPolygs;
  TObjArray*                         fBField;
  TObjArray*                         fGenerators;
  BFLWingedEdge*                     fVor;
  BFLVorOperator*                    fOp;  
  BFLCache*                          fCache;
  TIntList*                          fPolygons;
  TObjArray*                         fVertices;

  ClassDef(BFLInterpolation,0)   // Field Vectors interpolation class  
};

#endif

---