LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
Public Types | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
noisefilteralg::WireCellNoiseFilter Class Reference
Inheritance diagram for noisefilteralg::WireCellNoiseFilter:
art::EDProducer art::detail::Producer art::detail::LegacyModule art::Modifier art::ModuleBase art::ProductRegistryHelper

Public Types

using ModuleType = EDProducer
 
template<typename UserConfig , typename KeysToIgnore = void>
using Table = Modifier::Table< UserConfig, KeysToIgnore >
 

Public Member Functions

 WireCellNoiseFilter (fhicl::ParameterSet const &pset)
 
virtual ~WireCellNoiseFilter ()
 
void produce (art::Event &evt)
 
void reconfigure (fhicl::ParameterSet const &pset)
 
void doBeginJob (SharedResources const &resources)
 
void doEndJob ()
 
void doRespondToOpenInputFile (FileBlock const &fb)
 
void doRespondToCloseInputFile (FileBlock const &fb)
 
void doRespondToOpenOutputFiles (FileBlock const &fb)
 
void doRespondToCloseOutputFiles (FileBlock const &fb)
 
bool doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doEndRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
 
void fillProductDescriptions ()
 
void registerProducts (ProductDescriptions &productsToRegister)
 
ModuleDescription const & moduleDescription () const
 
void setModuleDescription (ModuleDescription const &)
 
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables () const
 
void sortConsumables (std::string const &current_process_name)
 
std::unique_ptr< Worker > makeWorker (WorkerParams const &wp)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Protected Member Functions

ConsumesCollector & consumesCollector ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 

Private Member Functions

void DoNoiseFilter (art::Event const &evt, const std::vector< raw::RawDigit > &, std::vector< raw::RawDigit > &) const
 

Private Attributes

IDFT::pointer m_dft
 
std::string fDigitModuleLabel
 
bool fDoNoiseFiltering
 
size_t fNumTicksToDropFront
 
size_t fWindowSize
 

Detailed Description

Definition at line 41 of file WireCellNoiseFilter_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer
inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>
using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>
inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

noisefilteralg::WireCellNoiseFilter::WireCellNoiseFilter ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 67 of file WireCellNoiseFilter_module.cc.

References reconfigure().

67  : EDProducer(pset)
68  {
69  this->reconfigure(pset);
70  produces<std::vector<raw::RawDigit>>();
71  }
art::EDProducer::EDProducer
EDProducer(fhicl::ParameterSet const &pset)
Definition: EDProducer.cc:6
noisefilteralg::WireCellNoiseFilter::reconfigure
void reconfigure(fhicl::ParameterSet const &pset)
Definition: WireCellNoiseFilter_module.cc:77
noisefilteralg::WireCellNoiseFilter::~WireCellNoiseFilter ( )
virtual

Definition at line 74 of file WireCellNoiseFilter_module.cc.

74 {}

Member Function Documentation

template<typename T , BranchType BT>
ProductToken< T > art::ModuleBase::consumes ( InputTag const &  tag)
protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

62  {
63  return collector_.consumes<T, BT>(tag);
64  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::consumes
ProductToken< T > consumes(InputTag const &)
Definition: ConsumesCollector.h:50
ConsumesCollector & art::ModuleBase::consumesCollector ( )
protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

58  {
59  return collector_;
60  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename T , BranchType BT>
void art::ModuleBase::consumesMany ( )
protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

76  {
77  collector_.consumesMany<T, BT>();
78  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::ModuleBase::consumesView ( InputTag const &  )
protectedinherited
template<typename T , BranchType BT>
ViewToken<T> art::ModuleBase::consumesView ( InputTag const &  tag)
inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

69  {
70  return collector_.consumesView<T, BT>(tag);
71  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::consumesView
ViewToken< Element > consumesView(InputTag const &)
void art::detail::Producer::doBeginJob ( SharedResources const &  resources)
inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

23  {
24  setupQueues(resources);
25  ProcessingFrame const frame{ScheduleID{}};
26  beginJobWithFrame(frame);
27  }
art::detail::Producer::setupQueues
virtual void setupQueues(SharedResources const &)=0
art::detail::Producer::beginJobWithFrame
virtual void beginJobWithFrame(ProcessingFrame const &)=0
bool art::detail::Producer::doBeginRun ( RunPrincipal &  rp,
ModuleContext const &  mc 
)
inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

66  {
67  auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
68  ProcessingFrame const frame{mc.scheduleID()};
69  beginRunWithFrame(r, frame);
70  r.commitProducts();
71  return true;
72  }
r
TRandom r
Definition: spectrum.C:23
art::detail::Producer::beginRunWithFrame
virtual void beginRunWithFrame(Run &, ProcessingFrame const &)=0
art::RangeSet::forRun
static RangeSet forRun(RunID)
Definition: RangeSet.cc:51
bool art::detail::Producer::doBeginSubRun ( SubRunPrincipal &  srp,
ModuleContext const &  mc 
)
inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

86  {
87  auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
88  ProcessingFrame const frame{mc.scheduleID()};
89  beginSubRunWithFrame(sr, frame);
90  sr.commitProducts();
91  return true;
92  }
art::detail::Producer::beginSubRunWithFrame
virtual void beginSubRunWithFrame(SubRun &, ProcessingFrame const &)=0
art::RangeSet::forSubRun
static RangeSet forSubRun(SubRunID)
Definition: RangeSet.cc:57
void art::detail::Producer::doEndJob ( )
inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

31  {
32  ProcessingFrame const frame{ScheduleID{}};
33  endJobWithFrame(frame);
34  }
art::detail::Producer::endJobWithFrame
virtual void endJobWithFrame(ProcessingFrame const &)=0
bool art::detail::Producer::doEndRun ( RunPrincipal &  rp,
ModuleContext const &  mc 
)
inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

76  {
77  auto r = rp.makeRun(mc, rp.seenRanges());
78  ProcessingFrame const frame{mc.scheduleID()};
79  endRunWithFrame(r, frame);
80  r.commitProducts();
81  return true;
82  }
r
TRandom r
Definition: spectrum.C:23
art::detail::Producer::endRunWithFrame
virtual void endRunWithFrame(Run &, ProcessingFrame const &)=0
bool art::detail::Producer::doEndSubRun ( SubRunPrincipal &  srp,
ModuleContext const &  mc 
)
inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

96  {
97  auto sr = srp.makeSubRun(mc, srp.seenRanges());
98  ProcessingFrame const frame{mc.scheduleID()};
99  endSubRunWithFrame(sr, frame);
100  sr.commitProducts();
101  return true;
102  }
art::detail::Producer::endSubRunWithFrame
virtual void endSubRunWithFrame(SubRun &, ProcessingFrame const &)=0
bool art::detail::Producer::doEvent ( EventPrincipal &  ep,
ModuleContext const &  mc,
std::atomic< std::size_t > &  counts_run,
std::atomic< std::size_t > &  counts_passed,
std::atomic< std::size_t > &  counts_failed 
)
inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

110  {
111  auto e = ep.makeEvent(mc);
112  ++counts_run;
113  ProcessingFrame const frame{mc.scheduleID()};
114  produceWithFrame(e, frame);
115  e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
116  ++counts_passed;
117  return true;
118  }
art::detail::Producer::checkPutProducts_
bool const checkPutProducts_
Definition: Producer.h:70
e
Float_t e
Definition: plot.C:35
art::detail::Producer::produceWithFrame
virtual void produceWithFrame(Event &, ProcessingFrame const &)=0
void noisefilteralg::WireCellNoiseFilter::DoNoiseFilter ( art::Event const &  evt,
const std::vector< raw::RawDigit > &  inputWaveforms,
std::vector< raw::RawDigit > &  outputWaveforms 
) const
private

Definition at line 141 of file WireCellNoiseFilter_module.cc.

References DEFINE_ART_MODULE, fNumTicksToDropFront, fWindowSize, raw::kNone, m_dft, noise(), geo::GeometryCore::Nwires(), art::Event::run(), detinfo::sampling_rate(), and lar::dump::vector().

Referenced by produce().

144  {
145  auto const runNum = e.run();
146 
147  // Recover services we will need
148  const lariov::ChannelStatusProvider& channelStatus =
149  art::ServiceHandle<lariov::ChannelStatusService const>()->GetProvider();
150  const lariov::DetPedestalProvider& pedestalValues =
151  art::ServiceHandle<lariov::DetPedestalService const>()->GetPedestalProvider();
152  const lariov::ElectronicsCalibProvider& elec_provider =
153  art::ServiceHandle<lariov::ElectronicsCalibService const>()->GetProvider();
154  const geo::GeometryCore& geometry = *lar::providerFrom<geo::Geometry>();
155  auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
156 
157  const unsigned int n_channels = inputWaveforms.size();
158 
159  // S&C microboone sampling parameter database
160  const double tick = sampling_rate(clock_data); // 0.5 * units::microsecond;
161  const size_t nsamples = inputWaveforms.at(0).NADC();
162  const size_t windowSize = std::min(fWindowSize, nsamples);
163 
164  // Q&D microboone channel map
165  geo::TPCID const tpcid{0, 0};
166  std::vector<int> uchans(geometry.Nwires({tpcid, 0})), vchans(geometry.Nwires({tpcid, 1})),
167  wchans(geometry.Nwires({tpcid, 2}));
168  const int nchans = uchans.size() + vchans.size() + wchans.size();
169  std::iota(uchans.begin(), uchans.end(), 0);
170  std::iota(vchans.begin(), vchans.end(), vchans.size());
171  std::iota(wchans.begin(), wchans.end(), vchans.size() + uchans.size());
172 
173  // Q&D nominal baseline
174  const double unombl = 2048.0, vnombl = 2048.0, wnombl = 400.0;
175 
176  // Q&D miss-configured channel database
177  std::vector<int> miscfgchan;
178  for (int channelIdx = 0; channelIdx < nchans; channelIdx++) {
179  if (elec_provider.ExtraInfo(channelIdx).GetBoolData("is_misconfigured")) {
180  miscfgchan.push_back(channelIdx);
181  }
182  }
183 
184  const double from_gain_mVfC = 4.7, to_gain_mVfC = 14.0, from_shaping = 1.0 * units::microsecond,
185  to_shaping = 2.0 * units::microsecond;
186 
187  // Recover bad channels from the database
188  std::vector<int> bad_channels;
189  for (int channelIdx = 0; channelIdx < nchans; channelIdx++)
190  if (channelStatus.IsBad(channelIdx)) bad_channels.push_back(channelIdx);
191 
192  // Q&D RC+RC time constant - all have same.
193  const double rcrc = 1.0 * units::millisecond;
194  std::vector<int> rcrcchans(nchans);
195  std::iota(rcrcchans.begin(), rcrcchans.end(), 0);
196 
197  //harmonic noises
198  std::vector<int> harmonicchans(uchans.size() + vchans.size());
199  std::iota(harmonicchans.begin(), harmonicchans.end(), 0);
200 
201  std::vector<int> special_chans;
202  special_chans.push_back(2240);
203 
204  WireCell::SigProc::SimpleChannelNoiseDB::mask_t h36kHz(0, 169, 173);
205  WireCell::SigProc::SimpleChannelNoiseDB::mask_t h108kHz(0, 513, 516);
206  WireCell::SigProc::SimpleChannelNoiseDB::mask_t hspkHz(0, 17, 19);
207  WireCell::SigProc::SimpleChannelNoiseDB::multimask_t hharmonic;
208  hharmonic.push_back(h36kHz);
209  hharmonic.push_back(h108kHz);
210  WireCell::SigProc::SimpleChannelNoiseDB::multimask_t hspecial;
211  hspecial.push_back(h36kHz);
212  hspecial.push_back(h108kHz);
213  hspecial.push_back(hspkHz);
214 
215  float u_resp_array[120] = {
216  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
217  0.0, 0.0, 0.0, 0.0, 0.364382, 0.387949,
218  0.411053, 0.433979, 0.456863, 0.479746, 0.502641, 0.52554,
219  0.548441, 0.57134, 0.591765, 0.609448, 0.626848, 0.644094,
220  0.661364, 0.678859, 0.695231, 0.710462, 0.726147, 0.742373,
221  0.761332, 0.783313, 0.806325, 0.830412, 0.857676, 0.888412,
222  0.920705, 0.954624, 0.990242, 1.02766, 1.06121, 1.09027,
223  1.12037, 1.15157, 1.18392, 1.21748, 1.25229, 1.28824,
224  1.32509, 1.36256, 1.40051, 1.43907, 1.47857, 1.51933,
225  1.56134, 1.60404, 1.72665, 1.94005, 2.16994, 2.42041,
226  2.69475, 3.07222, 3.67375, 4.60766, 5.91864, 7.30178,
227  8.3715, 8.94736, 8.93705, 8.40339, 7.2212, 5.76382,
228  3.8931, 1.07893, -3.52481, -11.4593, -20.4011, -29.1259,
229  -34.9544, -36.9358, -35.3303, -31.2068, -25.8614, -20.3613,
230  -15.3794, -11.2266, -7.96091, -5.50138, -3.71143, -2.44637,
231  -1.57662, -0.99733, -0.62554, -0.393562, -0.249715, -0.15914,
232  -0.100771, -0.062443, -0.037283, -0.0211508, -0.0112448, -0.00552085,
233  -0.00245133, -0.000957821, -0.000316912, -8.51679e-05, -2.21299e-05, -1.37496e-05,
234  -1.49806e-05, -1.36935e-05, -9.66758e-06, -5.20773e-06, -7.4787e-07, 3.71199e-06,
235  8.17184e-06, 1.26317e-05, 1.70916e-05, 2.15514e-05, 2.60113e-05, 3.04711e-05};
236  float v_resp_array[120] = {
237  0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
238  0.0, 0.0, 0.0, 0.0, 0.0865303, 0.0925559,
239  0.0983619, 0.104068, 0.109739, 0.115403, 0.121068, 0.126735,
240  0.132403, 0.138072, 0.143739, 0.149408, 0.155085, 0.160791,
241  0.166565, 0.172454, 0.178514, 0.184795, 0.191341, 0.198192,
242  0.205382, 0.212944, 0.220905, 0.229292, 0.238129, 0.247441,
243  0.257256, 0.267601, 0.278502, 0.28999, 0.298745, 0.304378,
244  0.310105, 0.315921, 0.321818, 0.327796, 0.333852, 0.339967,
245  0.346098, 0.352169, 0.358103, 0.363859, 0.36945, 0.374915,
246  0.380261, 0.385401, 0.39016, 0.394378, 0.39804, 0.401394,
247  0.405145, 0.410714, 0.4205, 0.437951, 0.467841, 0.516042,
248  0.587738, 0.694157, 0.840763, 1.01966, 1.22894, 1.5612,
249  2.12348, 3.31455, 5.59355, 9.10709, 14.1756, 18.4603,
250  19.9517, 17.4166, 10.6683, 1.40656, -10.0638, -19.034,
251  -23.654, -24.0558, -21.4418, -17.3229, -12.9485, -9.08912,
252  -6.05941, -3.86946, -2.38669, -1.43678, -0.853335, -0.503951,
253  -0.296551, -0.173029, -0.0990099, -0.0547172, -0.0287882, -0.0142758,
254  -0.00661815, -0.00284757, -0.00115702, -0.000456456, -0.000183439, -8.04214e-05,
255  -4.20533e-05, -2.62903e-05, -1.64098e-05, -6.68039e-06, 3.04903e-06, 1.27784e-05,
256  2.25079e-05, 3.22373e-05, 4.19667e-05, 5.16961e-05, 6.14255e-05, 7.11549e-05};
257  WireCell::Waveform::realseq_t u_resp(nsamples);
258  WireCell::Waveform::realseq_t v_resp(nsamples);
259  for (int i = 0; i != 120; i++) {
260  u_resp.at(i) = u_resp_array[i];
261  v_resp.at(i) = v_resp_array[i];
262  }
263  // WireCell::Waveform::compseq_t u_resp_freq = WireCell::Waveform::dft(u_resp);
264  // WireCell::Waveform::compseq_t v_resp_freq = WireCell::Waveform::dft(v_resp);
265  WireCell::Waveform::compseq_t u_resp_freq = fwd_r2c(m_dft, u_resp);
266  WireCell::Waveform::compseq_t v_resp_freq = fwd_r2c(m_dft, v_resp);
267 
268  int uplane_time_shift = 79;
269  int vplane_time_shift = 82;
270 
271  // do the coherent subtraction
272  std::vector<std::vector<int>> channel_groups;
273  for (unsigned int i = 0; i != 172; i++) {
274  //for (int i=150;i!=151;i++){
275  std::vector<int> channel_group;
276  for (int j = 0; j != 48; j++) {
277  channel_group.push_back(i * 48 + j);
278  }
279  channel_groups.push_back(channel_group);
280  }
281 
282  auto noise = new WireCell::SigProc::SimpleChannelNoiseDB;
283  // initialize
284  noise->set_sampling(tick, nsamples);
285  // set nominal baseline
286  noise->set_nominal_baseline(uchans, unombl);
287  noise->set_nominal_baseline(vchans, vnombl);
288  noise->set_nominal_baseline(wchans, wnombl);
289 
290  noise->set_response(uchans, u_resp_freq);
291  noise->set_response(vchans, v_resp_freq);
292 
293  noise->set_response_offset(uchans, uplane_time_shift);
294  noise->set_response_offset(vchans, vplane_time_shift);
295 
296  noise->set_pad_window_front(uchans, 20);
297  noise->set_pad_window_back(uchans, 10);
298  noise->set_pad_window_front(vchans, 10);
299  noise->set_pad_window_back(vchans, 10);
300  noise->set_pad_window_front(wchans, 10);
301  noise->set_pad_window_back(wchans, 10);
302 
303  // set misconfigured channels
304  noise->set_gains_shapings(miscfgchan, from_gain_mVfC, to_gain_mVfC, from_shaping, to_shaping);
305  // do the RCRC
306  noise->set_rcrc_constant(rcrcchans, rcrc);
307  // set initial bad channels
308  noise->set_bad_channels(bad_channels);
309  // set the harmonic filter
310  noise->set_filter(harmonicchans, hharmonic);
311  noise->set_filter(special_chans, hspecial);
312  noise->set_channel_groups(channel_groups);
313 
314  for (unsigned int i = 0; i != uchans.size(); i++) {
315  if (uchans.at(i) < 100) {
316  noise->set_min_rms_cut_one(uchans.at(i), 1);
317  noise->set_max_rms_cut_one(uchans.at(i), 5);
318  }
319  else if (uchans.at(i) < 2000) {
320  noise->set_min_rms_cut_one(uchans.at(i), 1.9);
321  noise->set_max_rms_cut_one(uchans.at(i), 11);
322  }
323  else {
324  noise->set_min_rms_cut_one(uchans.at(i), 0.9);
325  noise->set_max_rms_cut_one(uchans.at(i), 5);
326  }
327  }
328  for (unsigned int i = 0; i != vchans.size(); i++) {
329  if (vchans.at(i) < 290 + (int)uchans.size()) {
330  noise->set_min_rms_cut_one(vchans.at(i), 1);
331  noise->set_max_rms_cut_one(vchans.at(i), 5);
332  }
333  else if (vchans.at(i) < 2200 + (int)uchans.size()) {
334  noise->set_min_rms_cut_one(vchans.at(i), 1.9);
335  noise->set_max_rms_cut_one(vchans.at(i), 11);
336  }
337  else {
338  noise->set_min_rms_cut_one(vchans.at(i), 1);
339  noise->set_max_rms_cut_one(vchans.at(i), 5);
340  }
341  }
342 
343  // these are the one after the Hardware Fix
344  if (runNum > 8000) {
345  for (unsigned int i = 0; i != uchans.size(); i++) {
346  if (uchans.at(i) < 600) {
347  noise->set_min_rms_cut_one(uchans.at(i), 1 + (1.7 - 1) / 600. * i);
348  noise->set_max_rms_cut_one(uchans.at(i), 5);
349  }
350  else if (uchans.at(i) < 1800) {
351  noise->set_min_rms_cut_one(uchans.at(i), 1.7);
352  noise->set_max_rms_cut_one(uchans.at(i), 11);
353  }
354  else {
355  noise->set_min_rms_cut_one(uchans.at(i), 1 + (1.7 - 1) / 600. * (2399 - i));
356  noise->set_max_rms_cut_one(uchans.at(i), 5);
357  }
358  }
359  for (unsigned int i = 0; i != vchans.size(); i++) {
360  if (vchans.at(i) < 600 + (int)uchans.size()) {
361  noise->set_min_rms_cut_one(vchans.at(i), 0.8 + (1.7 - 0.8) / 600. * i);
362  noise->set_max_rms_cut_one(vchans.at(i), 5);
363  }
364  else if (vchans.at(i) < 1800 + (int)uchans.size()) {
365  noise->set_min_rms_cut_one(vchans.at(i), 1.7);
366  noise->set_max_rms_cut_one(vchans.at(i), 11);
367  }
368  else {
369  noise->set_min_rms_cut_one(vchans.at(i), 0.8 + (1.7 - 0.8) / 600. * (2399 - i));
370  noise->set_max_rms_cut_one(vchans.at(i), 5);
371  }
372  }
373  }
374 
375  noise->set_min_rms_cut(wchans, 1.3);
376  noise->set_max_rms_cut(wchans, 8.0);
377 
378  //Define database object
379  std::shared_ptr<WireCell::IChannelNoiseDatabase> noise_sp(noise);
380 
381  auto one = new WireCell::SigProc::Microboone::OneChannelNoise;
382  one->set_channel_noisedb(noise_sp);
383  std::shared_ptr<WireCell::IChannelFilter> one_sp(one);
384  auto many = new WireCell::SigProc::Microboone::CoherentNoiseSub;
385  many->set_channel_noisedb(noise_sp);
386  std::shared_ptr<WireCell::IChannelFilter> many_sp(many);
387 
388  //define noisefilter object
389  WireCell::SigProc::OmnibusNoiseFilter bus;
390  bus.set_channel_filters({one_sp});
391  bus.set_grouped_filters({many_sp});
392  bus.set_channel_noisedb(noise_sp);
393 
394  // Enable truncation
395  size_t startBin(fNumTicksToDropFront);
396  size_t stopBin(startBin + windowSize);
397 
398  //load waveforms into traces
399  WireCell::ITrace::vector traces;
400  for (unsigned int ich = 0; ich < n_channels; ich++) {
401  if (inputWaveforms.at(ich).NADC() < windowSize) continue;
402 
403  const raw::RawDigit::ADCvector_t& rawAdcVec = inputWaveforms.at(ich).ADCs();
404 
405  WireCell::ITrace::ChargeSequence charges;
406 
407  charges.resize(nsamples);
408 
409  std::transform(rawAdcVec.begin(), rawAdcVec.end(), charges.begin(), [](auto& adcVal) {
410  return float(adcVal);
411  });
412 
413  unsigned int chan = inputWaveforms.at(ich).Channel();
414  SimpleTrace* st = new SimpleTrace(chan, 0.0, charges);
415  traces.push_back(WireCell::ITrace::pointer(st));
416  }
417 
418  //Load traces into frame
419  SimpleFrame* sf = new SimpleFrame(0, 0, traces);
420  WireCell::IFrame::pointer frame = WireCell::IFrame::pointer(sf);
421  WireCell::IFrame::pointer quiet = NULL;
422 
423  //Do filtering
424  bus(frame, quiet);
425 
426  //std::cout << "HERE" << std::endl;
427  //return;
428  if (quiet == NULL) return;
429 
430  //Output results
431  std::vector<short> waveform(windowSize);
432 
433  auto quiet_traces = quiet->traces();
434  for (auto quiet_trace : *quiet_traces.get()) {
435  //int tbin = quiet_trace->tbin();
436  unsigned int channel = quiet_trace->channel();
437 
438  auto& quiet_charges = quiet_trace->charge();
439 
440  // Recover the database version of the pedestal, we'll offset the waveforms so it matches
441  float pedestal = pedestalValues.PedMean(channel);
442 
443  std::transform(quiet_charges.begin() + startBin,
444  quiet_charges.begin() + stopBin,
445  waveform.begin(),
446  [pedestal](auto charge) { return std::round(charge + pedestal); });
447 
448  outputWaveforms.emplace_back(raw::RawDigit(channel, waveform.size(), waveform, raw::kNone));
449  outputWaveforms.back().SetPedestal(pedestal, 1.75);
450  }
451 
452  return;
453  }
raw::RawDigit
Collection of charge vs time digitized from a single readout channel.
Definition: RawDigit.h:68
util::quantities::microsecond
microsecond_as<> microsecond
Type of time stored in microseconds, in double precision.
Definition: spacetime.h:116
raw::RawDigit::ADCvector_t
std::vector< short > ADCvector_t
Type representing a (compressed) vector of ADC counts.
Definition: RawDigit.h:72
raw::kNone
no compression
Definition: RawTypes.h:9
util::quantities::millisecond
millisecond_as<> millisecond
Type of time stored in milliseconds, in double precision.
Definition: spacetime.h:99
noise
unsigned int noise()
Definition: chem4.cc:261
lar::dump::vector
auto vector(Vector const &v)
Returns a manipulator which will print the specified array.
Definition: DumpUtils.h:289
util::quantities::tick
tick_as<> tick
Tick number, represented by std::ptrdiff_t.
Definition: electronics.h:73
geo::TPCID
The data type to uniquely identify a TPC.
Definition: geo_types.h:381
geo::GeometryCore
Description of geometry of one entire detector.
Definition: GeometryCore.h:119
geo::GeometryCore::Nwires
unsigned int Nwires(PlaneID const &planeid) const
Returns the total number of wires in the specified plane.
Definition: GeometryCore.h:1151
e
Float_t e
Definition: plot.C:35
detinfo::sampling_rate
double sampling_rate(DetectorClocksData const &data)
Returns the period of the TPC readout electronics clock.
Definition: DetectorClocksData.h:575
void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const &  fb)
inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

45  {
46  ProcessingFrame const frame{ScheduleID{}};
47  respondToCloseInputFileWithFrame(fb, frame);
48  }
art::detail::Producer::respondToCloseInputFileWithFrame
virtual void respondToCloseInputFileWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const &  fb)
inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

59  {
60  ProcessingFrame const frame{ScheduleID{}};
61  respondToCloseOutputFilesWithFrame(fb, frame);
62  }
art::detail::Producer::respondToCloseOutputFilesWithFrame
virtual void respondToCloseOutputFilesWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const &  fb)
inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

38  {
39  ProcessingFrame const frame{ScheduleID{}};
40  respondToOpenInputFileWithFrame(fb, frame);
41  }
art::detail::Producer::respondToOpenInputFileWithFrame
virtual void respondToOpenInputFileWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const &  fb)
inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

52  {
53  ProcessingFrame const frame{ScheduleID{}};
54  respondToOpenOutputFilesWithFrame(fb, frame);
55  }
art::detail::Producer::respondToOpenOutputFilesWithFrame
virtual void respondToOpenOutputFilesWithFrame(FileBlock const &, ProcessingFrame const &)=0
fb
TFile fb("Li6.root")
void art::Modifier::fillProductDescriptions ( )
inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

11  {
12  ProductRegistryHelper::fillDescriptions(moduleDescription());
13  }
art::ProductRegistryHelper::fillDescriptions
void fillDescriptions(ModuleDescription const &md)
Definition: ProductRegistryHelper.cc:48
art::ModuleBase::moduleDescription
ModuleDescription const & moduleDescription() const
Definition: ModuleBase.cc:13
std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const
inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

44  {
45  return collector_.getConsumables();
46  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::getConsumables
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables() const
Definition: ConsumesCollector.cc:13
std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const &  wp)
inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

38  {
39  return doMakeWorker(wp);
40  }
art::ModuleBase::doMakeWorker
virtual std::unique_ptr< Worker > doMakeWorker(WorkerParams const &wp)=0
template<typename T , BranchType BT>
ProductToken< T > art::ModuleBase::mayConsume ( InputTag const &  tag)
protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

83  {
84  return collector_.mayConsume<T, BT>(tag);
85  }
art::ConsumesCollector::mayConsume
ProductToken< T > mayConsume(InputTag const &)
Definition: ConsumesCollector.h:75
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename T , BranchType BT>
void art::ModuleBase::mayConsumeMany ( )
protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

97  {
98  collector_.mayConsumeMany<T, BT>();
99  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const &  )
protectedinherited
template<typename T , BranchType BT>
ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const &  tag)
inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

90  {
91  return collector_.mayConsumeView<T, BT>(tag);
92  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::mayConsumeView
ViewToken< Element > mayConsumeView(InputTag const &)
ModuleDescription const & art::ModuleBase::moduleDescription ( ) const
inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

14  {
15  if (md_.has_value()) {
16  return *md_;
17  }
18 
19  throw Exception{errors::LogicError,
20  "There was an error while calling moduleDescription().\n"}
21  << "The moduleDescription() base-class member function cannot be called\n"
22  "during module construction. To determine which module is "
23  "responsible\n"
24  "for calling it, find the '<module type>:<module "
25  "label>@Construction'\n"
26  "tag in the message prefix above. Please contact artists@fnal.gov\n"
27  "for guidance.\n";
28  }
art::Exception
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
art::ModuleBase::md_
std::optional< ModuleDescription > md_
Definition: ModuleBase.h:55
void noisefilteralg::WireCellNoiseFilter::produce ( art::Event &  evt)
virtual

Implements art::EDProducer.

Definition at line 88 of file WireCellNoiseFilter_module.cc.

References DoNoiseFilter(), fDigitModuleLabel, fDoNoiseFiltering, fNumTicksToDropFront, fWindowSize, art::ProductRetriever::getByLabel(), art::Handle< T >::isValid(), raw::kNone, and art::Event::put().

89  {
90  // Recover services we will need
91  const lariov::DetPedestalProvider& pedestalValues =
92  art::ServiceHandle<lariov::DetPedestalService const>()->GetPedestalProvider();
93 
94  art::Handle<std::vector<raw::RawDigit>> rawDigitHandle;
95  evt.getByLabel(fDigitModuleLabel, rawDigitHandle);
96 
97  // Define the output vector (in case we don't do anything)
98  std::unique_ptr<std::vector<raw::RawDigit>> filteredRawDigit(new std::vector<raw::RawDigit>);
99 
100  if (rawDigitHandle.isValid() && rawDigitHandle->size() > 0) {
101  const std::vector<raw::RawDigit>& rawDigitVector(*rawDigitHandle);
102 
103  // Make sure we have the correct window size (e.g. window size = 9600 but data is 9595)
104  size_t windowSize(std::min(fWindowSize, rawDigitVector.at(0).NADC()));
105 
106  if (fNumTicksToDropFront + windowSize > rawDigitVector.at(0).NADC())
107  throw cet::exception("WireCellNoiseFilter")
108  << "Ticks to drop + windowsize larger than input buffer\n";
109 
110  if (fDoNoiseFiltering)
111  DoNoiseFilter(evt, rawDigitVector, *filteredRawDigit);
112  else {
113  // Enable truncation
114  size_t startBin(fNumTicksToDropFront);
115  size_t stopBin(startBin + windowSize);
116 
117  raw::RawDigit::ADCvector_t outputVector(windowSize);
118 
119  for (const auto& rawDigit : rawDigitVector) {
120  if (rawDigit.NADC() < windowSize) continue;
121 
122  const raw::RawDigit::ADCvector_t& rawAdcVec = rawDigit.ADCs();
123 
124  unsigned int channel = rawDigit.Channel();
125  float pedestal = pedestalValues.PedMean(channel);
126 
127  std::copy(
128  rawAdcVec.begin() + startBin, rawAdcVec.begin() + stopBin, outputVector.begin());
129 
130  filteredRawDigit->emplace_back(
131  raw::RawDigit(channel, outputVector.size(), outputVector, raw::kNone));
132  filteredRawDigit->back().SetPedestal(pedestal, 2.0);
133  }
134  }
135  }
136 
137  //filtered raw digits
138  evt.put(std::move(filteredRawDigit));
139  }
raw::RawDigit
Collection of charge vs time digitized from a single readout channel.
Definition: RawDigit.h:68
noisefilteralg::WireCellNoiseFilter::DoNoiseFilter
void DoNoiseFilter(art::Event const &evt, const std::vector< raw::RawDigit > &, std::vector< raw::RawDigit > &) const
Definition: WireCellNoiseFilter_module.cc:141
raw::RawDigit::ADCvector_t
std::vector< short > ADCvector_t
Type representing a (compressed) vector of ADC counts.
Definition: RawDigit.h:72
raw::kNone
no compression
Definition: RawTypes.h:9
art::Handle::isValid
bool isValid() const noexcept
Definition: Handle.h:203
art::Event::put
PutHandle< PROD > put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
Definition: Event.h:77
art::ProductRetriever::getByLabel
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Definition: ProductRetriever.h:444
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
void noisefilteralg::WireCellNoiseFilter::reconfigure ( fhicl::ParameterSet const &  pset)

Definition at line 77 of file WireCellNoiseFilter_module.cc.

References fDigitModuleLabel, fDoNoiseFiltering, fNumTicksToDropFront, fWindowSize, fhicl::ParameterSet::get(), and m_dft.

Referenced by WireCellNoiseFilter().

78  {
79  auto dft_tn = pset.get<std::string>("dft", "FftwDFT");
80  m_dft = Factory::find_tn<IDFT>(dft_tn);
81  fDigitModuleLabel = pset.get<std::string>("DigitModuleLabel", "daq");
82  fDoNoiseFiltering = pset.get<bool>("DoNoiseFiltering", true);
83  fNumTicksToDropFront = pset.get<size_t>("NumTicksToDropFront", 2400);
84  fWindowSize = pset.get<size_t>("WindowSize", 6400);
85  }
void art::Modifier::registerProducts ( ProductDescriptions &  productsToRegister)
inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

17  {
18  ProductRegistryHelper::registerProducts(productsToRegister,
19  moduleDescription());
20  }
art::ProductRegistryHelper::registerProducts
void registerProducts(ProductDescriptions &productsToRegister, ModuleDescription const &md)
Definition: ProductRegistryHelper.cc:24
art::ModuleBase::moduleDescription
ModuleDescription const & moduleDescription() const
Definition: ModuleBase.cc:13
void art::ModuleBase::setModuleDescription ( ModuleDescription const &  md)
inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

32  {
33  md_ = md;
34  }
art::ModuleBase::md_
std::optional< ModuleDescription > md_
Definition: ModuleBase.h:55
void art::ModuleBase::sortConsumables ( std::string const &  current_process_name)
inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

50  {
51  // Now that we know we have seen all the consumes declarations,
52  // sort the results for fast lookup later.
53  collector_.sortConsumables(current_process_name);
54  }
art::ModuleBase::collector_
ConsumesCollector collector_
Definition: ModuleBase.h:56
art::ConsumesCollector::sortConsumables
void sortConsumables(std::string const &current_process_name)
Definition: ConsumesCollector.cc:19

Member Data Documentation

std::string noisefilteralg::WireCellNoiseFilter::fDigitModuleLabel
private

Definition at line 58 of file WireCellNoiseFilter_module.cc.

Referenced by produce(), and reconfigure().

bool noisefilteralg::WireCellNoiseFilter::fDoNoiseFiltering
private

Definition at line 59 of file WireCellNoiseFilter_module.cc.

Referenced by produce(), and reconfigure().

size_t noisefilteralg::WireCellNoiseFilter::fNumTicksToDropFront
private

Definition at line 60 of file WireCellNoiseFilter_module.cc.

Referenced by DoNoiseFilter(), produce(), and reconfigure().

size_t noisefilteralg::WireCellNoiseFilter::fWindowSize
private

Definition at line 61 of file WireCellNoiseFilter_module.cc.

Referenced by DoNoiseFilter(), produce(), and reconfigure().

IDFT::pointer noisefilteralg::WireCellNoiseFilter::m_dft
private

Definition at line 57 of file WireCellNoiseFilter_module.cc.

Referenced by DoNoiseFilter(), and reconfigure().

The documentation for this class was generated from the following file: