NuebarCCIncBins.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/Binning.h"
#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Core/HistAxis.h"

namespace nuebarccinc
{
  ////////////////////////////////////////////////////////////////////
  //Analysis Binning
  ////////////////////////////////////////////////////////////////////
  const std::vector<double> enuedges{0,1.0,1.50,1.8,2.1,2.4,
      2.7,3.0,3.5,4,4.5,5,5.5,6,120};
  const ana::Binning enubins   = ana::Binning::Custom(enuedges);

  //Electron Energy Bins
  const std::vector<double> eelecedges{0,1,1.4,1.65,2.0,2.5,3.0,3.5,
      4.1,4.7,6.0,120};
  const ana::Binning eelecbins = ana::Binning::Custom(eelecedges);

  //Cos Theta Bins
  const std::vector<double> angedges{-1,0.50,0.85,0.90,0.94,0.97,1.001};
  const ana::Binning costhetabins = ana::Binning::Custom(angedges);

  //Q2 Bins
  const std::vector<double> q2edges{0,0.075,0.15,0.225,0.3,0.4,0.5,0.75,
      1.00,2.0,120};
  const ana::Binning q2bins = ana::Binning::Custom(q2edges);

  const std::vector<double> pidedges{0,0.075,0.15,0.225,0.3,0.375,0.45,0.525,0.6,0.7,0.725,0.75,0.775,0.8,0.825,0.85,0.875,0.9,0.925,0.95,0.975,1.0};
  const ana::Binning pidbins = ana::Binning::Custom(pidedges);

  // electron ID
  const std::vector<double> electronid_edges = 
    {-0.6,-0.3,-0.2,-0.13,-0.07,0,0.07,0.14,0.21,0.28,0.35,0.42,0.6};
  const ana::Binning pidbins_electronid = ana::Binning::Custom(electronid_edges);

  //Custom Binning For Electron E Vs CosTheta, 2D
  const ana::Binning double_diff_bins = ana::Binning::Simple(eelecbins.NBins()*
                                                             costhetabins.NBins(),
                                                             0,
                                                             eelecbins.NBins()*
                                                             costhetabins.NBins());
  
  //Custom Binning For Neutrino Energy vs Electron E Vs CosTheta , 3D
  const ana::Binning three_d_bins = ana::Binning::Simple(eelecbins.NBins()*
                                                         costhetabins.NBins()*
                                                         enubins.NBins(),
                                                         0,
                                                         eelecbins.NBins()*
                                                         costhetabins.NBins()*
                                                         enubins.NBins());

}