Namespaces | Macros | Functions | Variables
NuSLoadMacroDefs.h File Reference
#include <stdio.h>
#include <string>
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Core/Cut.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Cuts/Cuts.h"
#include "NuXAna/Cuts/NusCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Vars/Vars.h"
#include "OscLib/IOscCalc.h"
#include "TVector3.h"

Go to the source code of this file.

Namespaces

 ana
 Cuts and Vars for the 2020 FD DiF Study.
 

Macros

#define MAXSPEC   5
 
#define SPECARR(label, bins, loaderbeam, loadercos, var, cut)   {Spectrum(label, bins, loaderbeam, var, cut && FlavSel[0]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[1]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[2]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[3]), Spectrum(label, bins, loadercos, var, cut)}
 
#define SPECARR2D(label, loaderbeam, loadercos, binsX, varX, binsY, varY, cut)   {Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[0]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[1]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[2]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[3]), Spectrum(label, loadercos, binsX, varX, binsY, varY, cut)}
 

Functions

void Preliminary ()
 Put NOvA Preliminary on plots. More...
 
void Simulation ()
 
void CenterTitles (TH1 *histo)
 
std::vector< std::stringana::MakeCosOutOfTimeList ()
 
const std::vector< std::stringana::fnameblind_concat (MakeCosOutOfTimeList())
 
const Cut ana::kNDTruthFull ([](const caf::SRProxy *sr){if(sr->hdr.ismc==false){return true;}if(sr->mc.nnu< 1){return false;}if(sr->mc.nu[0].vtx.X()< NDL){return false;}if(sr->mc.nu[0].vtx.X() > NDR){return false;}if(sr->mc.nu[0].vtx.Y()< NDB){return false;}if(sr->mc.nu[0].vtx.Y() > NDT){return false;}if(sr->mc.nu[0].vtx.Z()< NDF){return false;}double NDM=1587.;if(sr->mc.nu[0].vtx.Z() > NDM){return false;}if(sr->mc.nu[0].vtx.Z() > NDE){double MCT=(2./3.)*(NDT-NDB)+NDB;if(sr->mc.nu[0].vtx.Y() > MCT){return false;}}return true;})
 

Variables

const std::string ana::diskdir = "/pnfs/nova/persistent/users/gsdavies/nus/concat/S16-04-08/"
 
const std::string ana::diskdirFD = "/pnfs/nova/persistent/users/gsdavies/nus/concat/S16-05-20/"
 
const std::string ana::diskdir_fcs = "/pnfs/nova/persistent/production/concat/R16-11-12-feature_caf_size.a/"
 
const std::string ana::fnamenear = "defname: prod_decaf_S16-04-08_nd_genie_nonswap_genierw_fhc_nova_v08_full_nus_or_numu_contain_v1 with stride 2 with limit 10000"
 
const std::string ana::fnamefar = "prod_decaf_S16-04-08_fd_genie_nonswap_fhc_nova_v08_full_nus_contain_v1"
 
const std::string ana::fnameswap = "prod_decaf_S16-04-08_fd_genie_fluxswap_fhc_nova_v08_full_nus_contain_v1"
 
const std::string ana::fnametau = "prod_decaf_S16-04-08_fd_genie_tau_fhc_nova_v08_full_nus_contain_v1"
 
const std::string ana::fnamenearcaf = "defname: prod_caf_R16-03-03-prod2reco.a_nd_genie_nonswap_genierw_fhc_nova_v08_full_v1 with limit 10000"
 
const std::string ana::fnamefarcaf = "prod_caf_R16-03-03-prod2reco.a_fd_genie_nonswap_fhc_nova_v08_full_v1"
 
const std::string ana::fnameswapcaf = "prod_caf_R16-03-03-prod2reco.a_fd_genie_fluxswap_fhc_nova_v08_full_v1"
 
const std::string ana::fnametaucaf = "prod_caf_R16-03-03-prod2reco.a_fd_genie_tau_fhc_nova_v08_full_v1"
 
const std::string ana::fnamenear_concat = diskdir + "prod_decaf_S16-04-08_nd_genie_nonswap_genierw_fhc_nova_v08_full_nus_or_numu_contain_v1/*.root"
 
const std::string ana::fnamenear_concat_fcs = diskdir_fcs + "prod_decaf_R16-11-12-feature_decaf_size.a_nd_genie_nonswap_genierw_fhc_nova_v08_full_v1_nue_or_numu/*.root"
 
const std::vector< std::stringana::fFDMC_non
 
const std::vector< std::stringana::fFDMC_swp
 
const std::vector< std::stringana::fFDMC_tau
 
const std::string ana::fFDMC_non_3c = diskdirFD + "by_period/prod_decaf_S16-05-20_fd_genie_nonswap_fhc_nova_v08_epoch3c_nus_contain_v1_prod2-snapshot.root"
 
const std::string ana::fFDMC_swp_3c = diskdirFD + "by_period/prod_decaf_S16-05-20_fd_genie_fluxswap_fhc_nova_v08_epoch3c_nus_contain_v1_prod2-snapshot.root"
 
const std::string ana::fFDMC_tau_3c = diskdirFD + "by_period/prod_decaf_S16-05-20_fd_genie_tau_fhc_nova_v08_epoch3c_nus_contain_v1_prod2-snapshot.root"
 
const std::string ana::fFDMC_non_fcs = "defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_nonswap_fhc_nova_v08_period1_v1 OR defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_nonswap_fhc_nova_v08_period2_v1 OR defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_nonswap_fhc_nova_v08_epoch3b_v1"
 
const std::string ana::fFDMC_swp_fcs = "defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_fluxswap_fhc_nova_v08_period1_v1 OR defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_fluxswap_fhc_nova_v08_period2_v1 OR defname: bckhouse_newcaf_overlap_prod_caf_R16-11-12-feature_caf_size.a_fd_genie_fluxswap_fhc_nova_v08_epoch3b_v1"
 
const std::string ana::fFDMC_non_3c_fcs = "prod_caf_R16-11-12-feature_caf_size.a_fd_genie_nonswap_fhc_nova_v08_epoch3c_v1"
 
const std::string ana::fFDMC_swp_3c_fcs = "prod_caf_R16-11-12-feature_caf_size.a_fd_genie_fluxswap_fhc_nova_v08_epoch3c_v1"
 
const std::string ana::fnamefar_concat = diskdir + "prod_decaf_S16-04-08_fd_genie_nonswap_fhc_nova_v08_full_nus_contain_v1.root"
 
const std::string ana::fnameswap_concat = diskdir + "prod_decaf_S16-04-08_fd_genie_fluxswap_fhc_nova_v08_full_nus_contain_v1.root"
 
const std::string ana::fnametau_concat = diskdir + "prod_decaf_S16-04-08_fd_genie_tau_fhc_nova_v08_full_nus_contain_v1.root"
 
const std::string ana::fnamecos_concat = diskdir + "nus_decaf_S16-04-08_fd_data_cosmic/*.root"
 
const std::string ana::fnamecos_fcs = "prod_caf_R16-11-12-feature_caf_size.a_fd_cosmic_full_nueveto_v1_goodruns"
 
const std::string ana::fnamecosnueveto = "defname:nus_decaf_S16-04-08_fd_data_cosmic with limit 10000"
 
const std::string ana::fnameneardata_concat = diskdir + "prod_decaf_S16-04-08_nd_numi_fhc_full_nus_or_numu_contain_v1_goodruns/*.root"
 
const std::string ana::fnameneardata_concat_fcs = diskdir_fcs + "prod_decaf_R16-11-12-feature_decaf_size.a_nd_numi_fhc_full_v1_goodruns_nue_or_numu/*.root"
 
const std::string ana::fnamefardata = ""
 
const std::string ana::fnamenearsmall = "defname: prod_decaf_S16-04-08_nd_genie_nonswap_genierw_fhc_nova_v08_full_nus_or_numu_contain_v1 with stride 2 with limit 100"
 
const std::string ana::fnamefarsmall = "defname: prod_decaf_S16-04-08_fd_genie_nonswap_fhc_nova_v08_full_nus_contain_v1 with limit 100"
 
const std::string ana::fnameswapsmall = "defname: prod_decaf_S16-04-08_fd_genie_fluxswap_fhc_nova_v08_full_nus_contain_v1 with limit 100"
 
const std::string ana::fnametausmall = "defname: prod_decaf_S16-04-08_fd_genie_tau_fhc_nova_v08_full_nus_contain_v1 with limit 100"
 
const std::string ana::fnamecossmall = "defname: nus_decaf_S16-04-08_fd_data_cosmic with limit 100"
 
const std::string ana::fnameneardatasmall = "defname: prod_decaf_S16-04-08_nd_genie_nonswap_genierw_fhc_nova_v08_full_nus_contain_v1 with stride 2 with offset 1 with limit 100"
 

Macro Definition Documentation

#define MAXSPEC   5

Definition at line 19 of file NuSLoadMacroDefs.h.

#define SPECARR (   label,
  bins,
  loaderbeam,
  loadercos,
  var,
  cut 
)    {Spectrum(label, bins, loaderbeam, var, cut && FlavSel[0]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[1]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[2]), Spectrum(label, bins, loaderbeam, var, cut && FlavSel[3]), Spectrum(label, bins, loadercos, var, cut)}

Definition at line 174 of file NuSLoadMacroDefs.h.

#define SPECARR2D (   label,
  loaderbeam,
  loadercos,
  binsX,
  varX,
  binsY,
  varY,
  cut 
)    {Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[0]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[1]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[2]), Spectrum(label, loaderbeam, binsX, varX, binsY, varY, cut && FlavSel[3]), Spectrum(label, loadercos, binsX, varX, binsY, varY, cut)}

Definition at line 176 of file NuSLoadMacroDefs.h.

Function Documentation

void CenterTitles ( TH1 *  histo)

Definition at line 184 of file nue_nd_datastability.C.

Referenced by accum_nue_numu_equivalent(), accum_nue_numu_equivalent_fhc(), accum_nue_numu_equivalent_rhc(), DataMCPair::AestheticsArea(), DataMCPair::AestheticsPOT(), calcAlphaBetaEachBin(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), CenterTitles(), CompareDecompDataMC(), ana::CompareNDDataMC(), ana::ComparePredictions(), cosnumi_plot(), ana::DefaultFormatNue(), DrawContourCanvas(), DrawSliceCanvas(), FD_Data_PosComp(), febrates(), febshutoff(), fiducial_optimization_plots_2d(), FixMCHist(), FourBinAxis(), GetSpectToHist(), joint_fit_2017_contours(), main(), ana::Make2DFrom1D(), Make2DPlot(), make_plots(), MakePlot(), MakePlots(), MakeRatio(), MakeRawPlot(), MakeResultPlot(), MakeSplitCans(), MatrixElementContours(), monoprob(), MuonFit(), nhit_plot(), nova_official_data_release(), NusContoursPlusMINOS(), DataMCPair::OverlayDataMCSyst(), DataMCPair::OverlayDataMCSystNorm(), PID2DAxis(), Plot1DSlice(), plot_beam_peak(), plot_datapredictions(), plot_oscprob(), plot_pi0_xcheck(), plot_predictions(), plot_rationoosc(), plot_time(), ana::NumuCCIncAnalysis::PlotHelper(), PlotMCComponentsErrorBand(), ana::PlotMultiSyst(), plotNumuWrongSign(), PlotPionPlots(), PlotPurEff(), PlotStack(), PlotSyst(), ana::PlotSyst(), ana::PlotSystBand(), PlotSysts(), PlotSystShifts(), ana::PlotSystShifts(), Plotting_DataAndPrediction(), PlotVertices(), plotWrongSignFraction(), PositionComparison(), print_tables(), ProfileAxes(), ana::RatioPlot(), ana::NumuCCIncAnalysis::Result(), ana::BENDecomp::SavePlots(), ana::BENDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsPi(), SetHBarOptions(), ana::SetHistOptions(), Spec2DToHist(), Spec2DtoHist(), SpecToHist(), tbrecoplots(), time_wide_plot(), time_zoom_plot(), x1_plot(), y1_plot(), and z1_plot().

185 {
186  h->GetXaxis()->CenterTitle();
187  h->GetYaxis()->CenterTitle();
188 }
void Preliminary ( )

Put NOvA Preliminary on plots.

Definition at line 55 of file draw_decomp_plots.C.

References caf::Proxy< caf::SRVertexBranch >::elastic, caf::Proxy< caf::SRElastic >::fuzzyk, caf::Proxy< caf::SRElastic >::IsValid, kBlue, kNueShwContain, kNueVtxContain, plot_validation_datamc_2018::New(), caf::Proxy< caf::SRFuzzyK >::nshwlid, caf::Proxy< caf::SRFuzzyK >::png, prelim, sr, caf::Proxy< caf::SRElastic >::vtx, caf::Proxy< caf::StandardRecord >::vtx, caf::Proxy< caf::SRVector3D >::X(), caf::Proxy< caf::SRVector3D >::Y(), and caf::Proxy< caf::SRVector3D >::Z().

Referenced by accum_nue_numu_equivalent(), accum_nue_numu_equivalent_fhc(), accum_nue_numu_equivalent_rhc(), bicount(), calcAlphaBetaEachBin(), calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), CompareDecompDataMC(), CompareSlices(), cosnumi_plot(), demo6(), demoFitContours(), Draw1DProfile(), DrawCovMx(), FD_Data_PosComp(), FD_DataOnly_PosComp(), hadEFrac_nd_data_mc_systs(), joint_fit_2018_contours(), joint_fit_2019_contours(), main(), make_angle_plot(), make_plots(), MakeCanvasForData(), MakeCanvasForSim(), MakeHBar(), MakeNuEnergyPlots(), MakeRawPlot(), MakeResultPlot(), MakeSplitCans(), makeSystTable(), MatrixValidation(), median_contours(), monoprob(), NCPaperPlot(), ND_DataMC_Comp_Systs(), nhit_plot(), nova_t2k_comparison(), Nus17vsAn01_NuFACT_2D(), NusContoursPlusMINOS(), pi0_xcheck(), plot(), Plot1DSlice(), Plot2DSlice(), plot_beam_peak(), plot_DataMCComp_numu(), plot_joint_fit_2020_contours(), plot_ND_DataMC(), plot_ND_DataMC_energybinning(), plot_nd_spectra_2018(), plot_pi0_xcheck(), plot_pid(), plot_recoE_numu(), plot_time(), plotContProf(), plotDataPred(), plotmodulewise(), plotNumuWrongSign(), PlotNZ(), ana::PlotStack(), PlotSyst(), PlotSysts(), plotWrongSignFraction(), PlotXY(), PositionComparison(), print_tables(), ana::BENDecomp::SavePlotsKa(), ana::BENDecomp::SavePlotsPi(), selection_story_plots(), sensitivity_plot(), starPlot(), test_beam_errorband(), time_wide_plot(), time_zoom_plot(), Tutorial2019FitContours(), x1_plot(), y1_plot(), and z1_plot().

57 {
58  TLatex* prelim = new TLatex(.9, .95, "NOvA Preliminary");
59  prelim->SetTextColor(kBlue);
60  prelim->SetNDC();
61  prelim->SetTextSize(2/30.);
62  prelim->SetTextAlign(32);
63  prelim->Draw();
64 }
TLatex * prelim
Definition: Xsec_final.C:133
enum BeamMode kBlue
void Simulation ( )

Definition at line 16 of file tools.h.

17 {
18  TLatex* prelim = new TLatex(.9, .95, "NOvA Simulation");
19  prelim->SetTextColor(kGray+1);
20  prelim->SetNDC();
21  prelim->SetTextSize(2/30.);
22  prelim->SetTextAlign(32);
23  prelim->Draw();
24 }
TLatex * prelim
Definition: Xsec_final.C:133