get_data_histogram.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/NueCuts2018.h"
#include "CAFAna/Analysis/Exposures.h"

using namespace ana;

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

void get_data_histogram()
{
    /// File for making data histograms for fit
    /// Use it for extracting OLD data hitstograms from decafs made in 2019 release

    bool getfiles=0;
    if (getfiles){

    //std::string fnameData_rhc = "defname: prod4_sumrestricteddecaf_fd_numi_rhc_Ana2018_goodruns_nue2019";
    //std::string fnameData_fhc = "defname: prod4_sumrestricteddecaf_fd_numi_fhc_full_goodruns_nue2018";
    std::string fnameData_rhc = "defname: prod4_sumrestricteddecaf_fd_numi_rhc_full_goodruns_nue2019";
    std::string fnameData_fhc = "defname: prod4_sumrestricteddecaf_fd_numi_fhc_full_goodruns_nue2019";

    SpectrumLoader loaderData_rhc(fnameData_rhc);
    loaderData_rhc.SetSpillCut(kStandardSpillCuts);

    SpectrumLoader loaderData_fhc(fnameData_fhc);
    loaderData_fhc.SetSpillCut(kStandardSpillCuts);

    Spectrum data_rhc       (loaderData_rhc, kNue2018AxisMergedPeripheral, kNue2018FDAllSamples&&kInBeamSpill,  kNoShift);
    Spectrum data_fhc       (loaderData_fhc, kNue2018AxisMergedPeripheral, kNue2018FDAllSamples&&kInBeamSpill,  kNoShift);

    loaderData_rhc.Go();
    loaderData_fhc.Go();
    
    TFile* fout = new TFile("nue_data_2019.root", "RECREATE");
    data_rhc.SaveTo(fout, "spect_rhc");
    data_fhc.SaveTo(fout, "spect_fhc");

    delete fout;
    }
// in order to get spectra from *.root file do:
   if (!getfiles){

    TFile * file = new TFile ("nue_data_2019.root", "READ");
    auto data_rhc = Spectrum::LoadFrom(file, "spect_rhc");
    auto data_fhc = Spectrum::LoadFrom(file, "spect_fhc");
    TH1* hrhc = data_rhc->ToTH1(data_rhc->POT());
    TH1* hfhc = data_fhc->ToTH1(data_fhc->POT());
    cout<<" data rhc:  "<<hrhc->Integral()<<" POT "<<data_rhc->POT()<<endl;
    cout<<" data fhc:  "<<hfhc->Integral()<<" POT "<<data_fhc->POT()<<endl;
   }

}