make_michel.C

Go to the documentation of this file.

// Macro to generate 2018 nue FHC analysis root file for MichelDecomp plots


#include "OscLib/OscCalcDumb.h"

#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/PPFXWeights.h"

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "3FlavorAna/Cuts/NueCuts2018.h"
#include "3FlavorAna/Cuts/NueCuts2018Alt.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Cuts/NueCutsSecondAna.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"

#include "CAFAna/Analysis/SALoaders.h"
#include "3FlavorAna/Decomp/MichelDecomp.h"
#include "3FlavorAna/Decomp/BENDecomp.h"
#include "CAFAna/Decomp/ProportionalDecomp.h"
#include "CAFAna/Decomp/CheatDecomp.h"
#include "3FlavorAna/Decomp/NumuDecomp.h"
#include "CAFAna/Extrap/ModularExtrap.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "CAFAna/Systs/Systs.h"
#include "CAFAna/Core/SystShifts.h"


#include "Utilities/func/MathUtil.h"

#include "TCanvas.h"
#include "TH2.h"
#include "TFile.h"
#include "TLegend.h"
#include "TProfile.h"
#include "TRandom3.h"
#include "3FlavorAna/Vars/NumuVars.h"

#include <iostream>
#include <cmath>

using namespace ana;

void make_michel()
{  
  HistAxis NumuAxis("CCE", Binning::Simple(50,0,5), kCCE);
  const Var Wei = kXSecCVWgt2018 * kPPFXFluxCVWgt;
  
  // FHC
  SpectrumLoader lNDData("prod_caf_R17-09-05-prod4recopreview.f_nd_numi_fhc_full_v1_addShortSimpleCVN_goodruns");
  SpectrumLoader lNDMC("prod_caf_R17-11-14-prod4reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1");
  
  SpectrumLoader lFDnonswap("prod_caf_R17-11-14-prod4reco.d_fd_genie_nonswap_fhc_nova_v08_full_v1");
  SpectrumLoader lFDfluxswap("prod_caf_R17-11-14-prod4reco.d_fd_genie_fluxswap_fhc_nova_v08_full_v1");
  SpectrumLoader lFDtauswap("prod_caf_R17-11-14-prod4reco.d_fd_genie_tau_fhc_nova_v08_full_v1");
  
  lNDData.SetSpillCut(kStandardSpillCuts);
  lNDMC.SetSpillCut(kStandardDQCuts);
  
  lFDnonswap.SetSpillCut(kStandardSpillCuts);
  lFDfluxswap.SetSpillCut(kStandardDQCuts);
  lFDtauswap.SetSpillCut(kStandardSpillCuts);
  
  struct HistDef
  {
      std::string name;
      HistAxis axis;
  };

  const HistDef def = {"cvn2d", {"CVN 2D", kNue2018Binning, kNue2018AnaBin}};
  
  Cut NueND = kNue2018NDPresel && kNue2018CVNCut;

  IDecomp*            cheatDCMP; 
  IDecomp*            michelDCMP;
  IDecomp*            michelbenDCMP; 
  BENDecomp*          benDCMP; 
  ProportionalDecomp* propDCMP;
  
  const HistAxis& axis = def.axis;
  
  // Michel decomp with decomp as input
  cheatDCMP = new CheatDecomp(lNDMC, axis, NueND, kNoShift, Wei);
  michelDCMP = new MichelDecomp(lNDMC, lNDData, axis, NueND, cheatDCMP, kNoShift, kNoShift, Wei);
      
  // COMBO Decomp = Michel+BEN decomp
  benDCMP = new BENDecomp(lNDMC, lNDData, axis, NueND, kNoShift, kNoShift, Wei);
  michelbenDCMP = new MichelDecomp(lNDMC, lNDData, axis, NueND, benDCMP, kNoShift, kNoShift, Wei);
      
  // Proportional decomp
  propDCMP = new ProportionalDecomp(lNDMC, lNDData, axis, NueND, kNoShift, kNoShift, Wei);

  lNDMC.Go();
  lNDData.Go();
  
  lFDnonswap.Go();
  lFDfluxswap.Go();
  lFDtauswap.Go();
  
  TFile* file = new TFile("nue_dcmp.root","recreate");
  TDirectory* d = file->mkdir("nue_decomps");
  
  michelDCMP->SaveTo(d, "michelcheatDCMP");
  michelbenDCMP->SaveTo(d, "michelbenDCMP");
  propDCMP->SaveTo(d, "propDCMP");
      
  file->Close();  
}