icarus_tool::ROIFromDecoder Class Reference

Inheritance diagram for icarus_tool::ROIFromDecoder:

Public Member Functions
	ROIFromDecoder (const fhicl::ParameterSet &pset)
	~ROIFromDecoder ()
void	configure (const fhicl::ParameterSet &pset) override
void	initializeHistograms (art::TFileDirectory &) override
void	FindROIs (const art::Event &, const ArrayFloat &, const std::vector< raw::ChannelID_t > &, const geo::PlaneID &, ArrayFloat &, ArrayBool &) override
Public Member Functions inherited from icarus_tool::IROILocator
virtual	~IROILocator () noexcept=default

Private Types
using	TPCIDToLabelMap = std::map< geo::TPCID, std::string >

Private Attributes
TPCIDToLabelMap	fTPCIDToLabelMap
	Translate from TPCID to a substring. More...
std::vector< art::InputTag >	fROILabelVec
	List of input files to search. More...
const geo::GeometryCore *	fGeometry = lar::providerFrom<geo::Geometry>()

Additional Inherited Members
Public Types inherited from icarus_tool::IROILocator
using	VectorBool = std::vector< bool >
using	VectorFloat = std::vector< float >
using	ArrayBool = std::vector< VectorBool >
using	ArrayFloat = std::vector< VectorFloat >
using	PlaneIDVec = std::vector< geo::PlaneID >

Detailed Description

Definition at line 27 of file ROIFromDecoder_tool.cc.

Member Typedef Documentation

using icarus_tool::ROIFromDecoder::TPCIDToLabelMap = std::map<geo::TPCID,std::string>

private

Definition at line 41 of file ROIFromDecoder_tool.cc.

Constructor & Destructor Documentation

icarus_tool::ROIFromDecoder::ROIFromDecoder ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 53 of file ROIFromDecoder_tool.cc.

{
    std::cout << "*** In ROIFromDecoder constructor ***" << std::endl;
    configure(pset);
}

icarus_tool::ROIFromDecoder::~ROIFromDecoder

(

)

Definition at line 59 of file ROIFromDecoder_tool.cc.

{
}

Member Function Documentation

void icarus_tool::ROIFromDecoder::configure ( const fhicl::ParameterSet &  pset )

overridevirtual

Implements icarus_tool::IROILocator.

Definition at line 63 of file ROIFromDecoder_tool.cc.

{
    // Start by recovering the parameters
    fROILabelVec = pset.get<std::vector<art::InputTag>>("ROILabelVec", {""});

    // Build our map
    fTPCIDToLabelMap[geo::TPCID(0,0)] = "EE";
    fTPCIDToLabelMap[geo::TPCID(0,1)] = "EE";
    fTPCIDToLabelMap[geo::TPCID(0,2)] = "EW";
    fTPCIDToLabelMap[geo::TPCID(0,3)] = "EW";
    fTPCIDToLabelMap[geo::TPCID(1,0)] = "WE";
    fTPCIDToLabelMap[geo::TPCID(1,1)] = "WE";
    fTPCIDToLabelMap[geo::TPCID(1,2)] = "WW";
    fTPCIDToLabelMap[geo::TPCID(1,3)] = "WW";

    return;
}

void icarus_tool::ROIFromDecoder::FindROIs ( const art::Event &  event, const ArrayFloat &  inputImage, const std::vector< raw::ChannelID_t > &  channelVec, const geo::PlaneID &  planeID, ArrayFloat &  output, ArrayBool &  outputROIs  )

overridevirtual

Implements icarus_tool::IROILocator.

Definition at line 81 of file ROIFromDecoder_tool.cc.

{
    // First thing is find the correct data product to recover ROIs
    TPCIDToLabelMap::const_iterator tpcItr = fTPCIDToLabelMap.find(planeID.asTPCID());

    if (tpcItr == fTPCIDToLabelMap.end())
    {
        std::cout << "Can't happen? for planeID: " << planeID << std::endl;
        return;
    }

    art::InputTag roiLabel("");

    for(const auto& tag : fROILabelVec)
    {
        if (tag.instance().find(tpcItr->second) != std::string::npos)
        {
            roiLabel = tag;
            break;
        }
    }

    if (roiLabel != "")
    {
        // do stuff here
    
        // Read in the collection of full length deconvolved waveforms
        // Note we assume this list is sorted in increasing channel number!
        art::Handle< std::vector<recob::Wire>> wireVecHandle;
        
        event.getByLabel(roiLabel, wireVecHandle);

        // The first task is to build a mapping between channel and ROIs from the input data product. 
        // The input data product will contain information for every plane in a physical TPC... unfortunately, we only want a single plane within
        // a logical TPC... so we need to loop through to find what we want
        std::unordered_map<raw::ChannelID_t,const recob::Wire*> channelToWireMap;

        for(const auto& wireData : *wireVecHandle)
        {
            std::vector<geo::WireID> wireIDVec = fGeometry->ChannelToWire(wireData.Channel());

            for(const auto& wireID : wireIDVec)
            {
                if (wireID.asPlaneID() != planeID) continue;

                if (wireID.Wire > outputROIs.size())
                {
                    std::cout << "#################################### Wire out of bounds! Wire: " << wireID.Wire << ", max: " << outputROIs.size() << " #################" << std::endl;
                    continue;
                }

                channelToWireMap[wireData.Channel()] = &wireData;
            }
        }

        // Now we loop through the input image channels and make the new ROIs
        for(size_t channelIdx = 0; channelIdx < channelVec.size(); channelIdx++)
        {
            raw::ChannelID_t channel = channelVec[channelIdx];

            std::unordered_map<raw::ChannelID_t,const recob::Wire*>::const_iterator wireItr = channelToWireMap.find(channel);

            if (wireItr != channelToWireMap.end())
            {
                const recob::Wire& wireData = *(wireItr->second);

                // Get the WireIDs (again)
                std::vector<geo::WireID> wireIDVec = fGeometry->ChannelToWire(channel);

                for(const auto& wireID : wireIDVec)
                {
                    if (wireID.asPlaneID() != planeID) continue;

                    // Recover this wire's output vector
                    VectorBool& channelData = outputROIs[wireID.Wire];

                    // What we need to do is find the ROIs in the input wire data and then translate to the output
                    const recob::Wire::RegionsOfInterest_t signalROIs = wireData.SignalROI();

                    for(const auto& range : signalROIs.get_ranges())
                    {
                        size_t startTick = range.begin_index();
                        size_t roiLen    = range.data().size();
                        size_t stopTick  = startTick + roiLen;

                        if (startTick > channelData.size())
                        {
                            std::cout << "*** ROI decoder has start tick larger than output array, start: " << startTick << ", array size: " << channelData.size() << std::endl;
                            continue;
                        }

                        if (stopTick > channelData.size())
                        {
                            std::cout << "*** ROI decoder has ROI length larger than output array, start: " << startTick << ", end: " << stopTick << ", array size: " << channelData.size() << std::endl;
                            continue;
                        }

                        std::fill(channelData.begin() + startTick, channelData.begin() + stopTick, true);
                    }

                    if (planeID.Cryostat == 0 && planeID.TPC == 0 && planeID.Plane == 0) std::cout << std::endl;
                }
            }
            else
            {
                std::cout << "*** Not finding ROI information for channel: " << channel << std::endl;
            }
        }
    }

    return;
}

void icarus_tool::ROIFromDecoder::initializeHistograms ( art::TFileDirectory &  )

inlineoverridevirtual

Implements icarus_tool::IROILocator.

Definition at line 35 of file ROIFromDecoder_tool.cc.

{return;}

Member Data Documentation

const geo::GeometryCore* icarus_tool::ROIFromDecoder::fGeometry = lar::providerFrom<geo::Geometry>()

private

Definition at line 48 of file ROIFromDecoder_tool.cc.

std::vector<art::InputTag> icarus_tool::ROIFromDecoder::fROILabelVec

private

List of input files to search.

Definition at line 46 of file ROIFromDecoder_tool.cc.

TPCIDToLabelMap icarus_tool::ROIFromDecoder::fTPCIDToLabelMap

private

Translate from TPCID to a substring.

Definition at line 43 of file ROIFromDecoder_tool.cc.

---

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

• ROIFromDecoder_tool.cc