Inheritance diagram for evdb_tool::DrawRawHist:

Public Member Functions
	DrawRawHist (const fhicl::ParameterSet &pset)

void	configure (const fhicl::ParameterSet &pset) override

void	Fill (evdb::View2D &, raw::ChannelID_t &, float, float) override

void	Draw (const std::string &, float, float) override

float	getMaximum () const override

float	getMinimum () const override

Private Member Functions
void	BookHistogram (raw::ChannelID_t &, float, float)

Private Attributes
float	fMaximum

float	fMinimum

std::unique_ptr< TH1F >	fRawDigitHist

Detailed Description

Definition at line 27 of file DrawRawHist_tool.cc.

Constructor & Destructor Documentation

evdb_tool::DrawRawHist::DrawRawHist ( const fhicl::ParameterSet & pset )

explicit

Definition at line 49 of file DrawRawHist_tool.cc.

References configure().

   {
     configure(pset);
   }

Member Function Documentation

void evdb_tool::DrawRawHist::BookHistogram	(	raw::ChannelID_t &	channel,
		float	startTick,
		float	numTicks
	)

private

Definition at line 142 of file DrawRawHist_tool.cc.

References DEFINE_ART_CLASS_TOOL, fRawDigitHist, evd::ColorDrawingOptions::fRawQHigh, evd::ColorDrawingOptions::fRawQLow, and Get.

Referenced by Fill(), and getMinimum().

   {
     art::ServiceHandle<evd::ColorDrawingOptions const> cst;
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
 
     // Get rid of the previous histograms
     if (fRawDigitHist.get()) fRawDigitHist.reset();
 
     // figure out the signal type for this plane, assume that
     // plane n in each TPC/cryostat has the same type
     geo::SigType_t sigType = wireReadoutGeom.SignalType(channel);
     int numBins = numTicks;
 
     fRawDigitHist = std::make_unique<TH1F>(
       "fRAWQHisto", ";t [ticks];q [ADC]", numBins, startTick, startTick + numTicks);
 
     TH1F* histPtr = fRawDigitHist.get();
 
     histPtr->SetMaximum(cst->fRawQHigh[(size_t)sigType]);
     histPtr->SetMinimum(cst->fRawQLow[(size_t)sigType]);
 
     histPtr->SetLineColor(kBlack);
     histPtr->SetLineWidth(1);
 
     histPtr->GetXaxis()->SetLabelSize(0.10);   // was 0.15
     histPtr->GetXaxis()->SetLabelOffset(0.01); // was 0.00
     histPtr->GetXaxis()->SetTitleSize(0.10);   // was 0.15
     histPtr->GetXaxis()->SetTitleOffset(0.60); // was 0.80
 
     histPtr->GetYaxis()->SetLabelSize(0.10);    // was 0.15
     histPtr->GetYaxis()->SetLabelOffset(0.002); // was 0.00
     histPtr->GetYaxis()->SetTitleSize(0.10);    // was 0.15
     histPtr->GetYaxis()->SetTitleOffset(0.16);  // was 0.80
   }

void evdb_tool::DrawRawHist::configure ( const fhicl::ParameterSet & pset )

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 54 of file DrawRawHist_tool.cc.

Referenced by DrawRawHist().

54 {}

void evdb_tool::DrawRawHist::Draw	(	const std::string &	options,
		float	maxLowVal,
		float	maxHiVal
	)

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 130 of file DrawRawHist_tool.cc.

References fRawDigitHist.

   {
     TH1F* histPtr = fRawDigitHist.get();
 
     // Do we have valid limits to set?
     histPtr->SetMaximum(maxHiVal);
     histPtr->SetMinimum(maxLowVal);
 
     histPtr->Draw(options.c_str());
   }

void evdb_tool::DrawRawHist::Fill	(	evdb::View2D &	view2D,
		raw::ChannelID_t &	channel,
		float	lowBin,
		float	numTicks
	)

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 56 of file DrawRawHist_tool.cc.

References raw::RawDigit::ADCs(), BookHistogram(), raw::RawDigit::Channel(), raw::RawDigit::Compression(), fMaximum, fMinimum, evd::RawDrawingOptions::fPedestalOption, evd::RawDrawingOptions::fRawDataLabels, fRawDigitHist, evdb::EventHolder::GetEvent(), raw::RawDigit::GetPedestal(), evdb::EventHolder::Instance(), art::Handle< T >::isValid(), raw::RawDigit::Samples(), and raw::Uncompress().

   {
     art::ServiceHandle<evd::RawDrawingOptions const> rawOpt;
 
     //grab the singleton with the event
     const art::Event* event = evdb::EventHolder::Instance()->GetEvent();
     if (!event) return;
 
     // Handle histograms
     BookHistogram(channel, lowBin, numTicks);
 
     fMinimum = std::numeric_limits<float>::max();
     fMaximum = std::numeric_limits<float>::lowest();
 
     // Loop over the possible producers of RawDigits
     for (const auto& rawDataLabel : rawOpt->fRawDataLabels) {
       art::Handle<std::vector<raw::RawDigit>> rawDigitVecHandle;
       event->getByLabel(rawDataLabel, rawDigitVecHandle);
 
       if (!rawDigitVecHandle.isValid()) continue;
 
       for (size_t rawDigitIdx = 0; rawDigitIdx < rawDigitVecHandle->size(); rawDigitIdx++) {
         art::Ptr<raw::RawDigit> rawDigit(rawDigitVecHandle, rawDigitIdx);
 
         if (rawDigit->Channel() != channel) continue;
 
         // We will need the pedestal service...
         const lariov::DetPedestalProvider& pedestalRetrievalAlg =
           art::ServiceHandle<lariov::DetPedestalService const>()->GetPedestalProvider();
 
         // recover the pedestal
         float pedestal = 0;
 
         if (rawOpt->fPedestalOption == 0) { pedestal = pedestalRetrievalAlg.PedMean(channel); }
         else if (rawOpt->fPedestalOption == 1) {
           pedestal = rawDigit->GetPedestal();
         }
         else if (rawOpt->fPedestalOption == 2) {
           pedestal = 0;
         }
         else {
           mf::LogWarning("DrawRawHist")
             << " PedestalOption is not understood: " << rawOpt->fPedestalOption
             << ".  Pedestals not subtracted.";
         }
 
         std::vector<short> uncompressed(rawDigit->Samples());
         raw::Uncompress(rawDigit->ADCs(), uncompressed, rawDigit->Compression());
 
         TH1F* histPtr = fRawDigitHist.get();
 
         for (size_t idx = 0; idx < uncompressed.size(); idx++) {
           float signalVal = float(uncompressed[idx]) - pedestal;
 
           histPtr->Fill(float(idx) + 0.5, signalVal);
         }
 
         short minimumVal = *std::min_element(uncompressed.begin(), uncompressed.end());
         short maximumVal = *std::max_element(uncompressed.begin(), uncompressed.end());
 
         fMinimum = float(minimumVal) - pedestal;
         fMaximum = float(maximumVal) - pedestal;
 
         histPtr->SetLineColor(kBlack);
 
         // There is only one channel displayed so if here we are done
         break;
       }
     }
   }