599 std::map<std::string,const Cut> map_FurtherBenBreakdown;
602 map_FurtherBenBreakdown.emplace(
"over_fid_sig_ATE_rec_presel_minus_fidana",
604 map_FurtherBenBreakdown.emplace(
"over_fid_sig_ATE_rec_presel_minus_cont",
606 map_FurtherBenBreakdown.emplace(
"over_fid_sig_ATE_rec_presel_minus_nhit",
608 map_FurtherBenBreakdown.emplace(
"over_fid_sig_ATE_rec_presel_minus_energy",
610 map_FurtherBenBreakdown.emplace(
"over_fid_sig_ATE_rec_presel_minus_prongl",
613 for(
auto furthercut : map_FurtherBenBreakdown){
614 Cut kFurtherCut = furthercut.second;
615 const Cut kThisCut = kBaseCut && kFurtherCut;
616 map_Cuts.emplace(furthercut.first, kThisCut);
const Cut kNue2020NDNHits
const Cut kNue2020NDContain([](const caf::SRProxy *sr){for(unsigned int ix=0;ix< sr->vtx.elastic.fuzzyk.nshwlid;++ix){const caf::SRVector3DProxy &start=sr->vtx.elastic.fuzzyk.png[ix].shwlid.start;const caf::SRVector3DProxy &stop=sr->vtx.elastic.fuzzyk.png[ix].shwlid.stop;if(std::min(start.X(), stop.X())< -150.0) return false;if(std::max(start.X(), stop.X()) > 170.0) return false;if(std::min(start.Y(), stop.Y())< -170.0) return false;if(std::max(start.Y(), stop.Y()) > 150.0) return false;if(std::min(start.Z(), stop.Z())< 100.0) return false;if(std::max(start.Z(), stop.Z()) > 1225.0) return false;}if(sr->sel.contain.nplanestofront< 10) return false;return true;})
const Cut kNue2020NDProngLength
const Cut kNue2020NDEnergy
const Cut kNue2020RecoQuality([](const caf::SRProxy *sr){if(sr->sel.nuecosrej.hitsperplane >=8) return false;if(!sr->vtx.elastic.IsValid) return false;if(sr->vtx.elastic.fuzzyk.npng==0) return false;return true;})
const Cut kNue2020NDFiducial([](const caf::SRProxy *sr){assert(sr->vtx.elastic.IsValid &&"Must apply DQ cuts");if(sr->vtx.elastic.vtx.X()< -100.0) return false;if(sr->vtx.elastic.vtx.X() > 160.0) return false;if(sr->vtx.elastic.vtx.Y()< -160.0) return false;if(sr->vtx.elastic.vtx.Y() > 100.0) return false;if(sr->vtx.elastic.vtx.Z()< 150.0) return false;if(sr->vtx.elastic.vtx.Z() > 900.0) return false;return true;})