tracking::dqm::RawSiStripDQM Class Reference

DQM for silicon strip data read from the DAQ: raw hits, pedestal subtracted hits, the waveforms assembled from them and their pulse fits. More...

#include <RawSiStripDQM.h>

Public Member Functions
	RawSiStripDQM (const std::string &name, framework::Process &process)
void	configure (framework::config::Parameters &parameters) override
	Callback for the EventProcessor to configure itself from the given set of parameters.
void	analyze (const framework::Event &event) override
	Process the event and make histograms or summaries.
Public Member Functions inherited from framework::Analyzer
	Analyzer (const std::string &name, Process &process)
	Class constructor.
virtual void	process (Event &event) final
	Processing an event for an Analyzer is calling analyze.
virtual void	beforeNewRun (ldmx::RunHeader &run_header) final
	Don't allow Analyzers to add parameters to the run header.
Public Member Functions inherited from framework::EventProcessor
	DECLARE_FACTORY (EventProcessor, EventProcessor *, const std::string &, Process &)
	declare that we have a factory for this class
	EventProcessor (const std::string &name, Process &process)
	Class constructor.
virtual	~EventProcessor ()=default
	Class destructor.
virtual void	onNewRun (const ldmx::RunHeader &run_header)
	Callback for the EventProcessor to take any necessary action when the run being processed changes.
virtual void	onFileOpen (EventFile &event_file)
	Callback for the EventProcessor to take any necessary action when a new event input ROOT file is opened.
virtual void	onFileClose (EventFile &event_file)
	Callback for the EventProcessor to take any necessary action when a event input ROOT file is closed.
virtual void	onProcessStart ()
	Callback for the EventProcessor to take any necessary action when the processing of events starts, such as creating histograms.
virtual void	onProcessEnd ()
	Callback for the EventProcessor to take any necessary action when the processing of events finishes, such as calculating job-summary quantities.
template<class T >
const T &	getCondition (const std::string &condition_name)
	Access a conditions object for the current event.
TDirectory *	getHistoDirectory ()
	Access/create a directory in the histogram file for this event processor to create histograms and analysis tuples.
void	setStorageHint (framework::StorageControl::Hint hint)
	Mark the current event as having the given storage control hint from this module_.
void	setStorageHint (framework::StorageControl::Hint hint, const std::string &purposeString)
	Mark the current event as having the given storage control hint from this module and the given purpose string.
int	getLogFrequency () const
	Get the current logging frequency from the process.
int	getRunNumber () const
	Get the run number from the process.
std::string	getName () const
	Get the processor name.
void	createHistograms (const std::vector< framework::config::Parameters > &histos)
	Internal function which is used to create histograms passed from the python configuration @parma histos vector of Parameters that configure histograms to create.

Private Member Functions
std::vector< TH1 * >	getHists (const std::vector< std::string > &names)
	histograms by name, looked up once per event

Private Attributes
std::string	raw_hits_collection_
std::string	subtracted_hits_collection_
std::string	waveforms_collection_
std::string	fitted_hits_collection_
std::string	input_pass_name_
int	n_hybrids_ {0}
std::vector< std::string >	raw_adc_names_
	per-hybrid histogram names, built once
std::vector< std::string >	subtracted_adc_names_
std::vector< std::string >	pchannel_names_
std::vector< std::string >	peak_names_
std::map< int, std::string >	fit_amplitude_names_
	fit amplitude histogram name by layer id

Additional Inherited Members
Protected Member Functions inherited from framework::EventProcessor
void	abortEvent ()
	Abort the event immediately.
Protected Attributes inherited from framework::EventProcessor
HistogramPool	histograms_
	helper object for making and filling histograms
NtupleManager &	ntuple_ {NtupleManager::getInstance()}
	Manager for any ntuples.
logging::logger	the_log_
	The logger for this EventProcessor.

Detailed Description

DQM for silicon strip data read from the DAQ: raw hits, pedestal subtracted hits, the waveforms assembled from them and their pulse fits.

Definition at line 17 of file RawSiStripDQM.h.

Constructor & Destructor Documentation

RawSiStripDQM()

tracking::dqm::RawSiStripDQM::RawSiStripDQM ( const std::string & name, framework::Process & process )

inline

Definition at line 19 of file RawSiStripDQM.h.

      : framework::Analyzer(name, process) {};

Member Function Documentation

analyze()

void tracking::dqm::RawSiStripDQM::analyze ( const framework::Event & event )

overridevirtual

Process the event and make histograms or summaries.

Parameters

event

The Event to analyze

Implements framework::Analyzer.

Definition at line 39 of file RawSiStripDQM.cxx.

                                                       {
  // direct fills, ~1e5 samples per event
  if (event.exists(raw_hits_collection_, input_pass_name_)) {
    const auto& hits{event.getCollection<ldmx::RawSiStripHit>(
        raw_hits_collection_, input_pass_name_)};
    histograms_.fill("n_raw_hits", hits.size());
    auto* read_error{histograms_.get("raw_read_error")};
    const auto raw_adc{getHists(raw_adc_names_)};
    for (const auto& hit : hits) {
      read_error->Fill(hit.getReadError());
      if (hit.getHybridId() >= n_hybrids_) continue;
      for (short adc : hit.getSamples()) {
        raw_adc[hit.getHybridId()]->Fill(adc);
      }
    }
  }
 
  if (event.exists(subtracted_hits_collection_, input_pass_name_)) {
    const auto& hits{event.getCollection<ldmx::RawSiStripHit>(
        subtracted_hits_collection_, input_pass_name_)};
    const auto subtracted_adc{getHists(subtracted_adc_names_)};
    for (const auto& hit : hits) {
      if (hit.getHybridId() >= n_hybrids_) continue;
      for (short adc : hit.getSamples()) {
        subtracted_adc[hit.getHybridId()]->Fill(adc);
      }
    }
  }
 
  if (event.exists(waveforms_collection_, input_pass_name_)) {
    const auto& waveforms{event.getCollection<ldmx::SiStripWaveform>(
        waveforms_collection_, input_pass_name_)};
    histograms_.fill("n_waveforms", waveforms.size());
    for (const auto& waveform : waveforms) {
      histograms_.fill("waveform_n_triggers", waveform.getNTriggers());
      histograms_.fill("waveform_peak_trigger", waveform.peakTrigger());
      if (waveform.getHybridId() >= n_hybrids_) continue;
      histograms_.fill(pchannel_names_[waveform.getHybridId()],
                       waveform.getPchannel());
      histograms_.fill(peak_names_[waveform.getHybridId()],
                       waveform.peakAmplitude());
    }
  }
 
  // fit results live in their own collection, converged fits only
  if (!event.exists(fitted_hits_collection_, input_pass_name_)) return;
  const auto& fitted_hits{event.getCollection<ldmx::FittedSiStripHit>(
      fitted_hits_collection_, input_pass_name_)};
  histograms_.fill("n_fitted_hits", fitted_hits.size());
  for (const auto& hit : fitted_hits) {
    histograms_.fill("fit_t0", hit.getT0());
    if (hit.getNDF() > 0) {
      histograms_.fill("fit_chi2_ndf", hit.getReducedChi2());
    }
    auto name{fit_amplitude_names_.find(hit.getLayerID())};
    if (name != fit_amplitude_names_.end()) {
      histograms_.fill(name->second, hit.getAmplitude());
    }
  }
}

References framework::Event::exists(), framework::HistogramPool::fill(), fit_amplitude_names_, framework::HistogramPool::get(), getHists(), framework::EventProcessor::histograms_, and raw_adc_names_.

configure()

void tracking::dqm::RawSiStripDQM::configure ( framework::config::Parameters & parameters )

overridevirtual

Callback for the EventProcessor to configure itself from the given set of parameters.

The parameters a processor has access to are the member variables of the python class in the sequence that has class_name equal to the EventProcessor class name.

For an example, look at MyProcessor.

Parameters

parameters

Parameters for configuration.

Reimplemented from framework::EventProcessor.

Definition at line 9 of file RawSiStripDQM.cxx.

                                                                   {
  raw_hits_collection_ = parameters.get<std::string>("raw_hits_collection");
  subtracted_hits_collection_ =
      parameters.get<std::string>("subtracted_hits_collection");
  waveforms_collection_ = parameters.get<std::string>("waveforms_collection");
  fitted_hits_collection_ =
      parameters.get<std::string>("fitted_hits_collection");
  input_pass_name_ = parameters.get<std::string>("input_pass_name");
  n_hybrids_ = parameters.get<int>("n_hybrids");
  const auto layer_ids{parameters.get<std::vector<int>>("layer_ids")};
 
  for (int h{0}; h < n_hybrids_; ++h) {
    const std::string suffix{"_h" + std::to_string(h)};
    raw_adc_names_.push_back("raw_adc" + suffix);
    subtracted_adc_names_.push_back("subtracted_adc" + suffix);
    pchannel_names_.push_back("waveform_pchannel" + suffix);
    peak_names_.push_back("waveform_peak" + suffix);
  }
  for (int layer : layer_ids) {
    fit_amplitude_names_[layer] = "fit_amplitude_l" + std::to_string(layer);
  }
}

References fit_amplitude_names_, framework::config::Parameters::get(), and raw_adc_names_.

getHists()

std::vector< TH1 * > tracking::dqm::RawSiStripDQM::getHists ( const std::vector< std::string > & names )

private

histograms by name, looked up once per event

Definition at line 32 of file RawSiStripDQM.cxx.

                                       {
  std::vector<TH1*> hists;
  for (const auto& name : names) hists.push_back(histograms_.get(name));
  return hists;
}

References framework::HistogramPool::get(), and framework::EventProcessor::histograms_.

Referenced by analyze().

Member Data Documentation

fit_amplitude_names_

std::map<int, std::string> tracking::dqm::RawSiStripDQM::fit_amplitude_names_

private

fit amplitude histogram name by layer id

Definition at line 45 of file RawSiStripDQM.h.

Referenced by analyze(), and configure().

fitted_hits_collection_

std::string tracking::dqm::RawSiStripDQM::fitted_hits_collection_

private

Definition at line 35 of file RawSiStripDQM.h.

input_pass_name_

std::string tracking::dqm::RawSiStripDQM::input_pass_name_

private

Definition at line 36 of file RawSiStripDQM.h.

n_hybrids_

int tracking::dqm::RawSiStripDQM::n_hybrids_ {0}

private

Definition at line 37 of file RawSiStripDQM.h.

{0};

pchannel_names_

std::vector<std::string> tracking::dqm::RawSiStripDQM::pchannel_names_

private

Definition at line 42 of file RawSiStripDQM.h.

peak_names_

std::vector<std::string> tracking::dqm::RawSiStripDQM::peak_names_

private

Definition at line 43 of file RawSiStripDQM.h.

raw_adc_names_

std::vector<std::string> tracking::dqm::RawSiStripDQM::raw_adc_names_

private

per-hybrid histogram names, built once

Definition at line 40 of file RawSiStripDQM.h.

Referenced by analyze(), and configure().

raw_hits_collection_

std::string tracking::dqm::RawSiStripDQM::raw_hits_collection_

private

Definition at line 32 of file RawSiStripDQM.h.

subtracted_adc_names_

std::vector<std::string> tracking::dqm::RawSiStripDQM::subtracted_adc_names_

private

Definition at line 41 of file RawSiStripDQM.h.

subtracted_hits_collection_

std::string tracking::dqm::RawSiStripDQM::subtracted_hits_collection_

private

Definition at line 33 of file RawSiStripDQM.h.

waveforms_collection_

std::string tracking::dqm::RawSiStripDQM::waveforms_collection_

private

Definition at line 34 of file RawSiStripDQM.h.

---

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

• Tracking/include/Tracking/dqm/RawSiStripDQM.h
• Tracking/src/Tracking/dqm/RawSiStripDQM.cxx