icarus_tool::BaselineStandard Class Reference

Inheritance diagram for icarus_tool::BaselineStandard:

Public Member Functions
	BaselineStandard (const fhicl::ParameterSet &pset)
	~BaselineStandard ()
void	configure (const fhicl::ParameterSet &pset) override
void	outputHistograms (art::TFileDirectory &) const override
float	GetBaseline (std::vector< float > const &, raw::ChannelID_t, size_t, size_t) const override
Public Member Functions inherited from icarus_tool::IBaseline
virtual	~IBaseline () noexcept=default

Private Attributes
int	fNumBinsToAverage
art::ServiceHandle < icarusutil::SignalShapingICARUSService >	fSignalShaping

Detailed Description

Definition at line 24 of file BaselineStandard_tool.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 45 of file BaselineStandard_tool.cc.

{
    configure(pset);
}

icarus_tool::BaselineStandard::~BaselineStandard

(

)

Definition at line 50 of file BaselineStandard_tool.cc.

{
}

Member Function Documentation

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

overridevirtual

Implements icarus_tool::IBaseline.

Definition at line 54 of file BaselineStandard_tool.cc.

{
    // Start by recovering the parameters
    fNumBinsToAverage = pset.get<int>("NumBinsToAverage", 20);
    
    // Get signal shaping service.
    fSignalShaping = art::ServiceHandle<icarusutil::SignalShapingICARUSService>();
    
    return;
}

float icarus_tool::BaselineStandard::GetBaseline ( std::vector< float > const &  holder, raw::ChannelID_t  channel, size_t  roiStart, size_t  roiLen  ) const

overridevirtual

Implements icarus_tool::IBaseline.

Definition at line 66 of file BaselineStandard_tool.cc.

{
    float base=0;
    //1. Check Baseline match?
    // If not, include next ROI(if none, go to the end of signal)
    // If yes, proceed
    size_t nBinsToAve(fNumBinsToAverage);
    size_t roiStop(roiStart + roiLen);
    size_t newRoiStart(roiStart);
    size_t newRoiStop(roiStop);
    size_t nTries(4);
    
    // Calculate baslines from the very front of the deconvolution buffer and from the end
    float  basePre  = std::accumulate(holder.begin() + roiStart, holder.begin() + roiStart + nBinsToAve, 0.) / float(nBinsToAve);
    float  basePost = std::accumulate(holder.begin() + roiStop - nBinsToAve, holder.begin() + roiStop,   0.) / float(nBinsToAve);
    
    // emulate method for refining baseline from original version of CalWireROI
    float deconNoise = 1.26491 * fSignalShaping->GetDeconNoise(channel);    // 4./sqrt(10) * noise
    
    // If the estimated baseline from the front of the roi does not agree well with that from the end
    // of the roi then we'll extend the roi hoping for good agreement
    while(!(fabs(basePre - basePost) < deconNoise) && nTries++ < 3)
    {
        size_t nBinsToAdd(10);
        
        if (newRoiStart < nBinsToAdd)                newRoiStart  = 0;
        else                                         newRoiStart -= nBinsToAdd;
        if (newRoiStop + nBinsToAdd > holder.size()) newRoiStop   = holder.size();
        else                                         newRoiStop  += nBinsToAdd;
        
        basePre  = std::accumulate(holder.begin() + newRoiStart, holder.begin() + newRoiStart + nBinsToAve, 0.) / float(nBinsToAve);
        basePost = std::accumulate(holder.begin() + newRoiStop - nBinsToAve, holder.begin() + newRoiStop,   0.) / float(nBinsToAve);
    }
    
    // get spread in "pre" baseline
    float maxPre = *std::max_element(holder.begin() + roiStart, holder.begin() + roiStart + nBinsToAve);
    float minPre = *std::min_element(holder.begin() + roiStart, holder.begin() + roiStart + nBinsToAve);
    
    // Basically, we are hoping for a relatively smooth "front porch" for the waveform
    if ((maxPre - minPre) < 4.) base = basePre;
    else
    {
        // No success... see if the "back porch" is smooth
        size_t roiStop = roiStart + roiLen;
        
        float maxPost = *std::max_element(holder.begin() + roiStop - nBinsToAve, holder.begin() + roiStop);
        float minPost = *std::min_element(holder.begin() + roiStop - nBinsToAve, holder.begin() + roiStop);
        
        if ((maxPost - minPost) < 4.) base = basePost;
        // Starting to get desparate...
        else if ((maxPre  - minPre)  < 8. && std::fabs(basePre)  < std::fabs(basePost)) base = basePre;
        else if ((maxPost - minPost) < 8. && std::fabs(basePost) < std::fabs(basePre))  base = basePost;
        // Ok, apply brute force
        else
        {
            float min = *std::min_element(holder.begin()+roiStart,holder.begin()+roiStart+roiLen);
            float max = *std::max_element(holder.begin()+roiStart,holder.begin()+roiStart+roiLen);
            int nbin = std::ceil(max - min);
            if (nbin > 0){
                TH1F *h1 = new TH1F("h1","h1",nbin,min,max);
                for (unsigned int bin = roiStart; bin < roiStart+roiLen; bin++){
                    h1->Fill(holder[bin]);
                }
                int   pedBin = h1->GetMaximumBin();
                float ped    = h1->GetBinCenter(pedBin);
                float ave=0,ncount = 0;
                for (unsigned int bin = roiStart; bin < roiStart+roiLen; bin++){
                    if (fabs(holder[bin]-ped)<2){
                        ave +=holder[bin];
                        ncount ++;
                    }
                }
                if (ncount==0) ncount=1;
                ave = ave/ncount;
                h1->Delete();
                base = ave;
            }
        }
    }
    
    return base;
}

void icarus_tool::BaselineStandard::outputHistograms ( art::TFileDirectory &  histDir ) const

overridevirtual

Implements icarus_tool::IBaseline.

Definition at line 152 of file BaselineStandard_tool.cc.

{
    // It is assumed that the input TFileDirectory has been set up to group histograms into a common
    // folder at the calling routine's level. Here we create one more level of indirection to keep
    // histograms made by this tool separate.
/*
    std::string dirName = "BaselinePlane_" + std::to_string(fPlane);
    
    art::TFileDirectory dir = histDir.mkdir(dirName.c_str());
    
    auto const* detprop      = lar::providerFrom<detinfo::DetectorPropertiesService>();
    double      samplingRate = detprop->SamplingRate();
    double      numBins      = fBaselineVec.size();
    double      maxFreq      = 500. / samplingRate;
    std::string histName     = "BaselinePlane_" + std::to_string(fPlane);
    
    TH1D*       hist         = dir.make<TH1D>(histName.c_str(), "Baseline;Frequency(MHz)", numBins, 0., maxFreq);
    
    for(int bin = 0; bin < numBins; bin++)
    {
        double freq = maxFreq * double(bin + 0.5) / double(numBins);
        
        hist->Fill(freq, fBaselineVec.at(bin).Re());
    }
*/
    
    return;
}

Member Data Documentation

int icarus_tool::BaselineStandard::fNumBinsToAverage

private

Definition at line 38 of file BaselineStandard_tool.cc.

art::ServiceHandle<icarusutil::SignalShapingICARUSService> icarus_tool::BaselineStandard::fSignalShaping

private

Definition at line 40 of file BaselineStandard_tool.cc.

---

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

• BaselineStandard_tool.cc