recon::OverlayProducer Class Reference

Class to overlay in-time pile-up events from an overlay file. More...

#include <OverlayProducer.h>

Public Member Functions
	OverlayProducer (const std::string &name, framework::Process &process)
void	configure (framework::config::Parameters &parameters) override
	Configure the processor with input parameters from the python cofig.
void	onNewRun (const ldmx::RunHeader &) override
	At the start of the run, the pileup overlay file is set up, and the starting event number is chosen, using the RNSS.
void	produce (framework::Event &event) override
	Based on the list of collections to overlay, and the desired number of events, loop through all relevant collections and copy the sim event (once), and then add the corresponding collection from the pileup overlay file.
void	onProcessStart () override
	At the start of processing, the pileup overlay file is set up.
Public Member Functions inherited from framework::Producer
	Producer (const std::string &name, Process &process)
	Class constructor.
virtual void	process (Event &event) final
	Processing an event for a Producer is calling produce.
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	beforeNewRun (ldmx::RunHeader &run_header)
	Callback for Producers to add parameters to the run header before conditions are initialized.
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
framework::config::Parameters	params_
	The parameters used to configure this producer.
std::string	overlay_filename_
	Pileup overlay events input file name.
std::unique_ptr< framework::EventFile >	overlay_file_
	Pileup overlay events input file.
framework::Event	overlay_event_
	The overlay ldmx event bus.
std::vector< std::string >	calo_collections_
	List of SimCalorimeterHit collection(s) to loop over and add hits from, combining sim and pileup.
std::vector< std::string >	tracker_collections_
	List of SimTrackerHit collection(s) to loop over and add hits from, combining sim and pileup.
std::vector< std::string >	contrib_collections_
	List of SimCalorimeterHit collections which keep track of hit contribs.
std::string	overlay_passname_
	Pileup overlay events input pass name.
std::string	sim_passname_
	To use for finding the sim event bus passengers, mostly a disambiguation.
std::string	out_coll_postfix_
	Postfix to add to the collection name of the overlayed collections.
bool	do_poisson_in_time_ {false}
	Let the total number of in-time events be poisson distributed, or fix at the chosen value, poissonMu_.
bool	do_poisson_out_of_time_ {false}
	Let the total number of out-of-time events be poisson distributed, or fix at the chosen value, poissonMu_.
double	poisson_mu_ {0.}
	(average) total number of events
std::unique_ptr< TRandom2 >	rndm_
	Random number generator for number of overlaid events.
std::unique_ptr< TRandom2 >	rndm_time_
	Random number generator for pileup event time offset.
double	time_sigma_ {0.}
	Width of pileup bunch spread in time (in [ns]), specified as a sigma of a Gaussian distribution.
double	time_mean_ {0.}
	Average position in time (in [ns]) of pileup bunches, relative to the sim event.
double	bunch_spacing_ {0.}
	Spacing in time (in [ns]) between electron bunches.
int	n_earlier_ {0}
	Number of bunches before the sim event to pull pileup events from.
int	n_later_ {0}
	Number of bunches after the sim event to pull pileup events from.
int	start_event_min_ {1}
	Minimum event number to start overlaying from.
int	start_event_max_ {10000}
	Maximum event number to start overlaying from.
int	overlay_incident_id_ {-1000}
	For Ecal, overlay hits should be added as contribs.
int	overlay_track_id_ {-1000}
int	overlay_pdg_code_ {0}

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

Class to overlay in-time pile-up events from an overlay file.

Definition at line 29 of file OverlayProducer.h.

Constructor & Destructor Documentation

OverlayProducer()

recon::OverlayProducer::OverlayProducer ( const std::string & name, framework::Process & process )

inline

Definition at line 31 of file OverlayProducer.h.

      : framework::Producer(name, process), overlay_event_{"overlay"} {}

Member Function Documentation

configure()

void recon::OverlayProducer::configure ( framework::config::Parameters & parameters )

overridevirtual

Configure the processor with input parameters from the python cofig.

Print the parameters actually set. Helpful in case of typos.

Reimplemented from framework::EventProcessor.

Definition at line 5 of file OverlayProducer.cxx.

                                                                     {
  params_ = parameters;
 
  // name of file containing events to be overlaid, and a list of collections to
  // overlay
  overlay_filename_ = parameters.get<std::string>("overlay_filename");
  calo_collections_ =
      parameters.get<std::vector<std::string>>("calo_collections");
  tracker_collections_ =
      parameters.get<std::vector<std::string>>("tracker_collections");
  contrib_collections_ =
      parameters.get<std::vector<std::string>>("contrib_collections");
  sim_passname_ = parameters.get<std::string>("sim_passname");
  overlay_passname_ = parameters.get<std::string>("overlay_passname");
  out_coll_postfix_ = parameters.get<std::string>("out_coll_postfix");
  // overlay specifics:
  poisson_mu_ = parameters.get<double>("poisson_mu");
  do_poisson_in_time_ = parameters.get<bool>("do_poisson_in_time");
  do_poisson_out_of_time_ = parameters.get<bool>("do_poisson_out_of_time");
  time_sigma_ = parameters.get<double>("time_sigma");
  time_mean_ = parameters.get<double>("time_mean");
  n_earlier_ = parameters.get<int>("n_earlier");
  n_later_ = parameters.get<int>("n_later");
  bunch_spacing_ = parameters.get<double>("bunch_spacing");
  start_event_min_ = parameters.get<int>("start_event_min");
  start_event_max_ = parameters.get<int>("start_event_max");
 
  ldmx_log(debug) << "Got parameters \n \t overlayFileName = "
                  << overlay_filename_
                  << "\n\t sim pass name = " << sim_passname_
                  << "\n\t overlay pass name = " << overlay_passname_;
  ldmx_log(debug) << "\n\t overlayCaloHitCollections = ";
  for (const auto &coll : calo_collections_) {
    ldmx_log(debug) << coll << "; ";
  }
 
  ldmx_log(debug) << "\n\t overlayTrackerHitCollections = ";
  for (const std::string &coll : tracker_collections_) {
    ldmx_log(debug) << coll << "; ";
  }
 
  ldmx_log(trace) << "\n\t numberOverlaidInteractions = " << poisson_mu_
                  << "\n\t nEarlierBunchesToSample = " << n_earlier_
                  << "\n\t nLaterBunchesToSample = " << n_later_
                  << "\n\t bunchSpacing = " << bunch_spacing_
                  << "\n\t doPoissonIntime = " << do_poisson_in_time_
                  << "\n\t doPoissonOutoftime = " << do_poisson_out_of_time_
                  << "\n\t timeSpread = " << time_sigma_
                  << "\n\t timeMean = " << time_mean_
                  << "\n\t startEventMin = " << start_event_min_
                  << "\n\t startEventMax = " << start_event_max_;
 
  // given how the contrib collection specification is setup, it's possible for
  // spelling errors to cause doom here, so we want to flag those
  for (auto coll_name : contrib_collections_) {
    if (std::find(calo_collections_.begin(), calo_collections_.end(),
                  coll_name) == calo_collections_.end()) {
      EXCEPTION_RAISE("CollNameMismatch",
                      "The contrib-using collection " + coll_name +
                          " does not match any SimCalorimeterHit collection in "
                          "'calo_collections_'! Please check your spelling");
    }
  }
 
  return;
}  // end configure()

References bunch_spacing_, calo_collections_, contrib_collections_, do_poisson_in_time_, do_poisson_out_of_time_, framework::config::Parameters::get(), n_earlier_, n_later_, out_coll_postfix_, overlay_filename_, overlay_passname_, params_, poisson_mu_, sim_passname_, start_event_max_, start_event_min_, time_mean_, time_sigma_, and tracker_collections_.

onNewRun()

void recon::OverlayProducer::onNewRun ( const ldmx::RunHeader & )

overridevirtual

At the start of the run, the pileup overlay file is set up, and the starting event number is chosen, using the RNSS.

set up random seeds

Reimplemented from framework::EventProcessor.

Definition at line 73 of file OverlayProducer.cxx.

                                                    {
  const auto &rnss = getCondition<framework::RandomNumberSeedService>(
      framework::RandomNumberSeedService::CONDITIONS_OBJECT_NAME);
  rndm_ = std::make_unique<TRandom2>(rnss.getSeed("OverlayProducer::rndm"));
  rndm_time_ =
      std::make_unique<TRandom2>(rnss.getSeed("OverlayProducer::rndmTime"));
 
  // Pick a random event from the Pileup file
  int start_event = rndm_->Uniform(start_event_min_, start_event_max_);
  // EventFile::skipToEvent handles actual number of events in file
  int ev_number = overlay_file_->skipToEvent(start_event);
  if (ev_number < 0) {
    EXCEPTION_RAISE("BadRead", "Couldn't read to starting offset.");
  }
  overlay_event_.getEventHeader().setEventNumber(ev_number);
  ldmx_log(info) << "Starting overlay process with pileup event number "
                 << ev_number << " (random event number picked was "
                 << start_event << ").";
}  // end onNewRun()

References framework::RandomNumberSeedService::CONDITIONS_OBJECT_NAME, framework::EventProcessor::getCondition(), framework::Event::getEventHeader(), overlay_event_, overlay_file_, rndm_, rndm_time_, ldmx::EventHeader::setEventNumber(), start_event_max_, and start_event_min_.

onProcessStart()

void recon::OverlayProducer::onProcessStart ( )

overridevirtual

At the start of processing, the pileup overlay file is set up.

Reimplemented from framework::EventProcessor.

Definition at line 391 of file OverlayProducer.cxx.

                                     {
  // replace by this line once the corresponding tweak to EventFile is ready:
  overlay_file_ =
      std::make_unique<framework::EventFile>(params_, overlay_filename_, true);
  overlay_file_->setupEvent(&overlay_event_);
  // we update the iterator at the end of each event. so do this once here to
  // grab the first event in the processor
  // ldmx_log(trace) << "Used input file: "
  // << overlay_file_->getFileName() << " Got event info: " <<
  // overlay_file_->getEvent()->Print();
 
  return;
}  // end onProcessStart

References overlay_event_, overlay_file_, overlay_filename_, and params_.

produce()

void recon::OverlayProducer::produce ( framework::Event & event )

overridevirtual

Based on the list of collections to overlay, and the desired number of events, loop through all relevant collections and copy the sim event (once), and then add the corresponding collection from the pileup overlay file.

The event loop is the outer loop, and the inner loop is over the list of collections.

The collections have to be specified separately as a list of SimCalorimeterHit collections and a list of SimTrackerHit collections.

The collection name is parsed for "Ecal" to be flagged as a collection which needs overlay hits to be added as contribs. This is currently hardwired.

The resulting collections inherit the input collection name, with an appended string "Overlay". This name is also currently hardwired.

Go to next overlay event This overlay file has been configured to loop back to the beginning of the TTree when it reaches the end. This means nextEvent() will only return false if an error is occurred or if the overlay file is mis-configured.

Implements framework::Producer.

Definition at line 94 of file OverlayProducer.cxx.

                                                   {
  // using nextEvent to loop, we need to loop over overlay events and in an
  // inner loop, loop over collections, and store them. after all pileup events
  // have been added, the vector of collections is iterated over and added to
  // the event bus.
  std::map<std::string, std::vector<ldmx::SimCalorimeterHit>>
      calo_collection_map;
  std::map<std::string, std::vector<ldmx::SimTrackerHit>>
      tracker_collection_map;
  std::map<std::string, std::map<int, ldmx::SimCalorimeterHit>> hit_map;
 
  // start by copying over all the collections from the sim event
 
  /* ----------- first do the SimCalorimeterHits ----------- */
 
  // get the calo hits_ collections that we want to overlay, by looping over
  // the list of collections passed to the producer : calo_collections_
  for (const auto &coll_name : calo_collections_) {
    // search for the collection name in the list of collections that
    // need contribs to be added, contrib_collections_
    bool needs_contribs_added{std::find(contrib_collections_.begin(),
                                        contrib_collections_.end(),
                                        coll_name) != contrib_collections_.end()
                                  ? true
                                  : false};
 
    // start out by just copying the sim hits, unaltered.
    auto simhits_calo =
        event.getCollection<ldmx::SimCalorimeterHit>(coll_name, sim_passname_);
    // but don't copy contrib using hits immediately: for them, wait until
    // overlay contribs have been added. then add everything through the hit_map
    if (!needs_contribs_added) {
      calo_collection_map[coll_name + out_coll_postfix_] = simhits_calo;
    }
 
    ldmx_log(debug) << "in loop: start of collection " << coll_name
                    << "in loop: printing current sim event: ";
    ldmx_log(debug) << "in loop: size of sim hits_ vector " << coll_name
                    << " is " << simhits_calo.size();
 
    // we don't need to touch the hard process sim hits_, really... but we
    // might need the simhits in the hit map.
    if (needs_contribs_added) {
      ldmx_log(trace) << "Collection " << coll_name << " needs contribs added";
      for (const ldmx::SimCalorimeterHit &simhit : simhits_calo) {
        ldmx_log(trace) << simhit;
        // this copies the hit, its ID and its coordinates directly
        hit_map[coll_name + out_coll_postfix_][simhit.getID()] = simhit;
 
      }  // over calo simhit collection
    }  // if needContribs
 
  }  // over calo collections for sim event
 
  /* ----------- then do the same with SimTrackerHits! ----------- */
 
  // get the SimTrackerHit collections that we want to overlay, by looping
  // over the list of collections passed to the producer : tracker_collections_
  for (const auto &coll_name : tracker_collections_) {
    auto simhits_tracker =
        event.getCollection<ldmx::SimTrackerHit>(coll_name, sim_passname_);
    tracker_collection_map[coll_name + out_coll_postfix_] = simhits_tracker;
 
    // the rest is printouts for debugging
    ldmx_log(debug) << "in loop: size of sim hits_ vector " << coll_name
                    << " is " << simhits_tracker.size();
 
    ldmx_log(debug) << "in loop: start of collection " << coll_name
                    << "in loop: printing current sim event: ";
 
    for (const ldmx::SimTrackerHit &simhit : simhits_tracker) {
      ldmx_log(trace) << simhit;
    }
  }  // over tracker collections for sim event
 
  /* ----------- now do the pileup overlay ----------- */
 
  // we could shift these by a random number, effectively placing the
  // sim event at random positions in the interval, preserving the
  // overall interval length
  // int simBunch=  static_cast<int>(rndm_time_->Uniform(
  //                   -(n_earlier_+1) , n_later_+1));  // +1 to get
  // inclusive interval
  int start_bunch = -n_earlier_;
  int end_bunch = n_later_;
 
  // TODO -- figure out if we should also randomly shift the time of the sim
  // event (likely only needed if time bias gets picked up by BDT or ML by way
  // of pulse behaviour)
  for (int bunch_offset{start_bunch}; bunch_offset <= end_bunch;
       bunch_offset++) {
    // sample a poisson distribution, or use mu as fixed number of overlay
    // events
    int n_events_overlay = do_poisson_out_of_time_
                               ? rndm_->Poisson(poisson_mu_)
                               : static_cast<int>(poisson_mu_);
 
    // special case: in-time pileup at bunch 0
    if (bunch_offset == 0) {
      if (!do_poisson_in_time_) {
        // fix it to the average
        n_events_overlay = static_cast<int>(poisson_mu_);
      } else if (do_poisson_in_time_ && !do_poisson_out_of_time_) {
        // then we haven't set this yet
        n_events_overlay = rndm_->Poisson(poisson_mu_);
      }
 
      // paticularly useful in the poisson fluctuated case
      event.getEventHeader().setIntParameter("in_time_pu", n_events_overlay);
 
      // the total number of events is nPU + 1 (it includes the sim event)
      // in any case, subtract the sim event from nOverlay
      n_events_overlay -= 1;
 
    }  // end if bunch_offset == 0
 
    float bunchtime_offset = bunch_spacing_ * bunch_offset;
    ldmx_log(debug) << "Will overlay " << n_events_overlay
                    << " events on the simulated one";
 
    for (int i_ev = 0; i_ev < n_events_overlay; i_ev++) {
      if (!overlay_file_->nextEvent()) {
        ldmx_log(error) << "Couldn't read next overlay event!";
        return;
      }
 
      // a pileup event wide time offset to be applied to all its hits_.
      float time_offset = rndm_time_->Gaus(time_mean_, time_sigma_);
      time_offset += bunchtime_offset;
 
      ldmx_log(trace) << "in overlay loop: overlaying event " << "which is "
                      << i_ev + 1 << " out of " << n_events_overlay
                      << "\n\thit time offset is " << time_offset << " ns"
                      << "\n\tbunch position offset is " << bunch_offset
                      << ", leading to a total time offset of "
                      << bunchtime_offset << " ns";
 
      /* ----------- first do the SimCalorimeterHits overlay ----------- */
 
      // again get the calo hits_ collections that we want to overlay
      for (uint i_coll = 0; i_coll < calo_collections_.size(); i_coll++) {
        // search for the collection name in the list of collections that
        // need contribs to be added, contrib_collections_
        bool needs_contribs_added{
            std::find(contrib_collections_.begin(), contrib_collections_.end(),
                      calo_collections_[i_coll]) != contrib_collections_.end()
                ? true
                : false};
 
        std::vector<ldmx::SimCalorimeterHit> overlay_hits =
            overlay_event_.getCollection<ldmx::SimCalorimeterHit>(
                calo_collections_[i_coll], overlay_passname_);
 
        ldmx_log(debug) << "in loop: size of overlay hits_ vector is "
                        << overlay_hits.size();
 
        std::string out_coll_name =
            calo_collections_[i_coll] + out_coll_postfix_;
 
        ldmx_log(trace) << "in loop: printing overlay event: ";
 
        for (ldmx::SimCalorimeterHit &overlay_hit : overlay_hits) {
          ldmx_log(trace) << overlay_hit;
 
          const float overlay_time = overlay_hit.getTime() + time_offset;
          overlay_hit.setTime(overlay_time);
 
          if (needs_contribs_added) {  // special treatment for (for now only)
                                       // ecal
            auto &this_coll_hit_map{
                hit_map[calo_collections_[i_coll] + out_coll_postfix_]};
            int overlay_hit_id = overlay_hit.getID();
            if (this_coll_hit_map.find(overlay_hit_id) ==
                this_coll_hit_map
                    .end()) {  // there wasn't already a simhit in this id
              this_coll_hit_map[overlay_hit_id] = ldmx::SimCalorimeterHit();
              this_coll_hit_map[overlay_hit_id].setID(overlay_hit_id);
              std::vector<float> hit_pos = overlay_hit.getPosition();
              this_coll_hit_map[overlay_hit_id].setPosition(
                  hit_pos[0], hit_pos[1], hit_pos[2]);
            }
            // add the overlay hit (as a) contrib
            // incidentID = -1000, trackID = -1000, pdgCode = 0  <-- these are
            // set in the header for now but could be parameters
            this_coll_hit_map[overlay_hit_id].addContrib(
                overlay_incident_id_, overlay_track_id_, overlay_pdg_code_,
                overlay_hit.getEdep(), overlay_time);
          }  // if add overlay as contribs
          else {
            calo_collection_map[out_coll_name].push_back(overlay_hit);
 
            ldmx_log(trace) << "Adding non-Ecal overlay hit to outhit vector "
                            << out_coll_name;
          }  // end else !needs_contribs_added
        }  // over overlay calo simhit collection
 
        if (!needs_contribs_added)
          ldmx_log(debug) << "Nhits in overlay collection " << out_coll_name
                          << ": " << calo_collection_map[out_coll_name].size();
 
      }  // over calo_collections_
 
      /* ----------- now do simtracker hits_ overlay ----------- */
 
      // get the SimTrackerHit collections that we want to overlay
      for (const auto &coll : tracker_collections_) {
        auto overlay_tracker_hits{
            overlay_event_.getCollection<ldmx::SimTrackerHit>(
                coll, overlay_passname_)};
 
        ldmx_log(debug) << "in loop: size of overlay hits_ vector is "
                        << overlay_tracker_hits.size();
 
        std::string out_coll_name_tracker{coll + out_coll_postfix_};
 
        ldmx_log(trace) << "in loop: printing overlay event: ";
 
        for (auto &overlay_hit : overlay_tracker_hits) {
          auto overlay_time{overlay_hit.getTime() + time_offset};
          overlay_hit.setTime(overlay_time);
          auto overlay_track_id{overlay_hit.getTrackID()};
          overlay_hit.setTrackID(overlay_track_id);
          tracker_collection_map[out_coll_name_tracker].push_back(overlay_hit);
 
          ldmx_log(trace) << overlay_hit;
          ldmx_log(trace) << "Adding tracker overlay hit to outhit vector "
                          << out_coll_name_tracker;
        }  // over overlay tracker simhit collection
 
        ldmx_log(debug) << "Nhits in overlay collection "
                        << out_coll_name_tracker << ": "
                        << tracker_collection_map[out_coll_name_tracker].size();
 
      }  // over tracker_collections_
 
    }  // over overlay events
  }  // over bunches
 
  // after all events are done, the contrib-using collections' hit_maps are
  // final and can be written to the event output
  for (uint i_coll = 0; i_coll < contrib_collections_.size(); i_coll++) {
    // for each SimCalorimeterHit collection that uses contribs, add overlaid
    // hits_ as contribs/from hit_map rather than as copied simhits
    ldmx_log(trace) << "Hits in hit_map after overlay of "
                    << contrib_collections_[i_coll] << "Overlay :";
 
    for (auto &map_hit :
         hit_map[contrib_collections_[i_coll] + out_coll_postfix_]) {
      ldmx_log(trace) << map_hit.second;
 
      if (calo_collection_map.find(contrib_collections_[i_coll] +
                                   out_coll_postfix_) ==
          calo_collection_map.end()) {
        ldmx_log(debug) << "Adding first hit from hit map as first outhit "
                           "vector to calo_collection_map";
        calo_collection_map[contrib_collections_[i_coll] + out_coll_postfix_] =
            {map_hit.second};
      } else {
        calo_collection_map[contrib_collections_[i_coll] + out_coll_postfix_]
            .push_back(map_hit.second);
      }
    }  // over hit_map
  }  // second loop over contrib using collections, to collect hits_ from
     // hit_map
 
  // done collecting hits_.
 
  // now we can write the calo collections to the event bus
  for (auto &[name, coll] : calo_collection_map) {
    ldmx_log(debug) << "Writing " << name << " to event bus.";
 
    ldmx_log(trace) << "List of hits_ added: ";
    for (auto &hit : coll) {
      ldmx_log(trace) << hit;
    }
    event.add(name, coll);
  }
 
  // and now for the tracker hits_
  for (auto &[name, coll] : tracker_collection_map) {
    ldmx_log(debug) << "Writing " << name << " to event bus.";
    ldmx_log(trace) << "List of hits_ added: ";
    for (auto &hit : coll) {
      ldmx_log(trace) << hit;
    }
    event.add(name, coll);
  }
 
  return;
}  // end produce()

References bunch_spacing_, calo_collections_, contrib_collections_, do_poisson_in_time_, do_poisson_out_of_time_, framework::Event::getCollection(), n_earlier_, n_later_, out_coll_postfix_, overlay_event_, overlay_file_, overlay_incident_id_, overlay_passname_, poisson_mu_, rndm_, rndm_time_, sim_passname_, time_mean_, time_sigma_, and tracker_collections_.

Member Data Documentation

bunch_spacing_

double recon::OverlayProducer::bunch_spacing_ {0.}

private

Spacing in time (in [ns]) between electron bunches.

Definition at line 171 of file OverlayProducer.h.

{0.};

Referenced by configure(), and produce().

calo_collections_

std::vector<std::string> recon::OverlayProducer::calo_collections_

private

List of SimCalorimeterHit collection(s) to loop over and add hits from, combining sim and pileup.

Definition at line 95 of file OverlayProducer.h.

Referenced by configure(), and produce().

contrib_collections_

std::vector<std::string> recon::OverlayProducer::contrib_collections_

private

List of SimCalorimeterHit collections which keep track of hit contribs.

Definition at line 106 of file OverlayProducer.h.

Referenced by configure(), and produce().

do_poisson_in_time_

bool recon::OverlayProducer::do_poisson_in_time_ {false}

private

Let the total number of in-time events be poisson distributed, or fix at the chosen value, poissonMu_.

Definition at line 128 of file OverlayProducer.h.

{false};

Referenced by configure(), and produce().

do_poisson_out_of_time_

bool recon::OverlayProducer::do_poisson_out_of_time_ {false}

private

Let the total number of out-of-time events be poisson distributed, or fix at the chosen value, poissonMu_.

Definition at line 134 of file OverlayProducer.h.

{false};

Referenced by configure(), and produce().

n_earlier_

int recon::OverlayProducer::n_earlier_ {0}

private

Number of bunches before the sim event to pull pileup events from.

Defaults to 0 --> all events occur in the same bunch as the sim event.

Definition at line 178 of file OverlayProducer.h.

{0};

Referenced by configure(), and produce().

n_later_

int recon::OverlayProducer::n_later_ {0}

private

Number of bunches after the sim event to pull pileup events from.

Defaults to 0 --> all events occur in the same bunch as the sim event.

Definition at line 185 of file OverlayProducer.h.

{0};

Referenced by configure(), and produce().

out_coll_postfix_

std::string recon::OverlayProducer::out_coll_postfix_

private

Postfix to add to the collection name of the overlayed collections.

This is currently set to "Overlay".

Definition at line 122 of file OverlayProducer.h.

Referenced by configure(), and produce().

overlay_event_

framework::Event recon::OverlayProducer::overlay_event_

private

The overlay ldmx event bus.

Definition at line 89 of file OverlayProducer.h.

Referenced by onNewRun(), onProcessStart(), and produce().

overlay_file_

std::unique_ptr<framework::EventFile> recon::OverlayProducer::overlay_file_

private

Pileup overlay events input file.

Definition at line 84 of file OverlayProducer.h.

Referenced by onNewRun(), onProcessStart(), and produce().

overlay_filename_

std::string recon::OverlayProducer::overlay_filename_

private

Pileup overlay events input file name.

Definition at line 79 of file OverlayProducer.h.

Referenced by configure(), and onProcessStart().

overlay_incident_id_

int recon::OverlayProducer::overlay_incident_id_ {-1000}

private

For Ecal, overlay hits should be added as contribs.

But these are required to be unique, by the Ecal rconstruction code. So assign a nonsensical trackID, incidentID, and PDG ID to the contribs from overlay. These are hardwired right here.

Definition at line 205 of file OverlayProducer.h.

{-1000};

Referenced by produce().

overlay_passname_

std::string recon::OverlayProducer::overlay_passname_

private

Pileup overlay events input pass name.

Definition at line 111 of file OverlayProducer.h.

Referenced by configure(), and produce().

overlay_pdg_code_

int recon::OverlayProducer::overlay_pdg_code_ {0}

private

Definition at line 207 of file OverlayProducer.h.

{0};

overlay_track_id_

int recon::OverlayProducer::overlay_track_id_ {-1000}

private

Definition at line 206 of file OverlayProducer.h.

{-1000};

params_

framework::config::Parameters recon::OverlayProducer::params_

private

The parameters used to configure this producer.

Definition at line 74 of file OverlayProducer.h.

Referenced by configure(), and onProcessStart().

poisson_mu_

double recon::OverlayProducer::poisson_mu_ {0.}

private

(average) total number of events

Definition at line 139 of file OverlayProducer.h.

{0.};

Referenced by configure(), and produce().

rndm_

std::unique_ptr<TRandom2> recon::OverlayProducer::rndm_

private

Random number generator for number of overlaid events.

TRandom2 slightly (~10%) faster than TRandom3; shorter period but our input files will have way shorter period anyway.

Definition at line 146 of file OverlayProducer.h.

Referenced by onNewRun(), and produce().

rndm_time_

std::unique_ptr<TRandom2> recon::OverlayProducer::rndm_time_

private

Random number generator for pileup event time offset.

TRandom2 slightly (~10%) faster than TRandom3; shorter period but our input files will have way shorter period anyway.

Definition at line 153 of file OverlayProducer.h.

Referenced by onNewRun(), and produce().

sim_passname_

std::string recon::OverlayProducer::sim_passname_

private

To use for finding the sim event bus passengers, mostly a disambiguation.

Definition at line 116 of file OverlayProducer.h.

Referenced by configure(), and produce().

start_event_max_

int recon::OverlayProducer::start_event_max_ {10000}

private

Maximum event number to start overlaying from.

Inclusive.

Definition at line 197 of file OverlayProducer.h.

{10000};

Referenced by configure(), and onNewRun().

start_event_min_

int recon::OverlayProducer::start_event_min_ {1}

private

Minimum event number to start overlaying from.

Inclusive.

Definition at line 191 of file OverlayProducer.h.

{1};

Referenced by configure(), and onNewRun().

time_mean_

double recon::OverlayProducer::time_mean_ {0.}

private

Average position in time (in [ns]) of pileup bunches, relative to the sim event.

Should realistically be 0. Using a non-zero mean and sigma = 0 is however useful for validation.

Definition at line 166 of file OverlayProducer.h.

{0.};

Referenced by configure(), and produce().

time_sigma_

double recon::OverlayProducer::time_sigma_ {0.}

private

Width of pileup bunch spread in time (in [ns]), specified as a sigma of a Gaussian distribution.

Definition at line 159 of file OverlayProducer.h.

{0.};

Referenced by configure(), and produce().

tracker_collections_

std::vector<std::string> recon::OverlayProducer::tracker_collections_

private

List of SimTrackerHit collection(s) to loop over and add hits from, combining sim and pileup.

Definition at line 101 of file OverlayProducer.h.

Referenced by configure(), and produce().

---

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

• Recon/include/Recon/OverlayProducer.h
• Recon/src/Recon/OverlayProducer.cxx