NumubarCCpi0_Cuts.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/Cut.h"
#include "NDAna/numubarcc_pi0/NumubarCCpi0_Vars.h"
#include "StandardRecord/Proxy/SRProxy.h"

#include "TMath.h"

namespace ana
{
  namespace numubarccpi0
  {
    //====================//
    //     Truth cuts     //
    //====================//
    extern const Cut kTrueFiducial;
    extern const Cut kIsNumuCC; 
    extern const Cut kIsNC;     
    extern const Cut kIsNu;     
    extern const Cut kIsAntiNu; 

    const Cut kPrimaryPi0   = kNPrimaryPi0 > 0;
    const Cut kSecondaryPi0 = kNSecondaryPi0 > 0;

    //======================//
    //     Preselection     //
    //======================//
    extern const Cut kThreePngPlus;
    extern const Cut kIsFiducial;
    extern const Cut kNumuMyQuality;
    extern const Cut kNumuTightContainND;

    const Cut kMuonIDCut = kMuonID > 0.24;

    const Cut kAllNumuCCCuts =
              kMuonIDCut &&        // PID cut
              kIsFiducial &&       // fiducial cut
              kNumuMyQuality &&    // ensures reconstruction
              kNumuTightContainND; // containment

    //=====================//
    //     Signal cuts     //
    //=====================//
    const Cut kNumuCC = kTrueFiducial && kIsNumuCC;
    const Cut kNC = kTrueFiducial && kIsNC;

    const Cut kSignal = kNumuCC && kPrimaryPi0;
    const Cut kNuSig = kSignal && kIsNu;
    const Cut kAntiNuSig = kSignal && kIsAntiNu;
    const Cut kBkg = !kSignal;

    const Cut kPrimPi0 = kPrimaryPi0;
    const Cut kNoPrimPi0 = !kPrimaryPi0;
    const Cut kSecoPi0 = !kPrimaryPi0 && kSecondaryPi0;
    const Cut kNoPi0 = !kPrimaryPi0 && !kSecondaryPi0;
    const Cut kHasPi0 = kPrimaryPi0 || kSecondaryPi0;

    const Cut kIsOther = !kTrueFiducial || (kTrueFiducial && !kIsNC && !kIsNumuCC);

    const Cut kSignal_fixed = kNumuCC && kHasPi0;
    const Cut kBkg_fixed = !kSignal_fixed;

    //=========================//
    //     Prong NHit Cuts     //
    //=========================//
    const Cut kProng1NHitCut4_4view  = kProng1NHit_4view  >= 4;
    const Cut kProng1NHitCut4_emid   = kProng1NHit_emid   >= 4;
    const Cut kProng1NHitCut4_5label = kProng1NHit_5label >= 4;

    const Cut kProng2NHitCut4_4view  = kProng2NHit_4view  >= 4;
    const Cut kProng2NHitCut4_emid   = kProng2NHit_emid   >= 4;
    const Cut kProng2NHitCut4_5label = kProng2NHit_5label >= 4;

    const Cut kProng1NHitCut6_4view  = kProng1NHit_4view  >= 6;
    const Cut kProng1NHitCut6_emid   = kProng1NHit_emid   >= 6;
    const Cut kProng1NHitCut6_5label = kProng1NHit_5label >= 6;

    const Cut kProng2NHitCut6_4view  = kProng2NHit_4view  >= 6;
    const Cut kProng2NHitCut6_emid   = kProng2NHit_emid   >= 6;
    const Cut kProng2NHitCut6_5label = kProng2NHit_5label >= 6;

    const Cut kProngNHitCut4_4view  = kProng1NHitCut4_4view  && kProng2NHitCut4_4view;
    const Cut kProngNHitCut4_emid   = kProng1NHitCut4_emid   && kProng2NHitCut4_emid;
    const Cut kProngNHitCut4_5label = kProng1NHitCut4_5label && kProng2NHitCut4_5label;

    const Cut kProngNHitCut6_4view  = kProng1NHitCut6_4view  && kProng2NHitCut6_4view;
    const Cut kProngNHitCut6_emid   = kProng1NHitCut6_emid   && kProng2NHitCut6_emid;
    const Cut kProngNHitCut6_5label = kProng1NHitCut6_5label && kProng2NHitCut6_5label;

    //========================//
    //     Selection Cuts     //
    //========================//
    const Cut kPresel = kNumuMyQuality && kNumuTightContainND && kIsFiducial;
    const Cut kNumuCCsel = kPresel && kMuonIDCut;
    const Cut kMuonE40Cut = kRecoMuKE > 0.40;
    const Cut kPngPresel = kNumuCCsel && kMuonE40Cut && kThreePngPlus;

    const Cut kProngScoreCut_4view = kProng1Score_4view > 0.5 && kProng2Score_4view > 0.5;
    const Cut kProngScoreCut_emid_Cut4 = kProng2Score_emid > 0.22 && kProng1Score_emid > 0.32;
    const Cut kProngScoreCut_emid_Cut6 = kProng2Score_emid > 0.17 && kProng1Score_emid > 0.17;
    const Cut kProngScoreCut_5label_Cut4 = kProng2Score_5label > 0.21 && kProng1Score_5label > 0.27;
    const Cut kProngScoreCut_5label_Cut6 = kProng2Score_5label > 0.17 && kProng1Score_5label > 0.27;

  } // end of namespace numubarccpi0
} // end of namespace ana