CutFlow.h

Go to the documentation of this file.

// this holds the cut flow for both ND and FD selections
// the cuts are split into analysis cuts (FDCuts, NDCuts)
// and subspectra cuts (FDSubSpectra, NDSubSpectra), where 
// the subspectra define the total, signal, and out-of-det
// spectra

#pragma once

#include "CAFAna/Core/Cut.h"
#include "NuXAna/Cuts/NusCuts.h"
#include "NuXAna/macros/Nus20/Selection/Structs.h"

using namespace ana;

const ana::Cut kCCNC([](const caf::SRProxy* sr) { 
    if (sr->mc.nnu == 0) return false;          
    return (sr->mc.nu[0].iscc) ? false : true;
    });                                       

const ana::Cut kNDTruthSelection = kCCNC && !kNus20NDTrueVtxOutOfDet;
const ana::Cut kFDTruthSelection = kCCNC && !kNus20FDTrueVtxOutOfDet;

std::vector<CutContainer> cutsFD {
  CutContainer(kNus20FDPreSelection, "Presel"),
  CutContainer(kNus20FDPreSelection && kNus20QualityCuts, "PreselQuality"),
  CutContainer(kNus20FDPreSelection && kNus20QualityCuts && kNus20FDMinSliceDistCut, "PreselQualityFiducial"),
        CutContainer(kNus20FDPreSelection && kNus20QualityCuts && kNus20FDMinSliceDistCut && kNus20FDSelectionCuts, "PreselQualityFiducialSelection"),
        CutContainer(kNus20FDCuts,     "FullSelection"),
};


std::vector<CutContainer> cutsND {
  CutContainer(kNus20NDPreSelection, "Presel"),
        CutContainer(kNus20NDPreSelection && kNus20QualityCuts, "PreselQuality"),
        CutContainer(kNus20NDPreSelection && kNus20QualityCuts && kNus20NDVtxPosCut, "PreselQualityVtx"),
        CutContainer(kNus20NDPreSelection && kNus20QualityCuts && kNus20NDFiducialCuts, "PreselQualityFiducial"),
        CutContainer(kNus20NDPreSelection && kNus20NDCuts, "FullSelection")
};

std::vector<CutContainer> subCutsFD {
        CutContainer(kNoCut, "Total"),
        CutContainer(kFDTruthSelection, "Signal", false),
  //CutContainer(kNus20FDTrueVtxOutOfDet, "TrueVtxOutOfDet", false)
};

std::vector<CutContainer> subCutsND {
        CutContainer(kNoCut, "Total"),
        CutContainer(kNDTruthSelection, "Signal", false),
  //CutContainer(kNus20NDTrueVtxOutOfDet, "TrueVtxOutOfDet", false)
};