//void hstack(){
{
  TFile *f1=TFile::Open("Trimmer_FullData_final2016.root");
  TFile *f2=TFile::Open("Trimmer_FullData_final2016.root");
  TFile *f3=TFile::Open("Trimmer_FullData_final2016.root");
  TFile *f4=TFile::Open("Trimmer_FullData_final2016.root");
  TFile *f5=TFile::Open("Trimmer_FullData_final2016.root");
  TFile *f6=TFile::Open("Trimmer_FullData_final2016.root");

  TTree *T1 =(TTree*)f1->Get("monopoles");
  TTree *T2 =(TTree*)f2->Get("monopoles");
  TTree *T3 =(TTree*)f3->Get("monopoles");
  TTree *T4 =(TTree*)f4->Get("monopoles");
  TTree *T5 =(TTree*)f5->Get("monopoles");
  TTree *T6 =(TTree*)f6->Get("monopoles");

  TH1F *h1 = new TH1F("h1","",100,0.,30.0);
  TH1F *h2 = new TH1F("h2","",100,0.,30.0);
  TH1F *h3 = new TH1F("h3","",100,0.,30.0);
  TH1F *h4 = new TH1F("h4","",100,0.,30.0);
  TH1F *h5 = new TH1F("h5","",100,0.,30.0);
  TH1F *h6 = new TH1F("h6","",100,0.,30.0);

  // h500->Scale(1./h500->Integral("width"));
 
  //   h500->SetColorLine(kRed);
  //  h700->SetFillColor(46);
  T1->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h1", "HIso < 10 && Dist < 0.6 && abs(XYPar0) < 0.6 && abs(RZPar0) < 10 && abs(XYPar2) > 1000 && abs(RZPar1) < 999 && abs(RZPar2) < 0.005 && clust_E > 175 && 0 < seedFrac < 0.21");

  T2->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h2", "HIso < 10 && Dist < 0.6 && abs(XYPar0) < 0.6 && abs(RZPar0) < 10 && abs(XYPar2) > 1000 && abs(RZPar1) < 999 && abs(RZPar2) < 0.005 && clust_E > 175 && 0.21 <= seedFrac < 0.42");

  T3->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h3", "HIso < 10 && Dist < 0.6 && abs(XYPar0) < 0.6 && abs(RZPar0) < 10 && abs(XYPar2) > 1000 && abs(RZPar1) < 999 && abs(RZPar2) < 0.005 && clust_E > 175 && 0.42 <= seedFrac < 0.63");
 
  T4->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h4", "HIso < 10 && Dist < 0.6 && abs(XYPar0) < 0.6 && abs(RZPar0) < 10 && abs(XYPar2) > 1000 && abs(RZPar1) < 999 && abs(RZPar2) < 0.005 && clust_E > 175 && 0.63 <= seedFrac < 0.84");

  T5->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h5", "HIso < 10 && Dist < 0.6 && abs(XYPar0) < 0.6 && abs(RZPar0) < 10 && abs(XYPar2) > 1000 && abs(RZPar1) < 999 && abs(RZPar2) < 0.005 && clust_E > 175 && 0.84 <= seedFrac < 1.05");

    T6->Draw("sqrt(-TMath::Log(TMath::BinomialI(0.07,SubHits,SatSubHits))) >> h6");

 


  h1->Scale(1./h1->Integral());
  h2->Scale(1./h2->Integral());
  h3->Scale(1./h3->Integral());
  h4->Scale(1./h4->Integral());
  h5->Scale(1./h5->Integral());
  h6->Scale(1./h6->Integral());

  h1->SetLineColor(2);
  //  h500->SetLineWidth(2);

  h2->SetLineColor(3);
  h3->SetLineColor(4);
  h4->SetLineColor(1);
  h5->SetLineColor(2);
  h6->SetLineColor(6);


  h6->Draw();
  h1->Draw("E1 SAME");
  h2->Draw("E1 SAME");
  h3->Draw("E1 SAME");
  h4->Draw("E1 SAME");
  h5->Draw("E1 SAME");

  TLegend leg(0.3, 0.8, 0.6, 0.5);
  leg.AddEntry(h1, "0 < f51 < 0.21");
  leg.AddEntry(h2, "0.21 < f51 < 0.42 ");
  leg.AddEntry(h3, "0.42 < f51 < 0.63");
  leg.AddEntry(h4, "0.63 < f51 < 0.84");
  leg.AddEntry(h5, "0.84 < f51 < 1.05");
  leg.AddEntry(h6, "All Data");

  leg.SetFillStyle(0);
  leg.Draw();
  //  a->Add(h500);
   //   h700->SetColorLine(kGreen);
   //h700->Update();
  // a->Add(h700);
  // a->Draw();
  //  gStyle->SetErrorX();
  gPad-> SetLogy();
  gPad->Modified();
  gPad->Update();
  gPad->SaveAs("dedxsig_slice_FINAL.png");
  gPad->SaveAs("dedxsig_aLLsLICES_FINAL.root");



}