test: try it

JeffersonLab · Sep 28, 2024 · ae05cdf · ae05cdf
1 parent 6c9c250
commit ae05cdf
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 4 deletions.
diff --git a/src/iguana/algorithms/physics/InclusiveKinematics/Algorithm.cc b/src/iguana/algorithms/physics/InclusiveKinematics/Algorithm.cc
@@ -1,5 +1,4 @@
 #include "Algorithm.h"
-#include "iguana/services/RCDBReader.h"
 
 // ROOT
 #include <Math/Vector4D.h>
@@ -33,6 +32,9 @@ namespace iguana::physics {
     i_qz     = result_schema.getEntryOrder("qz");
     i_qE     = result_schema.getEntryOrder("qE");
 
+    // instantiate RCDB reader
+    m_rcdb = std::make_unique<RCDBReader>("RCDB|" + GetName());
+
     // parse config file
     ParseYAMLConfig();
     o_runnum         = ConcurrentParamFactory::Create<int>();
@@ -177,8 +179,12 @@ namespace iguana::physics {
     m_log->Trace("-> calling Reload({}, {})", runnum, key);
     o_runnum->Save(runnum, key);
 
+    // get the beam energy //////////////////////////////////////////////////////////
+    auto beam_energy = m_rcdb->GetBeamEnergy(runnum);
+    m_log->Error("Beam energy = {}", beam_energy);
+    /////////////////////////////////////////////////////////////////////////////////
+
     // parse config params
-    auto beam_energy     = GetOptionScalar<double>("beam_energy", {"initial_state", GetConfig()->InRange("runs", runnum), "beam_energy"});
     auto beam_direction  = GetOptionVector<double>("beam_direction", {"initial_state", GetConfig()->InRange("runs", runnum), "beam_direction"});
     auto target_particle = GetOptionScalar<std::string>("target_particle", {"initial_state", GetConfig()->InRange("runs", runnum), "target_particle"});
 

diff --git a/src/iguana/algorithms/physics/InclusiveKinematics/Algorithm.h b/src/iguana/algorithms/physics/InclusiveKinematics/Algorithm.h
@@ -3,6 +3,7 @@
 #include "iguana/algorithms/Algorithm.h"
 #include "iguana/algorithms/TypeDefs.h"
 #include "iguana/services/ConcurrentParam.h"
+#include "iguana/services/RCDBReader.h"
 
 namespace iguana::physics {
 
@@ -114,6 +115,8 @@ namespace iguana::physics {
       enum method_lepton_finder { highest_energy_FD_trigger };
       method_reconstruction o_method_reconstruction;
       method_lepton_finder o_method_lepton_finder;
+
+      std::unique_ptr<RCDBReader> m_rcdb;
   };
 
 }
diff --git a/src/iguana/algorithms/physics/InclusiveKinematics/Config.yaml b/src/iguana/algorithms/physics/InclusiveKinematics/Config.yaml
@@ -1,8 +1,7 @@
 physics::InclusiveKinematics:
 
-  initial_state: # FIXME: add run ranges
+  initial_state:
     - default:
-      beam_energy: 10.6
       beam_direction: [ 0.0, 0.0, 1.0 ]
       target_particle: proton