// -------------------------------------------------------------------------
// -----                    PNDMCTESTLAMBDALAMBDABAR header file               -----
// -----                    Created 19/03/13  by S.Esch                 -----
// -------------------------------------------------------------------------


/** PndAnaPidLambdaLambdaBar.h
 *@author S.Esch <s.esch@fz-juelich.de>
 **
 **
 **/


#ifndef PndAnaPidLambdaLambdaBar_H
#define PndAnaPidLambdaLambdaBar_H


// framework includes
#include "FairTask.h"
#include "PndMCMatch.h"
#include "PndDetectorList.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "PndGeoHandling.h"
#include "THStack.h"
#include "PndPidCandidate.h"

#include <vector>
#include <map>

class TClonesArray;

class PndAnaPidLambdaLambdaBar : public FairTask
{
 public:

  /** Default constructor **/
	PndAnaPidLambdaLambdaBar();

  /** Destructor **/
  virtual ~PndAnaPidLambdaLambdaBar();


  /** Virtual method Init **/
  virtual void SetParContainers();
  virtual InitStatus Init();


  /** Virtual method Exec **/
  virtual void Exec(Option_t* opt);

  virtual void Finish();

  void SetInputBranchName(TString inputbranchname){fInputBranchName = inputbranchname;};

  void SetNumberOfMVDHits(int numberofmvdhits){fNumberOfMVDHits = numberofmvdhits;};
  void SetNumberOfTotalHits(int numberoftotalhits){fNumberOfTotalHits = numberoftotalhits;};


 private:

void MVDStripCounter(FairMultiLinkedData MVDStriplinks, int &mvdstrip, int &mvdstriplambda);
int FillAnalysisHistoSingleDetectors(TH2I* SingleParticleHisto, TH2I* TotalHisto, TH2I* SingleParticleSelected, TH2I* TotalSelected, TH2I* SingleParticleRejected, TH2I* TotalRejected, int fts, int stt, int gem, int mvd);
int FillAnalysisHistoDoubleDetectors(TH2I* SingleParticleHisto, TH2I* TotalHisto, TH2I* SingleParticleSelected, TH2I* TotalSelected, TH2I* SingleParticleRejected, TH2I* TotalRejected, int fts, int stt, int gem, int mvd);
int FillAnalysisHistoTripleDetectors(TH2I* SingleParticleHisto, TH2I* TotalHisto, TH2I* SingleParticleSelected, TH2I* TotalSelected, TH2I* SingleParticleRejected, TH2I* TotalRejected, int fts, int stt, int gem, int mvd);
int FillAnalysisHistoQuatroDetectors(TH2I* SingleParticleHisto, TH2I* TotalHisto, TH2I* SingleParticleSelected, TH2I* TotalSelected, TH2I* SingleParticleRejected, TH2I* TotalRejected, int fts, int stt, int gem, int mvd);
bool CutCriteria(int fts, int stt, int gem, int mvd);
void SetHitAnalysisAxis(TH1* histo);
void FillAnalysisTracksPerEventSingleTrack(int proton, int antiproton, int pionminus, int pionplus, TH2I* histo);
void SetHitAnalysisEventAxis(TH1* histo);
int GEMCounter(FairMultiLinkedData GEMlinks);

   TString fInputBranchName;

   int fNumberOfMVDHits;
   int fNumberOfTotalHits;

   int eventcounter;

   PndMCMatch* fMCMatch;
   TClonesArray* fStripHit;
   TClonesArray* fPixelHit;
   TClonesArray* fMCPoint;
   TClonesArray* fMCTrack;
   TClonesArray* fEventHeader;
   TClonesArray* fSTTHit;
   TClonesArray* fGEMHit;
   TClonesArray* fGEMPoint;
   TClonesArray* fFTSHit;
   TClonesArray* fPID;
   int fErrorCounter;

   TH2I* fAnaPid_HitAnalysisTracksPerEvent2D;
   PndGeoHandling* fGeoH;

  void Register();

  void Reset();

  void ProduceHits();

  ClassDef(PndAnaPidLambdaLambdaBar,1);

};

#endif