void CerenkovTest3()
{

  const char *const filenames[5] = {
    "QweakSim_03_09_2000_kink.root",
    "QweakSim_03_11_3000_kink_tilt5.8.root",
    "QweakSim_03_09_2000_kink_tilt10.8.root",
    "QweakSim_03_09_2000_kink_tilt20.8.root", 
    "QweakSim_03_11_3000_kink_tilt25.8.root", 
//     "QweakSim_03_09_2000_kink_tilt15.8.root",
  }

  TCanvas *c1 = new TCanvas("c1","",600,400);
  TProfile *prfl[5];
  TProfile *prfr[5];
  
  for(int j = 0; j < 5; j++){
    prfl[j] = new TProfile("prfl","Left PMT Hits vs. X",30,-108,108);
    prfl[j]->SetDirectory(0);
    prfr[j] = new TProfile("prfr","Right PMT Hits vs. X",30,-108,108);
    prfr[j]->SetDirectory(0);
    TFile *f = new TFile(filenames[j],"READ");
    TTree *tree = (TTree *) gROOT->FindObject("Test_NT");
    Int_t NEvents     = tree->GetEntries(); 
    tree->Draw("Cerenkov.Detector.HitLocalPositionX:Cerenkov.PMT.PMTLeftNbOfHits:Cerenkov.PMT.PMTRightNbOfHits","","", 10000, 0);
    for(int i = 0; i < NEvents; i++){
      prfl[j]->Fill(tree->GetV1()[i],tree->GetV2()[i]);
      prfr[j]->Fill(tree->GetV1()[i],tree->GetV3()[i]);
    }
    f->Close();
    delete f;
  }  

  Int_t color[5] = {kBlack,kMagenta,kRed,kGreen,kBlue};
  TCanvas *c2 = new TCanvas("c1","",1200,250);
  gStyle->SetOptStat(kFALSE);
  c2->Divide(2,1);
  c2->SetFillColor(kWhite);
  prfl[0]->SetLineColor(color[0]);
  prfr[0]->SetLineColor(color[0]);
  prfl[0]->GetYaxis()->SetRangeUser(0,200);
  prfr[0]->GetYaxis()->SetRangeUser(0,200);
  c2->cd(1);
  gPad->SetFillColor(kWhite);
  prfl[0]->Draw();
  c2->cd(2);
  gPad->SetFillColor(kWhite);
  prfr[0]->Draw();

  for(int j = 1; j < 5; j++){
    prfl[j]->SetLineColor(color[j]);
    prfr[j]->SetLineColor(color[j]);
    c2->cd(1);
    prfl[j]->Draw("same");
    c2->cd(2);
    prfr[j]->Draw("same");
  }

  const char *const legendentry[5] = {
    "Tilt = 0 degrees",
    "Tilt = 5.8",
    "Tilt = 10.8",
    "Tilt = 20.8",
    "Tilt = 25.8",
  }

  c2->cd(1);
  TLegend *leg1 = new TLegend(0.58,0.82,0.98,0.98);
  for(int j = 0; j < 5; j++) leg1->AddEntry(prfl[j],legendentry[j],"l");
  leg1->SetTextFont(62);
  leg1->SetTextSize(0.030);
  leg1->Draw();
}