#ifndef __TPC_DEFINITION_H__
#define __TPC_DEFINITION_H__

#include "TApplication.h"

class TPC_Definition {
    int nbinsZ, nbinsRR, nbinsR;
	int LCE_min, LCE_max;
	double chamber_minZ, chamber_maxZ;
	double chamber_minRR, chamber_maxRR, chamber_maxR;
	double LXe_minZ, LXe_maxZ;
	double GXe_minZ, GXe_maxZ;
	double LXe_minRR, LXe_maxRR, LXe_maxR;
	double QE_top, QE_bottom;
	double PMTs_top, PMTs_bottom;
	
  public:
	TPC_Definition ();
	TPC_Definition (int,int,int);
    void TPC_Xe1T ();
	void TPC_MS ();
	void Init (TNamed*);
	void Set_Bins (int,int,int);
	void Set_LCE_max (int);
    int Get_nbinsZ (void) {return nbinsZ;}
	int Get_nbinsRR (void) {return nbinsRR;}
	int Get_nbinsR (void) {return nbinsR;}
	int Get_LCE_min (void) {return LCE_min;}
	int Get_LCE_max (void) {return LCE_max;}
	double Get_chamber_minZ (void) {return chamber_minZ;}
	double Get_chamber_maxZ (void) {return chamber_maxZ;}
	double Get_chamber_minRR (void) {return chamber_minRR;}
	double Get_chamber_maxRR (void) {return chamber_maxRR;}
	double Get_chamber_maxR (void) {return chamber_maxR;}
	double Get_LXe_minZ (void) {return LXe_minZ;}
	double Get_LXe_maxZ (void) {return LXe_maxZ;}
	double Get_GXe_minZ (void) {return GXe_minZ;}
	double Get_GXe_maxZ (void) {return GXe_maxZ;}
	double Get_LXe_minRR (void) {return LXe_minRR;}
	double Get_LXe_maxRR (void) {return LXe_maxRR;}
	double Get_LXe_maxR (void) {return LXe_maxR;}
	double Get_QE_top (void) {return QE_top;}
	double Get_QE_bottom (void) {return QE_bottom;}
	double Get_PMTs_top (void) {return PMTs_top;}
	double Get_PMTs_bottom (void) {return PMTs_bottom;}
};

TPC_Definition::TPC_Definition () {
	TPC_Xe1T ();
	nbinsZ  = 26;
	nbinsRR = 22;
	nbinsR  = 50;
	LCE_min = 0;
	LCE_max = 100;
}

TPC_Definition::TPC_Definition (int Z, int R, int RR) {
	TPC_Xe1T ();
	LCE_min = 0;
	LCE_max = 100;
	Set_Bins (Z, R, RR);
}

void TPC_Definition::TPC_Xe1T () {
  	chamber_minZ  = -140.;
	chamber_maxZ  = 6.;
	chamber_minRR = 0.;
	chamber_maxR  = 60.;
	chamber_maxRR = chamber_maxR*chamber_maxR;
	
	LXe_minZ  = -96.7;
	LXe_maxZ  = 0.;
	GXe_minZ  = 0.;
	GXe_maxZ  = 0.6;
	LXe_minRR = 0.;
	LXe_maxR  = 47.9;
	LXe_maxRR = LXe_maxR*LXe_maxR;
	
	QE_top = 0.314;
	QE_bottom = 0.366;
	PMTs_top = 127;
	PMTs_bottom = 121;
}

void TPC_Definition::TPC_MS () {
  	chamber_minZ  = -18.;
	chamber_maxZ  = 0.5;
	chamber_minRR = 0.;
	chamber_maxR  = 5.;
	chamber_maxRR = chamber_maxR*chamber_maxR;
	
	LXe_minZ  = -16.9;
	LXe_maxZ  = -0.2;
	GXe_minZ  = -0.2;
	GXe_maxZ  = 0.5;
	LXe_minRR = 0.;
	LXe_maxR  = 4.;
	LXe_maxRR = LXe_maxR*LXe_maxR;
	
	QE_top = 0.31;
	QE_bottom = 0.31;
	PMTs_top = 7;
	PMTs_bottom = 7;
}

void TPC_Definition::Init (TNamed *G4MCname) {
  	if ( strcmp(G4MCname->GetTitle(),"Xenon1t") == 0 ) {
		TPC_Xe1T();
		Set_LCE_max(50);
	}
	else if ( strcmp(G4MCname->GetTitle(),"TPC-MC") == 0 ){
		TPC_MS();
		Set_LCE_max(30);
	}
	else {
		cout << endl;
		cout << "x Error xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx" << endl;
		cout << "TPC name not known ..." << endl;
		cout << "xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx" << endl;
		cout << endl;
		gApplication->Terminate();
	}
}

void TPC_Definition::Set_Bins (int Z, int R, int RR) {
	if ( (Z>=200) || (R>=200) || (RR>=200) ) {
		cout << endl;
		cout << "x Error xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx" << endl;
		cout << "The amount of Bins can not be above 200x200:" << endl;
		cout << "-> Z: " << Z << " R: " << R << " RR: " << RR << endl;
		cout << "xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx" << endl;
		cout << endl;
		gApplication->Terminate();
	}
	nbinsZ  = Z;
	nbinsRR = RR;
	nbinsR  = R;
}

void TPC_Definition::Set_LCE_max (int LCE) {
	LCE_max = LCE;
}

#endif // __TPC_DEFINITION_H__