Namespaces | Functions | Variables
NumuCCIncVars.h File Reference
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/NumuEFxs.h"
#include "CAFAna/Core/HistAxis.h"
#include "CAFAna/Core/Utilities.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Vars/TruthVars.h"
#include "Utilities/func/MathUtil.h"
#include "NDAna/numucc_inc/NumuCCIncBins.h"
#include "NDAna/numucc_inc/NDXSecMuonPID.h"
#include "TH2.h"
#include "TMath.h"
#include "TH1.h"
#include "TAxis.h"
#include "TVector3.h"
#include "TParticlePDG.h"
#include "TDatabasePDG.h"
#include "TLorentzVector.h"

Go to the source code of this file.


 Cuts and Vars for the 2020 FD DiF Study.


float ana::MuonLLR (const float dedxll, const float scatll)
float ana::MuonEActandCat (double trklenactandcat)
float ana::MuonEAct (double trklenact)
float ana::MuonECat (double trklencat)
float ana::VisibleHadE (float vishadE)
double ana::MuonMass ()
const Var ana::kMuonMass ([](const caf::SRProxy *sr){return MuonMass();})
const Var ana::kMuonMass2 ([](const caf::SRProxy *sr){return util::sqr(MuonMass());})
double ana::NeutronMass ()
const Var ana::kNeutronMass ([](const caf::SRProxy *sr){return NeutronMass();})
const Var ana::kNeutronMass2 ([](const caf::SRProxy *sr){return util::sqr(NeutronMass());})
const HistAxis ana::kRecoMuKEVsCosStandardAxis ("Reconstructed T_{#mu} vs cos #{theta};T_{#mu} [GeV];cos #{theta}", angvsmukebins, kRecoMuKEVsCos)
const NuTruthHistAxis ana::kTrueMuKEVsCosStandardAxisST ("True T_{#mu} vs cos #{theta};T_{#mu} [GeV];cos #{theta}", angvsmukebins, kTrueMuKEVsCosST)
const HistAxis ana::kRecoEStandardAxis ("Reconstructed Neutrino Energy (GeV)", enubins, kRecoE)
const NuTruthHistAxis ana::kTrueEStandardAxisST ("True Neutrino Energy (GeV)", enubins, kTrueEST)
const NuTruthHistAxis ana::kTrueMuKEVsCosVsEnuStandardAxisST ("True T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevsebins, kTrueMuKEVsCosVsEnuST)
const NuTruthHistAxis ana::kTrueMuKEVsCosVsEavailStandardAxisST ("True T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevseavailbins, kTrueMuKEVsCosVsEavailST)
const HistAxis ana::kRecoMuKEVsCosVsEnuStandardAxis ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevsebins, kRecoMuKEVsCosVsEnu)
const HistAxis ana::kRecoMuKEVsCosVsEavailStandardAxis ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevseavailbins, kRecoMuKEVsCosVsEavail)
const HistAxis ana::kRecoMuKEVsCosVsEnuStandardAxis_onlyMuKEUp ("Reco. T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevsebins, kRecoMuKEVsCosVsEnu_onlyMuKEUp)
const HistAxis ana::kRecoMuKEVsCosVsEnuStandardAxis_onlyMuKEDw ("Reco. T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevsebins, kRecoMuKEVsCosVsEnu_onlyMuKEDw)
const HistAxis ana::kRecoMuKEVsCosVsEnuStandardAxisUp ("Reco. T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevsebins, kRecoMuKEVsCosVsEnuUp)
const HistAxis ana::kRecoMuKEVsCosVsEnuStandardAxisDw ("Reco. T_{#mu} vs cos #{theta} vs Neutrino Energy (GeV)", angvsmukevsebins, kRecoMuKEVsCosVsEnuDw)
const HistAxis ana::kRecoMuKEVsCosVsEavailStandardAxisUp ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevseavailbins, kRecoMuKEVsCosVsEavailUp)
const HistAxis ana::kRecoMuKEVsCosVsEavailStandardAxisDw ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevseavailbins, kRecoMuKEVsCosVsEavailDw)
const HistAxis ana::kRecoMuKEVsCosVsEavailStandardAxisAngleUp ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevseavailbins, kRecoMuKEVsCosVsEavailAngleUp)
const HistAxis ana::kRecoMuKEVsCosVsEavailStandardAxisAngleDw ("Reco. T_{#mu} vs cos #{theta} vs Available Energy (GeV)", angvsmukevseavailbins, kRecoMuKEVsCosVsEavailAngleDw)
const NuTruthHistAxis ana::kTrueQ2StandardAxisST ("True Q2 (GeV)", q2bins, kTrueQ2_NT)
const HistAxis ana::kRecoQ2StandardAxis ("Reco Q2 (GeV)", q2bins, kRecoq2)
const HistAxis ana::kRecoEStandardAxis_MuonEUp ("Reconstructed Neutrino Energy (GeV)", enubins, kRecoEUp)
const HistAxis ana::kRecoEStandardAxis_MuonEDw ("Reconstructed Neutrino Energy (GeV)", enubins, kRecoEDw)
const HistAxis ana::kRecoQ2StandardAxis_MuonEUp ("Reco Q2 (GeV)", q2bins, kRecoq2_MuonEUp)
const HistAxis ana::kRecoQ2StandardAxis_MuonEDw ("Reco Q2 (GeV)", q2bins, kRecoq2_MuonEDw)
const HistAxis ana::kRecoQ2StandardAxis_AngleUp ("Reco Q2 (GeV)", q2bins, kRecoq2_AngleUp)
const HistAxis ana::kRecoQ2StandardAxis_AngleDw ("Reco Q2 (GeV)", q2bins, kRecoq2_AngleDw)


const Var ana::kTrkLenActandCat = kTrkLenAct + kTrkLenCat
const Var ana::kIncXsecNuE = kIncXsecMuonE + kIncXsecHadE
const Var ana::kRecoE = kIncXsecNuE
const Var ana::kRecoMuKE = kIncXsecMuonE - kMuonMass
const Var ana::kThetaSystShift = Constant(0.0025)
const Var ana::kRecoMuThetaUp = kRecoMuTheta + kThetaSystShift
const Var ana::kRecoMuThetaDw = kRecoMuTheta - kThetaSystShift
const Var ana::kTrueNuE = VarFromNuTruthVar(kTrueEST)
const Var ana::kTrueMuE = VarFromNuTruthVar(kTrueMuEST)
const Var ana::kTrueMuKE = VarFromNuTruthVar(kTrueMuKEST)
const Var ana::kTrueMuCostheta = VarFromNuTruthVar(kTrueMuCosthetaST)
const Var ana::kRecoHadE = kIncXsecHadE
const NuTruthVar ana::kTrueEavailST = kTrueEavail_NT
const Var ana::kTrueMuKEVsCos = VarFromNuTruthVar(kTrueMuKEVsCosST)
const Var ana::kIncXsecNuEUp = kIncXsecMuonEUp + kIncXsecHadE
const Var ana::kIncXsecNuEDw = kIncXsecMuonEDw + kIncXsecHadE
const Var ana::kRecoEUp = kIncXsecNuEUp
const Var ana::kRecoEDw = kIncXsecNuEDw
const Var ana::kRecoMuKEUp = kIncXsecMuonEUp - kMuonMass
const Var ana::kRecoMuKEDw = kIncXsecMuonEDw - kMuonMass
const HistAxis ana::kTrueMuKEVsCosStandardAxis = HistAxisFromNuTruthHistAxis(kTrueMuKEVsCosStandardAxisST)
const HistAxis ana::kTrueEStandardAxis = HistAxisFromNuTruthHistAxis(kTrueEStandardAxisST)
const HistAxis ana::kTrueMuKEVsCosVsEnuStandardAxis = HistAxisFromNuTruthHistAxis(kTrueMuKEVsCosVsEnuStandardAxisST)
const HistAxis ana::kTrueMuKEVsCosVsEavailStandardAxis = HistAxisFromNuTruthHistAxis(kTrueMuKEVsCosVsEavailStandardAxisST)
const Var ana::kW = Sqrt(kWsqr)
const Var ana::kQ = Sqrt(kQsqr)
const Var ana::kRecoMuP = Sqrt(kRecoMuPsqr)
const Var ana::kRecoMuPUp = Sqrt(kRecoMuPsqrUp)
const Var ana::kRecoMuPDw = Sqrt(kRecoMuPsqrDw)
const Var ana::kResE = (kNumuCCXsecNuE - kTrueNuE)/kTrueNuE
const Var ana::kRecoMuE = kIncXsecMuonE
const Var ana::kRecoMuEUp = kIncXsecMuonEUp
const Var ana::kRecoMuEDw = kIncXsecMuonEDw
const Var ana::kTwo = Constant(2)
const Var ana::kRecoq2 = kTwo * kRecoE *(kRecoMuE- (kRecoMuP *kRecoMuCostheta ))-kMuonMass2
const Var ana::kRecoq2_MuonEUp = kTwo * kRecoEUp*(kRecoMuEUp-(kRecoMuPUp*kRecoMuCostheta ))-kMuonMass2
const Var ana::kRecoq2_MuonEDw = kTwo * kRecoEDw*(kRecoMuEDw-(kRecoMuPDw*kRecoMuCostheta ))-kMuonMass2
const Var ana::kRecoq2_AngleUp = kTwo * kRecoE *(kRecoMuE- (kRecoMuP *kRecoMuCosthetaUp))-kMuonMass2
const Var ana::kRecoq2_AngleDw = kTwo * kRecoE *(kRecoMuE- (kRecoMuP *kRecoMuCosthetaDw))-kMuonMass2
const Var ana::kRecow2 = kNeutronMass2+(kTwo*kNeutronMass*(kRecoE-kRecoMuE))-kRecoq2
const Var ana::kRecoq = Sqrt(kRecoq2)
const Var ana::kRecow = Sqrt(kRecow2)
const Var ana::kDeltatheta = kRecoMuCostheta - kTrueMuCostheta
const Var ana::kDeltamuKE = kRecoMuKE - kTrueMuKE
const HistAxis ana::kTrueQ2StandardAxis = HistAxisFromNuTruthHistAxis(kTrueQ2StandardAxisST)