//-----------------------------------------------------------
// File and Version Information:
// 
//
// Description:
//      Implementation of class TpcEFieldCylLs
//      see TpcEFieldCylLs.h for details
//
// Environment:
//      Software developed for the GEM-TPC Detector at FOPI.
//
// Author List:
//      Martin Berger    TUM            (original author)
//
//
//-----------------------------------------------------------

// This Class' Header ------------------
#include "TpcEFieldCylLs.h"
#include "TpcEFieldCyl.h"
#include "TMath.h"
#include <iostream>

ClassImp(TpcEFieldCylLs)

// Class Member definitions -----------

TpcEFieldCylLs::TpcEFieldCylLs()
:TpcEFieldCyl()
{
  lastStripField=200.;
  lastStripPos=3.5;
  driftField = TVector3(0.,0.,400.);
}

// TpcEFieldCylLs::TpcEFieldCylLs(TString name)
// :TpcEFieldCyl()
// {
//   fname=name;
//   lastStripField=200.;
//   lastStripPos=3.5;
//   driftField = TVector3(0.,0.,400.);
// }


TpcEFieldCylLs::~TpcEFieldCylLs()
{;}

TVector3 TpcEFieldCylLs::value(const TVector3& point) const
{
  return value2(point);
}

TVector3 TpcEFieldCylLs::value2(const TVector3& point) const
{
  TVector3 resultField;
  TVector3 resultLsField;
  double signFactor;
  double rad;
  double radfactor;
  double x_along_z, x_along_r;
  double width;

  resultField=driftField;
  rad=TMath::Sqrt(point[0]*point[0]+point[1]*point[1]);
  if (rad<5)
    return TVector3(0.,0.,0.);
  else
    rad-=5;

  signFactor=1;
  if (lastStripField<driftField.Mag())
    signFactor=-1;
  if (lastStripField==driftField.Mag())
    signFactor=0;

  radfactor=1  ;

  if (point[2]<lastStripPos)
    {
      //      resultField += TVector3(0.,0.,signFactor*lastStripField);
      //      resultField[2] +=lastStripField*radfactor*TMath::Gaus(0,lastStripPos,2)*TMath::Gaus(0,0,0.1/3);
      //      x_along_z=lastStripField*TMath::Gaus(point[2],lastStripPos,double(lastStripPos/3.));
      x_along_z=TMath::Exp(point[2]-lastStripPos);
    }  
  else
    {
      //x_along_z=0.;
      x_along_z=TMath::Gaus(point[2],lastStripPos,0.1/2.);
    }

  x_along_r=TMath::Exp(-rad);
  double angle=TMath::ATan(point[1]/point[0]);
  
  resultField[0]+=signFactor*lastStripField*x_along_z*x_along_r*TMath::Cos(angle);
  resultField[1]+=signFactor*lastStripField*x_along_z*x_along_r*TMath::Sin(angle);

  //  std::cout<<angle<<" "<<TMath::Cos(angle)<<" "<<TMath::Sin(angle)<<" "<<point[0]<<" "<<point[1]<<std::endl;

  //  if (rad < 5.25 || rad > 15)
  //    return 0.0;

  //  radfactor=TMath::Abs((10.-rad)/5.);


  //  resultLsField[2]=lastStripField*radfactor*TMath::Gaus(point[2],lastStripPos,2)*TMath::Gaus(rad,0,0.1/3);
  //  resultField[2]+=resultLsField[2]*signFactor;

  

  return resultField;
}

// void TpcEFlielCylLs::readFile()
// {

//  // try to find and open the file
//   std::ifstream infile(fname, std::fstream::in);
//   if (!infile.good()) 
//   {
//     Fatal("TpcEFieldCylLoader::load()","field-file could not be found.");
//     return;
//   }

//   ffieldmapR  = new std::vector<double>;
//   ffieldmapZ  = new std::vector<double>;
//   ffieldmapEr = new std::vector<double>;
//   ffieldmapEz = new std::vector<double>;

//   double R,Z,Er,Ez;

//   while (infile.good())
//     {
//       infile >> R >> Z >> Er >> Ez;
//     }
// }