{
	//script to be loaded in statesRatioFitting* scripts
	//

	//file with Be events only
	TCut cBe20 = "fBeIM>0 && fBeIM<20";
	TCut cBe3 = "fBeIM>0 && fBeIM<3";
	TCut cBeWork = "fBeIM>0 && fBeIM<10";
	TCut cBe0_14 = "fBeIM>0 && fBeIM<1.4";
//	TCut cBe0_14 = "fBeIM>0 && fBeIM<1.6";
//	TCut cBe0_14 = "fBeIM>0 && fBeIM<1.20";
	TCut cBe14_19 = "fBeIM>1.4 && fBeIM<1.9";
	TCut cBe19_25 = "fBeIM>1.9 && fBeIM<2.5";
	TCut cBe25_31 = "fBeIM>2.5 && fBeIM<3.1";
	TCut cBe31_37 = "fBeIM>3.1 && fBeIM<3.7";
	TCut cBeE[5] = {cBe0_14, cBe14_19, cBe19_25, cBe25_31, cBe31_37};

	TCut cEpsilonT = "fTpp/fBeIM<0.2";
//	TCut cEpsilonT = "fTpp/fBeIM<1.";
//	TCut cEpsilonY = "fTap/fBeIM<0.5";
	TCut cEpsilonY = "fTap/fBeIM>0.7";

	//raw files
	//energy cuts
	TCut crBe20 = "f6BeIM>0 && f6BeIM<20";
	TCut crBe3 = "f6BeIM>0 && f6BeIM<3";
	TCut crBeWork = "f6BeIM>0 && f6BeIM<10";
	TCut crBe0_14 = "f6BeIM>0 && f6BeIM<1.4";
//	TCut crBe0_14 = "f6BeIM>0 && f6BeIM<1.6";
//	TCut crBe0_14 = "f6BeIM>0.5 && f6BeIM<1.20";
	TCut crBe14_19 = "f6BeIM>1.4 && f6BeIM<1.9";
	TCut crBe19_25 = "f6BeIM>1.9 && f6BeIM<2.5";
	TCut crBe25_31 = "f6BeIM>2.5 && f6BeIM<3.1";
	TCut crBe31_37 = "f6BeIM>3.1 && f6BeIM<3.7";
	TCut crBeE[5] = {crBe0_14, crBe14_19, crBe19_25, crBe25_31, crBe31_37};
	//angular cuts


//	TCut crEpsilonT = "fTpp/f6BeIM>0.8";
	TCut crEpsilonT = "fTpp/f6BeIM<0.2";
//	TCut crEpsilonT = "fTpp/f6BeIM<1.";
//	TCut crEpsilonY = "fTap/f6BeIM<0.5";
	TCut crEpsilonY = "fTap/f6BeIM>0.7";

	//simulation input
	TCut ciBe0_14 = "E_IM>0 && E_IM<1.4";
//	TCut ciBe0_14 = "E_IM>0 && E_IM<1.6";
//	TCut ciBe0_14 = "E_IM>0.5 && E_IM<1.20";
	TCut ciBe14_19 = "E_IM>1.4 && E_IM<1.9";
	TCut ciBe19_25 = "E_IM>1.9 && E_IM<2.5";
	TCut ciBe25_31 = "E_IM>2.5 && E_IM<3.1";
	TCut ciBe31_37 = "E_IM>3.1 && E_IM<3.7";
	TCut ciBeE[5] = {ciBe0_14, ciBe14_19, ciBe19_25, ciBe25_31, ciBe31_37};

	TCut ciEpsilon = "sTpp/E_IM<0.2";
//	TCut ciEpsilon = "sTpp/E_IM<1.";
//	TCut ciEpsilonY = "sTap/E_IM<0.5";
	TCut ciEpsilonY = "sTap/E_IM>0.7";

	TCut cQ = "TMath::Abs(fQLiP)<10";
	TCut cProtons = "fP1Lab.fE-938.272<50 && fP2Lab.fE-938.272<50";
//	TCut cProtons = "fP1Lab.fE-938.272<34 && fP2Lab.fE-938.272<34 && fP1Lab.fE-938.272>1 && fP2Lab.fE-938.272>1";


}