NusCuts.cxx

Go to the documentation of this file.

#include "NuXAna/Cuts/NusCuts.h"

#include "StandardRecord/Proxy/SRProxy.h"

#include <cassert>

namespace ana
{
  const double NDL = -191.;
  const double NDR =  192.;
  const double NDB = -187.;
  const double NDT =  194.;
  const double NDF =    0.;
  const double NDE = 1270.;

  const double FDL = -758.;
  const double FDR =  765.;
  const double FDB = -749.;
  const double FDT =  765.;
  const double FDF =    0.;
  const double FDE = 5962.;

  const Cut kCosmicVeto([](const caf::SRProxy* sr)
                        {
                          return sr->sel.veto.keep;
                        });

  const Cut kNusEventQuality(
    [](const caf::SRProxy* sr)
    {
      if(!sr->vtx.elastic.IsValid)                  return false;
      if(sr->sel.nuecosrej.hitsperplane >= 8)    return false;
      if(sr->vtx.elastic.fuzzyk.nshwlid == 0) return false;
      if(sr->vtx.elastic.fuzzyk.png[0].shwlid.gap >= 100.)  return false;
      if(sr->slc.ncontplanes <= 2)               return false;
      return true;
    });

  const Cut kNusFDFiducial(
    [](const caf::SRProxy* sr)
    {
      // can't use the assertion for N-1 plots
      //assert(sr->vtx.elastic.IsValid && "Must apply DQ cuts");
      if(!sr->vtx.elastic.IsValid) return false;

      const TVector3 vtx = sr->vtx.elastic.vtx;
      if(vtx.X() < -680.0) return false;
      if(vtx.X() >  650.0) return false;
      if(vtx.Y() < -720.0) return false;
      if(vtx.Y() >  500.0) return false;
      if(vtx.Z() <   50.0) return false;
      if(vtx.Z() > 5450.0) return false;
      return true;
    });

  const Cut kNusFDFidLoose(
    [](const caf::SRProxy* sr)
    {
      double buf = 10.;
      if(!sr->vtx.elastic.IsValid) return false;

      const TVector3 vtx = sr->vtx.elastic.vtx;
      if(vtx.X() < (FDL + buf)) return false;
      if(vtx.X() > (FDR - buf)) return false;
      if(vtx.Y() < (FDB + buf)) return false;
      if(vtx.Y() > (FDT - buf)) return false;
      if(vtx.Z() < (FDF + buf)) return false;
      if(vtx.Z() > (FDE - buf)) return false;
      return true;
    });

  const Cut kNusFDContain(
    [](const caf::SRProxy* sr)
    {
      const caf::SRNueCosRejProxy& cr = sr->sel.nuecosrej;

      if(std::min(cr.starteast,   cr.stopeast)   < 10) return false;
      if(std::min(cr.startwest,   cr.stopwest)   < 10) return false;
      if(std::min(cr.starttop,    cr.stoptop)    < 10) return false;
      if(std::min(cr.startbottom, cr.stopbottom) < 10) return false;
      if(std::min(cr.startfront,  cr.stopfront)  < 10) return false;
      if(std::min(cr.startback,   cr.stopback)   < 10) return false;
      return true;
    });

  const Cut kNusNCSel(
    [](const caf::SRProxy* sr)
    {
      if(sr->slc.nhit >= 200)    return false;
      if(sr->slc.nhit <  20)     return false;
      if(sr->sel.cvn.ncid < 0.2) return false;

      return true;
    });

  const Cut kNusFDDecafCut(
    [](const caf::SRProxy* sr)
    {
      if(sr->sel.cosrej.numucontpid2019 <= 0.42) return false;
      return true;
    });

  const Cut kNusCosRej(
    [](const caf::SRProxy* sr)
    {
      double numucontpid2019 = sr->sel.cosrej.numucontpid2019;
      double partptp = sr->sel.nuecosrej.partptp;
      if(sr->vtx.elastic.fuzzyk.nshwlid == 0) return false;
      if(numucontpid2019 <= 0.5) return false;
      if(partptp >= 0.8) return false;
      if(sr->vtx.elastic.fuzzyk.png[0].shwlid.lid.ismuon == 1) return false;
      double nhit = (double)sr->slc.nhit;
      if(nhit <= 0.) return false;
      if(sr->slc.calE/nhit <= 0.018) return false;
      if(std::min(sr->sel.nuecosrej.starttop, sr->sel.nuecosrej.stoptop) < 480) return false;
      return true;
    });

  const Cut kNusNDFiducial(
    [](const caf::SRProxy* sr)
    {
      assert(sr->vtx.elastic.IsValid && "Must apply DQ cuts");
      if(sr->vtx.elastic.vtx.X() < -100.0) return false;
      if(sr->vtx.elastic.vtx.X() >  100.0) return false;
      if(sr->vtx.elastic.vtx.Y() < -100.0) return false;
      if(sr->vtx.elastic.vtx.Y() >  100.0) return false;
      if(sr->vtx.elastic.vtx.Z() <  200.0) return false;
      if(sr->vtx.elastic.vtx.Z() > 1000.0) return false;
      return true;
    });

  const Cut kNusNDFidLoose(
    [](const caf::SRProxy* sr)
    {
      double buf = 10.;
      if(!sr->vtx.elastic.IsValid) return false;

      const TVector3 vtx = sr->vtx.elastic.vtx;
      if(vtx.X() < (NDL + buf)) return false;
      if(vtx.X() > (NDR - buf)) return false;
      if(vtx.Y() < (NDB + buf)) return false;
      if(vtx.Y() > (NDT - buf)) return false;
      if(vtx.Z() < (NDF + buf)) return false;
      if(vtx.Z() > (NDE - buf)) return false;

      return true;
    });

  const Cut kNusNDContain(
    [](const caf::SRProxy* sr)
    {
      const caf::SRNueCosRejProxy& cr = sr->sel.nuecosrej;

      if(std::min(cr.starteast,   cr.stopeast)   < 25) return false;
      if(std::min(cr.startwest,   cr.stopwest)   < 25) return false;
      if(std::min(cr.starttop,    cr.stoptop)    < 25) return false;
      if(std::min(cr.startbottom, cr.stopbottom) < 25) return false;
      if(std::min(cr.startfront,  cr.stopfront)  < 25) return false;
      if(std::min(cr.startback,   cr.stopback)   < 25) return false;
      return true;
    });

  const Cut kNusNDHarshTrk(
    [](const caf::SRProxy* sr)
    {
      double buf = 25.;
      if(sr->trk.kalman.ntracks == 0) return false;

      for(unsigned int i = 0, n = sr->trk.kalman.ntracks; i < n; ++i) {
        TVector3 start = sr->trk.kalman.tracks[i].start;
        TVector3 stop  = sr->trk.kalman.tracks[i].stop;
        if(std::min(start.X(), stop.X()) < (NDL + buf)) return false;
        if(std::max(start.X(), stop.X()) > (NDR - buf)) return false;
        if(std::min(start.Y(), stop.Y()) < (NDB + buf)) return false;
        if(std::max(start.Y(), stop.Y()) > (NDT - buf)) return false;
        if(std::min(start.Z(), stop.Z()) < (NDF + buf)) return false;
        if(std::max(start.Z(), stop.Z()) > (NDE - buf)) return false;
      }

      return true;
    });

  const Cut kNusCosRejMod(
    [](const caf::SRProxy* sr)
    {
      double partptp = sr->sel.nuecosrej.partptp;
      if(sr->vtx.elastic.fuzzyk.nshwlid == 0) return false;
      if(partptp >= 0.8) return false;
      if(sr->vtx.elastic.fuzzyk.png[0].shwlid.lid.ismuon == 1) return false;
      double nhit = (double)sr->slc.nhit;
      if(nhit <= 0.) return false;
      if(sr->slc.calE/nhit <= 0.009) return false;
      return true;
    });
}