getInputList.C

Go to the documentation of this file.

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NumuCuts2017.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Core/EventList.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "3FlavorAna/Vars/NumuVars.h"

using namespace ana;

void getInputList()
{
  // files
  const std::string fd_MC = "prod_sumdecaf_R17-03-01-prod3reco.l_fd_genie_nonswap_fhc_nova_v08_full_v1_numu_sa_2017";

  // cuts
  const Cut kNumuCosmicRejSA(
                           [](const caf::SRProxy* sr)
                           { return ( sr->sel.cosrej.anglekal  > 0.5 &&
                                      sr->sel.cosrej.numuSAcontpid > 0.535 &&
                                      //sr->sel.cosrej.numucontpid > 0.535 &&
                                      sr->slc.nhit < 400);
                           }
                           );

  const Cut kNumuFDSA = kNumuQuality && kNumuContainFD && kNumuNCRej && kNumuCosmicRejSA;

  // Vars

  const Var kIsSA(
                      [kNumuFDSA](const caf::SRProxy *sr)
                      {
                        return kNumuFDSA(sr);
                      });

  const Var kIsA17(
                      [kNumuCutFD2017](const caf::SRProxy *sr)
                      {
                        return kNumuCutFD2017(sr);
                      });

  const Var kCCESA = SIMPLEVAR(energy.numu.trkccE);

  
  MakeTextListFile(fd_MC, {kNumuFDSA || kNumuCutFD2017}, {"Eventlist.txt"} , {&kTrueE, &kIsSA, &kIsA17, &kHadEFrac, &kCCE, &kCCESA, &kPPFXFluxCVWgt, &kXSecCVWgt2017}, &kStandardSpillCuts);

} // End of function