lbne::FilterWF Class Reference

Inheritance diagram for lbne::FilterWF:

Public Member Functions
	FilterWF (fhicl::ParameterSet const &pset)
virtual	~FilterWF ()
virtual void	produce (art::Event &evt) override
void	reconfigure (fhicl::ParameterSet const &pset)
void	beginRun (art::Run &run) override
void	beginJob () override
void	endRun (art::Run &run) override
Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDProducer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept
	Producer (fhicl::ParameterSet const &)
	Producer (Producer const &)=delete
	Producer (Producer &&)=delete
Producer &	operator= (Producer const &)=delete
Producer &	operator= (Producer &&)=delete
void	doBeginJob ()
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)
Public Member Functions inherited from art::Modifier
	~Modifier () noexcept
	Modifier ()
	Modifier (Modifier const &)=delete
	Modifier (Modifier &&)=delete
Modifier &	operator= (Modifier const &)=delete
Modifier &	operator= (Modifier &&)=delete
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
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)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fRawDigitModuleLabel
std::string	fRawDigitModuleInstance
bool	fSkipStuckCodes
bool	fSkipTicks
unsigned int	fLowerSkipTick
unsigned int	fUpperSkipTick
int	fMethod
art::ServiceHandle< geo::Geometry >	fGeom
art::ServiceHandle< lbne::ChannelMapService >	fChannelMap

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >
Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)
Protected Member Functions inherited from art::ModuleBase
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 ()

Detailed Description

Definition at line 65 of file FilterWF_module.cc.

Constructor & Destructor Documentation

lbne::FilterWF::FilterWF ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 96 of file FilterWF_module.cc.

                                                    : EDProducer{pset} {
  this->reconfigure(pset);
  produces<std::vector<raw::RawDigit> >();
}

lbne::FilterWF::~FilterWF ( )

virtual

Definition at line 102 of file FilterWF_module.cc.

                       {
}

Member Function Documentation

void lbne::FilterWF::beginJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 123 of file FilterWF_module.cc.

                            {
  art::ServiceHandle<art::TFileService> tfs;
  /*
  int NHisto =0;
  for (int Plane=0; Plane<3; Plane++) {
    for (int Reg=0; Reg < 128; ++Reg) {
      ++NHisto;
      std::stringstream oss;
      oss << "Channel_vs_tick_" << Plane << "_" << Reg;
      std::string Title = oss.str();
      std::cout << oss.str() << " " << Title << std::endl;
      Channel_vs_tick[Plane][Reg] = tfs->make<TH2D>(Title.c_str(),"Plot of Channel vs Tick; Tick; Channel", 6000, 0., 6000., 14, -0.5, 13.5);
    }
  }
  */
}

void lbne::FilterWF::beginRun ( art::Run &  run )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 119 of file FilterWF_module.cc.

                                        {
  return;
}

void lbne::FilterWF::endRun ( art::Run &  run )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 141 of file FilterWF_module.cc.

                                      {
  return;
}

void lbne::FilterWF::produce ( art::Event &  evt )

overridevirtual

Implements art::EDProducer.

Definition at line 146 of file FilterWF_module.cc.

                                        {

  art::Handle<std::vector<raw::RawDigit> > rawDigitHandle;
  evt.getByLabel(fRawDigitModuleLabel, fRawDigitModuleInstance, rawDigitHandle);
  std::vector<raw::RawDigit> const& rawDigitVector(*rawDigitHandle);
  
  // define temporary vector to hold filtered waveforms
  const unsigned int n_channels = fGeom->Nchannels();
  //const unsigned int n_channels = fGeom->rawDigitVector.size();
  std::vector<std::vector<short> > filterWf(n_channels); // MW: I had to make this a vector of shorts to make it compile (the RawDigit constructor does not accept ADC as float)
  
  // raw digit map
  std::map<int,raw::RawDigit> rawDigitMap;
  std::map<int,float> pedestalMap;
  unsigned int maxNumBins = 0;
  for (std::vector<raw::RawDigit>::const_iterator digitIt = rawDigitVector.begin(); digitIt != rawDigitVector.end(); ++digitIt) {
    rawDigitMap[digitIt->Channel()] = *digitIt;
    pedestalMap[digitIt->Channel()] = digitIt->GetPedestal();
    if (digitIt->NADC() > maxNumBins) maxNumBins = digitIt->NADC();
  }
  //std::cout<<"Method = "<<fMethod<<std::endl;
  // define set of induction and collection channels in each regulator group
  std::vector<std::vector<std::vector<unsigned int> > > Chs;  //Chs[plane id][group id][channel id]

  if (fMethod == 0){ //group by regulators
    Chs.resize(3);
    for (size_t i = 0; i<3; ++i){
      Chs[i].resize(32);
    }
    for (unsigned int i = 0; i<n_channels; ++i){//online channels
      unsigned int plane     = fChannelMap->PlaneFromOnlineChannel(i);
      unsigned int rce       = fChannelMap->RCEFromOnlineChannel(i);
      unsigned int regulator = fChannelMap->RegulatorFromOnlineChannel(i);
      Chs[plane][rce*2+regulator].push_back(i);
    }
  }
  else if (fMethod == 1){
    Chs.resize(3);
    for (size_t i = 0; i<3; ++i){
      Chs[i].resize(128);
    }
    for (unsigned int i = 0; i<n_channels; ++i){//online channels
      unsigned int plane     = fChannelMap->PlaneFromOnlineChannel(i);
      unsigned int rce       = fChannelMap->RCEFromOnlineChannel(i);
      unsigned int asic      = fChannelMap->ASICFromOnlineChannel(i);
      //std::cout<<i<<" "<<plane<<" "<<rce<<" "<<asic<<std::endl;
      Chs[plane][rce*8+asic].push_back(i);
    }
  }   
  // derive correction factors - require raw adc waveform and pedestal for each channel
  std::vector<Double_t> corrVals;
  // loop through time slices
  for (unsigned int s = 0; s < maxNumBins; s++) {
    for (size_t i = 0; i<Chs.size(); ++i){
      for (size_t j = 0; j<Chs[i].size(); ++j){
        //if (s==0) std::cout << "Looking at " << s << " " << i << " " << j << " " << Chs[i][j].size() << std::endl;
        corrVals.clear();
        for (size_t k = 0; k<Chs[i][j].size(); ++k){
          
          unsigned int offlineChan = fChannelMap->Offline(Chs[i][j][k]);
          if (rawDigitMap.count(offlineChan) == 0)
            continue;
          int adc = rawDigitMap.at(offlineChan).ADC(s);
          if ( fSkipStuckCodes && ( (adc & 0x3F) == 0x0 || (adc & 0x3F) == 0x3F ) ) 
            continue;
          if (adc < 10) //skip "sample dropping" problem
            continue;
          double mean = pedestalMap.at(offlineChan);
          if (mean < 10)
            continue;
          corrVals.push_back(adc - mean);
        }

        unsigned int corrValSize = corrVals.size();
        sort(corrVals.begin(),corrVals.end());
        double correction = 0;
        if (corrValSize < 2)
          correction = 0.0;
        else if ((corrValSize % 2) == 0)
          correction = (corrVals[corrValSize/2] + corrVals[(corrValSize/2)-1])/2.0;
        else
          correction = corrVals[(corrValSize-1)/2];

        for (size_t k = 0; k<Chs[i][j].size(); ++k){
                  
          unsigned int offlineChan = fChannelMap->Offline(Chs[i][j][k]);
          if (rawDigitMap.count(offlineChan) == 0)
            continue;
          int adc = rawDigitMap.at(offlineChan).ADC(s);
          double newAdc = adc - correction;
          if ( fSkipStuckCodes && ( (adc & 0x3F) == 0x0 || (adc & 0x3F) == 0x3F ) )
            newAdc = adc; //don't do anything about stuck code, will run stuck code removal later
          if ( fSkipTicks and s > fLowerSkipTick and s < fUpperSkipTick )
            newAdc = adc;
          // if the code unsticker is run first, then don't skip the stuck codes.
          //      if( adc < 10  ) //skip "sample dropping" problem
          //        newAdc = 0;
          filterWf.at(offlineChan).push_back(newAdc);
          //if (newAdc < 1000)
          //Channel_vs_tick[i][j] -> SetBinContent( Channel_vs_tick[i][j]->FindBin(s,k), newAdc);
          //if (s==0)
          //std::cout << "Set Channel_vs_tick["<<i<<"]["<<j<<"] bin (" << s << ", " << k << ") to " << Channel_vs_tick[i][j]->GetBinContent( Channel_vs_tick[i][j]->FindBin(s, k) ) << " it should have been " << newAdc << std::endl;
        }
      }//all groups
    }//all planes
  }// all ticks

  // loop over channels - save filtered waveforms into digits
  std::vector<raw::RawDigit> filterRawDigitVector;
  for (unsigned int ich = 0; ich < n_channels; ich++) {
    if (rawDigitMap.count(ich) == 0)
      continue;
    raw::RawDigit theRawDigit(ich, filterWf.at(ich).size(), filterWf.at(ich));
    theRawDigit.SetPedestal(rawDigitMap[ich].GetPedestal(),
                            rawDigitMap[ich].GetSigma());
    filterRawDigitVector.push_back(theRawDigit);
  }

  // save filtered waveforms to event
  evt.put(std::make_unique<decltype(filterRawDigitVector)>(std::move(filterRawDigitVector)));
}

void lbne::FilterWF::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 106 of file FilterWF_module.cc.

                                                            {
  fRawDigitModuleLabel = pset.get<std::string>("RawDigitModuleLabel");
  fRawDigitModuleInstance = pset.get<std::string>("RawDigitModuleInstance");
  //fOutputModuleLabel = "filterwf";
  fSkipStuckCodes = pset.get<bool>("SkipStuckCodes",true);
  fSkipTicks      = pset.get<bool>("SkipTicks",false); // MW: added option to ignore certain tick ranges -- turned off by default
  fLowerSkipTick  = pset.get<int>("LowerSkipTick",0);
  fUpperSkipTick  = pset.get<int>("UpperSkipTick",0);
  fMethod         = pset.get<int>("Method",0);  //0: group by regulator
                                                //1: group by ASIC
}

Member Data Documentation

art::ServiceHandle<lbne::ChannelMapService> lbne::FilterWF::fChannelMap

private

Definition at line 90 of file FilterWF_module.cc.

art::ServiceHandle<geo::Geometry> lbne::FilterWF::fGeom

private

Definition at line 89 of file FilterWF_module.cc.

unsigned int lbne::FilterWF::fLowerSkipTick

private

Definition at line 84 of file FilterWF_module.cc.

int lbne::FilterWF::fMethod

private

Definition at line 86 of file FilterWF_module.cc.

std::string lbne::FilterWF::fRawDigitModuleInstance

private

Definition at line 80 of file FilterWF_module.cc.

std::string lbne::FilterWF::fRawDigitModuleLabel

private

Definition at line 80 of file FilterWF_module.cc.

bool lbne::FilterWF::fSkipStuckCodes

private

Definition at line 82 of file FilterWF_module.cc.

bool lbne::FilterWF::fSkipTicks

private

Definition at line 83 of file FilterWF_module.cc.

unsigned int lbne::FilterWF::fUpperSkipTick

private

Definition at line 84 of file FilterWF_module.cc.

---

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

• dunetpc/dune/daqinput35t/FilterWF_module.cc