#ifndef PndScrutJpsiEtacTask_H
#define PndScrutJpsiEtacTask_H 1



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

class TClonesArray;
class TH1F;
class TH2F;
class TDatabasePDG;

class RhoMassParticleSelector;
class RhoEnergyParticleSelector;
class PndAnalysis;
class RhoCandList;
class RhoCandidate;
class RhoTuple;
class PndRhoTupleQA;



class PndScrutJpsiEtacTask : public FairTask
{

 public:

	// ** Default constructor
  PndScrutJpsiEtacTask(double pbarmom, TString outname, int decayMode=0, int recoilMode=0, int dataMode=0, TString pidQuality="All", Double_t chiqCut=50);

	// ** Destructor
	~PndScrutJpsiEtacTask();

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

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

	virtual void Finish();

 protected:


 private:
	// *** event counter
	int fEvtCount;

	// *** the initial 4-vector
	TLorentzVector fIni;

	// *** the PndAnalysis object
	PndAnalysis   *fAnalysis;
	TDatabasePDG  *fPdg;
	TFile 		  *fFile;
	TString		  fOutName;
    Double_t      fChiqCut;

	// *******
	// ******* DECLARE THE STUFF YOU NEED
	// *******
        int dMode, rMode, fMode;
	TString fPIDquality;

	//RhoTuple *ntp1;
	RhoTuple *ntp2;
	RhoTuple *nmc;

	RhoMassParticleSelector *jpsiMassSel;
	RhoMassParticleSelector *etaMassSel;
	RhoMassParticleSelector *etacMassSel;
	RhoMassParticleSelector *pi0MassSel;
	RhoEnergyParticleSelector *gammaSel;
	RhoMassParticleSelector *KsMassSel;
	RhoMassParticleSelector *pbarpMassSel;

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


	ClassDef(PndScrutJpsiEtacTask,1);

};

#endif