Added StrawGasStepFilter

Author: bonventre

Filters/fcl/SGSFilter.fcl

@@ -0,0 +1,40 @@
+#include "Offline/fcl/standardServices.fcl"
+process_name : SGSFilter
+
+source : {
+  module_type: RootInput
+}
+services : @local::Services.Core
+
+physics : {
+  filters : {
+    SGSFilter : {
+      StrawGasStepTag: "compressDetStepMCs"
+      MinSteps: 10
+      module_type: "StrawGasStepFilter"
+    }
+  }
+  FilterPath: [
+    SGSFilter
+  ]
+  e1        : [FilterOutput]
+  end_paths : [e1]
+}
+physics.trigger_paths : [ FilterPath ]
+
+outputs: {
+   FilterOutput: {
+      SelectEvents: [
+         "FilterPath"
+      ]
+      module_type: "RootOutput"
+      outputCommands: [
+         "keep *_*_*_*"
+      ]
+   }
+}
+
+services.ProditionsService.trackerStatus.Settings.useDb: true
+services.ProditionsService.trackerStatus.Settings.verbose: 2
+#services.DbService.textFile : ["TrackerStatus.txt" ]
+services.DbService.verbose: 2

Filters/src/SConscript

@@ -32,6 +32,8 @@ helper.make_plugins( [
   'mu2e_GeneralUtilities',
   'mu2e_DbService',
   'mu2e_DbTables',
+  'mu2e_TrackerConditions',
+  'mu2e_ConditionsService',
   babarlibs,
   'art_Framework_Core',
   'art_Framework_Principal',

Filters/src/StrawGasStepFilter_module.cc

@@ -0,0 +1,59 @@
+// requires a minimum number of steps in online straws
+
+#include "art/Framework/Core/EDFilter.h"
+#include "art/Framework/Principal/Event.h"
+#include "art/Framework/Principal/Handle.h"
+#include "fhiclcpp/ParameterSet.h"
+#include "Offline/MCDataProducts/inc/StrawGasStep.hh"
+#include "Offline/ProditionsService/inc/ProditionsHandle.hh"
+#include "Offline/TrackerConditions/inc/TrackerStatus.hh"
+
+#include <string>
+
+using namespace std;
+
+namespace mu2e{
+
+  class StrawGasStepFilter : public art::EDFilter {
+  public:
+    using Name=fhicl::Name;
+    using Comment=fhicl::Comment;
+
+    struct Config {
+      fhicl::Atom<art::InputTag> stepstag { Name("StrawGasStepTag"), Comment("StrawGasStep tag")};
+      fhicl::Atom<unsigned> minsteps { Name("MinSteps"), Comment("Minimum number of SGS")};
+    };
+    using Parameters = art::EDFilter::Table<Config>;
+    explicit StrawGasStepFilter(const Parameters& config);
+    virtual bool filter(art::Event& e);
+
+  private:
+    art::InputTag _stepsTag;
+    size_t        _minSteps;
+
+    ProditionsHandle<TrackerStatus> _trackerStatus_h;
+
+
+  };
+
+  StrawGasStepFilter::StrawGasStepFilter(const Parameters& config) :
+    art::EDFilter{config},
+    _stepsTag(config().stepstag()),
+    _minSteps(config().minsteps())
+    {
+    }
+
+  bool StrawGasStepFilter::filter(art::Event& event)  {
+    auto const& trackerStatus = _trackerStatus_h.getPtr(event.id());
+
+    size_t count = 0;
+    auto steps = event.getValidHandle<StrawGasStepCollection>(_stepsTag);
+    for (size_t i=0;i<steps->size();i++){
+      if (!trackerStatus->noSignal(steps->at(i).strawId()))
+        count += 1;
+    }
+    return count >= _minSteps;
+  }
+}
+
+DEFINE_ART_MODULE(mu2e::StrawGasStepFilter)