Add example calibrations factory

Author: nathanwbrei

src/examples/tutorial_with_podio_datamodel/11_calibrations/CMakeLists.txt

@@ -1,7 +1,7 @@
 
 
-#add_jana_plugin(podio_calibrations)
+add_jana_plugin(podio_calibrations)
 
-#target_link_libraries(podio_calibrations PUBLIC podio_datamodel)
+target_link_libraries(podio_calibrations PUBLIC jana2_tutorial_podio_datamodel)
 
 

src/examples/tutorial_with_podio_datamodel/11_calibrations/CalorimeterCluster_factory_filtered.h

@@ -0,0 +1,35 @@
+
+#include <jana2_tutorial_podio_datamodel/CalorimeterClusterCollection.h>
+
+#include <JANA/Components/JOmniFactory.h>
+#include <JANA/Calibrations/JCalibrationManager.h>
+
+
+struct CalorimeterCluster_factory_filtered:
+    public JOmniFactory<CalorimeterCluster_factory_filtered> {
+
+    PodioInput<CalorimeterCluster> m_clusters_in {this};
+    PodioOutput<CalorimeterCluster> m_clusters_out {this};
+
+    Service<JCalibrationManager> m_calib_manager {this};
+    double threshold = 100.0;
+
+    void Configure() {}
+
+    void ChangeRun(int32_t run_nr) {
+        m_clusters_out.SetCollectionName("filtered_clusters");
+        m_calib_manager->GetJCalibration(run_nr)
+                       ->Get("MyCAL/cluster_threshold", threshold);
+    }
+
+    void Execute(int32_t /*run_nr*/, uint64_t /*evt_nr*/) {
+        m_clusters_out()->setSubsetCollection(true);
+        for (auto cluster : *m_clusters_in()) {
+            if (cluster.getEnergy() >= threshold) {
+                m_clusters_out()->push_back(cluster);
+            }
+        }
+    }
+};
+
+

src/libraries/JANA/Topology/JEventPool.cc

@@ -70,7 +70,7 @@ void JEventPool::Ingest(JEvent* event, size_t location) {
     // This is necessary for interleaved events
     auto incoming_event_level = event->GetLevel();
     if (incoming_event_level != m_level) {
-        LOG << "Pool " << toString(m_level) << " forwarding event " << event->GetEventNumber() << " to parent pool " << toString(event->GetLevel());
+        //LOG << "Pool " << toString(m_level) << " forwarding event " << event->GetEventNumber() << " to parent pool " << toString(event->GetLevel());
         m_parent_pools.at(incoming_event_level)->Ingest(event, location);
         return;
     }
@@ -79,32 +79,32 @@ void JEventPool::Ingest(JEvent* event, size_t location) {
     auto finished_parents = event->ReleaseAllParents();
     // TODO: I'd prefer to not have to do an allocation each time, but this will work for now
     for (auto* parent : finished_parents) {
-        LOG << "JEventPool::Ingest: Found finished parent of level " << toString(parent->GetLevel());
+        //LOG << "JEventPool::Ingest: Found finished parent of level " << toString(parent->GetLevel());
         m_parent_pools.at(parent->GetLevel())->NotifyThatAllChildrenFinished(parent, location);
         // TODO: This is likely the wrong location. Obtain from parent event?
     }
 
     if (event->GetChildCount() == 0) {
         // There's no way for additional children to appear because Ingest takes the "original" parent
         //LOG << "JEventPool::Ingest: event at level " << toString(m_level) << " is PUSHED";
-        LOG << "JEventPool::Ingest: " << toString(m_level) << " event is pushed";
+        //LOG << "JEventPool::Ingest: " << toString(m_level) << " event is pushed";
         Push(event, location);
     }
     else {
         // We've received the original but we can't push it until all children have been pushed to 
         // _their_ pools, in which case we push once we receive the notification
-        LOG << "JEventPool::Ingest: " << toString(m_level) << " event is pending";
+        //LOG << "JEventPool::Ingest: " << toString(m_level) << " event is pending";
         m_pending.insert(event);
     }
 }
 
 
 void JEventPool::NotifyThatAllChildrenFinished(JEvent* event, size_t location) {
-    LOG << "JEventPool::Notify called for level " << toString(m_level);
+    //LOG << "JEventPool::Notify called for level " << toString(m_level);
     size_t was_present = m_pending.erase(event);
     if (was_present == 1) {
         Push(event, location);
-        LOG << "JEventPool at level " << toString(m_level) << " has pushed a parent event";
+        //LOG << "JEventPool at level " << toString(m_level) << " has pushed a parent event";
     }
 }
 

src/libraries/JANA/Topology/JMultilevelSourceArrow.cc

@@ -80,7 +80,7 @@ void JMultilevelSourceArrow::fire(JEvent* input, OutputData& outputs, size_t& ou
                     // Note that this only attaches parents that we already, so if the parents arrive in the wrong order they
                     // will just be missing. If this is expected behavior, you'll need to set your downstream parent inputs to be optional.
                     if (parent_pair.first != nullptr) {
-                        LOG_INFO(get_logger()) << "JMultilevelSourceArrow: Attaching parent: " << parent_pair.first->GetEventStamp() << " to event " << input->GetEventStamp();
+                        LOG_TRACE(get_logger()) << "JMultilevelSourceArrow: Attaching parent: " << parent_pair.first->GetEventStamp() << " to event " << input->GetEventStamp();
                         input->SetParent(parent_pair.first);
                     }
                 }