CheckTrueDate.C

Go to the documentation of this file.

#include <fstream>
#include <TH1F.h>
#include <TTree.h>
#include <TFile.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TDatime.h>

void CheckTrueDate(int idet, int timeOpt) {

  std::cout << "\n\nChecking Actual Dates for new files...\n\n";

  TDatime *TNowGMT = new TDatime;  // current GMT time
  TDatime *TNow    = new TDatime;  // current local time
  TDatime *SRtime  = 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 X4DyAgo  = XNow - 60*60*24*4;
  int XWekAgo  = XNow - 60*60*24*7;
  int X2HrAgo  = XNow - 60*60*2;
  int X8HrAgo  = XNow - 60*60*8;
  int Xsrtime   = 0;

  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);

  const unsigned int Npoint = 1440;
  double time = 0.0;
  
  //
  // Check date in files
  //
  std::ifstream inDay;
  std::ofstream out_file;
  if(idet==1){
    inDay.open("NewFarDetFiles.txt");
    out_file.open ("NewFarDetOnMonFiles.txt");
  }
  else if(idet==2){
    inDay.open("NewNearDetFiles.txt");
    out_file.open ("NewNearDetOnMonFiles.txt");
  }
  else cout << " invalid detector type in CheckTrueDate " << endl;
  while(inDay.good()) {

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

    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->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;

      if(timeOpt == 0 && Xsrtime > XDayAgo)  out_file << filenameD << std::endl;
      if(timeOpt == 1 && Xsrtime > X2HrAgo)  out_file << filenameD << std::endl;
      if(timeOpt == 2 && Xsrtime > X8HrAgo)  out_file << filenameD << std::endl;
      if(timeOpt == 3 && Xsrtime > XWekAgo)  out_file << filenameD << std::endl;
      if(timeOpt == 4 && Xsrtime > X4DyAgo)  out_file << filenameD << std::endl;
    }
  }

} // End making day plots.