NueCuts2017.h

Go to the documentation of this file.

#pragma once

// final summary for 2017 doc-db ?
//
// Blank for summer analysis cuts
// Only semi-finalized FD cuts for now
//

#include "CAFAna/Core/HistAxis.h"
#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "3FlavorAna/Vars/NueEnergy2017.h"

namespace ana
{

  ///////////////////////////////////////////////////////////////////////////
  //                                                                       //
  //  -------------------  Cuts Common to Both Detectors  ---------------  //
  //                                                                       //
  ///////////////////////////////////////////////////////////////////////////


  //--------------------------- DQ Cuts ---------------------------
  // CVN, LID doc db 19011, 19550

  //--------- CVN ---------
  extern const Cut kNueDQ2017CVN;

  //--------- LID ---------
  extern const Cut kNueDQ2017LID;

  /// Close to the s/sqrt(s+b) optimization.  TODO: DocDB reference?
  const Cut kNue2017CVNCut = kCVNe > 0.75;

  ///////////////////////////////////////////////////////////////////////////
  //                                                                       //
  //  -----------------------  FD specific  cuts  -----------------------  //
  //                                                                       //
  ///////////////////////////////////////////////////////////////////////////

  /// cosmic veto cut for FD
  extern const Cut kNue2017Veto;


  //------------------- Basic Cuts ----------------------
  // ALWAYS apply these
  bool kNue2017BasicPartFunc(const caf::SRProxy* sr);
  const Cut kNue2017BasicPart(kNue2017BasicPartFunc);


  //------------------- Preselection Cut -------------------
  // doc-db 17176
  const Cut kNue2017Presel = kNueEnergy2017 > 1 && kNueEnergy2017 < 4 &&
                             kNHit > 30 && kNHit < 150 &&
                             kLongestProng > 100 && kLongestProng < 500;

  //----------------------- Pt/P Cut -----------------------
  // doc-db 18105, 19550
  //--------- CVN ---------
  //
  /// "L-shaped" Cut in p_t/P vs maxy space to reject cosmic rays for the "core"
  /// samples, suitable for use with CVN.
  /// Where "X" is rejected and "." is selected.
  ///
  /// M ..XXXX
  /// a ..XXXX
  /// x ....XX
  /// y ....XX
  ///   ......
  ///  P_t / p
  const Cut kNue2017PtPCutCVN = ( kPtP < 0.58 ) || ( 0.58 <= kPtP && kPtP < 0.8 && kMaxY < 590 ) || ( 0.8 <= kPtP && kMaxY < 350 );

  //--------- LID ---------
  /// Somewhat tighter p_t/P vs maxy cuts suitable for use with LID
  const Cut kNue2017PtPCutLID = ( kPtP < 0.4 ) || ( 0.4 <= kPtP && kPtP < 0.65 && kMaxY < 600 );

  /// FD containment cuts, inherited from Second Analysis
  const Cut kNue2017ProngContainment =
    kDistAllTop    > 63  &&
    kDistAllBottom > 12  &&
    kDistAllEast   > 12  &&
    kDistAllWest   > 12  &&
    kDistAllFront  > 18  &&
    kDistAllBack   > 18;

  /// FD backward cosmic photon rejection cut, inherited from Second Analysis
  const Cut kNue2017BackwardCut =
    ( kDistAllBack < 200 && kSparsenessAsymm < -0.1 ) ||
    ( kDistAllBack >= 200 );


  /// Cuts which define the "core" sample and are reversed for the "peripheral" sample
  const Cut kNue2017CorePart =
    kNue2017Presel &&
    kNue2017ProngContainment &&
    kNue2017PtPCutCVN &&
    kNue2017BackwardCut;

  /// Cuts which define the "core" sample and are reversed for the "peripheral"
  /// sample. `ProngContainment`, `PtPCutCVN` and `BackwardCut` are unchanged
  /// from 2017 analysis at the moment.

  /// Full preselection for 2017 "core" sample with Cut-based CosRej
  const Cut kNue2017CorePresel = kNue2017BasicPart && kNue2017CorePart;

  /// Full preselection for 2017 "peripheral" sample
  /// To save us the headache of figuring out how to plot energy overflow
  /// on our FD distribution plot, only consider events with reco E < 4.5 GeV
  const Cut kNue2017PeripheralPresel =
    kNue2017BasicPart &&
    !kNue2017CorePart &&
    kNueEnergy2017 < 4.5 &&
    kIsFarDet;

  /// Lowest BDT value used in Matryoshka binning DocDB 19724.  Not finalized
  const Cut kNue2017BDTContainment = kCosPID > 0.56;

  /// Full FD nue preselection cuts for 2017 analysis with BDT-based CosRej
  const Cut kNue2017FullPreselBDTCosRej =
    kNue2017BasicPart &&
    kNue2017Presel &&
    kNue2017BDTContainment;

  // Preliminary nearest slice cuts for cosmic rejection based on delayed brem showers
  // TODO : Docdb reference
  extern const Cut kNearestSliceCosRejCore;

  /// Full FD nue selection for 2017 analysis
  const Cut kNue2017FD = kNue2017CorePresel && kNue2017CVNCut && kNearestSliceCosRejCore;

  /// Cosmic rejection cut in CVNe vs CosPIDContain space for peripheral sample
  ///
  /// This is a box of high CVN and high BDT space that is selected, plus a
  /// sliver of very-high CVN events that are selected regardless of BDT.
  ///
  /// Where "X" is rejected and "." selected
  ///
  /// C XXX...
  /// o XXX...
  /// s XXX...
  /// B XXXXX.
  /// D XXXXX.
  /// T XXXXX.
  ///    CVNe
  const Cut kNue2017CVNVsCosPID = (kCVNe>0.95 && kCosPIDContain>0.53) || kCVNe>0.99;

  // Nearest slice cuts for the peripheral bin
  extern const Cut kNearestSliceCosRejPeripheral;

  /// Definition of the FD peripheral sample around our core selected events
  const Cut kNue2017FDPeripheral = kNue2017PeripheralPresel && kNue2017CVNVsCosPID && kNearestSliceCosRejPeripheral;

  /// Our FD selection including all samples, for making predictions, etc.
  const Cut kNue2017FDAllSamples = kNue2017FD || kNue2017FDPeripheral;

  const Cut kNueFD2017DecafCut = kNue2017BasicPart && kCVNe > 0.5;


///////////////////////////////////////////////////////////////////////////
//
//  ---------------------  FD Peripheral Sidebands  ---------------------  //
//
///////////////////////////////////////////////////////////////////////////

  const Cut kNue2017PeripheralPreselHE =
    kNue2017BasicPart &&
    !kNue2017CorePart &&
    kNueEnergy2017 > 4.5 &&
    kNueEnergy2017 < 10 &&
    kIsFarDet;

  const Cut kNue2017CVNVsCosPIDPeripheralMP = (kCVNe > 0.87 && kCVNe<=0.95
                                               && kCosPIDContain>0.53);

  const Cut kNue2017CVNVsCosPIDPeripheralLP = (kCVNe > 0.75 && kCVNe<=0.95
                                               && kCosPIDContain>0.53);

  // Nearest slice cuts for the peripheral sidebands
  extern const Cut kNearestSliceCosRejPeripheralSideband;

  /// Definition of the FD peripheral High Energy Low PID
  const Cut kNue2017FDPeripheralHELPsideband =
    kNue2017PeripheralPreselHE &&
    kNue2017CVNVsCosPIDPeripheralLP &&
    kNearestSliceCosRejPeripheralSideband;

  /// Definition of the FD peripheral High Energy Mid PID
  const Cut kNue2017FDPeripheralHEMPsideband =
    kNue2017PeripheralPreselHE &&
    kNue2017CVNVsCosPIDPeripheralMP &&
    kNearestSliceCosRejPeripheralSideband;

  /// Definition of the FD peripheral Low PID
  const Cut kNue2017FDPeripheralLPsideband =
    kNue2017PeripheralPresel &&
    kNue2017CVNVsCosPIDPeripheralLP &&
    kNearestSliceCosRejPeripheralSideband;

  /// Definition of the FD peripheral Mid PID
  const Cut kNue2017FDPeripheralMPsideband =
    kNue2017PeripheralPresel &&
    kNue2017CVNVsCosPIDPeripheralMP &&
    kNearestSliceCosRejPeripheralSideband;

  /// Definition of the FD peripheral High Energy sideband
  const Cut kNue2017FDPeripheralHEsideband =
    kNue2017PeripheralPreselHE &&
    kNue2017CVNVsCosPID &&
    kNearestSliceCosRejPeripheralSideband;

  ///////////////////////////////////////////////////////////////////////////
  //
  //  ---------------------     FD Full Sidebands    ---------------------  //
  //
  ///////////////////////////////////////////////////////////////////////////

  //------------------- Preselection Cut -------------------
  const Cut kNue2017PreselHE =
    kNueEnergy2017 > 4.5 &&
    kNHit > 30 &&
    kLongestProng > 100 && kLongestProng < 500;

  /// Cuts which define the "core" HE sample and are reversed for the "peripheral" sample
  const Cut kNue2017CorePartHE =
    kNue2017PreselHE &&
    kNue2017ProngContainment &&
    kNue2017PtPCutCVN &&
    kNue2017BackwardCut;

  /// CosRej for the core sidebands
  extern const Cut kNearestSliceCosRejCoreSideband;

  /// Preselection for 2017 "core" HE sample with Cut-based CosRej
  const Cut kNue2017CorePreselHE = kNue2017BasicPart && kNue2017CorePartHE;

  /// FD nue selection for "core" HE sample
  const Cut kNue2017FDHE = kNue2017CorePreselHE && kNue2017CVNCut && kNearestSliceCosRejCoreSideband;

  /// All samples - High energy sideband
  const Cut kNue2017FDAllSamplesHEsidebands = kNue2017FDHE || kNue2017FDPeripheralHEsideband;

  /// Low PID cut for core sidebands
  const Cut kNue2017CVNCutLP = kCVNe > 0.50 && kCVNe <= 0.75;

  /// FD nue selection for 2017 low PID sideband
  const Cut kNue2017FDLP = kNue2017CorePresel && kNue2017CVNCutLP && kNearestSliceCosRejCoreSideband;

///////////////////////////////////////////////////////////////////////////
//
//  -----------------------  ND specific  cuts  -----------------------  //
//
///////////////////////////////////////////////////////////////////////////

  /// Cut on ND vertex position
  extern const Cut kNue2017NDFiducial;

  /// Loose containtment on start and end of all showers, docdb-12943
  extern const Cut kNue2017NDContain;

  /// No hits in the first 5 planes of the detector, doc 12879
  extern const Cut kNue2017NDFrontPlanes;

  // Make sure you cover extrap
  const Cut kNue2017NDNHits = kNHit >= 20 && kNHit <= 200;

  // Make sure you cover extrap in high energy sideband
  const Cut kNue2017NDNHitsHE = kNHit >= 20;

  // Make sure you cover extrap
  const Cut kNue2017NDEnergy = kNueEnergy2017 < 4.5;

  // Make sure you cover extrap in high energy sideband
  const Cut kNue2017NDEnergyHE = kNueEnergy2017 < 10;

  // Same as FA ND / FD cuts
  const Cut kNue2017NDProngLength = kLongestProng > 100 &&
                                    kLongestProng < 500;

  const Cut kNueND2017DecafCut =  kNueDQ2017CVN &&
                                  kNue2017NDFiducial &&
                                  kNue2017NDContain &&
                                  kNue2017NDFrontPlanes;

  const Cut kNue2017NDPresel = kNueDQ2017CVN &&
                               kNue2017NDFiducial &&
                               kNue2017NDContain &&
                               kNue2017NDFrontPlanes &&
                               kNue2017NDNHits &&
                               kNue2017NDEnergy &&
                               kNue2017NDProngLength;

  const Cut kNue2017NDPreselHE = kNueDQ2017CVN &&
                                 kNue2017NDFiducial &&
                                 kNue2017NDContain &&
                                 kNue2017NDFrontPlanes &&
                                 kNue2017NDNHitsHE &&
                                 kNue2017NDEnergyHE &&
                                 kNue2017NDProngLength;


  const Cut kNue2017NDCVNSsb = kNue2017NDPresel && kNue2017CVNCut;

  // Actually a Var, but it depends on a cut, so including it here.
  // This properly lays out our FD events into 3 CVN bins + peripheral events
  // In the ND, this should just fill the 3 core CVN bins without events in
  // the peripheral portion.  ModularExtrap will then just return the FD MC
  // for the background prediction for our background.  That may be good
  // enough, or we'll have to rethink.
  extern const Var kNue2017SelectionBin;
  extern const Var kNue2017AnaBin;
  extern const Binning kNue2017Binning;
  extern const HistAxis kNue2017Axis;

  // Similarly, for merging the peripheral sample into a single bin, we need
  // a new var / HistAxis.  Since the nominal analysis will need to extrap the
  // beam bkgd in reco-E, we need to extrap in energy and merge later, but for
  // simple studies, this is useful
  extern const Var kNue2017AnaBinMergedPeripheral;
  extern const Binning kNue2017BinningMergedPeripheral;
  extern const HistAxis kNue2017AxisMergedPeripheral;


  ///////////////////////////////////////////////////////////////////////////
  //                                                                       //
  //  -----------------------  MRE specific  cuts  ----------------------- //
  //                                                                       //
  ///////////////////////////////////////////////////////////////////////////


  // TODO: numu's probably changing their selection to include CVN, we'll
  // need to update this cut
  // This cut needs to be applied to the mre events, in addition to nue
  // selection cuts. It makes sure that the parent of the mre event was
  // a valid numu CC candidate
  extern const Cut kNue2017MRParentSliceCut;
}