GetDetectorType.cxx

Go to the documentation of this file.

#include "larpandora/LArPandoraInterface/Detectors/GetDetectorType.h"
#include "larcore/Geometry/WireReadout.h"
#include "larpandora/LArPandoraInterface/Detectors/DUNEFarDetVDThreeView.h"
#include "larpandora/LArPandoraInterface/Detectors/ICARUS.h"
#include "larpandora/LArPandoraInterface/Detectors/LArPandoraDetectorType.h"
#include "larpandora/LArPandoraInterface/Detectors/ProtoDUNEDualPhase.h"
#include "larpandora/LArPandoraInterface/Detectors/VintageLArTPCThreeView.h"

#include "cetlib_except/exception.h"

#include <limits>
#include <set>

namespace lar_pandora {

  LArPandoraDetectorType* detector_functions::GetDetectorType()
  {
    auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout>()->Get();

    const unsigned int nPlanes(wireReadoutGeom.MaxPlanes());
    std::set<geo::View_t> planeSet;
    for (auto const& plane : wireReadoutGeom.Iterate<geo::PlaneGeo>(geo::TPCID{0, 0}))
      planeSet.insert(plane.View());

    if (nPlanes == 3 && planeSet.count(geo::kU) && planeSet.count(geo::kY) &&
        planeSet.count(geo::kZ)) {
      return new DUNEFarDetVDThreeView; //TODO Address bare pointer
    }
    if (nPlanes == 3 && planeSet.count(geo::kU) && planeSet.count(geo::kV) &&
        planeSet.count(geo::kW)) {
      return new VintageLArTPCThreeView;
    }
    if (nPlanes == 3 && planeSet.count(geo::kU) && planeSet.count(geo::kV) &&
        planeSet.count(geo::kY)) {
      return new ICARUS;
    }
    if (nPlanes == 2 && planeSet.count(geo::kW) && planeSet.count(geo::kY)) {
      return new ProtoDUNEDualPhase;
    }

    throw cet::exception("LArPandora") << "LArPandoraDetectorType::GetDetectorType --- unable to "
                                          "determine the detector type from the geometry GDML";
  }

} // namespace lar_pandora