Add helper function to make trivial progagator

sfegan · Oct 10, 2023 · 982f29b · 982f29b
1 parent 69da781
commit 982f29b
Showing 1 changed file with 27 additions and 0 deletions.
diff --git a/include/simulation/vcl_iact_array.hpp b/include/simulation/vcl_iact_array.hpp
@@ -250,6 +250,7 @@ template<typename VCLArchitecture> class alignas(VCLArchitecture::vec_bytes) VCL
 
   CALIN_TYPEALIAS(FocalPlaneRayPropagator, calin::simulation::vcl_ray_propagator::VCLFocalPlaneRayPropagator<VCLArchitecture>);
   CALIN_TYPEALIAS(DaviesCottonVCLFocalPlaneRayPropagator, calin::simulation::vcl_ray_propagator::DaviesCottonVCLFocalPlaneRayPropagator<VCLArchitecture>);
+  CALIN_TYPEALIAS(TrivialVCLFocalPlaneRayPropagator, calin::simulation::vcl_ray_propagator::TrivialVCLFocalPlaneRayPropagator<VCLArchitecture>);
 
   VCLIACTArray(calin::simulation::atmosphere::LayeredRefractiveAtmosphere* atm,
     const calin::simulation::detector_efficiency::AtmosphericAbsorption& atm_abs,
@@ -288,6 +289,11 @@ template<typename VCLArchitecture> class alignas(VCLArchitecture::vec_bytes) VCL
     const std::string& propagator_name, SplinePEAmplitudeGenerator* pe_generator = nullptr,
     bool adopt_pe_processor = false, bool adopt_pe_generator = false);
 
+  TrivialVCLFocalPlaneRayPropagator* add_trivial_propagator(PEProcessor* pe_processor,
+    const DetectionEfficiency& detector_efficiency, const AngularEfficiency& fp_angular_efficiency,
+    const std::string& propagator_name, SplinePEAmplitudeGenerator* pe_generator = nullptr,
+    bool adopt_pe_processor = false, bool adopt_pe_generator = false);
+
   void point_telescope_az_el_phi_deg(unsigned iscope, double az_deg, double el_deg, double phi_deg);
   void point_telescope_az_el_deg(unsigned iscope, double az_deg, double el_deg);
 
@@ -565,6 +571,27 @@ VCLIACTArray<VCLArchitecture>::add_davies_cotton_propagator(
     adopt_pe_generator);
 }
 
+template<typename VCLArchitecture>
+calin::simulation::vcl_ray_propagator::TrivialVCLFocalPlaneRayPropagator<VCLArchitecture>*
+VCLIACTArray<VCLArchitecture>::add_trivial_propagator(PEProcessor* pe_processor,
+  const DetectionEfficiency& detector_efficiency, const AngularEfficiency& fp_angular_efficiency,
+  const std::string& propagator_name, SplinePEAmplitudeGenerator* pe_generator,
+  bool adopt_pe_processor, bool adopt_pe_generator)
+{
+  auto* propagator = new calin::simulation::vcl_ray_propagator::TrivialVCLFocalPlaneRayPropagator<VCLArchitecture>(
+    this->rng_, ref_index_, /* adopt_rng= */ false);
+
+  auto* bandwidth_manager = new VCLDCBandwidthManager<VCLArchitecture>(
+    &atm_abs_, detector_efficiency, fp_angular_efficiency, zobs_,
+    config_.detector_energy_lo(), config_.detector_energy_hi(),
+    config_.detector_energy_bin_width(), propagator_name);
+
+  add_propagator(propagator, pe_processor, bandwidth_manager, pe_generator, propagator_name,
+    /* adopt_propagator = */ true, adopt_pe_processor, adopt_pe_generator);
+
+  return propagator;
+}
+
 template<typename VCLArchitecture> void VCLIACTArray<VCLArchitecture>::
 point_telescope_az_el_phi_deg(unsigned iscope,
   double az_deg, double el_deg, double phi_deg)