butcher::EventButcher Class Reference

Inheritance diagram for butcher::EventButcher:

Public Types
using	Parameters = art::EDProducer::Table< EventButcherConfig >
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 >

Public Member Functions
	EventButcher (Parameters const &params)
virtual	~EventButcher ()
void	produce (art::Event &evt)
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 (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)
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
const EventButcherConfig	m_cfg
art::InputTag	m_rawtag
art::InputTag	m_sigtag

Additional Inherited Members
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 51 of file EventButcher_module.cc.

Member Typedef Documentation

using butcher::EventButcher::Parameters = art::EDProducer::Table<EventButcherConfig>

Definition at line 54 of file EventButcher_module.cc.

Constructor & Destructor Documentation

butcher::EventButcher::EventButcher ( Parameters const &  params )

explicit

Definition at line 76 of file EventButcher_module.cc.

    : EDProducer{params}
    , m_cfg(params())
    , m_rawtag{m_cfg.inRawTag()}
    , m_sigtag{m_cfg.inSigTag()}
//    , m_rawtok{consumes< std::vector<raw::RawDigit> >(m_rawtag)}
//    , m_sigtok{consumes< std::vector<recob::Wire> >(m_sigtag)}
{
    //cerr << "Producing: outraw:"<<m_cfg.outRawTag()<<" outsig:"<<m_cfg.outSigTag()<<" outras:" << m_cfg.outAssnTag() << endl;
    produces< std::vector<raw::RawDigit> >(m_cfg.outRawTag());
    produces< std::vector<recob::Wire> >(m_cfg.outSigTag());
    produces< art::Assns<raw::RawDigit,recob::Wire> >(m_cfg.outAssnTag());
}

butcher::EventButcher::~EventButcher ( )

virtual

Definition at line 90 of file EventButcher_module.cc.

{
}

Member Function Documentation

void butcher::EventButcher::produce ( art::Event &  evt )

virtual

Implements art::EDProducer.

Definition at line 95 of file EventButcher_module.cc.

{

    //cerr <<"In raw=" << m_rawtag << " in sig=" << m_sigtag << "\n";

    art::Handle< std::vector<raw::RawDigit> > raw;
    //event.getByToken(m_rawtok, raw);
    event.getByLabel(m_rawtag, raw);
    const size_t nraw = raw->size();

    art::Handle< std::vector<recob::Wire> > sig;
    //event.getByToken(m_sigtok, sig);
    event.getByLabel(m_sigtag, sig);
    const size_t nsig = sig->size();

    // https://cdcvs.fnal.gov/redmine/projects/art/wiki/The_PtrMaker_utility
    art::PtrMaker<raw::RawDigit> RawPtr(event, m_cfg.outRawTag());
    art::PtrMaker<recob::Wire> SigPtr(event, m_cfg.outSigTag());

    // keep track of raw digit Ptr by its channel id, for later look
    // up in making associations.
    std::unordered_map<raw::ChannelID_t, art::Ptr<raw::RawDigit> > chid2rawptr;

    // raw-signal association
    auto outrsa = std::make_unique< art::Assns<raw::RawDigit,recob::Wire> >();
    auto outraw = std::make_unique< std::vector<raw::RawDigit> >();
    auto outsig = std::make_unique< std::vector<recob::Wire> >();

    const int ndrop = m_cfg.ndrop();
    const int nkeep = m_cfg.nkeep();
    const double sigscale = m_cfg.sigscale();

    // Truncate the raw digits
    for (size_t iraw=0; iraw != nraw; ++iraw) {
        const auto& inrd = raw->at(iraw);
        const auto& inadcs = inrd.ADCs();
        const size_t inlen = inadcs.size();

        const int outlen = std::min(inlen-ndrop, nkeep < 0 ? inlen : nkeep);

        if (outlen <= 0) {
            continue;           // add a "keep_empty" option and save an empty place holder RawDigit here?
        }

        size_t outind = outraw->size();
        const raw::ChannelID_t chid = inrd.Channel();
        raw::RawDigit::ADCvector_t outadc(inadcs.begin()+ndrop, inadcs.begin()+ndrop+outlen);
        outraw->emplace_back(raw::RawDigit(chid, outlen, outadc, raw::kNone));
        // given the truncationn, this is technically a BOGUS thing to do
        auto& outrd = outraw->back();
        outrd.SetPedestal(inrd.GetPedestal(), inrd.GetSigma());

        chid2rawptr[chid] = RawPtr(outind);
    }


    // Truncate the signal and make assns
    for (size_t isig=0; isig != nsig; ++isig) {
        const auto& inw = sig->at(isig);
        std::vector<float> wave = inw.Signal();
        const size_t inlen = wave.size();

        const int outlen = std::min(inlen-ndrop, nkeep < 0 ? inlen : nkeep);

        if (outlen <= 0) {
            continue;
        }

        const auto chid = inw.Channel();
        const auto view = inw.View();

        // resparsify
        recob::Wire::RegionsOfInterest_t roi(outlen);
        auto first = wave.begin()+ndrop;
        auto done = wave.begin()+ndrop+outlen;
        auto beg = first;
        while (true) {
            beg = std::find_if(beg, done, [](float v){return v != 0.0;});
            if (beg == done) {
                break;
            }
            auto end = std::find_if(beg, done, [](float v){return v == 0.0;});

            std::vector<float> scaled(beg, end);
            for (int ind=0; ind<end-beg; ++ind) {
                scaled[ind] *= sigscale;
            }
            roi.add_range(beg-first, scaled.begin(), scaled.end());
            beg = end;
        }

        const size_t outind = outsig->size();
        outsig->emplace_back(recob::Wire(roi, chid, view));

        // associate
        auto rawit = chid2rawptr.find(chid);
        if (rawit == chid2rawptr.end()) {
            continue;           // emit warning about no accompaning raw digit?
        }
        auto const& rawptr = rawit->second;
        auto const sigptr = SigPtr(outind);
        outrsa->addSingle(rawptr, sigptr);
    }

    event.put(std::move(outraw), m_cfg.outRawTag());
    event.put(std::move(outsig), m_cfg.outSigTag());
    event.put(std::move(outrsa), m_cfg.outAssnTag());

}

Member Data Documentation

const EventButcherConfig butcher::EventButcher::m_cfg

private

Definition at line 63 of file EventButcher_module.cc.

art::InputTag butcher::EventButcher::m_rawtag

private

Definition at line 65 of file EventButcher_module.cc.

art::InputTag butcher::EventButcher::m_sigtag

private

Definition at line 65 of file EventButcher_module.cc.

---

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

• larwirecell/larwirecell/Modules/EventButcher_module.cc