dataprocess_numuccinc.C

Go to the documentation of this file.

#include "CAFAna/Vars/Vars.h"
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/SpectrumLoaderBase.h"
#include "3FlavorAna/Vars/NumuVars.h"

#include "NDAna/numucc_inc/NumuCCIncAnalysis.h"
#include "NDAna/numucc_inc/NumuCCIncBins.h"
#include "NDAna/numucc_inc/NumuCCIncVars.h"
#include "NDAna/numucc_inc/NumuCCIncCuts.h"
#include "NDAna/numucc_inc/NDXSecMuonPID.h"

#include "TFile.h"

#include <iostream>
#include <string>

#include "Utilities/func/MathUtil.h"

using namespace ana;

void dataprocess_numuccinc(const char* fOutname){

  //This script generates a data spectrum with the variables and binning of the NuMu CC Inclusive analysis. 
  std::cout<<"=> dataprocess_numucc_inc_xs() start "<<std::endl;

  //test run: cafe --numuccinc -s 100000 dataprocess_numuccinc.C test.root
  
  // FHC CAF prod4 ND respin 
  const std::string dataset = "defname: prod_caf_R17-11-14-prod4reco.remid-hotfix.b_nd_numi_fhc_full_v1_addShortSimpleCVN_goodruns";
  const std::string dataset_name = "Data";
  
  SpectrumLoader loader(dataset);
  loader.SetSpillCut(kStandardSpillCuts);
  
  Spectrum sDt(loader, kRecoMuKEVsCosVsEnuStandardAxis, kAllNumuCCCuts);
  loader.Go();
  
  //The structure of the file is Data/data to be compatible with the MC systematic ones. 
  TFile fOut(fOutname, "RECREATE"); 
  fOut.cd();
  TDirectory* dir = fOut.mkdir(Form("%s",dataset_name.c_str()));
  dir->cd();
  sDt.SaveTo(dir->mkdir("data"));
  
  std::cout<<"=> dataprocess_numucc_inc_xs() end "<<std::endl;
  
}