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

virtual	~EcalPreselectionSkimmer ()=default
	Destructor.

virtual void	configure (framework::config::Parameters &) final
	Callback for the EventProcessor to configure itself from the given set of parameters.

virtual void	produce (framework::Event &event) final
	Run the processor and select events that pass pre-selection in ECAL.

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	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 Attributes
std::string	ecal_veto_name_
	Collection Name for veto object.

std::string	ecal_veto_pass_
	Pass Name for veto object.

std::string	ecal_mip_name_
	Collection Name for mip result object.

std::string	ecal_mip_pass_
	Pass Name for mip result object.

double	summed_det_max_
	Max value for summed det.

double	summed_tight_iso_max_
	Max value for summed tigh iso.

double	ecal_back_energy_max_
	Max value for ecal back energy.

int	n_readout_hits_max_
	Max value for num readout hits_.

double	shower_rms_max_
	Max value for shower rms.

double	shower_y_std_max_
	Max value for shower rms in Y.

double	shower_x_std_max_
	Max value for shower rms in X.

double	max_cell_dep_max_
	Max value for maximal cell deposition.

int	std_layer_hit_max_
	Max value for std layer hits_.

int	n_straight_tracks_max_
	Max value for num straight tracks.

double	bdt_disc_min_
	Min value for the BDT disc variable.

int	fiducial_level_
	Level of interest in fiducial events 0: don't care if it's fiducial or not, 1: keep fiducial events only, 2: keep non-fid events only.

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

Definition at line 19 of file EcalPreselectionSkimmer.h.

Constructor & Destructor Documentation

◆ EcalPreselectionSkimmer()

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

Constructor.

Definition at line 11 of file EcalPreselectionSkimmer.cxx.

13 : framework::Producer(name, process) {}

framework::Producer

Base class for a module which produces a data product.

Definition EventProcessor.h:257

framework::Producer::process

virtual void process(Event &event) final

Processing an event for a Producer is calling produce.

Definition EventProcessor.h:278

Member Function Documentation

◆ configure()

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

finalvirtual

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 15 of file EcalPreselectionSkimmer.cxx.

                                                                     {
  ecal_veto_name_ = ps.get<std::string>("ecal_veto_name");
  ecal_veto_pass_ = ps.get<std::string>("ecal_veto_pass");
  ecal_mip_name_ = ps.get<std::string>("ecal_mip_name");
  ecal_mip_pass_ = ps.get<std::string>("ecal_mip_pass");
  summed_det_max_ = ps.get<double>("summed_det_max");              // MeV
  summed_tight_iso_max_ = ps.get<double>("summed_tight_iso_max");  // MeV
  ecal_back_energy_max_ = ps.get<double>("ecal_back_energy_max");  // MeV
  n_readout_hits_max_ = ps.get<int>("n_readout_hits_max");
  shower_rms_max_ = ps.get<double>("shower_rms_max");
  shower_y_std_max_ = ps.get<double>("shower_y_std_max");
  shower_x_std_max_ = ps.get<double>("shower_x_std_max");
  max_cell_dep_max_ = ps.get<double>("max_cell_dep_max");  // MeV
  std_layer_hit_max_ = ps.get<int>("std_layer_hit_max");
  n_straight_tracks_max_ = ps.get<int>("n_straight_tracks_max");
  bdt_disc_min_ = ps.get<double>("bdt_disc_min");
  fiducial_level_ = ps.get<int>("fiducial_level");
 
  return;
}

References bdt_disc_min_, ecal_back_energy_max_, ecal_mip_name_, ecal_mip_pass_, ecal_veto_name_, ecal_veto_pass_, fiducial_level_, framework::config::Parameters::get(), max_cell_dep_max_, n_readout_hits_max_, n_straight_tracks_max_, shower_rms_max_, shower_x_std_max_, shower_y_std_max_, std_layer_hit_max_, summed_det_max_, and summed_tight_iso_max_.

◆ produce()

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

finalvirtual

Run the processor and select events that pass pre-selection in ECAL.

Parameters

event The event to process.

Implements framework::Producer.

Definition at line 36 of file EcalPreselectionSkimmer.cxx.

                                                           {
  bool passed_preselection{false};
  bool fiducial_decision{true};
  const auto &ecal_veto{
      event.getObject<ldmx::EcalVetoResult>(ecal_veto_name_, ecal_veto_pass_)};
  const auto &mip_result{
      event.getObject<ldmx::EcalMipResult>(ecal_mip_name_, ecal_mip_pass_)};
  // Boolean to if we skim for fiducial / nonfiducial
  fiducial_decision = (fiducial_level_ == 0 ||
                       (fiducial_level_ == 1 && ecal_veto.getFiducial()) ||
                       (fiducial_level_ == 2 && !ecal_veto.getFiducial()));
 
  // Boolean to check if we pass preselection
  passed_preselection =
      (ecal_veto.getSummedDet() < summed_det_max_) &&
      (ecal_veto.getSummedTightIso() < summed_tight_iso_max_) &&
      (ecal_veto.getEcalBackEnergy() < ecal_back_energy_max_) &&
      (ecal_veto.getNReadoutHits() < n_readout_hits_max_) &&
      (ecal_veto.getShowerRMS() < shower_rms_max_) &&
      (ecal_veto.getYStd() < shower_y_std_max_) &&
      (ecal_veto.getXStd() < shower_x_std_max_) &&
      (ecal_veto.getMaxCellDep() < max_cell_dep_max_) &&
      (ecal_veto.getStdLayerHit() < std_layer_hit_max_) &&
      (mip_result.getNStraightTracks() < n_straight_tracks_max_) &&
      (ecal_veto.getDisc() > bdt_disc_min_) && fiducial_decision;
 
  // Tell the skimmer to keep or drop the event based on whether preselection
  // passed
  if (passed_preselection) {
    ldmx_log(debug) << "This event passed preselection!";
    setStorageHint(framework::HINT_SHOULD_KEEP);
  } else {
    setStorageHint(framework::HINT_SHOULD_DROP);
  }
  // Add the boolean to the event
  event.add("EcalPreselectionDecision", passed_preselection);
}

References bdt_disc_min_, ecal_back_energy_max_, ecal_mip_name_, ecal_mip_pass_, ecal_veto_name_, ecal_veto_pass_, fiducial_level_, framework::HINT_SHOULD_DROP, framework::HINT_SHOULD_KEEP, max_cell_dep_max_, n_readout_hits_max_, n_straight_tracks_max_, framework::EventProcessor::setStorageHint(), shower_rms_max_, shower_x_std_max_, shower_y_std_max_, std_layer_hit_max_, summed_det_max_, and summed_tight_iso_max_.

Member Data Documentation

◆ bdt_disc_min_

double recon::EcalPreselectionSkimmer::bdt_disc_min_

private

Min value for the BDT disc variable.

Definition at line 67 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ ecal_back_energy_max_

double recon::EcalPreselectionSkimmer::ecal_back_energy_max_

private

Max value for ecal back energy.

Definition at line 51 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ ecal_mip_name_

std::string recon::EcalPreselectionSkimmer::ecal_mip_name_

private

Collection Name for mip result object.

Definition at line 43 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ ecal_mip_pass_

std::string recon::EcalPreselectionSkimmer::ecal_mip_pass_

private

Pass Name for mip result object.

Definition at line 45 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ ecal_veto_name_

std::string recon::EcalPreselectionSkimmer::ecal_veto_name_

private

Collection Name for veto object.

Definition at line 39 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ ecal_veto_pass_

std::string recon::EcalPreselectionSkimmer::ecal_veto_pass_

private

Pass Name for veto object.

Definition at line 41 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ fiducial_level_

int recon::EcalPreselectionSkimmer::fiducial_level_

private

Level of interest in fiducial events 0: don't care if it's fiducial or not, 1: keep fiducial events only, 2: keep non-fid events only.

Definition at line 74 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ max_cell_dep_max_

double recon::EcalPreselectionSkimmer::max_cell_dep_max_

private

Max value for maximal cell deposition.

Definition at line 61 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ n_readout_hits_max_

int recon::EcalPreselectionSkimmer::n_readout_hits_max_

private

Max value for num readout hits_.

Definition at line 53 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ n_straight_tracks_max_

int recon::EcalPreselectionSkimmer::n_straight_tracks_max_

private

Max value for num straight tracks.

Definition at line 65 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ shower_rms_max_

double recon::EcalPreselectionSkimmer::shower_rms_max_

private

Max value for shower rms.

Definition at line 55 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ shower_x_std_max_

double recon::EcalPreselectionSkimmer::shower_x_std_max_

private

Max value for shower rms in X.

Definition at line 59 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ shower_y_std_max_

double recon::EcalPreselectionSkimmer::shower_y_std_max_

private

Max value for shower rms in Y.

Definition at line 57 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ std_layer_hit_max_

int recon::EcalPreselectionSkimmer::std_layer_hit_max_

private

Max value for std layer hits_.

Definition at line 63 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ summed_det_max_

double recon::EcalPreselectionSkimmer::summed_det_max_

private

Max value for summed det.

Definition at line 47 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

◆ summed_tight_iso_max_

double recon::EcalPreselectionSkimmer::summed_tight_iso_max_

private

Max value for summed tigh iso.

Definition at line 49 of file EcalPreselectionSkimmer.h.

Referenced by configure(), and produce().

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

Recon/include/Recon/Skims/EcalPreselectionSkimmer.h
Recon/src/Recon/Skims/EcalPreselectionSkimmer.cxx

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EcalPreselectionSkimmer()

Member Function Documentation

◆ configure()

◆ produce()

Member Data Documentation

◆ bdt_disc_min_

◆ ecal_back_energy_max_

◆ ecal_mip_name_

◆ ecal_mip_pass_

◆ ecal_veto_name_

◆ ecal_veto_pass_

◆ fiducial_level_

◆ max_cell_dep_max_

◆ n_readout_hits_max_

◆ n_straight_tracks_max_

◆ shower_rms_max_

◆ shower_x_std_max_

◆ shower_y_std_max_

◆ std_layer_hit_max_

◆ summed_det_max_

◆ summed_tight_iso_max_