AnalyzeSysts.C

Go to the documentation of this file.

//////////////////////////////////////////////////////////////
// AnalyzeSysts.C
// Jeremy Hewes (jhewes15@fnal.gov)
//
// Make plots from SystMakers for 2020 FHC NC analysis
//////////////////////////////////////////////////////////////

#include "TFile.h"

#include "CAFAna/Core/LoadFromFile.h"
#include "CAFAna/Analysis/Exposures.h"
#include "NuXAna/Analysis/SystMaker.h"

using namespace ana;

// -------------------------------------------------------------------------------------
void AnalyzeSysts(std::string ndFile="/nova/data/users/jhewes15/nus20/systmakers/systsmaker_nus_fhc_neardet.root",
                  std::string fdFile="/nova/data/users/jhewes15/nus20/systmakers/systsmaker_nus_fhc_fardet.root") {

  DontAddDirectory guard;

  std::string ofName = "syst_ana.root";

  TFile* inFile  = TFile::Open(ndFile.c_str(), "read");
  TFile* outFile = TFile::Open(ofName.c_str(), "recreate");

  // Get the syst maker from the file
  SystematicsMaker* systsMaker = LoadFromFile<SystematicsMaker>(inFile, "nus_fhc_neardet").release();

  for (auto systMaker : systsMaker->GetSystMakers()) {
    std::cout << "Name is " << systMaker.second->GetName() << std::endl;
  }

  // Run the syst makers
  systsMaker->Go(kAna2018SensitivityFHCNDPOT);

  // Draw plots
  systsMaker->PrintSysts();
  systsMaker->DrawSysts(outFile->mkdir("nus_fhc_neardet"));

  delete systsMaker;
  delete inFile;

  // Repeat for far detector
  inFile  = TFile::Open(fdFile.c_str(), "read");
  systsMaker = LoadFromFile<SystematicsMaker>(inFile, "nus_fhc_fardet").release();
  systsMaker->Go(kAna2018FHCPOT);
  systsMaker->PrintSysts();
  systsMaker->DrawSysts(outFile->mkdir("nus_fhc_fardet"));

  delete systsMaker;
  delete inFile;
  delete outFile;

}