#ifndef PNDANALYSIS_H
#define PNDANALYSIS_H 1


#include "FairTask.h"
#include <map>
#include <string>

class TClonesArray;
class TObjectArray;
class TH1F;

class TPidChargedSelector;
class TPidNeutralSelector;
class TPidPlusSelector;
class TPidMinusSelector;
class TPidMassSelector;
class TPidSimpleKaonSelector;
class TPidSimplePionSelector;

class PndAnalysis : public FairTask
{

 public:
	typedef std::map<Int_t, Float_t> mapper;

  /** Default constructor **/  
  PndAnalysis();


  /** Destructor **/
  ~PndAnalysis();


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


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

  virtual void Finish();
  //void CreateStructure();

 protected:

  int evcount;

  

  // **** some basic selectors
  //
  //TPidChargedSelector *chargedSel;
  TPidNeutralSelector *neutralSel;
  TPidPlusSelector    *plusSel;
  TPidMinusSelector   *minusSel;
  
  // **** mass selectors for the resonances/composites
  //
  TPidMassSelector *phiMSel;
  TPidMassSelector *pi0MSel;
  TPidMassSelector *dsMSel;
  TPidSimpleKaonSelector *kSel;
  TPidSimplePionSelector *piSel;

  /**book all the histograms**/
  TH1F *phimass;
  TH1F *pi0mass;
  TH1F *dsmass; 
  TH1F *ds0mass;
  TH1F *ppmass;

  TH1F *nmult;

 private: 
   
  /** Input array of LheTrack **/
  TClonesArray* fChargedArray;
  TClonesArray* fNeutralArray;
 
 
  /** Geo file to use **/

	/** Get parameter containers **/
	virtual void SetParContainers();

 
  ClassDef(PndAnalysis,1);
  
};

#endif