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 "larcore/CoreUtils/ServiceUtil.h"
#include "larcorealg/CoreUtils/counter.h"
#include "larcorealg/CoreUtils/enumerate.h"
#include "larcorealg/Geometry/CryostatGeo.h"
#include "larcorealg/Geometry/GeometryCore.h"
#include "larsim/IonizationScintillation/ISTPC.h"

#include "messagefacility/MessageLogger/MessageLogger.h"


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] : util::enumerate(fActiveVolumes)) {
        log << "\n - C:" << iCryo << ": " << box.Min() << " -- " << box.Max() << " cm";
      }
    } // local scope
  }


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


 bool
 ISTPC::isScintInActiveVolume(geo::Point_t const& ScintPoint)
 {
   return fActiveVolumes[0].ContainsPosition(ScintPoint);
 }
  //----------------------------------------------------------------------------


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.IterateCryostats()) {

    // 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;
}

}