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.

Namespaces

 ana
 Cuts and Vars for the 2020 FD DiF Study.
 

Functions

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)
 

Variables

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)