#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::string > | ana::MakeCosOutOfTimeList () |
const std::vector< std::string > | ana::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::string > | ana::fFDMC_non |
const std::vector< std::string > | ana::fFDMC_swp |
const std::vector< std::string > | ana::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" |
#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.
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().
void Preliminary | ( | ) |
Put NOvA Preliminary on plots.
Definition at line 55 of file draw_decomp_plots.C.
References ana::assert(), caf::Proxy< caf::SRVertexBranch >::elastic, stan::math::fabs(), caf::Proxy< caf::SRElastic >::fuzzyk, caf::Proxy< caf::SRElastic >::IsValid, kBlue, ana::kIsDytmanMEC, caf::kMEC, kNueShwContain, kNueVtxContain, caf::StandardRecord::mc, plot_validation_datamc_2018::New(), caf::SRTruthBranch::nnu, caf::Proxy< caf::SRFuzzyK >::nshwlid, caf::SRTruthBranch::nu, 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().