//
// Cut definitions for E06007 test analysis
//

{
//cout << "at the top 1" << endl;

TCut OneTrackOnly =  "L.tr.n==1&&R.tr.n==1" ;

TCut xR   = "abs(R.tr.r_x)<0.4" ;
TCut yR   = "R.tr.r_y>-0.005 && R.tr.r_y<0.03" ;
TCut thR   = "R.tr.r_th>-0.015 && R.tr.r_th<0.015" ;
TCut phR   = "R.tr.r_ph>0.0 && R.tr.r_ph<0.015" ;

TCut xL   = "abs(L.tr.r_x)<0.4" ;
TCut yL   = "L.tr.r_y>-0.025 && L.tr.r_y<0.01" ;
TCut thL   = "L.tr.r_th>-0.015 && L.tr.r_th<0.02" ;
TCut phL   = "L.tr.r_ph>-0.01 && L.tr.r_ph<0.0" ;

TCut centfpR = xR && yR && thR && phR ;
TCut centfpL = xL && yL && thL && phL ;
TCut centfpLR = centfpR && centfpL ;

// TCut YtargR  = "abs(R.gold.y-0.00416)<0.006" ;
TCut YtargR   = "abs(R.gold.y-0.004)<0.006" ;
TCut YtargL   = "abs(L.gold.y-0.003)<0.005" ;

TCut deltaL =  "abs(L.gold.dp)<0.05" ;
TCut deltaR =  "abs(R.gold.dp)<0.05" ;
TCut smalldeltaL = "abs(L.gold.dp)<0.045" ;
TCut smalldeltaR = "abs(R.gold.dp)<0.045" ;

TCut anglesL = "(L.gold.ph>-0.30 && L.gold.ph<0.030) && (L.gold.th>-0.06 && L.gold.th<0.060)" ;
TCut anglesR = "(R.gold.ph>-0.030 && R.gold.ph<0.030) && (R.gold.th>-0.060 && R.gold.th<0.060)";

//the best result

//TCut smallangL = "(L.gold.ph>-0.012 && L.gold.ph<0.012) && (L.gold.th>-0.024 && L.gold.th<0.004)";
//TCut smallangR = "(R.gold.ph>-0.012&& R.gold.ph<0.012) && (R.gold.th>-0.024 && R.gold.th<0.004)";

//to improve

//TCut smallangL = "(L.gold.ph>-0.015 && L.gold.ph<0.015 && L.gold.th>-0.028 && L.gold.th<0.004)";
//TCut smallangR = "(R.gold.ph>-0.015&& R.gold.ph<0.015 && R.gold.th>-0.028 && R.gold.th<0.004)";
TCut smallangL = "(L.gold.ph>-0.03 && L.gold.ph<0.03 && L.gold.th>-0.032 && L.gold.th<0.012)";
TCut smallangR = "(R.gold.ph>-0.026&& R.gold.ph<0.026 && R.gold.th>-0.032 && R.gold.th<0.016)";

//TCut normalangL = "(L.gold.ph>-0.03 && L.gold.ph<0.03 && L.gold.th>-0.07 && L.gold.th<0.07)";
//TCut normalangR = "(R.gold.ph>-0.03 && R.gold.ph<0.03 && R.gold.th>-0.06 && R.gold.th<0.08)";

//TCut normalangL = "(L.gold.ph>-0.01 && L.gold.ph<0.01) && (L.gold.th>-0.02 && L.gold.th<0.02)";
//TCut normalangR = "(R.gold.ph>-0.01&& R.gold.ph<0.01 && R.gold.th>-0.02 && R.gold.th<0.02)";

TCut normalangL = "(L.gold.ph>-0.035 && L.gold.ph<0.035) && (L.gold.th>-0.05 && L.gold.th<0.05)";
TCut normalangR = "(R.gold.ph>-0.035&& R.gold.ph<0.035 && R.gold.th>-0.05 && R.gold.th<0.05)";

TCut normalAcc = normalangL && normalangR && smalldeltaL && smalldeltaR;
TCut smallAcc = smallangL && smallangR && smalldeltaL && smalldeltaR;
TCut coinc = "abs(CT_LvR.ct_1by2*1000000000-200)<6 && abs(L.gold.y-R.gold.y+0.0058)<1000.10";
TCut ReactPt_L_y = "ReactPt_L.y>-0.0015&&ReactPt_L.y<0.0015"; 
TCut ReactPt_L_x = "ReactPt_L.x>0.00&&ReactPt_L.x<0.009";
TCut ReactPt_L =  ReactPt_L_x&&ReactPt_L_y   ;
TCut emiss = "abs(SKxceb.emiss*1000.-20.)<30";

TCut GoodPathlL = "abs(L.tr.pathl-24)<1." ;
TCut GoodPathlR = "abs(R.tr.pathl-24)<1." ;

TCut GoodLeft  = GoodPathlL && deltaL && anglesL && YtargL;
TCut GoodRight = GoodPathlR && deltaR && anglesR && YtargR ;

TCut GoodLR = GoodLeft && GoodRight ;

TCut T1 = "(D.evtypebits&0x02)==0x02" ;
TCut T3 = "(D.evtypebits&0x08)==0x08" ;
TCut T5 = "(D.evtypebits&0x20)==0x20" ;
TCut T6 = "(D.evtypebits&0x40)==0x40" ;
TCut T7 = "(D.evtypebits&0x80)==0x80" ;

TCut T1good = T1 && GoodRight ;
TCut T3good = T3 && GoodLeft ;
TCut T5good = T5 && GoodLR ;
TCut T6good = T6 && GoodLeft ;
TCut T7good = T7 && GoodLR ;

TCut dp1 = "abs(L.tr.tg_dp)<0.01" ;

TCut leftSing = (T3good||(T5&&GoodLeft));
TCut TrigType = T5good ;

TCut par_ct  = "abs(CT_S2LvR.ct_1by2*1.e9-258)<7.5" ;
TCut par_ctBck = "abs(CT_S2LvR.ct_1by2*1.e9-170)<14." ;
TCut fullct  = par_ct && GoodLR ;

// TCut u1oneclust = "L.vdc.u1.nclust==1 && L.vdc.u1.nhit>=5" ;
// TCut u2oneclust = "L.vdc.u2.nclust==1 && L.vdc.u2.nhit>=5 " ;
// TCut v1oneclust = "L.vdc.v1.nclust==1 && L.vdc.v1.nhit>=5" ;
// TCut v2oneclust = "L.vdc.v2.nclust==1 && L.vdc.v2.nhit>=5 " ;
// TCut vdconeclust = u1oneclust &&  u2oneclust && v1oneclust && v2oneclust ;

TCut u1oneclust = "L.vdc.u1.nclust==1 ";
TCut u2oneclust = "L.vdc.u2.nclust==1 " ;
TCut v1oneclust = "L.vdc.v1.nclust==1 " ;
TCut v2oneclust = "L.vdc.v2.nclust==1 " ;
TCut vdconeclust = u1oneclust &&  u2oneclust && v1oneclust && v2oneclust ;


// TCut par_ctPi = "abs(CT_S2LvR.ct_1by2*1.e9-170.2)<1." ; 
// TCut par_ctK = "abs(CT_S2LvR.ct_1by2*1.e9-172.8)<1." ;
// TCut par_ctP = "abs(CT_S2LvR.ct_1by2*1.e9-180.4)<1." ;

// TCut par_ctPi = "abs(CT_S2LvR.ct_1by2*1.e9-154.3)<1.2" ; 
// TCut par_ctK = "abs(CT_S2LvR.ct_1by2*1.e9-156.8)<1.2" ;
// TCut par_ctP = "abs(CT_S2LvR.ct_1by2*1.e9-164.4)<1.2" ;


TCut par_ctPi = "abs(CT_S2LvR.ct_1by2*1.e9-138.5)<1.4" ; 
TCut par_ctK = "abs(CT_S2LvR.ct_1by2*1.e9-141.3)<1.4" ;
TCut par_ctP = "abs(CT_S2LvR.ct_1by2*1.e9-148.6)<1.4" ;

TCut ctPi = par_ctPi && GoodLR  ;
TCut ctK  = par_ctK  && GoodLR  ;
TCut ctP  = par_ctP  && GoodLR  ;

TCut ctBck = par_ctBck && GoodLR  ;

TCut coinPi  = ctPi  ;
TCut coinK   = ctK  ;
TCut coinP   = ctP  ;
TCut coinBck = ctBck ;

TCut Lcer_e  = "L.cer.asum_c>200" ;
TCut Lcer_pi = "L.cer.asum_c<100" ;

TCut cutPi  =  fullct && coinPi ;
TCut cutK   =  fullct && coinK ;
TCut cutP   =  fullct && coinP;

TCut GoodFrancesco = YtargR && YtargL && anglesR && anglesL && deltaL && deltaR && par_ct;
TCut GoodSmallFrancesco = YtargR && YtargL && smallangR && smallangL && smalldeltaL && smalldeltaR && par_ct;
 /////////////////////////////////////////////////////////


	//TCut BeamAxmin = "urb.BPMA.x>-0.0025";
	//TCut BeamAxmax = "urb.BPMA.x<-0.00207";

		//TCut BeamAxmax = "urb.BPMA.x<-0.00212";
		//TCut BeamAxmin = "urb.BPMA.x>-0.00222";

		TCut BeamAxmax = "urb.BPMA.x<-0.00197";
		TCut BeamAxmin = "urb.BPMA.x>-0.00237";	

TCut BeamAxProve = "(urb.BPMA.x > -0.00202 || urb.BPMA.x < -0.00232)";
TCut AxProvem = "(urb.BPMA.x < -0.00220505 && urb.BPMA.x > -0.00235201)";

TCut BeamAx = BeamAxmin && BeamAxmax;

//TCut BeamAymin = "urb.BPMA.y>-0.00005";
//TCut BeamAymax = "urb.BPMA.y<0.00035";

TCut BeamAymin = "urb.BPMA.y>-0.00005";
TCut BeamAymax = "urb.BPMA.y<0.000382";

TCut BeamAy = BeamAymin && BeamAymax;

TCut BeamA = BeamAx && BeamAy;

	TCut BeamBxmin = "urb.BPMB.x>0.0015";
	TCut BeamBxmax = "urb.BPMB.x<0.00175";

//TCut BeamBxmin = "urb.BPMB.x>0.00146";
//TCut BeamBxmax = "urb.BPMB.x<0.00177";



TCut BeamBx = BeamBxmin && BeamBxmax;


TCut BeamBymin = "urb.BPMB.y>0.";
TCut BeamBymax = "urb.BPMB.y<0.000358";

TCut BeamBy = BeamBymin && BeamBymax;

TCut BeamB = BeamBx && BeamBy;

TCut BeamAB = BeamA && BeamB; 

//TCut Beamph = "(urb.ph > - 0.00001 &&  urb.ph < 0.000012)";

TCut Beamph = "(urb.ph > - 0.00004 &&  urb.ph < 0.00004)";

//TCut BeamC12 = "(IPM1C12.XPOS > - 0.04 && IPM1C12.XPOS < 0.04 && IPM1C12.XPOS !=0)";

TCut BeamC12 = "(IPM1C12.XPOS > - 0.04 && IPM1C12.XPOS < 0.04)";

TCut BeamC12Prove = "(IPM1C12.XPOS >  -0.015  && IPM1C12.XPOS < 0.015)";

TCut BeamC12tight = "(IPM1C12.XPOS > - 0.01 && IPM1C12.XPOS < 0.01 && IPM1C12.XPOS != 0.)";

TCut BeamC12medium = "((IPM1C12.XPOS > - 0.04 && IPM1C12.XPOS < - 0.015) ||(IPM1C12.XPOS > 0.015 && IPM1C12.XPOS < 0.040))";

TCut BeamC12N = "(IPM1C12.XPOS < -0.04.)";

TCut BeamC12Null = "(IPM1C12.XPOS  == 0)";

TCut BeamC12loose = "(IPM1C12.XPOS > -0.1 && IPM1C12.XPOS < 0.1 && IPM1C12.XPOS!=0.)";
TCut BeamC12loosenot = "(IPM1C12.XPOS < -0.1 || IPM1C12.XPOS > 0.1 || IPM1C12.XPOS == 0.)";

TCut BeamC12P = "(IPM1C12.XPOS > 0.04)";

TCut Beamhalla_MeV = "(halla_MeV > 3774. && halla_MeV <3778.)";
TCut Beamhalla_MeV_Normal = "(halla_MeV > 3000. && halla_MeV <4000.)";
TCut Beamhalla_MeVP = "(halla_MeV > -2000.)";

TCut Slia = "slia_MV_AO_SPREADS < 2.5 && slia_MV_AO_SPREADS > 1. ";
 

TCut RunCut = "(fEvtHdr.fRun>2579) && (fEvtHdr.fRun<2586)";
TCut RunCut1 = "(fEvtHdr.fRun>2715) && (fEvtHdr.fRun<2722)";
TCut RunCut2 = "(fEvtHdr.fRun>2630)";               

}