Feature/acastill barycenter fm upgrade

sbndcode/EvtsFiltering/CMakeLists.txt

@@ -0,0 +1,54 @@
+
+cet_build_plugin(FilterEventID art::module 
+        lardataobj::RecoBase
+        ROOT::Core
+        ROOT::Tree
+        art::Framework_Core
+        art::Framework_Principal
+        art::Utilities
+        art_root_io::tfile_support
+        art_root_io::TFileService_service
+        canvas::canvas
+        fhiclcpp::fhiclcpp
+        cetlib::cetlib
+        cetlib_except::cetlib_except
+        messagefacility::MF_MessageLogger
+        lardataobj::RawData
+        larcore::Geometry_Geometry_service
+        fhiclcpp::fhiclcpp
+        cetlib::cetlib
+        CLHEP::Random
+        ROOT::Geom
+        ROOT::XMLIO
+        ROOT::Gdml
+        ROOT::FFTW
+)
+
+cet_build_plugin(NoFilter art::module 
+        lardataobj::RecoBase
+        ROOT::Core
+        ROOT::Tree
+        art::Framework_Core
+        art::Framework_Principal
+        art::Utilities
+        art_root_io::tfile_support
+        art_root_io::TFileService_service
+        canvas::canvas
+        fhiclcpp::fhiclcpp
+        cetlib::cetlib
+        cetlib_except::cetlib_except
+        messagefacility::MF_MessageLogger
+        lardataobj::RawData
+        larcore::Geometry_Geometry_service
+        fhiclcpp::fhiclcpp
+        cetlib::cetlib
+        CLHEP::Random
+        ROOT::Geom
+        ROOT::XMLIO
+        ROOT::Gdml
+        ROOT::FFTW
+)
+
+install_headers()
+install_source()
+install_fhicl()

sbndcode/EvtsFiltering/FilterEventID_module.cc

@@ -0,0 +1,101 @@
+#include <iostream>
+#include <stdlib.h>
+#include <string>
+#include <vector>
+#include <numeric>
+#include <getopt.h>
+#include <chrono>
+#include <float.h>
+
+
+//some ROOT includes
+#include "TInterpreter.h"
+#include "TROOT.h"
+#include "TH1F.h"
+#include "TH2D.h"
+#include "TCanvas.h"
+#include "TTree.h"
+#include "TFile.h"
+#include "TStyle.h"
+#include "TSystem.h"
+#include "TGraph.h"
+#include "TFFTReal.h"
+#include "TTimeStamp.h"
+#include "TLatex.h"
+#include "TLine.h"
+
+
+// art includes
+#include "canvas/Utilities/InputTag.h"
+#include "canvas/Persistency/Common/FindOneP.h"
+#include "canvas/Persistency/Common/FindManyP.h"
+#include "canvas/Persistency/Common/Ptr.h"
+#include "art/Framework/Core/EDFilter.h" 
+#include "art/Framework/Core/ModuleMacros.h"
+#include "art/Framework/Principal/Event.h"
+#include "art/Framework/Principal/Handle.h"
+#include "art/Framework/Principal/Run.h" 
+#include "art/Framework/Principal/SubRun.h" 
+#include "fhiclcpp/ParameterSet.h" 
+#include "messagefacility/MessageLogger/MessageLogger.h"  
+#include "lardataobj/RecoBase/Hit.h"
+#include "lardataobj/RawData/RawDigit.h"
+#include "lardataobj/RawData/RDTimeStamp.h"
+#include "lardataobj/RawData/raw.h"
+#include "larcore/Geometry/Geometry.h"
+#include "canvas/Persistency/Provenance/Timestamp.h"
+
+#include "art_root_io/TFileService.h"
+
+#include "sbndcode/ChannelMaps/TPC/TPCChannelMapService.h"
+
+
+namespace filt{
+
+  class FilterEventID : public art::EDFilter {
+    public:
+
+      explicit FilterEventID(fhicl::ParameterSet const & pset);
+      virtual bool filter(art::Event& evt) override;
+      void reconfigure(fhicl::ParameterSet const& pset);
+      virtual void beginJob() override;
+
+    private:
+
+      int run;
+      int subrun;
+      int event;
+      std::vector<int> filter_evts;
+
+  };
+
+
+  FilterEventID::FilterEventID(fhicl::ParameterSet const & pset)
+  : EDFilter(pset)
+  {
+    this->reconfigure(pset);
+    filter_evts = pset.get<std::vector<int> >("filterevts");
+  }
+
+  void FilterEventID::reconfigure(fhicl::ParameterSet const& pset){
+  }
+
+
+  bool FilterEventID::filter(art::Event & evt){
+
+    run = evt.id().run();
+    subrun = evt.id().subRun();
+    event = evt.id().event();
+
+    bool is_in_list = std::find(filter_evts.begin(), filter_evts.end(), event) != filter_evts.end();
+    std::cout << "event: " << event << " is_in_list: " << is_in_list << std::endl;
+    return is_in_list;
+
+  }
+
+  void FilterEventID::beginJob() {
+  }
+
+  DEFINE_ART_MODULE(FilterEventID)
+
+}

sbndcode/EvtsFiltering/NoFilter_module.cc

@@ -0,0 +1,87 @@
+#include <iostream>
+#include <stdlib.h>
+#include <string>
+#include <vector>
+#include <numeric>
+#include <getopt.h>
+#include <chrono>
+#include <float.h>
+
+
+//some ROOT includes
+#include "TInterpreter.h"
+#include "TROOT.h"
+#include "TH1F.h"
+#include "TH2D.h"
+#include "TCanvas.h"
+#include "TTree.h"
+#include "TFile.h"
+#include "TStyle.h"
+#include "TSystem.h"
+#include "TGraph.h"
+#include "TFFTReal.h"
+#include "TTimeStamp.h"
+#include "TLatex.h"
+#include "TLine.h"
+
+
+// art includes
+#include "canvas/Utilities/InputTag.h"
+#include "canvas/Persistency/Common/FindOneP.h"
+#include "canvas/Persistency/Common/FindManyP.h"
+#include "canvas/Persistency/Common/Ptr.h"
+#include "art/Framework/Core/EDFilter.h" 
+#include "art/Framework/Core/ModuleMacros.h"
+#include "art/Framework/Principal/Event.h"
+#include "art/Framework/Principal/Handle.h"
+#include "art/Framework/Principal/Run.h" 
+#include "art/Framework/Principal/SubRun.h" 
+#include "fhiclcpp/ParameterSet.h" 
+#include "messagefacility/MessageLogger/MessageLogger.h"  
+#include "lardataobj/RecoBase/Hit.h"
+#include "lardataobj/RawData/RawDigit.h"
+#include "lardataobj/RawData/RDTimeStamp.h"
+#include "lardataobj/RawData/raw.h"
+#include "larcore/Geometry/Geometry.h"
+#include "canvas/Persistency/Provenance/Timestamp.h"
+
+#include "art_root_io/TFileService.h"
+
+#include "sbndcode/ChannelMaps/TPC/TPCChannelMapService.h"
+
+
+namespace filt{
+
+  class NoFilter : public art::EDFilter {
+    public:
+
+      explicit NoFilter(fhicl::ParameterSet const & pset);
+      virtual bool filter(art::Event& evt) override;
+      void reconfigure(fhicl::ParameterSet const& pset);
+      virtual void beginJob() override;
+
+    private:
+
+  };
+
+
+  NoFilter::NoFilter(fhicl::ParameterSet const & pset)
+  : EDFilter(pset)
+  {
+    this->reconfigure(pset);
+  }
+
+  void NoFilter::reconfigure(fhicl::ParameterSet const& pset){
+  }
+
+
+  bool NoFilter::filter(art::Event & evt){
+    return true;
+  }
+
+  void NoFilter::beginJob() {
+  }
+
+  DEFINE_ART_MODULE(NoFilter)
+
+}

sbndcode/EvtsFiltering/filtereventid.fcl

@@ -0,0 +1,51 @@
+#include "services_sbnd.fcl"
+
+process_name: FilterEventID
+
+services:
+{
+  TFileService: { fileName: "reco_hist.root"}
+  MemoryTracker:     {}
+  TimeTracker:       {}
+  RandomNumberGenerator: {} 
+  message:              @local::sbnd_message_services_prod_debug
+  FileCatalogMetadata:  @local::art_file_catalog_mc
+  GeometryConfigurationWriter:    {}
+  @table::sbnd_basic_services
+  @table::sbnd_random_services
+}
+
+source:
+{
+  module_type: RootInput
+  maxEvents:  -1
+}
+
+outputs:
+{
+ out1:
+ {
+   module_type: RootOutput
+   fileName:    "%ifb_eventidfiltered.root"
+   dataTier:    "filtered"
+   fastCloning: true
+   SelectEvents: [ filter ]
+ }
+}
+
+physics:
+{
+  filters:
+  {
+    filtereventid: {
+      module_type: FilterEventID
+      filterevts: [2223, 2295]
+    }
+  }
+
+  filter: [ filtereventid ]
+  stream1:  [ out1 ]
+  trigger_paths: [ filter ]
+  end_paths:     [ stream1 ]
+}
+

sbndcode/EvtsFiltering/nofilter.fcl

@@ -0,0 +1,50 @@
+#include "services_sbnd.fcl"
+
+process_name: NoFilter
+
+services:
+{
+  TFileService: { fileName: "reco_hist.root"}
+  MemoryTracker:     {}
+  TimeTracker:       {}
+  RandomNumberGenerator: {} 
+  message:              @local::sbnd_message_services_prod_debug
+  FileCatalogMetadata:  @local::art_file_catalog_mc
+  GeometryConfigurationWriter:    {}
+  @table::sbnd_basic_services
+  @table::sbnd_random_services
+}
+
+source:
+{
+  module_type: RootInput
+  maxEvents:  -1
+}
+
+outputs:
+{
+ out1:
+ {
+   module_type: RootOutput
+   fileName:    "merged.root"
+   dataTier:    "filtered"
+   fastCloning: true
+   SelectEvents: [ filter ]
+ }
+}
+
+physics:
+{
+  filters:
+  {
+    nofilter: {
+      module_type: NoFilter
+    }
+  }
+
+  filter: [ nofilter ]
+  stream1:  [ out1 ]
+  trigger_paths: [ filter ]
+  end_paths:     [ stream1 ]
+}
+

sbndcode/OpDetReco/OpFlash/FlashTools/DriftEstimatorPMTRatio_tool.cc

@@ -218,7 +218,7 @@ namespace lightana
         drift_distance=Interpolate(pmtratio);
       return drift_distance;
     }
-    else return fDriftCal[fNCalBins-1]; 
+    else return -999999.; 
   }
 
   double DriftEstimatorPMTRatio::GetPropagationTime(double drift){