AnalyzeNus18Systs.C

Go to the documentation of this file.

//////////////////////////////////////////////////////////////
// AnalyzeNus18Systs.C
// Mike Wallbank (wallbank@fnal.gov)
//
// Analyze all systematics for Nus18 analyses.
//////////////////////////////////////////////////////////////

#include "CAFAna/Analysis/Exposures.h"
#include "NuXAna/Analysis/NusSystsMaker.h"

using namespace ana;

// -------------------------------------------------------------------------------------
void AnalyzeNus18Systs(std::string inFileName, TString opt) {

  TH1::AddDirectory(0);

  // Running options
  bool fhc = opt.Contains("fhc", TString::kIgnoreCase);
  bool rhc = opt.Contains("rhc", TString::kIgnoreCase);
  assert(fhc or rhc);

  // Open the input file and an output file
  std::string fileName = inFileName.substr(0, inFileName.find(".root"));
  std::string saveName = fileName + "Ana.root";
  TFile* inFile = new TFile(inFileName.c_str(), "READ");
  TFile* outFile = new TFile(saveName.c_str(), "RECREATE");

  // Get the syst maker from the file
  NusSystematicsMaker* nus_systs = NusSystematicsMaker::LoadFrom(inFile, "Nus18Systs").release();
  //NusSystematics* nus_systs = NusSystematics::LoadFrom(inFile, "Nus18Systs");

  // Get POT
  double pot = fhc ? kAna2018FHCPOT : kAna2018RHCPOT;

  // Run the syst makers
  nus_systs->Go(kAna2018RHCPOT);

  // Print the systematics
  nus_systs->PrintSysts();

  // Draw the shifted systematics
  nus_systs->DrawSysts(outFile);

  // Get the ISysts
  std::vector<NusISyst*> systs = nus_systs->MakeISysts();
  for (const auto& syst : systs)
    std::cout << "Syst " << syst->ShortName() << std::endl;

  TFile* systfile = TFile::Open("isysts.root", "recreate");
  for (const auto& syst : systs)
    syst->SaveTo(systfile, syst->ShortName().c_str());

  return;

}