cellShifts.C File Reference

#include "TCanvas.h"
#include "TChain.h"
#include "TCut.h"
#include "TLine.h"
#include "TFile.h"
#include "TF2.h"
#include "TF1.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TStyle.h"
#include <cmath>

Go to the source code of this file.

Functions
void	FitForMinimum (TH2 *h2, double &x0, double &y0)
void	cellPos ()
void	initialise (bool IsData)
void	plot_residPerViewV (bool IsData)
void	plot_sumcellsV (int cell)
void	plot_sumcellsZ (int cell)
void	plot_residPerViewZ (bool IsData)
void	plot_resid (bool IsData, int p, int c)
void	plot_residV (bool IsData, int p, int c)
void	plot_residZ (bool IsData, int p, int c)
void	plot_CellPos (bool IsData, int p, int c)
void	plot_residPerPlane (int Plane)
void	perView_labelV (bool IsData)
void	perView_labelZ (bool IsData)
void	plane_label (bool IsData, int p, int c)
void	SetLimits (double minmax)
void	save (bool IsData, bool Bad)

Variables
const int	pmax = 200
const int	cmax = 96
const int	nbins = 200
TFile *	file = 0
TChain *	al = 0
TH1D *	fResidV [pmax][cmax] = {}
TH1D *	fResidZ [pmax][cmax] = {}
TH2D *	hCellPos [pmax][cmax] = {}
TH2D *	p1X [pmax] = {}

Function Documentation

void cellPos

(

)

Definition at line 26 of file cellShifts.C.

References initialise(), plot_residPerViewV(), and plot_residPerViewZ().

Referenced by align::SplitTracks::analyze(), slid::ParticleIDAlg::CalcPlaneCentroid(), jmshower::RecoJMShower::GetCellTransDedx(), and jmshower::RecoJMShower::GetTrkCPlaneCell().

{
  initialise(0);        // 0 = MC, 1 = Data
  plot_residPerViewV(0); // 0 = MC, 1 = Data
  plot_residPerViewZ(0);
}

void FitForMinimum	(	TH2 *	h2,
		double &	x0,
		double &	y0
	)

Definition at line 472 of file cellShifts.C.

References func().

Referenced by plot_residPerViewV(), and plot_residPerViewZ().

                                                    {

  // for speed, only fit right around minimum

  int binX0;
  int binY0;
  int binZ0;
  h2->GetMinimumBin(binX0,binY0,binZ0);

  int fitbinX1 = binX0 - 2;
  int fitbinX2 = binX0 + 2;
  int fitbinY1 = binY0 - 2;
  int fitbinY2 = binY0 + 2;

  if (fitbinX1<1) {
    fitbinX2 += (1-fitbinX1);
    fitbinX1 += (1-fitbinX1);
  }
  if (fitbinX2>h2->GetNbinsX()) {
    fitbinX1 -= (fitbinX2-h2->GetNbinsX());
    fitbinX2 -= (fitbinX2-h2->GetNbinsX());
  }
  if (fitbinY1<1) {
    fitbinY2 += (1-fitbinY1);
    fitbinY1 += (1-fitbinY1);
  }
  if (fitbinY2>h2->GetNbinsY()) {
    fitbinX1 += (fitbinX2-h2->GetNbinsY());
    fitbinX2 += (fitbinX2-h2->GetNbinsY());
  }

  TF2 func("func","[2]*(x-[0])*(x-[0])+[3]*(x-[0])*(y-[1])+[4]*(y-[1])*(y-[1])+[5]",
           h2->GetXaxis()->GetBinLowEdge(fitbinX1),
           h2->GetXaxis()->GetBinUpEdge (fitbinX2),
           h2->GetYaxis()->GetBinLowEdge(fitbinY1),
           h2->GetYaxis()->GetBinUpEdge (fitbinY2));

  func.SetParameters(
                     h2->GetXaxis()->GetBinCenter(binX0),
                     h2->GetYaxis()->GetBinCenter(binY0),
                     1000,-10,1000,h2->GetBinContent(binX0,binY0));

  h2->Fit(&func,"q0RN","goff");
  x0 = func.GetParameter(0);
  y0 = func.GetParameter(1);
}

void initialise

(

bool

IsData

)

Definition at line 33 of file cellShifts.C.

References getBrightness::cell, cmax, om::cout, allTimeWatchdog::endl, file, fResidV, fResidZ, hCellPos, findDuplicateFiles::key, NDAPDHVSetting::plane, pmax, and PandAna.Demos.pi0_spectra::tot.

Referenced by cellPos().

{
  file = TFile::Open("/nova/app/users/gsdavies/art/align/align_1X_cosmicsMC.root");
  file->cd("alignment");

  int tot = 0;

  // Read in histograms from file
  for (int plane = 0; plane < pmax; plane++){
    for (int cell =0; cell < cmax; cell++){
      TString hname;
      TString hname2;
      TString hname3;
      if(IsData) { 
        hname = Form("residualmcV_%03d_%03d",plane,cell);
        hname2 = Form("residualmcZ_%03d_%03d",plane,cell); 
        hname3 = Form("hCellPos_%03d_%03d",plane,cell);  
      }
      else { 
        hname = Form("residualmcV_%03d_%03d",plane,cell);
        hname2 = Form("residualmcZ_%03d_%03d",plane,cell); 
        hname3 = Form("hCellPos_%03d_%03d",plane,cell);  
      }
      TKey *key = gDirectory->FindKey(hname);
      TKey *key2 = gDirectory->FindKey(hname2);
      TKey *key3 = gDirectory->FindKey(hname3);
      if(key==0 || key2 == 0 || key3 == 0) { 
        //std::cerr  << "!!Histogram does not exist!!" << endl;
        tot++;
        continue;
      }
      else {
        fResidV[plane][cell] = (TH1D *)gDirectory->Get(hname);
        fResidZ[plane][cell] = (TH1D *)gDirectory->Get(hname2);
        hCellPos[plane][cell] = (TH2D *)gDirectory->Get(hname3);
      }
    }
    std::cout << "Processing plane : " << plane << std::endl;
  }
  std::cout << "Total histograms that do not exist: " << tot << std::endl;
}

void perView_labelV

(

bool

IsData

)

Definition at line 360 of file cellShifts.C.

References kBlue.

Referenced by plot_residPerViewV().

{
  
  mX->cd();
  TLatex* xz;
  if(IsData) { xz = new TLatex(.9, .95, "XZ-View [V](Data)"); }
  else { xz = new TLatex(.9, .95, "XZ-View [V](MC)"); }
  xz->SetTextColor(kBlue);
  xz->SetNDC();
  xz->SetTextSize(2/30.);
  xz->SetTextAlign(32);
  xz->Draw();
  
  mY->cd();
  TLatex* yz;
  if(IsData) { yz = new TLatex(.9, .95, "YZ-View [V](Data)"); }
  else { yz = new TLatex(.9, .95, "YZ-View [V](MC)"); }
  yz->SetTextColor(kBlue);
  yz->SetNDC();
  yz->SetTextSize(2/30.);
  yz->SetTextAlign(32);
  yz->Draw();
}

void perView_labelZ

(

bool

IsData

)

Definition at line 383 of file cellShifts.C.

References kBlue.

Referenced by plot_residPerViewZ().

{
  
  mXz->cd();
  TLatex* xz;
  if(IsData) { xz = new TLatex(.9, .95, "XZ-View [Z](Data)"); }
  else { xz = new TLatex(.9, .95, "XZ-View [Z](MC)"); }
  xz->SetTextColor(kBlue);
  xz->SetNDC();
  xz->SetTextSize(2/30.);
  xz->SetTextAlign(32);
  xz->Draw();
  
  mYz->cd();
  TLatex* yz;
  if(IsData) { yz = new TLatex(.9, .95, "YZ-View [Z](Data)"); }
  else { yz = new TLatex(.9, .95, "YZ-View [Z](MC)"); }
  yz->SetTextColor(kBlue);
  yz->SetNDC();
  yz->SetTextSize(2/30.);
  yz->SetTextAlign(32);
  yz->Draw();
}

void plane_label	(	bool	IsData,
		int	p,
		int	c
	)

Definition at line 407 of file cellShifts.C.

References make_syst_table_plots::h, hh, kBlue, and plotROC::title.

Referenced by plot_resid().

{
  h->cd();
  if(IsData) { 
    TString title = Form("Plane %03d, Cell %03d (Data)",p,c);
    TLatex* hh = new TLatex(.9, .95, title);
  }  
  else { 
    TString title = Form("Plane %03d, Cell %03d (MC)",p,c);
    TLatex* hh = new TLatex(.9, .95, title);
  }  
  hh->SetTextColor(kBlue);
  hh->SetNDC();
  hh->SetTextSize(2/30.);
  hh->SetTextAlign(32);
  hh->Draw();

}

void plot_CellPos	(	bool	IsData,
		int	p,
		int	c
	)

Definition at line 293 of file cellShifts.C.

References plot_validation_datamc::c, h3, and hCellPos.

{
  TCanvas* h3 = new TCanvas("h3","h3");
  h3->SetBorderMode(0);
  h3->Clear();
  h3->cd();
  hCellPos[p][c]->Draw("colz");// bin 111,9 -> plane 110, cell 8
  
  //plane_label(IsData,p,c);

}  

void plot_resid	(	bool	IsData,
		int	p,
		int	c
	)

Definition at line 257 of file cellShifts.C.

References plot_validation_datamc::c, make_syst_table_plots::h, and plane_label().

{
  TCanvas* h = new TCanvas("h","h");
  h->SetBorderMode(0);
  h->Clear();
  h->cd();
  fResid[p][c]->Draw();// bin 111,9 -> plane 110, cell 8
  
  plane_label(IsData,p,c);

}

void plot_residPerPlane

(

int

Plane

)

Definition at line 304 of file cellShifts.C.

References getBrightness::cell, cmax, fResidV, fResidZ, calib::j, plot_validation_datamc::p1, p1X, plot_validation_datamc::p2, NDAPDHVSetting::plane, and pmax.

                                  {
  TCanvas* p1 = new TCanvas("p1","p1");
  p1->SetBorderMode(0);
  p1->Divide(10,20);

  TCanvas* p2 = new TCanvas("p1","p1");
  p2->SetBorderMode(0);
  p2->Divide(10,20);
  for (int plane = 0; plane < pmax; plane++){
    TString name = Form("plane_%03d",plane);
    p1X[plane] = new TH2D(name,";cell;residual v (cm)",96,0.,95,1000,-50,50);
    p2X[plane] = new TH2D(name,";cell;residual z (cm)",96,0.,95,1000,-50,50);
    double resid;
    for (int cell = 0; cell < cmax; cell++){  
      if(fResidV[plane][cell]){  
        if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
          
          for (int j = 0; j < 1001; j++){
            resid = (j/5)-100; 
            p1X[plane]->Fill(cell,resid,fResidV[plane][cell]->GetBinContent(j));
          }  
        }
        else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
          for (int j = 0; j < 1001; j++){
            resid = (j/5)-100; 
            p1X[plane]->Fill(cell,resid,fResidV[plane][cell]->GetBinContent(j));
          }  
        }
      }
      if(fResidZ[plane][cell]){
        if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
          
          for (int j = 0; j < 1001; j++){
            resid = (j/5)-100; 
            p2X[plane]->Fill(cell,resid,fResidZ[plane][cell]->GetBinContent(j));
          }  
        }
        else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
          for (int j = 0; j < 1001; j++){
            resid = (j/5)-100; 
            p2X[plane]->Fill(cell,resid,fResidZ[plane][cell]->GetBinContent(j));
          }  
        }
      }  
      else {
        //std::cerr  << "!!Histogram not filled or does not exist!!" << endl;
      }
    }
  }
  p1->cd(); 
  p1X[Plane]->Draw("colz");
  p1->cd();
  p2X[Plane]->Draw("colz");

}

void plot_residPerViewV

(

bool

IsData

)

Definition at line 76 of file cellShifts.C.

References getBrightness::cell, cmax, om::cout, allTimeWatchdog::endl, FitForMinimum(), hCellPos, MECModelEnuComparisons::i, calib::j, nbins, perView_labelV(), NDAPDHVSetting::plane, and pmax.

Referenced by cellPos().

{
  TCanvas* mX = new TCanvas("mX","mX");
  
  mX->SetBorderMode(0);
  TCanvas* mY = new TCanvas("mY","mY");
  mY->SetBorderMode(0);
  
  TH2D *meanX = new TH2D("meanX",";plane;cell",nbins,0.,199,65,0.,64); 
  TH2D *meanY = new TH2D("meanY",";plane;cell",nbins,0.,199,96,0.,95); 
  for (int plane = 0; plane < pmax; plane++){
    cout << "Current plane in v mode: " << plane << endl;
    for (int cell = 0; cell < cmax; cell++){
      
      if(hCellPos[plane][cell]){

        double v0, z0;
        
        FitForMinimum(hCellPos[plane][cell],v0,z0);

        // hack for boundary of two vertical planes in NDOS
        if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
          //std::cout << "Plane/cell " << plane << ", " << cell << std::endl;
          meanX->Fill(plane,cell,v0);
        }
        else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
          meanY->Fill(plane,cell,v0);
        }
      }
      else {
      }
    }
  }
  
  // Clean up channels that don't exist in each view so they show up blank!
  for(int i =0; i< nbins+1; i++){
    for(int j=0; j< 66;j++){
      if(meanX->GetBinContent(i,j) == 0){
        meanX->SetBinContent(i,j,-999);
       }
     }
    for(int k=0; k< 97;k++){
      if(meanY->GetBinContent(i,k) == 0){
        meanY->SetBinContent(i,k,-999);
       }
     } 
  }
  
  mX->Clear();
  mX->cd();      

  meanX->Draw("colz");
  mY->Clear();
  mY->cd();

  meanY->Draw("colz");
  perView_labelV(IsData);  
}

void plot_residPerViewZ

(

bool

IsData

)

Definition at line 198 of file cellShifts.C.

References getBrightness::cell, cmax, om::cout, allTimeWatchdog::endl, FitForMinimum(), hCellPos, MECModelEnuComparisons::i, calib::j, nbins, perView_labelZ(), NDAPDHVSetting::plane, and pmax.

Referenced by cellPos().

{
  TCanvas* mXz = new TCanvas("mXz","mXz");
  
  mXz->SetBorderMode(0);
  TCanvas* mYz = new TCanvas("mYz","mYz");
  mYz->SetBorderMode(0);
  
  TH2D *meanXz = new TH2D("meanXz",";plane;cell",nbins,0.,199,65,0.,64); 
  TH2D *meanYz = new TH2D("meanYz",";plane;cell",nbins,0.,199,96,0.,95);

  for (int plane = 0; plane < pmax; plane++){
    cout << "Current plane in z mode: " << plane << endl;
    for (int cell = 0; cell < cmax; cell++){
      
      if(hCellPos[plane][cell]){

        double v0, z0;
        FitForMinimum(hCellPos[plane][cell],v0,z0);

        // hack for boundary of two vertical planes in NDOS
        if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
          //std::cout << "Plane/cell " << plane << ", " << cell << std::endl;
          meanXz->Fill(plane,cell,z0);
        }
        else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
          meanYz->Fill(plane,cell,z0);
        } 
      }
      else {
      }
    }
  }
  
  // Clean up channels that don't exist in each view so they show up blank!
  for(int i =0; i< nbins+1; i++){
    for(int j=0; j< 66;j++){
      if(meanXz->GetBinContent(i,j) == 0){
        meanXz->SetBinContent(i,j,-999);
       }
     }
    for(int k=0; k< 97;k++){
      if(meanYz->GetBinContent(i,k) == 0){
        meanYz->SetBinContent(i,k,-999);
       }
     } 
  }
  mXz->Clear();
  mXz->cd();      

  meanXz->Draw("colz");
  mYz->Clear();
  mYz->cd();

  meanYz->Draw("colz");
  perView_labelZ(IsData);
  
}

void plot_residV	(	bool	IsData,
		int	p,
		int	c
	)

Definition at line 270 of file cellShifts.C.

References plot_validation_datamc::c, fResidV, and h1.

{
  TCanvas* h1 = new TCanvas("h1","h1");
  h1->SetBorderMode(0);
  h1->Clear();
  h1->cd();
  fResidV[p][c]->Draw();// bin 111,9 -> plane 110, cell 8
  
  //plane_label(IsData,p,c);

}

void plot_residZ	(	bool	IsData,
		int	p,
		int	c
	)

Definition at line 282 of file cellShifts.C.

References plot_validation_datamc::c, fResidZ, and h2.

{
  TCanvas* h2 = new TCanvas("h2","h2");
  h2->SetBorderMode(0);
  h2->Clear();
  h2->cd();
  fResidZ[p][c]->Draw();// bin 111,9 -> plane 110, cell 8
  
  //plane_label(IsData,p,c);

}

void plot_sumcellsV

(

int

cell

)

Definition at line 135 of file cellShifts.C.

References fResidV, NDAPDHVSetting::plane, and pmax.

{
  TCanvas* sumXV = new TCanvas("sumXV","sumXV");
  sumXV->SetBorderMode(0);
  TCanvas* sumYV = new TCanvas("sumYV","sumYV");
  sumYV->SetBorderMode(0);
  TString nameXV = Form("rsumXV_cell_%03d",cell);
  TString nameYV = Form("rsumYV_cell_%03d",cell);
  TString strV = Form(" sum cell %03d",cell);
  TH1D *SumXV = new TH1D(nameXV,";Residual (cm)"+strV,1000,-50,50);
  TH1D *SumYV = new TH1D(nameYV,";Residual (cm)"+strV,1000,-50,50);

  for (int plane = 0; plane < pmax; plane++){
    if(fResidV[plane][cell]){
      // hack for boundary of two vertical planes in NDOS
      if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
       SumXV->Add(fResidV[plane][cell]);
      }
      else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
       SumYV->Add(fResidV[plane][cell]);
      } 
    }
    else {
    }
  }
  sumXV->cd();
  sumXV->Clear();
  SumXV->Draw();
  
  sumYV->cd();
  sumYV->Clear();
  SumYV->Draw();    
}

void plot_sumcellsZ

(

int

cell

)

Definition at line 169 of file cellShifts.C.

References fResidZ, NDAPDHVSetting::plane, and pmax.

{
  TCanvas* sumXZ = new TCanvas("sumXZ","sumXZ");
  sumXZ->SetBorderMode(0);
  TCanvas* sumYZ = new TCanvas("sumYZ","sumYZ");
  sumYZ->SetBorderMode(0);
  TString nameXZ = Form("rsumXZ_cell_%03d",cell);
  TString nameYZ = Form("rsumYZ_cell_%03d",cell);
  TString strZ = Form(" sum cell %03d",cell);
  TH1D *SumXZ = new TH1D(nameXZ,";Residual (cm)"+strZ,1000,-50,50);
  TH1D *SumYZ = new TH1D(nameYZ,";Residual (cm)"+strZ,1000,-50,50);

  for (int plane = 0; plane < pmax; plane++){
    if(fResidZ[plane][cell]){
      // hack for boundary of two vertical planes in NDOS
      if ((plane < 155 && plane % 2 == 0) || (plane >= 155 && plane % 2 != 0)){
       SumXZ->Fill(fResidZ[plane][cell]);
      }
      else if ((plane <= 155 && plane % 2 != 0) || (plane > 155 && plane % 2 == 0)){
       SumYZ->Fill(fResidZ[plane][cell]);
      } 
    }
    else {
    }
  }
  SumXZ->Draw();
  SumYZ->Draw();    
}

void save	(	bool	IsData,
		bool	Bad
	)

Definition at line 453 of file cellShifts.C.

Referenced by median_contours(), and ana::process_list().

{
  if(IsData) {
    mX->SaveAs("xzView-resid_data.png");
    mY->SaveAs("yzView-residual_data.png");
  }
  else { 
    if(Bad) {
      mX->SaveAs("xzView-resid_mcbad.png");
      mY->SaveAs("yzView-residual_mcbad.png");
    }
    else {
      mX->SaveAs("xzView-resid_mc.png");
      mY->SaveAs("yzView-residual_mc.png");
    }  
  }
}  

void SetLimits

(

double

minmax

)

Definition at line 429 of file cellShifts.C.

{
  mX->cd();
  meanX->SetMaximum(minmax);
  meanX->SetMinimum(-minmax);
  meanX->Draw("colz");
  
  mY->cd();
  meanY->SetMaximum(minmax);
  meanY->SetMinimum(-minmax);
  meanY->Draw("colz");
  
  mXz->cd();
  meanXz->SetMaximum(minmax);
  meanXz->SetMinimum(-minmax);
  meanXz->Draw("colz");
  
  mYz->cd();
  meanYz->SetMaximum(minmax);
  meanYz->SetMinimum(-minmax);
  meanYz->Draw("colz");
}

Variable Documentation

TChain* al = 0

Definition at line 18 of file cellShifts.C.

const int cmax = 96

Definition at line 14 of file cellShifts.C.

Referenced by initialise(), dt::Cand::MinMaxLines(), dt::Cand::NSegs(), plot_residPerPlane(), plot_residPerViewV(), plot_residPerViewZ(), dt::Cand::PopSegHelper(), novaddt::calib::BufferManager< bufferType, Nbuffers >::SetCounterMax(), slicer_production_validation(), and dt::Cand::TryAddSeg().

TFile* file = 0

Definition at line 17 of file cellShifts.C.

Referenced by compareEvents(), countEvents(), d4sigma_plot(), DCMOffsetCalculator(), doUnfolding(), fake_future_data(), fiducial_accounting(), fill_simple(), bp::FindFiles(), mcmc::FullFilename(), test::FullFilename(), futureSig_reach_singlePOTcombo_syst(), generate_sk_flux_histograms(), genie_preds_make(), get_data_and_cosmic(), get_data_histogram(), get_fdmc_spectra(), get_files(), get_numu_data_histogram(), get_spectra_2dplots(), GetManualSubruns(), ana::GetNueData2020(), GetNumuData2018(), GetNumuData2019(), ana::GetNumuData2020(), GetPartList(), GetPartList_4DB(), getPlotsFromNtupleFD(), getPlotsFromNtupleND(), GetProcessed(), genie::utils::xml::GetTFile(), genie::utils::xml::GetTH1D(), genie::utils::xml::GetTH1F(), genie::utils::xml::GetTH2D(), getTimePeakTotal(), genie::utils::xml::GetTVectorD(), goodness_of_fit(), HadFit2sp(), HadFit4sp(), HadFitCC4sp(), ana::MultiHistContainer::Hist(), initialise(), LoadPPFX(), main(), make_calib_hists(), make_decomp(), make_estimate_energy(), make_esyst_hists(), make_extrap_figure_hists(), make_michel(), make_nueFDprediction_kinematics(), make_nueFDprediction_kinematics_FHC(), make_nueFDprediction_kinematics_FHC_pTExtrap(), make_nueFDprediction_kinematics_FHC_REW(), make_nueFDprediction_kinematics_FHC_REW_pTExtrap(), make_nueFDprediction_kinematics_REW(), make_nueFDprediction_kinematics_RHC(), make_nueFDprediction_kinematics_RHC_pTExtrap(), make_nueFDprediction_kinematics_RHC_REW(), make_nueFDprediction_kinematics_RHC_REW_pTExtrap(), make_prediction(), make_prediction_rhc(), make_predictions_systs(), make_predictions_systs_fakedata(), make_quantile_cuts_predictions(), make_starplots(), make_surfprof(), ana::MakeCosOutOfTimeList(), MakeDetCfgCSV(), MakeDQValidityCSV(), MakeDQValidityCSVND(), ana::MakeNumiUnblindList(), makePlots(), MakePlots(), MakeProfileHadEFD(), MakeProfileHadEND(), MakeProfileMuEFD(), MakeProfileMuEND(), MakeRunList(), ana::MakeUnblindList(), MakeUnblindList(), MakingProfile(), MakingProfileActCatcherND(), MakingProfileAllCatcherND(), MakingProfileHadCC(), MakingProfileHadCCND(), MakingProfileHadNonQE(), MakingProfileHadNonQEND(), MakingProfileHadQE(), MakingProfileHadQEND(), MakingProfileND(), mf::MaybeLogger_< SEV, VERBATIM >::MaybeLogger_(), mec_tuning_preds_2020(), modularextrap_demo_nue(), modularextrap_demo_numu(), mrbrem_get_initial_spectra(), mrbrem_get_reweighted_spectra(), MRDiFStudy_FHC_Step1(), MRDiFStudy_RHC_Step1(), mre_blessed(), mre_comp_split(), MuonFit4sp(), MuonFitLin(), neut_code_from_rootracker(), om::HistoTable::NewLineChar(), nuebar_signif(), nus17_fiducial_accounting(), NuSTreeMaker(), TMVA::TMVAGlob::OpenFile(), INukeNucleonCorr::OutputFiles(), PhotonSim(), PhotonSim_mp(), plot_cont_expfriends(), plot_contprof(), plot_datapredictions(), plot_NDvsFD_weights(), plot_NDvsFD_weights_FHC(), plot_NDvsFD_weights_RHC(), plot_nueFD_Signal_REWvsNOM_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_FHC(), plot_nueFD_Signal_REWvsNOM_pTExtrap_RHC(), plot_nueFD_Signal_REWvsNOM_RHC(), plot_prediction_systs(), plot_predictions(), plot_rationoosc(), plotContProf(), plotDataPred(), PlotDeltaFractionVsSignificance(), PlotDeltaFractionVsYear(), PlotSignificanceVsYear_FixDelta(), PlotSignificanceVsYear_FixTh23(), plotTime(), calib::CDPStorage::postCloseOutputFile(), calib::CDPStorage::postOpenFile(), ana::PredictionSyst3Flavor2020::PredictionSyst3Flavor2020(), nutools::dbi::Table::PrintPMUsed(), nova::dbi::Table::PrintPMUsed(), nutools::dbi::Table::PrintVMUsed(), nova::dbi::Table::PrintVMUsed(), reach_2018_dCPfractions(), read_file(), read_t2k_rootracker(), daq2raw::TestBeamRawInputDriver::readBeamline(), om::HistoTable::ReadFile(), red(), redND(), resolution2018(), resPlot(), resPlotND(), rock(), om::GUIModelData::SetCSVFile(), test_beam_errorband(), test_fluxhadr_prod_weights_Flux(), test_genieweights(), test_nue2018_surface(), stan::io::program_reader::trace(), train(), train_tmva_cosrej_nonc(), trim_caf(), ValidateRootFile(), mcdatadaq::Configuration::writeToFile(), xsec_extrap_plots(), and zerosstudy().

TH1D* fResidV[pmax][cmax] = {}

Definition at line 19 of file cellShifts.C.

Referenced by initialise(), plot_residPerPlane(), plot_residV(), and plot_sumcellsV().

TH1D* fResidZ[pmax][cmax] = {}

Definition at line 20 of file cellShifts.C.

Referenced by initialise(), plot_residPerPlane(), plot_residZ(), and plot_sumcellsZ().

TH2D* hCellPos[pmax][cmax] = {}

Definition at line 21 of file cellShifts.C.

Referenced by align::Alignment::GetChannelHists(), initialise(), plot_CellPos(), plot_residPerViewV(), and plot_residPerViewZ().

const int nbins = 200

Definition at line 15 of file cellShifts.C.

Referenced by genie::HadronizationModelBase::ApplyRijk(), nueccinc_test::NueCCIncGlobalFitter::CheckSystsSize(), cmf::CMFDecorrelator::CMFDecorrelator(), CombineHistograms(), ana::CompareNDDataMCFromVector(), ana::CompareNDDataOneMC(), ana::CompareNDDataTwoMC(), ana::CompareOneShiftPred(), ana::CompareTwoShiftPred(), genie::HadronizationModelBase::CreateMultProbHist(), datamc_ND_numu_kinematics_FHC_REW(), datamc_ND_numu_kinematics_FHC_REW_pTBins(), datamc_ND_numu_kinematics_REW(), datamc_ND_numu_kinematics_RHC_REW(), datamc_ND_numu_kinematics_RHC_REW_pTBins(), ana::PredictionInterp::DebugPlot(), ana::PredictionInterp::DebugPlotColz(), DiagonalErrors(), DrawBins(), DrawCompBarChart(), DrawCovMx(), DrawPulls(), rsim::ExcessNoiseMaker::DrawSmearedPE(), murem::MuonRemoveAna::endJob(), chaninfo::BadChanValidate::endSubRun(), ana::ErrorBarChart(), FCTutorial2020(), fd_plot(), osc::_OscCalcDMP< T >::FillCache(), FOMCalc(), ana::Hist::FromDirectory(), Gaus(), generate_sk_flux_histograms(), genie_syst_pca(), GetEdges(), GetMedian(), ana::GetNDComponents(), ana::GetNDDecompsFromDecomp(), RandHisto::GetRandom(), GoodDataSelDate(), GoodDataSelDuration(), GoodDataSelEmptySpillFrac(), GoodDataSelGoodDB(), GoodDataSelLiveTime(), GoodDataSelMipRate(), GoodDataSelNumSlices(), GoodDataSelRuns(), GoodDataSelTimingPeakEnd(), GoodDataSelTimingPeakStart(), GoodDataSelTrkFrac(), histCompare(), ana::Hist::Initialized(), levy_func_problem(), LinBins(), ana::FrequentistSurface::LoadFromMulti(), LogBins(), LOverEPlot(), make_nueFDprediction_kinematics_FHC_REW(), make_nueFDprediction_kinematics_FHC_REW_pTExtrap(), make_nueFDprediction_kinematics_REW(), make_nueFDprediction_kinematics_RHC_REW(), make_nueFDprediction_kinematics_RHC_REW_pTExtrap(), make_plots(), MakeLegend(), rsim::ExcessNoiseMaker::MakeTemplate(), MatrixValidation(), mix::MixSimEvents::MixSimEvents(), genie::KNOHadronization::MultiplicityProb(), myBarChart(), Nearline_Watchdog_Plots(), NominalRat(), ana::RecoReweight::OptimalBinningHelper(), plot_multiplicity_prob(), plot_residPerViewV(), plot_residPerViewZ(), PlotEfficiency(), ana::NumuCCIncAnalysis::PlotHelper(), ppfx_syst_pca_fn(), ana::PrintTableHeader(), ana::SystMaker::ProcessSyst(), ResolutionBars(), resolutionplotter(), resolutionplotterbias(), SavePCAForCAFAna(), set_labels_1D(), SetBinsOnOffShift(), SetBinsPlusMinusOne(), osc::_OscCalcDMP< T >::SetCachedEnergies(), SmoothWithFourierFit(), test_mc_rejection_method_lin(), test_mc_rejection_method_log(), cmf::CovarianceTools::TH2ToMat(), TransverseEff(), cmf::CovarianceMatrixMaker::writeResults(), ana::SystMakerShift::WriteTable(), and ana::SystMaker::WriteTable().

TH2D* p1X[pmax] = {}

Definition at line 22 of file cellShifts.C.

Referenced by plot_residPerPlane().

const int pmax = 200

Definition at line 13 of file cellShifts.C.

Referenced by genie::GMCJDriver::ComputeEventProbability(), genie::GMCJDriver::ComputeInteractionProbabilities(), genie::GMCJDriver::ComputeProbScales(), genie::utils::nuclear::FmArea(), initialise(), main(), plot_residPerPlane(), plot_residPerViewV(), plot_residPerViewZ(), plot_sumcellsV(), plot_sumcellsZ(), slicer_production_validation(), and genie::GMCJDriver::~GMCJDriver().