NumuCuts2018.cxx
Go to the documentation of this file.
2 
4 
5 #include <iostream>
6 
7 namespace ana
8 {
9  const Cut kLooseKFC = SIMPLEVAR(sel.contain.kalfwdcell) > 2;
10  const Cut kLooseKBC = SIMPLEVAR(sel.contain.kalbakcell) > 2;
11  // --------------------------------------------------------------------------
12  Cut k2018PIDs(double bdtCut, double remidCut, double cvnCut)
13  {
14  std::cout << "ERROR::k2018PIDs. Looking for cvnProd3Train. Branch no longer exists." << std::endl;
15  abort();
16  Cut cut1 = SIMPLEVAR(sel.cosrej.anglekal) > 0.5;
17  return cut1;
18  /*
19  Cut cut2 = SIMPLEVAR(slc.nhit) < 400;
20  Cut cut3 = SIMPLEVAR(sel.nuecosrej.pngptp) < 0.9;
21  Cut cut4 = SIMPLEVAR(sel.cosrej.numucontpid2020) > bdtCut;
22  Cut cut5 = SIMPLEVAR(sel.remid.pid) > remidCut;
23  Cut cut6 = SIMPLEVAR(sel.cvnProd3Train.numuid) > cvnCut;
24 
25  return cut1 && cut2 && cut3 && cut4 && cut5 && cut6;
26  */
27  }
28 
29  // --------------------------------------------------------------------------
30  const Cut kNumuDecafPresel2018([](const caf::SRProxy* sr)
31  {
32  std::cout << "ERROR::kNumuPIDMiniprod5, cutting on both cvnProd3Train and cvn2017."
33  << " Neither exist." << std::endl;
34  abort();
35  return false;
36  /*
37  return (sr->slc.nhit > 20
38  && sr->sel.remid.pid > 0
39  && (sr->sel.cvn2017.numuid > 0 || sr->sel.cvnProd3Train.numuid > 0) );
40  */
41  }
42  );
43 
44  // --------------------------------------------------------------------------
46  [](const caf::SRProxy* sr)
47  { return ( //sr->sel.cosrej.numucontpid > 0.48 || sr->sel.cosrej.numucontpidoldcvn > 0.48 );
48  sr->sel.cosrej.numucontpid2020 > 0.48 || sr->sel.cosrej.numucontpid2019 > 0.48 );
49  }
50  );
51 
52  // --------------------------------------------------------------------------
53  const Cut kNumuCosmicRej2018(
54  [](const caf::SRProxy* sr)
55  { return ( sr->sel.cosrej.anglekal > 0.5 &&
56  sr->sel.cosrej.numucontpid2019 > 0.53 &&
57  sr->slc.nhit < 400 &&
58  sr->sel.nuecosrej.pngptp < 0.9
59  );
60  }
61  );
62 
63  // --------------------------------------------------------------------------
64  const Cut kNumuPID2018(
65  [](const caf::SRProxy* sr)
66  {
67  std::cout << "ERROR::kNumuPID2018, cutting on both cvnProd3Train and cvn2017."
68  << " Neither branch exists anymore. Returning False." << std::endl;
69  abort();
70  return false;
71  /*
72  return (sr->sel.remid.pid > 0.7
73  && sr->sel.cvnProd3Train.numuid > 0.7
74  && sr->sel.cvn2017.numuid > 0.1);
75  */
76  }
77  );
78 }
Cuts and Vars for the 2020 FD DiF Study.
Definition: vars.h:6
Proxy for caf::StandardRecord.
Definition: SRProxy.h:2126
const Cut kNumuDecafPresel2018([](const caf::SRProxy *sr){std::cout<< "ERROR::kNumuPIDMiniprod5, cutting on both cvnProd3Train and cvn2017."<< " Neither exist."<< std::endl;abort();return false;})
Definition: NumuCuts2018.h:13
caf::Proxy< caf::SRCosRej > cosrej
Definition: SRProxy.h:1252
const Cut kNumuCosmicRej2018([](const caf::SRProxy *sr){return(sr->sel.cosrej.anglekal > 0.5 && sr->sel.cosrej.numucontpid2019 > 0.53 && sr->slc.nhit< 400 && sr->sel.nuecosrej.pngptp< 0.9 );})
Definition: NumuCuts2018.h:19
caf::Proxy< caf::SRNueCosRej > nuecosrej
Definition: SRProxy.h:1265
Cut k2018PIDs(double bdtCut, double remidCut, double cvnCut)
caf::Proxy< unsigned int > nhit
Definition: SRProxy.h:1315
const Cut kLooseKBC
caf::StandardRecord * sr
caf::Proxy< float > numucontpid2020
Definition: SRProxy.h:884
OStream cout
Definition: OStream.cxx:6
caf::Proxy< caf::SRSlice > slc
Definition: SRProxy.h:2142
const Cut kNumuPID2018([](const caf::SRProxy *sr){std::cout<< "ERROR::kNumuPID2018, cutting on both cvnProd3Train and cvn2017."<< " Neither branch exists anymore. Returning False."<< std::endl;abort();return false;})
Definition: NumuCuts2018.h:22
#define SIMPLEVAR(CAFNAME)
For Vars where literally all you need is a single CAF variable.
Definition: Var.h:88
caf::Proxy< float > anglekal
Definition: SRProxy.h:859
caf::Proxy< float > pngptp
Definition: SRProxy.h:1061
const Cut kNumuCosmicRejLoose2018([](const caf::SRProxy *sr){return( sr->sel.cosrej.numucontpid2020 > 0.48||sr->sel.cosrej.numucontpid2019 > 0.48);})
Definition: NumuCuts2018.h:16
caf::Proxy< caf::SRIDBranch > sel
Definition: SRProxy.h:2141
caf::Proxy< float > numucontpid2019
Definition: SRProxy.h:883
const Cut kLooseKFC
Definition: NumuCuts2018.cxx:9