DCSPlotMaker.C

Go to the documentation of this file.

#include <TFile.h>
#include <TCanvas.h>
#include <TTree.h>
#include <TGraph.h>
#include <TDatime.h>
#include <TPaveText.h>
#include <TAxis.h>
#include <TGaxis.h>
#include <TStyle.h>
#include <TObject.h>
#include <TROOT.h>

#include <algorithm>
#include <string>
#include <iostream>
     
//TCanvas *c0 = new TCanvas("c0");
TPaveText *UpdateText = new TPaveText(0.1, 0.0, 0.5, 0.05, "NDC");
TPaveText *WarningText = new TPaveText(0.15, 0.15, 0.9, 0.3, "NDC");
int nwarnings = 0;
//Plots a single channel onto the canvas
void plot_with_zoom(TTree *t, char *branchname, int color, int marker, TCanvas *c1, char *title, float *ylimit,
                    float *zoom, UInt_t timeago, UInt_t XNow, int mode, char *saveme, int timeperiod)
{
  //cout << "mode = " << mode << endl;
  // Find timezone offset
  TDatime now,nowlocal,tstart,tend;
  now.Set(now.Convert(kTRUE));
  nowlocal.Set(nowlocal.Convert(kFALSE));
  int houroffset = now.GetHour()-nowlocal.GetHour();
  if (houroffset < -12 )   houroffset+=24;
  cout << "using houroffset " << houroffset << endl ;

  // Some lazy idiot got tired of scope problems amd moved the axes here.  It all needs a rewrite anyway.
  TGaxis *axis,*zoomaxis;

  //cout << " UTC offset is "  << houroffset << endl ;
  // The following constructor tries to open a file named "g", that's a problem.  LMM
  //TGraph *g = new TGraph("g");
  // Try with empty arguments or wait until we have something to plot. LMM
  TGraph *g = new TGraph();
  //char canvasname[20];
  //sprintf(canvasname,"c0_%s",branchname) ; 
  //if (gROOT->FindObject("c0")) 
  //    TCanvas *c0 = (TCanvas*)gROOT->FindObject("c0");
  // else 
  //    TCanvas *c0 = new TCanvas("c0");
  TCanvas *c0 = new TCanvas("c0","c0",1200,800);
  //TCanvas *c0 = new TCanvas(canvasname);
  cout << "Plotting with_zoom branchname = " << branchname << endl;
  char drawme[50];
  sprintf(drawme,"%s:Time_%s.Convert(kFALSE)-2*60*60*%d", branchname, branchname,houroffset);
  //for wierd Time based root axis reasons the offset is subracted twice, because the axis is also converted.
  //cout << "t->Draw(" << drawme << ");" << endl;
  t->Draw(drawme);
  if ( t->GetBranch(branchname) == 0 ) 
        {
        cout << "Branch " << branchname << " does not exist " << endl;
        return;
        }
  int nentries = t->GetBranch(branchname)->GetEntries();
  //cout << "Plotting " << nentries << " Entries " << endl;
  double *v = t->GetV1();
  //cout << "min,max = " << zoom[0] << ", " << zoom[1] << endl;
  if(*(std::min_element(v, v+nentries)) < zoom[0]) zoom[0] = *(std::min_element(v, v+nentries));
  if(*(std::max_element(v, v+nentries)) > zoom[1]) zoom[1] = *(std::max_element(v, v+nentries));
  //cout << "min,max = " << zoom[0] << ", " << zoom[1] << endl;
  double *times = t->GetV2();
  g = new TGraph(nentries,times,v);
  c1->cd();
  if(color==10)color=95; //10 is white. Useless. Make it orange.
  if(color==12)color=50; //12 is dark grey, too similar to black. Make it an ugly pink.
  g->SetMarkerColor(color);
  g->SetLineColor(color);
  g->SetMarkerStyle(marker);

  if(mode==0 || mode==3) //If it's the first channel plotted on this canvas (mode 0) or the only channel (mode 3), we need axes and titles
    {
      //cout << "Drawing axes" << endl;
      g->SetMinimum(ylimit[0]);
      g->SetMaximum(ylimit[1]);
      g->Draw("ALP");
      g->SetTitle(title);
      g->GetXaxis()->SetTimeDisplay(1);
      g->GetXaxis()->SetTimeOffset(0,"gmt");  // make X axis think it is GMT
      g->GetXaxis()->SetLabelSize(0.03);      
      // Adding information of hour and minute helped sort some confusing root date plotting issues.
      //if(timeperiod == 0) g->GetXaxis()->SetTimeFormat("%m/%d %H:%M");//for day plots, use hour:minute
      //else g->GetXaxis()->SetTimeFormat("%m/%d %H");//for week and month plots, use month/day
      if(timeperiod == 0) g->GetXaxis()->SetTimeFormat("%H:%M");//for day plots, use hour:minute
      else g->GetXaxis()->SetTimeFormat("%m/%d");//for week and month plots, use month/day
      g->GetXaxis()->SetLimits(timeago-houroffset*60*60,XNow-houroffset*60*60); //Convert limits to a UTC time
      g->GetYaxis()->SetTitleOffset(1.5);
      //tstart.Set(timeago,kFALSE);
      //tend.Set(XNow,kFALSE);
      //cout << " creating axes from " << tstart.AsString() << endl;
      //cout << " to " << tend.AsString() << endl ;
      //cout << " From timeago=" << timeago << " to XNow=" << XNow << endl;

    }
  else//Not first
    {
      g->Draw("same LP");
    }
  //There should be 1 entry per hour (i.e. 24 entries) but sometimes a reading is minutes (or even seconds) too early/too late so there are
  //only 23 entries. Anyway, if there are fewer than that, we need a warning due to missing entries
  //LMM Only worry about empty plots, goodNentries >0.
  int goodNEntries = 0;
  // Standard warning trigger
  if(timeperiod == 0)      goodNEntries = 23;
  else if(timeperiod == 1) goodNEntries = 23*7;
  else if(timeperiod == 2) goodNEntries = 23*7*30;
  // LMM 20180813 warning trigger at least 2/day
  if(timeperiod == 0)      goodNEntries = 1;
  else if(timeperiod == 1) goodNEntries = 1;
  else if(timeperiod == 2) goodNEntries = 1;
  if(nentries < goodNEntries && mode != 3) 
    {
      WarningText->SetLineColor(0);
      WarningText->SetFillColor(0);
      WarningText->SetBorderSize(0);
      WarningText->SetMargin(0.0);
      WarningText->SetTextAlign(11);
      WarningText->SetTextColor(color);
      char buff1[256];
      char buff2[10];
      if(nwarnings==0){
        sprintf(buff1, "Warning: check for missing entries in:");// %s", branchname);
        WarningText->AddText(buff1);
      }
      nwarnings+=1;
      sprintf(buff2, "%s", branchname);
      WarningText->AddText(buff2);
      //WarningText->Draw();
    }
  
  if(mode==2 || mode==3)//If it's the last or only plot - save and clear
    {
      char savename[100];
      UpdateText->Draw();
      if(nwarnings>0) WarningText->Draw();
      c1->Modified();
      c1->Update();
      if(timeperiod == 0) sprintf(savename,"%sDay.png", saveme);
      if(timeperiod == 1) sprintf(savename,"%sWeek.png", saveme);
      if(timeperiod == 2) sprintf(savename,"%sMonth.png", saveme);

      //TGraph *g1 = new TGraph("g1");
      TGraph *g1 = new TGraph();
      g1 = (TGraph*)gPad->GetPrimitive("Graph");
      gPad->SetGridx(1);
      char* newtitle = strrchr(title,';');
      //If it's a temperature plot, we need an extra axis so that it can be plotted in terms of Farenheit and Celsius
      if(title[strlen(title)-2] == 'C' || title[strlen(title)-2] == 'F')
      //if(temperature)
        {
          //cout << "Temperature plot" << endl;
          //char* newtitle = strrchr(title,';');
          newtitle = newtitle+1;
          //cout << newtitle << ", length = " << strlen(newtitle) << "--->";
          //TGaxis *axis;
          if(title[strlen(title)-2] == 'F')//Farenheit
            {
              axis = new TGaxis(gPad->GetUxmax(),gPad->GetUymin(),gPad->GetUxmax(),gPad->GetUymax(),(gPad->GetUymin()-32)/1.8,(gPad->GetUymax()-32)/1.8,510,"+L");
              //cout << "making 2nd axis with " << gPad->GetUxmax() << " " << gPad->GetUymin() << " " << gPad->GetUxmax() << " " 
              //<< gPad->GetUymax() << " " <<  (gPad->GetUymin()-32)/1.8 << " " << (gPad->GetUymax()-32)/1.8 << endl ; 
              newtitle[strlen(newtitle)-2] = 'C';
            }
          else if(title[strlen(title)-2] == 'C')//Celcius
            {
              axis = new TGaxis(gPad->GetUxmax(),gPad->GetUymin(),gPad->GetUxmax(),gPad->GetUymax(),(gPad->GetUymin()*1.8)+32,(gPad->GetUymax()*1.8)+32,510,"+L");
              //cout << "making 2nd axis with " << gPad->GetUxmax() << " " << 
              //gPad->GetUymin() << " " << gPad->GetUxmax() << " " << gPad->GetUymax() << " " 
              //<< (gPad->GetUymin()*1.8)+32 << " " << (gPad->GetUymax()*1.8)+32 << endl;
              newtitle[strlen(newtitle)-2] = 'F';
            }
          cout << newtitle << endl;
          axis->SetTitle(newtitle);
          axis->SetTitleSize(g->GetYaxis()->GetTitleSize());
          axis->SetLabelSize(g->GetYaxis()->GetLabelSize());
          axis->SetTitleFont(g->GetYaxis()->GetTitleFont());
          axis->SetLabelFont(g->GetYaxis()->GetTitleFont());
          axis->SetTitleOffset(1.5);
          axis->Draw();
        }

      c1->SaveAs(savename);
      //cout << "nonzoom file: " << saveme << endl;
      
      //For zoomed plot
      g1->SetMinimum(zoom[0] - (zoom[1]-zoom[0])*0.10);
      g1->SetMaximum(zoom[1] + (zoom[1]-zoom[0])*0.10);

      //cout << "Min: " << g1->GetMinimum() << " Max: " << g1->GetMaximum() << endl;
      //cout << "C Min: " << (g1->GetMinimum()-32)/1.8 << " C Max: " << (g1->GetMaximum()-32)/1.8 << endl;
      //cout << "title " << title << endl;
      
      if(title[strlen(title)-2] == 'C' || title[strlen(title)-2] == 'F')
        {
          //cout << "Temperature plot" << endl;
          axis->Delete();  // removes TGaxis from the unzoomed plot Why not just draw with new parameters? LMM
          if(newtitle[strlen(newtitle)-2] == 'C')
            zoomaxis = new TGaxis( gPad->GetUxmax(),g1->GetMinimum(),gPad->GetUxmax(),g1->GetMaximum(),(g1->GetMinimum()-32)/1.8,(g1->GetMaximum()-32)/1.8,510,"+L");
          else if(newtitle[strlen(newtitle)-2] == 'F')
            zoomaxis = new TGaxis( gPad->GetUxmax(),g1->GetMinimum(),gPad->GetUxmax(), g1->GetMaximum(),(g1->GetMinimum()*1.8)+32,(g1->GetMaximum()*1.8)+32,510,"+L");
          zoomaxis->SetTitle(newtitle);
          zoomaxis->SetTitleSize(g->GetYaxis()->GetTitleSize());
          zoomaxis->SetLabelSize(g->GetYaxis()->GetLabelSize());
          zoomaxis->SetTitleFont(g->GetYaxis()->GetTitleFont());
          zoomaxis->SetLabelFont(g->GetYaxis()->GetLabelFont());
          zoomaxis->Draw();
          zoomaxis->SetTitleOffset(1.5);
        }

      //cout << "modifying and updating canvas" << endl;
      c1->Modified();
      c1->Update();
      
      //cout << saveme << ":" << strlen(saveme);
      //  int pos = saveme.length();
      //cout <<"changing name..." << endl;
      if(timeperiod == 0)      sprintf(savename,"%szoomDay.png", saveme);
      else if(timeperiod == 1) sprintf(savename,"%szoomWeek.png", saveme);
      else if(timeperiod == 2) sprintf(savename,"%szoomMonth.png", saveme);
      //cout << "saving zoom..." << endl;
      c1->SaveAs(savename);
      //cout << "zoom file: " << saveme << endl;  
      c1->Clear();
      zoom[0] = 99999;
      zoom[1] =-99999;
      nwarnings=0;
      WarningText->Clear();
      //cout << "all clear, start again" << endl;
      //if(title[strlen(title)-2] == 'C' || title[strlen(title)-2] == 'F') zoomaxis->Delete();
    }
}

//Like plot_with_zoom, but calculates the strain from the 3 strain gauges
void plot_strain_with_zoom(TTree *t, char *branchname, int x, int y, int marker, TCanvas *c1, char *title, float *ylimit, float *zoom, int timeago, int XNow, int mode, char *saveme, int timeperiod)
{
  //TGraph *g = new TGraph("g");
  TDatime now,nowlocal;
  now.Set(now.Convert(kTRUE));
  nowlocal.Set(nowlocal.Convert(kFALSE));
  int houroffset = now.GetHour()-nowlocal.GetHour();
  if (houroffset < -12 )   houroffset+=24;
  //cout << " UTC offset is "  << houroffset << endl ;

  //char canvasname[20];
  //sprintf(canvasname,"c0_%s",branchname) ; 
  //if (gROOT->FindObject("c0"))
  //    TCanvas *c0 = (TCanvas*)gROOT->FindObject("c0");
  //else 
  //    TCanvas *c0 = new TCanvas("c0");
  TCanvas *c0 = new TCanvas("c0");
  //TCanvas *c0 = new TCanvas(canvasname);
  cout << "Plotting strain branchname = " << branchname << endl;
  char drawme[50];

  sprintf(drawme,"%s1:Time_%s1.Convert(kFALSE)-60*60*%d", branchname, branchname,houroffset);
  //cout << "t->Draw(" << drawme << ");" << endl;
  t->Draw(drawme);
  char realbranchname[20];
  sprintf(realbranchname,"%s1",branchname);
  int nentries = t->GetBranch(realbranchname)->GetEntries();
  double *v1 = t->GetV1();
  double *times = t->GetV2();

  sprintf(drawme,"%s2:Time_%s2.Convert(kFALSE)-60*60*%d", branchname, branchname,houroffset);
  //cout << "t->Draw(" << drawme << ");" << endl;
  t->Draw(drawme);
  double *v2 = t->GetV1();

  sprintf(drawme,"%s3:Time_%s3.Convert(kFALSE)-60*60*%d", branchname, branchname,houroffset);
  //cout << "t->Draw(" << drawme << ");" << endl;
  t->Draw(drawme);
  double *v3 = t->GetV1();

  for(int entry=0;entry<nentries;entry++)
    {
      v1[entry] = 0.5*(v1[entry]+v3[entry]) + (1/sqrt(2))*sqrt(pow(v1[entry]-v2[entry],2)+(pow(v2[entry]-v3[entry],2)));
      //cout << "Strain in " << branchname << " = " << v1[entry];
      if(v1[entry] < zoom[0]) zoom[0] = v1[entry];
      if(v1[entry] > zoom[1]) zoom[1] = v1[entry];
      //cout << "min, max: " << zoom[0] << ", " << zoom[1] << endl;
    }

  TGraph *g = new TGraph(nentries, times, v1);

  c1->cd();
  if(branchname[2] == 'N')
    {
      if(x==1)g->SetMarkerColor(kGreen+y-1);
      if(x==2)g->SetMarkerColor(kRed+y-1);
      if(x==3)g->SetMarkerColor(kGray+y-1);
      if(x==4)g->SetMarkerColor(kCyan+y-1);
      if(x==5)g->SetMarkerColor(kBlue+y-1);
      if(x==6)g->SetMarkerColor(kMagenta+y-1);
    }
  else
    {
      if(y==1)g->SetMarkerColor(kRed+3*(x-1));
      if(y==2)g->SetMarkerColor(kGreen+3*(x-1));
      if(y==3)g->SetMarkerColor(kBlue+3*(x-1));
    }
  g->SetLineColor(g->GetMarkerColor());  
  g->SetMarkerStyle(marker);
  //cout << "x,y,marker, marker size: " << x << "," << y << "," << marker << "," << g->GetMarkerSize() << endl;
  if(mode==0)//x==1 && y==1 && marker == 20 && g->GetMarkerSize()==1)
    {
      //cout << "Drawing axes" << endl;
      g->SetMinimum(ylimit[0]);
      g->SetMaximum(ylimit[1]);
      g->Draw("ALP");
      g->SetTitle(title);
      g->GetXaxis()->SetTimeDisplay(1);
      g->GetXaxis()->SetLabelSize(0.03);
      if(timeperiod == 0) g->GetXaxis()->SetTimeFormat("%H:%M");
      else                g->GetXaxis()->SetTimeFormat("%m/%d");
      g->GetXaxis()->SetLimits(timeago,XNow);
      g->GetYaxis()->SetTitleOffset(1.5);
    }
  else
    {
      g->Draw("same LP");
    }

  if(mode==2)
    {
      char savename[100];
      UpdateText->Draw();
      c1->Modified();
      c1->Update();
      if(timeperiod == 0)      sprintf(savename,"%sDay.png", saveme);
      else if(timeperiod == 1) sprintf(savename,"%sWeek.png", saveme);
      else if(timeperiod == 2) sprintf(savename,"%sMonth.png", saveme);


      //TGraph *g1 = new TGraph("g1");
      TGraph *g1 = new TGraph();
      g1 = (TGraph*)gPad->GetPrimitive("Graph");

      c1->SaveAs(savename);
      //cout << "nonzoom file: " << saveme << endl;
      
      //For zoomed plot
      g1->SetMinimum(zoom[0] - (zoom[1]-zoom[0])*0.1);
      g1->SetMaximum(zoom[1] + (zoom[1]-zoom[0])*0.1);

      //cout << "Min: " << g1->GetMinimum() << " Max: " << g1->GetMaximum() << endl;
      c1->Modified();
      c1->Update();
      
      //cout << saveme << ":" << strlen(saveme);
      //  int pos = saveme.length();
      if(timeperiod == 0)      sprintf(savename,"%szoomDay.png", saveme);
      else if(timeperiod == 1) sprintf(savename,"%szoomWeek.png", saveme);
      else if(timeperiod == 2) sprintf(savename,"%szoomMonth.png", saveme);
      c1->SaveAs(savename);
      //cout << "zoom file: " << saveme << endl;  
      c1->Clear();
      zoom[0] = 99999;
      zoom[1] =-99999;

    }

}

//Main function
// void DCSPlotMaker(int timeperiod, int detector)
void DCSPlotMaker(string timeperiodname, string detectorname)
{
  int timeperiod = 0;
  int detector = 0;

  if(timeperiodname.compare("day") == 0) {
    timeperiod = 0;
    std::cout << "Making day plots\n";
  } else if(timeperiodname.compare("week") == 0) {
    timeperiod = 1;
    std::cout << "Making week plots\n";
  } else if(timeperiodname.compare("month") == 0) {
      timeperiod = 2;
      std::cout << "Making month plots\n";
  } else {
    std::cout << "Macro should take input of 'day', 'week', or 'month'. \n"
              << "You have input '" << timeperiodname
              << "'\nSomething's wrong. Abort.";
    return;
  }

  if(detectorname.compare("ND") == 0) {
    detector = 1;
    std::cout << "Making plots for ND\n";
  } else if(detectorname.compare("FD") == 0) {
    detector = 2;
    std::cout << "Making plots for FD\n";
  } else {
    std::cout << "Macro should take input of 1 for ND, 2 for FD"
              << "You have input " << detector
              << "... Something's wrong. Abort.";
    return;
  }

  gStyle->SetPalette(55);
  gStyle->SetMarkerStyle(7);
  gStyle->SetFillStyle(0);

  //TFile *f = new TFile("/nova/app/users/mcampbel/dcs_db/automated/temp_dcs_db_readout_multichan_day.root");
  TFile *f;
  if(detector == 1){
    //    if(timeperiod == 0)      f = new TFile("dcs_db_readout_multichan_nd_day.root",   "read");
    if(timeperiod == 0)      f = new TFile("trees/dcs_db_readout_multichan_day_ND.root",   "read");
    else if(timeperiod == 1) f = new TFile("trees/dcs_db_readout_multichan_week_ND.root",  "read");
    else if(timeperiod == 2) f = new TFile("trees/dcs_db_readout_multichan_month_ND.root", "read");
  } else if(detector == 2) {
    if(timeperiod == 0)      f = new TFile("trees/dcs_db_readout_multichan_day_FD.root",   "read");
    else if(timeperiod == 1) f = new TFile("trees/dcs_db_readout_multichan_week_FD.root",  "read");
    else if(timeperiod == 2) f = new TFile("trees/dcs_db_readout_multichan_month_FD.root", "read");
  }

  TTree *t = (TTree*)f->Get("dcs_db_tree");
  TCanvas *c1 = new TCanvas("c1","c1",1200,800);
  TCanvas *c2 = new TCanvas("c2","c2",1200,800);

  //TLegend *leg = new TLegend(0.9,0.1,1.0,0.9);

  char name[10];
  char saveme[50];
  char saveme2[50];
  char title[100];
  char title2[100];
  char branchname[10];
  char* filepath;
  //If ND, call the plots NearDet-something
  if(detector == 1) {
    if(timeperiod == 0)      filepath = "dayplots/NearDet-";
    else if(timeperiod == 1) filepath = "weekplots/NearDet-";
    else if(timeperiod == 2) filepath = "monthplots/NearDet-";
  } else if(detector == 2){
    //If FD, call them FarDet-something
    if(timeperiod == 0)      filepath = "dayplots/FarDet-";
    else if(timeperiod == 1) filepath = "weekplots/FarDet-";
    else if(timeperiod == 2) filepath = "monthplots/FarDet-";
  }
  //I don't know why, but there's a segmentation violation after the first plot any time I try to define a full filepath
  //To get around this, the shell script directs the user to the correct directory and works from there
  //then copies the contents of dayplots to /nusoft/app/web/htdoc/nova/datacheck/nearline/plots/
  //Sorry for the hackiness...

  float ylimit[2];//For setting axis limits
  float zoom[2] = {99999, -99999};//For setting axis limits for zoomed plots
  float zoom2[2] = {99999, -99999};
  int mode = 0;
  int mode2 = 0;
  t->UseCurrentStyle();

  //FOR X AXIS
  //TDatime * TNowGMT = new TDatime;
  //TNowGMT->Set();
  //cout << TNowGMT->AsString() << endl;
  TDatime *TNow = new TDatime;
  TDatime tstart,tend;
  TNow->Set();  // sets time o now
  //cout << " Running at " << TNow->AsString() << endl;
  TNow->Set(TNow->GetYear(),TNow->GetMonth(),TNow->GetDay(),TNow->GetHour()+1,0,0);  //LMM added one hour
  //cout << " Dropped Min/Sec? " << TNow->AsString() << endl;  cout << " Running at " << TNow->AsString() << endl;
  // unused ? TDatime *SRtime = new TDatime;
  //int XNowGMT = TNowGMT->Convert(kTRUE);
  UInt_t XNow = TNow->Convert(kFALSE);  // Gets current time
  //cout << XNowGMT << "   " << XNow << endl;
  //int GMToffset = XNowGMT - XNow;
  //XNow = XNow - 30*60;//To make up for the time the python code takes to read the database

  //gStyle->SetTimeOffset(GMToffset);

  UInt_t timeago = 0;
  if(timeperiod == 0)      timeago = XNow - 60*60*24;
  else if(timeperiod == 1) timeago = XNow - 60*60*24*7;
  else if(timeperiod == 2) timeago = XNow - 60*60*24*30;

  //tstart.Set(timeago,kFALSE);
  //cout << " Tstart " << tstart.AsString() << endl;
  //tend.Set(XNow,kFALSE);
  //cout << " Tend " << tend.AsString() << endl;
  //cout << " will create axes from " << tstart.AsString() << endl ;
  //cout << " to " << tend.AsString() << endl ;
  //cout << " From timeago=" << timeago << " to XNow=" << XNow << endl;

  //FOR TIMESTAMP
  TDatime *Tfinish = new TDatime;
  Tfinish->Set();
  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);

  //POWER SUPPLY ---------------------------------------------------------------------------------------------------
  int ndiblocks = 0; //Number of diblocks
  if(detector == 1) ndiblocks = 4;
  else if(detector == 2) ndiblocks = 14;

  char lmh[4] = "LMH"; //Measure low, medium, and high voltage
  char measstr[3];
  int ndevices = 0; //Number of devices per diblock
  if(detector == 1) ndevices = 4;
  else if(detector == 2) ndevices = 12;
  int nps = 3; //Number of power supplies (i.e. low, med, high)
  int nmeas = 4; //Number of measurements per device
  int col = 800; //Colour of graph (I like 800-910 and 920 onwards)
  int col2 = 800;
  int mark = 20; //Marker style
  
  for(int meas = 0;meas<nmeas;meas++)
    {
      if(meas == 0)strcpy(measstr,"");//Voltage at device
      if(meas == 1)strcpy(measstr,"MC");//Current
      if(meas == 2)strcpy(measstr,"MV");//Voltage at source
      if(meas == 3)strcpy(measstr,"TP");//Temperature
      
      for(int ps = 0;ps<nps;ps++)
        {
          for(int diblock = 1;diblock<ndiblocks+1;diblock++)
            {
              //cout << "diblock = " << diblock << ", detector = " << detector;
              if(detector == 1 && diblock == 4)
                {
                  sprintf(name,"%c0M", lmh[ps], diblock);//muon catcher
                  //cout << endl << "CHANGE NAME TO ACCOUNT FOR MUON CATCHER" << endl;
                  ndevices = 2;
                }
              else if(diblock<10) sprintf(name,"%c0%i", lmh[ps], diblock);
              else sprintf(name,"%c%i", lmh[ps], diblock);
              //cout << name << endl;
              
              //Setting titles and y axis range Probably need adjustment/rethink.  LMM
              if(meas == 0)
                {
                  if(ps ==0){
                    sprintf(title, "Low V power supply, diblock %i: Voltage at device;;Voltage (V)", diblock);
                    sprintf(title2, "Low V power supply: Voltage at device;;Voltage (V)", diblock);
                    ylimit[0] = 2.5;
                    ylimit[1] = 4.5;//These limits are set by me looking at the plots and guessing a useful range, they are not set in stone
                  }
                  if(ps ==1){
                    sprintf(title, "Medium V power supply, diblock %i: Voltage at device;;Voltage (V)", diblock);
                    sprintf(title2, "Medium V power supply: Voltage at device;;Voltage (V)", diblock);
                    ylimit[0] = 20;
                    ylimit[1] = 26;
                  }
                  if(ps ==2){
                    sprintf(title, "High V power supply, diblock %i: Voltage at device;;Voltage (V)", diblock);
                    sprintf(title2, "High V power supply: Voltage at device;;Voltage (V)", diblock);
                    ylimit[0] = 422;
                    ylimit[1] = 428;
                  }
                }

              if(meas == 1)
                {
                  if(ps ==0){
                    sprintf(title, "Low V power supply, diblock %i: Current;;Current (mA)", diblock);
                    sprintf(title2, "Low V power supply: Current;;Current (mA)", diblock);
                    ylimit[0] = 30000;
                    ylimit[1] = 31000;
                  }
                  if(ps ==1){
                    sprintf(title, "Medium V power supply, diblock %i: Current;;Current (mA)", diblock);
                    sprintf(title2, "Medium V power supply: Current;;Current (mA)", diblock);
                    ylimit[0] = 8000;
                    ylimit[1] = 9900;
                  }                  
                  if(ps ==2){
                    sprintf(title, "High V power supply, diblock %i: Current;;Current (A)", diblock);
                    sprintf(title2, "High V power supply: Current;;Current (A)", diblock);
                    ylimit[0] = 3.;
                    ylimit[1] = 4.5;
                  }
                }

              if(meas == 2)
                {
                  if(ps ==0){
                    sprintf(title, "Low V power supply, diblock %i: Voltage at source;;Voltage (V)", diblock);
                    sprintf(title2, "Low V power supply: Voltage at source;;Voltage (V)", diblock);
                    ylimit[0] = 3.5;
                    ylimit[1] = 5.;
                  }                  
                  if(ps ==1){
                    sprintf(title, "Medium V power supply, diblock %i: Voltage at source;;Voltage (V)", diblock);
                    sprintf(title2, "Medium V power supply: Voltage at source;;Voltage (V)", diblock);
                    ylimit[0] = 22.;
                    ylimit[1] = 28.;
                  }                  
                  if(ps ==2){
                    sprintf(title, "High V power supply, diblock %i: Voltage at source;;Voltage (V)", diblock);
                    sprintf(title2, "High V power supply: Voltage at source;;Voltage (V)", diblock);
                    ylimit[0] = 422.;
                    ylimit[1] = 428.;
                  }
                }

              if(meas == 3)
                {
                  if(ps ==0){
                    sprintf(title, "Low V power supply, diblock %i: Temperature;;Temperature (C)", diblock);
                    sprintf(title2, "Low V power supply: Temperature;;Temperature (C)", diblock);
                    ylimit[0] = 20.;
                    ylimit[1] = 35.;
                  }                  
                  if(ps ==1){
                    sprintf(title, "Medium V power supply, diblock %i: Temperature;;Temperature (C)", diblock);
                    sprintf(title2, "Medium V power supply: Temperature;;Temperature (C)", diblock);
                    ylimit[0] = 20;
                    ylimit[1] = 30;
                  }                  
                  if(ps ==2){
                    sprintf(title, "High V power supply, diblock %i: Temperature;;Temperature (C)", diblock);
                    sprintf(title2, "High V power supply: Temperature;;Temperature (C)", diblock);
                    ylimit[0] = 25;
                    ylimit[1] = 35;
                  }
                }

              sprintf(saveme,"%s%s%s", filepath, name, measstr);
              sprintf(saveme2,"%s%c%s", filepath, lmh[ps], measstr);
          
              for(int device=1;device<ndevices+1;device++)
                {
                  if(device==ndevices)
                    {
                      mode=2;//If it's the last device, save
                      if(diblock==ndiblocks)mode2=2;
                    }
                  if(device<10) sprintf(branchname,"%s0%i%s", name, device, measstr);
                  else sprintf(branchname, "%s%i%s", name, device, measstr);
                  // plot once for diblock
                  plot_with_zoom(t, branchname, device, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
                  // plot again for detector -- maybe combine to access the data once?
                  plot_with_zoom(t, branchname, device, 19+diblock, c2, title2, ylimit, zoom2, timeago, XNow, mode2, saveme2, timeperiod);
                  mode = 1;
                  mode2 = 1;
                }
              //saveplot(c1, leg, UpdateText, saveme, zoom, ylimit);
              mode = 0;
            }//end diblock<ndiblocks
          mode2 = 0;
          // Reset number of devices/diblock 
          if(detector == 1) ndevices = 4;
          else if(detector == 2) ndevices = 12;
        }//end ps<nps
    }//end meas<nmeas
  mode = 0;
  mode2 = 0;
   //*/

   //UTLOPC-------------------------------------------------------------------------------------------------------------------------------
   ///*

   //----------BUILDING FLOW----------
   /*
   ylimit[0] = 0.;
   ylimit[1] = 5.;
   sprintf(name,"EBF");
   //cout << name << endl;

   sprintf(title,"Building flow;;");
   sprintf(branchname,"%s1",name);
   sprintf(legendentry,"BF1");
   plot_with_zoom(t, branchname, 1, c1, title, g1, legendentry, ylimit, zoom, timeago, XNow);
   sprintf(branchname,"%s2",name);
   sprintf(legendentry,"BF2");
   plot_with_zoom(t, branchname, 2, c1, title, g2, legendentry, ylimit, zoom, timeago, XNow);
   sprintf(branchname,"%s3",name);
   sprintf(legendentry,"BF3");
   plot_with_zoom(t, branchname, 3, c1, title, g3, legendentry, ylimit, zoom, timeago, XNow);

   c1->cd();
   leg->Draw();
   c1->Modified();
   c1->Update();
   sprintf(saveme,"%sEBFDay.png", filepath);

   //For zoomed plot
   g1 = (TGraph*)gPad->GetPrimitive("Graph");
   g1->SetMinimum(zoom[0] - (zoom[1]-zoom[0])*0.1);
   g1->SetMaximum(zoom[1] + (zoom[1]-zoom[0])*0.1);
   c1->Modified();
   c1->Update();
   sprintf(saveme,"%sEBFzoomDay.png", filepath, name);
   c1->SaveAs(saveme);
          
   c1->Clear();
   leg->Clear();
   zoom[0] = 99999;
   zoom[1] =-99999;
 */

  //AIR FLOW-----------------------------------------------------                                                                                                                 
  if(detector == 1)
    {
      //cout << "starting airflow" << endl;
      sprintf(title, "Rack Air Flow");
      ylimit[0] = 0;
      ylimit[1] = 300;
      sprintf(saveme,"%sAirFlow", filepath);
      
      sprintf(branchname,"ENNALV1");
      //sprintf(legendentry,"LV1");
      plot_with_zoom(t, branchname, 1, 20, c1, title, ylimit, zoom, timeago, XNow, 0, saveme, timeperiod);
      sprintf(branchname,"ENNALV2");
      //sprintf(legendentry,"LV2");
      plot_with_zoom(t, branchname, 2, 20, c1, title, ylimit, zoom, timeago, XNow, 1, saveme, timeperiod);
      sprintf(branchname,"ENNAHV1");
      //sprintf(legendentry,"HV");
      plot_with_zoom(t, branchname, 3, 20, c1, title, ylimit, zoom, timeago, XNow, 2, saveme, timeperiod);
    }

  //----------DEWPOINT & TEMPERATURE----------

  if(detector == 1){
    nmeas = 2;
    char dh[3] = "DH";

    zoom[0]=99999.;
    zoom[1]=-99999.;

    for(int iter=0;iter<nmeas;iter++)
      {
        if(iter==0)
          {
            sprintf(title, "Dewpoint;;Dewpoint (C)");
            ylimit[0] = -60;
            ylimit[1] = 30;
            sprintf(saveme,"%sDewPoint", filepath);
          }
        if(iter==1)
          {
            sprintf(title, "Humidity;;Humidity (%)");
            ylimit[0] = 0;
            ylimit[1] = 90;
            sprintf(saveme,"%sHumidity", filepath);
          }

        sprintf(branchname,"ENN%cDB1", dh[iter]);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;
        plot_with_zoom(t, branchname, 1, 20, c1, title, ylimit, zoom, timeago, XNow, 0, saveme, timeperiod);
        sprintf(branchname,"ENN%cDB2", dh[iter]);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;
        plot_with_zoom(t, branchname, 2, 20, c1, title, ylimit, zoom, timeago, XNow, 1, saveme, timeperiod);
        sprintf(branchname,"ENN%cDB3", dh[iter]);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;
        plot_with_zoom(t, branchname, 3, 20, c1, title, ylimit, zoom, timeago, XNow, 1, saveme, timeperiod);
        sprintf(branchname,"ENN%cDME", dh[iter]);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;
        plot_with_zoom(t, branchname, 4, 20, c1, title, ylimit, zoom, timeago, XNow, 1, saveme, timeperiod);
        sprintf(branchname,"ENN%cDMN", dh[iter]);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;
        plot_with_zoom(t, branchname, 5, 20, c1, title, ylimit, zoom, timeago, XNow, 2, saveme, timeperiod);
        //cout << " zoom params " << zoom[0] << " and " << zoom[1] << endl;

      }
  }//end if (ND)

  else if(detector == 2){
    sprintf(measstr,"HT");
    nmeas = 2;
    
    for(int iter = 0;iter<nmeas;iter++)
      {
        if(iter == 0)
          {
            ylimit[0] = 0.;
            ylimit[1] = 70.;
            sprintf(title,"Dewpoint;;Dewpoint (F)");
          }
        if(iter == 1)
          {
            ylimit[0] = 60;
            ylimit[1] = 90;
            sprintf(title,"Temperature;;Temperature (F)");
          }
        sprintf(name,"E%cL", measstr[iter]);
        //cout << name << endl;
        if(iter==0)sprintf(saveme,"%sDewPoint", filepath);
        else sprintf(saveme,"%sTemperature", filepath);
      
        sprintf(branchname,"%s101", name);//Weirdly ordered colors to match a previous code, where this order made sense!
        plot_with_zoom(t, branchname, 1, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        mode = 1;
        sprintf(branchname,"%s102", name);
        plot_with_zoom(t, branchname, 2, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        sprintf(branchname,"%s103", name);
        plot_with_zoom(t, branchname, 3, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
       
        if(iter == 1)
          {        
            sprintf(branchname,"%sL0ME", name);
            plot_with_zoom(t, branchname, 95, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
            sprintf(branchname,"%sL0NE", name);
            plot_with_zoom(t, branchname, 11, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
            sprintf(branchname,"%sL0SE", name);
            plot_with_zoom(t, branchname, 14, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
          }
      
        sprintf(branchname,"%sL0MW", name);
        plot_with_zoom(t, branchname, 4, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        sprintf(branchname,"%sL0NW", name);
        plot_with_zoom(t, branchname, 5, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        sprintf(branchname,"%sL0SW", name);
        plot_with_zoom(t, branchname, 6, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        /*
          sprintf(branchname,"%sL1MW", name);
          sprintf(legendentry,"L1MW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL1NW", name);
          sprintf(legendentry,"L1 NW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL1SW", name);
          sprintf(legendentry,"L1SW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL2MW", name);
          sprintf(legendentry,"L2 MW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL2NW", name);
          sprintf(legendentry,"L2 NW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL2SW", name);
          sprintf(legendentry,"L2 SW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL3MW", name);
          sprintf(legendentry,"L3MW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL3NW", name);
          sprintf(legendentry,"L3 NW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
          sprintf(branchname,"%sL3SW", name);
          sprintf(legendentry,"L3 SW");
          plot_with_zoom(t, branchname, 1, 20, c1, title, leg, legendentry, ylimit, zoom, timeago, XNow);
        */
        sprintf(branchname,"%sL4MW", name);
        plot_with_zoom(t, branchname, 7, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        sprintf(branchname,"%sL4NW", name);
        plot_with_zoom(t, branchname, 8, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        sprintf(branchname,"%sL4SW", name);
        mode = 2;
        plot_with_zoom(t, branchname, 9, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        mode = 0;
      }//end iter over h vs t
  }//end if (FD)

  //----------------STATUSES-------------------                                                                                                                                   
  if(detector == 1){
    int ndastatuses = 7;//7 dry air statuses                                                                                                                                        
    sprintf(title,"Status of systems");
    ylimit[0] = -1;
    ylimit[1] = 11;
    sprintf(saveme,"%sStatuses", filepath);
    mode = 0;
    
    for(int dastatus=1;dastatus<ndastatuses+1;dastatus++)
      {
        //sprintf(legendentry,"Dry air status %i",dastatus);                                                                                                                        
        sprintf(branchname,"ENNDAS%i",dastatus);
        if(dastatus==ndastatuses)mode = 2;
        plot_with_zoom(t, branchname, dastatus, 20, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
        mode = 1;
      }
    /*  sprintf(legendentry,"EPICS Reset status");
        sprintf(branchname,"ENNEPIC");                                                                                                                                                  plot_with_zoom(t, branchname, ndastatuses+1, c1, title, ylimit, zoom, timeago, XNow);*/
  }

  //PHIDGT VARIABLES--------------------------------------------------------------------------------------------------------------

  if(detector == 2){  
    char side[4] = "WNS";
    int nblocks = 29;//28 plus south bookend
    int nsides = 2;//North and south, or north and west
    int xmax;
    int ymax;
    col=800;
    int marker=20;
  
    for(int iter=0;iter<nsides;iter++)
      {

        //----------STRAIN----------------
        nblocks = 8;      
        for(int block=0;block < nblocks;block++)
          {
            col = 800;   
            ylimit[0] = -500;
            ylimit[1] = 1000;
          
            char realblock[3];
            if(block == 0)sprintf(realblock, "00");
            if(block == 1)sprintf(realblock, "01");
            if(block == 2)sprintf(realblock, "08");
            if(block == 3)sprintf(realblock, "09");
            if(block == 4)sprintf(realblock, "18");
            if(block == 5)sprintf(realblock, "19");
            if(block == 6)sprintf(realblock, "26");
            if(block == 7)sprintf(realblock, "27");
          
            if(iter==0)sprintf(title,"Strain measured on West side of blocks;;Micro strain");
            if(iter==1)sprintf(title,"Strain measured on North side of blocks;;Micro strain");
          
            //sprintf(name,"ES%i%c", block, side[iter]);
            sprintf(name,"ES");
            //cout << name << endl;
          
            if(iter==1)
              {
                xmax = 6;
                ymax = 4;
              }
            else if(iter==0)
              {
                xmax = 2;
                ymax = 3;
              }
            //cout << "xmax = " << xmax << endl;
          
            for(int x=1;x<xmax+1;x++)
              {
                //cout << "x = " << x << endl;
                for(int y=1;y<ymax+1;y++)
                  {
                    if(x==xmax && y==ymax) mode = 2;
                    //cout << "x,y,block: " << x << "," << y << "," << block << endl;
                    //print_strings("strain", block, x, y, name, branchname, legendentry);
                    sprintf(branchname,"%s%i%c%i%i",name,block,side[iter],x,y);
                    sprintf(saveme,"%s%s%c",filepath, name,side[iter]);
                    //cout << "col = " << col << endl;
                    plot_strain_with_zoom(t, branchname, x, y, marker, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
                    col += 1;
                    mode = 1;
                    if(col==910)col=920;//910-920 is white 
                  }
              }
            marker += 1;
            if(marker > 34)
              {
                marker = 20;
                gStyle->SetMarkerSize(2);
              }
          }//end iteration over blocks
      
        //sprintf(saveme,"%sFarDet-Strain%c",filepath,side[iter]);
        //saveplot(c1, leg, UpdateText, saveme, zoom, ylimit);
      
        gStyle->SetMarkerSize(1);
        marker = 20;
        mode = 0;

        //*//*
        //----------DISTANCE----------
        nblocks = 29;
        for(int block=0;block < nblocks;block++)
          {
            sprintf(saveme,"%sDist%c", filepath, side[iter+1]);

            col2 = 1;
            ylimit[0] = -0.2;
            ylimit[1] = 1.;
            if((block==0 && iter==1) || (block==27 && iter==0) || (block==28 && iter==1)) continue;

            if(iter==0) sprintf(title,"Distance from original position (North side of blocks);;Distance (inches)");
            if(iter==1) sprintf(title,"Distance from original position (South side of blocks);;Distance (inches)");
       
            if(block < 10) sprintf(name,"ED0%i",block);
            else if (block < 28) sprintf(name,"ED%i",block);
            else sprintf(name,"EDSB");
          
            //cout << name << endl;

            xmax = 2;

            for(int x=1; x<xmax+1; x++)
              {
                if(iter == 1 && x == 2)
                  {
                    x=6;
                    if(block==nblocks-2)mode = 2;
                  }
                //int nwchar = iter +1;
                sprintf(branchname,"%s%c%i1", name, side[iter+1], x);
                plot_with_zoom(t, branchname, col2, marker, c2, title, ylimit, zoom2, timeago, XNow, mode, saveme, timeperiod);
                mode = 1;
                col2 += 1;
                if(col2==910)col2=920;//910-920 is white

                if(iter==0)
                  {
                    if(block!=28)sprintf(branchname,"%s%c%i9",name, side[iter+1],x);
                    else
                      {
                        sprintf(branchname,"%s%c%i6",name, side[iter+1],x);
                        if(x==xmax && block==nblocks-1)mode = 2;
                      }
                    plot_with_zoom(t, branchname, col2, marker, c2, title, ylimit, zoom2, timeago, XNow, mode, saveme, timeperiod);
                    col2 += 1;
                    if(col2==910)col2=920;//910-920 is white
                  }
              }
            marker += 1;
            if(marker > 34)
              {
                marker = 20;
                gStyle->SetMarkerSize(2);
              }
          
          } //end iteration over blocks

        marker = 20;
        gStyle->SetMarkerSize(1);

        mode = 0;
      }  //end iteration over sides
    gStyle->SetMarkerSize(1);  

    //----------TEMPERATURE----------- 
   
    nsides = 2;
    nblocks = 8;
    ylimit[0] = 10;
    ylimit[1] = 25;
    xmax = 2;
    sprintf(title,"Internal temperature;;Temperature (C)");
    sprintf(saveme,"%sIntTemp", filepath);
  
    for(int s=0;s<nsides;s++)
      {
        if(s==0) sprintf(name,"ETN");
        else if(s==1) sprintf(name, "ETW");

        for(int block=0;block<nblocks;block++)
          {
            for(int x=1;x<xmax+1;x++)
              {
                if(s==nsides-1 && block==nblocks-1 && x==xmax) mode = 2;
                col = x + xmax*s;
                sprintf(branchname,"ET%i%c%i1",block, side[s],x);
                plot_with_zoom(t, branchname, col, 20+block, c1, title, ylimit, zoom, timeago, XNow, mode, saveme, timeperiod);
                mode = 1;
              }
          }
      }
    //*/

    //---------OAT, GLYCOT, HTWARN-------

    sprintf(title,"Outside air temperature;;Temperature (F)");
    sprintf(saveme, "%sOAT", filepath);
    sprintf(branchname,"EOAT");
    ylimit[0] = -40.;
    ylimit[1] = 90.;
    plot_with_zoom(t, branchname, 1, 7, c1, title, ylimit, zoom, timeago, XNow, 3, saveme, timeperiod);

    sprintf(title,"Cooling liquid temperature;;Temperature (F)");
    sprintf(saveme, "%sGLYCOT", filepath);
    sprintf(branchname,"EGLYCOT");
    ylimit[0] = 20.;
    ylimit[1] = 120.;
    plot_with_zoom(t, branchname, 1, 7, c1, title, ylimit, zoom, timeago, XNow, 3, saveme, timeperiod);

    sprintf(title,"High Temperature warning status");
    sprintf(saveme, "%sHTWarn", filepath);
    sprintf(branchname,"EHTWARN");
    ylimit[0] = -1;
    ylimit[1] = 2;
    plot_with_zoom(t, branchname, 1, 7, c1, title, ylimit, zoom, timeago, XNow, 3, saveme, timeperiod);

  }//End FD only plots
}//End EVERYTHING. We're done here.