NueCuts2020.h

Go to the documentation of this file.

#pragma once
#include "CAFAna/Cuts/AnalysisMasks.h"
#include "3FlavorAna/Cuts/NueCuts2018.h"
#include "3FlavorAna/Vars/NueEnergy2020.h"
#include "3FlavorAna/Cuts/3flavorCuts.h"

//VERY PRELIMINARY

namespace ana
{

//*************************************************************
//********************   Reco quality   ***********************
//*************************************************************

  extern const Cut kNue2020RecoQuality;

//*************************************************************
//**************************  FD  *****************************
//*************************************************************

  //extern const Cut kNue2020FDVeto;

  extern const Cut kNue2020FDNearestSlice;

  const Cut kNue2020FDBasicQuality =    kApplySecondAnalysisMask && 
                                        //kNue2020FDVeto &&
                                        k3flavor2020FDVeto &&
                                        kNue2020RecoQuality;

  const Cut kNue2020FDDecafCut = kNue2020FDBasicQuality && kCVNe_looseptp >= 0.5;
  const Cut kNue2020CoreBasicEventCut = kNueEnergy2020 > 1 && kNueEnergy2020 < 4 &&
                                        kNHit > 30 && kNHit < 150 &&
                                        kLongestProng > 100 && kLongestProng < 500;

  // Inherited from second ana, re-tuned/validated docDB 42840/42964
  const Cut kNue2020ProngContain = kDistAllTop   > 63 && kDistAllBottom > 12 &&
                                   kDistAllEast  > 12 && kDistAllWest   > 12 &&
                                   kDistAllFront > 18 && kDistAllBack   > 18;

  const Cut kNue2020CorePart = kNue2020CoreBasicEventCut &&
                               kNue2020ProngContain;
 
  const Cut kNue2020CorePresel = kNue2020FDBasicQuality &&  kNue2020CorePart;

  // CVN FHC boundaries
  const double kNue2020CVNFHCLowEdge = 0.84;
  const double kNue2020CVNFHCHighEdge = 0.97;
  const double kNue2020CosRejFHCEdge = 0.49;

  // CVN RHC boundaries
  const double kNue2020CVNRHCLowEdge = 0.85;
  const double kNue2020CVNRHCHighEdge = 0.97;
  const double kNue2020CosRejRHCEdge = 0.47;

  const Cut kNue2020PIDFHC = kCVNe_looseptp >= kNue2020CVNFHCLowEdge;
  const Cut kNue2020PIDRHC = kCVNe_looseptp >= kNue2020CVNRHCLowEdge;
  const Cut kNue2020CosRejFHC = kCosPIDCoreBDT >= kNue2020CosRejFHCEdge;
  const Cut kNue2020CosRejRHC = kCosPIDCoreBDT >= kNue2020CosRejRHCEdge;

  extern const Cut kNue2020PID;
  extern const Cut kNue2020CosRej;
  
  // TODO update when there is new cosrej cuts
  const Cut kNue2020FD = kNue2020CorePresel && kNue2020CosRej && kNue2020PID && kNue2020FDNearestSlice;

  //***********************  Peripheral  ***********************

  const Cut kNue2020PeripheralBasicEventCut = kIsFarDet &&  
                                              kNue2020FDBasicQuality && 
                                              kNueEnergy2020 < 4.5 && 
                                              kNueEnergy2020 >= 0; 

  const Cut kNue2020CVN_VsCosPIDRHC = ((kCVNe_looseptp >= kNue2020CVNRHCHighEdge && kCosPIDPeriBDT > 0.6) || (kCVNe_looseptp>=0.995 && kCosPIDPeriBDT > 0.56));
  const Cut kNue2020CVN_VsCosPIDFHC = ((kCVNe_looseptp >= kNue2020CVNFHCHighEdge && kCosPIDPeriBDT > 0.61) || (kCVNe_looseptp>=0.995 && kCosPIDPeriBDT > 0.57));
  
  bool kNue2020CVN_VsCosPIDFunc(const caf::SRProxy* sr);
  const Cut kNue2020CVN_VsCosPID(kNue2020CVN_VsCosPIDFunc);
  
  /// Use this cut for the full FD Peripheral sample, the same for RHC and FHC
  bool kNue2020FDPeripheralFunc(const caf::SRProxy* sr);
  const Cut kNue2020FDPeripheral(kNue2020FDPeripheralFunc);

  const Cut kNue2020FDAllSamples = kNue2020FD || kNue2020FDPeripheral;
  
  //  ---------------------  FD Sidebands  ---------------------  //
  const Cut kNue2020FHCCVNCutLP = kCVNe_looseptp >= 0.5 && kCVNe_looseptp < 0.65;
  const Cut kNue2020RHCCVNCutLP = kCVNe_looseptp >= 0.5 && kCVNe_looseptp < 0.65;
  extern const Cut kNue2020CVNCutLP;

  ///////////////////////////////////////////////////////////////////////////
  //
  //  ---------------------  FD Peripheral Sidebands  ---------------------  //
  //
  ///////////////////////////////////////////////////////////////////////////
  const Cut kNue2020PeripheralBasicEventCutHE = kIsFarDet &&  
                                                kNue2020FDBasicQuality && 
                                                kNueEnergy2020 > 4.5 && 
                                                kNueEnergy2020 <= 10; 

  const Cut kNue2020CVNVsCosPIDRHCLP = kCVNe_looseptp >= 0.5 && kCosPIDPeriBDT>0.56;
  const Cut kNue2020CVNVsCosPIDFHCLP = (kCVNe_looseptp >= 0.5 && kCosPIDPeriBDT>0.57);
  extern const Cut kNue2020CVNVsCosPIDLP;

  /// Use this cut for the sideband FD Peripheral samples, the same for RHC and FHC
  const Cut kNue2020FDPeripheralHE =  kNue2020PeripheralBasicEventCutHE && kNue2020CVN_VsCosPID && kNue2020FDNearestSlice && !kNue2020FD;
  const Cut kNue2020FDPeripheralLP =  kNue2020PeripheralBasicEventCut && kNue2020CVNVsCosPIDLP && kNue2020FDNearestSlice && !kNue2020FD;


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

  //------------------- Preselection Cut -------------------
  const Cut kNue2020CoreBasicEventCutHE =
    kNueEnergy2020 > 4.5 &&
    kNHit > 30 &&
    kLongestProng > 100 && kLongestProng < 500;
  const Cut kNue2020CorePartHE = kNue2020CoreBasicEventCutHE &&
                                 kNue2020ProngContain;
 
  const Cut kNue2020CorePreselHE = kNue2020FDBasicQuality &&  kNue2020CorePartHE;

  /// Use this cut for the sideband FD Core selections, apply for both RHC and FHC
  const Cut kNue2020FDHE = kNue2020CorePreselHE && kNue2020CosRej && kNue2020PID && kNue2020FDNearestSlice;
  const Cut kNue2020FDLP = kNue2020CorePresel && kNue2020CVNCutLP && kNue2020CosRej && kNue2020FDNearestSlice;

  const Cut kNue2020FDAllSamplesHE = kNue2020FDHE || kNue2020FDPeripheralHE;
  const Cut kNue2020FDAllSamplesLP = kNue2020FDLP || kNue2020FDPeripheralLP;

  ///////////////////////////////////////////////////////////////////////////
  //
  //  ---------------------     MRBrem cuts         ---------------------  //
  //
  ///////////////////////////////////////////////////////////////////////////
  // Core
  const Cut kNue2020FDBasicQuality_MRBrem = kApplySecondAnalysisMask && 
                                            //kCosRejVeto && // CNN veto is not applicable for mrbrem files
                                            kNue2020RecoQuality;

  const Cut kNue2020CorePresel_MRBrem = kNue2020FDBasicQuality_MRBrem &&  kNue2020CorePart;
  const Cut kNue2020FD_MRBrem = kNue2020CorePresel_MRBrem && kNue2020CosRej && kNue2020PID && kNue2020FDNearestSlice;

  // Peripheral
  const Cut kNue2020PeripheralBasicEventCut_MRBrem = kIsFarDet &&  
                                                     kNue2020FDBasicQuality_MRBrem && 
                                                     kNueEnergy2020 < 4.5 && 
                                                     kNueEnergy2020 >= 0;

  bool kNue2020FDPeripheralFunc_MRBrem(const caf::SRProxy* sr);
  const Cut kNue2020FDPeripheral_MRBrem(kNue2020FDPeripheralFunc_MRBrem);

//*************************************************************
//**************************  ND  *****************************
//*************************************************************

  extern const Cut kNue2020NDFrontPlanes;
  const Cut kNue2020NDNHits =       kNHit >= 20 && kNHit <= 200;
  const Cut kNue2020NDEnergy =      kNueEnergy2020 < 4.5 && kNueEnergy2020 >= 0;
  const Cut kNue2020NDProngLength = kLongestProng > 100 &&
                                    kLongestProng < 500;

  // Inherited from second ana, re-tuned/validated docDB 42840/42964
  extern const Cut kNue2020NDFiducial;
  extern const Cut kNue2020NDContain;

  const Cut kNue2020NDPresel = kNue2020RecoQuality &&
                               kNue2020NDFiducial &&
                               kNue2020NDContain &&
                               kNue2020NDNHits &&
                               kNue2020NDEnergy &&
                               kNue2020NDProngLength;

  const Cut kNue2020NDDecafCut = kNue2020RecoQuality && kNue2020NDFiducial && 
                                 kNue2020NDContain;
  const Cut kNue2020ND  = kNue2020NDPresel && kNue2020PID;

//*************************************************************
//**************  Energy/CVN Binning for 2020 ****************
//*************************************************************

//will be updated once we have our cvn tuning finalized
//TODO: looks weird for ND test... should be ok for ND/FD binning

  // 2bin layout for 2018, ana2017 style (the third bin is for peripheral sample)
  extern const Var kNue2020SelectionBin;

  /// Use this Analysis Binning for Ana2020, official Binning
  extern const Var kNue2020AnaBin;

  extern const Binning kNue2020Binning;
  /// Use this Axis for Ana2020, official Axis
  extern const HistAxis kNue2020Axis;

  extern const Var kNue2020AnaBinMergedPeripheral;

  extern const Binning kNue2020BinningMergedPeripheral;
  extern const HistAxis kNue2020AxisMergedPeripheral;
 
}