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

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 <TH2D.h>
#include <TClonesArray.h>
#include <TRandom.h>
#include <fstream>
#include <set>
//#include "PndHypGeGammaAnaStorage.h"
#include "PndHypGePoint.h"
#include "../../pnddata/HypData/PndHypPoint.h"
#include "PndMCTrack.h"
#include <TSpectrum.h>
#include <TF1.h>
#include "TGeoManager.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();
		void SetParContainers();
		void SetTarget(Bool_t  TargetAvailable = 1);
		void SetOmegaQuadrupolMode(Bool_t useOmegaQuadrupolMode_ext = 1);
		void SetPeakWidtchStrechFactor(Double_t PeakWidtchStrechFactor_ext = 1);

	protected:
	
		Double_t 									Compton(Double_t E, Double_t Th);
		Double_t									invCompton(Double_t E, Double_t Eprime);
		Double_t 									SmearHit(Double_t Energy);
		Double_t 									CalculatePeakWidth(Double_t Energy);
		void 											FitSpectrum();
		void 											FittingOutput();
		void 											FindPeaks();
		void 											GetPeakContent();
		void 											ActivateComptonCalculations(Bool_t isActivated = 1);
		void 											CalculateCompton();
		void 											WriteHistogramsToFile();
		void 											HistogramCosmetics(TH1D* histo, TString XTitle, Double_t XTitleSize, Double_t XTitleOffset, Double_t XLabelSize, TString YTitle, Double_t YTitleSize, Double_t YTitleOffset, Double_t YLabelSize);
		void 											HistogramCosmetics2D(TH2D* histo, TString XTitle, Double_t XTitleSize, Double_t XTitleOffset, Double_t XLabelSize, TString YTitle, Double_t YTitleSize, Double_t YTitleOffset, Double_t YLabelSize);

		bool											InTargetAbsorbedGammaWouldHitGermanium();
		TString 									TxtFileName;
		ofstream									TxtFile;
		
		Bool_t 										useOmegaQuadrupolMode;

		Double_t 									ResConstPar;
		Double_t 									ResGradientPar;
		Double_t									PeakWidtchStrechFactor;
		Double_t									PeakWidth;

		Double_t 									bufferEnergy;
		
		TH1D*											hGamEnergy;
		TH1D*											hGamEnergyTargetBackground;
		TH1D*											hNumberOfHits;
		TH1D*											hThetaCheck;
		TH1D*											hThetaCheckPrimary;
		TH1D*											hAbsorption;
		TH2D*											h2AbsorptionDistanceAngle;
		
		TClonesArray*							fMcTr;
		TClonesArray*							fMc;
		TClonesArray*							fGe;
		TClonesArray*							fGeAl;
		TClonesArray*							fSilicon;
		TClonesArray*							fTargetOther;

		PndHypGePoint 						*fHitGe;;
		PndHypPoint 							*fHitSi;
		PndHypPoint 							*fHitTargetOther;
		PndMCTrack 								*fMCGam;
		PndMCTrack 								*fMCGamPrimary;

		Long_t										EvtCount;
		Long_t 										NumberOfEvents;
		Double_t 									GammaEnergy;
		Double_t									iHistoUpperLimit;
		TString										fName;
		TRandom*									Rng;
		
		//PndHypGeGammaAnaStorage*	fStorage;
		
		FairRootManager* 					ioman;
		// Exec
		std::set<int> 						SetOfCrystalHit;
		TVector3 									fStartVertex;
		TString 									fVertexVolumeName;

		//Analysis
		Int_t 										nPeaks;
		Int_t 										PeakToLook;
		Double_t 									PeakX;
		Double_t 									PeakY; 
		Double_t 									DPeakY; 
		Double_t 									SumPeak;
		Double_t 									DSumPeak;
		Bool_t										DoComptonCalculations;
		Double_t 									lowComptonAngle;
		Double_t 									highComptonAngle;
		Double_t									SumCompton;
		Double_t									DSumCompton;
		Double_t 									PeakToCompton;
		Double_t 									DPeakToCompton;
		
		TF1*											GausBG;
		
		Float_t 									*xpeaks;	//!
		Float_t 									*ypeaks;	//!

		Bool_t										fTargetAvailable;



	ClassDef(PndHypGeGammaAna,1);
};

#endif /* PNDHYPGEGAMMAANA_H */