cut_L_dp_old()
{
	
	/////////////////////////////////////////////////////////////////
	TString CutName = "m4";
	UInt_t run = 2528;
	
// 	TString CutName = "m2";
// 	UInt_t run = 2529;
	
// 	TString CutName = "0";
// 	UInt_t run = 2531;
	
// 	TString CutName = "p2";
// 	UInt_t run = 2532;
	
// 	TString CutName = "p4";
// 	UInt_t run = 2533;
	/////////////////////////////////////////////////////////////////

	
	TFile* fdata = new TFile(Form("ROOTfiles/e06010_det_L_%d.root",run));
	TTree* t=T;
	
	gSystem->Exec("cp -vf ./cutfiles/gcut_L_dp.root ./cutfiles/gcut_L_dp.root.old");
	TFile* fcut = new TFile("./cutfiles/gcut_L_dp.root","UPDATE");fcut->cd();
	
	TCanvas * c1 = new TCanvas("asdf","asdf",1900,1100);
	t->Draw(Form("L.tr.r_ph:L.tr.tg_dp>>hDpCut%d(500,-.07,.07,500,-.07,.07)",run)
			,"L.tr.n==1 && abs(L.tr.tg_th)<.1 && abs(L.tr.tg_ph)<.05","COLZ"); 
	c1->Update();
	
	cutg = (TCutG*) (gPad->WaitPrimitive("CUTG","CutG")); // making cut, store to CUTG
	c1->Update();
	
	cutg->SetName(Form("L_dp_%s",CutName.Data())); //
  // set axises' name
	cutg->SetVarX("L.tr.tg_dp");
	cutg->SetVarY("L.tr.r_ph");
	cutg->Write("",TObject::kOverwrite); // Overwrite old cut
	
	fcut->Write();
	
}

cut_L_dp()
{
	
	/////////////////////////////////////////////////////////////////
// 	TString CutName = "m4";
// 	UInt_t run = 2528;
	
// 	TString CutName = "m2";
// 	UInt_t run = 2529;
	
// 	TString CutName = "0";
// 	UInt_t run = 2531;
	
// 	TString CutName = "p2";
// 	UInt_t run = 2532;
	
	TString CutName = "p4";
	UInt_t run = 2533;
	/////////////////////////////////////////////////////////////////

	
	
	TFile* fdata = new TFile(Form("ROOTfiles/e06010_det_L_%d.root",run));
	TTree* t=T;
	
	gSystem->Exec("cp -vf ./cutfiles/gcut_L_dpfull.root ./cutfiles/gcut_L_dpfull.root.old2");
	TFile* fcut = new TFile("./cutfiles/gcut_L_dpfull.root","UPDATE");fcut->cd();

	
	TCanvas * c1 = new TCanvas("asdf","asdf",1900,1100);

	for(UInt_t FoilID = 1;FoilID<=6;FoilID++)
	{

		TString VertexCut = Form("fcut_L_%d",FoilID);
		cout<<"Working on Kine "<<CutName<<", Vertex Cut "<<VertexCut<<endl;
		
		TCutG *cVertexCut = (TCutG*)gROOT->FindObject(VertexCut);
		assert(cVertexCut);
	
		t->Draw(Form("L.tr.tg_th:sqrt(PriKine_C.W2)>>hDpCut%d_Foil%d(1000,11.15,11.25,1000,-.1,.1)",run,FoilID)
				,VertexCut+"&& L.tr.n==1 && abs(L.tr.tg_th)<.1 && abs(L.tr.tg_ph)<.05","COLZ");
		c1->Update();
		
		cutg = (TCutG*) (gPad->WaitPrimitive("CUTG","CutG")); // making cut, store to CUTG
		c1->Update();
		
		cutg->SetName(Form("L_dp_%s_%d",CutName.Data(),FoilID)); //
	  // set axises' name
		cutg->SetVarX("sqrt(PriKine_C.W2)");
		cutg->SetVarY("L.tr.tg_th");
		cutg->Write("",TObject::kOverwrite); // Overwrite old cut
	}
	
	fcut->Write();
	
}