void hRates()
{
//=========Macro generated from canvas: c/c
//=========  (Sun Jan 13 14:40:43 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,-16018.89,67.5,144170);
   c_1->SetFillColor(0);
   c_1->SetBorderMode(0);
   c_1->SetBorderSize(2);
   c_1->SetFrameBorderMode(0);
   c_1->SetFrameBorderMode(0);
   
   TH1D *h_SL1L172 = new TH1D("h_SL1L172","Occupancy removing noisy cells",60,0,60);
   h_SL1L172->SetBinContent(1,20045.71);
   h_SL1L172->SetBinContent(2,20315.71);
   h_SL1L172->SetBinContent(3,19997.14);
   h_SL1L172->SetBinContent(4,19872.86);
   h_SL1L172->SetBinContent(5,19722.86);
   h_SL1L172->SetBinContent(6,20725.71);
   h_SL1L172->SetBinContent(7,20492.86);
   h_SL1L172->SetBinContent(8,20501.43);
   h_SL1L172->SetBinContent(9,19377.14);
   h_SL1L172->SetBinContent(10,17520);
   h_SL1L172->SetBinContent(11,17121.43);
   h_SL1L172->SetBinContent(12,16925.71);
   h_SL1L172->SetBinContent(13,17254.29);
   h_SL1L172->SetBinContent(14,17925.71);
   h_SL1L172->SetBinContent(15,17752.86);
   h_SL1L172->SetBinContent(16,16515.71);
   h_SL1L172->SetBinContent(17,17187.14);
   h_SL1L172->SetBinContent(18,17681.43);
   h_SL1L172->SetBinContent(19,18485.71);
   h_SL1L172->SetBinContent(20,18462.86);
   h_SL1L172->SetBinContent(21,19422.86);
   h_SL1L172->SetBinContent(22,19724.29);
   h_SL1L172->SetBinContent(23,19915.71);
   h_SL1L172->SetBinContent(24,19550);
   h_SL1L172->SetBinContent(25,20271.43);
   h_SL1L172->SetBinContent(26,18791.43);
   h_SL1L172->SetBinContent(27,18977.14);
   h_SL1L172->SetBinContent(28,20071.43);
   h_SL1L172->SetBinContent(29,20474.29);
   h_SL1L172->SetBinContent(30,21114.29);
   h_SL1L172->SetBinContent(31,22172.86);
   h_SL1L172->SetBinContent(32,24978.57);
   h_SL1L172->SetBinContent(33,48585.71);
   h_SL1L172->SetBinContent(34,62555.71);
   h_SL1L172->SetBinContent(35,71777.15);
   h_SL1L172->SetBinContent(36,82542.86);
   h_SL1L172->SetBinContent(37,92795.72);
   h_SL1L172->SetBinContent(38,86264.29);
   h_SL1L172->SetBinContent(39,62450);
   h_SL1L172->SetBinContent(40,25580);
   h_SL1L172->SetBinContent(41,21280);
   h_SL1L172->SetBinContent(42,20054.29);
   h_SL1L172->SetBinContent(43,18758.57);
   h_SL1L172->SetBinContent(44,18321.43);
   h_SL1L172->SetBinContent(45,17628.57);
   h_SL1L172->SetBinContent(46,17472.86);
   h_SL1L172->SetBinContent(47,16281.43);
   h_SL1L172->SetBinContent(48,15997.14);
   h_SL1L172->SetBinContent(49,15981.43);
   h_SL1L172->SetBinContent(50,15881.43);
   h_SL1L172->SetBinContent(51,15012.86);
   h_SL1L172->SetBinContent(52,14980);
   h_SL1L172->SetBinContent(53,14852.86);
   h_SL1L172->SetBinContent(54,14215.71);
   h_SL1L172->SetBinContent(55,14257.14);
   h_SL1L172->SetBinContent(56,13745.71);
   h_SL1L172->SetBinContent(57,13457.14);
   h_SL1L172->SetBinContent(58,12902.86);
   h_SL1L172->SetBinContent(59,12882.86);
   h_SL1L172->SetMaximum(128151.1);
   h_SL1L172->SetEntries(60);
   h_SL1L172->SetDirectory(0);
   h_SL1L172->SetStats(0);
   h_SL1L172->SetLineWidth(2);
   h_SL1L172->GetXaxis()->SetTitle("Cell");
   h_SL1L172->GetXaxis()->SetLabelFont(42);
   h_SL1L172->GetXaxis()->SetLabelSize(0.035);
   h_SL1L172->GetXaxis()->SetTitleSize(0.035);
   h_SL1L172->GetXaxis()->SetTitleFont(42);
   h_SL1L172->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L172->GetYaxis()->SetLabelFont(42);
   h_SL1L172->GetYaxis()->SetLabelSize(0.035);
   h_SL1L172->GetYaxis()->SetTitleSize(0.035);
   h_SL1L172->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L172->GetYaxis()->SetTitleFont(42);
   h_SL1L172->GetZaxis()->SetLabelFont(42);
   h_SL1L172->GetZaxis()->SetLabelSize(0.035);
   h_SL1L172->GetZaxis()->SetTitleSize(0.035);
   h_SL1L172->GetZaxis()->SetTitleFont(42);
   h_SL1L172->Draw("hist, ");
   
   TH1D *h_SL1L273 = new TH1D("h_SL1L273","Occupancy removing noisy cells",60,0,60);
   h_SL1L273->SetBinContent(4,4.285714);
   h_SL1L273->SetBinContent(7,20);
   h_SL1L273->SetBinContent(8,1.428571);
   h_SL1L273->SetBinContent(9,40);
   h_SL1L273->SetBinContent(15,1.428571);
   h_SL1L273->SetBinContent(23,7.142858);
   h_SL1L273->SetBinContent(24,1.428571);
   h_SL1L273->SetBinContent(28,2.857143);
   h_SL1L273->SetBinContent(29,1.428571);
   h_SL1L273->SetBinContent(31,5.714286);
   h_SL1L273->SetBinContent(32,10);
   h_SL1L273->SetBinContent(36,1.428571);
   h_SL1L273->SetBinContent(37,1.428571);
   h_SL1L273->SetBinContent(38,1.428571);
   h_SL1L273->SetBinContent(39,10);
   h_SL1L273->SetBinContent(40,2.857143);
   h_SL1L273->SetBinContent(41,22.85714);
   h_SL1L273->SetBinContent(44,14.28572);
   h_SL1L273->SetBinContent(48,5.714286);
   h_SL1L273->SetBinContent(49,111435.7);
   h_SL1L273->SetBinContent(52,4981.429);
   h_SL1L273->SetBinContent(53,7245.714);
   h_SL1L273->SetBinContent(56,32715.71);
   h_SL1L273->SetMaximum(128151.1);
   h_SL1L273->SetEntries(60);
   h_SL1L273->SetDirectory(0);
   h_SL1L273->SetStats(0);

   Int_t ci;      // for color index setting
   TColor *color; // for color definition with alpha
   ci = TColor::GetColor("#cc0000");
   h_SL1L273->SetLineColor(ci);
   h_SL1L273->SetLineWidth(2);
   h_SL1L273->GetXaxis()->SetTitle("Cell");
   h_SL1L273->GetXaxis()->SetLabelFont(42);
   h_SL1L273->GetXaxis()->SetLabelSize(0.035);
   h_SL1L273->GetXaxis()->SetTitleSize(0.035);
   h_SL1L273->GetXaxis()->SetTitleFont(42);
   h_SL1L273->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L273->GetYaxis()->SetLabelFont(42);
   h_SL1L273->GetYaxis()->SetLabelSize(0.035);
   h_SL1L273->GetYaxis()->SetTitleSize(0.035);
   h_SL1L273->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L273->GetYaxis()->SetTitleFont(42);
   h_SL1L273->GetZaxis()->SetLabelFont(42);
   h_SL1L273->GetZaxis()->SetLabelSize(0.035);
   h_SL1L273->GetZaxis()->SetTitleSize(0.035);
   h_SL1L273->GetZaxis()->SetTitleFont(42);
   h_SL1L273->Draw("hist, same");
   
   TH1D *h_SL1L374 = new TH1D("h_SL1L374","Occupancy removing noisy cells",60,0,60);
   h_SL1L374->SetBinContent(2,1.428571);
   h_SL1L374->SetBinContent(3,1.428571);
   h_SL1L374->SetBinContent(4,4.285714);
   h_SL1L374->SetBinContent(7,8.571428);
   h_SL1L374->SetBinContent(12,1.428571);
   h_SL1L374->SetBinContent(15,5.714286);
   h_SL1L374->SetBinContent(17,1.428571);
   h_SL1L374->SetBinContent(20,2.857143);
   h_SL1L374->SetBinContent(21,1.428571);
   h_SL1L374->SetBinContent(23,7.142858);
   h_SL1L374->SetBinContent(28,1.428571);
   h_SL1L374->SetBinContent(31,14.28572);
   h_SL1L374->SetBinContent(36,5.714286);
   h_SL1L374->SetBinContent(37,2.857143);
   h_SL1L374->SetBinContent(39,11.42857);
   h_SL1L374->SetBinContent(43,1.428571);
   h_SL1L374->SetBinContent(45,4.285714);
   h_SL1L374->SetBinContent(47,15.71429);
   h_SL1L374->SetBinContent(48,4.285714);
   h_SL1L374->SetBinContent(49,1871.429);
   h_SL1L374->SetBinContent(51,1.428571);
   h_SL1L374->SetBinContent(52,4.285714);
   h_SL1L374->SetBinContent(53,42354.29);
   h_SL1L374->SetBinContent(56,46981.43);
   h_SL1L374->SetMaximum(128151.1);
   h_SL1L374->SetEntries(60);
   h_SL1L374->SetDirectory(0);
   h_SL1L374->SetStats(0);

   ci = TColor::GetColor("#009900");
   h_SL1L374->SetLineColor(ci);
   h_SL1L374->SetLineWidth(2);
   h_SL1L374->GetXaxis()->SetTitle("Cell");
   h_SL1L374->GetXaxis()->SetLabelFont(42);
   h_SL1L374->GetXaxis()->SetLabelSize(0.035);
   h_SL1L374->GetXaxis()->SetTitleSize(0.035);
   h_SL1L374->GetXaxis()->SetTitleFont(42);
   h_SL1L374->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L374->GetYaxis()->SetLabelFont(42);
   h_SL1L374->GetYaxis()->SetLabelSize(0.035);
   h_SL1L374->GetYaxis()->SetTitleSize(0.035);
   h_SL1L374->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L374->GetYaxis()->SetTitleFont(42);
   h_SL1L374->GetZaxis()->SetLabelFont(42);
   h_SL1L374->GetZaxis()->SetLabelSize(0.035);
   h_SL1L374->GetZaxis()->SetTitleSize(0.035);
   h_SL1L374->GetZaxis()->SetTitleFont(42);
   h_SL1L374->Draw("hist, same");
   
   TH1D *h_SL1L475 = new TH1D("h_SL1L475","Occupancy removing noisy cells",60,0,60);
   h_SL1L475->SetBinContent(2,18512.86);
   h_SL1L475->SetBinContent(3,18058.57);
   h_SL1L475->SetBinContent(4,18238.57);
   h_SL1L475->SetBinContent(5,17994.29);
   h_SL1L475->SetBinContent(6,18768.57);
   h_SL1L475->SetBinContent(7,18747.14);
   h_SL1L475->SetBinContent(8,18800);
   h_SL1L475->SetBinContent(9,17991.43);
   h_SL1L475->SetBinContent(10,16837.14);
   h_SL1L475->SetBinContent(11,15851.43);
   h_SL1L475->SetBinContent(12,15335.71);
   h_SL1L475->SetBinContent(13,15347.14);
   h_SL1L475->SetBinContent(14,15697.14);
   h_SL1L475->SetBinContent(15,16062.86);
   h_SL1L475->SetBinContent(16,15578.57);
   h_SL1L475->SetBinContent(17,15062.86);
   h_SL1L475->SetBinContent(18,15575.71);
   h_SL1L475->SetBinContent(19,16251.43);
   h_SL1L475->SetBinContent(20,16628.57);
   h_SL1L475->SetBinContent(21,16972.86);
   h_SL1L475->SetBinContent(22,17570);
   h_SL1L475->SetBinContent(23,17702.86);
   h_SL1L475->SetBinContent(24,17781.43);
   h_SL1L475->SetBinContent(25,18327.14);
   h_SL1L475->SetBinContent(26,17795.71);
   h_SL1L475->SetBinContent(27,17274.29);
   h_SL1L475->SetBinContent(28,17627.14);
   h_SL1L475->SetBinContent(29,18790);
   h_SL1L475->SetBinContent(30,19320);
   h_SL1L475->SetBinContent(31,19972.86);
   h_SL1L475->SetBinContent(32,21907.14);
   h_SL1L475->SetBinContent(33,31925.71);
   h_SL1L475->SetBinContent(34,56208.57);
   h_SL1L475->SetBinContent(35,66882.86);
   h_SL1L475->SetBinContent(36,75797.15);
   h_SL1L475->SetBinContent(37,87365.72);
   h_SL1L475->SetBinContent(38,90308.57);
   h_SL1L475->SetBinContent(39,79202.86);
   h_SL1L475->SetBinContent(40,36058.57);
   h_SL1L475->SetBinContent(41,22424.29);
   h_SL1L475->SetBinContent(42,19414.29);
   h_SL1L475->SetBinContent(43,17981.43);
   h_SL1L475->SetBinContent(44,17261.43);
   h_SL1L475->SetBinContent(45,16608.57);
   h_SL1L475->SetBinContent(46,16545.71);
   h_SL1L475->SetBinContent(47,15798.57);
   h_SL1L475->SetBinContent(48,14785.71);
   h_SL1L475->SetBinContent(49,14290);
   h_SL1L475->SetBinContent(50,14642.86);
   h_SL1L475->SetBinContent(51,14384.29);
   h_SL1L475->SetBinContent(52,13502.86);
   h_SL1L475->SetBinContent(53,13475.71);
   h_SL1L475->SetBinContent(54,13035.71);
   h_SL1L475->SetBinContent(56,12780);
   h_SL1L475->SetBinContent(57,12475.71);
   h_SL1L475->SetBinContent(58,12220);
   h_SL1L475->SetBinContent(59,11880);
   h_SL1L475->SetMaximum(128151.1);
   h_SL1L475->SetEntries(60);
   h_SL1L475->SetDirectory(0);
   h_SL1L475->SetStats(0);

   ci = TColor::GetColor("#0000cc");
   h_SL1L475->SetLineColor(ci);
   h_SL1L475->SetLineWidth(2);
   h_SL1L475->GetXaxis()->SetTitle("Cell");
   h_SL1L475->GetXaxis()->SetLabelFont(42);
   h_SL1L475->GetXaxis()->SetLabelSize(0.035);
   h_SL1L475->GetXaxis()->SetTitleSize(0.035);
   h_SL1L475->GetXaxis()->SetTitleFont(42);
   h_SL1L475->GetYaxis()->SetTitle("Hz per wire");
   h_SL1L475->GetYaxis()->SetLabelFont(42);
   h_SL1L475->GetYaxis()->SetLabelSize(0.035);
   h_SL1L475->GetYaxis()->SetTitleSize(0.035);
   h_SL1L475->GetYaxis()->SetTitleOffset(1.45);
   h_SL1L475->GetYaxis()->SetTitleFont(42);
   h_SL1L475->GetZaxis()->SetLabelFont(42);
   h_SL1L475->GetZaxis()->SetLabelSize(0.035);
   h_SL1L475->GetZaxis()->SetTitleSize(0.035);
   h_SL1L475->GetZaxis()->SetTitleFont(42);
   h_SL1L475->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,-12.73214,67.5,114.5893);
   c_2->SetFillColor(0);
   c_2->SetBorderMode(0);
   c_2->SetBorderSize(2);
   c_2->SetFrameBorderMode(0);
   c_2->SetFrameBorderMode(0);
   
   TH1D *h_SL2L176 = new TH1D("h_SL2L176","Occupancy removing noisy cells",60,0,60);
   h_SL2L176->SetMaximum(101.8571);
   h_SL2L176->SetEntries(60);
   h_SL2L176->SetDirectory(0);
   h_SL2L176->SetStats(0);
   h_SL2L176->SetLineWidth(2);
   h_SL2L176->GetXaxis()->SetTitle("Cell");
   h_SL2L176->GetXaxis()->SetLabelFont(42);
   h_SL2L176->GetXaxis()->SetLabelSize(0.035);
   h_SL2L176->GetXaxis()->SetTitleSize(0.035);
   h_SL2L176->GetXaxis()->SetTitleFont(42);
   h_SL2L176->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L176->GetYaxis()->SetLabelFont(42);
   h_SL2L176->GetYaxis()->SetLabelSize(0.035);
   h_SL2L176->GetYaxis()->SetTitleSize(0.035);
   h_SL2L176->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L176->GetYaxis()->SetTitleFont(42);
   h_SL2L176->GetZaxis()->SetLabelFont(42);
   h_SL2L176->GetZaxis()->SetLabelSize(0.035);
   h_SL2L176->GetZaxis()->SetTitleSize(0.035);
   h_SL2L176->GetZaxis()->SetTitleFont(42);
   h_SL2L176->Draw("hist, ");
   
   TH1D *h_SL2L277 = new TH1D("h_SL2L277","Occupancy removing noisy cells",60,0,60);
   h_SL2L277->SetMaximum(101.8571);
   h_SL2L277->SetEntries(60);
   h_SL2L277->SetDirectory(0);
   h_SL2L277->SetStats(0);

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

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

   ci = TColor::GetColor("#0000cc");
   h_SL2L479->SetLineColor(ci);
   h_SL2L479->SetLineWidth(2);
   h_SL2L479->GetXaxis()->SetTitle("Cell");
   h_SL2L479->GetXaxis()->SetLabelFont(42);
   h_SL2L479->GetXaxis()->SetLabelSize(0.035);
   h_SL2L479->GetXaxis()->SetTitleSize(0.035);
   h_SL2L479->GetXaxis()->SetTitleFont(42);
   h_SL2L479->GetYaxis()->SetTitle("Hz per wire");
   h_SL2L479->GetYaxis()->SetLabelFont(42);
   h_SL2L479->GetYaxis()->SetLabelSize(0.035);
   h_SL2L479->GetYaxis()->SetTitleSize(0.035);
   h_SL2L479->GetYaxis()->SetTitleOffset(1.45);
   h_SL2L479->GetYaxis()->SetTitleFont(42);
   h_SL2L479->GetZaxis()->SetLabelFont(42);
   h_SL2L479->GetZaxis()->SetLabelSize(0.035);
   h_SL2L479->GetZaxis()->SetTitleSize(0.035);
   h_SL2L479->GetZaxis()->SetTitleFont(42);
   h_SL2L479->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,-272.509,67.5,2452.58);
   c_3->SetFillColor(0);
   c_3->SetBorderMode(0);
   c_3->SetBorderSize(2);
   c_3->SetFrameBorderMode(0);
   c_3->SetFrameBorderMode(0);
   
   TH1D *h_SL3L180 = new TH1D("h_SL3L180","Occupancy removing noisy cells",60,0,60);
   h_SL3L180->SetMaximum(2180.072);
   h_SL3L180->SetEntries(60);
   h_SL3L180->SetDirectory(0);
   h_SL3L180->SetStats(0);
   h_SL3L180->SetLineWidth(2);
   h_SL3L180->GetXaxis()->SetTitle("Cell");
   h_SL3L180->GetXaxis()->SetLabelFont(42);
   h_SL3L180->GetXaxis()->SetLabelSize(0.035);
   h_SL3L180->GetXaxis()->SetTitleSize(0.035);
   h_SL3L180->GetXaxis()->SetTitleFont(42);
   h_SL3L180->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L180->GetYaxis()->SetLabelFont(42);
   h_SL3L180->GetYaxis()->SetLabelSize(0.035);
   h_SL3L180->GetYaxis()->SetTitleSize(0.035);
   h_SL3L180->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L180->GetYaxis()->SetTitleFont(42);
   h_SL3L180->GetZaxis()->SetLabelFont(42);
   h_SL3L180->GetZaxis()->SetLabelSize(0.035);
   h_SL3L180->GetZaxis()->SetTitleSize(0.035);
   h_SL3L180->GetZaxis()->SetTitleFont(42);
   h_SL3L180->Draw("hist, ");
   
   TH1D *h_SL3L281 = new TH1D("h_SL3L281","Occupancy removing noisy cells",60,0,60);
   h_SL3L281->SetBinContent(41,8.571428);
   h_SL3L281->SetMaximum(2180.072);
   h_SL3L281->SetEntries(60);
   h_SL3L281->SetDirectory(0);
   h_SL3L281->SetStats(0);

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

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

   ci = TColor::GetColor("#0000cc");
   h_SL3L483->SetLineColor(ci);
   h_SL3L483->SetLineWidth(2);
   h_SL3L483->GetXaxis()->SetTitle("Cell");
   h_SL3L483->GetXaxis()->SetLabelFont(42);
   h_SL3L483->GetXaxis()->SetLabelSize(0.035);
   h_SL3L483->GetXaxis()->SetTitleSize(0.035);
   h_SL3L483->GetXaxis()->SetTitleFont(42);
   h_SL3L483->GetYaxis()->SetTitle("Hz per wire");
   h_SL3L483->GetYaxis()->SetLabelFont(42);
   h_SL3L483->GetYaxis()->SetLabelSize(0.035);
   h_SL3L483->GetYaxis()->SetTitleSize(0.035);
   h_SL3L483->GetYaxis()->SetTitleOffset(1.45);
   h_SL3L483->GetYaxis()->SetTitleFont(42);
   h_SL3L483->GetZaxis()->SetLabelFont(42);
   h_SL3L483->GetZaxis()->SetLabelSize(0.035);
   h_SL3L483->GetZaxis()->SetTitleSize(0.035);
   h_SL3L483->GetZaxis()->SetTitleFont(42);
   h_SL3L483->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,-3547.127,14,31924.14);
   c_6->SetFillColor(0);
   c_6->SetBorderMode(0);
   c_6->SetBorderSize(2);
   c_6->SetFrameBorderMode(0);
   c_6->SetFrameBorderMode(0);
   
   TH1D *hRateSub84 = new TH1D("hRateSub84","Average rate per layer for run 5119 (no noisy cells)",12,0.5,12.5);
   hRateSub84->SetBinContent(1,24675.67);
   hRateSub84->SetBinContent(2,12339.45);
   hRateSub84->SetBinContent(3,9131);
   hRateSub84->SetBinContent(4,23703.96);
   hRateSub84->SetBinContent(10,1.714286);
   hRateSub84->SetMinimum(0);
   hRateSub84->SetMaximum(28377.02);
   hRateSub84->SetEntries(2461165);
   hRateSub84->SetDirectory(0);
   hRateSub84->SetStats(0);

   ci = TColor::GetColor("#330066");
   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.06);
   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");
   
   TH1D *hRateSub85 = new TH1D("hRateSub85","Average rate per layer for run 5119 (no noisy cells)",12,0.5,12.5);
   hRateSub85->SetBinContent(1,16769.63);
   hRateSub85->SetBinContent(2,12355.38);
   hRateSub85->SetBinContent(3,9148);
   hRateSub85->SetBinContent(4,15694.95);
   hRateSub85->SetBinContent(10,1.6);
   hRateSub85->SetEntries(2461165);
   hRateSub85->SetDirectory(0);
   hRateSub85->SetStats(0);

   ci = TColor::GetColor("#cc6600");
   hRateSub85->SetLineColor(ci);
   hRateSub85->SetLineWidth(2);
   hRateSub85->GetXaxis()->SetBinLabel(1,"SL1L1");
   hRateSub85->GetXaxis()->SetBinLabel(2,"SL1L2");
   hRateSub85->GetXaxis()->SetBinLabel(3,"SL1L3");
   hRateSub85->GetXaxis()->SetBinLabel(4,"SL1L4");
   hRateSub85->GetXaxis()->SetBinLabel(5,"SL2L1");
   hRateSub85->GetXaxis()->SetBinLabel(6,"SL2L2");
   hRateSub85->GetXaxis()->SetBinLabel(7,"SL2L3");
   hRateSub85->GetXaxis()->SetBinLabel(8,"SL2L4");
   hRateSub85->GetXaxis()->SetBinLabel(9,"SL3L1");
   hRateSub85->GetXaxis()->SetBinLabel(10,"SL3L2");
   hRateSub85->GetXaxis()->SetBinLabel(11,"SL3L3");
   hRateSub85->GetXaxis()->SetBinLabel(12,"SL3L4");
   hRateSub85->GetXaxis()->SetLabelFont(42);
   hRateSub85->GetXaxis()->SetLabelSize(0.035);
   hRateSub85->GetXaxis()->SetTitleSize(0.035);
   hRateSub85->GetXaxis()->SetTitleFont(42);
   hRateSub85->GetYaxis()->SetTitle("Mean Hz/Cell");
   hRateSub85->GetYaxis()->SetLabelFont(42);
   hRateSub85->GetYaxis()->SetLabelSize(0.035);
   hRateSub85->GetYaxis()->SetTitleSize(0.035);
   hRateSub85->GetYaxis()->SetTitleOffset(1.3);
   hRateSub85->GetYaxis()->SetTitleFont(42);
   hRateSub85->GetZaxis()->SetLabelFont(42);
   hRateSub85->GetZaxis()->SetLabelSize(0.035);
   hRateSub85->GetZaxis()->SetTitleSize(0.035);
   hRateSub85->GetZaxis()->SetTitleFont(42);
   hRateSub85->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 5119");
   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 5119 (no noisy cells)");
   pt->Draw();
   c_6->Modified();
   c->cd();
   c->Modified();
   c->cd();
   c->SetSelected(c);
}