AnalyzeSysts.C File Reference

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

Go to the source code of this file.

Functions
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")

Function Documentation

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"
	)

Definition at line 17 of file AnalyzeSysts.C.

References om::cout, allTimeWatchdog::endl, inFile, ana::kAna2018FHCPOT, ana::kAna2018SensitivityFHCNDPOT, outFile, runNovaSAM::release, and string.

                                                                                                                {

  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;

}