diff --git a/macros/cer_ped.C b/macros/cer_ped.C
new file mode 100644
index 0000000000000000000000000000000000000000..6182a4d3e4d393a2725f9be070ad42fb3a70da1a
--- /dev/null
+++ b/macros/cer_ped.C
@@ -0,0 +1,133 @@
+void cer_ped() {
+
+  // Declare constants
+  static const Double_t clk2tdc = 0.100;
+  static const Double_t clk2adc = 0.0625;
+
+  static const UInt_t nhgc_pmts = 4;
+  static const UInt_t nngc_pmts = 4;
+
+  static const UInt_t maxTdcHits = 128;
+  static const UInt_t maxAdcHits = 4;
+
+  // Heavy gas Cherenkov ADC variables
+  Int_t phgc_hits;
+  Double_t phgc_pmt[maxAdcHits*nhgc_pmts], phgc_pulsePedRaw[maxAdcHits*nhgc_pmts];
+  Double_t phgc_pulseIntRaw[maxAdcHits*nhgc_pmts], phgc_pulseAmpRaw[maxAdcHits*nhgc_pmts];
+  Double_t phgc_pulsePed[maxAdcHits*nhgc_pmts], phgc_pulseInt[maxAdcHits*nhgc_pmts];
+  Double_t phgc_pulseAmp[maxAdcHits*nhgc_pmts], phgc_pulseTime[maxAdcHits*nhgc_pmts];
+
+  // Noble gas Cherenkov ADC variables
+  Int_t pngc_hits;
+  Double_t pngc_pmt[maxAdcHits*nngc_pmts], pngc_pulsePedRaw[maxAdcHits*nngc_pmts];
+  Double_t pngc_pulseIntRaw[maxAdcHits*nngc_pmts], pngc_pulseAmpRaw[maxAdcHits*nngc_pmts];
+  Double_t pngc_pulsePed[maxAdcHits*nngc_pmts], pngc_pulseInt[maxAdcHits*nngc_pmts];
+  Double_t pngc_pulseAmp[maxAdcHits*nngc_pmts], pngc_pulseTime[maxAdcHits*nngc_pmts];
+
+  // Trigger time
+  Double_t pT1_tdcTime;
+
+  // Declare Histos
+  TH2F *h2_phgc_pulsePedRaw_vs_pmt, *h2_phgc_pulseIntRaw_vs_pmt, *h2_phgc_pulseAmpRaw_vs_pmt;
+  TH2F *h2_phgc_pulsePed_vs_pmt, *h2_phgc_pulseInt_vs_pmt, *h2_phgc_pulseAmp_vs_pmt, *h2_phgc_pulseTimeCorr_vs_pmt;
+
+  TH2F *h2_pngc_pulsePedRaw_vs_pmt, *h2_pngc_pulseIntRaw_vs_pmt, *h2_pngc_pulseAmpRaw_vs_pmt;
+  TH2F *h2_pngc_pulsePed_vs_pmt, *h2_pngc_pulseInt_vs_pmt, *h2_pngc_pulseAmp_vs_pmt, *h2_pngc_pulseTimeCorr_vs_pmt;
+
+  // Declare ROOT files
+  TFile *rif = new TFile("ROOTfiles/shms_replay_488_-1.root");
+  TFile *rof = new TFile("cer_ped.root", "RECREATE");
+
+  // Declare tree
+  TTree *T = rif->Get("T");
+  Long64_t nentries; 
+    
+  // Acquire the number of entries
+  //nentries = T->GetEntries();
+  nentries = 100000; 
+
+  // Acquire the branches
+  T->SetBranchAddress("Ndata.P.hgcer.adcCounter", &phgc_hits);
+  T->SetBranchAddress("P.hgcer.adcCounter", phgc_pmt);
+  T->SetBranchAddress("P.hgcer.adcPedRaw", phgc_pulsePedRaw);
+  T->SetBranchAddress("P.hgcer.adcPed", phgc_pulsePed);
+  T->SetBranchAddress("P.hgcer.adcPulseIntRaw", phgc_pulseIntRaw);
+  T->SetBranchAddress("P.hgcer.adcPulseInt", phgc_pulseInt);
+  T->SetBranchAddress("P.hgcer.adcPulseAmpRaw", phgc_pulseAmpRaw);
+  T->SetBranchAddress("P.hgcer.adcPulseAmp", phgc_pulseAmp);
+  T->SetBranchAddress("P.hgcer.adcPulseTimeRaw", phgc_pulseTime);
+
+  T->SetBranchAddress("Ndata.P.ngcer.adcCounter", &pngc_hits);
+  T->SetBranchAddress("P.ngcer.adcCounter", pngc_pmt);
+  T->SetBranchAddress("P.ngcer.adcPedRaw", pngc_pulsePedRaw);
+  T->SetBranchAddress("P.ngcer.adcPed", pngc_pulsePed);
+  T->SetBranchAddress("P.ngcer.adcPulseIntRaw", pngc_pulseIntRaw);
+  T->SetBranchAddress("P.ngcer.adcPulseInt", pngc_pulseInt);
+  T->SetBranchAddress("P.ngcer.adcPulseAmpRaw", pngc_pulseAmpRaw);
+  T->SetBranchAddress("P.ngcer.adcPulseAmp", pngc_pulseAmp);
+  T->SetBranchAddress("P.ngcer.adcPulseTimeRaw", pngc_pulseTime);
+
+  T->SetBranchAddress("T.shms.pT1_tdcTime", &pT1_tdcTime);
+
+  // Declare histos
+  h2_phgc_pulsePedRaw_vs_pmt = new TH2F("h2_phgc_pulsePedRaw_vs_pmt", "SHMS HGC Raw Pulse Pedestal; PMT Number; Raw Pulse Pedestal / 1 ADC Count", 4, 0.5, 4.5, 8000, -4000, 4000);
+  h2_phgc_pulseIntRaw_vs_pmt = new TH2F("h2_phgc_pulseIntRaw_vs_pmt", "SHMS HGC Raw Pulse Integral; PMT Number; Raw Pulse Integral / 10 ADC Count", 4, 0.5, 4.5, 8000, -40000, 40000);
+  h2_phgc_pulseAmpRaw_vs_pmt = new TH2F("h2_phgc_pulseAmpRaw_vs_pmt", "SHMS HGC Raw Pulse Amplitude; PMT Number; Raw Pulse Amplitude / 1 ADC Count", 4, 0.5, 4.5, 8100, -4100, 4100);
+
+  h2_phgc_pulsePed_vs_pmt = new TH2F("h2_phgc_pulsePed_vs_pmt", "SHMS HGC Pulse Pedestal; PMT Number; Pulse Pedestal / 1 ADC Count", 4, 0.5, 4.5, 8000, -4000, 4000);
+  h2_phgc_pulseInt_vs_pmt = new TH2F("h2_phgc_pulseInt_vs_pmt", "SHMS HGC Pulse Integral; PMT Number; Pulse Integral / 10 ADC Count", 4, 0.5, 4.5, 8000, -40000, 40000);
+  h2_phgc_pulseAmp_vs_pmt = new TH2F("h2_phgc_pulseAmp_vs_pmt", "SHMS HGC Pulse Amplitude; PMT Number; Pulse Amplitude / 1 ADC Count", 4, 0.5, 4.5, 8100, -4100, 4100);
+  h2_phgc_pulseTimeCorr_vs_pmt = new TH2F("h2_phgc_pulseTimeCorr_vs_pmt", "SHMS HGC Pulse Time Corrected; PMT Number; Pulse Time / 1 ns", 4, 0.5, 4.5, 8100, -4100, 4100);
+
+  h2_pngc_pulsePedRaw_vs_pmt = new TH2F("h2_pngc_pulsePedRaw_vs_pmt", "SHMS NGC Raw Pulse Pedestal; PMT Number; Raw Pulse Pedestal / 1 ADC Count", 4, 0.5, 4.5, 8000, -4000, 4000);
+  h2_pngc_pulseIntRaw_vs_pmt = new TH2F("h2_pngc_pulseIntRaw_vs_pmt", "SHMS NGC Raw Pulse Integral; PMT Number; Raw Pulse Integral / 10 ADC Count", 4, 0.5, 4.5, 8000, -40000, 40000);
+  h2_pngc_pulseAmpRaw_vs_pmt = new TH2F("h2_pngc_pulseAmpRaw_vs_pmt", "SHMS NGC Raw Pulse Amplitude; PMT Number; Raw Pulse Amplitude / 1 ADC Count", 4, 0.5, 4.5, 8100, -4100, 4100);
+
+  h2_pngc_pulsePed_vs_pmt = new TH2F("h2_pngc_pulsePed_vs_pmt", "SHMS NGC Pulse Pedestal; PMT Number; Pulse Pedestal / 1 ADC Count", 4, 0.5, 4.5, 8000, -4000, 4000);
+  h2_pngc_pulseInt_vs_pmt = new TH2F("h2_pngc_pulseInt_vs_pmt", "SHMS NGC Pulse Integral; PMT Number; Pulse Integral / 10 ADC Count", 4, 0.5, 4.5, 8000, -40000, 40000);
+  h2_pngc_pulseAmp_vs_pmt = new TH2F("h2_pngc_pulseAmp_vs_pmt", "SHMS NGC Pulse Amplitude; PMT Number; Pulse Amplitude / 1 ADC Count", 4, 0.5, 4.5, 8100, -4100, 4100);
+  h2_pngc_pulseTimeCorr_vs_pmt = new TH2F("h2_pngc_pulseTimeCorr_vs_pmt", "SHMS NGC Pulse Time Corrected; PMT Number; Pulse Time / 1 ns", 4, 0.5, 4.5, 8100, -4100, 4100);
+
+    // Loop of entries in tree
+  for(UInt_t ievent = 0; ievent < nentries; ievent++) {
+
+    T->GetEntry(ievent);
+    if ((ievent)%1000 == 0) cout << "ievent = " << ievent << endl;
+
+    for (UInt_t ihgchit = 0; ihgchit < phgc_hits; ihgchit++) {
+
+      if (phgc_pulseAmpRaw[ihgchit] != 0.0) continue;
+      
+      h2_phgc_pulsePedRaw_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulsePedRaw[ihgchit]);
+      h2_phgc_pulseIntRaw_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulseIntRaw[ihgchit]);      
+      h2_phgc_pulseAmpRaw_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulseAmpRaw[ihgchit]);
+
+      h2_phgc_pulsePed_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulsePed[ihgchit]);
+      h2_phgc_pulseInt_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulseInt[ihgchit]);       
+      h2_phgc_pulseAmp_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulseAmp[ihgchit]);
+
+      h2_phgc_pulseTimeCorr_vs_pmt->Fill(phgc_pmt[ihgchit], phgc_pulseTime[ihgchit]*clk2adc - pT1_tdcTime*clk2tdc);
+
+    } // HGC loop
+
+    for (UInt_t ingchit = 0; ingchit < pngc_hits; ingchit++) {
+
+      if (pngc_pulseAmpRaw[ingchit] != 0.0) continue;
+      
+      h2_pngc_pulsePedRaw_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulsePedRaw[ingchit]);
+      h2_pngc_pulseIntRaw_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulseIntRaw[ingchit]);      
+      h2_pngc_pulseAmpRaw_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulseAmpRaw[ingchit]);
+
+      h2_pngc_pulsePed_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulsePed[ingchit]);
+      h2_pngc_pulseInt_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulseInt[ingchit]);      
+      h2_pngc_pulseAmp_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulseAmp[ingchit]);
+
+      h2_pngc_pulseTimeCorr_vs_pmt->Fill(pngc_pmt[ingchit], pngc_pulseTime[ingchit]*clk2adc - pT1_tdcTime*clk2tdc);
+
+    } // NGC loop
+
+  }  // Entry loop
+
+  rof->Write();
+
+}  // End of cer_ped()
diff --git a/macros/cer_ped.root b/macros/cer_ped.root
new file mode 100644
index 0000000000000000000000000000000000000000..8d38f0151dadac1c4349977750afbd6d2db25cc2
Binary files /dev/null and b/macros/cer_ped.root differ
diff --git a/macros/cer_ped_nocut.root b/macros/cer_ped_nocut.root
new file mode 100644
index 0000000000000000000000000000000000000000..16b8d990e3990a86ff83bba32a4676a5d12dba88
Binary files /dev/null and b/macros/cer_ped_nocut.root differ
diff --git a/macros/kpp_plots.C b/macros/kpp_plots.C
new file mode 100644
index 0000000000000000000000000000000000000000..42ea8f1f9effc02e88082b323353c8cfb122d822
--- /dev/null
+++ b/macros/kpp_plots.C
@@ -0,0 +1,59 @@
+#include <TH1.h>
+#include <TFile.h>
+
+void kpp_plots() {
+
+  TFile *run_487_kpp_df = new TFile("./hallc-online_HISTOGRAMS/kpp_487.root");
+
+  TH1D *h_p1X_plTime = new TH1D();
+  TH1D *h_p1Y_plTime = new TH1D();
+  TH1D *h_p2X_plTime = new TH1D();
+  TH1D *h_p2Y_plTime = new TH1D();
+
+  TH2F *h2_hgc_pulseTime = new TH2F();
+  TH1D *h_hgc_pulseTime = new TH1D();
+
+  h_p1X_plTime = dynamic_cast <TH1D*> (run_487_kpp_df->Get("h_p1X_plTime"));
+  h_p1Y_plTime = dynamic_cast <TH1D*> (run_487_kpp_df->Get("h_p1Y_plTime"));
+  h_p2X_plTime = dynamic_cast <TH1D*> (run_487_kpp_df->Get("h_p2X_plTime"));
+  h_p2Y_plTime = dynamic_cast <TH1D*> (run_487_kpp_df->Get("h_p2Y_plTime"));
+
+  h2_hgc_pulseTime = dynamic_cast <TH2F*> (run_487_kpp_df->Get("h2_pngc_pulseTime_pT1_diff"));
+  h_hgc_pulseTime = h2_hgc_pulseTime->ProjectionY("h_hgc_pulseTime", 3, 3);
+  
+  h_p1X_plTime->SetLineColor(8);
+  h_p1X_plTime->SetFillColor(8);
+  h_p1X_plTime->Rebin(2);
+  h_p1X_plTime->GetXaxis()->SetRangeUser(20,60);
+  h_p1X_plTime->GetYaxis()->SetTitle("Counts / 2 ns");
+  		     
+  h_p1Y_plTime->SetLineColor(50);
+  h_p1Y_plTime->SetFillColor(50);
+  h_p1Y_plTime->Rebin(2);
+  h_p1Y_plTime->GetXaxis()->SetRangeUser(20,60);
+  h_p1Y_plTime->GetYaxis()->SetTitle("Counts / 2 ns");
+  
+  h_p2X_plTime->SetLineColor(9);
+  h_p2X_plTime->SetFillColor(9);
+  h_p2X_plTime->Rebin(2);
+  h_p2X_plTime->GetXaxis()->SetRangeUser(20,60);
+  h_p2X_plTime->GetYaxis()->SetTitle("Counts / 2 ns");
+  
+  h_p2Y_plTime->SetLineColor(52);
+  h_p2Y_plTime->SetFillColor(52);
+  h_p2Y_plTime->Rebin(2);
+  h_p2Y_plTime->GetXaxis()->SetRangeUser(20,60);
+  h_p2Y_plTime->GetYaxis()->SetTitle("Counts / 2 ns");
+  
+  TCanvas *c_hodo_plTime = new TCanvas("c_hodo_plTime", "", 800, 800);
+  c_hodo_plTime->Divide(2, 2);
+  c_hodo_plTime->cd(1);
+  h_p1X_plTime->Draw();
+  c_hodo_plTime->cd(2);
+  h_p1Y_plTime->Draw();
+  c_hodo_plTime->cd(3);
+  h_p2X_plTime->Draw();
+  c_hodo_plTime->cd(4);
+  h_p2Y_plTime->Draw();
+
+}
diff --git a/macros/ngcGain.C b/macros/ngcGain.C
new file mode 100644
index 0000000000000000000000000000000000000000..02dd3e5e56c3e558903481f538ef785b5537543c
--- /dev/null
+++ b/macros/ngcGain.C
@@ -0,0 +1,79 @@
+{
+
+  TCanvas *can = new TCanvas("can", "", 1800, 800);
+  can->Divide(2, 2);
+ 
+ TCanvas *can2 = new TCanvas("can2", "", 1800, 800);
+  can2->Divide(2, 2);
+
+ TCanvas *can3 = new TCanvas("can3", "", 1800, 800);
+  can3->Divide(2, 2);
+
+ TCanvas *can4 = new TCanvas("can4", "", 1800, 800);
+  can4->Divide(2, 2);
+ TCanvas *can5 = new TCanvas("can5", "", 1800, 800);
+  can5->Divide(2, 2);
+  
+  can->cd(1);
+  T->Draw("P.ngcer.goodAdcPulseInt[0]>>pmt1(100, 0, 200)", "P.ngcer.goodAdcPulseInt[0]>0.0&&P.ngcer.numTracksMatched[0]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[0]>-40.&&P.ngcer.goodAdcTdcDiffTime[0]<-10.");
+
+  can->cd(2);
+  T->Draw("P.ngcer.goodAdcPulseInt[1]>>pmt2(100, 0, 200)", "P.ngcer.goodAdcPulseInt[1]>0.0&&P.ngcer.numTracksMatched[1]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[1]>-40.&&P.ngcer.goodAdcTdcDiffTime[1]<-10.");
+
+  can->cd(3);
+  T->Draw("P.ngcer.goodAdcPulseInt[2]>>pmt3(100, 0, 200)", "P.ngcer.goodAdcPulseInt[2]>0.0&&P.ngcer.numTracksMatched[2]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[2]>-40.&&P.ngcer.goodAdcTdcDiffTime[2]<-10.");
+
+  can->cd(4);
+  T->Draw("P.ngcer.goodAdcPulseInt[3]>>pmt4(100, 0, 200)", "P.ngcer.goodAdcPulseInt[3]>0.0&&P.ngcer.numTracksMatched[3]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcTdcDiffTime[3]>-40.&&P.ngcer.goodAdcTdcDiffTime[3]<0.");
+
+  can5->cd(1);
+  T->Draw("P.ngcer.npe[0]>>npe1(100, 0, 100)", "P.ngcer.goodAdcPulseInt[0]>0.0&&P.ngcer.numTracksMatched[0]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[0]>-40.&&P.ngcer.goodAdcTdcDiffTime[0]<-10.");
+
+  can5->cd(2);
+  T->Draw("P.ngcer.npe[1]>>npe2(100, 0, 100)", "P.ngcer.goodAdcPulseInt[1]>0.0&&P.ngcer.numTracksMatched[1]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[1]>-40.&&P.ngcer.goodAdcTdcDiffTime[1]<-10.");
+
+  can5->cd(3);
+  T->Draw("P.ngcer.npe[2]>>npe3(100, 0, 100)", "P.ngcer.goodAdcPulseInt[2]>0.0&&P.ngcer.numTracksMatched[2]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[3]<10.&&P.ngcer.goodAdcTdcDiffTime[2]>-40.&&P.ngcer.goodAdcTdcDiffTime[2]<-10.");
+
+  can5->cd(4);
+  T->Draw("P.ngcer.npe[3]>>npe4(100, 0, 100)", "P.ngcer.goodAdcPulseInt[3]>0.0&&P.ngcer.numTracksMatched[3]>0.0&&P.cal.etracknorm>0.7&&P.ngcer.goodAdcPulseAmp[0]<10.&&P.ngcer.goodAdcPulseAmp[1]<10.&&P.ngcer.goodAdcPulseAmp[2]<10.&&P.ngcer.goodAdcTdcDiffTime[3]>-40.&&P.ngcer.goodAdcTdcDiffTime[3]<0.");
+  
+
+  can2->cd(1);
+  T->Draw("P.ngcer.goodAdcPulseAmp[0]>>pmt11(500, 0, 100)", "P.ngcer.goodAdcPulseInt[0]>0.0&&P.cal.etracknorm<0.4");
+
+  can2->cd(2);
+  T->Draw("P.ngcer.goodAdcPulseAmp[1]>>pmt21(500, 0, 100)", "P.ngcer.goodAdcPulseInt[1]>0.0&&P.cal.etracknorm<0.4");
+
+  can2->cd(3);
+  T->Draw("P.ngcer.goodAdcPulseAmp[2]>>pmt31(500, 0, 100)", "P.ngcer.goodAdcPulseInt[2]>0.0&&P.cal.etracknorm<0.4");
+
+  can2->cd(4);
+  T->Draw("P.ngcer.goodAdcPulseAmp[3]>>pmt41(500, 0, 100)", "P.ngcer.goodAdcPulseInt[3]>0.0&&P.cal.etracknorm<0.4");
+
+
+  //Raw TDC Times
+  can3->cd(1);
+  T->Draw("P.ngcer.goodAdcTdcDiffTime[0]>>pmt111(200, -100, 100)", "P.ngcer.goodAdcPulseInt[0]>0.0&&P.ngcer.numTracksMatched[1]>0.0&&P.cal.etracknorm>0.7");
+
+  can3->cd(2);
+  T->Draw("P.ngcer.goodAdcTdcDiffTime[1]>>pmt211(200, -100, 100)", "P.ngcer.goodAdcPulseInt[1]>0.0&&P.ngcer.numTracksMatched[0]>0.0&&P.cal.etracknorm>0.7");
+
+  can3->cd(3);
+  T->Draw("P.ngcer.goodAdcTdcDiffTime[2]>>pmt311(200, -100, 100)", "P.ngcer.goodAdcPulseInt[2]>0.0&&P.ngcer.numTracksMatched[3]>0.0&&P.cal.etracknorm>0.7");
+
+  can3->cd(4);
+  T->Draw("P.ngcer.goodAdcTdcDiffTime[3]>>pmt411(200, -100, 100)", "P.ngcer.goodAdcPulseInt[3]>0.0&&P.ngcer.numTracksMatched[2]>0.0&&P.cal.etracknorm>0.7");
+
+  can4->cd(1);
+  T->Draw("P.ngcer.goodAdcPulseAmp[0]:10*T.shms.pNGCER_adcPulseAmp","T.shms.pNGCER_adcPulseAmp>0.&&P.ngcer.goodAdcPulseAmp[0]>0.0&&P.ngcer.numTracksMatched[0]>0.0&&P.cal.etracknorm>0.7","colz");
+  can4->cd(2);
+  T->Draw("P.ngcer.goodAdcPulseAmp[1]:10*T.shms.pNGCER_adcPulseAmp","T.shms.pNGCER_adcPulseAmp>0.&&P.ngcer.goodAdcPulseAmp[1]>0.0&&P.ngcer.numTracksMatched[1]>0.0&&P.cal.etracknorm>0.7","colz");
+  can4->cd(3);
+  T->Draw("P.ngcer.goodAdcPulseAmp[2]:10*T.shms.pNGCER_adcPulseAmp","T.shms.pNGCER_adcPulseAmp>0.&&P.ngcer.goodAdcPulseAmp[2]>0.0&&P.ngcer.numTracksMatched[2]>0.0&&P.cal.etracknorm>0.7","colz");
+  can4->cd(4);
+  T->Draw("P.ngcer.goodAdcPulseAmp[3]:10*T.shms.pNGCER_adcPulseAmp","T.shms.pNGCER_adcPulseAmp>0.&&P.ngcer.goodAdcPulseAmp[3]>0.0&&P.ngcer.numTracksMatched[3]>0.0&&P.cal.etracknorm>0.7","colz");
+
+  // T.ngcer_adcPulseTime
+
+}
diff --git a/macros/ngcGainMatch.C b/macros/ngcGainMatch.C
new file mode 100644
index 0000000000000000000000000000000000000000..17e6cdd31f7e0937ef3d4a43e78931f8079ff8cb
--- /dev/null
+++ b/macros/ngcGainMatch.C
@@ -0,0 +1,76 @@
+#define nruns 2
+#define ntubes 4
+{
+
+ 
+  //Int_t nruns  = 1;
+  //Int_t ntubes = 3;
+  Double_t tube1_mean[nruns] = {249.3, 305.4};
+  Double_t tube2_mean[nruns] = {197.6, 257.7};
+  Double_t tube3_mean[nruns] = {198.4, 220.2};
+  Double_t tube4_mean[nruns] = {206.0, 295.1};
+
+  Double_t tube1_mean_err[nruns] = {1.4, 4.1};
+  Double_t tube2_mean_err[nruns] = {0.9, 1.6};
+  Double_t tube3_mean_err[nruns] = {0.6, 0.9};
+  Double_t tube4_mean_err[nruns] = {1.0, 1.9};
+
+  Double_t tube1_hv[nruns] = {1950, 2000}; 
+  Double_t tube2_hv[nruns] = {2150, 2200};
+  Double_t tube3_hv[nruns] = {1965, 1990};
+  Double_t tube4_hv[nruns] = {1887, 1917};
+
+  Double_t zero_err[ntubes] = {0.0, 0.0, 0.0, 0.0};
+
+  TF1 *tube1_fit = new TF1("tube1_fit", "pol1", 2000, 2400);
+  TF1 *tube2_fit = new TF1("tube2_fit", "pol1", 2000, 2400);
+  TF1 *tube3_fit = new TF1("tube3_fit", "pol1", 2000, 2400);
+  TF1 *tube4_fit = new TF1("tube4_fit", "pol1", 2000, 2400);
+
+  tube1_fit->SetLineColor(2);
+  tube2_fit->SetLineColor(2);
+  tube3_fit->SetLineColor(2);  
+  tube4_fit->SetLineColor(2);
+
+  TGraphErrors *tube1_gr = new TGraphErrors(nruns, tube1_hv, tube1_mean, zero_err, tube1_mean_err);
+  tube1_gr->SetMarkerSize(1);
+  tube1_gr->SetMarkerColor(4);
+  tube1_gr->SetMarkerStyle(21);
+  tube1_gr->SetTitle("PMT 1");
+  //tube1_gr->Fit(tube1_fit, "R");
+  TGraphErrors *tube2_gr = new TGraphErrors(nruns, tube2_hv, tube2_mean, zero_err, tube2_mean_err);
+  tube2_gr->SetMarkerSize(1);
+  tube2_gr->SetMarkerColor(4);
+  tube2_gr->SetMarkerStyle(21);
+  tube2_gr->SetTitle("PMT 2");
+  //tube2_gr->Fit(tube2_fit, "R");
+  TGraphErrors *tube3_gr = new TGraphErrors(nruns, tube3_hv, tube3_mean, zero_err, tube3_mean_err);
+  tube3_gr->SetMarkerSize(1);
+  tube3_gr->SetMarkerColor(4);
+  tube3_gr->SetMarkerStyle(21);
+  tube3_gr->SetTitle("PMT 3");
+  //tube3_gr->Fit(tube3_fit, "R");
+  TGraphErrors *tube4_gr = new TGraphErrors(nruns, tube4_hv, tube4_mean, zero_err, tube4_mean_err);
+  tube4_gr->SetMarkerSize(1);
+  tube4_gr->SetMarkerColor(4);
+  tube4_gr->SetMarkerStyle(21);
+  tube4_gr->SetTitle("PMT 4");
+  //tube4_gr->Fit(tube4_fit, "R");
+
+  TCanvas *can = new TCanvas("can", "can", 1800, 600);
+  can->Divide(2, 2);
+  can->cd(1);
+  //tube1_gr->Fit(tube1_fit, "R");
+  tube1_gr->Draw("AP");
+  can->cd(2);
+  tube2_gr->Draw("AP");
+  //tube2_gr->Fit(tube2_fit, "R");
+  can->cd(3);
+  tube3_gr->Draw("AP");
+  //tube3_gr->Fit(tube3_fit, "R");
+  can->cd(4);
+  tube4_gr->Draw("AP");
+  //tube4_gr->Fit(tube4_fit, "R");
+
+
+}
diff --git a/macros/ngc_gain_KPP.C b/macros/ngc_gain_KPP.C
new file mode 100644
index 0000000000000000000000000000000000000000..8f0c2d6204d34e3d2ddd1f1ae8b0c4735064bb1d
--- /dev/null
+++ b/macros/ngc_gain_KPP.C
@@ -0,0 +1,71 @@
+#include "TMath.h"
+
+void ngc_gain() {
+
+  static const Int_t ntubes = 4;
+
+  Double_t tube1_mean[ntubes] = {1.68060e+03, 6.23210e+03, 3.41938e+03, 8.28936e+03};
+  Double_t tube2_mean[ntubes] = {9.45180e+02, 4.93875e+03, 2.34050e+03, 7.64824e+03};
+  Double_t tube3_mean[ntubes] = {1.44629e+03, 7.79568e+03, 3.51568e+03, 1.16708e+04};
+  Double_t tube4_mean[ntubes] = {2.06980e+03, 1.03371e+04, 3.74562e+03, 1.49880e+04};
+
+  Double_t tube1_mean_err[ntubes] = {1.30291e+01, 4.75684e+01, 2.56718e+01, 4.84467e+01};
+  Double_t tube2_mean_err[ntubes] = {5.42784e+00, 2.24776e+01, 1.51154e+01, 3.22023e+01};
+  Double_t tube3_mean_err[ntubes] = {6.76353e+00, 3.58100e+01, 1.72918e+01, 3.75558e+01};
+  Double_t tube4_mean_err[ntubes] = {2.46633e+01, 8.98439e+01, 7.70273e+01, 7.62538e+01};
+
+  Double_t tube1_hv[ntubes] = {1900.0, 2200.0, 2050.0, 2300.0}; 
+  Double_t tube2_hv[ntubes] = {1900.0, 2200.0, 2050.0, 2300.0};
+  Double_t tube3_hv[ntubes] = {1900.0, 2200.0, 2050.0, 2300.0};
+  Double_t tube4_hv[ntubes] = {1900.0, 2200.0, 2000.0, 2300.0};
+
+  Double_t zero_err[ntubes] = {0.0, 0.0, 0.0, 0.0};
+
+  TF1 *tube1_fit = new TF1("tube1_fit", "pol1", 1800, 2400);
+  TF1 *tube2_fit = new TF1("tube2_fit", "pol1", 1800, 2400);
+  TF1 *tube3_fit = new TF1("tube3_fit", "pol1", 1800, 2400);
+  TF1 *tube4_fit = new TF1("tube4_fit", "pol1", 1800, 2400);
+
+  tube1_fit->SetLineColor(2);
+  tube2_fit->SetLineColor(2);
+  tube3_fit->SetLineColor(2);  
+  tube4_fit->SetLineColor(2);
+
+  TGraphErrors *tube1_gr = new TGraphErrors(ntubes, tube1_hv, tube1_mean, zero_err, tube1_mean_err);
+  tube1_gr->SetMarkerSize(1);
+  tube1_gr->SetMarkerColor(4);
+  tube1_gr->SetMarkerStyle(21);
+  tube1_gr->SetTitle("PMT 1");
+  //tube1_gr->Fit(tube1_fit, "R");
+  TGraphErrors *tube2_gr = new TGraphErrors(ntubes, tube2_hv, tube2_mean, zero_err, tube2_mean_err);
+  tube2_gr->SetMarkerSize(1);
+  tube2_gr->SetMarkerColor(4);
+  tube2_gr->SetMarkerStyle(21);
+  tube2_gr->SetTitle("PMT 2");
+  //tube2_gr->Fit(tube2_fit, "R");
+  TGraphErrors *tube3_gr = new TGraphErrors(ntubes, tube3_hv, tube3_mean, zero_err, tube3_mean_err);
+  tube3_gr->SetMarkerSize(1);
+  tube3_gr->SetMarkerColor(4);
+  tube3_gr->SetMarkerStyle(21);
+  tube3_gr->SetTitle("PMT 3");
+  //tube3_gr->Fit(tube3_fit, "R");
+  TGraphErrors *tube4_gr = new TGraphErrors(ntubes, tube4_hv, tube4_mean, zero_err, tube4_mean_err);
+  tube4_gr->SetMarkerSize(1);
+  tube4_gr->SetMarkerColor(4);
+  tube4_gr->SetMarkerStyle(21);
+  tube4_gr->SetTitle("PMT 4");
+  //tube4_gr->Fit(tube4_fit, "R");
+
+  TCanvas *can = new TCanvas("can", "can", 800, 800);
+  can->Divide(2, 2);
+  can->cd(1);
+  tube1_gr->Draw("AP");
+  can->cd(2);
+  tube2_gr->Draw("AP");
+  can->cd(3);
+  tube3_gr->Draw("AP");
+  can->cd(4);
+  tube4_gr->Draw("AP");
+  
+
+}
diff --git a/macros/trig_eff_hms_hl.C b/macros/trig_eff_hms_hl.C
new file mode 100644
index 0000000000000000000000000000000000000000..cf5066ed5c17b83e534975785417cdd53398001d
--- /dev/null
+++ b/macros/trig_eff_hms_hl.C
@@ -0,0 +1,322 @@
+{
+  gStyle->SetOptStat(0);
+
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  
+  can1->Divide(3, 2);
+  can1->cd(1);
+  gPad->SetLogy();
+  T->Draw("H.cer.npeSum>>h1(100,0,50)", "H.cer.npeSum>0.0");
+  T->Draw("H.cer.npeSum>>h2(100,0,50)", "H.cer.npeSum>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npeSum>>h3(100,0,50)", "H.cer.npeSum>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  h3->SetLineColor(8);
+  h3->SetFillColor(8);
+  h3->SetFillStyle(3006);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum 3/4");
+  leg1->AddEntry(h2, "NPE Sum EL-REAL");
+  leg1->AddEntry(h3, "NPE Sum EL-CLEAN");
+  leg1->Draw();
+
+  can1->cd(4);
+  TH1F *h22 = (TH1F*)h1->Clone("h22");
+  h22->Sumw2();
+  h22->Divide(h2);
+  h22->SetLineColor(2);
+  h22->Draw();
+
+  TH1F *h23 = (TH1F*)h1->Clone("h23");
+  h23->Sumw2();
+  h23->Divide(h3);
+  h23->SetLineColor(4);
+  h23->Draw("SAME");
+
+  TH1F *h24 = (TH1F*)h2->Clone("h24");
+  h24->Sumw2();
+  h24->Divide(h3);
+  h24->SetLineColor(8);
+  h24->Draw("SAME");
+
+  TLegend *leg8 = new TLegend(0.7,0.7,0.9,0.9);
+  leg8->AddEntry(h22, "3/4 / EL-REAL", "lep");
+  leg8->AddEntry(h23, "3/4 / EL-CLEAN", "lep");
+  leg8->AddEntry(h24, "EL-REAL / EL-CLEAN", "lep");
+  leg8->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  gPad->SetLogy();
+  T->Draw("H.cer.npe[0]>>h4(100,0,50)", "H.cer.npe[0]>0.0");
+  T->Draw("H.cer.npe[0]>>h5(100,0,50)", "H.cer.npe[0]>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npe[0]>>h6(100,0,50)", "H.cer.npe[0]>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+
+  h6->SetLineColor(8);
+  h6->SetFillColor(8);
+  h6->SetFillStyle(3006);
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h4, "PMT 1 3/4");
+  leg2->AddEntry(h5, "PMT 1 EL-REAL");
+  leg2->AddEntry(h6, "PMT 1 EL-CLEAN");
+  leg2->Draw();
+
+  can1->cd(5);
+  TH1F *h25 = (TH1F*)h4->Clone("h25");
+  h25->Sumw2();
+  h25->Divide(h5);
+  h25->SetLineColor(2);
+  h25->Draw();
+
+  TH1F *h26 = (TH1F*)h4->Clone("h26");
+  h26->Sumw2();
+  h26->Divide(h6);
+  h26->SetLineColor(4);
+  h26->Draw("SAME");
+
+  TH1F *h27 = (TH1F*)h5->Clone("h27");
+  h27->Sumw2();
+  h27->Divide(h6);
+  h27->SetLineColor(8);
+  h27->Draw("SAME");
+
+  TLegend *leg9 = new TLegend(0.7,0.7,0.9,0.9);
+  leg9->AddEntry(h25, "3/4 / EL-REAL", "lep");
+  leg9->AddEntry(h26, "3/4 / EL-CLEAN", "lep");
+  leg9->AddEntry(h27, "EL-REAL / EL-CLEAN", "lep");
+  leg9->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  gPad->SetLogy();
+  T->Draw("H.cer.npe[1]>>h7(100,0,50)", "H.cer.npe[1]>0.0");
+  T->Draw("H.cer.npe[1]>>h8(100,0,50)", "H.cer.npe[1]>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npe[1]>>h9(100,0,50)", "H.cer.npe[1]>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+
+  h9->SetLineColor(8);
+  h9->SetFillColor(8);
+  h9->SetFillStyle(3006);
+
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h7, "PMT 2 3/4");
+  leg3->AddEntry(h8, "PMT 2 EL-REAL");
+  leg3->AddEntry(h9, "PMT 2 EL-CLEAN");
+  leg3->Draw();
+
+  can1->cd(6);
+  TH1F *h28 = (TH1F*)h7->Clone("h28");
+  h28->Sumw2();
+  h28->Divide(h8);
+  h28->SetLineColor(2);
+  h28->Draw();
+
+  TH1F *h29 = (TH1F*)h7->Clone("h29");
+  h29->Sumw2();
+  h29->Divide(h9);
+  h29->SetLineColor(4);
+  h29->Draw("SAME");
+
+  TH1F *h30 = (TH1F*)h8->Clone("h30");
+  h30->Sumw2();
+  h30->Divide(h9);
+  h30->SetLineColor(8);
+  h30->Draw("SAME");
+
+  TLegend *leg11 = new TLegend(0.7,0.7,0.9,0.9);
+  leg11->AddEntry(h28, "3/4 / EL-REAL", "lep");
+  leg11->AddEntry(h29, "3/4 / EL-CLEAN", "lep");
+  leg11->AddEntry(h30, "EL-REAL / EL-CLEAN", "lep");
+  leg11->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+  // /-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+
+  TCanvas *can2 = new TCanvas("can2", "can2", 1600, 800);
+  
+  can2->Divide(3, 2);
+  can2->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.eprtracknorm>>h10(75, 0, 1.5)", "H.cal.eprtracknorm>0.0");
+  T->Draw("H.cal.eprtracknorm>>h11(75, 0, 1.5)", "H.cal.eprtracknorm>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.eprtracknorm>>h12(75, 0, 1.5)", "H.cal.eprtracknorm>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h10->SetLineColor(2);
+  h10->SetFillColor(2);
+  h10->SetFillStyle(3004);
+
+  h11->SetLineColor(4);
+  h11->SetFillColor(4);
+  h11->SetFillStyle(3005);
+
+  h12->SetLineColor(8);
+  h12->SetFillColor(8);
+  h12->SetFillStyle(3006);
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h10, "PrSh E-Norm 3/4");
+  leg4->AddEntry(h11, "PrSh E-Norm EL-REAL");
+  leg4->AddEntry(h12, "PrSh E-Norm EL-CLEAN");
+  leg4->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.etracknorm>>h13(75, 0, 1.5)", "H.cal.etracknorm>0.0");
+  T->Draw("H.cal.etracknorm>>h14(75, 0, 1.5)", "H.cal.etracknorm>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.etracknorm>>h15(75, 0, 1.5)", "H.cal.etracknorm>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h13->SetLineColor(2);
+  h13->SetFillColor(2);
+  h13->SetFillStyle(3004);
+
+  h14->SetLineColor(4);
+  h14->SetFillColor(4);
+  h14->SetFillStyle(3005);
+
+  h15->SetLineColor(8);
+  h15->SetFillColor(8);
+  h15->SetFillStyle(3006);
+
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h13, "Sh E-Norm 3/4");
+  leg5->AddEntry(h14, "Sh E-Norm EL-REAL");
+  leg5->AddEntry(h15, "Sh E-Norm EL-CLEAN");
+  leg5->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.etottracknorm>>h31(75, 0, 1.5)", "H.cal.etottracknorm>0.0");
+  T->Draw("H.cal.etottracknorm>>h32(75, 0, 1.5)", "H.cal.etottracknorm>0.0&&T.hms.hEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.etottracknorm>>h33(75, 0, 1.5)", "H.cal.etottracknorm>0.0&&T.hms.hEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h31->SetLineColor(2);
+  h31->SetFillColor(2);
+  h31->SetFillStyle(3004);
+
+  h32->SetLineColor(4);
+  h32->SetFillColor(4);
+  h32->SetFillStyle(3005);
+
+  h33->SetLineColor(8);
+  h33->SetFillColor(8);
+  h33->SetFillStyle(3006);
+
+  TLegend *leg12 = new TLegend(0.7,0.7,0.9,0.9);
+  leg12->AddEntry(h31, "Total E-Norm 3/4");
+  leg12->AddEntry(h32, "Total E-Norm EL-REAL");
+  leg12->AddEntry(h33, "Total E-Norm EL-CLEAN");
+  leg12->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(4);
+  TH1F *h16 = (TH1F*)h10->Clone("h16");
+  h16->Sumw2();
+  h16->Divide(h11);
+  h16->SetLineColor(2);
+  h16->Draw();
+
+  TH1F *h17 = (TH1F*)h10->Clone("h17");
+  h17->Sumw2();
+  h17->Divide(h12);
+  h17->SetLineColor(4);
+  h17->Draw("SAME");
+
+  TH1F *h18 = (TH1F*)h11->Clone("h18");
+  h18->Sumw2();
+  h18->Divide(h12);
+  h18->SetLineColor(8);
+  h18->Draw("SAME");
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h16, "3/4 / EL-REAL", "lep");
+  leg6->AddEntry(h17, "3/4 / EL-CLEAN", "lep");
+  leg6->AddEntry(h18, "EL-REAL / EL-CLEAN", "lep");
+  leg6->Draw();
+
+  can2->cd(5);
+  TH1F *h19 = (TH1F*)h13->Clone("h19");
+  h19->Sumw2();
+  h19->Divide(h14);
+  h19->SetLineColor(2);
+  h19->Draw();
+
+  TH1F *h20 = (TH1F*)h13->Clone("h20");
+  h20->Sumw2();
+  h20->Divide(h15);
+  h20->SetLineColor(4);
+  h20->Draw("SAME");
+
+  TH1F *h21 = (TH1F*)h14->Clone("h21");
+  h21->Sumw2();
+  h21->Divide(h15);
+  h21->SetLineColor(8);
+  h21->Draw("SAME");
+
+  TLegend *leg7 = new TLegend(0.7,0.7,0.9,0.9);
+  leg7->AddEntry(h19, "3/4 / EL-REAL", "lep");
+  leg7->AddEntry(h20, "3/4 / EL-CLEAN", "lep");
+  leg7->AddEntry(h21, "EL-REAL / EL-CLEAN", "lep");
+  leg7->Draw();
+
+  can2->cd(6);
+  TH1F *h34 = (TH1F*)h31->Clone("h34");
+  h34->Sumw2();
+  h34->Divide(h32);
+  h34->SetLineColor(2);
+  h34->Draw();
+
+  TH1F *h35 = (TH1F*)h31->Clone("h35");
+  h35->Sumw2();
+  h35->Divide(h33);
+  h35->SetLineColor(4);
+  h35->Draw("SAME");
+
+  TH1F *h36 = (TH1F*)h32->Clone("h36");
+  h36->Sumw2();
+  h36->Divide(h34);
+  h36->SetLineColor(8);
+  h36->Draw("SAME");
+
+  TLegend *leg13 = new TLegend(0.7,0.7,0.9,0.9);
+  leg13->AddEntry(h34, "3/4 / EL-REAL", "lep");
+  leg13->AddEntry(h35, "3/4 / EL-CLEAN", "lep");
+  leg13->AddEntry(h36, "EL-REAL / EL-CLEAN", "lep");
+  leg13->Draw();
+  
+}
diff --git a/macros/trig_eff_hms_ll.C b/macros/trig_eff_hms_ll.C
new file mode 100644
index 0000000000000000000000000000000000000000..674c10140f86a89969f70469436f7380ea0ba0c5
--- /dev/null
+++ b/macros/trig_eff_hms_ll.C
@@ -0,0 +1,165 @@
+{
+  gStyle->SetOptStat(0);
+  
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  can1->Divide(4, 2);
+  
+  can1->cd(1);
+  //gPad->SetLogy();
+
+  T->Draw("H.cer.npeSum>>h1(100,0,25)", "H.cer.npeSum>0.0");
+  T->Draw("H.cer.npeSum>>h2(100,0,25)", "H.cer.npeSum>0.0&&T.hms.hCERSUM_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum");
+  leg1->AddEntry(h2, "CER");
+  leg1->Draw();
+
+  can1->cd(5);
+  TH1F *h3 = (TH1F*)h1->Clone("h3");
+  h3->Sumw2();
+  h3->Divide(h2);
+  h3->SetLineColor(9);
+  h3->SetMarkerStyle(21);
+  h3->SetMarkerSize(0.5);
+  h3->Draw();
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h3, "CER", "lep");
+  leg2->Draw();
+
+  TLine *line1 = new TLine(0.0, 1.0, 5.0, 1.0);
+  line1->SetLineWidth(3);
+  line1->SetLineStyle(2);
+  line1->SetLineColor(6);
+  line1->Draw();
+
+  //=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  //gPad->SetLogy();
+
+  T->Draw("H.cal.etottracknorm>>h4(75,0,1.5)", "H.cal.etottracknorm>0.0");
+  T->Draw("H.cal.etottracknorm>>h5(75,0,1.5)", "H.cal.etottracknorm>0.0&&T.hms.hPSHWRLO_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+  
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+  
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h4, "ETOT NORM");
+  leg3->AddEntry(h5, "PRSHWRLO");
+  leg3->Draw();
+
+  can1->cd(6);
+  TH1F *h6 = (TH1F*)h4->Clone("h6");
+  h6->Sumw2();
+  h6->Divide(h5);
+  h6->SetLineColor(9);
+  h6->SetMarkerStyle(21);
+  h6->SetMarkerSize(0.5);
+  h6->Draw();
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h6, "PRSHWRLO", "lep");
+  leg4->Draw();
+
+  TLine *line2 = new TLine(0.0, 1.0, 1.0, 1.0);
+  line2->SetLineWidth(3);
+  line2->SetLineStyle(2);
+  line2->SetLineColor(6);
+  line2->Draw();
+
+ //=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  //gPad->SetLogy();
+
+  T->Draw("H.cal.etottracknorm>>h7(75,0,1.5)", "H.cal.etottracknorm>0.0");
+  T->Draw("H.cal.etottracknorm>>h8(75,0,1.5)", "H.cal.etottracknorm>0.0&&T.hms.hPSHWRHI_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+  
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+  
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h7, "ETOT NORM");
+  leg5->AddEntry(h8, "PRSHWRHI");
+  leg5->Draw();
+
+  can1->cd(7);
+  TH1F *h9 = (TH1F*)h7->Clone("h9");
+  h9->Sumw2();
+  h9->Divide(h8);
+  h9->SetLineColor(9);
+  h9->SetMarkerStyle(21);
+  h9->SetMarkerSize(0.5);
+  h9->Draw();
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h9, "PRSHWRHI", "lep");
+  leg6->Draw();
+
+  TLine *line3 = new TLine(0.0, 1.0, 1.0, 1.0);
+  line3->SetLineWidth(3);
+  line3->SetLineStyle(2);
+  line3->SetLineColor(6);
+  line3->Draw();
+
+ //=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(4);
+  //gPad->SetLogy();
+
+  T->Draw("H.cal.etottracknorm>>h10(75,0,1.5)", "H.cal.etottracknorm>0.0");
+  T->Draw("H.cal.etottracknorm>>h11(75,0,1.5)", "H.cal.etottracknorm>0.0&&T.hms.hSHWR_tdcTimeRaw>0.0", "SAME");
+
+  h10->SetLineColor(2);
+  h10->SetFillColor(2);
+  h10->SetFillStyle(3004);
+  
+  h11->SetLineColor(4);
+  h11->SetFillColor(4);
+  h11->SetFillStyle(3005);
+  
+  TLegend *leg7 = new TLegend(0.7,0.7,0.9,0.9);
+  leg7->AddEntry(h10, "ETOT NORM");
+  leg7->AddEntry(h11, "SHWRLO");
+  leg7->Draw();
+
+  can1->cd(8);
+  TH1F *h12 = (TH1F*)h10->Clone("h12");
+  h12->Sumw2();
+  h12->Divide(h11);
+  h12->SetLineColor(9);
+  h12->SetMarkerStyle(21);
+  h12->SetMarkerSize(0.5);
+  h12->Draw();
+
+  TLegend *leg8 = new TLegend(0.7,0.7,0.9,0.9);
+  leg8->AddEntry(h12, "SHWRLO", "lep");
+  leg8->Draw();
+
+  TLine *line4 = new TLine(0.0, 1.0, 0.5, 1.0);
+  line4->SetLineWidth(3);
+  line4->SetLineStyle(2);
+  line4->SetLineColor(6);
+  line4->Draw();
+
+}
diff --git a/macros/trig_eff_hms_ml.C b/macros/trig_eff_hms_ml.C
new file mode 100644
index 0000000000000000000000000000000000000000..66ee267164bd59d2209500a31ce2059b99dab391
--- /dev/null
+++ b/macros/trig_eff_hms_ml.C
@@ -0,0 +1,325 @@
+{
+  gStyle->SetOptStat(0);
+
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  
+  can1->Divide(3, 2);
+  can1->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("H.cer.npeSum>>h1(100,0,50)", "H.cer.npeSum>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cer.npeSum>>h2(100,0,50)", "H.cer.npeSum>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npeSum>>h3(100,0,50)", "H.cer.npeSum>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  h3->SetLineColor(8);
+  h3->SetFillColor(8);
+  h3->SetFillStyle(3006);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum EL-LO");
+  leg1->AddEntry(h2, "NPE Sum EL-LO-CER");
+  leg1->AddEntry(h3, "NPE Sum EL-HI");
+  leg1->Draw();
+
+  can1->cd(4);
+  TH1F *h22 = (TH1F*)h1->Clone("h22");
+  h22->Sumw2();
+  h22->Divide(h2);
+  h22->SetLineColor(2);
+  h22->Draw();
+
+  TH1F *h23 = (TH1F*)h1->Clone("h23");
+  h23->Sumw2();
+  h23->Divide(h3);
+  h23->SetLineColor(4);
+  h23->Draw("SAME");
+
+  TH1F *h24 = (TH1F*)h2->Clone("h24");
+  h24->Sumw2();
+  h24->Divide(h3);
+  h24->SetLineColor(8);
+  h24->Draw("SAME");
+
+  TLegend *leg8 = new TLegend(0.7,0.7,0.9,0.9);
+  leg8->AddEntry(h22, "EL-LO / EL-LO-CER", "lep");
+  leg8->AddEntry(h23, "EL-LO / EL-HI", "lep");
+  leg8->AddEntry(h24, "EL-LO-CER / EL-HI", "lep");
+  leg8->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("H.cer.npe[0]>>h4(100,0,50)", "H.cer.npe[0]>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cer.npe[0]>>h5(100,0,50)", "H.cer.npe[0]>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npe[0]>>h6(100,0,50)", "H.cer.npe[0]>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+
+  h6->SetLineColor(8);
+  h6->SetFillColor(8);
+  h6->SetFillStyle(3006);
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h4, "PMT 1 EL-LO");
+  leg2->AddEntry(h5, "PMT 1 EL-LO-CER");
+  leg2->AddEntry(h6, "PMT 1 EL-HI");
+  leg2->Draw();
+
+  can1->cd(5);
+  TH1F *h25 = (TH1F*)h4->Clone("h25");
+  h25->Sumw2();
+  h25->Divide(h5);
+  h25->SetLineColor(2);
+  h25->Draw();
+
+  TH1F *h26 = (TH1F*)h4->Clone("h26");
+  h26->Sumw2();
+  h26->Divide(h6);
+  h26->SetLineColor(4);
+  h26->Draw("SAME");
+
+  TH1F *h27 = (TH1F*)h5->Clone("h27");
+  h27->Sumw2();
+  h27->Divide(h6);
+  h27->SetLineColor(8);
+  h27->Draw("SAME");
+
+  TLegend *leg9 = new TLegend(0.7,0.7,0.9,0.9);
+  leg9->AddEntry(h25, "EL-LO / EL-LO-CER", "lep");
+  leg9->AddEntry(h26, "EL-LO / EL-HI", "lep");
+  leg9->AddEntry(h27, "EL-LO-CER / EL-HI", "lep");
+  leg9->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("H.cer.npe[1]>>h7(100,0,50)", "H.cer.npe[1]>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cer.npe[1]>>h8(100,0,50)", "H.cer.npe[1]>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cer.npe[1]>>h9(100,0,50)", "H.cer.npe[1]>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+
+  h9->SetLineColor(8);
+  h9->SetFillColor(8);
+  h9->SetFillStyle(3006);
+
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h7, "PMT 2 EL-LO");
+  leg3->AddEntry(h8, "PMT 2 EL-LO-CER");
+  leg3->AddEntry(h9, "PMT 2 EL-HI");
+  leg3->Draw();
+
+  can1->cd(6);
+  TH1F *h28 = (TH1F*)h7->Clone("h28");
+  h28->Sumw2();
+  h28->Divide(h8);
+  h28->SetLineColor(2);
+  h28->Draw();
+
+  TH1F *h29 = (TH1F*)h7->Clone("h29");
+  h29->Sumw2();
+  h29->Divide(h9);
+  h29->SetLineColor(4);
+  h29->Draw("SAME");
+
+  TH1F *h30 = (TH1F*)h8->Clone("h30");
+  h30->Sumw2();
+  h30->Divide(h9);
+  h30->SetLineColor(8);
+  h30->Draw("SAME");
+
+  TLegend *leg11 = new TLegend(0.7,0.7,0.9,0.9);
+  leg11->AddEntry(h28, "EL-LO / EL-LO-CER", "lep");
+  leg11->AddEntry(h29, "EL-LO / EL-HI", "lep");
+  leg11->AddEntry(h30, "EL-LO-CER / EL-HI", "lep");
+  leg11->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+  // /-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+
+  TCanvas *can2 = new TCanvas("can2", "can2", 1600, 800);
+  
+  can2->Divide(3, 2);
+  can2->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.eprtracknorm>>h10(75, 0, 1.5)", "H.cal.eprtracknorm>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cal.eprtracknorm>>h11(75, 0, 1.5)", "H.cal.eprtracknorm>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.eprtracknorm>>h12(75, 0, 1.5)", "H.cal.eprtracknorm>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h10->SetLineColor(2);
+  h10->SetFillColor(2);
+  h10->SetFillStyle(3004);
+
+  h11->SetLineColor(4);
+  h11->SetFillColor(4);
+  h11->SetFillStyle(3005);
+
+  h12->SetLineColor(8);
+  h12->SetFillColor(8);
+  h12->SetFillStyle(3006);
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h10, "PrSh E-Norm EL-LO");
+  leg4->AddEntry(h11, "PrSh E-Norm EL-LO-CER");
+  leg4->AddEntry(h12, "PrSh E-Norm EL-HI");
+  leg4->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.etracknorm>>h13(75, 0, 1.5)", "H.cal.etracknorm>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cal.etracknorm>>h14(75, 0, 1.5)", "H.cal.etracknorm>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.etracknorm>>h15(75, 0, 1.5)", "H.cal.etracknorm>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h13->SetLineColor(2);
+  h13->SetFillColor(2);
+  h13->SetFillStyle(3004);
+
+  h14->SetLineColor(4);
+  h14->SetFillColor(4);
+  h14->SetFillStyle(3005);
+
+  h15->SetLineColor(8);
+  h15->SetFillColor(8);
+  h15->SetFillStyle(3006);
+
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h13, "Sh E-Norm EL-LO");
+  leg5->AddEntry(h14, "Sh E-Norm EL-LO-CER");
+  leg5->AddEntry(h15, "Sh E-Norm EL-HI");
+  leg5->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("H.cal.etottracknorm>>h31(75, 0, 1.5)", "H.cal.etottracknorm>0.0&&T.hms.hEL_LO_tdcTimeRaw>0.0");
+  T->Draw("H.cal.etottracknorm>>h32(75, 0, 1.5)", "H.cal.etottracknorm>0.0&&T.hms.hEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("H.cal.etottracknorm>>h33(75, 0, 1.5)", "H.cal.etottracknorm>0.0&&T.hms.hEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h31->SetLineColor(2);
+  h31->SetFillColor(2);
+  h31->SetFillStyle(3004);
+
+  h32->SetLineColor(4);
+  h32->SetFillColor(4);
+  h32->SetFillStyle(3005);
+
+  h33->SetLineColor(8);
+  h33->SetFillColor(8);
+  h33->SetFillStyle(3006);
+
+  TLegend *leg12 = new TLegend(0.7,0.7,0.9,0.9);
+  leg12->AddEntry(h31, "Total E-Norm EL-LO");
+  leg12->AddEntry(h32, "Total E-Norm EL-LO-CER");
+  leg12->AddEntry(h33, "Total E-Norm EL-HI");
+  leg12->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(4);
+  TH1F *h16 = (TH1F*)h10->Clone("h16");
+  h16->Sumw2();
+  h16->Divide(h11);
+  h16->SetLineColor(2);
+  h16->Draw();
+
+  TH1F *h17 = (TH1F*)h10->Clone("h17");
+  h17->Sumw2();
+  h17->Divide(h12);
+  h17->SetLineColor(4);
+  h17->Draw("SAME");
+
+  TH1F *h18 = (TH1F*)h11->Clone("h18");
+  h18->Sumw2();
+  h18->Divide(h12);
+  h18->SetLineColor(8);
+  h18->Draw("SAME");
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h16, "EL-LO / EL-LO-CER", "lep");
+  leg6->AddEntry(h17, "EL-LO / EL-HI", "lep");
+  leg6->AddEntry(h18, "EL-LO-CER / EL-HI", "lep");
+  leg6->Draw();
+
+  can2->cd(5);
+  TH1F *h19 = (TH1F*)h13->Clone("h19");
+  h19->Sumw2();
+  h19->Divide(h14);
+  h19->SetLineColor(2);
+  h19->Draw();
+
+  TH1F *h20 = (TH1F*)h13->Clone("h20");
+  h20->Sumw2();
+  h20->Divide(h15);
+  h20->SetLineColor(4);
+  h20->Draw("SAME");
+
+  TH1F *h21 = (TH1F*)h14->Clone("h21");
+  h21->Sumw2();
+  h21->Divide(h15);
+  h21->SetLineColor(8);
+  h21->Draw("SAME");
+
+  TLegend *leg7 = new TLegend(0.7,0.7,0.9,0.9);
+  leg7->AddEntry(h19, "EL-LO / EL-LO-CER", "lep");
+  leg7->AddEntry(h20, "EL-LO / EL-HI", "lep");
+  leg7->AddEntry(h21, "EL-LO-CER / EL-HI", "lep");
+  leg7->Draw();
+
+  can2->cd(6);
+  TH1F *h34 = (TH1F*)h31->Clone("h34");
+  h34->Sumw2();
+  h34->Divide(h32);
+  h34->SetLineColor(2);
+  h34->Draw();
+
+  TH1F *h35 = (TH1F*)h31->Clone("h35");
+  h35->Sumw2();
+  h35->Divide(h33);
+  h35->SetLineColor(4);
+  h35->Draw("SAME");
+
+  TH1F *h36 = (TH1F*)h32->Clone("h36");
+  h36->Sumw2();
+  h36->Divide(h34);
+  h36->SetLineColor(8);
+  h36->Draw("SAME");
+
+  TLegend *leg13 = new TLegend(0.7,0.7,0.9,0.9);
+  leg13->AddEntry(h34, "EL-LO / EL-LO-CER", "lep");
+  leg13->AddEntry(h35, "EL-LO / EL-HI", "lep");
+  leg13->AddEntry(h36, "EL-LO-CER / EL-HI", "lep");
+  leg13->Draw();
+
+}
diff --git a/macros/trig_eff_shms_hl.C b/macros/trig_eff_shms_hl.C
new file mode 100644
index 0000000000000000000000000000000000000000..91200e8900dd91cf629c830764f7c9bfa92e9096
--- /dev/null
+++ b/macros/trig_eff_shms_hl.C
@@ -0,0 +1,322 @@
+{
+  gStyle->SetOptStat(0);
+
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  
+  can1->Divide(3, 2);
+  can1->cd(1);
+  gPad->SetLogy();
+  T->Draw("P.hgcer.npeSum>>h1(100,0,50)", "P.hgcer.npeSum>0.0");
+  T->Draw("P.hgcer.npeSum>>h2(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npeSum>>h3(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  h3->SetLineColor(8);
+  h3->SetFillColor(8);
+  h3->SetFillStyle(3006);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum 3/4");
+  leg1->AddEntry(h2, "NPE Sum EL-REAL");
+  leg1->AddEntry(h3, "NPE Sum EL-CLEAN");
+  leg1->Draw();
+
+  can1->cd(4);
+  TH1F *h22 = (TH1F*)h1->Clone("h22");
+  h22->Sumw2();
+  h22->Divide(h2);
+  h22->SetLineColor(2);
+  h22->Draw();
+
+  TH1F *h23 = (TH1F*)h1->Clone("h23");
+  h23->Sumw2();
+  h23->Divide(h3);
+  h23->SetLineColor(4);
+  h23->Draw("SAME");
+
+  TH1F *h24 = (TH1F*)h2->Clone("h24");
+  h24->Sumw2();
+  h24->Divide(h3);
+  h24->SetLineColor(8);
+  h24->Draw("SAME");
+
+  TLegend *leg8 = new TLegend(0.7,0.7,0.9,0.9);
+  leg8->AddEntry(h22, "3/4 / EL-REAL", "lep");
+  leg8->AddEntry(h23, "3/4 / EL-CLEAN", "lep");
+  leg8->AddEntry(h24, "EL-REAL / EL-CLEAN", "lep");
+  leg8->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  gPad->SetLogy();
+  T->Draw("P.hgcer.npe[0]>>h4(100,0,50)", "P.hgcer.npe[0]>0.0");
+  T->Draw("P.hgcer.npe[0]>>h5(100,0,50)", "P.hgcer.npe[0]>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npe[0]>>h6(100,0,50)", "P.hgcer.npe[0]>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+
+  h6->SetLineColor(8);
+  h6->SetFillColor(8);
+  h6->SetFillStyle(3006);
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h4, "PMT 1 3/4");
+  leg2->AddEntry(h5, "PMT 1 EL-REAL");
+  leg2->AddEntry(h6, "PMT 1 EL-CLEAN");
+  leg2->Draw();
+
+  can1->cd(5);
+  TH1F *h25 = (TH1F*)h4->Clone("h25");
+  h25->Sumw2();
+  h25->Divide(h5);
+  h25->SetLineColor(2);
+  h25->Draw();
+
+  TH1F *h26 = (TH1F*)h4->Clone("h26");
+  h26->Sumw2();
+  h26->Divide(h6);
+  h26->SetLineColor(4);
+  h26->Draw("SAME");
+
+  TH1F *h27 = (TH1F*)h5->Clone("h27");
+  h27->Sumw2();
+  h27->Divide(h6);
+  h27->SetLineColor(8);
+  h27->Draw("SAME");
+
+  TLegend *leg9 = new TLegend(0.7,0.7,0.9,0.9);
+  leg9->AddEntry(h25, "3/4 / EL-REAL", "lep");
+  leg9->AddEntry(h26, "3/4 / EL-CLEAN", "lep");
+  leg9->AddEntry(h27, "EL-REAL / EL-CLEAN", "lep");
+  leg9->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  gPad->SetLogy();
+  T->Draw("P.hgcer.npe[1]>>h7(100,0,50)", "P.hgcer.npe[1]>0.0");
+  T->Draw("P.hgcer.npe[1]>>h8(100,0,50)", "P.hgcer.npe[1]>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npe[1]>>h9(100,0,50)", "P.hgcer.npe[1]>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+
+  h9->SetLineColor(8);
+  h9->SetFillColor(8);
+  h9->SetFillStyle(3006);
+
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h7, "PMT 2 3/4");
+  leg3->AddEntry(h8, "PMT 2 EL-REAL");
+  leg3->AddEntry(h9, "PMT 2 EL-CLEAN");
+  leg3->Draw();
+
+  can1->cd(6);
+  TH1F *h28 = (TH1F*)h7->Clone("h28");
+  h28->Sumw2();
+  h28->Divide(h8);
+  h28->SetLineColor(2);
+  h28->Draw();
+
+  TH1F *h29 = (TH1F*)h7->Clone("h29");
+  h29->Sumw2();
+  h29->Divide(h9);
+  h29->SetLineColor(4);
+  h29->Draw("SAME");
+
+  TH1F *h30 = (TH1F*)h8->Clone("h30");
+  h30->Sumw2();
+  h30->Divide(h9);
+  h30->SetLineColor(8);
+  h30->Draw("SAME");
+
+  TLegend *leg11 = new TLegend(0.7,0.7,0.9,0.9);
+  leg11->AddEntry(h28, "3/4 / EL-REAL", "lep");
+  leg11->AddEntry(h29, "3/4 / EL-CLEAN", "lep");
+  leg11->AddEntry(h30, "EL-REAL / EL-CLEAN", "lep");
+  leg11->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+  // /-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+
+  TCanvas *can2 = new TCanvas("can2", "can2", 1600, 800);
+  
+  can2->Divide(3, 2);
+  can2->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.eprtracknorm>>h10(75, 0, 1.5)", "P.cal.eprtracknorm>0.0");
+  T->Draw("P.cal.eprtracknorm>>h11(75, 0, 1.5)", "P.cal.eprtracknorm>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.eprtracknorm>>h12(75, 0, 1.5)", "P.cal.eprtracknorm>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h10->SetLineColor(2);
+  h10->SetFillColor(2);
+  h10->SetFillStyle(3004);
+
+  h11->SetLineColor(4);
+  h11->SetFillColor(4);
+  h11->SetFillStyle(3005);
+
+  h12->SetLineColor(8);
+  h12->SetFillColor(8);
+  h12->SetFillStyle(3006);
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h10, "PrSh E-Norm 3/4");
+  leg4->AddEntry(h11, "PrSh E-Norm EL-REAL");
+  leg4->AddEntry(h12, "PrSh E-Norm EL-CLEAN");
+  leg4->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.etracknorm>>h13(75, 0, 1.5)", "P.cal.etracknorm>0.0");
+  T->Draw("P.cal.etracknorm>>h14(75, 0, 1.5)", "P.cal.etracknorm>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.etracknorm>>h15(75, 0, 1.5)", "P.cal.etracknorm>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h13->SetLineColor(2);
+  h13->SetFillColor(2);
+  h13->SetFillStyle(3004);
+
+  h14->SetLineColor(4);
+  h14->SetFillColor(4);
+  h14->SetFillStyle(3005);
+
+  h15->SetLineColor(8);
+  h15->SetFillColor(8);
+  h15->SetFillStyle(3006);
+
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h13, "Sh E-Norm 3/4");
+  leg5->AddEntry(h14, "Sh E-Norm EL-REAL");
+  leg5->AddEntry(h15, "Sh E-Norm EL-CLEAN");
+  leg5->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.etottracknorm>>h31(75, 0, 1.5)", "P.cal.etottracknorm>0.0");
+  T->Draw("P.cal.etottracknorm>>h32(75, 0, 1.5)", "P.cal.etottracknorm>0.0&&T.shms.pEL_REAL_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.etottracknorm>>h33(75, 0, 1.5)", "P.cal.etottracknorm>0.0&&T.shms.pEL_CLEAN_tdcTimeRaw>0.0", "SAME");
+
+  h31->SetLineColor(2);
+  h31->SetFillColor(2);
+  h31->SetFillStyle(3004);
+
+  h32->SetLineColor(4);
+  h32->SetFillColor(4);
+  h32->SetFillStyle(3005);
+
+  h33->SetLineColor(8);
+  h33->SetFillColor(8);
+  h33->SetFillStyle(3006);
+
+  TLegend *leg12 = new TLegend(0.7,0.7,0.9,0.9);
+  leg12->AddEntry(h31, "Total E-Norm 3/4");
+  leg12->AddEntry(h32, "Total E-Norm EL-REAL");
+  leg12->AddEntry(h33, "Total E-Norm EL-CLEAN");
+  leg12->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(4);
+  TH1F *h16 = (TH1F*)h10->Clone("h16");
+  h16->Sumw2();
+  h16->Divide(h11);
+  h16->SetLineColor(2);
+  h16->Draw();
+
+  TH1F *h17 = (TH1F*)h10->Clone("h17");
+  h17->Sumw2();
+  h17->Divide(h12);
+  h17->SetLineColor(4);
+  h17->Draw("SAME");
+
+  TH1F *h18 = (TH1F*)h11->Clone("h18");
+  h18->Sumw2();
+  h18->Divide(h12);
+  h18->SetLineColor(8);
+  h18->Draw("SAME");
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h16, "3/4 / EL-REAL", "lep");
+  leg6->AddEntry(h17, "3/4 / EL-CLEAN", "lep");
+  leg6->AddEntry(h18, "EL-REAL / EL-CLEAN", "lep");
+  leg6->Draw();
+
+  can2->cd(5);
+  TH1F *h19 = (TH1F*)h13->Clone("h19");
+  h19->Sumw2();
+  h19->Divide(h14);
+  h19->SetLineColor(2);
+  h19->Draw();
+
+  TH1F *h20 = (TH1F*)h13->Clone("h20");
+  h20->Sumw2();
+  h20->Divide(h15);
+  h20->SetLineColor(4);
+  h20->Draw("SAME");
+
+  TH1F *h21 = (TH1F*)h14->Clone("h21");
+  h21->Sumw2();
+  h21->Divide(h15);
+  h21->SetLineColor(8);
+  h21->Draw("SAME");
+
+  TLegend *leg7 = new TLegend(0.7,0.7,0.9,0.9);
+  leg7->AddEntry(h19, "3/4 / EL-REAL", "lep");
+  leg7->AddEntry(h20, "3/4 / EL-CLEAN", "lep");
+  leg7->AddEntry(h21, "EL-REAL / EL-CLEAN", "lep");
+  leg7->Draw();
+
+  can2->cd(6);
+  TH1F *h34 = (TH1F*)h31->Clone("h34");
+  h34->Sumw2();
+  h34->Divide(h32);
+  h34->SetLineColor(2);
+  h34->Draw();
+
+  TH1F *h35 = (TH1F*)h31->Clone("h35");
+  h35->Sumw2();
+  h35->Divide(h33);
+  h35->SetLineColor(4);
+  h35->Draw("SAME");
+
+  TH1F *h36 = (TH1F*)h32->Clone("h36");
+  h36->Sumw2();
+  h36->Divide(h34);
+  h36->SetLineColor(8);
+  h36->Draw("SAME");
+
+  TLegend *leg13 = new TLegend(0.7,0.7,0.9,0.9);
+  leg13->AddEntry(h34, "3/4 / EL-REAL", "lep");
+  leg13->AddEntry(h35, "3/4 / EL-CLEAN", "lep");
+  leg13->AddEntry(h36, "EL-REAL / EL-CLEAN", "lep");
+  leg13->Draw();
+  
+}
diff --git a/macros/trig_eff_shms_ll.C b/macros/trig_eff_shms_ll.C
new file mode 100644
index 0000000000000000000000000000000000000000..3af2d6cb365a491978d1f603837008c184ddcc2f
--- /dev/null
+++ b/macros/trig_eff_shms_ll.C
@@ -0,0 +1,101 @@
+{
+  gStyle->SetOptStat(0);
+
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  can1->Divide(3, 2);
+  
+  can1->cd(1);
+  //gPad->SetLogy();
+
+  T->Draw("P.hgcer.npeSum>>h1(100,0,50)", "P.hgcer.npeSum>0.0");
+  T->Draw("P.hgcer.npeSum>>h2(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pHGCSUM_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum");
+  leg1->AddEntry(h2, "CER");
+  leg1->Draw();
+
+  can1->cd(4);
+  TH1F *h3 = (TH1F*)h1->Clone("h3");
+  h3->Sumw2();
+  h3->Divide(h2);
+  h3->SetLineColor(2);
+  h3->Draw();
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h3, "CER", "lep");
+  leg2->Draw();
+
+  //=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  //gPad->SetLogy();
+
+  T->Draw("P.cal.etottracknorm>>h4(150,0,1.5)", "P.cal.etottracknorm>0.0");
+  T->Draw("P.cal.etottracknorm>>h5(150,0,1.5)", "P.cal.etottracknorm>0.0&&T.shms.pPSHWRLO_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+  
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+  
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h4, "ETOT NORM");
+  leg3->AddEntry(h5, "PRSHWRLO");
+  leg3->Draw();
+
+  can1->cd(5);
+  TH1F *h6 = (TH1F*)h4->Clone("h6");
+  h6->Sumw2();
+  h6->Divide(h5);
+  h6->SetLineColor(2);
+  h6->Draw();
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h6, "PRSHWRLO", "lep");
+  leg4->Draw();
+
+ //=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  //gPad->SetLogy();
+
+  T->Draw("P.cal.etottracknorm>>h7(150,0,1.5)", "P.cal.etottracknorm>0.0");
+  T->Draw("P.cal.etottracknorm>>h8(150,0,1.5)", "P.cal.etottracknorm>0.0&&T.shms.pPSHWRHI_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+  
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+  
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h7, "ETOT NORM");
+  leg5->AddEntry(h8, "PRSHWRHI");
+  leg5->Draw();
+
+  can1->cd(6);
+  TH1F *h9 = (TH1F*)h7->Clone("h9");
+  h9->Sumw2();
+  h9->Divide(h8);
+  h9->SetLineColor(2);
+  h9->Draw();
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h9, "PRSHWRHI", "lep");
+  leg6->Draw();
+
+}
diff --git a/macros/trig_eff_shms_ml.C b/macros/trig_eff_shms_ml.C
new file mode 100644
index 0000000000000000000000000000000000000000..49c18865d6e831b09f99db2941838830f1462139
--- /dev/null
+++ b/macros/trig_eff_shms_ml.C
@@ -0,0 +1,325 @@
+{
+  gStyle->SetOptStat(0);
+
+  TCanvas *can1 = new TCanvas("can1", "can1", 1600, 800);
+  
+  can1->Divide(3, 2);
+  can1->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("P.hgcer.npeSum>>h1(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.hgcer.npeSum>>h2(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npeSum>>h3(100,0,50)", "P.hgcer.npeSum>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h1->SetLineColor(2);
+  h1->SetFillColor(2);
+  h1->SetFillStyle(3004);
+  
+  h2->SetLineColor(4);
+  h2->SetFillColor(4);
+  h2->SetFillStyle(3005);
+  
+  h3->SetLineColor(8);
+  h3->SetFillColor(8);
+  h3->SetFillStyle(3006);
+  
+  TLegend *leg1 = new TLegend(0.7,0.7,0.9,0.9);
+  leg1->AddEntry(h1, "NPE Sum EL-LO");
+  leg1->AddEntry(h2, "NPE Sum EL-LO-CER");
+  leg1->AddEntry(h3, "NPE Sum EL-HI");
+  leg1->Draw();
+
+  can1->cd(4);
+  TH1F *h22 = (TH1F*)h1->Clone("h22");
+  h22->Sumw2();
+  h22->Divide(h2);
+  h22->SetLineColor(2);
+  h22->Draw();
+
+  TH1F *h23 = (TH1F*)h1->Clone("h23");
+  h23->Sumw2();
+  h23->Divide(h3);
+  h23->SetLineColor(4);
+  h23->Draw("SAME");
+
+  TH1F *h24 = (TH1F*)h2->Clone("h24");
+  h24->Sumw2();
+  h24->Divide(h3);
+  h24->SetLineColor(8);
+  h24->Draw("SAME");
+
+  TLegend *leg8 = new TLegend(0.7,0.7,0.9,0.9);
+  leg8->AddEntry(h22, "EL-LO / EL-LO-CER", "lep");
+  leg8->AddEntry(h23, "EL-LO / EL-HI", "lep");
+  leg8->AddEntry(h24, "EL-LO-CER / EL-HI", "lep");
+  leg8->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("P.hgcer.npe[0]>>h4(100,0,50)", "P.hgcer.npe[0]>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.hgcer.npe[0]>>h5(100,0,50)", "P.hgcer.npe[0]>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npe[0]>>h6(100,0,50)", "P.hgcer.npe[0]>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h4->SetLineColor(2);
+  h4->SetFillColor(2);
+  h4->SetFillStyle(3004);
+
+  h5->SetLineColor(4);
+  h5->SetFillColor(4);
+  h5->SetFillStyle(3005);
+
+  h6->SetLineColor(8);
+  h6->SetFillColor(8);
+  h6->SetFillStyle(3006);
+
+  TLegend *leg2 = new TLegend(0.7,0.7,0.9,0.9);
+  leg2->AddEntry(h4, "PMT 1 EL-LO");
+  leg2->AddEntry(h5, "PMT 1 EL-LO-CER");
+  leg2->AddEntry(h6, "PMT 1 EL-HI");
+  leg2->Draw();
+
+  can1->cd(5);
+  TH1F *h25 = (TH1F*)h4->Clone("h25");
+  h25->Sumw2();
+  h25->Divide(h5);
+  h25->SetLineColor(2);
+  h25->Draw();
+
+  TH1F *h26 = (TH1F*)h4->Clone("h26");
+  h26->Sumw2();
+  h26->Divide(h6);
+  h26->SetLineColor(4);
+  h26->Draw("SAME");
+
+  TH1F *h27 = (TH1F*)h5->Clone("h27");
+  h27->Sumw2();
+  h27->Divide(h6);
+  h27->SetLineColor(8);
+  h27->Draw("SAME");
+
+  TLegend *leg9 = new TLegend(0.7,0.7,0.9,0.9);
+  leg9->AddEntry(h25, "EL-LO / EL-LO-CER", "lep");
+  leg9->AddEntry(h26, "EL-LO / EL-HI", "lep");
+  leg9->AddEntry(h27, "EL-LO-CER / EL-HI", "lep");
+  leg9->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can1->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("P.hgcer.npe[1]>>h7(100,0,50)", "P.hgcer.npe[1]>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.hgcer.npe[1]>>h8(100,0,50)", "P.hgcer.npe[1]>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.hgcer.npe[1]>>h9(100,0,50)", "P.hgcer.npe[1]>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h7->SetLineColor(2);
+  h7->SetFillColor(2);
+  h7->SetFillStyle(3004);
+
+  h8->SetLineColor(4);
+  h8->SetFillColor(4);
+  h8->SetFillStyle(3005);
+
+  h9->SetLineColor(8);
+  h9->SetFillColor(8);
+  h9->SetFillStyle(3006);
+
+  TLegend *leg3 = new TLegend(0.7,0.7,0.9,0.9);
+  leg3->AddEntry(h7, "PMT 2 EL-LO");
+  leg3->AddEntry(h8, "PMT 2 EL-LO-CER");
+  leg3->AddEntry(h9, "PMT 2 EL-HI");
+  leg3->Draw();
+
+  can1->cd(6);
+  TH1F *h28 = (TH1F*)h7->Clone("h28");
+  h28->Sumw2();
+  h28->Divide(h8);
+  h28->SetLineColor(2);
+  h28->Draw();
+
+  TH1F *h29 = (TH1F*)h7->Clone("h29");
+  h29->Sumw2();
+  h29->Divide(h9);
+  h29->SetLineColor(4);
+  h29->Draw("SAME");
+
+  TH1F *h30 = (TH1F*)h8->Clone("h30");
+  h30->Sumw2();
+  h30->Divide(h9);
+  h30->SetLineColor(8);
+  h30->Draw("SAME");
+
+  TLegend *leg11 = new TLegend(0.7,0.7,0.9,0.9);
+  leg11->AddEntry(h28, "EL-LO / EL-LO-CER", "lep");
+  leg11->AddEntry(h29, "EL-LO / EL-HI", "lep");
+  leg11->AddEntry(h30, "EL-LO-CER / EL-HI", "lep");
+  leg11->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+  // /-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+
+  TCanvas *can2 = new TCanvas("can2", "can2", 1600, 800);
+  
+  can2->Divide(3, 2);
+  can2->cd(1);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.eprtracknorm>>h10(75, 0, 1.5)", "P.cal.eprtracknorm>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.cal.eprtracknorm>>h11(75, 0, 1.5)", "P.cal.eprtracknorm>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.eprtracknorm>>h12(75, 0, 1.5)", "P.cal.eprtracknorm>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h10->SetLineColor(2);
+  h10->SetFillColor(2);
+  h10->SetFillStyle(3004);
+
+  h11->SetLineColor(4);
+  h11->SetFillColor(4);
+  h11->SetFillStyle(3005);
+
+  h12->SetLineColor(8);
+  h12->SetFillColor(8);
+  h12->SetFillStyle(3006);
+
+  TLegend *leg4 = new TLegend(0.7,0.7,0.9,0.9);
+  leg4->AddEntry(h10, "PrSh E-Norm EL-LO");
+  leg4->AddEntry(h11, "PrSh E-Norm EL-LO-CER");
+  leg4->AddEntry(h12, "PrSh E-Norm EL-HI");
+  leg4->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(2);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.etracknorm>>h13(75, 0, 1.5)", "P.cal.etracknorm>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.cal.etracknorm>>h14(75, 0, 1.5)", "P.cal.etracknorm>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.etracknorm>>h15(75, 0, 1.5)", "P.cal.etracknorm>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h13->SetLineColor(2);
+  h13->SetFillColor(2);
+  h13->SetFillStyle(3004);
+
+  h14->SetLineColor(4);
+  h14->SetFillColor(4);
+  h14->SetFillStyle(3005);
+
+  h15->SetLineColor(8);
+  h15->SetFillColor(8);
+  h15->SetFillStyle(3006);
+
+  TLegend *leg5 = new TLegend(0.7,0.7,0.9,0.9);
+  leg5->AddEntry(h13, "Sh E-Norm EL-LO");
+  leg5->AddEntry(h14, "Sh E-Norm EL-LO-CER");
+  leg5->AddEntry(h15, "Sh E-Norm EL-HI");
+  leg5->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(3);
+  gPad->SetLogy();
+
+  T->Draw("P.cal.etottracknorm>>h31(75, 0, 1.5)", "P.cal.etottracknorm>0.0&&T.shms.pEL_LO_tdcTimeRaw>0.0");
+  T->Draw("P.cal.etottracknorm>>h32(75, 0, 1.5)", "P.cal.etottracknorm>0.0&&T.shms.pEL_LO_CER_tdcTimeRaw>0.0", "SAME");
+  T->Draw("P.cal.etottracknorm>>h33(75, 0, 1.5)", "P.cal.etottracknorm>0.0&&T.shms.pEL_HI_tdcTimeRaw>0.0", "SAME");
+
+  h31->SetLineColor(2);
+  h31->SetFillColor(2);
+  h31->SetFillStyle(3004);
+
+  h32->SetLineColor(4);
+  h32->SetFillColor(4);
+  h32->SetFillStyle(3005);
+
+  h33->SetLineColor(8);
+  h33->SetFillColor(8);
+  h33->SetFillStyle(3006);
+
+  TLegend *leg12 = new TLegend(0.7,0.7,0.9,0.9);
+  leg12->AddEntry(h31, "Total E-Norm EL-LO");
+  leg12->AddEntry(h32, "Total E-Norm EL-LO-CER");
+  leg12->AddEntry(h33, "Total E-Norm EL-HI");
+  leg12->Draw();
+
+  // =:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:=:
+
+  can2->cd(4);
+  TH1F *h16 = (TH1F*)h10->Clone("h16");
+  h16->Sumw2();
+  h16->Divide(h11);
+  h16->SetLineColor(2);
+  h16->Draw();
+
+  TH1F *h17 = (TH1F*)h10->Clone("h17");
+  h17->Sumw2();
+  h17->Divide(h12);
+  h17->SetLineColor(4);
+  h17->Draw("SAME");
+
+  TH1F *h18 = (TH1F*)h11->Clone("h18");
+  h18->Sumw2();
+  h18->Divide(h12);
+  h18->SetLineColor(8);
+  h18->Draw("SAME");
+
+  TLegend *leg6 = new TLegend(0.7,0.7,0.9,0.9);
+  leg6->AddEntry(h16, "EL-LO / EL-LO-CER", "lep");
+  leg6->AddEntry(h17, "EL-LO / EL-HI", "lep");
+  leg6->AddEntry(h18, "EL-LO-CER / EL-HI", "lep");
+  leg6->Draw();
+
+  can2->cd(5);
+  TH1F *h19 = (TH1F*)h13->Clone("h19");
+  h19->Sumw2();
+  h19->Divide(h14);
+  h19->SetLineColor(2);
+  h19->Draw();
+
+  TH1F *h20 = (TH1F*)h13->Clone("h20");
+  h20->Sumw2();
+  h20->Divide(h15);
+  h20->SetLineColor(4);
+  h20->Draw("SAME");
+
+  TH1F *h21 = (TH1F*)h14->Clone("h21");
+  h21->Sumw2();
+  h21->Divide(h15);
+  h21->SetLineColor(8);
+  h21->Draw("SAME");
+
+  TLegend *leg7 = new TLegend(0.7,0.7,0.9,0.9);
+  leg7->AddEntry(h19, "EL-LO / EL-LO-CER", "lep");
+  leg7->AddEntry(h20, "EL-LO / EL-HI", "lep");
+  leg7->AddEntry(h21, "EL-LO-CER / EL-HI", "lep");
+  leg7->Draw();
+
+  can2->cd(6);
+  TH1F *h34 = (TH1F*)h31->Clone("h34");
+  h34->Sumw2();
+  h34->Divide(h32);
+  h34->SetLineColor(2);
+  h34->Draw();
+
+  TH1F *h35 = (TH1F*)h31->Clone("h35");
+  h35->Sumw2();
+  h35->Divide(h33);
+  h35->SetLineColor(4);
+  h35->Draw("SAME");
+
+  TH1F *h36 = (TH1F*)h32->Clone("h36");
+  h36->Sumw2();
+  h36->Divide(h34);
+  h36->SetLineColor(8);
+  h36->Draw("SAME");
+
+  TLegend *leg13 = new TLegend(0.7,0.7,0.9,0.9);
+  leg13->AddEntry(h34, "EL-LO / EL-LO-CER", "lep");
+  leg13->AddEntry(h35, "EL-LO / EL-HI", "lep");
+  leg13->AddEntry(h36, "EL-LO-CER / EL-HI", "lep");
+  leg13->Draw();
+
+}