//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      Data object carrying infomration about a fitted track
//      (residuals, chi2,...)
//      This class is meant to serve as a container which can be
//      used witht the ROOT TreeViewer to produce plots
//
// Environment:
//      Software developed for the GEMTPC Detector
//
// Author List:
//      Martin Berger    TUM            (original author)
//
//
//-----------------------------------------------------------

#ifndef PARTICLE_HH
#define PARTICLE_HH

// Base Class Headers ----------------
#include "../../base/TpcdEdx.h"
#include "TObject.h"
#include "TVector3.h"

// Collaborating Class Headers -------
#include "GFDetPlane.h"

class Particle : public TObject
{

public:
	Particle();
	//Particle(Particle& other);
	//virtual ~Particle();
	void SetFitFlag					 (Int_t flag)	   {fFitFlag=flag;}
	void SetTrackMom				 (TVector3 mom)    {fTrackMom=mom;}
	void SetTrackPos				 (TVector3 pos)	   {fTrackPos=pos;}
	void SetHasTpc					 (Bool_t opt=kTRUE){fHasTpc=opt;}
	void SetHasCdc					 (Bool_t opt=kTRUE){fHasCdc=opt;}
	void SetdEdx					 (TpcdEdx dedx)	   {fdedx=dedx;}


	void appendDetID				 (Int_t id)		   {fDetID.push_back(id);}
	void appendDetName 				 (TString name)    {fDetName.push_back(name);}

	void appendPosXYZ				 (TVector3 vec)    {fPosXYZ.push_back(vec);}
	void appendTrackPosXYZ			 (TVector3 vec)    {fTrackPosXYZ.push_back(vec);}
	void appendTrackPosXYZErr		 (TVector3 vec)    {fTrackPosXYZErr.push_back(vec);}
	void appendTrackMom				 (TVector3 vec)    {fTrackPosMom.push_back(vec);}
	void appendTrackMomErr			 (TVector3 vec)    {fTrackPosMomErr.push_back(vec);}
	void appendResXYZ                (TVector3 vec)    {fResXYZ.push_back(vec);}

	void appendResUV                 (TVector3 vec)    {fResUV.push_back(vec);}
	void appendResA                  (TVector3 vec)    {fResA.push_back(vec);}

	void appendSigXYZ                (TVector3 vec)    {fSigXYZ.push_back(vec);}
	void appendSigA                  (TVector3 vec)    {fSigA.push_back(vec);}
	void appendSigUV                 (TVector3 vec)    {fSigUV.push_back(vec);}

	void appendAxis1                 (TVector3 vec)    {fAxis1.push_back(vec);}
	void appendAxis2                 (TVector3 vec)    {fAxis2.push_back(vec);}
	void appendAxis3                 (TVector3 vec)    {fAxis3.push_back(vec);}

	void appendTrackPlane            (GFDetPlane plane){fTrackPlane.push_back(plane);}
	void appendMidPlane              (GFDetPlane plane){fMidPlane.push_back(plane);}

	void appendClusterSize           (Int_t x)         {fClusterSize.push_back(x);}
	void appendClusterSizeShare      (Double_t x)      {fClusterSizeShare.push_back(x);}
	void append2DClusterSize         (Int_t x)         {f2DClusterSize.push_back(x);}
	void appendFullClusterSize       (TVector3 vec)    {fFullClusterSize.push_back(vec);}
	void appendFullClusterSizeTrack  (TVector3 vec)    {fFullClusterSizeTrack.push_back(vec);}
	void appendFullClusterSizeAxis   (TVector3 vec)    {fFullClusterSizeAxis.push_back(vec);}

	void appendClusterAmp            (Double_t x)      {fClusterAmp.push_back(x);}

//getters
	Int_t   getFitFlag				 ()	   	 	{return fFitFlag				;}
	Int_t   getNumHits				 ()			{return fDetID.size()			;}
	Int_t   nhits				 	 ()			{return fDetID.size()			;}
	TVector3 getTrackPos			 ()			{return fTrackPos				;}
	TVector3 getTrackMom		     () 		{return fTrackMom				;}
	TpcdEdx getdEdx					 ()			{return fdedx					;}
	Bool_t hasTpc					 ()         {return fHasTpc				  	;}
	Bool_t hasCdc					 ()         {return fHasCdc				   	;}

	Int_t   getDetID				 (Int_t i)	  {return fDetID			 [i];}
	TString getDetName 			     (Int_t i)    {return fDetName			 [i];}

	TVector3 getPosXYZ				 (Int_t i)    {return fPosXYZ			    [i];}
	TVector3 getTrackPosXYZ			 (Int_t i)    {if(i<fTrackPosXYZ.size())	return fTrackPosXYZ		[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 getTrackPosXYZErr		 (Int_t i)    {if(i<fTrackPosXYZErr.size())	return fTrackPosXYZErr	[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 getTrackMom			 (Int_t i)    {if(i<fTrackPosMom.size())	return fTrackPosMom		[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 getTrackMomErr			 (Int_t i)    {if(i<fTrackPosMomErr.size())	return fTrackPosMomErr	[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 getResXYZ               (Int_t i)    {if(i<fResXYZ.size()) 		return fResXYZ 			[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 getResUV                (Int_t i)    {return fResUV			    [i];}
	TVector3 getResA                 (Int_t i)    {return fResA			        [i];}

	TVector3 getSigXYZ               (Int_t i)    {return fSigXYZ			    [i];}
	TVector3 getSigA                 (Int_t i)    {return fSigA			        [i];}
	TVector3 getSigUV                (Int_t i)    {return fSigUV			    [i];}

	TVector3 getAxis1                (Int_t i)    {return fAxis1			    [i];}
	TVector3 getAxis2                (Int_t i)    {return fAxis2			    [i];}
	TVector3 getAxis3                (Int_t i)    {return fAxis3			    [i];}

	GFDetPlane getTrackPlane         (Int_t i)    {return fTrackPlane		    [i];}
	GFDetPlane getMidPlane           (Int_t i)    {return fMidPlane		    	[i];}

	Double_t getClusterSizeShare     (Int_t i)    {return fClusterSizeShare   	[i];}
	Int_t    getClusterSize          (Int_t i)    {return fClusterSize		    [i];}
	Int_t    get2DClusterSize        (Int_t i)    {return f2DClusterSize	    [i];}
	TVector3 getFullClusterSize      (Int_t i)    {return fFullClusterSize	    [i];}
	TVector3 getFullClusterSizeTrack (Int_t i)    {return fFullClusterSizeTrack	[i];}
	TVector3 getFullClusterSizeAxis  (Int_t i)    {return fFullClusterSizeAxis	[i];}

	Double_t getClusterAmp           (Int_t i)    {return fClusterAmp			[i];}

	//for compatabalibility
	TVector3 GetRes          	     (Int_t i)    {if(i<fResXYZ.size()) 		return fResXYZ 			[i]; else return TVector3(-9999,-9999,-9999);}
	TVector3 GetResA                 (Int_t i)    {return fResA			        [i];}
	TVector3 GetResUV                (Int_t i)    {return fResUV			    [i];}
	TVector3 GetSigXYZ               (Int_t i)    {return fSigXYZ			    [i];}
	TVector3 GetSigA                 (Int_t i)    {return fSigA			        [i];}
	TVector3 GetSigUV                (Int_t i)    {return fSigUV			    [i];}
	TVector3 GetFullClusterSize      (Int_t i)    {return fFullClusterSize	    [i];}
	TVector3 GetFullClusterSizeAxis  (Int_t i)    {return fFullClusterSizeAxis	[i];}
	TVector3 GetHitPositionXYZ		 (Int_t i)    {return fTrackPosXYZ			[i];}

private:
	// Private Data Members ------------
	//track wise values
	Int_t fFitFlag;
	TVector3 fTrackMom;
	TVector3 fTrackPos;
	Bool_t fHasTpc;
	Bool_t fHasCdc;
	TpcdEdx fdedx;

	//hitpoint wise values
	std::vector<Int_t>		fDetID;
	std::vector<TString>	fDetName;

	std::vector<TVector3>	fPosXYZ;
	std::vector<TVector3>   fTrackPosXYZ;
	std::vector<TVector3>   fTrackPosMom;
	std::vector<TVector3>   fTrackPosXYZErr;
	std::vector<TVector3>   fTrackPosMomErr;


	std::vector<TVector3>   fResXYZ;
	std::vector<TVector3>   fResUV;
	std::vector<TVector3>   fResA;

	std::vector<TVector3>   fSigXYZ;
	std::vector<TVector3>   fSigA;
	std::vector<TVector3>   fSigUV;

	std::vector<TVector3>   fAxis1;
	std::vector<TVector3>   fAxis2;
	std::vector<TVector3>   fAxis3;

	std::vector<GFDetPlane> fTrackPlane;
	std::vector<GFDetPlane> fMidPlane;

	std::vector<Int_t>  	fClusterSize;
	std::vector<Double_t>   fClusterSizeShare;
	std::vector<Int_t> 		f2DClusterSize;
	std::vector<TVector3> 	fFullClusterSize;
	std::vector<TVector3> 	fFullClusterSizeTrack;
	std::vector<TVector3> 	fFullClusterSizeAxis;

	std::vector<Double_t>   fClusterAmp;


public:
	ClassDef(Particle,1)

};

#endif