ISTPC.cxx

Go to the documentation of this file.

//Class: ISTPC
//File: ISTPC.h and ISTPC.cxx
//Description: Class including common functions needed for using the
//Hybrid Model for Photon Propagation
//UGR, 2021

#include "larsim/IonizationScintillation/ISTPC.h"
#include "larcore/CoreUtils/ServiceUtil.h"
#include "larcorealg/CoreUtils/counter.h"
#include "larcorealg/Geometry/CryostatGeo.h"
#include "larcorealg/Geometry/GeometryCore.h"

#include "messagefacility/MessageLogger/MessageLogger.h"

#include "range/v3/view/enumerate.hpp"

namespace larg4 {
  //----------------------------------------------------------------------------
  ISTPC::ISTPC(geo::GeometryCore const& geom) : fActiveVolumes{extractActiveLArVolume(geom)}
  {
    mf::LogTrace("IonAndScint") << "IonizationAndScintillation/ISTPC Initialize.\n"
                                << "Initializing the geometry of the detector.";

    {
      auto log = mf::LogTrace("IonAndScint") << "IonAndScint: active volume boundaries from "
                                             << fActiveVolumes.size() << " volumes:";
      for (auto const& [iCryo, box] : fActiveVolumes | ranges::views::enumerate) {
        log << "\n - C:" << iCryo << ": " << box.Min() << " -- " << box.Max() << " cm";
      }
    } // local scope
  }

  //----------------------------------------------------------------------------

  bool ISTPC::isScintInActiveVolume(geo::Point_t const& ScintPoint)
  {
    for (auto const& box : fActiveVolumes) {
      if (box.ContainsPosition(ScintPoint)) { return true; }
    }
    return false;
  }
  //----------------------------------------------------------------------------

  std::vector<geo::BoxBoundedGeo> ISTPC::extractActiveLArVolume(geo::GeometryCore const& geom)
  {
    std::vector<geo::BoxBoundedGeo> activeVolumes;
    activeVolumes.reserve(geom.Ncryostats());

    for (geo::CryostatGeo const& cryo : geom.Iterate<geo::CryostatGeo>()) {

      // can't use it default-constructed since it would always include origin

      geo::BoxBoundedGeo box{cryo.TPC(0).ActiveBoundingBox()};

      for (geo::TPCGeo const& TPC : cryo.IterateTPCs())
        box.ExtendToInclude(TPC.ActiveBoundingBox());

      activeVolumes.push_back(std::move(box));

    } // for cryostats

    return activeVolumes;
  }

}