#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Vars/XsecTunes.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/XSec/Multiverse.h"
#include "CAFAna/XSec/SystematicDef.h"
#include "CAFAna/XSec/GenieMultiverseSyst.h"
#include "CAFAna/Systs/XSecSystLists.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Analysis/Style.h"
#include "CAFAna/Systs/BeamSysts.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TPaletteAxis.h"
#include "THStack.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncCuts.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncVars.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncBins.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncExtra.h"
#include "NDAna/muonid/NDXSecMuonPID.h"
Go to the source code of this file.
Namespaces | |
nuebarccinc | |
Functions | |
void | load_libs_muonid () |
TH2 * | ToTH2 (const TH1 *h, const HistAxis &xaxis, const HistAxis &yaxis) |
TH1 * | MostConservative (TH1 *up, TH1 *down, TH1 *nominal) |
void | PlotMultiverse (Multiverse &mv, TH1 *nominal, std::string title, std::string name, TH1 *error=0) |
void | Plot1D (TH1 *h1, std::string h1_label, std::map< std::string, TH1 * > systs, std::vector< std::string > to_plot, std::string name, double xmax=-1) |
void | Plot2D (TH1 *total, std::map< std::string, TH1 * > systs, std::vector< std::string > to_plot, std::string basename) |
void | abs (TH1 *hist) |
void | MovePalette (TH1 *hist) |
void | wrong_sign_uncertainty (bool make_plots=true, std::string input_file_name="/nova/ana/users/ddoyle/NuebarCCInc/WrongSign/wrong_sign_uncertainty.root", std::string plot_dump="/nova/ana/users/ddoyle/NuebarCCInc/WrongSign/plots") |
Variables | |
std::vector< std::string > | all_systs = {"ppfx", "genie", "calib_shape", "calib_shift", "ll", "cherenkov", "beam_focusing", "stat"} |
std::vector< std::string > | beam_systs = {"ppfx", "beam_focusing"} |
std::map< std::string, nuebarccinc::Style > | syst_styles |
std::map< std::string, double > | xmax |
std::map< string, const ana::BeamSyst * > | beamSystematics |
void abs | ( | TH1 * | hist | ) |
Definition at line 716 of file wrong_sign_uncertainty.C.
References std::abs(), febshutoff_auto::val, submit_syst::x, submit_syst::y, and test::z.
Referenced by fuzz::FuzzyKVertex::AddProng(), AddUpDwUncerts(), AddUpDwUncerts2D(), remid::ReMIdTrain::analyze(), mcchk::ShowerAnaCheck::analyze(), slicer::S4DParamCalc::analyze(), beamlinereco::ToFSingleCounterAnalysis::analyze(), febstat::EventFEBStatus::analyze(), mcchk::CosmicAna::analyze(), novaddt::SNAna::analyze(), murem::MuonRemoveAna::analyze(), mcchk::NeutrinoAna::analyze(), zcl::MMCheater::analyze(), fuzz::FuzzyKValidate::analyze(), bpfit::BPFTmvaTrainer::analyze(), ncs::Xeff::analyze(), slid::LIDTraining::analyze(), ncs::NCAna::analyze(), comi::Leana::AnaNHit(), lem::MakePID::AvgAndFit(), slid::ParticleIDAlg::CalcCellPlaneTransverseDedx(), ana::CalcPrimE(), cheat::BackTracker::CalcTotalEscapingEnergy(), Calculate2DChi2(), trk::KalmanTrackMerge::CanJoinTracks(), osc::_OscCalcDMP< T >::ChannelCacheIdx(), osc::OscCalcPMNSOptEigen::ChannelCacheIdx(), trk::KalmanTrackAna::CheckRecoTracks(), checkCalibrationCSV.NovaCalibTest::checkRunCoverage(), genie::NievesQELCCPXSec::CNCTCLimUcalc(), om::PlotViewer::ContentProjectionTH2F(), gibuu::GiBUURegen::CopyGenieEvent(), genie::NievesQELCCPXSec::deltaLindhard(), novaddt::HoughTrackMaker::do_like_zukai(), plot_xsec_1d::do_single_shift_1bin(), plot_xsec_1d::do_up_down_1bin(), evd::HitTower::Draw(), ana::SystMakerShift::DrawBkg(), ana::SystMaker::DrawBkg(), ana::SystMakerShift::DrawShift(), ana::SystMakerShift::DrawSig(), ana::SystMaker::DrawSig(), evgen::CosmicsGen::Enhance(), beamlineutil::BeamlineMagneticField::FieldFromCurrent(), calib::StopperThreshold::FillHist(), slidt::FillHists(), std_candles.StdCandlesPlotter::FillMichelEPlots(), caf::FillNumuInfo(), lem::FillPotential(), calib::StopperThreshold::FillTree(), calib::EnergyStandardCandles::FillTruthInfo(), murem::MuonRemove::FillTruthInfo(), comi::CanaFilter::filter(), util::TimeFilter::filter(), mono::MonopoleTrack::filter(), novaddt::MichelETrigger::filter(), trk::KalmanTrack::FilterTracker(), caf::FindAndAddMichels(), generic.tools.HistogramTools::findMaxAsymmetry(), common_tools::findMaxAsymmetry(), generic.tools.HistogramTools::findMaxSkew(), common_tools::findMaxSkew(), bsf::BremShowerFilter::findShowerByTruth(), trk::KalmanTrack::FindTracks(), nuesand::FillNueSandbox::FracAngChanges(), lem::FracChargeMatched(), GetAbsolute(), GetAbsolute2D(), GetBinByEdge(), nuesand::FillNueSandbox::GetECF(), gibuu::GiBUURegen::GetEvent(), cvn::GetInteractionType(), cvn::GetInteractionTypeFromSlice(), GetMaxUpDown(), GetMaxUpDown2D(), numusand::NumuSandFxs::getMissingE(), cvn::GetMultiplicityMap(), getPlots(), getPlotsFromNtupleFD(), getPlotsFromNtupleND(), gibuu::GiBUURegen::GetRecordList(), ana::PredictionWriter::GetShiftName(), GetShifts2017(), ana::GetShifts2018(), ana::GetShifts2020(), ana::MINERvA_MEC_Wgt_Var::GetWeight(), mono::Cluster::good_track(), novaddt::smt::Trigger::good_track(), DecTreeLoader::HandleRecord(), ana::hasBPFE(), genie::BaryonResonanceDecayer::HasEvolvedBRs(), ana::MultiHistContainer::Hist(), lem::DistanceMap::InvDist(), genie::BaryonResonanceDecayer::IsDelta(), ana::IsElectronPng(), ana::IsMuonPng(), ana::IsPionPng(), ana::IsProtonPng(), lem::EventSummary::IsSig(), caf::IsTrueDiF(), calib::RockMuonStopperSelection::IsTrulyContainedStopper(), calib::RockMuonStopperSelection::IsTrulyInFiducial(), calib::ThroughgoingSelection::IsTrulyUncontained(), skim::EvaluatorNumuCCpi::KeepEventSlices(), skim::EvaluatorNumuCCpi::KeepSlice(), kinematics(), ana::kMinPngNhit(), ana::kNue2017BasicPartFunc(), load_libs_muonid(), readFlux::Loop(), mono::Cluster::make(), make_dst_cosrejbdttrain(), make_fc_mass_and_oct_nersc_2018(), make_fc_mass_and_oct_nersc_2019(), make_mass_and_oct_fc_input_2019(), cvn::MakeParticleLabel(), airshower::AirSlicer::MakeTrackSlices(), ana::FluxDecomp::MakeWeightsNumuFromKaon(), ana::BENDecomp::MakeWeightsNumuFromKaon(), MakingProfileHadCC(), hough::MultiHough2P::Map(), dt::ViewMerger::MatchScore(), lem::FindLEMMatches::MatchToVertex(), lem::LEM::MatchToVertex(), genie::geometry::ROOTGeomAnalyzer::MaxPathLengthsFluxMethod(), evd::SimulationDrawer::MCTruthTrajectoriesAnyD(), ana::MuonECat(), novaddt::DBSlicer::NeighborScore(), NuMu2020_TrimCAFs(), ana::GetBPFMuonID::operator()(), ana::CCFlavSel::operator()(), ana::GetBestBPFTrack::operator()(), ana::NCFlavSel::operator()(), oscW(), osc::OscCalcDumb::P(), osc::OscCalcPMNS_CPT::P(), osc::OscCalcPMNS_NSI::P(), osc::_OscCalcPMNS< T >::P(), osc::_OscCalcPMNSOpt< T >::P(), osc::OscCalcSterile::P(), osc::OscCalcGeneral::P(), osc::OscCalcPMNSOptEigen::P(), cvn::ParticlesSliceClassify(), PlotSystErrorBand(), ppfx_make_systs(), ppfx_smooth_weights_save(), Print(), plot_xsec_2d::process_ckv_calibshape(), plot_xsec_2d::process_focusing(), plot_xsec_2d::process_light_calibration_mue(), demo.ScoreSD::ProcessHits(), ana::SystMakerShift::ProcessShift(), ana::Prod3AbsCalibLoaders::Prod3AbsCalibLoaders(), ana::Prod3LightLevelLoaders::Prod3LightLevelLoaders(), ana::Prod4AbsCalibLoaders::Prod4AbsCalibLoaders(), ana::Prod4LightLevelLoaders::Prod4LightLevelLoaders(), ana::Prod5AbsCalibLoaders::Prod5AbsCalibLoaders(), ana::Prod5LightLevelLoaders::Prod5LightLevelLoaders(), cheat::ShowerCheater::produce(), cheat::ProngCheater::produce(), slid::Recluster::produce(), numusand::FillSandbox::produce(), calib::FiberCalibration::produce(), lem::AlignLibToVtx::produce(), zcl::SPCluster::produce(), evgen::CosmicPionGen::produce(), beamlinereco::CFDHitReco::produce(), dif::DiFShowerFinder::produce(), evgen::CORSIKAGen::produce(), cheat::MCCheater::produce(), g4rwgt::MakeG4WeightTable::produce(), murem::MuonRemove::produce(), caf::CAFMaker::produce(), evgen::CosmicsGen::ProjectMuonsToDetectorBigBox(), ana::ProngTruthCut(), PlotUtilityNew::pull_NormData(), murem::MuonRemove::RemoveByEfficiency(), murem::MuonRemove::RemoveByProngCVN(), murem::MuonRemove::RemoveByTruth(), resolution_event_tree(), RockFilter_D(), novaddt::smt::Trigger::run_algorithm(), evgen::SingleGen::Sample(), screen_to_activecell(), ana::ReinteractionSyst::SetKnobIndices(), simb::MCTruth::SetNeutrino(), ana::NueExtrapSystBkg2017::Shift(), ana::NueAcceptSystBkg2018FHC::Shift(), ana::NueExtrapSystSignalKin2017::Shift(), ana::NuWROSyst::Shift(), ana::NueAcceptSystSignalKin2018FHC::Shift(), ana::NueAcceptSystSignalKin2020FHC::Shift(), ana::NumuSummedSmallGENIESyst::Shift(), ana::NueAcceptSystBkg2018RHC::Shift(), ana::BeamNueScaleSyst::Shift(), ana::NueAcceptSystSignalKin2018RHC::Shift(), ana::NumuCCScaleSyst::Shift(), ana::NueAcceptSystSignalKin2020RHC::Shift(), ana::NumuTauContaminationSyst::Shift(), signFunc(), ncs::NCAna::SimpleOscProb(), lem::Library::SimpleSurvivalProb(), lem::SimpleSurvivalProb(), murem::MuonRemove::SliceSelection(), murem::MRCCAna::SliceSelection(), cheat::BackTracker::SlicesToMCTruthsTable(), lem::LEMSummarizer::SliceToLEMInput(), osc::_PMNS< T >::SortEigenvalues(), test_flux(), test_fluxhadr_prod_weights_Flux(), test_genieweights(), caf_tools::tightBeamQualityCut(), ana::TopHatPriorInFitVar(), mono::Track3D::Track3D(), cheat::BackTracker::TrackIDToMCTruth(), cheat::BackTracker::TrackIDToMotherParticle(), cheat::BackTracker::TrackIDToParticle(), train_dectree_caf(), TrimCAFs(), trimmubarid(), trimvar(), nuebarccinc::true_electron_png(), slid::NuEEnergyAlg::VertexEnergy(), HandyFuncs::VerticalRange(), ana::VisibleHadE(), and wrong_sign_uncertainty().
void load_libs_muonid | ( | ) |
Definition at line 32 of file wrong_sign_uncertainty.C.
References abs(), exit(), make_syst_table_plots::h, h1, analysePickle::hist, MostConservative(), MovePalette(), Plot1D(), Plot2D(), PlotMultiverse(), runNovaSAM::ret, string, systs, plotROC::title, produceInfoGainMatrix::total, ToTH2(), xaxis, xmax, and yaxis.
Referenced by wrong_sign_uncertainty().
TH1 * MostConservative | ( | TH1 * | up, |
TH1 * | down, | ||
TH1 * | nominal | ||
) |
Definition at line 821 of file wrong_sign_uncertainty.C.
References std::abs(), std::max(), runNovaSAM::ret, ana::UniqueName(), submit_syst::x, submit_syst::y, and test::z.
Referenced by load_libs_muonid(), and wrong_sign_uncertainty().
void MovePalette | ( | TH1 * | hist | ) |
Definition at line 841 of file wrong_sign_uncertainty.C.
Referenced by load_libs_muonid(), Plot2D(), and wrong_sign_uncertainty().
void Plot1D | ( | TH1 * | h1, |
std::string | h1_label, | ||
std::map< std::string, TH1 * > | systs, | ||
std::vector< std::string > | to_plot, | ||
std::string | name, | ||
double | xmax = -1 |
||
) |
Definition at line 790 of file wrong_sign_uncertainty.C.
References plot_validation_datamc::c, MECModelEnuComparisons::leg, std::max(), maxy, and syst_styles.
Referenced by load_libs_muonid(), and wrong_sign_uncertainty().
void Plot2D | ( | TH1 * | total, |
std::map< std::string, TH1 * > | systs, | ||
std::vector< std::string > | to_plot, | ||
std::string | basename | ||
) |
Definition at line 760 of file wrong_sign_uncertainty.C.
References plot_validation_datamc::c, nuebarccinc::costhetabins, nuebarccinc::eelecbins, MovePalette(), ana::ToTH2Helper(), and ana::UniqueName().
Referenced by estimate_energy(), GetHistsFD(), GetHistsND(), load_libs_muonid(), and wrong_sign_uncertainty().
void PlotMultiverse | ( | Multiverse & | mv, |
TH1 * | nominal, | ||
std::string | title, | ||
std::string | name, | ||
TH1 * | error = 0 |
||
) |
Definition at line 728 of file wrong_sign_uncertainty.C.
References std::abs(), plot_validation_datamc::c, Draw(), make_mec_shifts_plots::GetMaximum(), ana::Multiverse::GetUniversesHist(), hists, kBlue, std::max(), maxy, submit_hadd::u, and submit_syst::x.
Referenced by load_libs_muonid(), and wrong_sign_uncertainty().
Referenced by ana::GenieMultiverseSpectra::GenieMultiverseSpectra(), and load_libs_muonid().
void wrong_sign_uncertainty | ( | bool | make_plots = true , |
std::string | input_file_name = "/nova/ana/users/ddoyle/NuebarCCInc/WrongSign/wrong_sign_uncertainty.root" , |
||
std::string | plot_dump = "/nova/ana/users/ddoyle/NuebarCCInc/WrongSign/plots" |
||
) |
Definition at line 97 of file wrong_sign_uncertainty.C.
References abs(), all_systs, file_size_ana::axes, allInOneTrainingPlots::axis, beam_systs, beamSystematics, plot_validation_datamc::c, nuebarccinc::costhetabins, Cut(), update_sam_good_runs_metadata::cuts, den, ana::Multiverse::Divide(), nuebarccinc::eelecbins, genie::utils::style::Format(), ana::getAllXsecSysts_2020(), ana::GetkPPFXFluxUnivWgt(), ana::Multiverse::GetMinusOneSigmaShift(), ana::Multiverse::GetPlusOneSigmaShift(), ana::GenieMultiverseParameters::GetSystShifts(), hists, MECModelEnuComparisons::i, ana::kAna2020RHCPOT, kGreen, nuebarccinc::kMuonIDProd4Cut, ana::kNoShift, nuebarccinc::kNuebarCC, nuebarccinc::kNueCC, ana::kPPFXFluxCVWgt, nuebarccinc::kPreselectionLoose, ana::kStandardSpillCuts, nuebarccinc::kTrueElectronEVsCosStandardAxis, nuebarccinc::kTrueNeutrinoEnergyStandardAxis, nuebarccinc::kTrueQ2StandardAxis, ana::kXSecCVWgt2020, MECModelEnuComparisons::leg, load_libs_muonid(), loaders, ana::LoadFrom(), make_plots(), MostConservative(), MovePalette(), getGoodRuns4SAM::n, num, output, Plot1D(), Plot2D(), PlotMultiverse(), pot, pots, cet::pow(), nuebarccinc::PROD5_MC_RHC_CALIB_DOWN, nuebarccinc::PROD5_MC_RHC_CALIB_SHAPE, nuebarccinc::PROD5_MC_RHC_CALIB_UP, nuebarccinc::PROD5_MC_RHC_CHERENKOV, nuebarccinc::PROD5_MC_RHC_LIGHT_DOWN, nuebarccinc::PROD5_MC_RHC_LIGHT_UP, nuebarccinc::PROD5_MC_RHC_NOMINAL, std::sqrt(), string, ana::ToTH2Helper(), ana::UniqueName(), febshutoff_auto::val, submit_syst::x, xmax, submit_syst::y, and test::z.
std::vector<std::string> all_systs = {"ppfx", "genie", "calib_shape", "calib_shift", "ll", "cherenkov", "beam_focusing", "stat"} |
Definition at line 47 of file wrong_sign_uncertainty.C.
Referenced by make_cm_pullterm_pdf(), MakeCovarSim(), and wrong_sign_uncertainty().
std::vector<std::string> beam_systs = {"ppfx", "beam_focusing"} |
Definition at line 48 of file wrong_sign_uncertainty.C.
Referenced by wrong_sign_uncertainty().
std::map<string, const ana::BeamSyst*> beamSystematics |
Definition at line 83 of file wrong_sign_uncertainty.C.
Referenced by wrong_sign_uncertainty().
std::map<std::string, nuebarccinc::Style> syst_styles |
Definition at line 50 of file wrong_sign_uncertainty.C.
Referenced by Plot1D().
std::map<std::string, double> xmax |
Definition at line 79 of file wrong_sign_uncertainty.C.
Referenced by genie::Spline::Add(), calib::HitEfficiency::analyze(), showere::ShowerEnergyAna::analyze(), evd::TZProjPad::AutoZoom(), evd::TZProjPad::AutoZoomBox(), evd::TZProjPad::AutoZoomHelper(), evd::TZProjPad::AutoZoomSlice(), evd::TZProjPad::AutoZoomTruth(), ana::BayesianSurface::BayesianSurface(), novaddt::DAQHeaderAna::beginJob(), comi::NumiFilteringAna::beginJob(), genie::BLI2DNonUnifGrid::BLI2DNonUnifGrid(), genie::BLI2DUnifGrid::BLI2DUnifGrid(), ana::BayesianSurface::BuildHist(), genie::Spline::BuildSpline(), CalcFiberLoopCorr(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), genie::DISKinematicsGenerator::ComputeMaxXSec(), genie::SKKinematicsGenerator::ComputeMaxXSec(), genie::DMDISKinematicsGenerator::ComputeMaxXSec(), trk::KalmanTrackAna::Containment(), DCSPlotMaker(), demo_CPT(), genie::utils::kinematics::DISImportanceSamplingEnvelope(), genie::Spline::Divide(), Draw1DProfile(), DrawFrame(), DrawMedianAndUniverses(), DrawSliceCanvas(), ana::CovMxSurface::DrawSurfacePoint(), DrawSurfacePoint(), ana::ExtendGraphToTop(), showere::ShowerEnergyFilterMC::filter(), calib::AttenFit::FitQuality(), calib::AttenuationFit::FitQuality(), evd::XZYZProjectionsView::FixAspectRatio(), genie::Spline::GetAsTGraph(), novaddt::MoonShadow::GetHisto(), moonshadowana::MoonShadowAna::GetHisto(), moonshadowana::MoonShadowAnaHough::GetHisto(), GoodDataSelDate(), GoodDataSelDuration(), GoodDataSelEmptySpillFrac(), GoodDataSelGoodDB(), GoodDataSelLiveTime(), GoodDataSelMipRate(), GoodDataSelNumSlices(), GoodDataSelRuns(), GoodDataSelTimingPeakEnd(), GoodDataSelTimingPeakStart(), GoodDataSelTrkFrac(), genie::BLI2DUnifGrid::Init(), genie::BLI2DNonUnifGrid::Init(), remid::ReMIdValidate::IsContained(), remid::RecoMuon::IsTrackContained(), jointsensitivity(), load_libs_muonid(), ana::LogLikelihood(), main(), trk::CosmicTrackAna::Make1DRecoHistograms(), trk::CosmicTrackAna::Make2DRecoHistograms(), MakeBkgEMPlot(), ana::MakeHistCanvasReady_Quant(), MakeLegend(), genie::Spline::Multiply(), PID2DAxis(), PlotEfficiency(), plotTrainTest(), genie::Spline::Print(), PrintPlot(), numusand::FillSandbox::produce(), request_edarea_size(), resolutionplotter(), resolutionplotterbias(), sensitivity2020(), evd::TZProjPad::ShowFull(), evd::TZProjPad::ShowPartial(), SliceLegend(), ana::SurfaceKrige::SurfaceKrige(), test_mc_rejection_method_lin(), test_mc_rejection_method_log(), TrimGraph(), wrong_sign_uncertainty(), and genie::BLI2DGrid::ZMax().