/**
 * \file ExtTrack.h
 * \author Elena Vannuccini
 */
#ifndef exttrack_h
#define exttrack_h

#include <TrkLevel2.h>

class ExtTrack : public TObject {

private:

public:

    	
    Float_t al[5];           ///<TRACK STATE VECTOR 
    Float_t coval[5][5];     ///<covariance matrix 
    Float_t chi2;          
    Int_t  nstep;           
    Int_t nplanes;
    Int_t*   xgood ;        //[nplanes]       
    Int_t*   ygood ;        //[nplanes]        
    Float_t* xm ;           //[nplanes]
    Float_t* ym ;           //[nplanes]
    Float_t* zm ;           //[nplanes]
    Float_t* resx ;         //[nplanes]
    Float_t* resy ;         //[nplanes]
    Float_t* xv ;           //[nplanes]
    Float_t* yv ;           //[nplanes]
    Float_t* zv ;           //[nplanes]
    Float_t* axv ;          //[nplanes]
    Float_t* ayv ;          //[nplanes]
    Float_t* dedx_x ;       //[nplanes]
    Float_t* dedx_y ;       //[nplanes]


    ExtTrack();
    ExtTrack(Int_t dim);
    ExtTrack(const ExtTrack&);

    ~ExtTrack(){ Delete(); };

    void Clear(Option_t* option="");
    void Delete();
    void Dump();

    void FitReset();
    void SetDimension(Int_t dim);
    void SetMeasure(double *xmeas, double *ymeas, double *zmeas);
    void SetResolution(double *rx, double *ry);
    void SetGood(Int_t *xg, Int_t *yg);

    void FillMiniStruct(cMini2track& track);
    void SetFromMiniStruct(cMini2track *track);


    Bool_t XGood(int ip){ return xgood[ip]; };
    Bool_t YGood(int ip){ return ygood[ip]; };

    Int_t GetNX(){Int_t n=0; for(Int_t i=0; i<nplanes; i++)n+=(Int_t)XGood(i); return n;}; 
    Int_t GetNY(){Int_t n=0; for(Int_t i=0; i<nplanes; i++)n+=(Int_t)YGood(i); return n;};
    Int_t GetNXY(){Int_t n=0; for(Int_t i=0; i<nplanes; i++)n+=(Int_t)YGood(i)*XGood(i); return n;};
    Int_t GetNtot(){return GetNX()+GetNY();};

    ClassDef(ExtTrack,1);

};

#endif