EpsilonT() implemented.

Be/BeAnalysis.cpp

@@ -22,6 +22,17 @@ BeAnalysis::BeAnalysis() : che(0), spectra(1) {
 	SetCuts();
 	SetChainsToDraw();
 
+	noIntervals = 5;
+	epsilonT = new Bool_t[noIntervals];
+	for (Int_t i = 0; i < noIntervals; i++) {
+		epsilonT[i] = 0;
+	}
+	epsilonT[0] = 1;
+	epsilonT[2] = 1;
+
+	kVerbose = 1;
+	kRangeProportionIn = 2.0;
+
 }
 
 BeAnalysis::~BeAnalysis() {
@@ -371,35 +382,13 @@ void BeAnalysis::Spectra() {
 
 void BeAnalysis::EpsilonT() {
 
-	const Int_t nointervals = 5;
-	const Bool_t epsilonT[nointervals] = {1,0,0,0,0};
-	const Bool_t kVerbose = 1;
-	const Float_t kRangeProportionIn = 2.0;
-
-	Long64_t eEventsECuts[5][6] = {3950000, 3950000, 3800000, 3850000, 3950000, 3900000,
-									eMaxEvents, eMaxEvents, 4100000, eMaxEvents, eMaxEvents, eMaxEvents,
-									3050000, 3150000, 3200000, 3200000, 3250000, 3350000,
-									3850000, 3950000, 3900000, 3900000, 4000000, 4000000,
-									4100000, 4100000, 4000000, 4100000, eMaxEvents, eMaxEvents};
-
-	Long64_t sEventsECuts[5][6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
-									sMaxEvents, 2630000, sMaxEvents, sMaxEvents, 2660000, 2660000,
-									sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
-									sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
-									sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents};
 
-	const Double_t eTRange[5][6] = {300, 300, 300, 300, 300, 300,
-									380, 380, 380, 380, 380, 380,
-									300, 300, 300, 300, 300, 300,
-									500, 500, 500, 500, 500, 500,
-									500, 500, 500, 500, 500, 500};
+	TCanvas *canEpsilonT[noIntervals];
+	for (Int_t i = 0; i<noIntervals; i++) {	//energy intervals (canvases)
 
-	TCanvas *canEpsilonT[nointervals];
-	for (Int_t i = 0; i<nointervals; i++) {	//energy intervals (canvases)
-
-		TH1F *hseT[nointervals][6];
-		TH1F *heeT[nointervals][6];
-		TH1F *hieT[nointervals][6];
+		TH1F *hseT[noIntervals][6];
+		TH1F *heeT[noIntervals][6];
+		TH1F *hieT[noIntervals][6];
 
 		if (!epsilonT[i]) continue;
 //		{
@@ -491,3 +480,83 @@ void BeAnalysis::EpsilonT() {
 	}//for i
 
 }
+
+void BeAnalysis::CosThetaTk() {
+
+}
+
+void BeAnalysis::EpsilonY() {
+}
+
+void BeAnalysis::CosThetaYk() {
+}
+
+void BeAnalysis::ThetaP1T() {
+}
+
+void BeAnalysis::ExpEventsECuts() {
+	eEventsECuts[0][0] = 3950000;
+	eEventsECuts[0][1] = 3950000;
+	eEventsECuts[0][2] = 3800000;
+	eEventsECuts[0][3] = 3850000;
+	eEventsECuts[0][4] = 3950000;
+	eEventsECuts[0][5] = 3900000;
+
+	eEventsECuts[1][0] = eMaxEvents;
+	eEventsECuts[1][1] = eMaxEvents;
+	eEventsECuts[1][2] = 4100000;
+	eEventsECuts[1][3] = eMaxEvents;
+	eEventsECuts[1][4] = eMaxEvents;
+	eEventsECuts[1][5] = eMaxEvents;
+
+	eEventsECuts[2][0] = 3050000;
+	eEventsECuts[2][1] = 3150000;
+	eEventsECuts[2][2] = 3200000;
+	eEventsECuts[2][3] = 3200000;
+	eEventsECuts[2][4] = 3250000;
+	eEventsECuts[2][5] = 3350000;
+
+	eEventsECuts[3][0] = 3850000;
+	eEventsECuts[3][1] = 3950000;
+	eEventsECuts[3][2] = 3900000;
+	eEventsECuts[3][3] = 3900000;
+	eEventsECuts[3][4] = 4000000;
+	eEventsECuts[3][5] = 4000000;
+
+	eEventsECuts[4][0]	= 4100000;
+	eEventsECuts[4][1]	= 4100000;
+	eEventsECuts[4][2]	= 4000000;
+	eEventsECuts[4][3]	= 4100000;
+	eEventsECuts[4][4]	= eMaxEvents;
+	eEventsECuts[4][5]	= eMaxEvents;
+}
+
+void BeAnalysis::SimEventsECuts() {
+
+	for (Int_t i = 0; i < 5; i++) {
+		for (Int_t j = 0; j < 6; j++) {
+			sEventsECuts[i][j] = sMaxEvents;
+		}
+	}
+
+	sEventsECuts[1][1] = 2630000;
+	sEventsECuts[1][4] = 2660000;
+	sEventsECuts[1][5] = 2660000;
+
+}
+
+void BeAnalysis::EpsilonTRange() {
+
+	for (Int_t i = 0; i < 6; i++) {
+		eTRange[0][i] = 300;
+		eTRange[1][i] = 380;
+		eTRange[2][i] = 300;
+		eTRange[3][i] = 500;
+		eTRange[4][i] = 500;
+	}
+
+
+}
+
+void BeAnalysis::ThetaAT() {
+}

Be/BeAnalysis.h

@@ -40,9 +40,18 @@ public:
 	void SetNoExpEvents();
 	void SetNoSimEvents();
 	void SetRatiosGStoEX();
+	void ExpEventsECuts();
+	void SimEventsECuts();
+	void EpsilonTRange();
 
 	void Spectra();
 	void EpsilonT();
+	void CosThetaTk();
+	void EpsilonY();
+	void CosThetaYk();
+	void ThetaP1T();
+	void ThetaAT();
+
 //
 private:
 	TChain *che;		//chain with experimental trees
@@ -56,9 +65,13 @@ private:
 	TCut crAngles;
 
 	Bool_t kChains[6];
+	Int_t noIntervals;
+	Bool_t *epsilonT;	//!
 
 	Bool_t kAutoRange;
 	Float_t kRangeProportion;
+	Bool_t kVerbose;
+	Float_t kRangeProportionIn;
 
 	Long64_t eMaxEvents;
 	Long64_t eEventsAl0;
@@ -86,6 +99,10 @@ private:
 	TCut sRatioNoAl90;
 	TCut sRatio[6];
 
+	Long64_t eEventsECuts[5][6];
+	Long64_t sEventsECuts[5][6];
+	Double_t eTRange[5][6];
+
 	//cuts
 	TCut cBe20;
 	TCut cBe3;

macros/BeCorrPRC/lib_test.cxx

@@ -11,6 +11,10 @@ void lib_test() {
 	ana.SetNoSimEvents();
 	ana.SetRatiosGStoEX();
 
-	ana.Spectra();
+//	ana.Spectra();
+
+	ana.ExpEventsECuts();
+	ana.SimEventsECuts();
+	ana.EpsilonTRange();
 	ana.EpsilonT();
 }