Functions | |
def | GetMaximum (hist) |
def | PrepHist (hist, scale, ytitle, line_style, color) |
Variables | |
int | line_width = 4 |
float | top_margin = 0.08 |
float | bottom_margin = 0.12 |
float | left_margin = 0.15 |
float | right_margin = 0.05 |
string | default_outdir = "./plots/mec_shifts/" |
Parse Options. More... | |
parser = OptionParser() | |
dest | |
default | |
None | |
help | |
action | |
False | |
opts | |
args | |
string | indir = "/pnfs/nova/persistent/users/mislivec/xsec_tuning_paper/" |
hists_file_name_data = os.path.join( indir, "hists_2018_%s_data_full.root" % opts.beam ) | |
hists_file_name_mc = os.path.join( indir, "hists_2018_%s_mc_full.root" % opts.beam ) | |
hists_file_data = TFile( hists_file_name_data ) | |
hists_file_mc = TFile( hists_file_name_mc ) | |
tree_pot_data = hists_file_data.Get( "pot" ) | |
float | pot_data = 0.0 |
tree_pot_mc = hists_file_mc.Get( "pot" ) | |
float | pot_mc = 0.0 |
mc_pot_scale = pot_data/pot_mc | |
int | event_scale_power = 4 |
float | event_scale = 1.0 |
string | text_events = "10^{%d} Events" |
string | isRHC = "RHC" |
plots_dir = opts.outdir | |
float | label_text_size = 0.04 |
int | text_font = 42 |
float | text_angle = 0.0 |
bool | setNDC = True |
float | ypos_top = 1.0 |
string | text_beam = "Neutrino Beam" |
label_beam = root_style.GetLabel( text_beam, left_margin + label_text_size, ypos_top, label_text_size, kBlack, text_font, 13, text_angle, setNDC ) | |
string | text_nova = "NOvA" |
string | text_minerva = "MINERvA" |
label_prelim = root_style.GetLabel( text_nova + " Preliminary", 1.0 - right_margin, 1.0 - top_margin + 0.02, 0.06, kBlue, text_font, 31, text_angle, setNDC ) | |
hist_data = hists_file_data.Get( "%s_Data" % var ) | |
hist_cv_mec = hists_file_mc.Get( "%s_MEC" % var ) | |
hist_cv_tot_mc = hists_file_mc.Get( "%s_TotMC" % var ) | |
hist_cv_non_mec = hist_cv_tot_mc.Clone( "%s_NonMEC" % var ) | |
hist_neg_shift_tot_mc = hists_file_mc.Get( "%s_MECShapeDown_TotMC" % var ) | |
hist_pos_shift_tot_mc = hists_file_mc.Get( "%s_MECShapeUp_TotMC" % var ) | |
hist_minerva_tot_mc = hists_file_mc.Get( "%s_MEC_MINERvA" % var ) | |
mc_scale = event_scale*mc_pot_scale | |
legend = root_style.GetLegend() | |
xmax = hist_data.GetXaxis().GetXmax() | |
float | xmin = 0.0 |
xndiv = hist_data.GetXaxis().GetNdivisions() | |
float | ymax_scale = 1.5 |
ymax = GetMaximum( hist_data ) | |
string | canvas_name = "Canvas_" |
canvas = TCanvas( canvas_name, canvas_name, 300, 300 ) | |
hist_data_ratio = hist_data.Clone( "%s_Data_Ratio" % var ) | |
hist_pos_shift_tot_mc_ratio = hist_pos_shift_tot_mc.Clone( "%s_ShiftsQELike_TotMC_Ratio" % var ) | |
hist_neg_shift_tot_mc_ratio = hist_neg_shift_tot_mc.Clone( "%s_ShiftsRESLike_TotMC_Ratio" % var ) | |
hist_cv_tot_mc_ratio = hist_cv_tot_mc.Clone( "%s_TotMC_Ratio" % var ) | |
hist_minerva_tot_mc_ratio = hist_minerva_tot_mc.Clone( "%s_TotMC_MINERvA_Ratio" % var ) | |
float | lower_fraction = 0.35 |
float | upper_fraction = 1.0 |
float | title_size = 0.05 |
float | label_size = 0.035 |
pad_upper = canvas.GetPad( 1 ) | |
pad_lower = canvas.GetPad( 2 ) | |
string | plot_name = "%s_%s_MEC_Systs" |
def make_mec_shifts_plots.GetMaximum | ( | hist | ) |
Definition at line 20 of file make_mec_shifts_plots.py.
Referenced by calculateWrongSignNue(), calculateWrongSignNumuQ1(), calculateWrongSignNumuQ2(), calculateWrongSignNumuQ3(), calculateWrongSignNumuQ4(), data_over_mc_profile(), DrawSurfacePoint(), lem::EventSummary.FillHists(), FindLimitY(), Landau(), MakeComparisonPlots(), MakeComparisonPlotsWithRatio(), MakeMePlot(), makeNuECanvas(), makeNuMuCanvas(), MakeSplitCans(), mre_comp_split(), muonID_plots(), Plot(), plot_DataMCComp_numu(), plot_nd_spectra_2018(), plot_pid(), plot_predictions(), plot_recoE_numu(), ana.PlotAll(), ana::CutOptimization.PlotDebug(), ana.PlotDebug(), PlotFitResults(), PlotMCComponentsErrorBand(), ana.PlotMCComponentsErrorBand(), PlotMultiverse(), ana.PlotNDDataTotalMCComparison(), PlotNDDataTotalMCComparison(), PlotRat(), Plotting_DataAndPrediction(), PlotWithRatio2D(), and PlotXSec().
def make_mec_shifts_plots.PrepHist | ( | hist, | |
scale, | |||
ytitle, | |||
line_style, | |||
color | |||
) |
Definition at line 24 of file make_mec_shifts_plots.py.
References SetTitle().
make_mec_shifts_plots.action |
Definition at line 46 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.args |
Definition at line 51 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.bottom_margin = 0.12 |
Definition at line 15 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.canvas = TCanvas( canvas_name, canvas_name, 300, 300 ) |
Definition at line 196 of file make_mec_shifts_plots.py.
string make_mec_shifts_plots.canvas_name = "Canvas_" |
Definition at line 195 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.default |
Definition at line 44 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.default_outdir = "./plots/mec_shifts/" |
Parse Options.
Definition at line 41 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.dest |
Definition at line 44 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.event_scale = 1.0 |
Definition at line 88 of file make_mec_shifts_plots.py.
int make_mec_shifts_plots.event_scale_power = 4 |
Definition at line 87 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.False |
Definition at line 46 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.help |
Definition at line 44 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_cv_mec = hists_file_mc.Get( "%s_MEC" % var ) |
Definition at line 125 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_cv_non_mec = hist_cv_tot_mc.Clone( "%s_NonMEC" % var ) |
Definition at line 127 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_cv_tot_mc = hists_file_mc.Get( "%s_TotMC" % var ) |
Definition at line 126 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_cv_tot_mc_ratio = hist_cv_tot_mc.Clone( "%s_TotMC_Ratio" % var ) |
Definition at line 211 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_data = hists_file_data.Get( "%s_Data" % var ) |
Definition at line 123 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_data_ratio = hist_data.Clone( "%s_Data_Ratio" % var ) |
Definition at line 208 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_minerva_tot_mc = hists_file_mc.Get( "%s_MEC_MINERvA" % var ) |
Definition at line 132 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_minerva_tot_mc_ratio = hist_minerva_tot_mc.Clone( "%s_TotMC_MINERvA_Ratio" % var ) |
Definition at line 212 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_neg_shift_tot_mc = hists_file_mc.Get( "%s_MECShapeDown_TotMC" % var ) |
Definition at line 130 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_neg_shift_tot_mc_ratio = hist_neg_shift_tot_mc.Clone( "%s_ShiftsRESLike_TotMC_Ratio" % var ) |
Definition at line 210 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_pos_shift_tot_mc = hists_file_mc.Get( "%s_MECShapeUp_TotMC" % var ) |
Definition at line 131 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hist_pos_shift_tot_mc_ratio = hist_pos_shift_tot_mc.Clone( "%s_ShiftsQELike_TotMC_Ratio" % var ) |
Definition at line 209 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hists_file_data = TFile( hists_file_name_data ) |
Definition at line 66 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hists_file_mc = TFile( hists_file_name_mc ) |
Definition at line 67 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hists_file_name_data = os.path.join( indir, "hists_2018_%s_data_full.root" % opts.beam ) |
Definition at line 59 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.hists_file_name_mc = os.path.join( indir, "hists_2018_%s_mc_full.root" % opts.beam ) |
Definition at line 60 of file make_mec_shifts_plots.py.
Definition at line 57 of file make_mec_shifts_plots.py.
Referenced by ConsolidateSystMaker().
string make_mec_shifts_plots.isRHC = "RHC" |
Definition at line 100 of file make_mec_shifts_plots.py.
Referenced by comi::NearlineAna.analyze(), Beam(), caf::Proxy< caf::SRNeutrino >.CheckEquals(), caf::Proxy< caf::SRSpill >.CheckEquals(), flat::Flat< caf::SRNeutrino >.Clear(), flat::Flat< caf::SRSpill >.Clear(), ana::CompNormSyst.CompNormSyst(), datamc_ND_numu_kinematics_FHC_pTBins(), datamc_ND_numu_kinematics_FHC_REW_pTBins(), datamc_ND_numu_kinematics_RHC_pTBins(), datamc_ND_numu_kinematics_RHC_REW_pTBins(), flat::Flat< caf::SRNeutrino >.Fill(), flat::Flat< caf::SRSpill >.Fill(), FillGraphs(), GetExtrapolationDefs(), ana.GetJointFitSystematicList(), GetJointFitSystematicList(), remid::RecoMuon.Init(), ana.initNumuHadE2018Var(), ana.initNumuHadE2020Var(), ana.initNumuMuE2018Var(), ana.initNumuMuE2020Var(), make_fc_mass_and_oct_nersc_2018(), make_fc_mass_and_oct_nersc_2019(), make_fc_mh_nersc_2018(), make_fc_mh_nersc_2019(), make_fc_oct_nersc_2018(), make_fc_oct_nersc_2019(), make_fc_surfaces_2020(), make_fc_surfaces_2020_validation(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), make_fc_surfaces_nersc_2019(), make_nueFDprediction_kinematics_FHC_pTExtrap(), make_nueFDprediction_kinematics_FHC_REW_pTExtrap(), make_nueFDprediction_kinematics_RHC_pTExtrap(), make_nueFDprediction_kinematics_RHC_REW_pTExtrap(), make_predictions_systs(), mec_nux_tester_2020(), mec_tuning(), mec_tuning_fitter_2020(), mec_tuning_preds_2020(), ND_DataMC(), cmf::ShifterAndWeighter.NormSystWeight(), caf::Proxy< caf::SRNeutrino >.operator=(), caf::Proxy< caf::SRSpill >.operator=(), NumuEnergyFunc.predict_prod4_nd_had_energy_p(), NumuEnergyFunc.predict_prod5_nd_had_energy_p(), util::RewriteSpillData.produce(), cosrej::CosRej.produce(), ndreco::NDReco.produce(), resolution2018(), SliceLID::FillSliceLID.selectModel(), LSTME::FillLSTME.selectModel(), and ana.TAHadEND().
make_mec_shifts_plots.label_beam = root_style.GetLabel( text_beam, left_margin + label_text_size, ypos_top, label_text_size, kBlack, text_font, 13, text_angle, setNDC ) |
Definition at line 114 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.label_prelim = root_style.GetLabel( text_nova + " Preliminary", 1.0 - right_margin, 1.0 - top_margin + 0.02, 0.06, kBlue, text_font, 31, text_angle, setNDC ) |
Definition at line 119 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.label_size = 0.035 |
Definition at line 248 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.label_text_size = 0.04 |
Definition at line 105 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.left_margin = 0.15 |
Definition at line 16 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.legend = root_style.GetLegend() |
Definition at line 148 of file make_mec_shifts_plots.py.
Referenced by AnaPlotMaker(), calibrationconstants(), compare_fits(), ana::ModularExtrapComponent.ComparisonPlot(), ContourComparisonPlot(), corrected_PE(), CreateLegend(), Plotter.DataMCComparison(), DrawLegendBF(), DrawLegendBF2(), DrawLegendBF3(), DrawLegendBFNull(), example_plot(), main(), mlpPIDelecallEetot(), mlpPIDelecallEtot(), OnMonAnaHistos(), plot_contprof(), PandAna.reco_validation.prod5_pid_validation_data_mc.plot_data_mc(), plot_datapredictions(), plot_ehade_quantbound(), plot_predictions(), plot_quantbound_twosamples(), plot_quantile_boundaries_2020(), plot_quantile_boundaries_twosamples_2020(), PandAna.reco_validation.prod5_pid_validation.plot_ratio(), plot_rationoosc(), plotContProf(), plotDataPred(), plotDataPred_SplitCanvas(), TMVAEIDTraining.runDiagnostics(), SaveToPsFile(), tbrecoplots(), thresholdshadowing(), Plotter.TruthPlots(), uncorrected_PE(), and Plotter.WeightPlots().
int make_mec_shifts_plots.line_width = 4 |
Definition at line 13 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.lower_fraction = 0.35 |
Definition at line 244 of file make_mec_shifts_plots.py.
Definition at line 81 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.mc_scale = event_scale*mc_pot_scale |
Definition at line 137 of file make_mec_shifts_plots.py.
Referenced by llike_hist().
make_mec_shifts_plots.None |
Definition at line 44 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.opts |
Definition at line 51 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.pad_lower = canvas.GetPad( 2 ) |
Definition at line 264 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.pad_upper = canvas.GetPad( 1 ) |
Definition at line 250 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.parser = OptionParser() |
Definition at line 43 of file make_mec_shifts_plots.py.
string make_mec_shifts_plots.plot_name = "%s_%s_MEC_Systs" |
Definition at line 292 of file make_mec_shifts_plots.py.
Referenced by DrawSurface(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_contours(), joint_fit_2019_slices(), MatrixValidation(), plot_joint_fit_2020_contours(), plot_joint_fit_2020_slices(), and xsec_extrap_plots().
make_mec_shifts_plots.plots_dir = opts.outdir |
Definition at line 102 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.pot_data = 0.0 |
Definition at line 70 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.pot_mc = 0.0 |
Definition at line 76 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.right_margin = 0.05 |
Definition at line 17 of file make_mec_shifts_plots.py.
bool make_mec_shifts_plots.setNDC = True |
Definition at line 108 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.text_angle = 0.0 |
Definition at line 107 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.text_beam = "Neutrino Beam" |
Definition at line 112 of file make_mec_shifts_plots.py.
string make_mec_shifts_plots.text_events = "10^{%d} Events" |
Definition at line 89 of file make_mec_shifts_plots.py.
int make_mec_shifts_plots.text_font = 42 |
Definition at line 106 of file make_mec_shifts_plots.py.
string make_mec_shifts_plots.text_minerva = "MINERvA" |
Definition at line 117 of file make_mec_shifts_plots.py.
string make_mec_shifts_plots.text_nova = "NOvA" |
Definition at line 116 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.title_size = 0.05 |
Definition at line 247 of file make_mec_shifts_plots.py.
Referenced by ana.Nue2017FourBinAxis(), ana.Nue2018ThreeBinAxis(), and Nue2018ThreeBinAxis_Reduced().
float make_mec_shifts_plots.top_margin = 0.08 |
Definition at line 14 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.tree_pot_data = hists_file_data.Get( "pot" ) |
Definition at line 69 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.tree_pot_mc = hists_file_mc.Get( "pot" ) |
Definition at line 75 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.upper_fraction = 1.0 |
Definition at line 245 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.xmax = hist_data.GetXaxis().GetXmax() |
Definition at line 163 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.xmin = 0.0 |
Definition at line 164 of file make_mec_shifts_plots.py.
Referenced by genie::Spline.Add(), showere::ShowerEnergyAna.analyze(), evd::TZProjPad.AutoZoom(), evd::TZProjPad.AutoZoomBox(), evd::TZProjPad.AutoZoomSlice(), evd::TZProjPad.AutoZoomTruth(), ana::BayesianSurface.BayesianSurface(), 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(), demo_CPT(), genie::utils::kinematics.DISImportanceSamplingEnvelope(), genie::Spline.Divide(), Draw1DProfile(), draw_background(), DrawFrame(), DrawMedianAndUniverses(), DrawSliceCanvas(), DrawSurfacePoint(), ana.ExpandedHistogram(), ana.ExtendGraphToTop(), showere::ShowerEnergyFilterMC.filter(), FindSignPoint(), ana.FindValley(), calib::AttenFit.FitQuality(), calib::AttenuationFit.FitQuality(), evd::XZYZProjectionsView.FixAspectRatio(), genie::Spline.GetAsTGraph(), evd::SliceNavigator.GetBox(), GoodDataSelDate(), GoodDataSelDuration(), GoodDataSelEmptySpillFrac(), GoodDataSelGoodDB(), GoodDataSelLiveTime(), GoodDataSelMipRate(), GoodDataSelNumSlices(), GoodDataSelRuns(), GoodDataSelTimingPeakEnd(), GoodDataSelTimingPeakStart(), GoodDataSelTrkFrac(), genie::BLI2DUnifGrid.Init(), genie::BLI2DNonUnifGrid.Init(), remid::ReMIdValidate.IsContained(), remid::RecoMuon.IsTrackContained(), jointsensitivity(), ana.LogLikelihood(), main(), ana.MakeHistCanvasReady_Quant(), MakeLegend(), genie::Spline.Multiply(), ReadBDTG.NormVariable(), ReadBDTA.NormVariable(), PID2DAxis(), PlotEfficiency(), plotTrainTest(), genie::Spline.Print(), PrintPlot(), numusand::FillSandbox.produce(), resolutionplotter(), resolutionplotterbias(), sensitivity2020(), evd::TZProjPad.ShowFull(), evd::TZProjPad.ShowPartial(), SliceLegend(), ana::SurfaceKrige.SurfaceKrige(), test_mc_rejection_method_lin(), test_mc_rejection_method_log(), TrimGraph(), and genie::BLI2DGrid.ZMax().
int make_mec_shifts_plots.xndiv = hist_data.GetXaxis().GetNdivisions() |
Definition at line 165 of file make_mec_shifts_plots.py.
make_mec_shifts_plots.ymax = GetMaximum( hist_data ) |
Definition at line 186 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.ymax_scale = 1.5 |
Definition at line 166 of file make_mec_shifts_plots.py.
float make_mec_shifts_plots.ypos_top = 1.0 |
Definition at line 110 of file make_mec_shifts_plots.py.