TH2D Systematics

[jtenavidal]

ROOT treats them as uncorrelated but they are not

Option 1. Add New function
void plotting::PlotProjectionWithRatio( const std::vector<TH2D*>& mcHists, const std::vector<TH2D*>& breakdown, const TH2D* data_raw, const std::vector<TH2D*>& acceptance, const TH2D* radcorr, const std::string& xobservable, const std::string& yobservable, TPad* topPad, TPad* bottomPad, bool logScale, const std::string& axis, double y_cut_min, double y_cut_max )

Notice the data is before corrections!

We have to first, integrate, then compute the sys. for each bin.

In this space, the uncertainties are uncorrelated

[jtenavidal]

Uncorrelated data available in Plot2DXSec

hist_data_uncorr

And also per sector
need to add a call after

          input_data_location, output_location, output_file_name + "_with_breakdown", systematic_map, true, analysis_id, store_root, log_scale);```
For the slices

plotting::PlotTotalSlicedXSec(mc_hists_xsec, breakdown_xsec, hist_data, acceptance, radcorr, x_observable, y_observable, y_cuts, title, data_name, model, input_MC_location,
          input_data_location, output_location, output_file_name + "_with_breakdown", systematic_map, true, analysis_id, store_root, log_scale)


It is more moving this out of ```PlotTotal2DXSec``` So it is completely decoupled.

```      std::string output_name_2 = output_file_name + "_dxsec_d" + x_observable + "_vs_" + y_observable +"_projectionX";
      plotting::PlotSlicesGeneralized(mc_hists, breakdown, data, x_observable, y_observable, y_cuts, title, output_location, output_name_2, store_root, log_scale );

[jtenavidal]

Modify PlotSlicesGeneralized, PlotProjectionWithRatio and PlotProjectionsStack to add acceptance and rad. corr.