#ifndef _TCAD
#define _TCAD

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// TCAD                                                                 //
//                                                                      //
// Class for importation of TCAD objects                                //
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include "TMath.h"
#include "nrutil.h"
#include "math.h"
#include <stdio.h>
#include <stdlib.h>
#include "TArray.h"
#include "TArrayI.h"
#include "TArrayF.h"
#include "TArrayD.h"
#include "TString.h"
#include "TH2.h"
#include "TH1.h"
#include <string.h>
#include <TObject.h>

class TCAD : public TObject {

  //private:
 

public: 
  Int_t num;
  TArrayF Gx;
  TArrayF Gy;
  TArrayF Gz;
  TArrayF U;
  TArrayF X;
  TArrayF Y;
  TArrayF Z;
  Float_t MaxX;
  Float_t MinX;
  Float_t MaxY;
  Float_t MinY;
  Float_t MaxZ;
  Float_t MinZ;

  Int_t Index[50][300];
  Int_t IndexNum1,IndexNum2;
  TCAD(Char_t *,Char_t *,Int_t=0);
  ~TCAD();
 void Sort2(Int_t , Int_t, Char_t *);
 void SortAll();
 void GetMaxMin();
  void Print(Int_t=0, Int_t=-1);
  Float_t FieldValue(Float_t ,Float_t , Float_t );
  Int_t FindPlane(Int_t ,Float_t ,Float_t ,Int_t *);
  TH1F *Slice1D(Option_t *, Axis_t,Axis_t);
  TH2F *Slice2D(Option_t *, Axis_t );
  //  void GetField();
  ClassDef(TCAD,1) 
};

#endif