void hRates()
{
//=========Macro generated from canvas: c/c
//=========  (Sun Jan 13 14:40:53 2019) by ROOT version6.02/13
   TCanvas *c = new TCanvas("c", "c",0,0,1500,800);
   gStyle->SetOptStat(0);
   c->Range(0,0,1,1);
   c->SetFillColor(0);
   c->SetBorderMode(0);
   c->SetBorderSize(2);
   c->SetFrameBorderMode(0);
  
// ------------>Primitives in pad: c_1
   TPad *c_1 = new TPad("c_1", "c_1",0.01,0.51,0.3233333,0.99);
   c_1->Draw();
   c_1->cd();
   c_1->Range(-7.500001,-15182.62,67.5,136643.6);
   c_1->SetFillColor(0);
   c_1->SetBorderMode(0);
   c_1->SetBorderSize(2);
   c_1->SetFrameBorderMode(0);
   c_1->SetFrameBorderMode(0);
   
   TH1D *h_SL1L171 = new TH1D("h_SL1L171","Occupancy removing noisy cells",60,0,60);
   h_SL1L171->SetBinContent(1,47785.26);
   h_SL1L171->SetBinContent(2,47985.12);
   h_SL1L171->SetBinContent(3,47355.57);
   h_SL1L171->SetBinContent(4,47268.14);
   h_SL1L171->SetBinContent(5,46423.75);
   h_SL1L171->SetBinContent(6,49232.96);
   h_SL1L171->SetBinContent(7,48404.81);
   h_SL1L171->SetBinContent(8,48809.52);
   h_SL1L171->SetBinContent(9,45056);
   h_SL1L171->SetBinContent(10,41677.2);
   h_SL1L171->SetBinContent(11,39275.2);
   h_SL1L171->SetBinContent(12,38387.09);
   h_SL1L171->SetBinContent(13,39411.35);
   h_SL1L171->SetBinContent(14,41147.59);
   h_SL1L171->SetBinContent(15,42163.1);
   h_SL1L171->SetBinContent(16,38353.37);
   h_SL1L171->SetBinContent(17,39282.7);
   h_SL1L171->SetBinContent(18,41230.03);
   h_SL1L171->SetBinContent(19,42551.57);
   h_SL1L171->SetBinContent(20,42993.75);
   h_SL1L171->SetBinContent(21,44744.97);
   h_SL1L171->SetBinContent(22,45731.75);
   h_SL1L171->SetBinContent(23,45801.7);
   h_SL1L171->SetBinContent(24,45811.7);
   h_SL1L171->SetBinContent(25,46032.79);
   h_SL1L171->SetBinContent(26,42490.36);
   h_SL1L171->SetBinContent(27,43144.89);
   h_SL1L171->SetBinContent(28,44966.06);
   h_SL1L171->SetBinContent(29,46251.38);
   h_SL1L171->SetBinContent(30,46601.12);
   h_SL1L171->SetBinContent(31,48072.55);
   h_SL1L171->SetBinContent(32,52972.74);
   h_SL1L171->SetBinContent(33,102750.3);
   h_SL1L171->SetBinContent(34,105074.9);
   h_SL1L171->SetBinContent(35,101701.1);
   h_SL1L171->SetBinContent(36,100480.7);
   h_SL1L171->SetBinContent(37,95630.48);
   h_SL1L171->SetBinContent(38,83496.8);
   h_SL1L171->SetBinContent(39,63197.81);
   h_SL1L171->SetBinContent(40,43764.44);
   h_SL1L171->SetBinContent(41,41268.75);
   h_SL1L171->SetBinContent(42,40513.05);
   h_SL1L171->SetBinContent(43,39215.24);
   h_SL1L171->SetBinContent(44,38674.39);
   h_SL1L171->SetBinContent(45,37380.33);
   h_SL1L171->SetBinContent(46,37277.9);
   h_SL1L171->SetBinContent(47,35219.4);
   h_SL1L171->SetBinContent(48,33792.93);
   h_SL1L171->SetBinContent(49,34676.04);
   h_SL1L171->SetBinContent(50,34852.16);
   h_SL1L171->SetBinContent(51,33906.6);
   h_SL1L171->SetBinContent(52,33488.16);
   h_SL1L171->SetBinContent(53,33298.29);
   h_SL1L171->SetBinContent(54,32808.65);
   h_SL1L171->SetBinContent(55,32632.53);
   h_SL1L171->SetBinContent(56,31789.39);
   h_SL1L171->SetBinContent(57,31272.27);
   h_SL1L171->SetBinContent(58,30240.52);
   h_SL1L171->SetBinContent(59,30527.81);
   h_SL1L171->SetMaximum(121461);
   h_SL1L171->SetEntries(60);
   h_SL1L171->SetDirectory(0);
   h_SL1L171->SetStats(0);
   h_SL1L171->SetLineWidth(2);
   h_SL1L171->GetXaxis()->SetTitle("Cell");
   h_SL1L171->GetXaxis()->SetLabelFont(42);
   h_SL1L171->GetXaxis()->SetLabelSize(0.035);
   h_SL1L171->GetXaxis()->SetTitleSize(0.035);
   h_SL1L171->GetXaxis()->SetTitleFont(42);
   h_SL1L171->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L171->GetYaxis()->SetLabelFont(42);
   h_SL1L171->GetYaxis()->SetLabelSize(0.035);
   h_SL1L171->GetYaxis()->SetTitleSize(0.035);
   h_SL1L171->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L171->GetYaxis()->SetTitleFont(42);
   h_SL1L171->GetZaxis()->SetLabelFont(42);
   h_SL1L171->GetZaxis()->SetLabelSize(0.035);
   h_SL1L171->GetZaxis()->SetTitleSize(0.035);
   h_SL1L171->GetZaxis()->SetTitleFont(42);
   h_SL1L171->Draw("hist, ");
   
   TH1D *h_SL1L272 = new TH1D("h_SL1L272","Occupancy removing noisy cells",60,0,60);
   h_SL1L272->SetBinContent(5,1.249092);
   h_SL1L272->SetBinContent(21,1.249092);
   h_SL1L272->SetBinContent(22,1.249092);
   h_SL1L272->SetBinContent(25,1.249092);
   h_SL1L272->SetBinContent(31,1.249092);
   h_SL1L272->SetBinContent(32,2.498184);
   h_SL1L272->SetBinContent(34,1.249092);
   h_SL1L272->SetBinContent(46,1.249092);
   h_SL1L272->SetBinContent(49,1461.438);
   h_SL1L272->SetBinContent(54,1.249092);
   h_SL1L272->SetBinContent(55,1.249092);
   h_SL1L272->SetBinContent(56,11.24183);
   h_SL1L272->SetMaximum(121461);
   h_SL1L272->SetEntries(60);
   h_SL1L272->SetDirectory(0);
   h_SL1L272->SetStats(0);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#cc0000");
   h_SL1L272->SetLineColor(ci);
   h_SL1L272->SetLineWidth(2);
   h_SL1L272->GetXaxis()->SetTitle("Cell");
   h_SL1L272->GetXaxis()->SetLabelFont(42);
   h_SL1L272->GetXaxis()->SetLabelSize(0.035);
   h_SL1L272->GetXaxis()->SetTitleSize(0.035);
   h_SL1L272->GetXaxis()->SetTitleFont(42);
   h_SL1L272->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L272->GetYaxis()->SetLabelFont(42);
   h_SL1L272->GetYaxis()->SetLabelSize(0.035);
   h_SL1L272->GetYaxis()->SetTitleSize(0.035);
   h_SL1L272->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L272->GetYaxis()->SetTitleFont(42);
   h_SL1L272->GetZaxis()->SetLabelFont(42);
   h_SL1L272->GetZaxis()->SetLabelSize(0.035);
   h_SL1L272->GetZaxis()->SetTitleSize(0.035);
   h_SL1L272->GetZaxis()->SetTitleFont(42);
   h_SL1L272->Draw("hist, same");
   
   TH1D *h_SL1L373 = new TH1D("h_SL1L373","Occupancy removing noisy cells",60,0,60);
   h_SL1L373->SetBinContent(5,1.249092);
   h_SL1L373->SetBinContent(12,1.249092);
   h_SL1L373->SetBinContent(21,1.249092);
   h_SL1L373->SetBinContent(26,2.498184);
   h_SL1L373->SetBinContent(27,2.498184);
   h_SL1L373->SetBinContent(31,1.249092);
   h_SL1L373->SetBinContent(53,43.71822);
   h_SL1L373->SetBinContent(56,81.19098);
   h_SL1L373->SetBinContent(57,1.249092);
   h_SL1L373->SetBinContent(58,1.249092);
   h_SL1L373->SetMaximum(121461);
   h_SL1L373->SetEntries(60);
   h_SL1L373->SetDirectory(0);
   h_SL1L373->SetStats(0);

   ci = TColor::GetColor("#009900");
   h_SL1L373->SetLineColor(ci);
   h_SL1L373->SetLineWidth(2);
   h_SL1L373->GetXaxis()->SetTitle("Cell");
   h_SL1L373->GetXaxis()->SetLabelFont(42);
   h_SL1L373->GetXaxis()->SetLabelSize(0.035);
   h_SL1L373->GetXaxis()->SetTitleSize(0.035);
   h_SL1L373->GetXaxis()->SetTitleFont(42);
   h_SL1L373->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L373->GetYaxis()->SetLabelFont(42);
   h_SL1L373->GetYaxis()->SetLabelSize(0.035);
   h_SL1L373->GetYaxis()->SetTitleSize(0.035);
   h_SL1L373->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L373->GetYaxis()->SetTitleFont(42);
   h_SL1L373->GetZaxis()->SetLabelFont(42);
   h_SL1L373->GetZaxis()->SetLabelSize(0.035);
   h_SL1L373->GetZaxis()->SetTitleSize(0.035);
   h_SL1L373->GetZaxis()->SetTitleFont(42);
   h_SL1L373->Draw("hist, same");
   
   TH1D *h_SL1L474 = new TH1D("h_SL1L474","Occupancy removing noisy cells",60,0,60);
   h_SL1L474->SetBinContent(2,44129.17);
   h_SL1L474->SetBinContent(3,42727.69);
   h_SL1L474->SetBinContent(4,43281.04);
   h_SL1L474->SetBinContent(5,42164.35);
   h_SL1L474->SetBinContent(6,44551.36);
   h_SL1L474->SetBinContent(7,44269.07);
   h_SL1L474->SetBinContent(8,44361.5);
   h_SL1L474->SetBinContent(9,43037.46);
   h_SL1L474->SetBinContent(10,39702.39);
   h_SL1L474->SetBinContent(11,36533.44);
   h_SL1L474->SetBinContent(12,35104.48);
   h_SL1L474->SetBinContent(13,34984.57);
   h_SL1L474->SetBinContent(14,37190.46);
   h_SL1L474->SetBinContent(15,37305.38);
   h_SL1L474->SetBinContent(16,36467.24);
   h_SL1L474->SetBinContent(17,35236.89);
   h_SL1L474->SetBinContent(18,36597.14);
   h_SL1L474->SetBinContent(19,38057.34);
   h_SL1L474->SetBinContent(20,39210.25);
   h_SL1L474->SetBinContent(21,39788.57);
   h_SL1L474->SetBinContent(22,41163.83);
   h_SL1L474->SetBinContent(23,41128.85);
   h_SL1L474->SetBinContent(24,41317.46);
   h_SL1L474->SetBinContent(25,41774.63);
   h_SL1L474->SetBinContent(26,40599.24);
   h_SL1L474->SetBinContent(27,38989.16);
   h_SL1L474->SetBinContent(28,39818.55);
   h_SL1L474->SetBinContent(29,42010.71);
   h_SL1L474->SetBinContent(30,42288.01);
   h_SL1L474->SetBinContent(31,42762.66);
   h_SL1L474->SetBinContent(32,46468.72);
   h_SL1L474->SetBinContent(33,73019.42);
   h_SL1L474->SetBinContent(34,105618.2);
   h_SL1L474->SetBinContent(35,100369.5);
   h_SL1L474->SetBinContent(36,99635.07);
   h_SL1L474->SetBinContent(37,95113.36);
   h_SL1L474->SetBinContent(38,86699.48);
   h_SL1L474->SetBinContent(39,73066.88);
   h_SL1L474->SetBinContent(40,45654.31);
   h_SL1L474->SetBinContent(41,39490.04);
   h_SL1L474->SetBinContent(42,38117.29);
   h_SL1L474->SetBinContent(43,36490.97);
   h_SL1L474->SetBinContent(44,36147.47);
   h_SL1L474->SetBinContent(45,35073.25);
   h_SL1L474->SetBinContent(46,35279.36);
   h_SL1L474->SetBinContent(47,33476.91);
   h_SL1L474->SetBinContent(48,31088.65);
   h_SL1L474->SetBinContent(49,31002.46);
   h_SL1L474->SetBinContent(50,31997.99);
   h_SL1L474->SetBinContent(51,31868.08);
   h_SL1L474->SetBinContent(52,30596.51);
   h_SL1L474->SetBinContent(53,30678.95);
   h_SL1L474->SetBinContent(54,29714.65);
   h_SL1L474->SetBinContent(55,1.249092);
   h_SL1L474->SetBinContent(56,30089.38);
   h_SL1L474->SetBinContent(57,29349.91);
   h_SL1L474->SetBinContent(58,28786.57);
   h_SL1L474->SetBinContent(59,27350.12);
   h_SL1L474->SetMaximum(121461);
   h_SL1L474->SetEntries(60);
   h_SL1L474->SetDirectory(0);
   h_SL1L474->SetStats(0);

   ci = TColor::GetColor("#0000cc");
   h_SL1L474->SetLineColor(ci);
   h_SL1L474->SetLineWidth(2);
   h_SL1L474->GetXaxis()->SetTitle("Cell");
   h_SL1L474->GetXaxis()->SetLabelFont(42);
   h_SL1L474->GetXaxis()->SetLabelSize(0.035);
   h_SL1L474->GetXaxis()->SetTitleSize(0.035);
   h_SL1L474->GetXaxis()->SetTitleFont(42);
   h_SL1L474->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L474->GetYaxis()->SetLabelFont(42);
   h_SL1L474->GetYaxis()->SetLabelSize(0.035);
   h_SL1L474->GetYaxis()->SetTitleSize(0.035);
   h_SL1L474->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L474->GetYaxis()->SetTitleFont(42);
   h_SL1L474->GetZaxis()->SetLabelFont(42);
   h_SL1L474->GetZaxis()->SetLabelSize(0.035);
   h_SL1L474->GetZaxis()->SetTitleSize(0.035);
   h_SL1L474->GetZaxis()->SetTitleFont(42);
   h_SL1L474->Draw("hist, same");
   
   TPaveText *pt = new TPaveText(0.1753311,0.9302332,0.8246689,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   TText *AText = pt->AddText("Occupancy removing noisy cells");
   pt->Draw();
   c_1->Modified();
   c->cd();
  
// ------------>Primitives in pad: c_2
   TPad *c_2 = new TPad("c_2", "c_2",0.3433333,0.51,0.6566667,0.99);
   c_2->Draw();
   c_2->cd();
   c_2->Range(-7.500001,-0.179557,67.5,1.616013);
   c_2->SetFillColor(0);
   c_2->SetBorderMode(0);
   c_2->SetBorderSize(2);
   c_2->SetFrameBorderMode(0);
   c_2->SetFrameBorderMode(0);
   
   TH1D *h_SL2L175 = new TH1D("h_SL2L175","Occupancy removing noisy cells",60,0,60);
   h_SL2L175->SetMaximum(1.436456);
   h_SL2L175->SetEntries(60);
   h_SL2L175->SetDirectory(0);
   h_SL2L175->SetStats(0);
   h_SL2L175->SetLineWidth(2);
   h_SL2L175->GetXaxis()->SetTitle("Cell");
   h_SL2L175->GetXaxis()->SetLabelFont(42);
   h_SL2L175->GetXaxis()->SetLabelSize(0.035);
   h_SL2L175->GetXaxis()->SetTitleSize(0.035);
   h_SL2L175->GetXaxis()->SetTitleFont(42);
   h_SL2L175->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L175->GetYaxis()->SetLabelFont(42);
   h_SL2L175->GetYaxis()->SetLabelSize(0.035);
   h_SL2L175->GetYaxis()->SetTitleSize(0.035);
   h_SL2L175->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L175->GetYaxis()->SetTitleFont(42);
   h_SL2L175->GetZaxis()->SetLabelFont(42);
   h_SL2L175->GetZaxis()->SetLabelSize(0.035);
   h_SL2L175->GetZaxis()->SetTitleSize(0.035);
   h_SL2L175->GetZaxis()->SetTitleFont(42);
   h_SL2L175->Draw("hist, ");
   
   TH1D *h_SL2L276 = new TH1D("h_SL2L276","Occupancy removing noisy cells",60,0,60);
   h_SL2L276->SetMaximum(1.436456);
   h_SL2L276->SetEntries(60);
   h_SL2L276->SetDirectory(0);
   h_SL2L276->SetStats(0);

   ci = TColor::GetColor("#cc0000");
   h_SL2L276->SetLineColor(ci);
   h_SL2L276->SetLineWidth(2);
   h_SL2L276->GetXaxis()->SetTitle("Cell");
   h_SL2L276->GetXaxis()->SetLabelFont(42);
   h_SL2L276->GetXaxis()->SetLabelSize(0.035);
   h_SL2L276->GetXaxis()->SetTitleSize(0.035);
   h_SL2L276->GetXaxis()->SetTitleFont(42);
   h_SL2L276->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L276->GetYaxis()->SetLabelFont(42);
   h_SL2L276->GetYaxis()->SetLabelSize(0.035);
   h_SL2L276->GetYaxis()->SetTitleSize(0.035);
   h_SL2L276->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L276->GetYaxis()->SetTitleFont(42);
   h_SL2L276->GetZaxis()->SetLabelFont(42);
   h_SL2L276->GetZaxis()->SetLabelSize(0.035);
   h_SL2L276->GetZaxis()->SetTitleSize(0.035);
   h_SL2L276->GetZaxis()->SetTitleFont(42);
   h_SL2L276->Draw("hist, same");
   
   TH1D *h_SL2L377 = new TH1D("h_SL2L377","Occupancy removing noisy cells",60,0,60);
   h_SL2L377->SetBinContent(9,1.249092);
   h_SL2L377->SetMaximum(1.436456);
   h_SL2L377->SetEntries(60);
   h_SL2L377->SetDirectory(0);
   h_SL2L377->SetStats(0);

   ci = TColor::GetColor("#009900");
   h_SL2L377->SetLineColor(ci);
   h_SL2L377->SetLineWidth(2);
   h_SL2L377->GetXaxis()->SetTitle("Cell");
   h_SL2L377->GetXaxis()->SetLabelFont(42);
   h_SL2L377->GetXaxis()->SetLabelSize(0.035);
   h_SL2L377->GetXaxis()->SetTitleSize(0.035);
   h_SL2L377->GetXaxis()->SetTitleFont(42);
   h_SL2L377->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L377->GetYaxis()->SetLabelFont(42);
   h_SL2L377->GetYaxis()->SetLabelSize(0.035);
   h_SL2L377->GetYaxis()->SetTitleSize(0.035);
   h_SL2L377->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L377->GetYaxis()->SetTitleFont(42);
   h_SL2L377->GetZaxis()->SetLabelFont(42);
   h_SL2L377->GetZaxis()->SetLabelSize(0.035);
   h_SL2L377->GetZaxis()->SetTitleSize(0.035);
   h_SL2L377->GetZaxis()->SetTitleFont(42);
   h_SL2L377->Draw("hist, same");
   
   TH1D *h_SL2L478 = new TH1D("h_SL2L478","Occupancy removing noisy cells",60,0,60);
   h_SL2L478->SetMaximum(1.436456);
   h_SL2L478->SetEntries(60);
   h_SL2L478->SetDirectory(0);
   h_SL2L478->SetStats(0);

   ci = TColor::GetColor("#0000cc");
   h_SL2L478->SetLineColor(ci);
   h_SL2L478->SetLineWidth(2);
   h_SL2L478->GetXaxis()->SetTitle("Cell");
   h_SL2L478->GetXaxis()->SetLabelFont(42);
   h_SL2L478->GetXaxis()->SetLabelSize(0.035);
   h_SL2L478->GetXaxis()->SetTitleSize(0.035);
   h_SL2L478->GetXaxis()->SetTitleFont(42);
   h_SL2L478->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L478->GetYaxis()->SetLabelFont(42);
   h_SL2L478->GetYaxis()->SetLabelSize(0.035);
   h_SL2L478->GetYaxis()->SetTitleSize(0.035);
   h_SL2L478->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L478->GetYaxis()->SetTitleFont(42);
   h_SL2L478->GetZaxis()->SetLabelFont(42);
   h_SL2L478->GetZaxis()->SetLabelSize(0.035);
   h_SL2L478->GetZaxis()->SetTitleSize(0.035);
   h_SL2L478->GetZaxis()->SetTitleFont(42);
   h_SL2L478->Draw("hist, same");
   
   pt = new TPaveText(0.1753311,0.9302332,0.8246689,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   AText = pt->AddText("Occupancy removing noisy cells");
   pt->Draw();
   c_2->Modified();
   c->cd();
  
// ------------>Primitives in pad: c_3
   TPad *c_3 = new TPad("c_3", "c_3",0.6766667,0.51,0.99,0.99);
   c_3->Draw();
   c_3->cd();
   c_3->Range(-7.500001,-0.7182279,67.5,6.464051);
   c_3->SetFillColor(0);
   c_3->SetBorderMode(0);
   c_3->SetBorderSize(2);
   c_3->SetFrameBorderMode(0);
   c_3->SetFrameBorderMode(0);
   
   TH1D *h_SL3L179 = new TH1D("h_SL3L179","Occupancy removing noisy cells",60,0,60);
   h_SL3L179->SetBinContent(9,4.996368);
   h_SL3L179->SetBinContent(10,4.996368);
   h_SL3L179->SetBinContent(11,4.996368);
   h_SL3L179->SetBinContent(12,3.747276);
   h_SL3L179->SetMaximum(5.745823);
   h_SL3L179->SetEntries(60);
   h_SL3L179->SetDirectory(0);
   h_SL3L179->SetStats(0);
   h_SL3L179->SetLineWidth(2);
   h_SL3L179->GetXaxis()->SetTitle("Cell");
   h_SL3L179->GetXaxis()->SetLabelFont(42);
   h_SL3L179->GetXaxis()->SetLabelSize(0.035);
   h_SL3L179->GetXaxis()->SetTitleSize(0.035);
   h_SL3L179->GetXaxis()->SetTitleFont(42);
   h_SL3L179->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L179->GetYaxis()->SetLabelFont(42);
   h_SL3L179->GetYaxis()->SetLabelSize(0.035);
   h_SL3L179->GetYaxis()->SetTitleSize(0.035);
   h_SL3L179->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L179->GetYaxis()->SetTitleFont(42);
   h_SL3L179->GetZaxis()->SetLabelFont(42);
   h_SL3L179->GetZaxis()->SetLabelSize(0.035);
   h_SL3L179->GetZaxis()->SetTitleSize(0.035);
   h_SL3L179->GetZaxis()->SetTitleFont(42);
   h_SL3L179->Draw("hist, ");
   
   TH1D *h_SL3L280 = new TH1D("h_SL3L280","Occupancy removing noisy cells",60,0,60);
   h_SL3L280->SetBinContent(49,4.996368);
   h_SL3L280->SetMaximum(5.745823);
   h_SL3L280->SetEntries(60);
   h_SL3L280->SetDirectory(0);
   h_SL3L280->SetStats(0);

   ci = TColor::GetColor("#cc0000");
   h_SL3L280->SetLineColor(ci);
   h_SL3L280->SetLineWidth(2);
   h_SL3L280->GetXaxis()->SetTitle("Cell");
   h_SL3L280->GetXaxis()->SetLabelFont(42);
   h_SL3L280->GetXaxis()->SetLabelSize(0.035);
   h_SL3L280->GetXaxis()->SetTitleSize(0.035);
   h_SL3L280->GetXaxis()->SetTitleFont(42);
   h_SL3L280->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L280->GetYaxis()->SetLabelFont(42);
   h_SL3L280->GetYaxis()->SetLabelSize(0.035);
   h_SL3L280->GetYaxis()->SetTitleSize(0.035);
   h_SL3L280->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L280->GetYaxis()->SetTitleFont(42);
   h_SL3L280->GetZaxis()->SetLabelFont(42);
   h_SL3L280->GetZaxis()->SetLabelSize(0.035);
   h_SL3L280->GetZaxis()->SetTitleSize(0.035);
   h_SL3L280->GetZaxis()->SetTitleFont(42);
   h_SL3L280->Draw("hist, same");
   
   TH1D *h_SL3L381 = new TH1D("h_SL3L381","Occupancy removing noisy cells",60,0,60);
   h_SL3L381->SetBinContent(1,2.498184);
   h_SL3L381->SetBinContent(2,2.498184);
   h_SL3L381->SetBinContent(3,2.498184);
   h_SL3L381->SetBinContent(4,3.747276);
   h_SL3L381->SetBinContent(31,1.249092);
   h_SL3L381->SetMaximum(5.745823);
   h_SL3L381->SetEntries(60);
   h_SL3L381->SetDirectory(0);
   h_SL3L381->SetStats(0);

   ci = TColor::GetColor("#009900");
   h_SL3L381->SetLineColor(ci);
   h_SL3L381->SetLineWidth(2);
   h_SL3L381->GetXaxis()->SetTitle("Cell");
   h_SL3L381->GetXaxis()->SetLabelFont(42);
   h_SL3L381->GetXaxis()->SetLabelSize(0.035);
   h_SL3L381->GetXaxis()->SetTitleSize(0.035);
   h_SL3L381->GetXaxis()->SetTitleFont(42);
   h_SL3L381->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L381->GetYaxis()->SetLabelFont(42);
   h_SL3L381->GetYaxis()->SetLabelSize(0.035);
   h_SL3L381->GetYaxis()->SetTitleSize(0.035);
   h_SL3L381->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L381->GetYaxis()->SetTitleFont(42);
   h_SL3L381->GetZaxis()->SetLabelFont(42);
   h_SL3L381->GetZaxis()->SetLabelSize(0.035);
   h_SL3L381->GetZaxis()->SetTitleSize(0.035);
   h_SL3L381->GetZaxis()->SetTitleFont(42);
   h_SL3L381->Draw("hist, same");
   
   TH1D *h_SL3L482 = new TH1D("h_SL3L482","Occupancy removing noisy cells",60,0,60);
   h_SL3L482->SetMaximum(5.745823);
   h_SL3L482->SetEntries(60);
   h_SL3L482->SetDirectory(0);
   h_SL3L482->SetStats(0);

   ci = TColor::GetColor("#0000cc");
   h_SL3L482->SetLineColor(ci);
   h_SL3L482->SetLineWidth(2);
   h_SL3L482->GetXaxis()->SetTitle("Cell");
   h_SL3L482->GetXaxis()->SetLabelFont(42);
   h_SL3L482->GetXaxis()->SetLabelSize(0.035);
   h_SL3L482->GetXaxis()->SetTitleSize(0.035);
   h_SL3L482->GetXaxis()->SetTitleFont(42);
   h_SL3L482->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L482->GetYaxis()->SetLabelFont(42);
   h_SL3L482->GetYaxis()->SetLabelSize(0.035);
   h_SL3L482->GetYaxis()->SetTitleSize(0.035);
   h_SL3L482->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L482->GetYaxis()->SetTitleFont(42);
   h_SL3L482->GetZaxis()->SetLabelFont(42);
   h_SL3L482->GetZaxis()->SetLabelSize(0.035);
   h_SL3L482->GetZaxis()->SetTitleSize(0.035);
   h_SL3L482->GetZaxis()->SetTitleFont(42);
   h_SL3L482->Draw("hist, same");
   
   pt = new TPaveText(0.1753311,0.9302332,0.8246689,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   AText = pt->AddText("Occupancy removing noisy cells");
   pt->Draw();
   c_3->Modified();
   c->cd();
  
// ------------>Primitives in pad: c_4
   TPad *c_4 = new TPad("c_4", "c_4",0.01,0.01,0.3233333,0.49);
   c_4->Draw();
   c_4->cd();
   c_4->Range(0,0,1,1);
   c_4->SetFillColor(0);
   c_4->SetBorderMode(0);
   c_4->SetBorderSize(2);
   c_4->SetFrameBorderMode(0);
   c_4->Modified();
   c->cd();
  
// ------------>Primitives in pad: c_5
   TPad *c_5 = new TPad("c_5", "c_5",0.3433333,0.01,0.6566667,0.49);
   c_5->Draw();
   c_5->cd();
   c_5->Range(0,0,1,1);
   c_5->SetFillColor(0);
   c_5->SetBorderMode(0);
   c_5->SetBorderSize(2);
   c_5->SetFrameBorderMode(0);
   c_5->Modified();
   c->cd();
  
// ------------>Primitives in pad: c_6
   TPad *c_6 = new TPad("c_6", "c_6",0.05,0,0.95,0.5);
   c_6->Draw();
   c_6->cd();
   c_6->Range(-1,-6788.77,14,61098.93);
   c_6->SetFillColor(0);
   c_6->SetBorderMode(0);
   c_6->SetBorderSize(2);
   c_6->SetFrameBorderMode(0);
   c_6->SetFrameBorderMode(0);
   
   TH1D *hRateSub83 = new TH1D("hRateSub83","Average rate per layer for run 5072 (no noisy cells)",12,0.5,12.5);
   hRateSub83->SetBinContent(1,47226.22);
   hRateSub83->SetBinContent(4,44670.2);
   hRateSub83->SetMinimum(0);
   hRateSub83->SetMaximum(54310.16);
   hRateSub83->SetEntries(4287300);
   hRateSub83->SetDirectory(0);
   hRateSub83->SetStats(0);

   ci = TColor::GetColor("#330066");
   hRateSub83->SetLineColor(ci);
   hRateSub83->SetLineWidth(2);
   hRateSub83->GetXaxis()->SetBinLabel(1,"SL1L1");
   hRateSub83->GetXaxis()->SetBinLabel(2,"SL1L2");
   hRateSub83->GetXaxis()->SetBinLabel(3,"SL1L3");
   hRateSub83->GetXaxis()->SetBinLabel(4,"SL1L4");
   hRateSub83->GetXaxis()->SetBinLabel(5,"SL2L1");
   hRateSub83->GetXaxis()->SetBinLabel(6,"SL2L2");
   hRateSub83->GetXaxis()->SetBinLabel(7,"SL2L3");
   hRateSub83->GetXaxis()->SetBinLabel(8,"SL2L4");
   hRateSub83->GetXaxis()->SetBinLabel(9,"SL3L1");
   hRateSub83->GetXaxis()->SetBinLabel(10,"SL3L2");
   hRateSub83->GetXaxis()->SetBinLabel(11,"SL3L3");
   hRateSub83->GetXaxis()->SetBinLabel(12,"SL3L4");
   hRateSub83->GetXaxis()->SetLabelFont(42);
   hRateSub83->GetXaxis()->SetLabelSize(0.035);
   hRateSub83->GetXaxis()->SetTitleSize(0.035);
   hRateSub83->GetXaxis()->SetTitleFont(42);
   hRateSub83->GetYaxis()->SetTitle("Mean Hz/Cell");
   hRateSub83->GetYaxis()->SetLabelFont(42);
   hRateSub83->GetYaxis()->SetLabelSize(0.06);
   hRateSub83->GetYaxis()->SetTitleSize(0.035);
   hRateSub83->GetYaxis()->SetTitleOffset(1.3);
   hRateSub83->GetYaxis()->SetTitleFont(42);
   hRateSub83->GetZaxis()->SetLabelFont(42);
   hRateSub83->GetZaxis()->SetLabelSize(0.035);
   hRateSub83->GetZaxis()->SetTitleSize(0.035);
   hRateSub83->GetZaxis()->SetTitleFont(42);
   hRateSub83->Draw("hist");
   
   TH1D *hRateSub84 = new TH1D("hRateSub84","Average rate per layer for run 5072 (no noisy cells)",12,0.5,12.5);
   hRateSub84->SetBinContent(1,39696.68);
   hRateSub84->SetBinContent(4,36978.31);
   hRateSub84->SetEntries(4287300);
   hRateSub84->SetDirectory(0);
   hRateSub84->SetStats(0);

   ci = TColor::GetColor("#cc6600");
   hRateSub84->SetLineColor(ci);
   hRateSub84->SetLineWidth(2);
   hRateSub84->GetXaxis()->SetBinLabel(1,"SL1L1");
   hRateSub84->GetXaxis()->SetBinLabel(2,"SL1L2");
   hRateSub84->GetXaxis()->SetBinLabel(3,"SL1L3");
   hRateSub84->GetXaxis()->SetBinLabel(4,"SL1L4");
   hRateSub84->GetXaxis()->SetBinLabel(5,"SL2L1");
   hRateSub84->GetXaxis()->SetBinLabel(6,"SL2L2");
   hRateSub84->GetXaxis()->SetBinLabel(7,"SL2L3");
   hRateSub84->GetXaxis()->SetBinLabel(8,"SL2L4");
   hRateSub84->GetXaxis()->SetBinLabel(9,"SL3L1");
   hRateSub84->GetXaxis()->SetBinLabel(10,"SL3L2");
   hRateSub84->GetXaxis()->SetBinLabel(11,"SL3L3");
   hRateSub84->GetXaxis()->SetBinLabel(12,"SL3L4");
   hRateSub84->GetXaxis()->SetLabelFont(42);
   hRateSub84->GetXaxis()->SetLabelSize(0.035);
   hRateSub84->GetXaxis()->SetTitleSize(0.035);
   hRateSub84->GetXaxis()->SetTitleFont(42);
   hRateSub84->GetYaxis()->SetTitle("Mean Hz/Cell");
   hRateSub84->GetYaxis()->SetLabelFont(42);
   hRateSub84->GetYaxis()->SetLabelSize(0.035);
   hRateSub84->GetYaxis()->SetTitleSize(0.035);
   hRateSub84->GetYaxis()->SetTitleOffset(1.3);
   hRateSub84->GetYaxis()->SetTitleFont(42);
   hRateSub84->GetZaxis()->SetLabelFont(42);
   hRateSub84->GetZaxis()->SetLabelSize(0.035);
   hRateSub84->GetZaxis()->SetTitleSize(0.035);
   hRateSub84->GetZaxis()->SetTitleFont(42);
   hRateSub84->Draw("hist, same");
   
   TLegend *leg = new TLegend(0.7,0.92,0.9,0.97,NULL,"brNDC");
   leg->SetBorderSize(1);
   leg->SetLineColor(1);
   leg->SetLineStyle(1);
   leg->SetLineWidth(1);
   leg->SetFillColor(0);
   leg->SetFillStyle(1001);
   TLegendEntry *entry=leg->AddEntry("hRateSub","Occupancy","l");

   ci = TColor::GetColor("#330066");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry->SetTextFont(42);
   entry=leg->AddEntry("hRateSub","Background","l");

   ci = TColor::GetColor("#cc6600");
   entry->SetLineColor(ci);
   entry->SetLineStyle(1);
   entry->SetLineWidth(2);
   entry->SetMarkerColor(1);
   entry->SetMarkerStyle(21);
   entry->SetMarkerSize(1);
   entry->SetTextFont(42);
   leg->Draw();
   
   pt = new TPaveText(0.91,0.45,1,0.55,"brNDC");
   pt->SetFillColor(0);
   pt->SetTextAlign(12);
   pt->SetTextFont(42);
   pt->SetTextSize(0.07);
   AText = pt->AddText("Run 5072");
   pt->Draw();
   
   pt = new TPaveText(0.3113755,0.9328238,0.6886245,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   AText = pt->AddText("Average rate per layer for run 5072 (no noisy cells)");
   pt->Draw();
   c_6->Modified();
   c->cd();
   c->Modified();
   c->cd();
   c->SetSelected(c);
}