{
	gROOT->Reset();
	gStyle->SetOptStat(1111111);
	gStyle->SetOptDate(1);



        TH1F *TDC0=new TH1F("TDC0","Neutron stop F1(left)",800,0.5,1000);
	TH1F *TDC1=new TH1F("TDC1","Neutron stop F2(right)",800,0.5,1000);

	TH1F *TDC2=new TH1F("TDC2","Neutron stop G1(left)",800,0.5,1000);
        TH1F *TDC3=new TH1F("TDC3","Neutron stop G2(right)",800,0.5,1000);

	TH1F *TDC4=new TH1F("TDC4","Neutron stop H1(left)",800,0.5,1000);
	TH1F *TDC5=new TH1F("TDC5","Neutron stop H2(right)",800,0.5,1000);

        TH1F *TDC6=new TH1F("TDC6","Neutron stop I1(left)",800,0.5,1000);
	TH1F *TDC7=new TH1F("TDC7","Neutron stop I2(right)",800,0.5,1000);
        

        TH1F *TDC10=new TH1F("TDC10","Neutron stop M1(left)",800,0.5,1000);
        TH1F *TDC11=new TH1F("TDC11","Neutron stop M2(right)",800,0.5,1000);
   
        TH1F *TDC14=new TH1F("TDC14","Neutron stop E1(left)",800,0.5,1000);
        TH1F *TDC15=new TH1F("TDC15","Neutron stop E2(right)",800,0.5,1000);

        TH1F *TDC16=new TH1F("TDC16","Neutron stop K1(left)",800,0.5,1000);
        TH1F *TDC17=new TH1F("TDC17","Neutron stop K2(right)",800,0.5,1000);



  TH1F *Fsubtr1=new TH1F("Fsubtr1","Subtraction of the stops for E",3000,-1000,1000);
  TH1F *Fsubtr2=new TH1F("Fsubtr2","Subtraction of the stops for M",3000,-1000,1000);
  TH1F *Fsubtr3=new TH1F("Fsubtr3","Subtraction of the stops for F",3000,-1000,1000);
  TH1F *Fsubtr4=new TH1F("Fsubtr4","Subtraction of the stops for G",3000,-1000,1000);
  TH1F *Fsubtr5=new TH1F("Fsubtr5","Subtraction of the stops for H",3000,-1000,1000);
  TH1F *Fsubtr6=new TH1F("Fsubtr6","Subtraction of the stops for K",3000,-1000,1000);
  TH1F *Fsubtr7=new TH1F("Fsubtr7","Subtraction of the stops for I",3000,-1000,1000);


  TH1F *Fsum1=new TH1F("Fsum1","Sum of the stops for E",3000,20,1500);
  TH1F *Fsum2=new TH1F("Fsum2","Sum of the stops for M",3000,20,1500);
  TH1F *Fsum3=new TH1F("Fsum3","Sum of the stops for F",3000,20,1500);
  TH1F *Fsum4=new TH1F("Fsum4","Sum of the stops for G",3000,20,1500);
  TH1F *Fsum5=new TH1F("Fsum5","Sum of the stops for H",3000,20,1500);
  TH1F *Fsum6=new TH1F("Fsum6","Sum of the stops for K",3000,20,1500);
  TH1F *Fsum7=new TH1F("Fsum7","Sum of the stops for I",3000,20,1500);


	TCanvas *c = new TCanvas("c", "TDC", 900, 1000);

c->Divide(2,4);
	c->Update();

c->cd(1);
ntuple->Draw("TDC0*0.223>>TDC0");
	c->cd(1)->SetLogy(1);
	c->cd(1)->SetGrid(1,1);
	TDC0->SetXTitle("TDC ns");
	TDC0->SetYTitle("Events");
	TDC0->SetLineColor(4);
	c->Update();

	

c->cd(2);
ntuple->Draw("TDC1*0.223>>TDC1");
	c->cd(2)->SetLogy(1);
	c->cd(2)->SetGrid(1,1);
	TDC1->SetXTitle("TDC ns");
	TDC1->SetYTitle("Events");
	TDC1->SetLineColor(4);
	c->Update();


c->cd(3);
ntuple->Draw("TDC2*0.223>>TDC2");
	c->cd(3)->SetLogy(1);
	c->cd(3)->SetGrid(1,1);
	TDC2->SetXTitle("TDC ns");
	TDC2->SetYTitle("Events");
	TDC2->SetLineColor(4);
	c->Update();


c->cd(4);
ntuple->Draw("TDC3*0.223>>TDC3");
	c->cd(4)->SetLogy(1);
	c->cd(4)->SetGrid(1,1);
	TDC3->SetXTitle("TDC ns");
	TDC3->SetYTitle("Events");
	TDC3->SetLineColor(4);
	c->Update();

c->cd(5);
ntuple->Draw("TDC4*0.223>>TDC4");
	c->cd(5)->SetLogy(1);
	c->cd(5)->SetGrid(1,1);
	TDC4->SetXTitle("TDC ns");
	TDC4->SetYTitle("Events");
	TDC4->SetLineColor(4);
	c->Update();


c->cd(6);
ntuple->Draw("TDC5*0.223>>TDC5");
	c->cd(6)->SetLogy(1);
	c->cd(6)->SetGrid(1,1);
	TDC5->SetXTitle("TDC ns");
	TDC5->SetYTitle("Events");
	TDC5->SetLineColor(4);

c->cd(7);
ntuple->Draw("TDC6*0.223>>TDC6");
	c->cd(7)->SetLogy(1);
	c->cd(7)->SetGrid(1,1);
	TDC6->SetXTitle("TDC ns");
	TDC6->SetYTitle("Events");
	TDC6->SetLineColor(4);

c->cd(8);
ntuple->Draw("TDC7*0.223>>TDC7");
	c->cd(8)->SetLogy(1);
	c->cd(8)->SetGrid(1,1);
	TDC7->SetXTitle("TDC ns");
	TDC7->SetYTitle("Events");
	TDC7->SetLineColor(4);


	c->Update();
	


TCanvas *c7 = new TCanvas("c7", "TDC1", 900, 1000);

c7->Divide(2,3);
	c7->Update();

c7->cd(1);
ntuple->Draw("TDC14*0.223>>TDC14");
	c7->cd(1)->SetLogy(1);
	c7->cd(1)->SetGrid(1,1);
	TDC14->SetXTitle("TDC ns");
	TDC14->SetYTitle("Events");
	TDC14->SetLineColor(4);
	c7->Update();

	

c7->cd(2);
ntuple->Draw("TDC15*0.223>>TDC15");
	c7->cd(2)->SetLogy(1);
	c7->cd(2)->SetGrid(1,1);
	TDC15->SetXTitle("TDC ns");
	TDC15->SetYTitle("Events");
	TDC15->SetLineColor(4);
	c7->Update();


c7->cd(3);
ntuple->Draw("TDC10*0.223>>TDC10");
	c7->cd(3)->SetLogy(1);
	c7->cd(3)->SetGrid(1,1);
	TDC10->SetXTitle("TDC ns");
	TDC10->SetYTitle("Events");
	TDC10->SetLineColor(4);
	c7->Update();

c7->cd(4);
ntuple->Draw("TDC11*0.223>>TDC11");
	c7->cd(4)->SetLogy(1);
	c7->cd(4)->SetGrid(1,1);
	TDC11->SetXTitle("TDC ns");
	TDC11->SetYTitle("Events");
	TDC11->SetLineColor(4);
	c7->Update();


c7->cd(5);
ntuple->Draw("TDC16*0.223>>TDC16");
	c7->cd(5)->SetLogy(1);
	c7->cd(5)->SetGrid(1,1);
	TDC16->SetXTitle("TDC ns");
	TDC16->SetYTitle("Events");
	TDC16->SetLineColor(4);
	c7->Update();


c7->cd(6);
ntuple->Draw("TDC17*0.223>>TDC17");
	c7->cd(6)->SetLogy(1);
	c7->cd(6)->SetGrid(1,1);
	TDC17->SetXTitle("TDC ns");
	TDC17->SetYTitle("Events");
	TDC17->SetLineColor(4);
	c7->Update();



TCanvas *c2 = new TCanvas("c2", "Delta_TDC", 900, 800);

c2->Divide(3,3);
	c2->Update();
//Det F
c2->cd(1);
c2->cd(1)->SetLogy(1);
c2->cd(1)->SetGrid(1,1);
ntuple -> Draw("((((TDC0)-(TDC1))*0.223)+18)*7.37>>Fsubtr3","(((((TDC0)-(TDC1))*0.223)+18)*7.37>=0 && ((((TDC0)-(TDC1))*0.223)+18)*7.37<=220) && ((TDC0)-(TDC1))*0.223 !=0");
Fsubtr3->SetXTitle("TDC [ns]");

//Det G
c2->cd(2);
c2->cd(2)->SetLogy(1);
c2->cd(2)->SetGrid(1,1);
ntuple -> Draw("((((TDC2)-(TDC3))*0.223)+14)*7.24>>Fsubtr4"," ( ((((TDC2)-(TDC3))*0.223)+14)*7.24>=0 && ((((TDC2)-(TDC3))*0.223)+14)*7.24<=230 ) && ((TDC2)-(TDC3))*0.223 !=0");
Fsubtr4->SetXTitle("TDC [ns]");

//Det H
c2->cd(3);
c2->cd(3)->SetLogy(1);
c2->cd(3)->SetGrid(1,1);
ntuple -> Draw("((((TDC4)-(TDC5))*0.223)+10)*7.75>>Fsubtr5"," (((((TDC4)-(TDC5))*0.223)+10)*7.75>=0 &&((((TDC4)-(TDC5))*0.223)+10)*7.75<=190) && ((TDC4)-(TDC5))*0.223 !=0");
Fsubtr5->SetXTitle("TDC [ns]");

//Det I
c2->cd(4);
c2->cd(4)->SetLogy(1);
c2->cd(4)->SetGrid(1,1);
ntuple -> Draw("((((TDC6)-(TDC7))*0.223)+16)*5.48>>Fsubtr7"," ((((TDC6)-(TDC7))*0.223)+16)*5.48>=0 && ((((TDC6)-(TDC7))*0.223)+16)*5.48<=120 && ((TDC6)-(TDC7))*0.223 !=0");
Fsubtr7->SetXTitle("TDC [ns]");

//Det M
c2->cd(5);
c2->cd(5)->SetLogy(1);
c2->cd(5)->SetGrid(1,1);
ntuple -> Draw("((((TDC10)-(TDC11))*0.223)+54)*6.82>>Fsubtr2","(((((TDC10)-(TDC11))*0.223)+54)*6.82>=0 && ((((TDC10)-(TDC11))*0.223)+54)*6.82<=200) && ((TDC10)-(TDC11))*0.223 !=0");
Fsubtr2->SetXTitle("TDC [ns]");

//Det E
c2->cd(6);
c2->cd(6)->SetLogy(1);
c2->cd(6)->SetGrid(1,1);
ntuple -> Draw("((((TDC14)-(TDC15))*0.223)-3)*7.41>>Fsubtr1", "( ((((TDC14)-(TDC15))*0.223)-3)*7.41>=0 && ((((TDC14)-(TDC15))*0.223)-3)*7.41<=120.0 ) && ((TDC14)-(TDC15))*0.223 !=0");
Fsubtr1->SetXTitle("TDC [ns]");

//Det K
c2->cd(7);
c2->cd(7)->SetLogy(1);
c2->cd(7)->SetGrid(1,1);
ntuple -> Draw("((((TDC16)-(TDC17))*0.223)+28)*6.71>>Fsubtr6", "((((TDC16)-(TDC17))*0.223)+28)*6.71>=0 && ((((TDC16)-(TDC17))*0.223)+28)*6.71<=230  &&((TDC16)-(TDC17))*0.223 !=0");
Fsubtr6->SetXTitle("TDC [ns]");

c2->Update();


TCanvas *c3 = new TCanvas("c3", "Sum_TDC", 900, 800);

c3->Divide(3,3);
	c3->Update();

c3->cd(1);
c3->cd(1)->SetLogy(1);
c3->cd(1)->SetGrid(1,1);
ntuple -> Draw("((TDC0)+(TDC1))*0.223>>Fsum3","((TDC0)+(TDC1))*0.223 !=0");
Fsum3->SetXTitle("TDC [ns]");

c3->cd(2);
c3->cd(2)->SetLogy(1);
c3->cd(2)->SetGrid(1,1);
ntuple -> Draw("((TDC2)+(TDC3))*0.223>>Fsum4","((TDC2)+(TDC3))*0.223 !=0");
Fsum4->SetXTitle("TDC [ns]");

c3->cd(3);
c3->cd(3)->SetLogy(1);
c3->cd(3)->SetGrid(1,1);
ntuple -> Draw("((TDC4)+(TDC5))*0.223>>Fsum5","((TDC4)+(TDC5))*0.223 !=0");
Fsum5->SetXTitle("TDC [ns]");


c3->cd(4);
c3->cd(4)->SetLogy(1);
c3->cd(4)->SetGrid(1,1);
ntuple -> Draw("((TDC6)+(TDC7))*0.223>>Fsum7","((TDC6)+(TDC7))*0.223 !=0");
Fsum7->SetXTitle("TDC [ns]");


c3->cd(5);
c3->cd(5)->SetLogy(1);
c3->cd(5)->SetGrid(1,1);
ntuple -> Draw("((TDC10)+(TDC11))*0.223>>Fsum2","((TDC10)+(TDC11))*0.223 !=0");
Fsum2->SetXTitle("TDC [ns]");

c3->cd(6);
c3->cd(6)->SetLogy(1);
c3->cd(6)->SetGrid(1,1);
ntuple -> Draw("((TDC14)+(TDC15))*0.223>>Fsum1", "((TDC14)+(TDC15))*0.223 !=0");
Fsum1->SetXTitle("TDC [ns]");

c3->cd(7);
c3->cd(7)->SetLogy(1);
c3->cd(7)->SetGrid(1,1);
ntuple -> Draw("((TDC16)+(TDC17))*0.223>>Fsum6", "((TDC16)+(TDC17))*0.223 !=0");
Fsum6->SetXTitle("TDC [ns]");

c3->Update();

}