AnaPlotMaker.h File Reference

#include <fstream>
#include <TH1F.h>
#include <TH2F.h>
#include <TGraph.h>
#include <TGraphErrors.h>
#include <TTree.h>
#include <TCanvas.h>
#include <TFile.h>
#include <TPaveText.h>
#include <TLegend.h>
#include <TStyle.h>
#include <TSystem.h>
#include <TDatime.h>
#include <TMultiGraph.h>
#include <TGaxis.h>
#include <TPaveStats.h>
#include <TSQLServer.h>
#include <TSQLResult.h>
#include <TSQLStatement.h>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <string>
#include <cstring>

Go to the source code of this file.

Functions
void	eff_plots (unsigned int, string, string, string, string)

Variables
ifstream	inFile
const unsigned int	Npar = 5
TDatime *	TNowGMT
TDatime *	TNow
TDatime *	SRtime
Int_t	GMToffset
unsigned int	run
unsigned int	subrun
UInt_t	par
UInt_t	LastRun
UInt_t	LastSR
UInt_t	Last_Run [Npar]
UInt_t	Last_SR [Npar]
TH2F *	hHitsVwcellDay [Npar]
TH2F *	hCellsVwcellDay [Npar]

Function Documentation

void eff_plots	(	unsigned	int,
		string	,
		string	,
		string	,
		string
	)

Definition at line 3332 of file AnaPlotMaker.C.

References demo5::c1, demo5::c2, chisquared::c3, chisquared::c4, getBrightness::cell, nd_projection_maker::eff, allTimeWatchdog::endl, copy_new_pkl_files::f_name, GMToffset, hCellsVwcellDay, hHitsVwcellDay, Last_Run, Last_SR, new_module::module, getGoodRuns4SAM::n, NDAPDHVSetting::plane, PLOT_DIR, std::sqrt(), UpdateText, PandAna.Demos.tute_pid_validation::var, POTSpillRate::view, and w.

Referenced by AnaPlotMaker().

                                                                                             {
  

  // Output file list of low efficiency modules.
  ofstream modules;

  int ncellbins_local = 344065;
  double cellmin_local = -0.5;
  double cellmax_local = 344064.5;
  int nwbins_local = 32;
  double wmin_local = -800;
  double wmax_local = 800;

  double wavglength = 200.0; // low w length to average cell/module efficiencies over
  int wbinlow = 1;
  int wbinhigh = 5;

  char f_name[56];
  char hmod_name[56];

  // Defing binning for histograms that will be made here
  int neffbins = 101;
  double effmin = -0.05;
  double effmax = 1.05;
  int nplbins = 897;
  double plmin = -0.5;
  double plmax = nplbins-0.5;
  int nmodbins = 12;
  double modmin = 0;
  double modmax = 12;

//   TH1F* mod_Eff[12];
//   for(int imod=0;imod<12;imod++){
//     sprintf(hmod_name,"mod%ieffperplane",imod);
//     mod_Eff[imod] = new TH1F(hmod_name,hmod_name,nplbins,0,nplbins);
//   }

  TH2F* moduleEfficiencyAverage = new TH2F("moduleEfficiencyAverage",";Plane;Module",
                                           nplbins,plmin,plmax,nmodbins,modmin,modmax);
  TH2F* moduleEfficiencyRMS = new TH2F("moduleEfficiencyRMS",";Plane;Module",
                                       nplbins,plmin,plmax,nmodbins,modmin,modmax);
  TH2F* moduleEfficiencyAverageXView = new TH2F("moduleEfficiencyAverageXView",";Plane;Module",
                                                nplbins,plmin,plmax,nmodbins,modmin,modmax);
  TH2F* moduleEfficiencyRMSXView = new TH2F("moduleEfficiencyRMSXView",";Plane;Module",
                                            nplbins,plmin,plmax,nmodbins,modmin,modmax);

  TH2F* moduleEfficiencyAverageYView = new TH2F("moduleEfficiencyAverageYView",";Plane;Module",
                                                nplbins,plmin,plmax,nmodbins,modmin,modmax);
  TH2F* moduleEfficiencyRMSYView = new TH2F("moduleEfficiencyRMSYView",";Plane;Module",
                                            nplbins,plmin,plmax,nmodbins,modmin,modmax);

  // Histogram of efficiency for a module
  TH1F* modeff = new TH1F("modeff",";Efficiency;NCells",neffbins,effmin,effmax);

  // Histogram of efficiency vs w for one cell
  TH1F* eff = new TH1F("eff","",nwbins_local,wmin_local,wmax_local);

  int currmodule = 0;
  int currplane = 0;
  int currview = 0; // View convention y = 0, x = 1

  // loop over every cell in each partition
  if(par_eff==1 && period=="Month"){
    char lowEffFile[150];
    sprintf(lowEffFile, "%s/low_eff_modules.txt", PLOT_DIR.c_str());
    modules.open(lowEffFile);
    modules<<"Partition "<<par_eff<<std::endl;
    modules<<std::endl;
  }
  for(int icell = 1; icell <= ncellbins_local; ++icell){

    // convert this to a plane and cell number
    int cell = (icell-1)%384;
    int plane = (icell-1-cell)/384;
    int modcell = cell%32;
    int module = (cell-modcell)/32;
    int view = plane%2; // view = 0 y, view = 1 x

    if(plane != currplane || module != currmodule){
      // fill in the histogram of average efficiency and rms of average for the last module
      double effavg = modeff->GetMean();
      double effrms = modeff->GetRMS();
      //double effavg_far = modeff_far->GetMean();
      int planebin = moduleEfficiencyAverage->GetXaxis()->FindFixBin(currplane);
      int modulebin = moduleEfficiencyAverage->GetYaxis()->FindFixBin((double)currmodule + 0.5);

      if(par_eff==1 && period=="Month"){
        if(effavg<0.3){
          modules<<"module:  "<<modulebin<<"  plane:  "<<planebin<<"  efficiency:  "<<effavg<<std::endl;
        }
      }
      //mod_Eff[modulebin-1]->Fill(planebin,effavg);
      moduleEfficiencyAverage->SetBinContent(planebin,modulebin,effavg);
      moduleEfficiencyRMS->SetBinContent(planebin,modulebin,effrms);
      if(currview){
        // x view
        moduleEfficiencyAverageXView->SetBinContent(planebin,modulebin,effavg);
        moduleEfficiencyRMSXView->SetBinContent(planebin,modulebin,effrms);
      }
      else{
        // y view
        moduleEfficiencyAverageYView->SetBinContent(planebin,modulebin,effavg);
        moduleEfficiencyRMSYView->SetBinContent(planebin,modulebin,effrms);
      }
      // reset things for the new module
      modeff->Reset("ICES");
      currmodule = module;
      currplane = plane;
      currview = view;
    }
    

    // need to switch the module in cell for x view
 
    if(view){ modcell = 31-modcell; }

    // loop over every w bin
    double totcell = 0.0;
    double tothit = 0.0;

    // Calculate the cell efficiency vs w
    for(int iw = 1; iw <= nwbins_local; ++iw){
      // get the efficiency at this cell/w, eff = k/n
      double k = hHitsVwcellDay[par_eff]->GetBinContent(icell,iw);
      double n = hCellsVwcellDay[par_eff]->GetBinContent(icell,iw);

      totcell+=n;
      tothit+=k;
      double ef = k/n;
      if(n == 0){ ef = 0; }

      // calculate the variance of this
      double var = (k+1)*(k+2)/((n+2)*(n+3)) - (k+1)*(k+1)/((n+2)*(n+2));
      
      double w = hHitsVwcellDay[par_eff]->GetYaxis()->GetBinCenter(iw);
      
      eff->SetBinContent(iw,ef);
      eff->SetBinError(iw,sqrt(var));
           
    }


    // Average over the last 200 cm of the cell
    double efftot = eff->Integral(wbinlow,wbinhigh);
    double navgbins = wbinhigh-wbinlow+1.0;

    modeff->Fill(efftot/navgbins);
           
  }
  if(par_eff==1 && period=="Month") modules.close();
  
  char title_buff[128];

  gStyle->SetTimeOffset(0);
  string t_title;
  if(period == "Day") t_title = "(past 24 hrs.)";
  if(period == "Week") t_title = "(past week)";
  if(period == "Month") t_title = "(past month)";
  if(period == "6Months") t_title = "(past 26 weeks)";
  if(period == "Month") t_title = "(past 52 weeks)";

  TDatime *Ttemp = new TDatime;  // finish time
  int Xfin = Ttemp->Convert() - GMToffset;
  TDatime *Tfinish = new TDatime(Xfin);  // finish time
  TPaveText *UpdateText = new TPaveText(0.1, 0.0, 0.5, 0.05, "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", Last_Run[par_eff], Last_SR[par_eff]);
  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);

  TCanvas* c1 = new TCanvas("c1","c1",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"Average Module Efficiency %s - partition %.1u",t_title.c_str(),par_eff);
  moduleEfficiencyAverage->SetTitle(title_buff);
  moduleEfficiencyAverage->Draw("colz");
  //moduleEfficiencyAverage->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats1 = (TPaveStats*)moduleEfficiencyAverage->FindObject("stats");
//   stats1->SetX1NDC(0.65);
//   stats1->SetY1NDC(0.65);
//   stats1->SetX2NDC(0.90);
//   stats1->SetY2NDC(0.90);  
  c1->Modified();
  c1->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uAvgModuleEff%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c1->Print(f_name);
  delete c1;
  c1 = 0;

  //////////////////////////////////////////////////////////////////////
  // Creating list of low efficiency modules and corresponding planes.//
  //////////////////////////////////////////////////////////////////////
  /* TFile *test_out = new TFile("testout.root","RECREATE");
  modules.open("low_eff_modules.txt");
  //int Nbinsy = moduleEfficiencyAverage->GetNbinsY();
  //int Nbinsx = moduleEfficiencyAverage->GetNbinsX();
  if(det=="FarDet" && period=="Day"){
    for(int i=1;i<257;i++){
      for(int j=1;j<13;j++){
        //Int_t bin_num = moduleEfficiencyAverage->FindFixBin(double(i)+0.5,j);
        //Int_t bin_num = moduleEfficiencyAverage->GetBin(i+1,j+1);
        Double_t eff_mod = moduleEfficiencyAverage->GetBinContent(i,j);
        //std::cout<<bin_num<<"      "<<eff_mod<<std::endl;
        //if(eff_mod<0.6){
        //modules<<" plane:  "<<i<<"   module:   "<<j<<"   efficiency:  "<<eff_mod<<std::endl;
        modules << moduleEfficiencyAverage->GetBinContent(i,j)<<std::endl;
          //}
      }
    }
  }
  modules.close();
  test_out->cd();
  for(int ii=0;ii<12;ii++){
  mod_Eff[ii]->Write();
  }
  //moduleEfficiencyAverage->Write();
  //test_out->Close();
  */

  TCanvas* c2 = new TCanvas("c2","c2",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"%s Average Module Efficiency RMS %s - partition %.1u",trigname.c_str(),t_title.c_str(),par_eff);
  moduleEfficiencyRMS->SetTitle(title_buff);
  moduleEfficiencyRMS->Draw("colz");
  //moduleEfficiencyRMS->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats2 = (TPaveStats*)moduleEfficiencyRMS->FindObject("stats");
//   stats2->SetX1NDC(0.65);
//   stats2->SetY1NDC(0.65);
//   stats2->SetX2NDC(0.90);
//   stats2->SetY2NDC(0.90); 
  c2->Modified();
  c2->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uModuleEffRMS%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c2->Print(f_name);
  delete c2;
  c2 = 0;

  TCanvas* c3 = new TCanvas("c3","c3",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"%s Average Module Efficiency In X-view %s - partition %.1u",trigname.c_str(),t_title.c_str(),par_eff);
  moduleEfficiencyAverageXView->SetTitle(title_buff);
  moduleEfficiencyAverageXView->Draw("colz");
  moduleEfficiencyAverageXView->RebinX(2);
  //moduleEfficiencyAverageXView->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats3 = (TPaveStats*)moduleEfficiencyAverageXView->FindObject("stats");
//   stats3->SetX1NDC(0.65);
//   stats3->SetY1NDC(0.65);
//   stats3->SetX2NDC(0.90);
//   stats3->SetY2NDC(0.90); 
  c3->Modified();
  c3->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uAvgModuleXEff%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c3->Print(f_name);
  delete c3;
  c3 = 0;

  TCanvas* c4 = new TCanvas("c4","c4",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"%s Average Module Efficiency In X-view RMS %s - partition %.1u",trigname.c_str(),t_title.c_str(),par_eff);
  moduleEfficiencyRMSXView->SetTitle(title_buff);
  moduleEfficiencyRMSXView->Draw("colz");
  moduleEfficiencyRMSXView->RebinX(2);
  //moduleEfficiencyRMSXView->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats4 = (TPaveStats*)moduleEfficiencyRMSXView->FindObject("stats");
//   stats4->SetX1NDC(0.65);
//   stats4->SetY1NDC(0.65);
//   stats4->SetX2NDC(0.90);
//   stats4->SetY2NDC(0.90); 
  c4->Modified();
  c4->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uModuleXEffRMS%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c4->Print(f_name);
  delete c4;
  c4 = 0;

  TCanvas* c5 = new TCanvas("c5","c5",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"%s Average Module Efficiency In Y-view %s - partition %.1u",trigname.c_str(),t_title.c_str(),par_eff);
  moduleEfficiencyAverageYView->SetTitle(title_buff);
  moduleEfficiencyAverageYView->Draw("colz");
  moduleEfficiencyAverageYView->RebinX(2);
  //moduleEfficiencyAverageYView->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats5 = (TPaveStats*)moduleEfficiencyAverageYView->FindObject("stats");
//   stats5->SetX1NDC(0.65);
//   stats5->SetY1NDC(0.65);
//   stats5->SetX2NDC(0.90);
//   stats5->SetY2NDC(0.90); 
  c5->Modified();
  c5->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uAvgModuleYEff%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c5->Print(f_name);
  delete c5;
  c5 = 0;

  TCanvas* c6 = new TCanvas("c6","c6",1200,800);
  gStyle->SetOptStat(0);
  sprintf(title_buff,"%s Average Module Efficiency In Y-view RMS %s - partition %.1u",trigname.c_str(),t_title.c_str(),par_eff);
  moduleEfficiencyRMSYView->SetTitle(title_buff);
  moduleEfficiencyRMSYView->Draw("colz");
  moduleEfficiencyRMSYView->RebinX(2);
  //moduleEfficiencyRMSYView->GetXaxis()->SetRangeUser(0.,300.);
  gPad->Update();
//   TPaveStats *stats6 = (TPaveStats*)moduleEfficiencyRMSYView->FindObject("stats");
//   stats6->SetX1NDC(0.65);
//   stats6->SetY1NDC(0.65);
//   stats6->SetX2NDC(0.90);
//   stats6->SetY2NDC(0.90); 
  c6->Modified();
  c6->Update();
  UpdateText->Draw();
  sprintf(f_name,"%s/%s-%s-P%.1uModuleYEffRMS%s.png",PLOT_DIR.c_str(),det.c_str(),trig.c_str(),par_eff,period.c_str());
  c6->Print(f_name);
  delete c6;
  c6 = 0;
            
}

Variable Documentation

Int_t GMToffset

Definition at line 43 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), CheckTrueDate(), DDTOnMonPlotMaker(), eff_plots(), GoodDataSelDate(), GoodDataSelDuration(), GoodDataSelEmptySpillFrac(), GoodDataSelGoodDB(), GoodDataSelLiveTime(), GoodDataSelMipRate(), GoodDataSelNumSlices(), GoodDataSelRuns(), GoodDataSelTimingPeakEnd(), GoodDataSelTimingPeakStart(), GoodDataSelTrkFrac(), Nearline_Watchdog_Plots(), OnMonPlotMaker(), PlotDB_Web(), PlotEfficiency(), plotTime(), TB_Nearline_Watchdog_Plots(), and UpdatedText().

TH2F* hCellsVwcellDay[Npar]

Definition at line 71 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), and eff_plots().

TH2F* hHitsVwcellDay[Npar]

Definition at line 70 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), and eff_plots().

ifstream inFile

Definition at line 34 of file AnaPlotMaker.h.

Referenced by AnalyzeNus18Systs(), AnalyzeSysts(), AnaPlotMaker(), append(), BestFit(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), Compare_NoExtrap(), Compare_Spectra(), ComparisonPlots_Data(), ComparisonPlots_MC(), CutFlow_NearDet(), DCMOffsetCalculator(), DDTOnMonPlotMaker(), demo2p5b(), DoTraining(), DrawCovMx(), DrawExtrapSurface(), EnergyCont_macro(), example_plot(), getData(), GetExtrapolationDefs(), ana::GetNueQuantCuts2020(), ana::GetNumuEhadFracQuantCuts2020(), ana::GetNumuPtQuantCuts2020(), getPredictions(), GetQuantileCuts(), hadEFrac_nd_data_mc_systs(), HtmlCosmic(), HtmlGenie(), hyperon_macro(), hyperon_nom_macro(), daq2raw::TestBeamRawInputDriver::loadBeamlineEvent(), ana::LoadSystsFromFile(), make_dst_cosrejbdttrain(), make_eff_plots_areaNorm(), make_prediction_extrap(), make_prediction_noextrap(), make_starplots(), make_xsec_tuning_hists_mp5(), make_xsec_wgts_2018_hists(), MakeCovMx(), ana::MakeCutVec(), MakeISysts(), meanWeight_plot(), mec_nux_tester_2020(), mec_tuning_fitter_2020(), MergePredictions(), mlpPIDelecallEetot(), mlpPIDelecallEtot(), multiverse_efficiency_plot(), ND_DataMC_Comp_Systs(), neutKEsyst(), neutronE_macro(), NoExtrap(), NuSCalculateCorr(), NuSTreeMaker(), OnMonPlotMaker(), plot_DataMCComp_numu(), plot_pi0_xcheck(), plot_pid(), daq2raw::TestBeamRawInputDriver::readBeamline(), resolution2018(), SystematicComp(), and TrimCAFs().

UInt_t Last_Run[Npar]

Definition at line 53 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), and eff_plots().

UInt_t Last_SR[Npar]

Definition at line 54 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), and eff_plots().

UInt_t LastRun

Definition at line 49 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), CheckTrueDate(), DDTOnMonPlotMaker(), Nearline_Watchdog_Plots(), NewNearlinePlots(), OnMonPlotMaker(), and TB_Nearline_Watchdog_Plots().

UInt_t LastSR

Definition at line 50 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), CheckTrueDate(), DDTOnMonPlotMaker(), Nearline_Watchdog_Plots(), OnMonPlotMaker(), and TB_Nearline_Watchdog_Plots().

const unsigned int Npar = 5

Definition at line 38 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), and OnMonPlotMaker().

UInt_t par

Definition at line 47 of file AnaPlotMaker.h.

unsigned int run

Definition at line 45 of file AnaPlotMaker.h.

TDatime* SRtime

Definition at line 42 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), CheckTrueDate(), DDTOnMonPlotMaker(), Nearline_Watchdog_Plots(), OnMonPlotMaker(), and TB_Nearline_Watchdog_Plots().

unsigned int subrun

Definition at line 46 of file AnaPlotMaker.h.

TDatime* TNow

Definition at line 41 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), APDPlots(), ChanPlots(), CheckTrueDate(), DCSPlotMaker(), DDTOnMonPlotMaker(), Nearline_Watchdog_Plots(), NewNearlinePlots(), OnMonPlotMaker(), and TB_Nearline_Watchdog_Plots().

TDatime* TNowGMT

Definition at line 40 of file AnaPlotMaker.h.

Referenced by AnaPlotMaker(), CheckTrueDate(), DDTOnMonPlotMaker(), Nearline_Watchdog_Plots(), OnMonPlotMaker(), and TB_Nearline_Watchdog_Plots().