/**
 * \file CbmAnaConversionTomography.h
 *
 * \author Sascha Reinecke <reinecke@uni-wuppertal.de>
 * \date 2014
 * Using the true-MC data to do some tomography studies,
 * via conversion of gammas (gamma -> e+ e-)
 * somewhere in the detector material and plotting
 * the point of conversion.
 **/
 
 
#ifndef CBM_ANA_CONVERSION_TOMOGRAPHY
#define CBM_ANA_CONVERSION_TOMOGRAPHY

// included from ROOT
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TStopwatch.h>
#include <TClonesArray.h>

// included from CbmRoot
#include "CbmMCTrack.h"




using namespace std;


class CbmAnaConversionTomography
{

public:
	CbmAnaConversionTomography();
	virtual ~CbmAnaConversionTomography();

	void Init();
	void InitHistos();
	void Finish();

	void DoTomography();

	void TomographyMC(CbmMCTrack* mctrack);
	void TomographyReco(CbmMCTrack* mctrack);







private:
	TClonesArray* fMcTracks;

	vector<TH1*> fHistoList_tomography;	// list of all histograms


	TH1D * fhGammaZ;
	TH3D * fhTomography;
	TH2D * fhTomography_XZ;
	TH2D * fhTomography_YZ;
	TH2D * fhTomography_uptoRICH;
	TH2D * fhTomography_RICH_complete;
	TH2D * fhTomography_RICH_beampipe;
	TH2D * fhTomography_STS_end;
	TH2D * fhTomography_STS_lastStation;
	TH2D * fhTomography_RICH_frontplate;
	TH2D * fhTomography_RICH_backplate;

	TH1D * fhConversion;
	TH1D * fhConversion_inSTS;
	TH1D * fhConversion_prob;
	TH1D * fhConversion_energy;
	TH1D * fhConversion_p;


	TH3D * fhTomography_reco;
	TH2D * fhTomography_reco_XZ;
	TH2D * fhTomography_reco_YZ;


	// timer
	TStopwatch timer;
	Double_t fTime;

	CbmAnaConversionTomography(const CbmAnaConversionTomography&);
	CbmAnaConversionTomography operator=(const CbmAnaConversionTomography&);
     
	ClassDef(CbmAnaConversionTomography,1)
};

#endif