Inheritance diagram for evdb_tool::DrawWireHist:

Public Member Functions
	DrawWireHist (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::vector< int >	fColorMap

std::unordered_map< std::string, std::unique_ptr< TH1F > >	fRecoHistMap

Detailed Description

Definition at line 24 of file DrawWireHist_tool.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 47 of file DrawWireHist_tool.cc.

References configure().

   {
     configure(pset);
   }

Member Function Documentation

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

private

Definition at line 130 of file DrawWireHist_tool.cc.

References DEFINE_ART_CLASS_TOOL, fRecoHistMap, evd::ColorDrawingOptions::fRecoQHigh, evd::ColorDrawingOptions::fRecoQLow, evd::RecoDrawingOptions::fWireLabels, and Get.

Referenced by Fill(), and getMinimum().

   {
     art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
     art::ServiceHandle<evd::ColorDrawingOptions const> cst;
     art::ServiceHandle<evd::RawDrawingOptions const> drawopt;
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
 
     // Get rid of the previous histograms
     fRecoHistMap.clear();
 
     // Now add a histogram for each of the wire labels
     for (auto& tag : recoOpt->fWireLabels) {
       // 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);
       std::string tagString(tag.encode());
       int numBins = numTicks;
 
       fRecoHistMap[tagString] = std::make_unique<TH1F>(
         "fCALTQHisto", ";t [ticks];q [ADC]", numBins, startTick, startTick + numTicks);
 
       TH1F* histPtr = fRecoHistMap.at(tagString).get();
 
       histPtr->SetMaximum(cst->fRecoQHigh[(size_t)sigType]);
       histPtr->SetMinimum(cst->fRecoQLow[(size_t)sigType]);
 
       histPtr->SetLineColor(kBlue);
       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::DrawWireHist::configure ( const fhicl::ParameterSet & pset )

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 52 of file DrawWireHist_tool.cc.

References fColorMap, and fRecoHistMap.

Referenced by DrawWireHist().

   {
     fColorMap.push_back(kBlue);
     fColorMap.push_back(kMagenta);
     fColorMap.push_back(kBlack);
     fColorMap.push_back(kRed);
 
     fRecoHistMap.clear();
   }

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

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 116 of file DrawWireHist_tool.cc.

References fRecoHistMap.

   {
     for (const auto& histMap : fRecoHistMap) {
       TH1F* histPtr = histMap.second.get();
 
       // Set the limits
       histPtr->SetMaximum(maxHiVal);
       histPtr->SetMinimum(maxLowVal);
 
       histPtr->Draw(options.c_str());
     }
   }

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

overridevirtual

Implements evdb_tool::IWaveformDrawer.

Definition at line 62 of file DrawWireHist_tool.cc.

References BookHistogram(), recob::Wire::Channel(), art::InputTag::encode(), fColorMap, evd::RawDrawingOptions::fDrawRawDataOrCalibWires, fMaximum, fMinimum, fRecoHistMap, evd::RecoDrawingOptions::fWireLabels, evdb::EventHolder::GetEvent(), evdb::EventHolder::Instance(), and recob::Wire::Signal().

   {
     art::ServiceHandle<evd::RawDrawingOptions const> rawOpt;
     art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
 
     // Check if we're supposed to draw raw hits at all
     if (rawOpt->fDrawRawDataOrCalibWires == 0) return;
 
     //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();
 
     int nWireLabels = 0;
     for (size_t imod = 0; imod < recoOpt->fWireLabels.size(); ++imod) {
       // Step one is to recover the hits for this label that match the input channel
       art::InputTag const which = recoOpt->fWireLabels[imod];
 
       art::Handle<std::vector<recob::Wire>> wireVecHandle;
       if (!event->getByLabel(which, wireVecHandle)) continue;
       ++nWireLabels;
 
       for (size_t wireIdx = 0; wireIdx < wireVecHandle->size(); wireIdx++) {
         art::Ptr<recob::Wire> wire(wireVecHandle, wireIdx);
 
         if (wire->Channel() != channel) continue;
 
         const std::vector<float>& signalVec = wire->Signal();
 
         TH1F* histPtr = fRecoHistMap.at(which.encode()).get();
 
         for (size_t idx = 0; idx < signalVec.size(); idx++) {
           histPtr->Fill(float(idx) + 0.5, signalVec[idx]);
 
           fMinimum = std::min(fMinimum, signalVec[idx]);
           fMaximum = std::max(fMaximum, signalVec[idx]);
         }
 
         histPtr->SetLineColor(fColorMap.at((nWireLabels - 1) % recoOpt->fWireLabels.size()));
 
         // There is only one channel displayed so if here we are done
         break;
       }
     } //end loop over HitFinding modules
   }