#ifndef _SiCapCal
#define _SiCapCal

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// SiCapCal                                                             //
//                                                                      //
// Description of the scenario class                                    //
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include "TObject.h"
#include "TF1.h"
#include <TMath.h>
#include "TArray.h"
#include "TArrayI.h"
#include "TArrayF.h"
#include "TGraph.h"
#include <stdio.h>
#include <stdlib.h>
#include "TObject.h"
#include <nrutil.h>


void nrerror(char error_text[]);
void free_ivector(int *v, long nl, long nh);
int *ivector(long nl, long nh);
float *vector(long nl, long nh);
float **matrix(long nrl, long nrh, long ncl, long nch);
void free_vector(float *v, long nl, long nh);

class SiCapCal {

 private:
  void lubksb(float **, int , int *, float []);
  void ludcmp(float **, int , int *, float *);
  public:
  Int_t Nread;          //Number of read-out strips
  Int_t Niter;          //Number of interpolation strips
  Double_t Pitch;       //read-out pitch
  Double_t Width;       //strip width
  Double_t Thick;       //detector thichness
  Double_t Lenght;      //strip lenght
  Int_t Nwir;           //number of wire used in the approximation
  Int_t Sigstr;         //strip with the signal
  Int_t ReturnCap;      //Do you want the field of Capatcitance

   SiCapCal(Int_t=15, Int_t=2, Double_t=200, Double_t=10,Double_t=280,Double_t=1,Int_t=10,Int_t=1);
  ~SiCapCal() {};
  Double_t *CapCal();
  Double_t *CapCal(int,int,double,double,double,double,int,int);
  ClassDef(SiCapCal,1) 
    };
#endif