Producer making a collection based on some truth info cuts. More...

#include <TruthHitProducer.h>

Public Member Functions
	TruthHitProducer (const std::string &name, framework::Process &process)
	Constructor.

	~TruthHitProducer ()=default
	Destructor.

void	configure (framework::config::Parameters &parameters) override
	Configure the processor using the given user specified parameters.

void	produce (framework::Event &event) override
	Process the event and put new data products into it.

Public Member Functions inherited from framework::Producer
	Producer (const std::string &name, Process &process)
	Class constructor.

virtual void	beforeNewRun (ldmx::RunHeader &header)
	Handle allowing producers to modify run headers before the run begins.

Public Member Functions inherited from framework::EventProcessor
	EventProcessor (const std::string &name, Process &process)
	Class constructor.

virtual	~EventProcessor ()
	Class destructor.

virtual void	onNewRun (const ldmx::RunHeader &runHeader)
	Callback for the EventProcessor to take any necessary action when the run being processed changes.

virtual void	onFileOpen (EventFile &eventFile)
	Callback for the EventProcessor to take any necessary action when a new event input ROOT file is opened.

virtual void	onFileClose (EventFile &eventFile)
	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.

Public Attributes
bool	verbose_ {false}
	Class to set the verbosity level.

std::string	inputCollection_
	Name of the input collection containing the sim hits.

std::string	inputPassName_
	Name of the pass that the input collection is on (empty string means take any pass)

std::string	outputCollection_
	Name of the output collection that will be used to store the selected sim hits.

Additional Inherited Members
Static Public Member Functions inherited from framework::EventProcessor
static void	declare (const std::string &classname, int classtype, EventProcessorMaker *)
	Internal function which is part of the PluginFactory machinery.

Static Public Attributes inherited from framework::Producer
static const int	CLASSTYPE {1}
	Constant used to track EventProcessor types by the PluginFactory.

Protected Member Functions inherited from framework::EventProcessor
void	abortEvent ()
	Abort the event immediately.

Protected Attributes inherited from framework::EventProcessor
HistogramHelper	histograms_
	Interface class for making and filling histograms.

NtupleManager &	ntuple_ {NtupleManager::getInstance()}
	Manager for any ntuples.

logging::logger	theLog_
	The logger for this EventProcessor.

Detailed Description

Producer making a collection based on some truth info cuts.

Like, skimming out only the beam electrons. This producer creates a new hit collection, that can then be fed to digi etc like anything else.

Definition at line 23 of file TruthHitProducer.h.

Constructor & Destructor Documentation

◆ TruthHitProducer()

trigscint::TruthHitProducer::TruthHitProducer	(	const std::string &	name,
		framework::Process &	process
	)

Constructor.

Parameters

name	Name for this instance of the class.
process	The Process class associated with EventProcessor, provided by the framework.

Definition at line 6 of file TruthHitProducer.cxx.

8 : Producer(name, process) {}

framework::Producer::Producer

Producer(const std::string &name, Process &process)

Class constructor.

Definition EventProcessor.cxx:68

Member Function Documentation

◆ configure()

void trigscint::TruthHitProducer::configure ( framework::config::Parameters & parameters )

overridevirtual

Configure the processor using the given user specified parameters.

The user specified parameters that are availabed are defined in the python configuration class. Look at the my_processor.py module of the EventProc python for the python structure.

Parameters

parameters Set of parameters used to configure this processor.

Reimplemented from framework::EventProcessor.

Definition at line 10 of file TruthHitProducer.cxx.

                                                                      {
  inputCollection_ = parameters.getParameter<std::string>("input_collection");
  inputPassName_ = parameters.getParameter<std::string>("input_pass_name");
  outputCollection_ = parameters.getParameter<std::string>("output_collection");
  verbose_ = parameters.getParameter<bool>("verbose");
 
  if (verbose_) {
    ldmx_log(info) << "In TruthHitProducer: configure done!";
    ldmx_log(info) << "Got parameters:  "
                   << "\nInput collection:     " << inputCollection_
                   << "\nInput pass name:     " << inputPassName_
                   << "\nOutput collection:    " << outputCollection_
                   << "\nVerbose: " << verbose_;
  }
}

References framework::config::Parameters::getParameter(), inputCollection_, inputPassName_, outputCollection_, and verbose_.

◆ produce()

void trigscint::TruthHitProducer::produce ( framework::Event & event )

overridevirtual

Process the event and put new data products into it.

Parameters

event The event to process.

Implements framework::Producer.

Definition at line 26 of file TruthHitProducer.cxx.

                                                    {
  // Check if the collection exists.  If not, don't bother processing the event.
  if (!event.exists(inputCollection_)) {
    ldmx_log(error) << "No input collection called " << inputCollection_
                    << " found; skipping!";
    return;
  }
  // looper over sim hits and aggregate energy depositions for each detID
  const auto simHits{event.getCollection<ldmx::SimCalorimeterHit>(
      inputCollection_, inputPassName_)};
  auto particleMap{event.getMap<int, ldmx::SimParticle>("SimParticles")};
 
  std::vector<ldmx::SimCalorimeterHit> truthBeamElectrons;
 
  // TODO: Convert this to using a for_each and lambda
  for (const auto &simHit : simHits) {
    bool keep{false};
    // check if hit is from beam electron and, if so, add to output collection
    for (int i = 0; i < simHit.getNumberOfContribs(); i++) {
      auto contrib = simHit.getContrib(i);
      if (verbose_) {
        ldmx_log(debug) << "contrib " << i << " trackID: " << contrib.trackID
                        << " pdgID: " << contrib.pdgCode
                        << " edep: " << contrib.edep;
        ldmx_log(debug) << "\t particle id: "
                        << particleMap[contrib.trackID].getPdgID()
                        << " particle status: "
                        << particleMap[contrib.trackID].getGenStatus();
      }
      // if the trackID is in the map
      if (particleMap.find(contrib.trackID) != particleMap.end()) {
        // beam electron (PDGID = 11, genStatus == 1)
        if (particleMap[contrib.trackID].getPdgID() == 11 &&
            particleMap[contrib.trackID].getGenStatus() == 1) {
          keep = true;
        }
      }
      if (keep) truthBeamElectrons.push_back(simHit);
    }
  }
  event.add(outputCollection_, truthBeamElectrons);
}

References framework::Event::exists(), inputCollection_, inputPassName_, outputCollection_, and verbose_.

Member Data Documentation

◆ inputCollection_

std::string trigscint::TruthHitProducer::inputCollection_

Name of the input collection containing the sim hits.

Definition at line 60 of file TruthHitProducer.h.

Referenced by configure(), and produce().

◆ inputPassName_

std::string trigscint::TruthHitProducer::inputPassName_

Name of the pass that the input collection is on (empty string means take any pass)

Definition at line 64 of file TruthHitProducer.h.

Referenced by configure(), and produce().

◆ outputCollection_

std::string trigscint::TruthHitProducer::outputCollection_

Name of the output collection that will be used to store the selected sim hits.

Definition at line 68 of file TruthHitProducer.h.

Referenced by configure(), and produce().

◆ verbose_

bool trigscint::TruthHitProducer::verbose_ {false}

Class to set the verbosity level.

Definition at line 57 of file TruthHitProducer.h.

57{false};

Referenced by configure(), and produce().

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

TrigScint/include/TrigScint/TruthHitProducer.h
TrigScint/src/TrigScint/TruthHitProducer.cxx

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TruthHitProducer()

Member Function Documentation

◆ configure()

◆ produce()

Member Data Documentation

◆ inputCollection_

◆ inputPassName_

◆ outputCollection_

◆ verbose_