NumuEFxs.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/Var.h"

#include "NumuEnergy/func/Numu2017Fits.h"
#include "NumuEnergy/func/Numu2018Fits.h"
#include "NumuEnergy/func/Numu2020Fits.h"

// Energy tunings can be placed here


namespace ana {

  float SAHadEND(double hadVisE);
  float SAHadEFD(double hadVisE, int run);
  float SAMuE(double trkLen);
  float SAGetNDCatEffLen(double trkLenCat);

  double SAMuEND(double lenact, double lencat,
                 double actE,  double tranE,
                 double catE);

  float TAHadEND(double HadVisE);
  float TAMuE(double trklen);
  float TAGetNDActLen(double trkLenact, double trkLencat);


  extern const Var kNumuCCXsecMuonE;
  extern const Var kNumuCCXsecHadE;


  const Var kNumuCCXsecNuE = kNumuCCXsecHadE + kNumuCCXsecMuonE;

  extern const Var kNumuMuE2017FDp1;
  extern const Var kNumuMuE2017FDp2;
  extern const Var kNumuMuE2017FDp3;
  extern const Var kNumuMuE2017FDp4;
  extern const Var kNumuMuE2017FDp5;
 
  extern const Var kNumuHadE2017FDp1;
  extern const Var kNumuHadE2017FDp2;
  extern const Var kNumuHadE2017FDp3;
  extern const Var kNumuHadE2017FDp4;
  extern const Var kNumuHadE2017FDp5;

  const Var kNumuE2017FDp1 = kNumuMuE2017FDp1 + kNumuHadE2017FDp1;
  const Var kNumuE2017FDp2 = kNumuMuE2017FDp2 + kNumuHadE2017FDp2;
  const Var kNumuE2017FDp3 = kNumuMuE2017FDp3 + kNumuHadE2017FDp3;
  const Var kNumuE2017FDp4 = kNumuMuE2017FDp4 + kNumuHadE2017FDp4;
  const Var kNumuE2017FDp5 = kNumuMuE2017FDp5 + kNumuHadE2017FDp5;

  extern const Var kNumuMuE2017ND;
  extern const Var kNumuHadE2017ND;

  const Var kNumuE2017ND = kNumuMuE2017ND + kNumuHadE2017ND;

  extern const Var kNumuMuE2017RHCND;
  extern const Var kNumuHadE2017RHCND;

  const Var kNumuE2017RHCND = kNumuMuE2017RHCND + kNumuHadE2017RHCND;


  extern const Var kNumuMuE2017;
  extern const Var kNumuHadE2017;

  const Var kNumuE2017 = kNumuMuE2017 + kNumuHadE2017;

  /* Numu 2018(aka prod4) Tunes */

  Var initNumuMuE2018FDpXVar(const std::function<NumuEnergyFunc::prod4_fd_energy_p_func_t> &f);
  Var initNumuHadE2018FDpXVar(const std::function<NumuEnergyFunc::prod4_fd_energy_p_func_t> &f);
  Var initNumuMuE2018NDpXVar(const std::function<NumuEnergyFunc::prod4_nd_energy_p_func_t> &f_act,
                             const std::function<NumuEnergyFunc::prod4_nd_energy_p_func_t> &f_cat);
  Var initNumuHadE2018NDpXVar(const std::function<NumuEnergyFunc::prod4_nd_energy_p_func_t> &f);
  Var initNumuMuE2018Var(const std::function<NumuEnergyFunc::prod4_fd_energy_func_t> &f_fd,
                         const std::function<NumuEnergyFunc::prod4_nd_energy_func_t> &f_nd_act,
                         const std::function<NumuEnergyFunc::prod4_nd_energy_func_t> &f_nd_cat);
  Var initNumuHadE2018Var(const std::function<NumuEnergyFunc::prod4_fd_energy_func_t> &f_fd,
                          const std::function<NumuEnergyFunc::prod4_nd_energy_func_t> &f_nd);

  /* Numu 2018(aka prod4) Tunes */

  const Var kNumuMuE2018FDp1 = initNumuMuE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_muon_energy_p<1>
  );

  const Var kNumuMuE2018FDp2 = initNumuMuE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_muon_energy_p<2>
  );

  const Var kNumuMuE2018FDp3 = initNumuMuE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_muon_energy_p<3>
  );

  const Var kNumuMuE2018FDp4 = initNumuMuE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_muon_energy_p<4>
  );

  const Var kNumuMuE2018FDp5 = initNumuMuE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_muon_energy_p<5>
  );

  const Var kNumuHadE2018FDp1 = initNumuHadE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_had_energy_p<1>
  );

  const Var kNumuHadE2018FDp2 = initNumuHadE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_had_energy_p<2>
  );

  const Var kNumuHadE2018FDp3 = initNumuHadE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_had_energy_p<3>
  );

  const Var kNumuHadE2018FDp4 = initNumuHadE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_had_energy_p<4>
  );

  const Var kNumuHadE2018FDp5 = initNumuHadE2018FDpXVar(
    NumuEnergyFunc::predict_prod4_fd_had_energy_p<5>
  );

  const Var kNumuMuE2018ND = initNumuMuE2018NDpXVar(
    NumuEnergyFunc::predict_prod4_nd_act_energy_p<3>,
    NumuEnergyFunc::predict_prod4_nd_cat_energy_p<3>
  );

  const Var kNumuMuE2018RHCND = initNumuMuE2018NDpXVar(
    NumuEnergyFunc::predict_prod4_nd_act_energy_p<4>,
    NumuEnergyFunc::predict_prod4_nd_cat_energy_p<4>
  );

  const Var kNumuHadE2018ND = initNumuHadE2018NDpXVar(
    NumuEnergyFunc::predict_prod4_nd_had_energy_p<3>
  );

  const Var kNumuHadE2018RHCND = initNumuHadE2018NDpXVar(
    NumuEnergyFunc::predict_prod4_nd_had_energy_p<4>
  );

  const Var kNumuMuE2018 = initNumuMuE2018Var(
    NumuEnergyFunc::predict_prod4_fd_muon_energy,
    NumuEnergyFunc::predict_prod4_nd_act_energy,
    NumuEnergyFunc::predict_prod4_nd_cat_energy
  );

  const Var kNumuHadE2018 = initNumuHadE2018Var(
    NumuEnergyFunc::predict_prod4_fd_had_energy,
    NumuEnergyFunc::predict_prod4_nd_had_energy
  );

  const Var kNumuE2018FDp1 = kNumuMuE2018FDp1 + kNumuHadE2018FDp1;
  const Var kNumuE2018FDp2 = kNumuMuE2018FDp2 + kNumuHadE2018FDp2;
  const Var kNumuE2018FDp3 = kNumuMuE2018FDp3 + kNumuHadE2018FDp3;
  const Var kNumuE2018FDp4 = kNumuMuE2018FDp4 + kNumuHadE2018FDp4;
  const Var kNumuE2018FDp5 = kNumuMuE2018FDp5 + kNumuHadE2018FDp5;


  extern const Var kNumuMuE2018ND;
  extern const Var kNumuMuE2018RHCND;
  extern const Var kNumuHadE2018ND;
  extern const Var kNumuHadE2018RHCND;

  const Var kNumuE2018ND    = kNumuMuE2018ND + kNumuHadE2018ND;
  const Var kNumuE2018RHCND = kNumuMuE2018RHCND + kNumuHadE2018RHCND;


  extern const Var kNumuMuE2018;
  extern const Var kNumuHadE2018;

  const Var kNumuE2018 = kNumuMuE2018 + kNumuHadE2018;


    /* Numu 2020 (prod5) Tunes */
    Var initNumuMuE2020FDpXVar(const std::function<NumuEnergyFunc::prod5_fd_energy_p_func_t> &f);
    Var initNumuHadE2020FDpXVar(const std::function<NumuEnergyFunc::prod5_fd_energy_p_func_t> &f);
    Var initNumuMuE2020NDpXVar(const std::function<NumuEnergyFunc::prod5_nd_energy_p_func_t> &f_act,          const std::function<NumuEnergyFunc::prod5_nd_energy_p_func_t> &f_cat);
    Var initNumuHadE2020NDpXVar(const std::function<NumuEnergyFunc::prod5_nd_energy_p_func_t> &f);
    Var initNumuMuE2020Var(const std::function<NumuEnergyFunc::prod5_fd_energy_func_t> &f_fd,        const std::function<NumuEnergyFunc::prod5_nd_energy_func_t> &f_nd_act,const std::function<NumuEnergyFunc::prod5_nd_energy_func_t> &f_nd_cat);
    Var initNumuHadE2020Var(const std::function<NumuEnergyFunc::prod5_fd_energy_func_t> &f_fd,        const std::function<NumuEnergyFunc::prod5_nd_energy_func_t> &f_nd);
        
    const Var kNumuMuE2020FDp1  = initNumuMuE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_muon_energy_p<1>);
    const Var kNumuMuE2020FDp2  = initNumuMuE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_muon_energy_p<2>);
    const Var kNumuMuE2020FDfhc = initNumuMuE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_muon_energy_p<3>);
    const Var kNumuMuE2020FDrhc = initNumuMuE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_muon_energy_p<4>);
        
    const Var kNumuHadE2020FDp1  = initNumuHadE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_had_energy_p<1>);
    const Var kNumuHadE2020FDp2  = initNumuHadE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_had_energy_p<2>);
    const Var kNumuHadE2020FDfhc = initNumuHadE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_had_energy_p<3>);
    const Var kNumuHadE2020FDrhc = initNumuHadE2020FDpXVar(NumuEnergyFunc::predict_prod5_fd_had_energy_p<4>);
        
    const Var kNumuMuE2020FHCND = initNumuMuE2020NDpXVar(NumuEnergyFunc::predict_prod5_nd_act_energy_p<3>,NumuEnergyFunc::predict_prod5_nd_cat_energy_p<3>);
    const Var kNumuMuE2020RHCND = initNumuMuE2020NDpXVar(NumuEnergyFunc::predict_prod5_nd_act_energy_p<4>, NumuEnergyFunc::predict_prod5_nd_cat_energy_p<4>);
    
    const Var kNumuHadE2020FHCND = initNumuHadE2020NDpXVar(NumuEnergyFunc::predict_prod5_nd_had_energy_p<3>);
    const Var kNumuHadE2020RHCND = initNumuHadE2020NDpXVar(NumuEnergyFunc::predict_prod5_nd_had_energy_p<4>);
    
    const Var kNumuMuE2020  = initNumuMuE2020Var (NumuEnergyFunc::predict_prod5_fd_muon_energy,NumuEnergyFunc::predict_prod5_nd_act_energy,NumuEnergyFunc::predict_prod5_nd_cat_energy);
    const Var kNumuHadE2020 = initNumuHadE2020Var(NumuEnergyFunc::predict_prod5_fd_had_energy, NumuEnergyFunc::predict_prod5_nd_had_energy);
    
    const Var kNumuE2020FDp1  = kNumuMuE2020FDp1  + kNumuHadE2020FDp1;
    const Var kNumuE2020FDp2  = kNumuMuE2020FDp2  + kNumuHadE2020FDp2;
    const Var kNumuE2020FDfhc = kNumuMuE2020FDfhc + kNumuHadE2020FDfhc;
    const Var kNumuE2020FDrhc = kNumuMuE2020FDrhc + kNumuHadE2020FDrhc;    
    
    extern const Var kNumuMuE2020FHCND;
    extern const Var kNumuMuE2020RHCND;
    extern const Var kNumuHadE2020FHCND;
    extern const Var kNumuHadE2020RHCND;
    
    const Var kNumuE2020FHCND = kNumuMuE2020FHCND + kNumuHadE2020FHCND;
    const Var kNumuE2020RHCND = kNumuMuE2020RHCND + kNumuHadE2020RHCND;
    
    extern const Var kNumuMuE2020;
    extern const Var kNumuHadE2020;
    
    const Var kNumuE2020 = kNumuMuE2020 + kNumuHadE2020;
    const Var kNumuHadEFrac2020 = kNumuHadE2020 / kNumuE2020;
}