make_training_sample.C

Go to the documentation of this file.

/*

 * Extract features for BDT training.

 * 

 * To run through only 10 files, do

 * $ cafe -l 10 make_training_sample.C

 */

#ifdef __CINT__

void make_training_sample(int stride = 5000, int offset = 0)

{

  std::cout << "Sorry, you must run in compiled mode" << std::endl;

}

#else

// CAFAna includes

#include "CAFAna/Core/EventList.h"

#include "CAFAna/Core/MultiVar.h"

//#include "3FlavorAna/Cuts/NumuCuts2018.h"

#include "CAFAna/Cuts/SpillCuts.h"

#include "CAFAna/Vars/Vars.h"




///////NCPi0 Headers//////
#include "NDAna/ncpi0_semi_inc_png_cvn/ncpi0Bins.h"
#include "NDAna/ncpi0_semi_inc_png_cvn/ncpi0Cuts.h"
#include "NDAna/ncpi0_semi_inc_png_cvn/ncpi0Vars.h"
/////////////////////////


using namespace ana;

// Define my own Vars and Cuts

const Var kRun    = SIMPLEVAR(hdr.run);

const Var kSubRun = SIMPLEVAR(hdr.subrun);

const Var kCycle  = SIMPLEVAR(hdr.cycle);

const Var kEvt    = SIMPLEVAR(hdr.evt);

const Var kSlice  = SIMPLEVAR(hdr.subevt);



void make_training_sample(int stride = 5000, int offset = 0)
//void make_training_sample(int stride = 496, int offset = 0)

{

  //string dataset =  "defname: prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1 with stride "+to_string(stride)+" offset "+to_string(offset);
  string dataset =  "defname: CVN_Prong_Prod4_ND_NuMI_MC_2view_Genielike_FlatFlux_CAFs with stride "+to_string(stride)+" offset "+to_string(offset);
  

  MakeTextListFile(dataset, {kIsPreSel}, {"nd_fhc_nc_selection_stride_nominaldataset"+to_string(stride)+"_offset"+to_string(offset)+".txt"},{&kRun, &kSubRun, &kCycle, &kEvt, &kSlice, &kNpCVN1gamma, &kNpCVN2gamma, &kNCGapPng1, &kNCGapPng2, &kNCCurrent, &kPdgPi0},0);

  MakeTextListFile(dataset, {kIsPreSel}, {"nd_fhc_nc_selection_stride_2viewdataset"+to_string(stride)+"_offset"+to_string(offset)+".txt"},{&kRun, &kSubRun, &kCycle, &kEvt, &kSlice, &k2ViewCVN1gamma, &k2ViewCVN2gamma, &kNCGapPng1, &kNCGapPng2, &kNCCurrent, &kPdgPi0},0);

  MakeTextListFile(dataset, {kIsPreSel}, {"nd_fhc_nc_selection_stride_FlatFluxdataset"+to_string(stride)+"_offset"+to_string(offset)+".txt"},{&kRun, &kSubRun, &kCycle, &kEvt, &kSlice, &k2FlatFluxCVN1gamma, &k2FlatFluxCVN2gamma, &kNCGapPng1, &kNCGapPng2, &kNCCurrent, &kPdgPi0},0);

  MakeTextListFile(dataset, {kIsPreSel}, {"nd_fhc_nc_selection_stride_GenieLikedataset"+to_string(stride)+"_offset"+to_string(offset)+".txt"},{&kRun, &kSubRun, &kCycle, &kEvt, &kSlice, &k2GenieLikeCVN1gamma, &k2GenieLikeCVN2gamma, &kNCGapPng1, &kNCGapPng2, &kNCCurrent, &kPdgPi0},0);


}

#endif