1 files changed, 106 insertions, 0 deletions
diff --git a/ext/dsent/model/optical/GatedLaserSource.cc b/ext/dsent/model/optical/GatedLaserSource.cc
new file mode 100644
index 000000000..e32474daf
--- /dev/null
+++ b/ext/dsent/model/optical/GatedLaserSource.cc
@@ -0,0 +1,106 @@
+#include "model/optical/GatedLaserSource.h"
+
+#include "model/optical_graph/OpticalWaveguide.h"
+#include "model/optical_graph/OpticalWavelength.h"
+#include "model/optical_graph/OpticalLaser.h"
+#include "model/optical_graph/OpticalGraph.h"
+
+namespace DSENT
+{
+    GatedLaserSource::GatedLaserSource(const String& instance_name_, const TechModel* tech_model_)
+        : OpticalModel(instance_name_, tech_model_)
+    {
+        initParameters();
+        initProperties();
+    }
+
+    GatedLaserSource::~GatedLaserSource()
+    {}
+
+    void GatedLaserSource::initParameters()
+    {
+        addParameterName("OutStart");
+        addParameterName("OutEnd");
+        addParameterName("MaxDetectors");
+        return;
+    }
+
+    void GatedLaserSource::initProperties()
+    {
+		addPropertyName("OptUtil", 1.0);
+        addPropertyName("LaserEventTime");
+        return;
+    }
+
+    void GatedLaserSource::constructModel()
+    {
+        // Create Area result
+        Result* area_result = new AtomicResult("Photonic");
+        addAreaResult(area_result);
+        // Create NDD power result
+        Result* energy_result = new AtomicResult("Laser");
+        addEventResult(energy_result);
+    
+        // Get parameters
+        WavelengthGroup laser_wavelengths = makeWavelengthGroup(getParameter("OutStart"), getParameter("OutEnd"));
+
+        // Create optical ports
+        createOpticalOutputPort(    "Out",      laser_wavelengths);
+        // Create the filter
+        createLaser(                "Laser",    laser_wavelengths);
+        OpticalLaser* laser = getLaser("Laser");
+        // Connect the laser to the output
+        laser->addDownstreamNode(getWaveguide("Out"));
+    }
+    
+    void GatedLaserSource::updateModel()
+    {
+        // Get properties
+        double laser_efficiency = getTechModel()->get("Laser->CW->Efficiency");
+        double laser_area = getTechModel()->get("Laser->CW->Area");
+        double laser_diode_loss = getTechModel()->get("Laser->CW->LaserDiodeLoss");
+
+        // Get parameters
+        WavelengthGroup laser_wavelengths = makeWavelengthGroup(getParameter("OutStart"), getParameter("OutEnd"));
+        unsigned int number_wavelengths = laser_wavelengths.second - laser_wavelengths.first + 1;
+        // Update losses
+        OpticalLaser* laser = getLaser("Laser");
+        laser->setLoss(laser_diode_loss);
+        laser->setEfficiency(laser_efficiency);
+        // Update area
+        getAreaResult("Photonic")->setValue(laser_area * number_wavelengths);
+    }
+    
+    void GatedLaserSource::evaluateModel()
+    {    
+        // Get parameters
+        unsigned int max_detectors = getParameter("MaxDetectors");
+        double laser_event_time = getProperty("LaserEventTime");
+        WavelengthGroup laser_wavelengths = makeWavelengthGroup(getParameter("OutStart"), getParameter("OutEnd"));
+                
+		// Get properties
+		double opt_util = getProperty("OptUtil");
+
+        // Create optical graph object
+        OpticalGraph* optical_graph = new OpticalGraph("LaserTrace", this);        
+        // Ask optical graph object to perform power optimization
+        bool success = optical_graph->performPowerOpt(getLaser("Laser"), laser_wavelengths, max_detectors, opt_util);
+        if (!success)
+        {
+            Log::printLine(std::cerr, "[Warning] " + getInstanceName() +
+                " -> Wavelengths contains data paths with no possible modulator configurations!");
+        }        
+        // Trace the wavelengths the laser is outputting to find the output
+        // power needed by the laser        
+        OpticalWavelength* wavelength = optical_graph->traceWavelength(laser_wavelengths, getLaser("Laser"));                
+        // Calculate the power needed by the wavelength
+        double laser_power = wavelength->getLaserPower(max_detectors);
+        // Calculate NDD power
+        getEventResult("Laser")->setValue(laser_power * laser_event_time);
+        
+        delete wavelength;
+        delete optical_graph;
+    }
+    
+} // namespace DSENT
+