get_cosmic_spectra.C

Go to the documentation of this file.

#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "3FlavorAna/Cuts/NueCuts2017.h"

using namespace ana;

#include "TFile.h"
#include <iostream>

void get_cosmic_spectra()
{


    std::string fnameData = "prod_restrictedcaf_R17-03-01-prod3reco.k_fd_numi_fhc_full_v1_goodruns";

    SpectrumLoader loaderData(fnameData);
    loaderData.SetSpillCut(kStandardSpillCuts);

    Spectrum data_cosmic       (loaderData, kNue2017Axis, kNue2017FDAllSamples&&kInTimingSideband,  kNoShift);
    //TH1 will be scaled to spill, but the livetime will still be high. Use with care
    Spectrum data_cosmic_scaled(loaderData, kNue2017Axis, kNue2017FDAllSamples&&kInTimingSideband,  kNoShift, kTimingSidebandWeight);

    loaderData.Go();
    
    TFile* fout = new TFile("cosmic_prediction.root", "RECREATE");
    data_cosmic.SaveTo(fout, "cosmic_spect");
    data_cosmic_scaled.SaveTo(fout, "cosmic_spect_weighted");

    delete fout;

// in order to get spectra from *.root file do:

/*  TFile * file = new TFile (predFileName.c_str(), "READ");
    auto data_cosmic = Spectrum::LoadFrom(file, "cosmic_spect");
    TH1* cos = data_cosmic->ToTH1(data_cosmic->Livetime(), kGray+1, kSolid, ana::EExposureType::kLivetime);
    cout<<" bkg prediction:  "<<cos->Integral()<<" with "<<data_cosmic->Livetime()<<" s livetime "<<endl;*/
}