Public Member Functions
	HcalSimpleDigiAndRecProducer (const std::string &name, framework::Process &process)

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

void	configure (framework::config::Parameters &ps) override
	Callback for the EventProcessor to configure itself from the given set of 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	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.

Private Attributes
std::string	input_coll_name_ {}

std::string	input_pass_name_ {}

std::string	output_coll_name_ {}

double	mev_per_mip_ {}

double	pe_per_mip_ {}

double	attenuation_length_ {}

double	mean_noise_ {}

std::mt19937	rng_ {}

std::unique_ptr< std::normal_distribution< double > >	position_resolution_smear_

std::unique_ptr< ldmx::NoiseGenerator >	noiseGenerator_ {nullptr}

int	readout_threshold_ {2}

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

Definition at line 14 of file HcalSimpleDigiAndRecProducer.h.

Constructor & Destructor Documentation

◆ HcalSimpleDigiAndRecProducer()

hcal::HcalSimpleDigiAndRecProducer::HcalSimpleDigiAndRecProducer	(	const std::string &	name,
		framework::Process &	process )

inline

Definition at line 16 of file HcalSimpleDigiAndRecProducer.h.

18 : framework::Producer{name, process} {}

framework::Producer

Base class for a module which produces a data product.

Definition EventProcessor.h:252

Member Function Documentation

◆ configure()

void hcal::HcalSimpleDigiAndRecProducer::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 className equal to the EventProcessor class name.

For an example, look at MyProcessor.

Parameters

parameters Parameters for configuration.

Reimplemented from framework::EventProcessor.

Definition at line 8 of file HcalSimpleDigiAndRecProducer.cxx.

                                   {
  input_coll_name_ = ps.getParameter<std::string>("input_coll_name");
  input_pass_name_ = ps.getParameter<std::string>("input_pass_name");
  output_coll_name_ = ps.getParameter<std::string>("output_coll_name");
  mev_per_mip_ = ps.getParameter<double>("mev_per_mip");
  pe_per_mip_ = ps.getParameter<double>("pe_per_mip");
  attenuation_length_ = ps.getParameter<double>("attenuation_length");
  readout_threshold_ = ps.getParameter<int>("readout_threshold");
  mean_noise_ = ps.getParameter<double>("mean_noise");
  position_resolution_smear_ =
      std::make_unique<std::normal_distribution<double>>(
          0.0, ps.getParameter<double>("position_resolution"));
}

◆ onNewRun()

void hcal::HcalSimpleDigiAndRecProducer::onNewRun ( const ldmx::RunHeader & runHeader )

overridevirtual

Callback for the EventProcessor to take any necessary action when the run being processed changes.

Parameters

runHeader The RunHeader containing run information.

Reimplemented from framework::EventProcessor.

Definition at line 23 of file HcalSimpleDigiAndRecProducer.cxx.

                                                                {
  noiseGenerator_ = std::make_unique<ldmx::NoiseGenerator>(mean_noise_, false);
  // hard-code this number, create noise hits for non-zero PEs!
  noiseGenerator_->setNoiseThreshold(1);
  const framework::RandomNumberSeedService& rseed =
      getCondition<framework::RandomNumberSeedService>(
          framework::RandomNumberSeedService::CONDITIONS_OBJECT_NAME);
  noiseGenerator_->seedGenerator(
      rseed.getSeed("HcalSimpleDigiAndRecProducer::NoiseGenerator"));
  rng_.seed(rseed.getSeed("HcalSimpleDigiAndRecProducer"));
}

References framework::RandomNumberSeedService::CONDITIONS_OBJECT_NAME, framework::EventProcessor::getCondition(), and framework::RandomNumberSeedService::getSeed().

◆ produce()

void hcal::HcalSimpleDigiAndRecProducer::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 35 of file HcalSimpleDigiAndRecProducer.cxx.

                                                                {
  const auto& hcalGeometry = getCondition<ldmx::HcalGeometry>(
      ldmx::HcalGeometry::CONDITIONS_OBJECT_NAME);
 
  std::vector<ldmx::HcalHit> hcalRecHits;
 
  auto simHits{event.getCollection<ldmx::SimCalorimeterHit>(input_coll_name_,
                                                            input_pass_name_)};
  std::unordered_map<unsigned int, std::vector<const ldmx::SimCalorimeterHit*>>
      hits_by_id{};
  // Important, has to be a reference so that we don't take the address of a
  // variable that goes out of scope!
  for (const auto& hit : simHits) {
    auto id{hit.getID()};
    auto found{hits_by_id.find(id)};
    if (found == hits_by_id.end()) {
      hits_by_id[id] = std::vector<const ldmx::SimCalorimeterHit*>{&hit};
    } else {
      hits_by_id[id].push_back(&hit);
    }
  }
  for (const auto& [barID, simhits_in_bar] : hits_by_id) {
    ldmx::HcalHit& recHit = hcalRecHits.emplace_back();
    double edep{};
    double time{};
    std::vector<double> pos{0, 0, 0};
    for (auto hit : simhits_in_bar) {
      edep += hit->getEdep();
      double edep_hit = hit->getEdep();
      time += hit->getTime() * edep_hit;
      auto hitPos{hit->getPosition()};
      pos[0] += hitPos[0] * edep_hit;
      pos[1] += hitPos[1] * edep_hit;
      pos[2] += hitPos[2] * edep_hit;
    }
    ldmx::HcalID hitID{barID};
 
    // Position smearing
    double mean_pe{(edep / mev_per_mip_) * pe_per_mip_};
    double xpos{pos[0] / edep};
    double ypos{pos[1] / edep};
    double zpos{pos[2] / edep};
    time /= edep;
 
    auto orientation{hcalGeometry.getScintillatorOrientation(barID)};
    double half_total_width{
        hcalGeometry.getHalfTotalWidth(hitID.section(), hitID.layer())};
    double scint_bar_length{hcalGeometry.getScintillatorLength(hitID)};
 
    auto stripCenter{hcalGeometry.getStripCenterPosition(hitID)};
    if (hitID.section() == ldmx::HcalID::HcalSection::BACK) {
      double distance_along_bar =
          (orientation ==
           ldmx::HcalGeometry::ScintillatorOrientation::horizontal)
              ? xpos
              : ypos;
      if (orientation ==
          ldmx::HcalGeometry::ScintillatorOrientation::horizontal) {
        ypos = stripCenter.y();
        xpos += (*position_resolution_smear_)(rng_);
      } else {
        xpos = stripCenter.x();
        ypos += (*position_resolution_smear_)(rng_);
      }
      zpos = stripCenter.z();
      // Attenuation
      mean_pe *= exp(1. / attenuation_length_);
      double mean_pe_close =
          mean_pe * exp(-1. *
                        ((half_total_width - distance_along_bar) /
                         (scint_bar_length * 0.5)) /
                        attenuation_length_);
      double mean_pe_far =
          mean_pe * exp(-1. *
                        ((half_total_width + distance_along_bar) /
                         (scint_bar_length * 0.5)) /
                        attenuation_length_);
      int PE_close{
          std::poisson_distribution<int>(mean_pe_close + mean_noise_)(rng_)};
      int PE_far{
          std::poisson_distribution<int>(mean_pe_far + mean_noise_)(rng_)};
      recHit.setPE(PE_close + PE_far);
      recHit.setMinPE(std::min(PE_close, PE_far));
    } else {
      // Side HCAL, no attenuation business since single ended readout
      int PE{std::poisson_distribution<int>(mean_pe + mean_noise_)(rng_)};
      recHit.setPE(PE);
      recHit.setMinPE(PE);
 
      // Checks orientation of side Hcal bars, sets center positions and add
      // smearing along bar orientation axis
      if (orientation ==
          ldmx::HcalGeometry::ScintillatorOrientation::horizontal) {
        xpos += (*position_resolution_smear_)(rng_);
        ypos = stripCenter.y();
        zpos = stripCenter.z();
      } else if (orientation ==
                 ldmx::HcalGeometry::ScintillatorOrientation::vertical) {
        xpos = stripCenter.x();
        ypos += (*position_resolution_smear_)(rng_);
        zpos = stripCenter.z();
      } else if (orientation ==
                 ldmx::HcalGeometry::ScintillatorOrientation::depth) {
        xpos = stripCenter.x();
        ypos = stripCenter.y();
        zpos += (*position_resolution_smear_)(rng_);
      } else {
        xpos = stripCenter.x();
        ypos = stripCenter.y();
        zpos = stripCenter.z();
        ldmx_log(warn) << "Bar orientation not found. Hit" << hitID.raw()
                       << "positioned at bar center.";
      }
    }
 
    recHit.setID(hitID.raw());
    recHit.setXPos(xpos);
    recHit.setNoise(false);
    recHit.setYPos(ypos);
    recHit.setZPos(zpos);
    recHit.setTime(time);
    recHit.setSection(hitID.section());
    recHit.setStrip(hitID.strip());
    recHit.setLayer(hitID.layer());
    recHit.setEnergy(edep);
    recHit.setOrientation(static_cast<int>(orientation));
  }
  event.add(output_coll_name_, hcalRecHits);
}