Functions | Variables
ana::xsec::numubarcc Namespace Reference

Functions

bool VtxInBounds (const caf::SRVector3DProxy *vec, const TVector3 *vmin, const TVector3 *vmax)
 Main analysis truth cuts. More...
 
const Var kBestMuonIDIndex ([](const caf::SRProxy *sr) -> int{int bestidx=sr->trk.kalman.idxmuonid;if(bestidx< 0||bestidx >=(int) sr->trk.kalman.ntracks) return-99;return bestidx;})
 
const Var kMuonID ([](const caf::SRProxy *sr) -> float{int bestidx=kBestMuonIDIndex(sr);if(bestidx< 0||bestidx >=(int) sr->trk.kalman.ntracks) return-99.0;return sr->trk.kalman.tracks[bestidx].muonid;})
 

Variables

const Binning kMuonEActBinning = Binning::Simple(150, 0, 5)
 
const Binning kHadEnergyBinning = Binning::Simple(150, 0, 5)
 
const Binning kTrkLenActBinning = Binning::Simple(160, 0, 16)
 
const Binning kTrkLenCatBinning = Binning::Simple(150, 0, 4)
 
const Binning kMuonECatBinning = Binning::Simple(150, 0, 3)
 
const Binning kVisHadEBinning = Binning::Simple(150., 0.,2.)
 
const Binning kMuonEnergyBinningFine = Binning::Simple(500, 0.0, 5.0)
 
const std::vector< double > mukeedges
 
const Binning muon_ke_bins = Binning::Custom(mukeedges)
 
const std::vector< double > angedges {-1.0, 0.5, 0.56, 0.62, 0.68, 0.74, 0.80, 0.85, 0.88, 0.91, 0.94, 0.96, 0.98, 0.99, 1.0}
 
const Binning angle_bins = Binning::Custom(angedges)
 
const std::vector< double > eavailedges
 
const Binning eavail_bins = Binning::Custom(eavailedges)
 
const bool testo = true
 
const NuTruthCut kIsNumubarCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc &&nu->pdg==-14);})
 
const Cut kIsNumubarCC = CutFromNuTruthCut(kIsNumubarCC_NT)
 
const NuTruthCut kIsNumuCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc &&nu->pdg==+14);})
 
const Cut kIsNumuCC = CutFromNuTruthCut(kIsNumuCC_NT)
 
const NuTruthCut kIsNumuNumubarCC_NT = kIsNumuCC_NT || kIsNumubarCC_NT
 
const Cut kIsNumuNumubarCC = CutFromNuTruthCut(kIsNumuNumubarCC_NT)
 
const NuTruthCut kIsNueCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc &&(nu->pdg==12));})
 Signal cuts. More...
 
const NuTruthCut kIsNuebarCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc &&(nu->pdg==-12));})
 
const NuTruthCut kIsNueorbarCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc &&(std::abs(nu->pdg)==12));})
 
const NuTruthCut kIsCC_NT ([](const caf::SRNeutrinoProxy *nu){return(nu->iscc);})
 
const NuTruthCut kIsNC_NT = !kIsCC_NT
 
const NuTruthCut kTrueMuon_NT ([](const caf::SRNeutrinoProxy *nu){for(int iprim=0;iprim< (int) nu->prim.size();++iprim) if(std::abs(nu->prim[iprim].pdg)==13) return true;return false;})
 
const Cut kIsNueCC = CutFromNuTruthCut(kIsNueCC_NT)
 
const Cut kIsNuebarCC = CutFromNuTruthCut(kIsNuebarCC_NT)
 
const Cut kIsNueorbarCC = CutFromNuTruthCut(kIsNueorbarCC_NT)
 
const Cut kIsCC = CutFromNuTruthCut(kIsCC_NT)
 
const Cut kIsNC = CutFromNuTruthCut(kIsNC_NT)
 
const Cut kTrueMuon = CutFromNuTruthCut(kTrueMuon_NT)
 
const Cut kTrueKalmanMuon ([](const caf::SRProxy *sr){ int ibesttrk=kBestMuonIDIndex(sr);if(ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return false;return std::abs(sr->trk.kalman.tracks[ibesttrk].truth.pdg)==13;})
 
const TVector3 * detector_vtx_min = new TVector3(-191, -187, 0)
 Geometry. More...
 
const TVector3 * detector_vtx_max = new TVector3( 192, 194, 1270)
 
const TVector3 * loose_vtx_min = new TVector3(-160, -160, 25)
 
const TVector3 * loose_vtx_max = new TVector3(160, 160, 1150)
 
const TVector3 * numucc_fiducial_min = new TVector3(-130, -130, 100)
 
const TVector3 * numucc_fiducial_max = new TVector3( 140, 140, 1000)
 
const TVector3 * containLow = new TVector3(-180, -180, 20)
 
const TVector3 * containHigh = new TVector3(180, 180, 1525)
 
const NuTruthCut kTrueVtxDetectorCut_NT ([](const caf::SRNeutrinoProxy *nu){return VtxInBounds(&nu->vtx, detector_vtx_min, detector_vtx_max);})
 
const Cut kTrueVtxDetectorCut = ana::CutFromNuTruthCut(kTrueVtxDetectorCut_NT)
 
const NuTruthCut kTrueVtxLooseCut_NT ([](const caf::SRNeutrinoProxy *nu){return VtxInBounds(&nu->vtx, loose_vtx_min, loose_vtx_max);})
 
const Cut kTrueVtxLooseCut = ana::CutFromNuTruthCut(kTrueVtxLooseCut_NT)
 
const NuTruthCut kTrueVtxCut_NT = kTrueVtxLooseCut_NT
 
const Cut kTrueVtxCut = ana::CutFromNuTruthCut(kTrueVtxCut_NT)
 
const NuTruthCut kTrueVtxNumuFiducialCut_NT ([](const caf::SRNeutrinoProxy *nu){return VtxInBounds(&nu->vtx, numucc_fiducial_min, numucc_fiducial_max);})
 
const Cut kTrueVtxNumuFiducialCut = ana::CutFromNuTruthCut(kTrueVtxNumuFiducialCut_NT)
 
const Cut kQualityCut ([](const caf::SRProxy *sr){return(sr->trk.kalman.ntracks > 0 &&sr->slc.nhit > 20 &&sr->slc.ncontplanes > 4);})
 Preselection. More...
 
const Cut kContainmentCut ([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return false;int ibesttrk=kBestMuonIDIndex(sr);if(ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return false;for(const caf::SRFuzzyKProngProxy &prong:sr->vtx.elastic.fuzzyk.png) if(!VtxInBounds(&prong.shwlid.start, containLow, containHigh)||!VtxInBounds(&prong.shwlid.stop, containLow, containHigh)) return false;if(sr->trk.kalman.ntracks< 1) return false;const unsigned short muon_catcher_edge=1275;for(unsigned int i=0;i< sr->trk.kalman.ntracks;++i){if(int(i)==ibesttrk) continue;if(sr->trk.kalman.tracks[i].start.Z() > muon_catcher_edge|| sr->trk.kalman.tracks[i].stop.Z() > muon_catcher_edge) return false;}const caf::SRKalmanTrackProxy &besttrack=sr->trk.kalman.tracks[ibesttrk];return((besttrack.stop.Z()< muon_catcher_edge||besttrack.trkyposattrans< 55) &&besttrack.trkfwdcellnd > 5 &&besttrack.trkbakcellnd > 10);})
 
const Cut kPreselect = kQualityCut && kContainmentCut
 
const Cut kRecoVtxDetectorCut ([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return false;return VtxInBounds(&sr->vtx.elastic.vtx, detector_vtx_min, detector_vtx_max);})
 Main analysis reconstruction cuts. More...
 
const Cut kRecoVtxLooseCut ([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return false;return VtxInBounds(&sr->vtx.elastic.vtx, loose_vtx_min, loose_vtx_max);})
 
const Cut kRecoVtxNumuFiducialCut ([](const caf::SRProxy *sr){if(!sr->vtx.elastic.IsValid) return false;return VtxInBounds(&sr->vtx.elastic.vtx, numucc_fiducial_min, numucc_fiducial_max);})
 
const Cut kRecoVtxCut = kRecoVtxLooseCut
 
const float kMuonIDCutVal = 0.0
 
const Cut kMuonIDCut ([](const caf::SRProxy *sr){return kMuonID(sr) >=kMuonIDCutVal;})
 
const Cut kActive ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(sr->trk.kalman.ntracks< 1||ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return false;const caf::SRKalmanTrackProxy &bestmuon=sr->trk.kalman.tracks[ibesttrk];return(bestmuon.leninact > 0 && bestmuon.lenincat< 0);})
 
const Cut kActiveAndCatcher ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(sr->trk.kalman.ntracks< 1||ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return false;const caf::SRKalmanTrackProxy &bestmuon=sr->trk.kalman.tracks[ibesttrk];return(bestmuon.leninact > 0 && bestmuon.lenincat > 0);})
 
const TVector3 * contain_low
 
const TVector3 * contain_high
 
const Var kBestMuonIDIndex
 
const Var kMuonID
 
const std::string nominal_prod5 = "nd_rhc_prod5_caf_minus_muonid_training_minus_fakedata"
 Dataset definitions. More...
 
const std::string nominal_prod5_parent = "prod_caf_R19-11-18-prod5reco.d_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_v1"
 
const std::string nominal_prod5_decaf = "nd_rhc_prod5_decaf_minus_muonid_training_minus_fakedata"
 
const std::string nominal_prod5_1 = "prod_caf_R20-11-25-prod5.1reco.a_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_v1"
 
const std::string nominal_prod5_1_decaf = "prod_sumdecaf_development_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_ndphysics_contain_v1"
 
const std::string nominal_prod5_1_flatcaf = "prod_flatsumdecaf_development_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_ndphysics_contain_v1"
 
const std::string nominal_dataset = nominal_prod5_1
 
const std::string nominal_dataset_parent = nominal_prod5_1
 
const std::string nominal_decaf = nominal_prod5_1_decaf
 
const std::string fhc_prod5_dataset = "nd_fhc_prod5_caf_minus_muonid_training_minus_fakedata"
 
const std::string fhc_dataset_parent = "prod_caf_R19-11-18-prod5reco.d.h.l_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"
 
const std::string fhc_prod5_decaf = "nd_fhc_prod5_decaf_minus_muonid_training_minus_fakedata"
 
const std::string fhc_parent_decaf = "prod_decaf_R19-11-18-prod5reco.d.h.l_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"
 
const std::string fhc_prod5_1 = "prod_caf_R20-11-25-prod5.1reco.a_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"
 
const std::string fhc_prod5_1_decaf = "prod_sumdecaf_development_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_ndphysics_contain_v1"
 
const std::string fhc_dataset = fhc_prod5_1
 
const std::string fhc_decaf = fhc_prod5_1_decaf
 
std::map< const std::string, const std::stringsyst_datasets
 
const double MUON_MASS = TDatabasePDG::Instance()->GetParticle(13)->Mass()
 
const Var kLongestTrack ([](const caf::SRProxy *sr){if(sr->trk.kalman.ntracks==0) return 0.f;auto idx=sr->trk.kalman.idxlongest;return float(sr->trk.kalman.tracks[idx].len);})
 
const NuTruthVar kTrueNuE_NT ([](const caf::SRNeutrinoProxy *nu){return nu->E;})
 
const Var kTrueNuE = VarFromNuTruthVar(kTrueNuE_NT)
 
const NuTruthVar kTrueMuonE_NT ([](const caf::SRNeutrinoProxy *nu){if(abs(nu->pdg)!=14||!nu->iscc) return(float)-5.0;for(const caf::SRTrueParticleProxy &primary:nu->prim){if(abs(primary.pdg)==13){return(float) primary.p.T();}}return(float)-5.0;})
 
const Var kTrueMuonE = VarFromNuTruthVar(kTrueMuonE_NT)
 
const NuTruthVar kTrueMuonKE_NT ([](const caf::SRNeutrinoProxy *nu){if(abs(nu->pdg)!=14||!nu->iscc) return-5.0;float muone=kTrueMuonE_NT(nu);if(muone< MUON_MASS) return-5.0;return muone-MUON_MASS;})
 
const Var kTrueMuonKE = VarFromNuTruthVar(kTrueMuonKE_NT)
 
const Var kTrueActiveE ([](const caf::SRProxy *sr){return kTrueMuonE(sr)-kTrueCatcherE(sr);})
 
const Var kTrueCatcherE = SIMPLEVAR(energy.numu.mc.truemuoncatcherE)
 
const Var kTrkLenAct ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(sr->trk.kalman.ntracks< 1||ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return-1000.f;const caf::SRKalmanTrackProxy &bestmuon=sr->trk.kalman.tracks[ibesttrk]; if(bestmuon.leninact > 0 && bestmuon.lenincat<=0) return float((bestmuon.leninact/100.) +(bestmuon.lenincat/100.)); if(bestmuon.leninact > 0 && bestmuon.lenincat > 0) return float(bestmuon.leninact/100.); if(bestmuon.leninact<=0 && bestmuon.lenincat > 0) return 0.f;return-1000.f;})
 
const Var kTrkLenCat ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(sr->trk.kalman.ntracks< 1||ibesttrk< 0||ibesttrk >=int(sr->trk.kalman.ntracks)) return-1000.f;const caf::SRKalmanTrackProxy &bestmuon=sr->trk.kalman.tracks[ibesttrk]; if(bestmuon.leninact > 0 && bestmuon.lenincat<=0) return 0.f; if(bestmuon.leninact > 0 && bestmuon.lenincat > 0) return float(bestmuon.lenincat/100.); if(bestmuon.leninact<=0 && bestmuon.lenincat > 0) return float((bestmuon.leninact/100.) +(bestmuon.lenincat/100.));return-1000.f;})
 
const Var kRecoMuonE ([](const caf::SRProxy *sr){double trklenact=kTrkLenAct(sr)*100.;double trklencat=kTrkLenCat(sr)*100.;if(trklenact > 0 &&trklencat<=0) return NumuEnergyFunc::predict_prod5_nd_act_energy(trklenact, sr->spill.isRHC);else if(trklenact<=0 &&trklencat > 0) return NumuEnergyFunc::predict_prod5_nd_cat_energy(trklencat, sr->spill.isRHC);else if(trklenact > 0 &&trklencat > 0) return NumuEnergyFunc::predict_prod5_nd_act_energy(trklenact, sr->spill.isRHC)+NumuEnergyFunc::predict_prod5_nd_cat_energy(trklencat, sr->spill.isRHC);return-1000.;})
 
const Var kRecoMuonKE ([](const caf::SRProxy *sr){return kRecoMuonE(sr)-MUON_MASS;})
 
const NuTruthVar kTrueMuCosTheta_NT ([](const caf::SRNeutrinoProxy *nu){if(abs(nu->pdg)!=14||!nu->iscc) return-5.0;TVector3 beamdir=NuMIBeamDirection(caf::kNEARDET);for(const caf::SRTrueParticleProxy &primary:nu->prim){if(abs(primary.pdg)==13){TVector3 mudir=primary.p.Vect();return mudir.Unit().Dot(beamdir.Unit());}}return-5.0;})
 
const Var kTrueMuCosTheta = VarFromNuTruthVar(kTrueMuCosTheta_NT)
 
const Var kRecoMuCosTheta ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(ibesttrk< 0||ibesttrk >=(int) sr->trk.kalman.ntracks) return-99.0;TVector3 dir=sr->trk.kalman.tracks[ibesttrk].dir;TVector3 beamdir=NuMIBeamDirection(sr->hdr.det);return dir.Dot(beamdir);})
 
const Var kRecoMuTheta ([](const caf::SRProxy *sr){int ibesttrk=kBestMuonIDIndex(sr);if(ibesttrk< 0||ibesttrk >=(int) sr->trk.kalman.ntracks) return-99.0;TVector3 dir=sr->trk.kalman.tracks[ibesttrk].dir;TVector3 beamdir=NuMIBeamDirection(sr->hdr.det);return dir.Angle(beamdir);})
 
const Var kRecoMuCosThetaUp ([](const caf::SRProxy *sr){return(float) std::cos(kRecoMuThetaUp(sr));})
 
const Var kRecoMuCosThetaDw ([](const caf::SRProxy *sr){return(float) std::cos(kRecoMuThetaDw(sr));})
 
const Var kRecoEAvail = kRecoEAvailRHC
 
const Var kRecoEAvailFHC ([](const caf::SRProxy *sr){return 1.68 *ana::kNumuHadVisE(sr)+0.02 *pow(ana::kNumuHadVisE(sr), 2);})
 
const Var kRecoEAvailRHC ([](const caf::SRProxy *sr){return 1.796 *ana::kNumuHadVisE(sr)-0.140;})
 
const HistAxis kTrueRecoMuonKEAxis ("reco_muon_energy", ana::xsec::numubarcc::kMuonEnergyBinningFine, kRecoMuonKE,"true_muon_energy", ana::xsec::numubarcc::kMuonEnergyBinningFine, kTrueMuonKE)
 
const ana::NuTruthVar xsec_wgt_ST = ana::kXSecCVWgt2020_NT
 
const ana::Var xsec_wgt = ana::VarFromNuTruthVar(xsec_wgt_ST, 1)
 
const ana::NuTruthVar flux_wgt_ST = ana::kPPFXFluxCVWgtST
 
const ana::Var flux_wgt = ana::VarFromNuTruthVar(flux_wgt_ST, 1)
 
const ana::NuTruthVar std_wgt_ST = flux_wgt_ST*xsec_wgt_ST
 
const ana::Var std_wgt = ana::VarFromNuTruthVar(std_wgt_ST, 1)
 
const ana::Var kRecoMuonActEnergy
 
const ana::Var kRecoMuonCatEnergy
 
const ana::Var kThetaSystShift = Constant(0.0)
 
const ana::Var kRecoMuThetaUp = kRecoMuTheta + kThetaSystShift
 
const ana::Var kRecoMuThetaDw = kRecoMuTheta - kThetaSystShift
 
const ana::NuTruthVar kTrueEAvail_NT = ana::kTrueEavail_NT
 
const ana::Var kTrueEAvail = ana::kTrueEavail
 

Function Documentation

const Var ana::xsec::numubarcc::kBestMuonIDIndex ( [] (const caf::SRProxy *sr) -> int{int bestidx=sr->trk.kalman.idxmuonid;if(bestidx< 0||bestidx >=(int) sr->trk.kalman.ntracks) return-99;return bestidx;}  )
const Var ana::xsec::numubarcc::kMuonID ( [] (const caf::SRProxy *sr) -> float{int bestidx=kBestMuonIDIndex(sr);if(bestidx< 0||bestidx >=(int) sr->trk.kalman.ntracks) return-99.0;return sr->trk.kalman.tracks[bestidx].muonid;}  )
bool ana::xsec::numubarcc::VtxInBounds ( const caf::SRVector3DProxy vec,
const TVector3 *  vmin,
const TVector3 *  vmax 
)

Main analysis truth cuts.

Definition at line 12 of file NumubarCCIncCuts.cxx.

References caf::Proxy< caf::SRNeutrino >::iscc, kIsNumubarCC_NT, caf::Proxy< caf::SRNeutrino >::pdg, caf::Proxy< caf::SRVector3D >::X(), caf::Proxy< caf::SRVector3D >::Y(), and caf::Proxy< caf::SRVector3D >::Z().

12  {
13  return (vec->X() >= vmin->X() && vec->Y() >= vmin->Y() && vec->Z() >= vmin->Z() &&
14  vec->X() <= vmax->X() && vec->Y() <= vmax->Y() && vec->Z() <= vmax->Z());
15  }

Variable Documentation

const std::vector< double > ana::xsec::numubarcc::angedges {-1.0, 0.5, 0.56, 0.62, 0.68, 0.74, 0.80, 0.85, 0.88, 0.91, 0.94, 0.96, 0.98, 0.99, 1.0}

Definition at line 25 of file NumubarCCIncBins.cxx.

const Binning ana::xsec::numubarcc::angle_bins = Binning::Custom(angedges)

Definition at line 26 of file NumubarCCIncBins.cxx.

const TVector3* ana::xsec::numubarcc::contain_high
const TVector3* ana::xsec::numubarcc::contain_low
const TVector3* ana::xsec::numubarcc::containHigh = new TVector3(180, 180, 1525)

Definition at line 92 of file NumubarCCIncCuts.cxx.

const TVector3* ana::xsec::numubarcc::containLow = new TVector3(-180, -180, 20)

Definition at line 91 of file NumubarCCIncCuts.cxx.

const TVector3* ana::xsec::numubarcc::detector_vtx_max = new TVector3( 192, 194, 1270)

Definition at line 82 of file NumubarCCIncCuts.cxx.

const TVector3 * ana::xsec::numubarcc::detector_vtx_min = new TVector3(-191, -187, 0)

Geometry.

Definition at line 81 of file NumubarCCIncCuts.cxx.

const Binning ana::xsec::numubarcc::eavail_bins = Binning::Custom(eavailedges)

Definition at line 33 of file NumubarCCIncBins.cxx.

const std::vector< double > ana::xsec::numubarcc::eavailedges
Initial value:
{
0.0, 0.25, 0.50, 0.75,
1.0, 1.25, 1.50, 1.75,
2.0, 2.25, 2.50, 120}

Definition at line 29 of file NumubarCCIncBins.cxx.

const std::string ana::xsec::numubarcc::fhc_dataset = fhc_prod5_1

Definition at line 26 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_dataset_parent = "prod_caf_R19-11-18-prod5reco.d.h.l_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"

Definition at line 19 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_decaf = fhc_prod5_1_decaf

Definition at line 27 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_parent_decaf = "prod_decaf_R19-11-18-prod5reco.d.h.l_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"

Definition at line 21 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_prod5_1 = "prod_caf_R20-11-25-prod5.1reco.a_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_v1"

Definition at line 23 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_prod5_1_decaf = "prod_sumdecaf_development_nd_genie_N1810j0211a_nonswap_fhc_nova_v08_full_ndphysics_contain_v1"

Definition at line 24 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_prod5_dataset = "nd_fhc_prod5_caf_minus_muonid_training_minus_fakedata"

Definition at line 18 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::fhc_prod5_decaf = "nd_fhc_prod5_decaf_minus_muonid_training_minus_fakedata"

Definition at line 20 of file NumubarCCIncSystDefs.cxx.

const ana::Var ana::xsec::numubarcc::flux_wgt = ana::VarFromNuTruthVar(flux_wgt_ST, 1)

Definition at line 28 of file NumubarCCIncVars.h.

Referenced by muonid_fhc_dist(), and muonid_rhc_dist().

const ana::NuTruthVar ana::xsec::numubarcc::flux_wgt_ST = ana::kPPFXFluxCVWgtST

Definition at line 27 of file NumubarCCIncVars.h.

const Cut ana::xsec::numubarcc::kActive

Definition at line 132 of file NumubarCCIncCuts.h.

Referenced by estimate_energy().

const Cut ana::xsec::numubarcc::kActiveAndCatcher

Definition at line 133 of file NumubarCCIncCuts.h.

Referenced by estimate_energy(), and GetHistsND().

const ana::Var ana::xsec::numubarcc::kBestMuonIDIndex

Definition at line 54 of file NumubarCCIncVars.h.

const Cut ana::xsec::numubarcc::kContainmentCut
const Binning ana::xsec::numubarcc::kHadEnergyBinning = Binning::Simple(150, 0, 5)

Definition at line 7 of file NumubarCCIncBins.cxx.

const Cut ana::xsec::numubarcc::kIsCC = CutFromNuTruthCut(kIsCC_NT)

Definition at line 58 of file NumubarCCIncCuts.cxx.

Referenced by resolution_event_tree().

const NuTruthCut ana::xsec::numubarcc::kIsCC_NT

Definition at line 74 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kIsNC = CutFromNuTruthCut(kIsNC_NT)

Definition at line 59 of file NumubarCCIncCuts.cxx.

Referenced by muonID(), muonid_fhc_dist(), muonid_rhc_dist(), and purity().

const NuTruthCut ana::xsec::numubarcc::kIsNC_NT = !kIsCC_NT

Definition at line 46 of file NumubarCCIncCuts.cxx.

const Cut ana::xsec::numubarcc::kIsNuebarCC = CutFromNuTruthCut(kIsNuebarCC_NT)

Definition at line 56 of file NumubarCCIncCuts.cxx.

const NuTruthCut ana::xsec::numubarcc::kIsNuebarCC_NT

Definition at line 70 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kIsNueCC = CutFromNuTruthCut(kIsNueCC_NT)

Definition at line 55 of file NumubarCCIncCuts.cxx.

const NuTruthCut ana::xsec::numubarcc::kIsNueCC_NT

Signal cuts.

Signal definitions.

Definition at line 68 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kIsNueorbarCC = CutFromNuTruthCut(kIsNueorbarCC_NT)

Definition at line 57 of file NumubarCCIncCuts.cxx.

Referenced by muonID(), muonid_fhc_dist(), muonid_rhc_dist(), and purity().

const NuTruthCut ana::xsec::numubarcc::kIsNueorbarCC_NT

Definition at line 72 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kIsNumubarCC = CutFromNuTruthCut(kIsNumubarCC_NT)

Definition at line 20 of file NumubarCCIncCuts.cxx.

Referenced by cafana_vars(), muonID(), muonid_fhc_dist(), muonid_rhc_dist(), and purity().

const NuTruthCut ana::xsec::numubarcc::kIsNumubarCC_NT
const Cut ana::xsec::numubarcc::kIsNumuCC = CutFromNuTruthCut(kIsNumuCC_NT)

Definition at line 26 of file NumubarCCIncCuts.cxx.

Referenced by muonID(), muonid_fhc_dist(), muonid_rhc_dist(), and purity().

const NuTruthCut ana::xsec::numubarcc::kIsNumuCC_NT

Definition at line 36 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kIsNumuNumubarCC = CutFromNuTruthCut(kIsNumuNumubarCC_NT)

Definition at line 29 of file NumubarCCIncCuts.cxx.

Referenced by estimate_energy().

const NuTruthCut ana::xsec::numubarcc::kIsNumuNumubarCC_NT = kIsNumuCC_NT || kIsNumubarCC_NT

Definition at line 28 of file NumubarCCIncCuts.cxx.

const Var ana::xsec::numubarcc::kLongestTrack

Definition at line 23 of file NumubarCCIncVars.h.

const Binning ana::xsec::numubarcc::kMuonEActBinning = Binning::Simple(150, 0, 5)

Definition at line 6 of file NumubarCCIncBins.cxx.

Referenced by estimate_energy().

const Binning ana::xsec::numubarcc::kMuonECatBinning = Binning::Simple(150, 0, 3)

Definition at line 10 of file NumubarCCIncBins.cxx.

Referenced by estimate_energy().

const Binning ana::xsec::numubarcc::kMuonEnergyBinningFine = Binning::Simple(500, 0.0, 5.0)

Definition at line 13 of file NumubarCCIncBins.cxx.

Referenced by energy_resolution().

const ana::Var ana::xsec::numubarcc::kMuonID
const Cut ana::xsec::numubarcc::kMuonIDCut

Definition at line 127 of file NumubarCCIncCuts.h.

Referenced by estimate_energy().

const float ana::xsec::numubarcc::kMuonIDCutVal = 0.0

Definition at line 204 of file NumubarCCIncCuts.cxx.

const Cut ana::xsec::numubarcc::kPreselect = kQualityCut && kContainmentCut

Definition at line 151 of file NumubarCCIncCuts.cxx.

Referenced by estimate_energy().

const Cut ana::xsec::numubarcc::kQualityCut

Preselection.

Preselection Cuts.

Definition at line 87 of file NumubarCCIncCuts.h.

Referenced by fiducial_opt(), muonID(), muonid_opt(), muonid_optimization(), purity(), and vertex_fiducial_optimization().

const ana::Var ana::xsec::numubarcc::kRecoEAvail = kRecoEAvailRHC

Definition at line 235 of file NumubarCCIncVars.cxx.

const ana::Var ana::xsec::numubarcc::kRecoEAvailFHC

Definition at line 70 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kRecoEAvailRHC

Definition at line 71 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kRecoMuCosTheta

Definition at line 60 of file NumubarCCIncVars.h.

const Var ana::xsec::numubarcc::kRecoMuCosThetaDw([](const caf::SRProxy *sr){return(float) std::cos(kRecoMuThetaDw(sr));})
const Var ana::xsec::numubarcc::kRecoMuCosThetaUp([](const caf::SRProxy *sr){return(float) std::cos(kRecoMuThetaUp(sr));})
const ana::Var ana::xsec::numubarcc::kRecoMuonActEnergy
const ana::Var ana::xsec::numubarcc::kRecoMuonCatEnergy
const ana::Var ana::xsec::numubarcc::kRecoMuonE

Definition at line 44 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kRecoMuonKE

Definition at line 45 of file NumubarCCIncVars.h.

Referenced by energy_resolution().

const ana::Var ana::xsec::numubarcc::kRecoMuTheta

Definition at line 61 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kRecoMuThetaDw = kRecoMuTheta - kThetaSystShift

Definition at line 64 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kRecoMuThetaUp = kRecoMuTheta + kThetaSystShift

Definition at line 63 of file NumubarCCIncVars.h.

const Cut ana::xsec::numubarcc::kRecoVtxCut = kRecoVtxLooseCut

Definition at line 202 of file NumubarCCIncCuts.cxx.

Referenced by muonID(), muonid_optimization(), and purity().

const Cut ana::xsec::numubarcc::kRecoVtxDetectorCut

Main analysis reconstruction cuts.

Definition at line 118 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kRecoVtxLooseCut

Definition at line 119 of file NumubarCCIncCuts.h.

Referenced by muonid_opt().

const Cut ana::xsec::numubarcc::kRecoVtxNumuFiducialCut

Definition at line 120 of file NumubarCCIncCuts.h.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kThetaSystShift = Constant(0.0)

Definition at line 62 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kTrkLenAct

Definition at line 50 of file NumubarCCIncVars.h.

Referenced by estimate_energy(), and muonID().

const Binning ana::xsec::numubarcc::kTrkLenActBinning = Binning::Simple(160, 0, 16)

Definition at line 8 of file NumubarCCIncBins.cxx.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kTrkLenCat

Definition at line 51 of file NumubarCCIncVars.h.

Referenced by estimate_energy(), and muonID().

const Binning ana::xsec::numubarcc::kTrkLenCatBinning = Binning::Simple(150, 0, 4)

Definition at line 9 of file NumubarCCIncBins.cxx.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kTrueActiveE

Definition at line 41 of file NumubarCCIncVars.h.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kTrueCatcherE = SIMPLEVAR(energy.numu.mc.truemuoncatcherE)

Definition at line 66 of file NumubarCCIncVars.cxx.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kTrueEAvail = ana::kTrueEavail

Definition at line 68 of file NumubarCCIncVars.h.

const ana::NuTruthVar ana::xsec::numubarcc::kTrueEAvail_NT = ana::kTrueEavail_NT

Definition at line 67 of file NumubarCCIncVars.h.

const Cut ana::xsec::numubarcc::kTrueKalmanMuon

Definition at line 82 of file NumubarCCIncCuts.h.

Referenced by estimate_energy().

const ana::Var ana::xsec::numubarcc::kTrueMuCosTheta = VarFromNuTruthVar(kTrueMuCosTheta_NT)

Definition at line 196 of file NumubarCCIncVars.cxx.

const ana::NuTruthVar ana::xsec::numubarcc::kTrueMuCosTheta_NT

Definition at line 58 of file NumubarCCIncVars.h.

const Cut ana::xsec::numubarcc::kTrueMuon = CutFromNuTruthCut(kTrueMuon_NT)

Definition at line 61 of file NumubarCCIncCuts.cxx.

const NuTruthCut ana::xsec::numubarcc::kTrueMuon_NT

Definition at line 79 of file NumubarCCIncCuts.h.

const ana::Var ana::xsec::numubarcc::kTrueMuonE = VarFromNuTruthVar(kTrueMuonE_NT)

Definition at line 47 of file NumubarCCIncVars.cxx.

Referenced by estimate_energy().

const ana::NuTruthVar ana::xsec::numubarcc::kTrueMuonE_NT

Definition at line 33 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kTrueMuonKE = VarFromNuTruthVar(kTrueMuonKE_NT)

Definition at line 60 of file NumubarCCIncVars.cxx.

Referenced by energy_resolution().

const ana::NuTruthVar ana::xsec::numubarcc::kTrueMuonKE_NT

Definition at line 35 of file NumubarCCIncVars.h.

const ana::Var ana::xsec::numubarcc::kTrueNuE = VarFromNuTruthVar(kTrueNuE_NT)

Definition at line 32 of file NumubarCCIncVars.cxx.

const ana::NuTruthVar ana::xsec::numubarcc::kTrueNuE_NT

Definition at line 38 of file NumubarCCIncVars.h.

const ana::HistAxis ana::xsec::numubarcc::kTrueRecoMuonKEAxis

Definition at line 74 of file NumubarCCIncVars.h.

Referenced by energy_resolution().

const Cut ana::xsec::numubarcc::kTrueVtxCut = ana::CutFromNuTruthCut(kTrueVtxCut_NT)

Definition at line 106 of file NumubarCCIncCuts.cxx.

Referenced by muonID(), and purity().

const NuTruthCut ana::xsec::numubarcc::kTrueVtxCut_NT = kTrueVtxLooseCut_NT

Definition at line 105 of file NumubarCCIncCuts.cxx.

Referenced by muonid_optimization().

const Cut ana::xsec::numubarcc::kTrueVtxDetectorCut = ana::CutFromNuTruthCut(kTrueVtxDetectorCut_NT)

Definition at line 98 of file NumubarCCIncCuts.cxx.

const NuTruthCut ana::xsec::numubarcc::kTrueVtxDetectorCut_NT

Definition at line 53 of file NumubarCCIncCuts.h.

const Cut ana::xsec::numubarcc::kTrueVtxLooseCut = ana::CutFromNuTruthCut(kTrueVtxLooseCut_NT)

Definition at line 104 of file NumubarCCIncCuts.cxx.

Referenced by muonid_fhc_dist(), and muonid_rhc_dist().

const NuTruthCut ana::xsec::numubarcc::kTrueVtxLooseCut_NT

Definition at line 56 of file NumubarCCIncCuts.h.

Referenced by fiducial_opt(), muonid_opt(), and vertex_fiducial_optimization().

const Cut ana::xsec::numubarcc::kTrueVtxNumuFiducialCut = ana::CutFromNuTruthCut(kTrueVtxNumuFiducialCut_NT)

Definition at line 112 of file NumubarCCIncCuts.cxx.

Referenced by estimate_energy().

const NuTruthCut ana::xsec::numubarcc::kTrueVtxNumuFiducialCut_NT

Definition at line 62 of file NumubarCCIncCuts.h.

const Binning ana::xsec::numubarcc::kVisHadEBinning = Binning::Simple(150., 0.,2.)

Definition at line 11 of file NumubarCCIncBins.cxx.

const TVector3 * ana::xsec::numubarcc::loose_vtx_max = new TVector3(160, 160, 1150)

Definition at line 85 of file NumubarCCIncCuts.cxx.

Referenced by vertex_fiducial_optimization_plots().

const TVector3 * ana::xsec::numubarcc::loose_vtx_min = new TVector3(-160, -160, 25)

Definition at line 84 of file NumubarCCIncCuts.cxx.

Referenced by vertex_fiducial_optimization_plots().

const std::vector< double > ana::xsec::numubarcc::mukeedges
Initial value:
{-10.0,
0.5, 0.6, 0.7, 0.8, 0.9,
1.0, 1.1, 1.2, 1.3, 1.4,
1.5, 1.6, 1.7, 1.8, 1.9,
2.0, 2.1, 2.2, 2.3, 2.4,
2.5, 120}

Definition at line 16 of file NumubarCCIncBins.cxx.

const Binning ana::xsec::numubarcc::muon_ke_bins = Binning::Custom(mukeedges)

Definition at line 22 of file NumubarCCIncBins.cxx.

const double ana::xsec::numubarcc::MUON_MASS = TDatabasePDG::Instance()->GetParticle(13)->Mass()

Definition at line 19 of file NumubarCCIncVars.cxx.

const std::string ana::xsec::numubarcc::nominal_dataset = nominal_prod5_1

Definition at line 14 of file NumubarCCIncSystDefs.cxx.

Referenced by fiducial_opt(), muonid_opt(), and muonid_rhc_dist().

const std::string ana::xsec::numubarcc::nominal_dataset_parent = nominal_prod5_1

Definition at line 15 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::nominal_decaf = nominal_prod5_1_decaf

Definition at line 16 of file NumubarCCIncSystDefs.cxx.

Referenced by energy_resolution().

const std::string ana::xsec::numubarcc::nominal_prod5 = "nd_rhc_prod5_caf_minus_muonid_training_minus_fakedata"

Dataset definitions.

Definition at line 6 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::nominal_prod5_1 = "prod_caf_R20-11-25-prod5.1reco.a_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_v1"

Definition at line 10 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::nominal_prod5_1_decaf = "prod_sumdecaf_development_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_ndphysics_contain_v1"

Definition at line 11 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::nominal_prod5_1_flatcaf = "prod_flatsumdecaf_development_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_ndphysics_contain_v1"

Definition at line 12 of file NumubarCCIncSystDefs.cxx.

Referenced by binningResolution(), and effPurityStudies().

const std::string ana::xsec::numubarcc::nominal_prod5_decaf = "nd_rhc_prod5_decaf_minus_muonid_training_minus_fakedata"

Definition at line 8 of file NumubarCCIncSystDefs.cxx.

const std::string ana::xsec::numubarcc::nominal_prod5_parent = "prod_caf_R19-11-18-prod5reco.d_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_v1"

Definition at line 7 of file NumubarCCIncSystDefs.cxx.

const TVector3 * ana::xsec::numubarcc::numucc_fiducial_max = new TVector3( 140, 140, 1000)

Definition at line 88 of file NumubarCCIncCuts.cxx.

const TVector3 * ana::xsec::numubarcc::numucc_fiducial_min = new TVector3(-130, -130, 100)

Definition at line 87 of file NumubarCCIncCuts.cxx.

const ana::Var ana::xsec::numubarcc::std_wgt = ana::VarFromNuTruthVar(std_wgt_ST, 1)
const ana::NuTruthVar ana::xsec::numubarcc::std_wgt_ST = flux_wgt_ST*xsec_wgt_ST

Definition at line 29 of file NumubarCCIncVars.h.

Referenced by binningResolution(), and effPurityStudies().

std::map< const std::string, const std::string > ana::xsec::numubarcc::syst_datasets
Initial value:
{
{"cherenkov", "prod_caf_R19-11-18-prod5reco.j_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_cherenkov_v1"},
{"lightup", "prod_caf_R19-11-18-prod5reco.j_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_llup_v1"},
{"lightdw", "prod_caf_R19-11-18-prod5reco.j_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_lldown_v1"},
{"calibup", "prod_caf_R19-11-18-prod5reco.l_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_calibup_v1"},
{"calibdw", "prod_caf_R19-11-18-prod5reco.l_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_calibdown_v1"},
{"calibshape", "prod_caf_R19-11-18-prod5reco.m_nd_genie_N1810j0211a_nonswap_rhc_nova_v08_full_calibshape_v1"}
}

Definition at line 30 of file NumubarCCIncSystDefs.cxx.

Referenced by fiducial_opt(), and muonid_opt().

const bool ana::xsec::numubarcc::testo = true

Definition at line 6 of file NumubarCCIncCuts.cxx.

const ana::Var ana::xsec::numubarcc::xsec_wgt = ana::VarFromNuTruthVar(xsec_wgt_ST, 1)

Definition at line 26 of file NumubarCCIncVars.h.

Referenced by fiducial_opt().

const ana::NuTruthVar ana::xsec::numubarcc::xsec_wgt_ST = ana::kXSecCVWgt2020_NT

Definition at line 25 of file NumubarCCIncVars.h.