dqm::GenieTruthDQM Class Reference

Generate histograms/ntuple to extract genie output info. More...

#include <GenieTruthDQM.h>

Public Member Functions
	GenieTruthDQM (const std::string &name, framework::Process &process)
	Constructor.
virtual void	configure (framework::config::Parameters &ps)
	Input python configuration parameters.
virtual void	onProcessStart ()
	Construct histograms/ntuples.
virtual void	onNewRun (const ldmx::RunHeader &runHeader)
	Grab the run number...
virtual void	analyze (const framework::Event &event)
	Fills histograms/ntuples.
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	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	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 Attributes
std::string	hepmc3_coll_name_
	Pass Name for genie objects.
std::string	hepmc3_pass_name_
int	run_number_

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

Generate histograms/ntuple to extract genie output info.

Definition at line 18 of file GenieTruthDQM.h.

Constructor & Destructor Documentation

GenieTruthDQM()

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

inline

Constructor.

Blank Analyzer constructor

Definition at line 25 of file GenieTruthDQM.h.

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

Member Function Documentation

analyze()

void dqm::GenieTruthDQM::analyze ( const framework::Event & event )

virtual

Fills histograms/ntuples.

Implements framework::Analyzer.

Definition at line 72 of file GenieTruthDQM.cxx.

                                                       {
  getHistoDirectory();
 
  ntuple_.clear();
  ntuple_.setVar<int>("run", run_number_);
  ntuple_.setVar<int>("event", event.getEventNumber());
 
  auto hepmc3_col = event.getObject<std::vector<ldmx::HepMC3GenEvent> >(
      hepmc3_coll_name_, hepmc3_pass_name_);
 
  if (hepmc3_col.size() < 1) {
    ntuple_.fill();
    return;
  }
 
  auto const& hepmc3_ev = hepmc3_col.at(0).getHepMCGenEvent();
 
  // set interaction/scattering codes
  auto interaction_type_ptr = hepmc3_ev.attribute<HepMC3::IntAttribute>(
      "GENIE.Interaction.InteractionType");
  if (interaction_type_ptr)
    ntuple_.setVar<int>("interaction_type", interaction_type_ptr->value());
  auto scattering_type_ptr = hepmc3_ev.attribute<HepMC3::IntAttribute>(
      "GENIE.Interaction.ScatteringType");
  if (scattering_type_ptr)
    ntuple_.setVar<int>("scattering_type", scattering_type_ptr->value());
  auto rescatter_code_ptr =
      hepmc3_ev.attribute<HepMC3::IntAttribute>("GENIE.RescatterCode");
  if (rescatter_code_ptr)
    ntuple_.setVar<int>("rescatter_code", rescatter_code_ptr->value());
 
  // get kinematic vars
  auto kvar_labels_ptr = hepmc3_ev.attribute<HepMC3::VectorIntAttribute>(
      "GENIE.Interaction.KineVarLabels");
  auto kvar_values_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.KineVarValues");
  if (kvar_labels_ptr && kvar_values_ptr) {
    auto kvar_labels = kvar_labels_ptr->value();
    auto kvar_values = kvar_values_ptr->value();
    for (size_t i = 0; i < kvar_labels.size(); ++i) {
      if (kvar_labels[i] == genie::EKineVar::kKVSelx)
        ntuple_.setVar<double>("x_bj", kvar_values[i]);
      else if (kvar_labels[i] == genie::EKineVar::kKVSely)
        ntuple_.setVar<double>("y_inel", kvar_values[i]);
      else if (kvar_labels[i] == genie::EKineVar::kKVSelQ2)
        ntuple_.setVar<double>("Q2", kvar_values[i]);
      else if (kvar_labels[i] == genie::EKineVar::kKVSelW)
        ntuple_.setVar<double>("W", kvar_values[i]);
    }
  }
 
  // electron info
  auto lep_pdg_ptr =
      hepmc3_ev.attribute<HepMC3::IntAttribute>("GENIE.Interaction.ProbePDG");
  if (lep_pdg_ptr) ntuple_.setVar<int>("lep_pdg", lep_pdg_ptr->value());
 
  auto lep_i_4vec_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.ProbeP4");
  if (lep_i_4vec_ptr) {
    auto lep_i_4vec = lep_i_4vec_ptr->value();
    ntuple_.setVar<double>("lep_i_px", lep_i_4vec[0]);
    ntuple_.setVar<double>("lep_i_py", lep_i_4vec[1]);
    ntuple_.setVar<double>("lep_i_pz", lep_i_4vec[2]);
    ntuple_.setVar<double>("lep_i_e", lep_i_4vec[3]);
  }
  auto lep_f_4vec_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.FSLeptonP4");
  if (lep_f_4vec_ptr) {
    auto lep_f_4vec = lep_f_4vec_ptr->value();
    ntuple_.setVar<double>("lep_f_px", lep_f_4vec[0]);
    ntuple_.setVar<double>("lep_f_py", lep_f_4vec[1]);
    ntuple_.setVar<double>("lep_f_pz", lep_f_4vec[2]);
    ntuple_.setVar<double>("lep_f_e", lep_f_4vec[3]);
  }
 
  // target info
  auto tgt_pdg_ptr =
      hepmc3_ev.attribute<HepMC3::IntAttribute>("GENIE.Interaction.TargetPDG");
  if (tgt_pdg_ptr) ntuple_.setVar<int>("tgt_pdg", tgt_pdg_ptr->value());
 
  auto tgt_4vec_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.TargetP4");
  if (tgt_4vec_ptr) {
    auto tgt_4vec = tgt_4vec_ptr->value();
    ntuple_.setVar<double>("tgt_px", tgt_4vec[0]);
    ntuple_.setVar<double>("tgt_py", tgt_4vec[1]);
    ntuple_.setVar<double>("tgt_pz", tgt_4vec[2]);
    ntuple_.setVar<double>("tgt_e", tgt_4vec[3]);
  }
 
  // hit nucleon info
  auto hnuc_pdg_ptr = hepmc3_ev.attribute<HepMC3::IntAttribute>(
      "GENIE.Interaction.HitNucleonPDG");
  if (hnuc_pdg_ptr) ntuple_.setVar<int>("hnuc_pdg", hnuc_pdg_ptr->value());
 
  auto hitnuc_4vec_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.HitNucleonP4");
  if (hitnuc_4vec_ptr) {
    auto hitnuc_4vec = hitnuc_4vec_ptr->value();
    ntuple_.setVar<double>("hnuc_px", hitnuc_4vec[0]);
    ntuple_.setVar<double>("hnuc_py", hitnuc_4vec[1]);
    ntuple_.setVar<double>("hnuc_pz", hitnuc_4vec[2]);
    ntuple_.setVar<double>("hnuc_e", hitnuc_4vec[3]);
  }
  // hit quark info
  // note: it's only there for some interaction types!
  auto hqrkpdg_ptr = hepmc3_ev.attribute<HepMC3::IntAttribute>(
      "GENIE.Interaction.HitQuarkPDG");
  if (hqrkpdg_ptr) ntuple_.setVar<int>("hqrk_pdg", hqrkpdg_ptr->value());
  auto hqrksea_ptr = hepmc3_ev.attribute<HepMC3::IntAttribute>(
      "GENIE.Interaction.HitSeaQuark");
  if (hqrksea_ptr) ntuple_.setVar<int>("hqrk_sea", hqrksea_ptr->value());
 
  auto hadsys_4vec_ptr = hepmc3_ev.attribute<HepMC3::VectorDoubleAttribute>(
      "GENIE.Interaction.HadSystP4");
  if (hadsys_4vec_ptr) {
    auto hadsys_4vec = hadsys_4vec_ptr->value();
    ntuple_.setVar<double>("hadsys_px", hadsys_4vec[0]);
    ntuple_.setVar<double>("hadsys_py", hadsys_4vec[1]);
    ntuple_.setVar<double>("hadsys_pz", hadsys_4vec[2]);
    ntuple_.setVar<double>("hadsys_e", hadsys_4vec[3]);
  }
  ntuple_.fill();
 
  return;
}

References framework::NtupleManager::clear(), framework::EventProcessor::getHistoDirectory(), hepmc3_coll_name_, framework::EventProcessor::ntuple_, and framework::NtupleManager::setVar().

configure()

void dqm::GenieTruthDQM::configure ( framework::config::Parameters & ps )

virtual

Input python configuration parameters.

Reimplemented from framework::EventProcessor.

Definition at line 12 of file GenieTruthDQM.cxx.

                                                           {
  hepmc3_coll_name_ = ps.get<std::string>("hepmc3_coll_name");
  hepmc3_pass_name_ = ps.get<std::string>("hepmc3_pass_name");
  return;
}

References framework::config::Parameters::get(), and hepmc3_coll_name_.

onNewRun()

void dqm::GenieTruthDQM::onNewRun ( const ldmx::RunHeader & runHeader )

virtual

Grab the run number...

Reimplemented from framework::EventProcessor.

Definition at line 68 of file GenieTruthDQM.cxx.

                                                           {
  run_number_ = runHeader.getRunNumber();
}

References ldmx::RunHeader::getRunNumber().

onProcessStart()

void dqm::GenieTruthDQM::onProcessStart ( )

virtual

Construct histograms/ntuples.

Reimplemented from framework::EventProcessor.

Definition at line 18 of file GenieTruthDQM.cxx.

                                   {
  getHistoDirectory();
 
  ntuple_.create("genie_truth");
 
  // event info
  ntuple_.addVar<int>("genie_truth", "run");
  ntuple_.addVar<int>("genie_truth", "event");
 
  ntuple_.addVar<int>("genie_truth", "interaction_type");
  ntuple_.addVar<int>("genie_truth", "scattering_type");
  ntuple_.addVar<int>("genie_truth", "rescatter_code");
 
  ntuple_.addVar<double>("genie_truth", "x_bj");
  ntuple_.addVar<double>("genie_truth", "y_inel");
  ntuple_.addVar<double>("genie_truth", "Q2");
  ntuple_.addVar<double>("genie_truth", "W");
 
  ntuple_.addVar<int>("genie_truth", "lep_pdg");
 
  ntuple_.addVar<double>("genie_truth", "lep_i_px");
  ntuple_.addVar<double>("genie_truth", "lep_i_py");
  ntuple_.addVar<double>("genie_truth", "lep_i_pz");
  ntuple_.addVar<double>("genie_truth", "lep_i_e");
 
  ntuple_.addVar<double>("genie_truth", "lep_f_px");
  ntuple_.addVar<double>("genie_truth", "lep_f_py");
  ntuple_.addVar<double>("genie_truth", "lep_f_pz");
  ntuple_.addVar<double>("genie_truth", "lep_f_e");
 
  ntuple_.addVar<int>("genie_truth", "tgt_pdg");
  ntuple_.addVar<double>("genie_truth", "tgt_px");
  ntuple_.addVar<double>("genie_truth", "tgt_py");
  ntuple_.addVar<double>("genie_truth", "tgt_pz");
  ntuple_.addVar<double>("genie_truth", "tgt_e");
 
  ntuple_.addVar<int>("genie_truth", "hnuc_pdg");
  ntuple_.addVar<double>("genie_truth", "hnuc_px");
  ntuple_.addVar<double>("genie_truth", "hnuc_py");
  ntuple_.addVar<double>("genie_truth", "hnuc_pz");
  ntuple_.addVar<double>("genie_truth", "hnuc_e");
 
  ntuple_.addVar<int>("genie_truth", "hqrk_pdg");
  ntuple_.addVar<int>("genie_truth", "hqrk_sea");
  ntuple_.addVar<double>("genie_truth", "hadsys_px");
  ntuple_.addVar<double>("genie_truth", "hadsys_py");
  ntuple_.addVar<double>("genie_truth", "hadsys_pz");
  ntuple_.addVar<double>("genie_truth", "hadsys_e");
}

References framework::NtupleManager::addVar(), framework::NtupleManager::create(), framework::EventProcessor::getHistoDirectory(), and framework::EventProcessor::ntuple_.

Member Data Documentation

hepmc3_coll_name_

std::string dqm::GenieTruthDQM::hepmc3_coll_name_

private

Pass Name for genie objects.

Definition at line 50 of file GenieTruthDQM.h.

Referenced by analyze(), and configure().

hepmc3_pass_name_

std::string dqm::GenieTruthDQM::hepmc3_pass_name_

private

Definition at line 51 of file GenieTruthDQM.h.

run_number_

int dqm::GenieTruthDQM::run_number_

private

Definition at line 53 of file GenieTruthDQM.h.

---

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

• DQM/include/DQM/GenieTruthDQM.h
• DQM/src/DQM/GenieTruthDQM.cxx