demo3.C File Reference

#include "OscLib/OscCalcPMNSOpt.h"
#include "CAFAna/Core/Binning.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Experiment/SingleSampleExperiment.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/Vars.h"
#include "StandardRecord/Proxy/SRProxy.h"
#include "TCanvas.h"
#include "TH2.h"

Go to the source code of this file.

Functions
void	demo3 ()

Function Documentation

void demo3

(

)

Definition at line 20 of file demo3.C.

References std::asin(), ana::bins, calc, e, plot_validation_datamc::fname, ana::SpectrumLoader::Go(), make_syst_table_plots::h, ana::kIsNumuCC(), ana::kNHit, kRed, demo0::loader, demo1::nhit, caf::Proxy< caf::SRSlice >::nhit, ana::OscillatableSpectrum::Oscillated(), pot, osc::_OscCalcPMNSOpt< T >::SetdCP(), osc::_OscCalcPMNSOpt< T >::SetDmsq21(), osc::_OscCalcPMNSOpt< T >::SetDmsq32(), osc::_OscCalcPMNSOpt< T >::SetL(), osc::_OscCalcPMNSOpt< T >::SetRho(), osc::_OscCalcPMNSOpt< T >::SetTh12(), osc::_OscCalcPMNSOpt< T >::SetTh13(), osc::_OscCalcPMNSOpt< T >::SetTh23(), ana::Binning::Simple(), caf::Proxy< caf::StandardRecord >::slc, demo2::sosc, demo3::soscd, std::sqrt(), sr, string, ana::Spectrum::ToTH1(), and ana::ReweightableSpectrum::ToTH2().

{
  const std::string fname = "prod_sumdecaf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1_numu2020";

  SpectrumLoader loader(fname);

  const Binning bins = Binning::Simple(100, 0, 1000);

  const Var kNHit([](const caf::SRProxy* sr)
                  {
                    return sr->slc.nhit;
                  });

  Spectrum nhit("Number of hits in slice", bins, loader, kNHit, kIsNumuCC);

  // With true energy axis too
  OscillatableSpectrum sosc("Number of hits in slice", bins, loader, kNHit, kIsNumuCC);

  // Do it!
  loader.Go();

  // How to scale histograms
  const double pot = 18e20;

  // 2D histograms
  new TCanvas;
  sosc.ToTH2(pot)->Draw("colz");

  // Make a calculator. This is the fastest variant
  osc::OscCalcPMNSOpt calc;

  calc.SetL(810);
  calc.SetRho(2.75);
  calc.SetDmsq21(7.6e-5);
  calc.SetDmsq32(2.35e-3);
  calc.SetTh12(asin(sqrt(.87))/2);
  calc.SetTh13(asin(sqrt(.10))/2);
  calc.SetTh23(TMath::Pi()/4);
  calc.SetdCP(0);

  // Can oscillate a spectrum (numu->numu)
  Spectrum soscd = sosc.Oscillated(&calc, 14, 14);
  TH1* h = soscd.ToTH1(pot);
  h->SetLineColor(kRed);
  new TCanvas;
  nhit.ToTH1(pot)->Draw("hist");
  h->Draw("hist same");
}