NusPlotLists.cxx

Go to the documentation of this file.

#include "NuXAna/Analysis/NusPlotLists.h"

#include "NuXAna/Vars/NusVars.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "3FlavorAna/Vars/NueVarsExtra.h"

#include "StandardRecord/Proxy/SRProxy.h"

namespace ana
{

  void AddNusHistDefAna01Basic(std::vector<HistDef>& hd)
  {
    hd.push_back({"cvnnc", {"CVN nc", Binning::Simple(100, 0, 1), kCVNnc}});
    hd.push_back({"nhits", {"Number of hits", Binning::Simple(50, 0, 250), kNHit}});
    hd.push_back({"calE", {"Calorimetric energy (GeV)", Binning::Simple(100, 0, 5), kCaloE}});
    hd.push_back({"longprong", {"Longest prong (cm)", Binning::Simple(60, 0, 600), kLongestProng}});
    hd.push_back({"numucontpid", {"BDT score (numucontpid)", Binning::Simple(100, 0, 1), kNumuContPID}});
    hd.push_back({"e_per_hit", {"Calorimetric energy / hit (GeV)", Binning::Simple(50, 0, .05), kEPerHit}});
    hd.push_back({"phot_ptp", {"Photon p_{T}/p", Binning::Simple(50, 0., 1.), kPtpMassless}});
    hd.push_back({"part_ptp", {"Particle p_{T}/p", Binning::Simple(50, 0., 1.), kPartPtp}});
    hd.push_back({"hits_per_plane", {"Hits per plane", Binning::Simple(10, 0, 10), kHitsPerPlane}});
    hd.push_back({"nmichel", {"Number of Michels", Binning::Simple(10, 0, 10), kNMichels}});
  }

  void AddNusHistDefTruth(std::vector<HistDef>& hd)
  {
    hd.push_back({"true_pt",{"True p_{T}",Binning::Simple(50,0,1),kTruetpT}});
    hd.push_back({"true_pz",{"True p_{z}",Binning::Simple(60,0,60),kTruetpz}});
    hd.push_back({"trueE",{"True E",Binning::Simple(100,0,5),kTrueE}});
    hd.push_back({"true_muE",{"True #mu E",Binning::Simple(100,0,5),kTrueMuonE}});
  }

  void AddNusHistDefFDExtra(std::vector<HistDef>& hd)
  {
    hd.push_back({"mindistface", {"Min. Distance to any Detector Face (cm)", Binning::Simple(160, 0., 800.), kDistTop}});
    hd.push_back({"dist_top", {"Distance to Detector Top (cm)", Binning::Simple(160, 0., 800.), kMinDistFace}});
  }

  void AddNusHistDefTrack(std::vector<HistDef>& hd)
  {
    hd.push_back({"trklength_kal", {"Kalman Track Length", Binning::Simple(60, 0, 600), kTrkLength}});
    hd.push_back({"trkEperHit_kal", {"Kalman Track Energy / hit (GeV)", Binning::Simple(50, 0., 0.05), kTrkCalEPerNHit}});
    hd.push_back({"trkNht_kal", {"Number of kalman track hits", Binning::Simple(50, 0., 250.), kTrkNhits}});
  }

  std::vector<HistDef> getAllNusAna01NDHistDefs()
  {
    std::vector<HistDef> ret;
    AddNusHistDefAna01Basic(ret);
    AddNueHistDefVertexND(ret);
    AddNueHistDefShowerND(ret);
    AddNueHistDefShower(ret);
    AddNusHistDefTrack(ret);
    AddNueHistDefPIDs(ret);
    AddNusHistDefTruth(ret);
    return ret;
  }

  std::vector<HistDef> getAllNusAna01FDHistDefs()
  {
    std::vector <HistDef> ret;
    AddNusHistDefAna01Basic(ret);
    AddNusHistDefFDExtra(ret);
    AddNueHistDefVertexFD(ret);
    AddNueHistDefShowerFD(ret);
    AddNueHistDefShower(ret);
    AddNusHistDefTrack(ret);
    AddNueHistDefPIDs(ret);
    return ret;
  }

}//namespace