/******************************************************

Analysis Task created by M.Steinen steinen@kph.uni-mainz.de
Analysis of Gamma Simulation with hypGe detectors
*******************************************************/ 


#ifndef PNDHYPGEGAMMAANA_H
#define PNDHYPGEGAMMAANA_H 1

#include <FairTask.h>
#include <TH1D.h>
#include <TClonesArray.h>
#include <TRandom.h>
#include <fstream>
#include "PndHypGeGammaAnaStorage.h"
#include <TSpectrum.h>
#include <TF1.h>


class PndHypGeGammaAna: public FairTask
{
	public:
		PndHypGeGammaAna();
		PndHypGeGammaAna(TString TxtFileName_Ext, Double_t GammaEnergy_Ext,Int_t nEvents = 0,Int_t nPeaks_Ext=1,Int_t PeakToLook_Ext = 1);
		~PndHypGeGammaAna();
	
		virtual InitStatus Init();
		virtual void Exec(Option_t* opt);
		void Finish();//TString cat);
		
		
		Double_t Compton(Double_t E, Double_t Th);
		Double_t invCompton(Double_t E, Double_t Eprime);
		
		
	protected:
	
		TString 									TxtFileName;
		ofstream									TxtFile;
		
		Double_t 									Resolution;
		
		Double_t 									bufferEnergy;
		
		TH1D*											hGamEnergy;
		TH1D*											hNumberOfHits;
		
		
		TClonesArray*							fMcTr;
		TClonesArray*							fMc;
		TClonesArray*							fGe;

		Long_t										EvtCount;
		Long_t 										NumberOfEvents;
		Double_t 									GammaEnergy;
		Double_t									iHistoUpperLimit;
		TString										fName;
		TRandom*									Rng;
		
		PndHypGeGammaAnaStorage*	fStorage;
		
		//Analysis
		Int_t 										nPeaks;
		Int_t 										PeakToLook;
		Double_t 									PeakX;
		Double_t 									PeakY; 
		Double_t 									DPeakY; 
		Double_t 									SumPeak;
		Double_t 									DSumPeak;
		Double_t									SumCompton;
		Double_t									DSumCompton;
		Double_t 									PeakToCompton;
		Double_t 									DPeakToCompton;
		
		TF1*											GausBG;
		
	ClassDef(PndHypGeGammaAna,1);
};

#endif /* PNDHYPGEGAMMAANA_H */