SBNAna/Cuts/TruthCuts.cxx

Go to the documentation of this file.

#include "sbnana/SBNAna/Cuts/TruthCuts.h"
#include "sbnana/SBNAna/Vars/TruthVars.h"
#include "sbnana/SBNAna/Cuts/VolumeDefinitions.h"

#include "sbnanaobj/StandardRecord/Proxy/SRProxy.h"

namespace ana{

  //---------------------------------------------------------------
  // Slice cuts 
  const Cut kTrueActiveVolumeND(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVolAbsX(slc->truth.position, avnd);
    }
    );

  const Cut kTrueFiducialVolumeND(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVolAbsX(slc->truth.position, fvndAbs);
    }
    );

  const Cut kTrueActiveVolumeFDCryo1(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVol(slc->truth.position, avfd_cryo1);
    }
    );

  const Cut kTrueActiveVolumeFDCryo2(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVol(slc->truth.position, avfd_cryo2);
    }
    );

  const Cut kTrueFiducialVolumeFDCryo1(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVol(slc->truth.position, fvfd_cryo1);
    }
    );

  const Cut kTrueFiducialVolumeFDCryo2(
    [](const caf::SRSliceProxy* slc){
      return kHasNu(slc) && PtInVol(slc->truth.position, fvfd_cryo2);
    }
    );

  const Cut kIsAntiNu([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return false;
            return slc->truth.pdg < 0;
  });

  const Cut kIsNu([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return false;
            return slc->truth.pdg > 0;
  });

  const Cut kHasNu([](const caf::SRSliceProxy* slc){
            return slc->truth.index >= 0;
  });

  const Cut kIsNue([](const caf::SRSliceProxy* slc){
            return slc->truth.index >= 0 && abs(slc->truth.pdg) == 12;
  });

  const Cut kIsNumu([](const caf::SRSliceProxy* slc){
            return slc->truth.index >= 0 && abs(slc->truth.pdg) == 14;
  });

  const Cut kIsNutau([](const caf::SRSliceProxy* slc){
            return slc->truth.index >= 0 && abs(slc->truth.pdg) == 16;
  });

  const Cut kIsCC([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return false;
            return (slc->truth.iscc==1);
  });

  const Cut kIsNC([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return false;
            return (slc->truth.isnc==1);
  });

  const Cut kVtxDistMagCut([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return true;
            return (kTruthVtxDistMag(slc) < 1);
  });

  const Cut kSlcCompletenessCut([](const caf::SRSliceProxy* slc){
            if(slc->truth.index < 0) return false;
            return (kCompletness(slc) > 0.5);
  });

  //---------------------------------------------------------------
  // Spill cuts 

  const SpillCut kIsCosmicSpill([](const caf::SRSpillProxy* sr) {
    return (sr->mc.nnu == 0);
  });

  const SpillCut kIsSingleNuSpill([](const caf::SRSpillProxy* sr) {
    return (sr->mc.nnu < 2);
  });

  const SpillCut kIsNueSpill([](const caf::SRSpillProxy* sr){
    if (kIsCosmicSpill(sr) || !kIsSingleNuSpill(sr)) return false;
    return (std::abs(sr->mc.nu[0].pdg) == 12);
  });

  const SpillCut kIsNumuSpill([](const caf::SRSpillProxy* sr){
    if (kIsCosmicSpill(sr) || !kIsSingleNuSpill(sr)) return false;
    return (std::abs(sr->mc.nu[0].pdg) == 14);
  });

  const SpillCut kIsNutauSpill([](const caf::SRSpillProxy* sr){
    if (kIsCosmicSpill(sr) || !kIsSingleNuSpill(sr)) return false;
    return (std::abs(sr->mc.nu[0].pdg) == 16);
  });

  const SpillCut kIsCCSpill([](const caf::SRSpillProxy* sr) {
    if (kIsCosmicSpill(sr) || !kIsSingleNuSpill(sr)) return false;
    return ((bool)sr->mc.nu[0].iscc);
  });

  const SpillCut kIsNCSpill([](const caf::SRSpillProxy* sr) {
    if (kIsCosmicSpill(sr) || !kIsSingleNuSpill(sr)) return false;
    return ((bool)sr->mc.nu[0].isnc);
  });



}