diff --git a/sbncode/CAFMaker/CAFMaker_module.cc b/sbncode/CAFMaker/CAFMaker_module.cc index 9f2c01590..3b0ea775c 100644 --- a/sbncode/CAFMaker/CAFMaker_module.cc +++ b/sbncode/CAFMaker/CAFMaker_module.cc @@ -664,6 +664,11 @@ void CAFMaker::InitVolumes() { double YMax = std::max_element(tpcs.begin(), tpcs.end(), [](auto &lhs, auto &rhs) { return lhs.MaxY() < rhs.MaxY(); })->MaxY(); double ZMax = std::max_element(tpcs.begin(), tpcs.end(), [](auto &lhs, auto &rhs) { return lhs.MaxZ() < rhs.MaxZ(); })->MaxZ(); + std::cout << "Starting new active volume!" << std::endl; + std::cout << "XMin: " << XMin << " XMax: " << XMax << std::endl; + std::cout << "YMin: " << YMin << " YMax: " << YMax << std::endl; + std::cout << "ZMin: " << ZMin << " ZMax: " << ZMax << std::endl; + fActiveVolumes.emplace_back(XMin, XMax, YMin, YMax, ZMin, ZMax); } } diff --git a/sbncode/CAFMaker/FillTrue.cxx b/sbncode/CAFMaker/FillTrue.cxx index ef1b4cc0c..9d44309b8 100644 --- a/sbncode/CAFMaker/FillTrue.cxx +++ b/sbncode/CAFMaker/FillTrue.cxx @@ -464,6 +464,7 @@ namespace caf { const std::vector &active_volumes) { srneutrino.index = i; + srneutrino.cryostat = -1001; for (int c = 0; c < 2; c++) { SRTrueInteractionPlaneInfo init; @@ -599,8 +600,10 @@ namespace caf { } // Set the cryostat of the position - for (int icryo = 0; icryo < 2; icryo++) { + for (int icryo = 0; icryo < (int)active_volumes.size(); icryo++) { + std::cout << "You have entered loop! " << icryo << std::endl; if (active_volumes[icryo].ContainsPosition(nu.Nu().Position().Vect())) { + std::cout << "Boolean passed!" << std::endl; srneutrino.cryostat = icryo; break; }