exporter_fd.C

Go to the documentation of this file.

#ifdef __CINT__
void exporter_fd()
{
  std::cout << "Sorry, you must run in compiled mode" << std::endl;
}
#else

#include <string>
#include <iostream>

#include "CAFAna/numu/SplineEE/exporters/var_defs.h"
#include "CAFAna/Analysis/CSVMaker.h"

#include "3FlavorAna/Cuts/NumuCuts2018.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Cuts/SpillCuts.h"

#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/XsecTunes.h"

using namespace ana;

const std::string data =
  "dataset_def_name_newest_snapshot "
  "prod_caf_R17-11-14-prod4reco.d_fd_genie_nonswap_fhc_nova_v08_period3_v1";

const std::string output = "data_fd_p3.csv";

const Cut cut    = kIsNumuCC && (kTrueE < 5) && (
     kNumuBasicQuality
  && kNumuContainFD2017 && kNumuPID2018 && kNumuCosmicRej2018
);
const Var weight = kPPFXFluxCVWgt * kXSecCVWgt2018;


void exporter_fd()
{
  CSVMaker maker(data, output);

  maker.setCut(cut);
  maker.SetSpillCut(kStandardSpillCuts);

  maker.addVars(VARS_TRUTH);
  maker.addVars(VARS_BASE);
  maker.addVars(VARS_FD);

  maker.addVar("weight", weight);

  maker.Go();
}

#endif