#ifndef PndSolmazAnaTask_H
#define PndSolmazAnaTask_H 1


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

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

class RhoMassParticleSelector;
class PndAnalysis;
class RhoCandList;
class RhoTuple;


class PndSolmazAnaTask : public FairTask
{

 public:
	
	// ** Default constructor   
        PndSolmazAnaTask(double pbarmom, TString outname, TString pidalg, TString pidcrit);
	
	// ** Destructor 
	~PndSolmazAnaTask();	
	
	// ** Virtual method Init 
	virtual InitStatus Init();
	
	// ** Virtual method Exec 
	virtual void Exec(Option_t* opt);
	
	
	virtual void Finish();

 protected:
	
	
 private: 
	// *** event counter
	int fEvtCount;	
	
	// *** mass selector for the pi0
	RhoMassParticleSelector *pi0MassSel;
	
	// *** a method 
	int  SelectTruePid(PndAnalysis *ana, RhoCandList &l);
	
	// *** declare some histograms

	TH1F *hPi0Mult;
	TH1F *hGamMult;

	// *** declare some ntuples
	
	RhoTuple *nmc;
	RhoTuple *npi0;

	// *** the initial 4-vector
	TLorentzVector fIni;
	
	// *** the PndAnalysis object
	PndAnalysis *fAnalysis;
        TDatabasePDG  *fPdg;
        TFile             *fFile;
        TString           fOutName;
        TString           fPidAlg;
        TString           fPidCrit;
	
	// *** Get parameter containers
	virtual void SetParContainers();

	ClassDef(PndSolmazAnaTask,1);
  
};

#endif