{
        gROOT->Reset();
        c1 = new TCanvas("c1","",10,100,1000,500);
        TChain ch1("h1000");
	ch1.Add("pl_et_must_cible_kin_r115a168.root");
	ch1.Add("pl_et_must_cible_kin_r169a197.root");
	ch1.Add("pl_et_must_cible_kin_r198a228.root");
// Cut graphique en Energie-Temps dans les strips :
//        TCutG *cutg = new TCutG("CUTG",11);
//   cutg->SetVarX("Tmuca2");
//   cutg->SetVarY("Ep");
//   cutg->SetPoint(0,2026.7,699.809);
//   cutg->SetPoint(1,2022.76,481.77);
//   cutg->SetPoint(2,2001.77,295.47);
//   cutg->SetPoint(3,1980.79,172.651);
//   cutg->SetPoint(4,1966.36,106.411);
//   cutg->SetPoint(5,2029.32,107.791);
//   cutg->SetPoint(6,2051.62,241.65);
//   cutg->SetPoint(7,2064.73,353.43);
//   cutg->SetPoint(8,2076.54,507.99);
//   cutg->SetPoint(9,2075.23,634.949);
//   cutg->SetPoint(10,2026.7,699.809);
   TCutG *cutg = new TCutG("CUTBE",11);
   cutg->SetVarX("Tmuca2");
   cutg->SetVarY("Ep");
   cutg->SetTitle("Graph");
   cutg->SetFillColor(1);
   cutg->SetPoint(0,2044.66,687.925);
   cutg->SetPoint(1,2080,687.925);
   cutg->SetPoint(2,2064.58,371.464);
   cutg->SetPoint(3,2033.73,145.421);
   cutg->SetPoint(4,2018.31,100.212);
   cutg->SetPoint(5,1970.75,100.212);
   cutg->SetPoint(6,1972.67,258.443);
   cutg->SetPoint(7,1999.67,288.582);
   cutg->SetPoint(8,2021.52,416.673);
   cutg->SetPoint(9,2044.66,650.251);
   cutg->SetPoint(10,2044.66,687.925);
   TCutG *cutg = new TCutG("CUTLI",11);
   cutg->SetVarX("Tmuca2");
   cutg->SetVarY("Ep");
   cutg->SetTitle("Graph");
   cutg->SetFillColor(1);
   cutg->SetPoint(0,1930.9,695.46);
   cutg->SetPoint(1,1976.96,695.46);
   cutg->SetPoint(2,1961.96,431.743);
   cutg->SetPoint(3,1922.33,273.512);
   cutg->SetPoint(4,1880.55,183.095);
   cutg->SetPoint(5,1851.63,115.282);
   cutg->SetPoint(6,1832.35,85.1429);
   cutg->SetPoint(7,1819.5,92.6777);
   cutg->SetPoint(8,1836.64,250.908);
   cutg->SetPoint(9,1911.62,492.021);
   cutg->SetPoint(10,1930.9,695.46);
// Selection runs :
	TCut rcut ="Num_run>114&&Num_run<156||Num_run>168&&Num_run<229";
//	TCut rcut ="Num_run>155&&Num_run<161";
//	TCut rcut ="Num_run>185&&Num_run<198";
// Suppression sili :
        TCut scut ="E_sili_1==0&&E_sili_2==0&&E_sili_3==0&&E_sili_4==0&&E_sili_5==0&&E_sili_6==0&&E_sili_7==0&&E_sili_8==0";
// Selections physiques	:
//	TCut pcut ="T_pl1m<3900&&T_pl3m>1150&&E_pl1m>460&&E_pl3m<270&&E_pl2m>200";
	TCut pcut ="T_pl1m>500&&T_pl1m<3900&&T_pl2m<1500&&T_pl3m>1150&&E_pl1m>460&&E_pl3m<270&&E_pl2m>300";
// Selection sur des angles labo :
        TCut acut ="Tetha>130";
// Et de ces dernieres selections : 
        cucu = rcut&&scut&&pcut&&acut;
// Au total 
	cut ="CUTBE"&&cucu; 	 
// histo 100 keV/cl 
//        TH1F h1("h1","11Be residual spectrum",400,-10.,30.);
        TH1F *h1 = new TH1F("h1","11Be residual spectrum",400,-10.,30.);

        ch1.Draw("Ex>>h1");
//        ch1->Draw("Ex",cut);
        ch1->Draw("Ex>>h1", cut);
        c1->Modified();
	c1->cd();
     }