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


/** PNDMCTESTLAMBDALAMBDABAR.h
 *@author S.Esch <s.esch@fz-juelich.de>
 **
 ** Analysis for Lambda Lambda Bar Decay
 **/

#ifndef PNDMCTESTLAMBDALAMBDABAR_H
#define PNDMCTESTLAMBDALAMBDABAR_H

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

#include <vector>
#include <map>

class TClonesArray;

class PndMCTestLambdaLambdaBar : public FairTask
{
 public:

  /** Default constructor **/
	PndMCTestLambdaLambdaBar();

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


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


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

  virtual void Finish();

 private:

  int eventcounter;

   PndMCMatch* fMCMatch;
   TClonesArray* fMCPoint;
   TClonesArray* fMCTrack;
   TClonesArray* fEventHeader;

   TH1F* fMeanLifeLambdaCMS;
   TH1F* fMeanLifeLambdaBarCMS;
   TH1F* fMeanLifeLambdaLAB;
   TH1F* fMeanLifeLambdaBarLAB;

   TH2F* fPositionBothLambdasDecayed;
   TH2F* fPositionDecayVertexLambdaLambdaBar2D;
   TH2F* fPositionDecayVertexLambda2D;
   TH2F* fPositionDecayVertexLambdaBar2D;

   TH2F* fArmenterosPodolanskiPlot;
   TH1F* fAngularDistributionLambdaCMSCosTheta;
   TH1F* fAngularDistributionLambdabarCMSCosTheta;
   TH1F* fAngularDistributionProtonCMSCosTheta;
   TH1F* fAngularDistributionAntiprotonCMSCosTheta;
   TH1F* fAngularDistributionPionminusCMSCosTheta;
   TH1F* fAngularDistributionPionplusCMSCosTheta;

   TH1F* fAngularDistributionLambdaCMSTheta;
   TH1F* fAngularDistributionLambdabarCMSTheta;
   TH1F* fAngularDistributionProtonCMSTheta;
   TH1F* fAngularDistributionAntiprotonCMSTheta;
   TH1F* fAngularDistributionPionminusCMSTheta;
   TH1F* fAngularDistributionPionplusCMSTheta;

   TH1F* fAngularDistributionLambdaLABCosTheta;
   TH1F* fAngularDistributionLambdabarLABCosTheta;
   TH1F* fAngularDistributionProtonLABCosTheta;
   TH1F* fAngularDistributionAntiprotonLABCosTheta;
   TH1F* fAngularDistributionPionminusLABCosTheta;
   TH1F* fAngularDistributionPionplusLABCosTheta;

   TH1F* fAngularDistributionLambdaLABTheta;
   TH1F* fAngularDistributionLambdabarLABTheta;
   TH1F* fAngularDistributionProtonLABTheta;
   TH1F* fAngularDistributionAntiprotonLABTheta;
   TH1F* fAngularDistributionPionminusLABTheta;
   TH1F* fAngularDistributionPionplusLABTheta;

   TH2F* fPzVsPtLambda;
   TH2F* fPzVsPtLambdabar;

   TH2F* fDetectorCoverageMVDMCPointsProton;
   TH2F* fDetectorCoverageMVDMCPointsAntiproton;
   TH2F* fDetectorCoverageMVDMCPointsPionplus;
   TH2F* fDetectorCoverageMVDMCPointsPionminus;

   TH2F* fPositionMvdMcPoint;

   TH1I* fMCTracksPerEvent;

   FairLogger* fLogger;

   PndGeoHandling* fGeoH;

  void Register();

  void Reset();

  void ProduceHits();

  ClassDef(PndMCTestLambdaLambdaBar,1);

};

#endif