MRDiFCutsAndVars.h

Go to the documentation of this file.

#pragma once

//#include "CAFAna/Core/Cut.h"
#include "CAFAna/Vars/Vars.h"
//#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/AnalysisMasks.h"
#include "3FlavorAna/Vars/NueEnergy2020.h"
#include "3FlavorAna/Cuts/NueCuts2020.h"
#include "3FlavorAna/Cuts/3flavorCuts.h"

#include "StandardRecord/Proxy/SRProxy.h"
#include <cassert>

/// Cuts and Vars for the 2020 FD DiF Study



namespace ana {


  //--------Vars-----------------------------------------
  
  //These are used to mimic the kNue2020ProngContain cut
  
  const Var kdistTop([](const caf::SRProxy* sr)
                     {
                       if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                       float y_val  = sr->vtx.elastic.vtx.y;
                       for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                         {
                           float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.y;
                           if (tmp > y_val) y_val = tmp;
                         }
                       float dist_to_top = 765 - y_val;
                       return dist_to_top;
                     });
  
  const Var kdistBottom([](const caf::SRProxy* sr)
                        {
                          if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                          float y_val  = sr->vtx.elastic.vtx.y;
                          for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                            {
                              float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.y;
                              if (tmp < y_val) y_val = tmp;
                            }
                          float dist_to_bot = y_val + 749;
                          return dist_to_bot;
                        });
  
  const Var kdistFront([](const caf::SRProxy* sr)
                       {
                         if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                         float z_val  = sr->vtx.elastic.vtx.z;
                         for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                           {
                             float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.z;
                             if (tmp < z_val) z_val = tmp;
                           }
                         return z_val;
                       });

  const Var kdistBack([](const caf::SRProxy* sr)
                      {
                        if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                        float z_val  = sr->vtx.elastic.vtx.z;
                        for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                          {
                            float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.z;
                            if (tmp > z_val) z_val = tmp;
                           }
                        return 5962 - z_val;
                      });
  
  const Var kdistWest([](const caf::SRProxy* sr)
                      {
                        if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                        float x_val  = sr->vtx.elastic.vtx.x;
                        for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                          {
                            float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.x;
                            if (tmp < x_val) x_val = tmp;
                          }
                        float dist_to_west = x_val + 758;
                        return dist_to_west;
                      });

  const Var kdistEast([](const caf::SRProxy* sr)
                      {
                        if (!sr->vtx.elastic.IsValid || !(sr->vtx.elastic.fuzzyk.png.size() >0)) return (float)-5;
                        float x_val  = sr->vtx.elastic.vtx.x;
                        for(unsigned int pngID = 0; pngID < sr->vtx.elastic.fuzzyk.png.size(); pngID ++)
                          {
                             float tmp = sr->vtx.elastic.fuzzyk.png[pngID].shwlid.stop.x;
                             if (tmp > x_val) x_val = tmp;
                          }
                        float dist_to_east = 765 - x_val;
                        return dist_to_east;
                      });
  
  /*
  extern const Var kdistTop;
  extern const Var kdistBottom;
  extern const Var kdistFront;
  extern const Var kdistBack;
  extern const Var kdistEast;
  extern const Var kdistWest;
  */
  
  //------------Cuts-------------------------------------
  
  const Cut kDiFContainment = (kdistTop > 53 && kdistBottom > 12
                               && kdistFront > 18 && kdistBack > 18
                               && kdistEast > 12 && kdistWest >12);
  
  const Cut kSingleSlice([](const caf::SRProxy* sr) { return (sr->hdr.subevt < 2); });
  
  const Cut kDiFBasicQuality([](const caf::SRProxy* sr) 
                             {
                               if(!sr->vtx.elastic.IsValid)
                                 return false;
                               if(sr->vtx.elastic.fuzzyk.npng == 0)
                                 return false;
                               return (kApplySecondAnalysisMask(sr) && kCosRejVeto(sr));
                             } );
  
  const Cut kLooseCVNeCut ([](const caf::SRProxy* sr) { return sr->sel.cvnloosepreselptp.nueid > 0.3;});
  
  /*
  extern const Cut kSingleSlice;
  extern const Cut kDiFBasicQuality;
  extern const Cut kLooseCVNeCut;
  */

  //Final preselection and selection cuts
  //Note: in TimingCuts.h, there is MR studies specific timing cut to use as well

  const Cut kDiFPreselectionCut = (kSingleSlice && kDiFBasicQuality 
                                   && kNue2020CoreBasicEventCut && kLooseCVNeCut 
                                   && kDiFContainment);
  const Cut kDiFSelectionCut = kDiFPreselectionCut && kNue2020PID;
  
}