OnMonPlotMaker.C File Reference

#include <fstream>
#include <TH1F.h>
#include <TH2F.h>
#include <TGraph.h>
#include <TTree.h>
#include <TCanvas.h>
#include <TFile.h>
#include <TPaveText.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TDatime.h>
#include <TLegend.h>
#include <TLegendEntry.h>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <string>
#include <cstring>
#include "OnlineMonitoring/viewer/HwDetLabel.cxx"

Go to the source code of this file.

Functions
void	OnMonPlotMaker ()
void	OnMonPlotMaker (TString det_type, TString t_prd, TString trigger)

Function Documentation

void OnMonPlotMaker

(

)

Definition at line 24 of file OnMonPlotMaker.C.

{}

void OnMonPlotMaker	(	TString	det_type,
		TString	t_prd,
		TString	trigger
	)

Definition at line 26 of file OnMonPlotMaker.C.

References a, std::abs(), readPickles::body, run_hadd::cmd, om::HwDetLabel::Config(), GenerateCanaryReport::content, om::cout, febshutoff_auto::day, fillBadChanDBTables::det, allInOneTrainingPlots::df, om::HwDetLabel::Draw(), E, e, allTimeWatchdog::endl, shutoffs::filename, GMToffset, submit_syst::header, febshutoff_auto::hour, MECModelEnuComparisons::i, inFile, makeTrainCVSamples::int, calib::j, kBlue, kGreen, kOrange, kRed, central_limit::l1, LastRun, LastSR, MECModelEnuComparisons::leg, livetime, std::log10(), cet::sqlite::max(), min(), plot_xsec_2d::NbinsY, norm, Npar, pad1, pad2, pad3, par, DCS_db_parser::period, PLOT_DIR, PandAna.Demos.demo1::range, updateRunHistoryTables::run, SRtime, string, getGoodRuns4SAM::subrun, system(), dump_to_simple_cpp::threshold, submit_nova_art::timestamp, plotROC::title, TNow, TNowGMT, produceInfoGainMatrix::total, UpdateText, WarningText, and Munits::year.

                                                                      {
  
  const std::string DATA_DIR = "/nusoft/app/web/htdoc/nova/datacheck/nearline/data";
  const std::string PLOT_DIR = "/nusoft/app/web/htdoc/nova/datacheck/nearline/plots";
  //const std::string PLOT_DIR = "/nusoft/app/web/htdoc/nova/datacheck/nearline/test_plots";

  // TODO:
  //   *  make the histograms smart enough to get their binnings from the onmon-histos file
  //      (or make them smart enough to figure out how big they should be...???)  In fact,
  //      make nothing hardcoded about the histos (get all info from onmon-histos.csv)
  //
  //   *  make hitmap style plot of DCM duty cycle (using the empty DCM histo)
  //
  //   *  fix this so that if one of the text files is empty, it doesn't crash
  //
  //   *  document the heck out of this code...
  //
  //   *  should I normalize the NFEBperSR histos so that each bin is the average of each subrun that contributed?
  //
  //   *  I think that I am doing the GMT to CST conversion correctly, but I also think that I'm doing it in a convoluted way.  Maybe I should make this better...
  //
  //   *  I don't think that I need to count all of the things for the TGraphs every time.  Maybe I only need to count them during the month loop and then
  //      just zoom each plot appropriately?
  //
  //   *  normalize Nhits per spill plot to be percentage???
  //
  //   *  figure out a better way to do the arrays (other than to have fixed sizes...)
  //
  //   *  make the rasta plot have 1 hour bins (since subruns are now ~20 minutes)???
  //
  //   * keep track of last run/subrun (for the text boxes) by partition




  TDatime *TNowGMT = new TDatime;  // current GMT time
  TDatime *TNow    = new TDatime;  // current local time
  TDatime *SRtime  = new TDatime;
  TDatime *SRtime_fine = new TDatime;

  int XNowGMT      = TNowGMT->Convert(kTRUE);
  int XNow         = TNow   ->Convert(kFALSE);

  int GMToffset    = XNowGMT - XNow;
  XNow             = XNow - GMToffset;

  int XDayAgo       = XNow - 60*60*24;
  int XWeekAgo      = XNow - 60*60*24*7;
  int XMonthAgo     = XNow - 60*60*24*30;
  int XSixMonthsAgo = XNow - 60*60*24*26*7;
  int XYearAgo      = XNow - 60*60*24*52*7;
  
  int Xsrtime   = 0;
  int Xsrtime_fine = 0;
  int time_ago  = 0;

  string det;
  string period;
  string trig;
  string trigname;

  if(det_type == "FarDet"){
    det      = "FarDet";
    trig     = "t02";
    trigname = " ";

    if(t_prd == "Day"){
      period   = "Day";
      time_ago = XDayAgo;
    }
    else if(t_prd == "Week"){
      period   = "Week";
      time_ago = XWeekAgo;
    }
    else if(t_prd == "Month"){
      period   = "Month";
      time_ago = XMonthAgo;
    }
    else if(t_prd == "6Months"){
      period = "6Months";
      time_ago = XSixMonthsAgo;
    }
    else if(t_prd == "Year"){
      period = "Year";
      time_ago = XYearAgo;
    }
  }

  if(det_type == "NearDet"){
    det = "NearDet";

    if (trigger == "t00"){
      trig     = "t00";
      trigname = "NuMI";
    }
    else if(trigger == "DDActivity1"){
      trig     = "DDActivity1";
      trigname = "DDActivity";
    }
    else if(trigger == "t02"){
      trig = "t02";
      trigname = "Cosmic";
    }
    if(t_prd == "Day"){
      period   = "Day";
      time_ago = XDayAgo;
    }
    else if(t_prd == "Week"){
      period   = "Week";
      time_ago = XWeekAgo;
    }
    else if(t_prd == "Month"){
      period   = "Month";
      time_ago = XMonthAgo;
    }
    else if(t_prd == "6Months"){
      period = "6Months";
      time_ago = XSixMonthsAgo;
    }
    else if(t_prd == "Year"){
      period = "Year";
      time_ago = XYearAgo;
    }
  }

  if(det_type == "TestBeam"){
    det      = "TestBeam";

    if (trigger == "Beamline"){
      trig     = "Beamline";
      trigname = "Beamline";
    }
    else if(trigger == "DDActivity1"){
      trig     = "DDActivity1";
      trigname = "DDActivity";
    }
    else if(trigger == "t02"){
      trig = "t02";
      trigname = "Cosmic";
    }
    else if(trigger == "Spill"){
      trig = "Spill";
      trigname = "Full Spill";
    }
    if(t_prd == "Day"){
      period   = "Day";
      time_ago = XDayAgo;
    }
    else if(t_prd == "Week"){
      period   = "Week";
      time_ago = XWeekAgo;
    }
    else if(t_prd == "Month"){
      period   = "Month";
      time_ago = XMonthAgo;
    }
    else if(t_prd == "6Months"){
      period = "6Months";
      time_ago = XSixMonthsAgo;
    }
    else if(t_prd == "Year"){
      period = "Year";
      time_ago = XYearAgo;
    }
  }

  unsigned int       year  = 0;
  unsigned short int month = 0;
  unsigned short int day   = 0;
  double Hour = 0.0;
  double Min  = 0.0;
  int hour = 0;
  int min  = 0;
  int sec  = 0;

  unsigned int run    = 0;
  unsigned int subrun = 0;
  UInt_t       par    = 0;

  unsigned int LastRun = 0;
  unsigned int LastSR  = 0;

  unsigned int       yearEnd  = 0;
  unsigned short int monthEnd = 0;
  unsigned short int dayEnd   = 0;
  double HourEnd = 0.0;
  double SRdur   = 0;

  // Root is stupid and you have to make it think that it is in a different time zone so that it will draw any
  // plots with time on the X-axis correctly.
  gSystem->Setenv("TZ","UTC");
  gStyle->SetTimeOffset(0);



  // Define the number of partitions in the data.  Make this number one more
  // than the actual number so that there can be a default unassigned value.
  // Define the maximum number of points to go into a TGraph (assuming a max
  // of one subrun every minute.)
  const unsigned int Npar   = 5;
  double livetime = 0.0;

  // Define the number of errors in the AE plot and alerts in the AA plot.
  const unsigned int Nerr = 10;
  const unsigned int Nalert = 5;

  int t_n = 1;
  if(t_prd == "Day") t_n = 1;
  if(t_prd == "Week") t_n = 7;
  if(t_prd == "Month") t_n = 30;
  if(t_prd == "6Months") t_n = 7*26;
  if(t_prd == "Year") t_n = 7*52;
  const unsigned int Npoint = 1440 * t_n;



  //
  // Make plots
  //
  // make titles for lists [in testing for automatic trigger naming of list files]
  // char filelist_title[128];
  
  // Spill files only exist for ND/FD
  char SpillInfoFile[128]="";
  //char BNBSpillInfoFile[128]="";
  if (det_type == "NearDet" || det_type == "FarDet"){
    
    if(t_prd == "Day")
      sprintf(SpillInfoFile, "%s/POTSpillRateDay.root",DATA_DIR.c_str());
    if(t_prd == "Week")
      sprintf(SpillInfoFile, "%s/POTSpillRateWeek.root",DATA_DIR.c_str());
    if(t_prd == "Month")
      sprintf(SpillInfoFile, "%s/POTSpillRateMonth.root",DATA_DIR.c_str());
    if(t_prd == "6Months")
      sprintf(SpillInfoFile, "%s/POTSpillRate6Months.root",DATA_DIR.c_str());
    if(t_prd == "Year")
      sprintf(SpillInfoFile, "%s/POTSpillRateYear.root",DATA_DIR.c_str());
    
    
    // if(t_prd == "Day")
    //sprintf(BNBSpillInfoFile, "%s/BNBPOTSpillRateDay.root",DATA_DIR.c_str());
    //if(t_prd == "Week")
    //sprintf(BNBSpillInfoFile, "%s/BNBPOTSpillRateWeek.root",DATA_DIR.c_str());
    //if(t_prd == "Month")
    //sprintf(BNBSpillInfoFile, "%s/BNBPOTSpillRateMonth.root",DATA_DIR.c_str());
    //if(t_prd == "6Months")
    //sprintf(BNBSpillInfoFile, "%s/BNBPOTSpillRate6Months.root",DATA_DIR.c_str());
    //if(t_prd == "Year")
    //sprintf(BNBSpillInfoFile, "%s/BNBPOTSpillRateYear.root",DATA_DIR.c_str());
  }
  
  std::ifstream inFile;
  if(det_type == "FarDet"){
    if(t_prd == "Day"){
      //      sprintf(filelist_title,"FDOnMonWeekFileList_t02-t00.txt");
      inFile.open("FDOnMonDayFileList_t02-t00.txt");
    }
    else if(t_prd == "Week"){
      //sprintf(filelist_title,"FDOnMonWeekFileList_%s.txt",trig.c_str());
      inFile.open("FDOnMonWeekFileList_t02-t00.txt");
    }
    else if(t_prd == "Month"){
      //sprintf(filelist_title,"FDOnMonMonthFileList_%s.txt",trig.c_str());
      inFile.open("FDOnMonMonthFileList_t02-t00.txt");
    }
    else if(t_prd == "6Months"){
     inFile.open("FDOnMon6MonthsFileList_t02-t00.txt");
    }
    else if(t_prd == "Year"){
     inFile.open("FDOnMonYearFileList_t02-t00.txt");
    }
  }
  if(det_type == "NearDet"){

    if(t_prd == "Day"){
     // sprintf(filelist_title,"NDOnMonDayFileList_%s.txt",trig.c_str());
     if(trig == "t00")
       inFile.open("NDOnMonDayFileList_t00.txt");
     if(trig == "t02")
       inFile.open("NDOnMonDayFileList_t02.txt");
     else if(trig == "DDActivity1")
       inFile.open("NDOnMonDayFileList_DDActivity1.txt");
    }
    else if(t_prd == "Week"){
     // sprintf(filelist_title,"NDOnMonWeekFileList_%s.txt",trig.c_str());
      if(trig == "t00")
       inFile.open("NDOnMonWeekFileList_t00.txt");
      if(trig == "t02")
       inFile.open("NDOnMonWeekFileList_t02.txt");
      else if(trig == "DDActivity1")
       inFile.open("NDOnMonWeekFileList_DDActivity1.txt");
    }
    else if(t_prd == "Month"){
     // sprintf(filelist_title,"NDOnMonMonthFileList_%s.txt",trig.c_str());
      if(trig == "t00")
       inFile.open("NDOnMonMonthFileList_t00.txt");
      if(trig == "t02")
       inFile.open("NDOnMonMonthFileList_t02.txt");
     else if(trig == "DDActivity1")
       inFile.open("NDOnMonMonthFileList_DDActivity1.txt");
    }
    else if(t_prd == "6Months"){
      if(trig== "t00")
        inFile.open("NDOnMon6MonthsFileList_t00.txt");
      if(trig == "t02")
        inFile.open("NDOnMon6MonthsFileList_t02.txt");
      else if(trig == "DDActivity1")
        inFile.open("NDOnMon6MonthsFileList_DDActivity1.txt");
    }
    else if(t_prd == "Year"){
      if(trig== "t00")
        inFile.open("NDOnMonYearFileList_t00.txt");
      if(trig== "t02")
        inFile.open("NDOnMonYearFileList_t02.txt");
      else if(trig == "DDActivity1")
        inFile.open("NDOnMonYearFileList_DDActivity1.txt");
    }
  }

  if(det_type == "TestBeam"){

    if(t_prd == "Day"){
     if(trig == "Beamline")
       inFile.open("TBOnMonDayFileList_Beamline.txt");
     else if(trig == "t02")
       inFile.open("TBOnMonDayFileList_t02.txt");
     else if(trig == "DDActivity1")
       inFile.open("TBOnMonDayFileList_DDActivity1.txt");
     else if(trig == "Spill")
       inFile.open("TBOnMonDayFileList_Spill.txt");
    }
    else if(t_prd == "Week"){
      if(trig == "Beamline")
       inFile.open("TBOnMonWeekFileList_Beamline.txt");
      else if(trig == "t02")
       inFile.open("TBOnMonWeekFileList_t02.txt");
      else if(trig == "DDActivity1")
       inFile.open("TBOnMonWeekFileList_DDActivity1.txt");
      else if(trig == "Spill")
        inFile.open("TBOnMonWeekFileList_Spill.txt");
    }
    else if(t_prd == "Month"){
      if(trig == "Beamline")
       inFile.open("TBOnMonMonthFileList_Beamline.txt");
      else if(trig == "t02")
       inFile.open("TBOnMonMonthFileList_t02.txt");
     else if(trig == "DDActivity1")
       inFile.open("TBOnMonMonthFileList_DDActivity1.txt");
     else if(trig == "Spill")
       inFile.open("TBOnMonMonthFileList_Spill.txt");
    }
    else if(t_prd == "6Months"){
      if(trig== "Beamline")
        inFile.open("TBOnMon6MonthsFileList_Beamline.txt");
      if(trig == "t02")
        inFile.open("TBOnMon6MonthsFileList_t02.txt");
      else if(trig == "DDActivity1")
        inFile.open("TBOnMon6MonthsFileList_DDActivity1.txt");
      else if(trig == "Spill")
        inFile.open("TBOnMon6MonthsFileList_Spill.txt");
    }
    else if(t_prd == "Year"){
      if(trig== "Beamline")
        inFile.open("TBOnMonYearFileList_Beamline.txt");
      else if(trig== "t02")
        inFile.open("TBOnMonYearFileList_t02.txt");
      else if(trig == "DDActivity1")
        inFile.open("TBOnMonYearFileList_DDActivity1.txt");
      else if(trig == "Spill")
        inFile.open("TBOnMonYearFileList_Spill.txt");
    }
  }

  // fix trig string
  if(trigger == "DDActivity1")
    trig = "Act";

  // Make array of histo names and titles for the all errors plots.
  // Histo titles taken from OnlineMonitoring/viewer/PlotOptions.cxx. If they ever change there (which they never will) they need to change here as well.
  // 19 July 2018: They changed. We were mislabelling 5 errors/alerts. This has now been fixed in OnMon.
  std::string errorHistoNames[Nerr];

  errorHistoNames[0] = "ErrorEventIncomplete";
  errorHistoNames[1] = "ErrorDataMissing";
  errorHistoNames[2] = "ErrorMilliSliceIncomplete"; // This is never actually filled in OnMon. Plot will always be blank.
  errorHistoNames[3] = "ErrorDCMDataPresent";
  errorHistoNames[4] = "ErrorComm";
  errorHistoNames[5] = "ErrorADC";
  errorHistoNames[6] = "ErrorTECC";
  errorHistoNames[7] = "ErrorDataDrop";
  errorHistoNames[8] = "ErrorLinkStatus";
  errorHistoNames[9] = "ErrorFEBDataPresent";

  std::string errorHistoTitles[Nerr];

  errorHistoTitles[0] = "Event Incomplete Error";
  errorHistoTitles[1] = "Data Missing Error";
  errorHistoTitles[2] = "Millislice Incomplete Error";
  errorHistoTitles[3] = "DCM Data Present Error";
  errorHistoTitles[4] = "Comm. Error";
  errorHistoTitles[5] = "ADC Error";
  errorHistoTitles[6] = "TECC Error";
  errorHistoTitles[7] = "Data Drop Error";
  errorHistoTitles[8] = "Link Status Error";
  errorHistoTitles[9] = "FEB Data Present Error";

  // Make array of histo names and titles for all alert plots.
  std::string alertHistoNames[Nalert];

  alertHistoNames[0] = "AlertSimulationsFlag";
  alertHistoNames[1] = "AlertDAQEnabled";
  alertHistoNames[2] = "AlertTECCEnable";
  alertHistoNames[3] = "AlertFIFOFull";
  alertHistoNames[4] = "AlertDebugMode";

  std::string alertHistoTitles[Nalert];

  alertHistoTitles[0] = "Simulations Flag";
  alertHistoTitles[1] = "DAQ Enabled";
  alertHistoTitles[2] = "TECC Enabled";
  alertHistoTitles[3] = "Event FIFO Full";
  alertHistoTitles[4] = "Debug Mode";


  // book histos, arrays for TGraphs, etc.
  TH1F *hErrorHistos[Npar][Nerr];
  TH1F *hAlertHistos[Npar][Nalert];

  TH1F *hQDay[Npar];
  TH2F *hFEBRateHitMapDay[Npar];
  TH2F *hFEBDropOutMapDay[Npar];
  TH2F *hFEBTimeHitMapDay[Npar];
  TH2F *hFEBaveADCHitMapDay[Npar];
  TH2F *hFEBhitsHitMapDay[Npar];

  TH1F *hFEBRateSpectrumDay[Npar];
  TH2F *hFEBRateSpectrumVsTimeDay[Npar];

  TH1F *hTPlotALLDay[Npar];
  TH2F *hTPlotALLVsTimeDay[Npar];

  TH1F *hTPlotZOOMDay[Npar];
  TH2F *hTPlotZOOMVsTimeDay[Npar];

  TH1F *hNormalizationVsTimeDay[Npar];
  TH1F *hTimingNormVsTimeDay[Npar];

  int *NDCMtimeDay[Npar];
  int *NDCMDay[Npar];
  int NDCMDayCount[Npar];
  int *NFEBtimeDay[Npar];
  int *NFEBDay[Npar];
  int NFEBDayCount[Npar];
  int *NFEBtimeDDay[Npar];
  int *NFEBDDay[Npar];
  int NFEBDDayCount[Npar];
  int *NPIXperSRtimeDay[Npar];
  int *NPIXperSRDay[Npar];
  int NPIXperSRDayCount[Npar];
  int *NPIXperEVTtimeDay[Npar];
  int *NPIXperEVTDay[Npar];
  int NPIXperEVTDayCount[Npar];

  double *Trig00timeDay[Npar];
  double *Trig00Day[Npar];
  double *Trig01timeDay[Npar];
  double *Trig01Day[Npar];
  double *Trig02timeDay[Npar];
  double *Trig02Day[Npar];
  double *TrigActtimeDay[Npar];
  double *TrigActDay[Npar];
  int    NTrigDayCount[Npar];

  int    NSpillEntryCount[Npar];
  double    *SpillCountValue[Npar];
  double    *SpillCountTime[Npar];

  //  int    BNBNSpillEntryCount[Npar];
  //double    *BNBSpillCountValue[Npar];
  //double    *BNBSpillCountTime[Npar];

  double *dcm0timeDay[Npar];
  double *Ndcm0Day[Npar];
  int    Ndcm0DayCount[Npar];



  // assigning for every partition an array of dimension Npoint.
  for(unsigned int i = 0; i < Npar; i++){
      NDCMtimeDay[i] = new int[Npoint];
      NDCMDay[i] = new int[Npoint];
      NFEBtimeDay[i] = new int[Npoint];
      NFEBtimeDDay[i] = new int[Npoint];
      NFEBDay[i] = new int[Npoint];
      NFEBDDay[i] = new int[Npoint];
      NPIXperSRtimeDay[i] = new int[Npoint];
      NPIXperSRDay[i] = new int[Npoint];
      NPIXperEVTtimeDay[i] = new int[Npoint];
      NPIXperEVTDay[i] = new int[Npoint];

      Trig00timeDay[i] = new double[Npoint];
      Trig00Day[i] = new double[Npoint];
      Trig01timeDay[i] = new double[Npoint];
      Trig01Day[i] = new double[Npoint];
      Trig02timeDay[i] = new double[Npoint];
      Trig02Day[i] = new double[Npoint];
      TrigActtimeDay[i] = new double[Npoint];
      TrigActDay[i] = new double[Npoint];

      SpillCountValue[i] = new double[Npoint];
      SpillCountTime[i] = new double[Npoint];

      //BNBSpillCountValue[i] = new double[Npoint];
      //BNBSpillCountTime[i] = new double[Npoint];

      dcm0timeDay[i] = new double[Npoint];
      Ndcm0Day[i] = new double[Npoint];
  }

  string time;
  if(period == "Day")   time = "past 24 hrs.";
  if(period == "Week")  time = "past 7 days";
  if(period == "Month") time = "past 30 days";
  if(period == "6Months") time = "past 26 week";
  if(period == "Year") time = "past 52 weeks";

  int FEBxbin_n   = 0;
  int FEBxbin_min = 0;
  int FEBxbin_max = 0;

  int FEBybin_n   = 0;
  int FEBybin_min = 0;
  int FEBybin_max = 0;

  int Timebin_n   = 72;
  // Since NearDet (and TestBeam) subruns are longer than 20 minutes, we will use fewer than 72 bins
  // for the 24 hour NearDet (and TestBeam) plots to avoid creating vertical white stripes.
  if( (det_type == "NearDet" || det_type == "TestBeam") && period == "Day") Timebin_n = 18;

  if(det_type == "FarDet"){
    FEBxbin_n   = 224;
    FEBxbin_min = 0;
    FEBxbin_max = 896;
    FEBybin_n   = 48;
    FEBybin_min = 0;
    FEBybin_max = 384;
  }
  if(det_type == "NearDet"){
    FEBxbin_n   = 64;
    FEBxbin_min = 0;
    FEBxbin_max = 256;
    FEBybin_n   = 16;
    FEBybin_min = 0;
    FEBybin_max = 128;
  }
  if(det_type == "TestBeam"){
    FEBxbin_n   = 16;
    FEBxbin_min = 0;
    FEBxbin_max = 64;
    FEBybin_n   = 12;
    FEBybin_min = 0;
    FEBybin_max = 96;
  }

  // initalize histos
  for(unsigned int p = 0; p < Npar; ++p) {
    char histoname[56];
    char histotitle[128];

    for(unsigned int e = 0; e < Nerr; ++e) {
      sprintf(histoname,"h%sP%.1u",errorHistoNames[e].c_str(),p);
      sprintf(histotitle,"%s %s (%s) - partition %.1u;(central time);",trigname.c_str(),errorHistoTitles[e].c_str(),time.c_str(),p);
      hErrorHistos[p][e] = new TH1F(histoname,histotitle,1440,time_ago,XNow);
      // NOTE:   Using 1440 bins since that will line up with the original AE histo X-axis bin size. This allows for easier comparisons to the original histogram.
    }

    for(unsigned int a = 0; a < Nalert; ++a) {
      sprintf(histoname,"h%sP%.1u",alertHistoNames[a].c_str(),p);
      sprintf(histotitle,"%s %s (%s) - partition %.1u;(central time);",trigname.c_str(),alertHistoTitles[a].c_str(),time.c_str(),p);
      hAlertHistos[p][a] = new TH1F(histoname,histotitle,1440,time_ago,XNow);
      // NOTE:   Using 1440 bins since that will line up with the original AA histo X-axis bin size. This allows for easier comparisons to the original histogram.
    }

    sprintf(histoname,"hQDayP%.1u",p);
    sprintf(histotitle,"%s ADC Spectra (%s) - partition %.1u;ADC;",trigname.c_str(),time.c_str(),p);
    hQDay[p] = new TH1F(histoname,histotitle,4100,0,4100);

    sprintf(histoname,"hFEBRateHitMapDayP%.1u",p);
    sprintf(histotitle,"%s FEB Hit Rates (%s) - partition %.1u",trigname.c_str(),time.c_str(),p);
    hFEBRateHitMapDay[p] = new TH2F(histoname,histotitle,FEBxbin_n,FEBxbin_min,FEBxbin_max,FEBybin_n,FEBybin_min,FEBybin_max);

    sprintf(histoname, "hFEBDropOutMapDayP%.1u", p);
    sprintf(histotitle,"%s FEB Number of Drop Outs (%s) - partition %.1u",trigname.c_str(),time.c_str(),p);
    hFEBDropOutMapDay[p] = new TH2F(histoname,histotitle,FEBxbin_n,FEBxbin_min,FEBxbin_max,FEBybin_n,FEBybin_min,FEBybin_max);

    sprintf(histoname,"hFEBTimeHitMapDayP%.1u",p);
    sprintf(histotitle,"%s FEB Live Time (%s) - partition %.1u",trigname.c_str(),time.c_str(),p);
    hFEBTimeHitMapDay[p] = new TH2F(histoname,histotitle,FEBxbin_n,FEBxbin_min,FEBxbin_max,FEBybin_n,FEBybin_min,FEBybin_max);

    sprintf(histoname,"hFEBaveADCHitMapDayP%.1u",p);
    sprintf(histotitle,"%s Average ADC per FEB (%s) - partition %.1u",trigname.c_str(),time.c_str(),p);
    hFEBaveADCHitMapDay[p] = new TH2F(histoname,histotitle,FEBxbin_n,FEBxbin_min,FEBxbin_max,FEBybin_n,FEBybin_min,FEBybin_max);

    sprintf(histoname,"hFEBhitsHitMapDayP%.1u",p);
    sprintf(histotitle,"%s FEB Number of Hits (%s) - partition %.1u",trigname.c_str(),time.c_str(),p);
    hFEBhitsHitMapDay[p] = new TH2F(histoname,histotitle,FEBxbin_n,FEBxbin_min,FEBxbin_max,FEBybin_n,FEBybin_min,FEBybin_max);

    sprintf(histoname,"hFEBRateSpectrumDayP%.1u",p);
    sprintf(histotitle,"%s FEB Hit Rate Spectrum (%s) - partition %.1u;Log_{10}(rate);",trigname.c_str(),time.c_str(),p);
    hFEBRateSpectrumDay[p] = new TH1F(histoname,histotitle,200,1,6);

    sprintf(histoname,"hFEBRateSpectrumVsTimeDayP%.1u",p);
    sprintf(histotitle,"%s FEB Hit Rate Spectra vs. Time - partition %.1u;(central time);Log_{10}(rate)",trigname.c_str(),p);
    hFEBRateSpectrumVsTimeDay[p] = new TH2F(histoname,histotitle,Timebin_n,time_ago,XNow,200,1,6);

    sprintf(histoname,"hTPlotALLDayP%.1u",p);
    sprintf(histotitle,"%s Time Distribution for All Hits (%s) - partition %.1u;Hit Time [#musec];",trigname.c_str(),time.c_str(),p);
    hTPlotALLDay[p] = new TH1F(histoname,histotitle,350,-100,600);
    sprintf(histoname,"hTPlotALLVsTimeDayP%.1u",p);
    sprintf(histotitle,"%s Time Distribution for All Hits vs. Time - partition %.1u;(central time);Hit Time [#musec]",trigname.c_str(),p);
    hTPlotALLVsTimeDay[p] = new TH2F(histoname,histotitle,Timebin_n,time_ago,XNow,350,-100,600);
    sprintf(histoname,"hTPlotZOOMDayP%.1u",p);
    sprintf(histotitle,"%s Time Distribution for All Hits (%s) - partition %.1u;Hit Time [#musec];",trigname.c_str(),time.c_str(),p);
    hTPlotZOOMDay[p] = new TH1F(histoname,histotitle,200,212,237);
    sprintf(histoname,"hTPlotZOOMVsTimeDayP%.1u",p);
    sprintf(histotitle,"%s Time Distribution for All Hits vs. Time - partition %.1u;(central time);Hit Time [#musec]",trigname.c_str(),p);
    hTPlotZOOMVsTimeDay[p] = new TH2F(histoname,histotitle,Timebin_n,time_ago,XNow,200,212,237);
    sprintf(histoname,"hNormalizationVsTimeDayP%.1u",p);
    sprintf(histotitle,"%s Day Normalization - partition %.1u",trigname.c_str(),p);
    hNormalizationVsTimeDay[p] = new TH1F(histoname,histotitle,Timebin_n,time_ago,XNow);
    sprintf(histoname,"hTimingNormVsTimeDayP%.1u",p);
    sprintf(histotitle,"%s Day Timing Normalization - partition %.1u",trigname.c_str(),p);
    hTimingNormVsTimeDay[p] = new TH1F(histoname,histotitle,Timebin_n,time_ago,XNow);
    NDCMDayCount[p]       = 0;
    NFEBDayCount[p]       = 0;
    NFEBDDayCount[p]      = 0;
    NPIXperSRDayCount[p]  = 0;
    NPIXperEVTDayCount[p] = 0;
    NTrigDayCount[p]      = 0;
    NSpillEntryCount[0]   = 0;
    //BNBNSpillEntryCount[0]= 0;
    Ndcm0DayCount[p]      = 0;
  }

  int NSpills = 0;
  //int BNBNSpills = 0;
  // Only get Beam files for FD/ND.
  if(det_type == "FarDet" || det_type == "NearDet"){
    //Get Spill Summary
    std::cout << "\n\nChecking for Beam Summary File...\n\n";
    //  TFile* fileBeam = TFile::Open("/home/novanearline/Nearline-test-releases/S15-02-05/Commissioning/Nearline/POTBeamDay.root");
    TFile* fileBeam = TFile::Open(SpillInfoFile);
    TTree* BeamTree = dynamic_cast<TTree*>(fileBeam->Get("BeamMetrics"));
    double timestamp;
    double spills;
    BeamTree->SetBranchAddress("Time", &timestamp);
    BeamTree->SetBranchAddress("SpillRate", &spills);
    
    int n_entries = BeamTree->GetEntries();
    
    std::cout<<n_entries<<std::endl;
    //loop over all entries
    NSpillEntryCount[0] =0;
    for (int entryid = 0; entryid < n_entries; entryid++){
      BeamTree->GetEntry(entryid);
      std::cout<<entryid<<"  t="<<timestamp<<"   S= "<<spills<<std::endl;
      
      timestamp = timestamp -GMToffset;
      if ( timestamp > time_ago && timestamp < XNow ){
        SpillCountTime[0][entryid]=timestamp;
        SpillCountValue[0][entryid]=spills;
        NSpillEntryCount[0] ++;
        NSpills++;
      }
    }

    //TFile* BNBfileBeam = TFile::Open(BNBSpillInfoFile);
    //TTree* BNBBeamTree = dynamic_cast<TTree*>(BNBfileBeam->Get("BeamMetrics"));
    //double BNBtimestamp;
    //double BNBspills;
    //BNBBeamTree->SetBranchAddress("Time", &BNBtimestamp);
    //BNBBeamTree->SetBranchAddress("SpillRate", &BNBspills);
    
    //int BNB_n_entries = BNBBeamTree->GetEntries();
    
    //std::cout<<BNB_n_entries<<std::endl;
    //loop over all entries
    //BNBNSpillEntryCount[0] =0;
    //for (int entryid = 0; entryid < BNB_n_entries; entryid++){
    //BNBBeamTree->GetEntry(entryid);
    //std::cout<<entryid<<"  t="<<BNBtimestamp<<"   S= "<<BNBspills<<std::endl;
      
    //BNBtimestamp = BNBtimestamp - GMToffset;
    //if ( BNBtimestamp > time_ago && BNBtimestamp < XNow ){
    //  BNBSpillCountTime[0][entryid]=BNBtimestamp;
    //  BNBSpillCountValue[0][entryid]=BNBspills;
    //  BNBNSpillEntryCount[0] ++;
    //  BNBNSpills++;
    // }
    //}
    
  } // end Beam Summary Files
  
 
  unsigned int nFilesProcessed = 0;

  std::cout << "\n\nChecking for file...\n\n";
  while(inFile.good()) {
    if ((++nFilesProcessed % 100) == 0) std::cout << "\nProcessed " << nFilesProcessed << " files.";

    char filenameD[256];
    inFile >> filenameD;
    TFile fileD(filenameD);

    // Get subrun info out of the header.
    TTree *header  = (TTree*)fileD.Get("Header");
    if(header != 0) {
      header->SetBranchAddress("StartYear",  &year);
      header->SetBranchAddress("StartMonth", &month);
      header->SetBranchAddress("StartDay",   &day);
      header->SetBranchAddress("StartHour",  &Hour);
      header->SetBranchAddress("EndYear",    &yearEnd);
      header->SetBranchAddress("EndMonth",   &monthEnd);
      header->SetBranchAddress("EndDay",     &dayEnd);
      header->SetBranchAddress("EndHour",    &HourEnd);
      header->SetBranchAddress("Run",        &run);
      header->SetBranchAddress("Subrun",     &subrun);
      header->SetBranchAddress("Partition",  &par);
      header->GetEntry(0);

      hour = Hour;
      Min  = (Hour-hour)*60.0;
      min  = (Hour-hour)*60.0;
      sec  = (Min-min)*60.0;

      // OnMon/DAQ keeps track of GMT time, so convert to local time
      SRtime->Set(year,month,day,hour,min,sec);
      Xsrtime = SRtime->Convert() - GMToffset;

      hour = HourEnd;
      Min  = (HourEnd-hour)*60.0;
      min  = (HourEnd-hour)*60.0;
      sec  = (Min-min)*60.0;

      // Get the end time and compute subrun duration.
      SRtime->Set(yearEnd,monthEnd,dayEnd,hour,min,sec);
      if(dayEnd == day)
        SRdur = 3600.0*(HourEnd-Hour);
      else
        SRdur = 3600.0*(HourEnd+24.0-Hour);

      if(Xsrtime < time_ago) continue;

      if(run >= LastRun) {
        if(run > LastRun) {
          LastRun = run;
          LastSR  = subrun;
        }
        else if(subrun > LastSR) {
          LastRun = run;
          LastSR  = subrun;
        }
      }

    }


    // Protection against subruns with only one event in them (yes, this did happen once...)
    if(SRdur == 0.0) continue;

    std::cout << "\nLOOP:   adding run " << run << ", subrun " << subrun;

    // Get the all error plot
    TH2F *hAEtemp = (TH2F*)fileD.FindObjectAny("AEVsHour");

    // loop over the errors and fill the error vs. time plots:
    if(hAEtemp != 0 && header != 0) {
      // loop over errors (Y-axis)
      for(unsigned int e = 0; e < Nerr; ++e) {
        // loop over time bins (X-axis)
        for(unsigned int xb = 61; xb < 1500; ++xb) {
          // Histo bins start at 1 AND in this histo, I put a spacer bin between each of the error bins (to make it look better.)
          // Therefore, the bins I am interested are actually enumerated by 2*e+1
          double content = hAEtemp->GetBinContent(hAEtemp->GetBin(xb,2*e+1));

          double eventTime = hAEtemp->GetXaxis()->GetBinCenter(xb);

          int hour_fine = (int)eventTime;
          int min_fine = (int) ((eventTime - (double)hour_fine)*60.0);
          int sec_fine = (int) ((eventTime - (double)hour_fine - (double)min_fine/60.0)*3600.0);

          SRtime_fine->Set(year,month,day,hour_fine,min_fine,sec_fine);
          Xsrtime_fine = SRtime_fine->Convert() - GMToffset;

          hErrorHistos[par][e]->Fill(Xsrtime_fine,content);
        } // end loop over Xbins
      } // end loop over number of errors
    } // end check if histo and header are okay

    // Get the all alerts plot
    TH2F *hAAtemp = (TH2F*)fileD.FindObjectAny("AAVsHour");

    // loop over the alerts and fill the alert vs. time plots:
    if(hAAtemp != 0 && header != 0) {
      // loop over alerts (Y-axis)
      for(unsigned int a = 0; a < Nalert; ++a) {
        // loop over time bins (X-axis)
        for(unsigned int xb = 61; xb < 1500; ++xb) {
          // Histo bins start at 1 AND in this histo, I put a spacer bin between each of the alert bins (to make it look better.)
          // Therefore, the bins I am interested are actually enumerated by 2*a+1
          double content = hAAtemp->GetBinContent(hAAtemp->GetBin(xb,2*a+1));

          double eventTime = hAAtemp->GetXaxis()->GetBinCenter(xb);

          int hour_fine = (int)eventTime;
          int min_fine = (int) ((eventTime - (double)hour_fine)*60.0);
          int sec_fine = (int) ((eventTime - (double)hour_fine - (double)min_fine/60.0)*3600.0);

          SRtime_fine->Set(year,month,day,hour_fine,min_fine,sec_fine);
          Xsrtime_fine = SRtime_fine->Convert() - GMToffset;

          hAlertHistos[par][a]->Fill(Xsrtime_fine,content);
        } // end loop over Xbins
      } // end loop over number of alerts
    } // end check if histo and header are okay

    TH1F *hQtemp = (TH1F*)fileD.FindObjectAny("QPlotALL");
    if(hQtemp != 0)
      hQDay[par]->Add(hQtemp,1.0);

    livetime = 0.0;
    TH2F *hFEBRateTemp = (TH2F*)fileD.FindObjectAny("FEBHitMap");
    TH1F *hTime = (TH1F*)fileD.FindObjectAny("RecordedTime");
    if(hFEBRateTemp != 0 && hTime != 0) {
      livetime = (hTime->Integral())*(5.0E-7);   // recorded time is in units of 500 ns
      hFEBRateHitMapDay[par]->Add(hFEBRateTemp,1.0);
      for(int i = 1; i <= hFEBRateTemp->GetNbinsX(); ++i) {
        for(int j = 1; j <= hFEBRateTemp->GetNbinsY(); ++j) {
          double bin1 = hFEBRateTemp->GetBinContent(hFEBRateTemp->GetBin(i,j));
          double bin2 = hFEBTimeHitMapDay[par]->GetBinContent(hFEBTimeHitMapDay[par]->GetBin(i,j));
          if(bin1 > 0.0) hFEBTimeHitMapDay[par]->SetBinContent(hFEBTimeHitMapDay[par]->GetBin(i,j),bin2+livetime);
        }
      }
    }

    TH2F *hFEBDropOut = (TH2F*)fileD.FindObjectAny("FEBShutOff");
    if(hFEBDropOut != 0 && header !=0){
      hFEBDropOutMapDay[par]->Add(hFEBDropOut,1.0);
    }

    TH2F *hFEBtotADCTemp = (TH2F*)fileD.FindObjectAny("TotADCFEBHitMap");
    if(hFEBtotADCTemp != 0 && hFEBRateTemp != 0) {
      hFEBaveADCHitMapDay[par]->Add(hFEBtotADCTemp,1.0);
      hFEBhitsHitMapDay[par]->Add(hFEBRateTemp,1.0);
    }

    TH2F *hDCMHit = (TH2F*)fileD.FindObjectAny("DCMHitMap");
    if(hDCMHit != 0 && header !=0) {
    int DCMhitcount = 0;
      for(int i = 1; i <= hDCMHit->GetNbinsX(); i++){
        for(int j = 1; j <= hDCMHit->GetNbinsY(); j++){
          if(hDCMHit->GetBinContent(hDCMHit->GetBin(i,j)) > 0)
            DCMhitcount++;
        }
      }
      NDCMtimeDay[par][NDCMDayCount[par]] = Xsrtime;
      NDCMDay    [par][NDCMDayCount[par]] = DCMhitcount;
      NDCMDayCount[par]++;
    }

    TH2F  *hFEBHit = (TH2F*)fileD.FindObjectAny("FEBHitMap");
    TH2F  *hFEBRate = (TH2F*)fileD.FindObjectAny("FEBHitRateMap");
    TH2F  *hFEBRateSpec = (TH2F*)fileD.FindObjectAny("FEBHitRateSpectrum");
    if(hFEBRate != 0 && header != 0) {
      hNormalizationVsTimeDay[par]->Fill(Xsrtime);
      hTimingNormVsTimeDay[par]->Fill(Xsrtime, SRdur);
      int FEBcount = 0;
      for(int i = 1; i <= hFEBRate->GetNbinsX(); ++i) {
        for(int j = 1; j <= hFEBRate->GetNbinsY(); ++j) {
          double rate = hFEBRate->GetBinContent(hFEBRate->GetBin(i,j));
          if(rate > 0.0) {
            rate = log10(rate);
            FEBcount++;
          }
          else
            rate = -100.0;
          hFEBRateSpectrumVsTimeDay[par]->Fill(Xsrtime,rate);
        }
      }
    }
    if(hFEBHit != 0 && header != 0){
      int FEBhitcount = 0;
      for(int i = 1; i <= hFEBHit->GetNbinsX(); i++){
        for(int j = 1; j <= hFEBHit->GetNbinsY(); j++){
          if(hFEBHit->GetBinContent(hFEBHit->GetBin(i,j)) > 0)
            FEBhitcount++;
        }
      }
      NFEBDay    [par][NFEBDayCount[par]] = FEBhitcount;
      NFEBDayCount[par]++;
      NFEBtimeDay[par][NFEBDayCount[par]] = Xsrtime;
    }
    if(hFEBDropOut != 0 && header != 0){
      int FEBdropoutcount = 0;
      for(int i = 1; i <= hFEBDropOut->GetNbinsX(); i++){
        for(int j = 1; j <= hFEBDropOut->GetNbinsY(); j++){
          FEBdropoutcount += hFEBDropOut->GetBinContent(hFEBDropOut->GetBin(i,j));
        }
      }
      NFEBDDay[par][NFEBDDayCount[par]] = FEBdropoutcount;
      NFEBDDayCount[par]++;
      NFEBtimeDDay[par][NFEBDDayCount[par]] = Xsrtime;
    }
    TH1F *hTPlotALLTemp = (TH1F*)fileD.FindObjectAny("TPlotALL");
    if(hTPlotALLTemp != 0 && header != 0) {
      hTPlotALLDay[par]->Add(hTPlotALLTemp);

      for(int i = 1; i <= hTPlotALLTemp->GetNbinsX(); ++i) {
        double content = hTPlotALLTemp->GetBinContent(hTPlotALLTemp->GetBin(i));
        hTPlotALLVsTimeDay[par]->Fill(Xsrtime,hTPlotALLTemp->GetBinCenter(i),content);
      }
    }

    TH1F *hTPlotZOOMTemp = (TH1F*)fileD.FindObjectAny("TPlotZOOM");
    if(hTPlotZOOMTemp != 0 && header != 0) {
      hTPlotZOOMDay[par]->Add(hTPlotZOOMTemp);

      for(int i = 1; i <= hTPlotZOOMTemp->GetNbinsX(); ++i) {
        double content = hTPlotZOOMTemp->GetBinContent(hTPlotZOOMTemp->GetBin(i));
        hTPlotZOOMVsTimeDay[par]->Fill(Xsrtime,hTPlotZOOMTemp->GetBinCenter(i),content);
      }
    }

    TH2F *hPixelRateSpec = (TH2F*)fileD.FindObjectAny("PixelHitRateSpectrum");
    if(hPixelRateSpec != 0 && header !=0) {
      NPIXperSRtimeDay[par][NPIXperSRDayCount[par]] = Xsrtime;
      NPIXperSRDay    [par][NPIXperSRDayCount[par]] = hPixelRateSpec->Integral();
      NPIXperSRDayCount[par]++;
    }

    TH2F *hNPIXvsTime = (TH2F*)fileD.FindObjectAny("NhitVsHour");
    if(hNPIXvsTime != 0 && header != 0) {
      NPIXperEVTtimeDay[par][NPIXperEVTDayCount[par]] = Xsrtime;
      NPIXperEVTDay    [par][NPIXperEVTDayCount[par]] = hNPIXvsTime->GetMean(2);
      NPIXperEVTDayCount[par]++;
    }

    TH2F *hTrig = (TH2F*)fileD.FindObjectAny("TriggerVsHourGeneral");
    if(hTrig != 0 && header != 0) {
      double Trig00 = hTrig->Integral(61,1500,1,1);
      double Trig01 = hTrig->Integral(61,1500,2,2); //BNB @ ND, Beamline @ TB
      double Trig02 = hTrig->Integral(61,1500,3,3);
      double TrigAct= hTrig->Integral(61,1500,22,22); //DDActivity
      Trig00timeDay[par][NTrigDayCount[par]] = Xsrtime;
      Trig01timeDay[par][NTrigDayCount[par]] = Xsrtime;
      Trig02timeDay[par][NTrigDayCount[par]] = Xsrtime;
      TrigActtimeDay[par][NTrigDayCount[par]]= Xsrtime;
      Trig00Day[par][NTrigDayCount[par]]     = Trig00/SRdur;
      Trig01Day[par][NTrigDayCount[par]]     = Trig01/SRdur;
      Trig02Day[par][NTrigDayCount[par]]     = Trig02/SRdur;
      TrigActDay[par][NTrigDayCount[par]]    = TrigAct/SRdur;
      NTrigDayCount[par]++;
    }

    TH2F *hNdcm = (TH2F*)fileD.FindObjectAny("NdcmsVsHour");
    double Zerodcm = 0.0;
    if(hNdcm != 0 && header != 0) {
      int    NbinsY = hNdcm->GetNbinsY();
      double total  = hNdcm->Integral(61,1500,1,NbinsY);
      if(total > 0.0)
        Zerodcm = ((double)(hNdcm->Integral(61,1500,1,1)))/(total)*100.0;
      dcm0timeDay[par][Ndcm0DayCount[par]] = Xsrtime;
      Ndcm0Day[par][Ndcm0DayCount[par]]    = Zerodcm;
      Ndcm0DayCount[par]++;
    }

    fileD.Close();
  }

  for(unsigned int p = 0; p < Npar; ++p) {
    hFEBaveADCHitMapDay[p]->Divide(hFEBhitsHitMapDay[p]);
    hFEBRateHitMapDay[p]->Divide(hFEBTimeHitMapDay[p]);

    for(int i = 1; i <= hFEBRateHitMapDay[p]->GetNbinsX(); ++i) {
      for(int j = 1; j <= hFEBRateHitMapDay[p]->GetNbinsY(); ++j) {
        double rate = hFEBRateHitMapDay[p]->GetBinContent(hFEBRateHitMapDay[p]->GetBin(i,j));
        if(rate > 0.0)
          rate = log10(rate);
        else
          rate = -100.0;
        hFEBRateSpectrumDay[p]->Fill(rate);
      }
    }

    // DIVIDE BY NORM HISTO
    for(int i = 1; i <= hFEBRateSpectrumVsTimeDay[p]->GetNbinsX(); ++i) {
      double norm = hNormalizationVsTimeDay[p]->GetBinContent(hNormalizationVsTimeDay[p]->GetBin(i));
      for(int j = 1; j <= hFEBRateSpectrumVsTimeDay[p]->GetNbinsY(); ++j) {
        double content = 0.0;
        if(norm > 0.0) content = (hFEBRateSpectrumVsTimeDay[p]->GetBinContent(hFEBRateSpectrumVsTimeDay[p]->GetBin(i,j)))/norm;
        hFEBRateSpectrumVsTimeDay[p]->SetBinContent(hFEBRateSpectrumVsTimeDay[p]->GetBin(i,j),content);
      }
    }

    // Different norm histo for timing peak. Normalize by subrun length instead of # subruns
    // This way shorter subruns don't cause weird features in the plot
    for(int i = 1; i <= hTPlotALLVsTimeDay[p]->GetNbinsX(); ++i) {
      double norm = hTimingNormVsTimeDay[p]->GetBinContent(hTimingNormVsTimeDay[p]->GetBin(i));
      for(int j = 1; j <= hTPlotALLVsTimeDay[p]->GetNbinsY(); ++j) {
        double content = 0.0;
        if(norm > 0.0) content = (hTPlotALLVsTimeDay[p]->GetBinContent(hTPlotALLVsTimeDay[p]->GetBin(i,j)))/norm;
        hTPlotALLVsTimeDay[p]->SetBinContent(hTPlotALLVsTimeDay[p]->GetBin(i,j),content);
      }
    }

    for(int i = 1; i <= hTPlotZOOMVsTimeDay[p]->GetNbinsX(); ++i) {
      double norm = hTimingNormVsTimeDay[p]->GetBinContent(hTimingNormVsTimeDay[p]->GetBin(i));
      for(int j = 1; j <= hTPlotZOOMVsTimeDay[p]->GetNbinsY(); ++j) {
        double content = 0.0;
        if(norm > 0.0) content = (hTPlotZOOMVsTimeDay[p]->GetBinContent(hTPlotZOOMVsTimeDay[p]->GetBin(i,j)))/norm;
        hTPlotZOOMVsTimeDay[p]->SetBinContent(hTPlotZOOMVsTimeDay[p]->GetBin(i,j),content);
      }
    }

  }





  //
  // Make time/date stamp for each plot
  //
  TDatime *Ttemp = new TDatime;  // finish time
  int Xfin = Ttemp->Convert() - GMToffset;
  TDatime *Tfinish = new TDatime(Xfin);  // finish time
  TPaveText *UpdateText = new TPaveText(0.05, 0.0, 0.5, 0.06, "NDC");
  UpdateText->SetLineColor(0);
  UpdateText->SetFillColor(0);
  UpdateText->SetBorderSize(1);
  UpdateText->SetMargin(0.0);
  UpdateText->SetTextAlign(11);
  char buff1[256];
  sprintf(buff1, "Last updated on:    %s (central time)", Tfinish->AsString());
  UpdateText->AddText(buff1);
  char buff2[256];
  sprintf(buff2, "Last run / subrun:   %d / %d", LastRun, LastSR);
  UpdateText->AddText(buff2);

  TPaveText *StatusText = new TPaveText(0.1, 0.1, 0.9, 0.9, "NDC");
  StatusText->SetLineColor(1);
  StatusText->SetFillColor(0);
  StatusText->SetBorderSize(1);
  // StatusText->SetMargin(0.0);
  StatusText->AddText(buff1);
  StatusText->AddText(buff2);


  //
  // Draw pretty canvases...
  //
  int currentTime = 0;
  int threshold = 0;
  int timeVariable = 0;
  double current  = 0.0, ave = 0.0, max = 0.0, max2 = 0.0, min1 = 99999, min2 = 99998;
  TPaveText *LastPoint = new TPaveText(0.3,0.88,0.93,0.93,"NDC");
  TPaveText *WarningText = new TPaveText(0,0,1,1,"NBNDC");
  //WarningText->SetHighLightColor(kRed);
  //WarningText->SetLineColor(kRed);
  WarningText->SetTextColor(kWhite);
  WarningText->SetFillColorAlpha(kRed,0);
  WarningText->SetTextSize(0.07);
  LastPoint->SetLineColor(1);
  LastPoint->SetFillColor(0);
  LastPoint->SetBorderSize(1);
  char lptext[128];
  char warntext[128];
  bool send_spill_email=false;
  bool send_feb_email=false;
  bool send_ddactrate_email=false;
  int nfeb=0;
  float ddactrate=0;

  // Make conditional standard time axis labels depending on period (Day, Week, Month)
  string taxis_labels;
  if(period == "Day"){
    taxis_labels = "%H:%M";
  }
  else{
    taxis_labels = "%m/%d";
  }

  // make the OnMon hardware label drawer
  om::HwDetLabel *hwdetlbl = new om::HwDetLabel();

  // loop over all partitions to make 24 hour plots
  for(unsigned int p = 0; p < Npar; ++p) {
    char filename[128];
    char title[128];
    
    for(unsigned int e = 0; e < Nerr; ++e) {
      TCanvas *cErr = new TCanvas("cErr","Errors",2000,800);
      cErr->cd();

      hErrorHistos[p][e]->SetLineWidth(2);
      hErrorHistos[p][e]->SetLineColor(kBlue);
      hErrorHistos[p][e]->Draw("hist");
      UpdateText->Draw();
      hErrorHistos[p][e]->GetXaxis()->SetTimeDisplay(1);
      hErrorHistos[p][e]->GetXaxis()->SetLabelSize(0.03);
      hErrorHistos[p][e]->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
      hErrorHistos[p][e]->SetStats(kFALSE);

      sprintf(filename,"%s/%s-%s-P%.1u%s%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,errorHistoNames[e].c_str(),period.c_str());
      cErr->Print(filename);

      delete cErr;
      cErr = 0;
    }

    for(unsigned int a = 0; a < Nalert; ++a) {
      TCanvas *cAlert = new TCanvas("cAlert","Alerts",2000,800);
      cAlert->cd();

      hAlertHistos[p][a]->SetLineWidth(2);
      hAlertHistos[p][a]->SetLineColor(kBlue);
      hAlertHistos[p][a]->Draw("hist");
      UpdateText->Draw();
      hAlertHistos[p][a]->GetXaxis()->SetTimeDisplay(1);
      hAlertHistos[p][a]->GetXaxis()->SetLabelSize(0.03);
      hAlertHistos[p][a]->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
      hAlertHistos[p][a]->SetStats(kFALSE);

      sprintf(filename,"%s/%s-%s-P%.1u%s%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,alertHistoNames[a].c_str(),period.c_str());
      cAlert->Print(filename);

      delete cAlert;
      cAlert = 0;
    }


    TCanvas *cQday = new TCanvas("cQday","Qplot Day",1200,800);
    cQday->cd();
    gPad->SetLogx();
    gPad->SetLogy();
    gStyle->SetOptStat(111111);
    hQDay[p]->SetAxisRange(5.0,5000.0,"X");
    hQDay[p]->SetLineWidth(2);
    hQDay[p]->SetLineColor(kRed);
    hQDay[p]->Draw();
    UpdateText->Draw();
    // Change this to normal directory but different plot name
    sprintf(filename,"%s/%s-%s-P%.1uQPlot%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cQday->Print(filename);
    delete cQday;
    cQday = 0;

    TCanvas *cHMRday = new TCanvas("cHMRday","FEB HitMap Rates Day",1200,800);
    cHMRday->cd();
    cHMRday->SetLogz();
    gStyle->SetPalette(1); //kRainbow
    hFEBRateHitMapDay[p]->SetAxisRange(100.0,40000.0,"Z");
    hFEBRateHitMapDay[p]->GetXaxis()->SetTickLength(0);
    hFEBRateHitMapDay[p]->GetXaxis()->SetLabelColor(0);
    hFEBRateHitMapDay[p]->GetYaxis()->SetTickLength(0);
    hFEBRateHitMapDay[p]->GetYaxis()->SetLabelColor(0);
    hFEBRateHitMapDay[p]->SetStats(kFALSE);
    hwdetlbl->Config(hFEBRateHitMapDay[p]);
    hFEBRateHitMapDay[p]->Draw("colz");
    hwdetlbl->Draw();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uFEBRateHitMap%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cHMRday->Print(filename);
    delete cHMRday;
    cHMRday = 0;

    TCanvas *cHMdropday = new TCanvas("cHMdropday", "FEB Dropouts Day",1200,800);
    cHMdropday->cd();
    cHMdropday->SetLogz();
    hFEBDropOutMapDay[p]->SetAxisRange(1,100.0,"Z");
    hFEBDropOutMapDay[p]->GetXaxis()->SetTickLength(0);
    hFEBDropOutMapDay[p]->GetXaxis()->SetLabelColor(0);
    hFEBDropOutMapDay[p]->GetYaxis()->SetTickLength(0);
    hFEBDropOutMapDay[p]->GetYaxis()->SetLabelColor(0);
    hFEBDropOutMapDay[p]->SetStats(kFALSE);
    hwdetlbl->Config(hFEBDropOutMapDay[p]);
    hFEBDropOutMapDay[p]->Draw("colz");
    hwdetlbl->Draw();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uFEBDropOutMap%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cHMdropday->Print(filename);
    delete cHMdropday;
    cHMdropday = 0;

    TCanvas *cHMadcday = new TCanvas("cHMadcday","FEB ave ADC HitMap Day",1200,800);
    cHMadcday->cd();
    cHMadcday->SetLogz();
    hFEBaveADCHitMapDay[p]->SetAxisRange(0.5,1000,"Z");
    hFEBaveADCHitMapDay[p]->GetXaxis()->SetTickLength(0);
    hFEBaveADCHitMapDay[p]->GetXaxis()->SetLabelColor(0);
    hFEBaveADCHitMapDay[p]->GetYaxis()->SetTickLength(0);
    hFEBaveADCHitMapDay[p]->GetYaxis()->SetLabelColor(0);
    hFEBaveADCHitMapDay[p]->SetStats(kFALSE);
    hwdetlbl->Config(hFEBaveADCHitMapDay[p]);
    hFEBaveADCHitMapDay[p]->Draw("colz");
    hwdetlbl->Draw();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uFEBaveADCHitMap%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cHMadcday->Print(filename);
    delete cHMadcday;
    cHMadcday = 0;

    TCanvas *cFRday = new TCanvas("cFRday","FEB Rate Day",1200,800);
    cFRday->cd();
    gPad->SetLogy();
    gStyle->SetOptStat(111111);
    hFEBRateSpectrumDay[p]->SetLineWidth(2);
    hFEBRateSpectrumDay[p]->SetLineColor(kGreen+2);
    hFEBRateSpectrumDay[p]->Draw();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uFEBRateSpectrum%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cFRday->Print(filename);
    delete cFRday;
    cFRday = 0;

    TCanvas *cFEBRTday = new TCanvas("cFEBRTday","FEB Rate vs Time day",1200,800);
    cFEBRTday->cd();
    gPad->SetGridx();
    gPad->SetLogz();
    hFEBRateSpectrumVsTimeDay[p]->GetXaxis()->SetTimeDisplay(1);
    hFEBRateSpectrumVsTimeDay[p]->GetXaxis()->SetLabelSize(0.03);
    hFEBRateSpectrumVsTimeDay[p]->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    hFEBRateSpectrumVsTimeDay[p]->SetStats(kFALSE);
    hFEBRateSpectrumVsTimeDay[p]->Draw("colz");
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uFEBRateSpectrumVsTime%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cFEBRTday->Print(filename);
    delete cFEBRTday;
    cFEBRTday = 0;

    TCanvas *cTPAday = new TCanvas("cTPAday","TPlotALL Day",1200,800);
    cTPAday->cd();
    gStyle->SetOptStat(111111);
    hTPlotALLDay[p]->SetLineWidth(2);
    hTPlotALLDay[p]->SetLineColor(kBlue);
    hTPlotALLDay[p]->Draw();
    hTPlotALLDay[p]->GetYaxis()->UnZoom();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uTPlotALL%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTPAday->Print(filename);
    delete cTPAday;
    cTPAday = 0;

    TCanvas *cTPATday = new TCanvas("cTPATday","TPlotALL vs Time day",1200,800);
    cTPATday->cd();
    gPad->SetGridx();
    hTPlotALLVsTimeDay[p]->GetXaxis()->SetTimeDisplay(1);
    hTPlotALLVsTimeDay[p]->GetXaxis()->SetLabelSize(0.03);
    hTPlotALLVsTimeDay[p]->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    hTPlotALLVsTimeDay[p]->SetStats(kFALSE);
    hTPlotALLVsTimeDay[p]->Draw("colz");
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uTPlotALLVsTime%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTPATday->Print(filename);
    delete cTPATday;
    cTPATday = 0;

    TCanvas *cTPZday = new TCanvas("cTPZday","TPlotZOOM Day",1200,800);
    cTPZday->cd();
    gStyle->SetOptStat(111111);
    hTPlotZOOMDay[p]->SetLineWidth(2);
    hTPlotZOOMDay[p]->SetLineColor(kBlue);
    hTPlotZOOMDay[p]->Draw();
    hTPlotZOOMDay[p]->GetYaxis()->UnZoom();
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uTPloZOOM%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTPZday->Print(filename);
    delete cTPZday;
    cTPZday = 0;

    TCanvas *cTPZTday = new TCanvas("cTPZTday","TPlotZOOM vs Time day",1200,800);
    cTPZTday->cd();
    gPad->SetGridx();
    hTPlotZOOMVsTimeDay[p]->GetXaxis()->SetTimeDisplay(1);
    hTPlotZOOMVsTimeDay[p]->GetXaxis()->SetLabelSize(0.03);
    hTPlotZOOMVsTimeDay[p]->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    hTPlotZOOMVsTimeDay[p]->SetStats(kFALSE);
    hTPlotZOOMVsTimeDay[p]->Draw("colz");
    UpdateText->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uTPlotZOOMVsTime%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTPZTday->Print(filename);
    delete cTPZTday;
    cTPZTday = 0;


    // ........................................................................
    // Number of Active DCMs per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NDCMDayCount[p]; ++i) {
      ave += (double)NDCMDay[p][i];
      if(NDCMtimeDay[p][i] > currentTime) {
        currentTime = NDCMtimeDay[p][i];
        current     = NDCMDay[p][i];
      }
    }
    if(NDCMDayCount[p] > 0) ave = ave/(double)NDCMDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNDCMday = new TCanvas("cNDCMday","# of DCMs day",1200,800);
    cNDCMday->cd();
    gPad->SetGridx();
    TGraph *gNDCMday = new TGraph(NDCMDayCount[p],NDCMtimeDay[p],NDCMDay[p]);
    sprintf(title,"Number of Active DCMs per Subrun - partition %.1u",p);
    gNDCMday->SetTitle(title);
    gNDCMday->SetMarkerColor(kBlue);
    gNDCMday->GetXaxis()->SetTimeDisplay(1);
    gNDCMday->GetXaxis()->SetLabelSize(0.03);
    gNDCMday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNDCMday->GetXaxis()->SetLimits(time_ago,XNow);
    gNDCMday->GetXaxis()->SetTitle("(central time)");
    gNDCMday->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNDCMperSR%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNDCMday->Print(filename);
    delete cNDCMday;
    cNDCMday = 0;


    // ........................................................................
    // Number of Active FEBs per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NFEBDayCount[p]; ++i) {
      ave += (double)NFEBDay[p][i];
      if(NFEBtimeDay[p][i] > currentTime) {
        currentTime = NFEBtimeDay[p][i];
        current     = NFEBDay[p][i];
      }
    }
    if(NFEBDayCount[p] > 0) ave = ave/(double)NFEBDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNFEBday = new TCanvas("cNFEBday","# of FEBs day",1200,800);
    cNFEBday->cd();
    gPad->SetGridx();
    TGraph *gNFEBday = new TGraph(NFEBDayCount[p],NFEBtimeDay[p],NFEBDay[p]);
    sprintf(title,"%s Number of Active FEBs per Subrun - partition %.1u",trigname.c_str(), p);
    gNFEBday->SetTitle(title);

    // Change plot color to red if latest point drops below threshold
    gNFEBday->SetMarkerColor(kBlue);
    gNFEBday->GetXaxis()->SetTimeDisplay(1);
    gNFEBday->GetXaxis()->SetLabelSize(0.03);
    gNFEBday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNFEBday->GetXaxis()->SetLimits(time_ago,XNow);
    gNFEBday->GetXaxis()->SetTitle("(central time)");
    if(det_type == "FarDet"){
      gNFEBday->SetMinimum(0);
      gNFEBday->SetMaximum(11700);
    }
    double lmax=gPad->GetUymax();
    int FEBmax = 0;
    if(det_type == "FarDet"){
      FEBmax = 10749;
    }
    if(det_type == "NearDet"){
      FEBmax = 631;
    }
    if(det_type =="TestBeam"){
      FEBmax = 126;
    }
    double range = FEBmax - 620;
    double rangeBuffer = 0.10 * range;
    double thisMax = ((rangeBuffer > 1) ? FEBmax + rangeBuffer : FEBmax + 6);
    gNFEBday->SetMaximum(thisMax);
    gNFEBday->Draw("A*");
    gPad->Update();
    TLine* lMaxFEBs = new TLine(gPad->GetUxmin(),FEBmax,gPad->GetUxmax(),FEBmax);
    lMaxFEBs->SetLineColor(kGreen+2);
    lMaxFEBs->SetLineWidth(3);
    lMaxFEBs->SetLineStyle(9);
    lMaxFEBs->Draw();
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNFEBperSR%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNFEBday->Print(filename);
    //delete cNFEBday;
    //cNFEBday = 0;

    // Number of Drop Outs v. Time Plot
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NFEBDDayCount[p]; ++i) {
      ave += (double)NFEBDDay[p][i];
      if(NFEBtimeDDay[p][i] > currentTime) {
        currentTime = NFEBtimeDDay[p][i];
        current     = NFEBDDay[p][i];
      }
    }
    if(NFEBDDayCount[p] > 0) ave = ave/(double)NFEBDDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNFEBDday = new TCanvas("cNFEBDday","# of FEB Drop Outs day",1200,800);
    cNFEBDday->cd();
    gPad->SetGridx();
    TGraph *gNFEBDday = new TGraph(NFEBDDayCount[p],NFEBtimeDDay[p],NFEBDDay[p]);
    sprintf(title,"%s Number of FEB Drop Outs per Subrun - partition %.1u",trigname.c_str(), p);
    gNFEBDday->SetTitle(title);

    gNFEBDday->SetMarkerColor(kBlue);
    gNFEBDday->GetXaxis()->SetTimeDisplay(1);
    gNFEBDday->GetXaxis()->SetLabelSize(0.03);
    gNFEBDday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNFEBDday->GetXaxis()->SetLimits(time_ago,XNow);
    gNFEBDday->GetXaxis()->SetTitle("(central time)");
    gNFEBDday->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNDropOutFEBperSR%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNFEBDday->Print(filename);
    //delete cNFEBDday;
    //cNFEBDday = 0;
    // ........................................................................
    // TODO: Number of Active FEBs per Subrun TEST TEST TEST TEST
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NFEBDayCount[p]; ++i) {
      ave += (double)NFEBDay[p][i];
      if(NFEBtimeDay[p][i] > currentTime) {
        currentTime = NFEBtimeDay[p][i];
        current     = NFEBDay[p][i];
      }
    }
    if(NFEBDayCount[p] > 0) ave = ave/(double)NFEBDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNFEBdayTest = new TCanvas("cNFEBdayTest","# of FEBs day",1200,800);
    cNFEBdayTest->SetFillColor(kRed - 4);
    cNFEBdayTest->cd();
    gPad->SetGridx();
    TGraph *gNFEBdayTest = new TGraph(NFEBDayCount[p],NFEBtimeDay[p],NFEBDay[p]);
    sprintf(title,"%s Number of Active FEBs per Subrun - partition %.1u",trigname.c_str(), p);
    gNFEBdayTest->SetTitle(title);

    // Change plot color to red if latest point drops below threshold
    if((det_type == "FarDet") && (current < 10690)) {
      //      gNFEBdayTest->SetMarkerColor(kRed);
      // gPad->SetFrameFillColor(kWhite);
      // cNFEBdayTest->SetFillColor(kRed);
      gNFEBdayTest->SetMarkerColor(kBlue);
    }
    else {
      //      gNFEBdayTest->SetMarkerColor(kWhite);
    }
    gNFEBdayTest->GetXaxis()->SetTimeDisplay(1);
    gNFEBdayTest->GetXaxis()->SetLabelSize(0.03);
    gNFEBdayTest->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNFEBdayTest->GetXaxis()->SetLimits(time_ago,XNow);
    gNFEBdayTest->GetXaxis()->SetTitle("(central time)");
    if(det_type == "FarDet"){
      gNFEBdayTest->SetMinimum(10650);
      gNFEBdayTest->SetMaximum(10720);
    }
    gNFEBdayTest->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNFEBperSR%s-TEST.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNFEBdayTest->Print(filename);
    delete cNFEBdayTest;
    cNFEBdayTest = 0;
    //////////////////////////////////////////////////////////////////////////////////



    // ........................................................................
    // Number of Active FEBs per Subrun (ZOOM)
    if(det_type != "NearDet"){
      timeVariable = 0;
      currentTime = 0;
      current = 0.0;
      ave = 0.0;
      max = 0.0;
      max2 = 0.0;
      min1 = 99999;
      min2 = 99998;
      for(int i = 0; i < NFEBDayCount[p]; ++i) {
        ave += (double)NFEBDay[p][i];
        if(i==0)
          timeVariable = NFEBtimeDay[p][i];
        if(NFEBtimeDay[p][i] > currentTime) {
          currentTime = NFEBtimeDay[p][i];
          current     = NFEBDay[p][i];
        }
        //std::cout << "i, Time: " << i << ", " << NFEBtimeDay[p][i] << std::endl;
        //find min 2 points
        if(NFEBDay[p][i] <= min1){
          //if zero data points are over x seconds old, ignore them (to avoid always zooming out)
          if((NFEBDay[p][i] == 0 ) && ((timeVariable - NFEBtimeDay[p][i]) > 7400))
            continue;
          min2 = min1;
          min1 = NFEBDay[p][i];
        }
        if(NFEBDay[p][i] > min1){
          if(NFEBDay[p][i] < min2)
            min2 = NFEBDay[p][i];
        }

        //find max 2 points
        if(NFEBDay[p][i] >= max){
          max2 = max;
          max = NFEBDay[p][i];
        }
        if(NFEBDay[p][i] < max){
          if(NFEBDay[p][i] > max2)
            max2 = NFEBDay[p][i];
        }
        //std::cout << "Current, Max, max2, min, min2: " << NFEBDay[p][i] << ", " << max << ", " << max2 << ", " << min1 << ", " << min2 << std::endl;
      }
      //threshold for showing warning error for NFEB plot
      threshold = max2 - 45;
      //max2 += (max2-min2)*0.1;
      //min2 -= (max2-min2)*0.1;
      if(NFEBDayCount[p] > 0) ave = ave/(double)NFEBDayCount[p];
      sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
      LastPoint->Clear();
      LastPoint->AddText(lptext);
      TCanvas *cNFEBday_zoom = new TCanvas("cNFEBday_zoom","# of FEBs day",1200,800);
      cNFEBday_zoom->cd();
      gPad->SetGridx();
      TGraph *gNFEBday_zoom = new TGraph(NFEBDayCount[p],NFEBtimeDay[p],NFEBDay[p]);
      sprintf(title,"%s Number of Active FEBs per Subrun - partition %.1u",trigname.c_str(), p);
      gNFEBday_zoom->SetTitle(title);
      //std::cout << "Current threshold: " << threshold << std::endl;
      // Change plot color to red if latest point drops below threshold
      int FEBmax = 0;
      if(det_type == "FarDet"){
        FEBmax = 10749;
      }
      if(det_type == "NearDet"){
        FEBmax = 631;
      }
      if(det_type == "TestBeam"){
        FEBmax = 126;
        threshold = 122;
      }
      if((det_type == "FarDet" || det_type == "TestBeam") && (current < threshold) && (period == "Day")) {
        cNFEBday_zoom->SetFillColor(kRed);
        sprintf(warntext,"WARNING! Low Number of FEBs! (%.f)",current);
        WarningText->Clear();
        WarningText->SetFillColorAlpha(kRed,0.5);
        WarningText->AddText(warntext);
        if (p==1 && det_type == "TestBeam") {
          send_feb_email=true;
          nfeb = current;
        }
      }
      else if((det_type == "FarDet") && (current < 10500)){
        gNFEBday_zoom->SetMarkerColor(kWhite);
        cNFEBday_zoom->SetFillColor(kRed);
        sprintf(warntext,"WARNING! Low Number of FEBs! (%.f)",current);
        WarningText->Clear();
        WarningText->AddText(warntext);
      }
      else {
        gNFEBday_zoom->SetMarkerColor(kBlue);
      }
      gNFEBday_zoom->GetXaxis()->SetTimeDisplay(1);
      gNFEBday_zoom->GetXaxis()->SetLabelSize(0.03);
      gNFEBday_zoom->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
      gNFEBday_zoom->GetXaxis()->SetLimits(time_ago,XNow);
      gNFEBday_zoom->GetXaxis()->SetTitle("(central time)");
      if(period == "Month" || period == "Week"){
        gNFEBday_zoom->SetMinimum(ave-105);
      }
      else{
        gNFEBday_zoom->SetMinimum(min2-(max2-min2)*0.1);
      }
      double range = FEBmax - gNFEBday_zoom->GetMinimum();
      double rangeBuffer = 0.10 * range;
      double thisMax = ((rangeBuffer > 1) ? FEBmax + rangeBuffer : FEBmax + 1);
      gNFEBday_zoom->SetMaximum(thisMax);
      gNFEBday_zoom->Draw("A*");
      cNFEBday_zoom->Update();
      TLine* lMaxFEBs = new TLine(gPad->GetUxmin(),FEBmax,gPad->GetUxmax(),FEBmax);
      lMaxFEBs->SetLineColor(kGreen+2);
      lMaxFEBs->SetLineWidth(3);
      lMaxFEBs->SetLineStyle(9);
      lMaxFEBs->Draw();
      UpdateText->Draw();
      LastPoint->Draw();
      if(((det_type == "FarDet" || det_type == "TestBeam") && (current < threshold) && (period == "Day")) || ((det_type == "FarDet") && (current < 10500) && (period == "Day") )){
        WarningText->Draw();
      }
      sprintf(filename,"%s/%s-%s-P%.1uNFEBperSRzoom%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
      cNFEBday_zoom->Print(filename);
      delete cNFEBday_zoom;
      cNFEBday_zoom = 0;
    }


    // ........................................................................
    // Number of Active Channels per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NPIXperSRDayCount[p]; ++i) {
      ave += (double)NPIXperSRDay[p][i];
      if(NPIXperSRtimeDay[p][i] > currentTime) {
        currentTime = NPIXperSRtimeDay[p][i];
        current     = NPIXperSRDay[p][i];
      }
    }
    if(NPIXperSRDayCount[p] > 0) ave = ave/(double)NPIXperSRDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNPIXperSRday = new TCanvas("cNPIXperSRday","# of channels per subrun day",1200,800);
    cNPIXperSRday->cd();
    gPad->SetGridx();
    TGraph *gNPIXperSRday = new TGraph(NPIXperSRDayCount[p],NPIXperSRtimeDay[p],NPIXperSRDay[p]);
    sprintf(title,"%s Number of Active Channels per Subrun - partition %.1u",trigname.c_str(), p);
    gNPIXperSRday->SetTitle(title);
    gNPIXperSRday->SetMarkerColor(kBlue);
    gNPIXperSRday->GetXaxis()->SetTimeDisplay(1);
    gNPIXperSRday->GetXaxis()->SetLabelSize(0.03);
    gNPIXperSRday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNPIXperSRday->GetXaxis()->SetLimits(time_ago,XNow);
    gNPIXperSRday->GetXaxis()->SetTitle("(central time)");
    gNPIXperSRday->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNPIXperSR%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNPIXperSRday->Print(filename);
    delete cNPIXperSRday;
    cNPIXperSRday = 0;


    // ........................................................................
    // Average Number of Hits per Spill
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NPIXperEVTDayCount[p]; ++i) {
      ave += (double)NPIXperEVTDay[p][i];
      if(NPIXperEVTtimeDay[p][i] > currentTime) {
        currentTime = NPIXperEVTtimeDay[p][i];
        current     = NPIXperEVTDay[p][i];
      }
    }
    if(NPIXperEVTDayCount[p] > 0) ave = ave/(double)NPIXperEVTDayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNPIXperEVTday = new TCanvas("cNPIXperEVTday","# of hits per spill day",1200,800);
    cNPIXperEVTday->cd();
    gPad->SetGridx();
    TGraph *gNPIXperEVTday = new TGraph(NPIXperEVTDayCount[p],NPIXperEVTtimeDay[p],NPIXperEVTDay[p]);
    sprintf(title,"%s Average Number of Hits per Spill - partition %.1u",trigname.c_str(),p);
    gNPIXperEVTday->SetTitle(title);
    gNPIXperEVTday->SetMarkerColor(kBlue);
    gNPIXperEVTday->GetXaxis()->SetTimeDisplay(1);
    gNPIXperEVTday->GetXaxis()->SetLabelSize(0.03);
    gNPIXperEVTday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNPIXperEVTday->GetXaxis()->SetLimits(time_ago,XNow);
    gNPIXperEVTday->GetXaxis()->SetTitle("(central time)");
    gNPIXperEVTday->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNPIXperEVT%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNPIXperEVTday->Print(filename);
    delete cNPIXperEVTday;
    cNPIXperEVTday = 0;


    // ........................................................................
    // Percentage of Empty Spills
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < Ndcm0DayCount[p]; ++i) {
      ave += (double)Ndcm0Day[p][i];
      if(dcm0timeDay[p][i] > currentTime) {
        currentTime = dcm0timeDay[p][i];
        current     = Ndcm0Day[p][i];
      }
    }
    if(Ndcm0DayCount[p] > 0) ave = ave/(double)Ndcm0DayCount[p];
    sprintf(lptext,"Last Point = %.f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cNdcm0day = new TCanvas("cNdcm0day","# of zero DCM spills day",1200,800);
    cNdcm0day->cd();
    gPad->SetGridx();
    TGraph *gNdcm0day = new TGraph(Ndcm0DayCount[p],dcm0timeDay[p],Ndcm0Day[p]);
    sprintf(title,"%s Percentage of Empty Spills (zero DCMs reporting) - partition %.1u",trigname.c_str(),p);
    gNdcm0day->SetTitle(title);
    gNdcm0day->SetMarkerColor(kBlue);
    gNdcm0day->GetXaxis()->SetTimeDisplay(1);
    gNdcm0day->GetXaxis()->SetLabelSize(0.03);
    gNdcm0day->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gNdcm0day->GetXaxis()->SetLimits(time_ago,XNow);
    gNdcm0day->GetYaxis()->SetTitle("percent");
    gNdcm0day->GetXaxis()->SetTitle("(central time)");
    gNdcm0day->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNdcm0%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cNdcm0day->Print(filename);
    delete cNdcm0day;
    cNdcm0day = 0;


    // ........................................................................
    // Average NuMI Trigger Rate Seen by OnMon per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NTrigDayCount[p]; ++i) {
      ave += (double)Trig00Day[p][i];
      if(Trig00timeDay[p][i] > currentTime) {
        currentTime = Trig00timeDay[p][i];
        current     = Trig00Day[p][i];
      }
    }
    if(NTrigDayCount[p] > 0) ave = ave/(double)NTrigDayCount[p];
    sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cTrig00day = new TCanvas("cTrig00day","trig00 rate day",1200,800);
    cTrig00day->cd();
    gPad->SetGridx();
    TGraph *gTrig00day = new TGraph(NTrigDayCount[p],Trig00timeDay[p],Trig00Day[p]);
    sprintf(title,"Average NuMI Trigger Rate (Hz) Seen by OnMon per Subrun - partition %.1u",p);
    gTrig00day->SetTitle(title);
    gTrig00day->SetMarkerColor(kBlue);
    gTrig00day->GetXaxis()->SetTimeDisplay(1);
    gTrig00day->GetXaxis()->SetLabelSize(0.03);
    gTrig00day->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gTrig00day->GetXaxis()->SetLimits(time_ago,XNow);
    gTrig00day->GetYaxis()->SetLimits(0,11);
    gTrig00day->GetXaxis()->SetTitle("(central time)");
    gTrig00day->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uNuMITriggerRate%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTrig00day->Print(filename);


    // ........................................................................
    // Average BNB/Beamline(TB) Trigger Rate Seen by OnMon per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NTrigDayCount[p]; ++i) {
      ave += (double)Trig01Day[p][i];
      if(Trig01timeDay[p][i] > currentTime) {
        currentTime = Trig01timeDay[p][i];
        current     = Trig01Day[p][i];
      }
    }
    if(NTrigDayCount[p] > 0) ave = ave/(double)NTrigDayCount[p];
    sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cTrig01day = new TCanvas("cTrig01day","trig01 rate day",1200,800);
    cTrig01day->cd();
    gPad->SetGridx();
    TGraph *gTrig01day = new TGraph(NTrigDayCount[p],Trig01timeDay[p],Trig01Day[p]);
    sprintf(title,"Average Beamline Trigger Rate (Hz) Seen by OnMon per Subrun - partition %.1u",p);
    gTrig01day->SetTitle(title);
    gTrig01day->SetMarkerColor(kBlue);
    gTrig01day->GetXaxis()->SetTimeDisplay(1);
    gTrig01day->GetXaxis()->SetLabelSize(0.03);
    gTrig01day->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gTrig01day->GetXaxis()->SetLimits(time_ago,XNow);
    gTrig01day->GetYaxis()->SetLimits(0,11);
    gTrig01day->GetXaxis()->SetTitle("(central time)");
    gTrig01day->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uBeamlineTriggerRate%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTrig01day->Print(filename);


    // ........................................................................
    // Average Cosmic Trigger Rate Seen by OnMon per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NTrigDayCount[p]; ++i) {
      ave += (double)Trig02Day[p][i];
      if(Trig02timeDay[p][i] > currentTime) {
        currentTime = Trig02timeDay[p][i];
        current     = Trig02Day[p][i];
      }
    }
    if(NTrigDayCount[p] > 0) ave = ave/(double)NTrigDayCount[p];
    sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cTrig02day = new TCanvas("cTrig02day","trig02 rate day",1200,800);
    cTrig02day->cd();
    gPad->SetGridx();
    TGraph *gTrig02day = new TGraph(NTrigDayCount[p],Trig02timeDay[p],Trig02Day[p]);
    sprintf(title,"Average Cosmic Trigger Rate (Hz) Seen by OnMon per Subrun - partition %.1u",p);
    gTrig02day->SetTitle(title);
    gTrig02day->SetMarkerColor(kBlue);
    gTrig02day->SetMarkerStyle(21);
    gTrig02day->GetXaxis()->SetTimeDisplay(1);
    gTrig02day->GetXaxis()->SetLabelSize(0.03);
    gTrig02day->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gTrig02day->GetXaxis()->SetLimits(time_ago,XNow);
    gTrig02day->GetYaxis()->SetLimits(0,1.1);
    gTrig02day->GetXaxis()->SetTitle("(central time)");
    gTrig02day->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uCosmicTriggerRate%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTrig02day->Print(filename);


    // ........................................................................
    // Average DDActivity Trigger Rate Seen by OnMon per Subrun
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NTrigDayCount[p]; ++i) {
      ave += (double)TrigActDay[p][i];
      if(TrigActtimeDay[p][i] > currentTime) {
        currentTime = TrigActtimeDay[p][i];
        current     = TrigActDay[p][i];
      }
    }
    if(NTrigDayCount[p] > 0) ave = ave/(double)NTrigDayCount[p];
    sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cTrigActday = new TCanvas("cTrigActday","trigact rate day",1200,800);
    cTrigActday->cd();
    gPad->SetGridx();
    TGraph *gTrigActday = new TGraph(NTrigDayCount[p],TrigActtimeDay[p],TrigActDay[p]);
    sprintf(title,"Average DDActivity Trigger Rate (Hz) Seen by OnMon per Subrun - partition %.1u",p);
    gTrigActday->SetTitle(title);
    
    // Change plot color to red if latest point drops below threshold
    if((det_type == "TestBeam") && (current < 10) && (period == "Day")) {
      cTrigActday->SetFillColor(kRed);
      sprintf(warntext,"WARNING! Trigger rate low! (%.2f Hz)",current);
      WarningText->Clear();
      WarningText->SetFillColorAlpha(kRed,0.5);
      WarningText->AddText(warntext);
      if (p==1) {
        send_ddactrate_email=true;
        ddactrate=current;
      }
    }
    gTrigActday->SetMarkerColor(kBlue);
    gTrigActday->GetXaxis()->SetTimeDisplay(1);
    gTrigActday->GetXaxis()->SetLabelSize(0.03);
    gTrigActday->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gTrigActday->GetXaxis()->SetLimits(time_ago,XNow);
    gTrigActday->GetYaxis()->SetLimits(0,11);
    gTrigActday->GetXaxis()->SetTitle("(central time)");
    gTrigActday->Draw("A*");
    gPad->Update();
    UpdateText->Draw();
    LastPoint->Draw();
    if((det_type == "TestBeam") && (current < 10) && (period == "Day")) {
      WarningText->Draw();
    }
    sprintf(filename,"%s/%s-%s-P%.1uActTriggerRate%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTrigActday->Print(filename);


    // ........................................................................
    // Spill Count per bin interval
    std::cout<<NSpillEntryCount[0]<<"......"<<std::endl;
    std::cout<<NSpills<<"......"<<std::endl;

    TGraph *gSpillCount = new TGraph(NSpillEntryCount[0],SpillCountTime[0],SpillCountValue[0]);

    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    for(int i = 0; i < NSpillEntryCount[0]; ++i) {
      ave += (double)SpillCountValue[0][i];
      if(SpillCountTime[0][i] > currentTime) {
        currentTime = SpillCountTime[0][i];
        current     = SpillCountValue[0][i];
      }
    }
    if(NSpillEntryCount[0] > 0) ave = ave/(double)NSpillEntryCount[0];
    sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    LastPoint->Clear();
    LastPoint->AddText(lptext);
    TCanvas *cSpill = new TCanvas("cSpill","Spill Count day",1200,800);
    cSpill->cd();
    gPad->SetGridx();
    //    TGraph *gTrig00day = new TGraph(NTrigDayCount[p],Trig00timeDay[p],Trig00Day[p]);
    sprintf(title,"Spill count per bin interval");
    gSpillCount->SetTitle(title);
    gSpillCount->SetMarkerColor(kMagenta);
    gSpillCount->SetMarkerStyle(22);
    gSpillCount->GetXaxis()->SetTimeDisplay(1);
    gSpillCount->GetXaxis()->SetLabelSize(0.03);
    gSpillCount->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    gSpillCount->GetXaxis()->SetLimits(time_ago,XNow);
    gSpillCount->GetXaxis()->SetTitle("(central time)");
    gSpillCount->Draw("A*");
    UpdateText->Draw();
    LastPoint->Draw();
    sprintf(filename,"%s/%s-%s-P%.1uSpillCount%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cSpill->Print(filename);
    
    // ........................................................................
    //     BNB Spill Count per bin interval
    //std::cout<<BNBNSpillEntryCount[0]<<"......"<<std::endl;
    //std::cout<<BNBNSpills<<"......"<<std::endl;

    //    TGraph *gBNBSpillCount = new TGraph(BNBNSpillEntryCount[0],BNBSpillCountTime[0],BNBSpillCountValue[0]);
    
    currentTime = 0;
    current = 0.0;
    ave = 0.0;
    //for(int i = 0; i < BNBNSpillEntryCount[0]; ++i) {
    //ave += (double)BNBSpillCountValue[0][i];
    //if(BNBSpillCountTime[0][i] > currentTime) {
    //  currentTime = BNBSpillCountTime[0][i];
    //  current     = BNBSpillCountValue[0][i];
    //}
    //}
    //if(BNBNSpillEntryCount[0] > 0) ave = ave/(double)BNBNSpillEntryCount[0];
    //sprintf(lptext,"Last Point = %f  /  Average = %f",current,ave);
    //LastPoint->Clear();
    //LastPoint->AddText(lptext);
    //TCanvas *cBNBSpill = new TCanvas("cBNBSpill","BNB Spill Count day",1200,800);
    //cBNBSpill->cd();
    //gPad->SetGridx();
    sprintf(title,"Spill count per bin interval");
    //gBNBSpillCount->SetTitle(title);
    //gBNBSpillCount->SetMarkerColor(kMagenta);
    //gBNBSpillCount->SetMarkerStyle(22);
    //gBNBSpillCount->GetXaxis()->SetTimeDisplay(1);
    //gBNBSpillCount->GetXaxis()->SetLabelSize(0.03);
    //gBNBSpillCount->GetXaxis()->SetTimeFormat(taxis_labels.c_str());
    //gBNBSpillCount->GetXaxis()->SetLimits(time_ago,XNow);
    //gBNBSpillCount->GetXaxis()->SetTitle("(central time)");
    //gBNBSpillCount->Draw("A*");
    //UpdateText->Draw();
    //LastPoint->Draw();
    //sprintf(filename,"%s/%s-%s-P%.1uBNBSpillCount%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    //    cBNBSpill->Print(filename);

    // Combined Active FEBs and FEB Drop Outs

    TCanvas *cCombo = new TCanvas("cCombo", "feb active/dropout day", 1200, 800);
    TPad *activepad = new TPad("activepad","",0,0,1,1);
    TPad *dropoutpad = new TPad("dropoutpad","",0,0,1,1);

    dropoutpad->SetFillStyle(4000); //transparent
    dropoutpad->SetFillColor(0); 
    dropoutpad->SetFrameFillStyle(4000);
    activepad->Draw();
    activepad->cd();
    activepad->SetGridx();

    sprintf(title,"Number of Active and Drop Out FEBs - partition %.1u",p);
    gNFEBday->SetTitle(title);
    gNFEBday->SetMarkerStyle(22);
    gNFEBday->SetMarkerColor(kBlue);
    gNFEBday->SetFillColor(kBlue);
    //gNFEBday->GetYaxis()->SetLabelOffset(1);
    gNFEBday->GetYaxis()->SetLabelColor(kBlue);
    //gNFEBday->GetYaxis()->SetNdivisions(0, kFALSE);
    if( det == "FarDet")
        gNFEBday->GetYaxis()->SetRangeUser(10740,10750);
    if( det == "NearDet")
        gNFEBday->GetYaxis()->SetRangeUser(625,635);
    gNFEBday->Draw("AP");

    activepad->Modified();
    activepad->Update();
 
    dropoutpad->Draw();
    dropoutpad->cd();
    dropoutpad->SetGridx();

    gNFEBDday->SetTitle("");
    gNFEBDday->SetMarkerStyle(23);
    gNFEBDday->SetMarkerColor(kOrange+7);
    gNFEBDday->SetFillColor(kOrange+7);
    gNFEBDday->GetYaxis()->SetLabelColor(kOrange+7);
    gNFEBDday->GetYaxis()->SetRangeUser(0,10);
    gNFEBDday->Draw("APY+");
    
    UpdateText->Draw();
    
    dropoutpad->Modified();
    dropoutpad->Update();

    TLegend *legc = new TLegend(0.62,0.69,0.87,0.78);
    legc->SetFillColor(0);
    TLegendEntry *af = legc->AddEntry(gNFEBday,"#bf{Active FEBs}","F");
    af->SetTextColor(kBlue);
    TLegendEntry *df = legc->AddEntry(gNFEBDday,"#bf{FEB Drop Outs}","F");
    df->SetTextColor(kOrange+7);
    legc->SetFillStyle(4000);
    legc->SetLineColor(0);
    legc->Draw();

    sprintf(filename,"%s/%s-%s-P%.1uNDropOutActiveFEBperSR%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cCombo->Print(filename);
    delete cCombo;
    delete cNFEBDday;
    delete cNFEBday;
    cCombo = 0;
    cNFEBDday = 0;
    cNFEBday  = 0;
    
    // ........................................................................
    // Average Trigger and Spill Rates
    gTrig00day->SetMarkerColor(kRed);
    TCanvas *cTrig = new TCanvas("cTrig", "trig rate day",1200, 800);
    
    float currentNuMITrigRate = Trig00Day[1][NTrigDayCount[1] - 1];
    float currentSpillRate    = SpillCountValue[0][NSpillEntryCount[0] - 1];
    
    //float spillRateThreshold = 0.8;
    float triggerRateThreshold = 0.9;
    float spillNumiDiffThreshold = 0.05;
    
    TPad *pad1 = new TPad("pad1","",0,0,1,1);
    TPad *pad2 = new TPad("pad2","",0,0,1,1);
    TPad *pad3 = new TPad("pad3","",0,0,1,1);
    TPad *pad4 = new TPad("pad4","",0,0,1,1);
    // Warning: spill rate too high
    // NOTE: This is not working right now, so commenting it out.
    // Hopefully someone will come back and fix this someday.
    //std::cout << "CHECKING FOR HIGH SPILL RATE......";
    //if (currentSpillRate > spillRateThreshold){
    //  std::cout << "BAD" << std::endl;
    //  pad1->SetFillColor(kRed - 4);
    //  
    //  // Email experts
    //  char cmd[1000];
    //  char to[1000] = "nova_watchdogs@fnal.gov nova_nearlineonmon_support@fnal.gov jvasel@indiana.edu";
    //  char subject[1000] = "Nearline--Current spill rate exceeds alarm thresholds";
    //  char body[1000] = "There may be an issue with the timing system.";
    //  
    //  sprintf(cmd, "echo %s | mail -s %s %s", body, subject, to);
    //  system(cmd);
    //}
    //else {
    //  std::cout << "OK" << std::endl;
    //}
    
    // Warning: spill and trigger rate differ
    //std::cout << "CHECKING FOR SPILL/TRIGGER RATE MISMATCH......";
    //if (std::abs(currentSpillRate - currentNuMITrigRate) > spillNumiDiffThreshold){
    //  std::cout << "BAD";
    //  pad1->SetFillColor(kRed - 4);
    //  
    //  // Email experts
    //  char cmd[1000];
    //  char to[1000] = "nova_watchdogs@fnal.gov nova_nearlineonmon_support@fnal.gov jvasel@indiana.edu";
    //  char subject[1000] = "Nearline--Current spill rate differs from NuMI trigger rate";
    //  char body[1000] = "There may be an issue with the spill server.";
    //  
    //  sprintf(cmd, "echo %s | mail -s %s %s", body, subject, to);
    //  system(cmd);
    //}
    //else {
    //  std::cout << "OK";
    //}
    //std::cout << " (SpillRate=" << currentSpillRate << " TrigRate=" << currentNuMITrigRate
    //          << " Diff=" << std::abs(currentSpillRate - currentNuMITrigRate) << ")" << std::endl;
    pad2->SetFillStyle(4000);  // pad2 will be transparent
    pad2->SetFillColor(0);
    pad2->SetFrameFillStyle(4000);
    pad3->SetFillStyle(4000);  // pad2 will be transparent
    pad3->SetFillColor(0);
    pad3->SetFrameFillStyle(4000);
    pad4->SetFillStyle(4000);  // pad2 will be transparent
    pad4->SetFillColor(0);
    pad4->SetFrameFillStyle(4000);
    pad1->Draw();
    pad1->cd();
    pad1->SetGridx();

    sprintf(title,"Average Trigger and Spill Rates (Hz) - partition %.1u",p);
    gTrig00day->SetTitle(title);
    gTrig00day->SetMarkerStyle(21);
    gTrig00day->SetMarkerColor(kRed+3);
    gTrig00day->SetFillColor(kRed+3);
    if (t_prd != "Day")     gTrig00day->SetMarkerStyle(7);
    gTrig00day->GetYaxis()->SetLabelOffset(1);
    gTrig00day->GetYaxis()->SetLabelColor(kRed+3);
    gTrig00day->GetYaxis()->SetNdivisions(0, kFALSE); // Disable tick marks
    gTrig00day->GetYaxis()->SetRangeUser(0,1.1);

    WarningText->Clear();
    std::cout << "CHECKING FOR HIGH NUMI TRIGGER RATE......";
    std::cout<<"Current NuMI trigger rate = "<<currentNuMITrigRate<<std::endl;
    if (currentNuMITrigRate > triggerRateThreshold){
      std::cout << "NUMI TRIGGER RATE TOO HIGH" << std::endl;
      pad1->SetFillColor(kRed);
      sprintf(warntext,"WARNING! NuMI trigger rate too high! Call DAQ expert!");
      WarningText->AddText(warntext);
      send_spill_email=true;
    }
    else {
      std::cout << "OK" << std::endl;
    }
    std::cout << " (SpillRate=" << currentSpillRate << " TrigRate=" << currentNuMITrigRate
              << " Diff=" << std::abs(currentSpillRate - currentNuMITrigRate) << ")" << std::endl;

    gTrig00day->Draw("AP");
    WarningText->Draw();

    pad1->Modified();
    pad1->Update();

    pad4->Draw();
    pad4->cd();
    pad4->SetGridx();

    //pad2->Draw();
    //pad2->cd();
    //pad2->SetGridx();

    //gBNBSpillCount->SetTitle("");
    //gBNBSpillCount->GetYaxis()->SetRangeUser(0,30);
    //gBNBSpillCount->GetYaxis()->SetLabelColor(kBlue+3);
    //gBNBSpillCount->SetMarkerStyle(34);
    //gBNBSpillCount->SetFillColor(kBlue+3);
    //gBNBSpillCount->SetMarkerColor(kBlue+3);
    //gBNBSpillCount->Draw("AP");

    if(det_type != "NearDet"){
      gTrig02day->SetTitle("");
      gTrig02day->SetMarkerStyle(21);
      if (t_prd != "Day")     gTrig02day->SetMarkerStyle(7);
      gTrig02day->SetFillColor(kBlue);
      gTrig02day->GetYaxis()->SetLabelColor(kBlue);
      gTrig02day->GetYaxis()->SetRangeUser(0,30);
      gTrig02day->Draw("AP");
    }
    pad4->Modified();
    pad4->Update();

    pad3->Draw();
    pad3->cd();
    pad3->SetGridx();

    gSpillCount->SetTitle("");
    gSpillCount->GetYaxis()->SetRangeUser(0,1.1);
    gSpillCount->GetYaxis()->SetLabelColor(kRed);
    gSpillCount->SetMarkerStyle(34);
    gSpillCount->SetFillColor(kRed);
    gSpillCount->SetMarkerColor(kRed);
    gSpillCount->Draw("APY+");
    UpdateText->Draw();

    pad3->Modified();
    pad3->Update();

    TLegend *leg = new TLegend(0.62,0.69,0.87,0.78);
    leg->SetFillColor(0);
    if(det_type != "NearDet"){
      TLegendEntry *l2 = leg->AddEntry(gTrig02day,"#bf{Cosmic Trigger Rate}","F");
      l2->SetTextColor(kBlue);
    }
    //    TLegendEntry *l4 = leg->AddEntry(gBNBSpillCount,"#bf{Delivered BNB Spill Rate}","F");
    //l4->SetTextColor(kBlue+3);
    TLegendEntry *l1 = leg->AddEntry(gTrig00day,"#bf{NuMI Trigger Rate}","F");
    l1->SetTextColor(kRed+3);
    TLegendEntry *l3 = leg->AddEntry(gSpillCount,"#bf{Delivered NuMI Spill Rate}","F");
    l3->SetTextColor(kRed);
    leg->SetFillStyle(4000);
    leg->SetLineColor(0);
    leg->Draw();

    sprintf(filename,"%s/%s-%s-P%.1uCosmicNuMITriggerRate%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),p,period.c_str());
    cTrig->Print(filename);
    delete cTrig;
    delete cTrig00day;
    delete cTrig02day;
    delete cSpill;
    cSpill=0;
    cTrig02day = 0;
    cTrig00day = 0;
    cTrig      = 0;
  }

  // If NuMI trigger rate is too high for any partition, send email.                                      
  if (send_spill_email) {
    char cmd[1000];
    char to[1000] = "tlackey@indiana.edu nova_runcoordinator@fnal.gov mcgroh@iu.edu jvasel@iu.edu nova_daq_support@fnal.gov";
    char subject[1000] = "NuMI Trigger rate too high.";
    char body[1000] = "NuMI Trigger rate is too high. There may be an issue with the spill server.";

    sprintf(cmd, "echo -e \"%s\" | mail -s \"%s\" %s",body,subject,to);
    system(cmd);
  }

  // If Number of active FEBs too low, send email.
  if (send_feb_email) {
    char cmd[1000];
    char to[1000] = "tlackey@indiana.edu";
    char subject[1000] = "NFEBs too low.";
    char body[1000];// = "Number of Active FEBs is too low. The TDU may be disconnected.";
    sprintf(body,"Number of Active FEBs (%d) is too low. The TDU may be disconnected.",nfeb);

    sprintf(cmd, "echo -e \"%s\" | mail -s \"%s\" %s",body,subject,to);
    system(cmd);
  }
  // If DDActivity rate is too low, send email.
  if (send_ddactrate_email) {
    char cmd[1000];
    char to[1000] = "tlackey@indiana.edu";
    char subject[1000] = "DDActivity rate too low.";
    char body[1000];// = "DDActivity rate is too low. The run may need restarting.";
    sprintf(body,"DDActivity rate is too low: %f Hz. The run may need restarting.",ddactrate);

    sprintf(cmd, "echo -e \"%s\" | mail -s \"%s\" %s",body,subject,to);
    system(cmd);
  }
}