Grid.h

Go to the documentation of this file.

#ifndef GeFiCa_Grid
#define GeFiCa_Grid
namespace GeFiCa { class Points; class FieldLine; class Grid; class Detector; }
#include <stdlib.h>
#include <vector>
class GeFiCa::Points
{
   public:
      std::vector<double> C1; 
      std::vector<double> C2; 
      std::vector<double> C3; 
      std::vector<double> Vp; 
      std::vector<double> Et; 
      std::vector<double> E1; 
      std::vector<double> E2; 
      std::vector<double> E3; 
      std::vector<double> dC1p; 
      std::vector<double> dC1m; 
      std::vector<double> dC2p; 
      std::vector<double> dC2m; 
      std::vector<double> dC3p; 
      std::vector<double> dC3m; 
      size_t GetN() { return C1.size(); } 
};
#include <TNamed.h>
class TGraph;
class GeFiCa::FieldLine : public Points, public TNamed
{
   public:
      FieldLine() : Points(), TNamed("fl", "electric field line"), fGl(0) {};
      TGraph* GetGraph();
   private:
      TGraph* fGl; 
      ClassDef(FieldLine,1);
};
class TTree;
class GeFiCa::Grid : public Points, public TNamed
{
   public:
      size_t N1; 
      size_t N2; 
      size_t N3; 
      std::vector<double> Src; 
      size_t MaxIterations;    
      double RelaxationFactor; 
      double Tolerance;        
      size_t Iterations;       

      Grid(size_t n1=0, size_t n2=0, size_t n3=0);
      virtual ~Grid() {};
      virtual void SetupWith(Detector &detector);
      void SuccessiveOverRelax();
      void SolveAnalytically();
      double GetV(double c1, double c2=0, double c3=0) const
      {return GetData(Vp,c1,c2,c3); }
      double GetE(double c1, double c2=0, double c3=0) const
      { return GetData(Et,c1,c2,c3); }
      double GetE1(double c1, double c2=0, double c3=0) const
      { return GetData(E1,c1,c2,c3); }
      double GetE2(double c1, double c2=0, double c3=0) const
      { return GetData(E2,c1,c2,c3); }
      double GetE3(double c1, double c2=0, double c3=0) const
      { return GetData(E3,c1,c2,c3); }
      double GetC();
      TTree* GetTree(bool createNew=false);
      bool IsDepleted();
      Grid& operator*=(double scale);
      Grid& operator+=(Grid& other);
      virtual FieldLine* GetFieldLineFrom(double c1, double c2, bool positive=true);
      virtual FieldLine* GetFieldLineFrom(double c1, double c2, double c3,bool positive=true);
   protected:
      std::vector<bool> fIsFixed; 
      std::vector<bool> fIsDepleted; 
      TTree* fTree; 
      Detector* fDetector; 

      virtual void OverRelaxAt(size_t idx) {};
      size_t GetIdxOfPointToTheRightOf(double c1,
            size_t begin, size_t end) const;
      size_t GetIdxOfPointToTheRightOf(double c1, double c2,
            size_t begin, size_t end) const;
      size_t GetIdxOfPointToTheRightOf(double c1, double c2, double c3,
            size_t begin, size_t end) const;
      virtual double GetData(const std::vector<double> &data,
            double c1, double c2, double c3) const;
      size_t GetIdxOfMaxV();
      size_t GetIdxOfMinV();
      virtual void CalculateE();
      double twopoint(double dataset[2],double tarlocationset,double pointxset[2])const;
      double threepoint(double dataset[3],double tarlocationset[2],double pointxset[3],double pointyset[3])const;
      double fourpoint(double dataset[4],double tarlocationset[2],double pointxset[4],double pointyset[4])const;
      ClassDef(Grid,1);
};
#endif