reduce_numu_ana2018_cosmics.C

Go to the documentation of this file.

#include "CAFAna/Cuts/Cuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "3FlavorAna/Cuts/NumuCuts2018.h"
#include "CAFAna/Core/FileReducer.h"
#include "CAFAna/Vars/Vars.h"

#include "StandardRecord/Proxy/SRProxy.h"
#include "StandardRecord/StandardRecord.h"
using namespace ana;

void Clear2DRecoInfo(caf::StandardRecord* sr) // the stuff we never use
{
  sr->trk.kalman.tracks2d.clear();
  sr->trk.kalman.ntracks2d = 0;
  for(caf::SRElastic& elastic: sr->vtx.elastic){
    elastic.fuzzyk.png2d.clear();
    elastic.fuzzyk.npng2d = 0;
  }
  sr->sel.cvnProd3Train.output.clear();
}

void reduce_numu_ana2018_cosmics(std::string fname, std::string outname)
{
  FileReducer reduce(fname, outname);

  reduce.SetMetadata("nova.decaf_skim","numu_loose_contain");

  const Cut kIsND([](const caf::SRProxy* sr)
                  {
                    return sr->hdr.det==caf::kNEARDET;
                  });

  const Cut kIsFD([](const caf::SRProxy* sr)
                  {
                    return sr->hdr.det==caf::kFARDET;
                  });

  const Cut kNumuDecafCut = (kIsND && kNumuDecafCutND2018) || (kIsFD && kCosmicDecafCut2019);

  reduce.AddEventCut(kNumuDecafCut);

  // reduce.AddReductionStep(ClearWindowTrackInfo);
  // reduce.AddReductionStep(ClearDiscreteTrackInfo);
  // reduce.AddReductionStep(Clear2DRecoInfo);

  reduce.Go();
}