void ES1_4560()
{

	TStopwatch sw;
	sw.Start();

	BeAnalysis ana;

	ana.SetCMAngularRange(45, 60);

	ana.SetFigures("figures/4560/", ".eps");

	//	ana.SetChainsToDraw(0, 0, 0, 0, 1, 1);
	ana.SetChainsToDraw();

	Bool_t epsilonT[6] = {0,0,1,0,0};
	Bool_t cosThetaT[6] = {0,0,1,0,0};
	Bool_t epsilonY[6] = {0,0,1,0,0};
	Bool_t cosThetaY[6] = {0,0,1,0,0};
	Bool_t thetaAT[6] = {0,0,1,0,0};

	ana.SetSpectra();
	ana.SetEpsilonTintervals(epsilonT);
	ana.SetEpsilonTintervals();
	ana.SetCosThetaTkIntervals(cosThetaT);
	ana.SetCosThetaTkIntervals();
	ana.SetEpsilonYintervals(epsilonY);
	ana.SetEpsilonYintervals();
	ana.SetCosThetaYkIntervals(cosThetaY);
	ana.SetCosThetaYkIntervals();
	ana.SetThetaATintervals(thetaAT);
	ana.SetThetaATintervals();



	ana.SetExpChain("../../../be/rootdata/correlations/v5_exp/Be.", 0, 40);

	const Long64_t eMaxEvents = 5000000;
	Long64_t eEvents[6] = {3150000, 3150000, 3150000, 3150000, 3150000, 3150000};
	ana.SetNoExpEvents(eEvents);

	const Long64_t sMaxEvents = 3000000;
	Long64_t sEvents[6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents};
	ana.SetNoSimEvents(sEvents);

	Int_t lowSimFile[6] = {0, 0, 0, 0, 0, 0};
	Int_t upSimFile[6] = {15, 15, 15, 15, 15, 15};
	ana.SetNoSimFiles(lowSimFile, upSimFile);

	ana.SetSimChains("../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_0_35-70_",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_0_35-70_",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_180_35-70_",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_180_35-70_",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_90_35-70_",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_90_35-70_");

	ana.SetInputChains("../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_0_35-70_000.root",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_0_35-70_000.root",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_180_35-70_000.root",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_180_35-70_000.root",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_al_90_35-70_000.root",
					"../../../be/rootdata/correlations/v5_4560/Sim_mix_br_noal_90_35-70_000.root");


	ana.OpenExpChain();
	ana.OpenSimChains();


	ana.SetRatiosGStoEX("sRatio>0.05 && sRatio<0.095",
			"sRatio>0.05 && sRatio<0.095",
			"sRatio>0.05 && sRatio<0.095",
			"sRatio>0.05 && sRatio<0.095",
			"sRatio>0.05 && sRatio<0.095",
			"sRatio>0.05 && sRatio<0.095");

//	ana.SetRatiosGStoEX("sRatio>0.051 && sRatio<0.101",
//			"sRatio>0.050 && sRatio<0.100",
//			"sRatio>0.047 && sRatio<0.097",
//			"sRatio>0.048 && sRatio<0.098",
//			"sRatio>0.051 && sRatio<0.101",
//			"sRatio>0.051 && sRatio<0.101");
	ana.SetCuts();

	Long64_t eEventsECuts[5][6] = {2900000, 2900000, 2900000, 2900000, 2900000, 2900000,
										3400000, 3400000, 3400000, 3500000, 3400000, 3400000,
										2400000, 2400000, 2400000, 2400000, 2400000, 2400000,
										2900000, 3000000, 2900000, 3000000, 2900000, 3000000,
										3200000, 3300000, 3300000, 3400000, 3300000, 3400000};
	ana.ExpEventsECuts(eEventsECuts);

	Long64_t sEventsECuts[5][6] = {sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents,
										sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents, sMaxEvents};
	ana.SimEventsECuts(sEventsECuts);

	ana.SetRangeProportion(1.1, 1);

	///////////////////////////////////////////////////////////////////////////////////

	ana.Spectra();

	ana.EpsilonT();

	sw.Print();
	sw.Continue();

	ana.CosThetaTk();

	sw.Print();
	sw.Continue();

	ana.EpsilonY();

	sw.Print();
	sw.Continue();

	ana.CosThetaYk();

	sw.Print();
	sw.Continue();

	ana.ThetaAT();

	sw.Print();
	sw.Continue();

	Info("ES1_4560.cxx", "Finished.");

	return;

}