recon::TriggerProcessor Class Reference

Provides a trigger decision for recon using a TriggerResult object. More...

#include <TriggerProcessor.h>

Public Member Functions
	TriggerProcessor (const std::string &name, framework::Process &process)
	Class constructor.
virtual	~TriggerProcessor ()
	Class destructor.
void	configure (framework::config::Parameters &parameters) override
	Configure the processor using the given user specified parameters.
void	produce (framework::Event &event) override
	Run the trigger algorithm and create a TriggerResult object to contain info about the trigger decision such as pass/fail, number of saved variables, etc.
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.

Private Attributes
std::vector< double >	layerESumCuts_
	The energy sum to make cut on.
double	beamEnergy_
	The Beam energy [MeV].
int	mode_ {0}
	The trigger mode to run in.
int	startLayer_ {0}
	The first layer of layer sum.
int	endLayer_ {0}
	The endpoint layer of layer sum.
TString	algoName_
	The name of the trigger algorithm used.
std::string	inputColl_
	The name of the input collection (the Ecal hits).
std::string	inputPass_
	The pass name of the input (the Ecal hits).
std::string	outputColl_
	The name of the output collection (the trigger decision).

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

Provides a trigger decision for recon using a TriggerResult object.

Note

TriggerProcessor takes in a set of parameters to be used in defining the trigger algorithm. An event is passed to the processor and the relevant algorithms are then run on the event (ECAL layer sum). A trigger decision is executed and the decision along with the algorithm name and relevant variables are stored in a TriggerResult object which is added to the collection.

Definition at line 32 of file TriggerProcessor.h.

Constructor & Destructor Documentation

TriggerProcessor()

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

inline

Class constructor.

Definition at line 37 of file TriggerProcessor.h.

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

~TriggerProcessor()

virtual recon::TriggerProcessor::~TriggerProcessor ( )

inlinevirtual

Class destructor.

Definition at line 43 of file TriggerProcessor.h.

{}

Member Function Documentation

configure()

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

overridevirtual

Configure the processor using the given user specified parameters.

Parameters

parameters

Set of parameters used to configure this processor.

Reimplemented from framework::EventProcessor.

Definition at line 8 of file TriggerProcessor.cxx.

                                                                      {
  layerESumCuts_ = parameters.getParameter<std::vector<double>>("thresholds");
  beamEnergy_ = parameters.getParameter<double>("beamEnergy");
  mode_ = parameters.getParameter<int>("mode");
  startLayer_ = parameters.getParameter<int>("start_layer");
  endLayer_ = parameters.getParameter<int>("end_layer");
  inputColl_ = parameters.getParameter<std::string>("input_collection");
  inputPass_ = parameters.getParameter<std::string>("input_pass");
  outputColl_ = parameters.getParameter<std::string>("trigger_collection");
 
  if (mode_ == 0) {
    algoName_ = "LayerSumTrig";
  } else if (mode_ == 1) {
    algoName_ = "CenterTower";
  }
}

References algoName_, beamEnergy_, endLayer_, framework::config::Parameters::getParameter(), inputColl_, inputPass_, layerESumCuts_, mode_, outputColl_, and startLayer_.

produce()

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

overridevirtual

Run the trigger algorithm and create a TriggerResult object to contain info about the trigger decision such as pass/fail, number of saved variables, etc.

param event The event to run trigger algorithm on.

Grab the Ecal hit collection for the given event

There are three scenarios:

1. No Incoming Electrons - If the electron count is 0 or negative (undetermined), then we set the sum-energy cut to zero so the event always fails.
2. Num Electrons Listed in Thresholds - Pull cut from list
3. More electrons than listed - Set threshold as 'threshold_for_1e + nExtraElectrons*beamEnergy' Note that the "overflow" formula here is too naive. It should be a fct( nElectrons, 1e_thr, beamE), taking how sigma evolves with multiplicity into account. a simple scaling might suffice there too assuming energy cuts are listed as [ Ecut_1e, Ecut_2e, ... ]

Loop over all ecal hits in the given event

Implements framework::Producer.

Definition at line 25 of file TriggerProcessor.cxx.

                                                    {
  const std::vector<ldmx::EcalHit> ecalRecHits =
      event.getCollection<ldmx::EcalHit>(inputColl_, inputPass_);
 
  // number of electrons in this event
  const int nElectrons{event.getElectronCount()};
 
  double layerESumCut{0.};
  if (nElectrons <= 0)
    layerESumCut = 0.;  // always fail if no electrons
  else if (nElectrons <= layerESumCuts_.size())
    layerESumCut = layerESumCuts_.at(nElectrons - 1);
  else
    layerESumCut = layerESumCuts_.at(0) + (nElectrons - 1) * beamEnergy_;
 
  ldmx_log(debug) << "Got trigger energy cut " << layerESumCut << " for "
                  << nElectrons << " electrons counted in the event.";
 
  std::vector<double> layerDigiE(100, 0.0);  // big empty vector..
 
  for (const ldmx::EcalHit& hit : ecalRecHits) {
    ldmx::EcalID id(hit.getID());
    if (id.layer() < layerDigiE.size()) {  // just to be safe...
      if (mode_ == 0) {  // Sum over all cells in a given layer
        layerDigiE[id.layer()] += hit.getEnergy();
      } else if (mode_ == 1) {  // Sum over cells in central tower only
                                // std::pair<float, float> xyPos =
        // hit->getCellCentroidXYPair(hit->getID()); float cellRadius =
        // sqrt(pow(xyPos.first, 2) + pow(xyPos.second, 2)); if (cellRadius <
        // MAGICNUMBERHERE) {
        //    layerDigiE[hit->getLayer()] += hit->getEnergy();
        //}
      }
    }
  }
 
  float layerSum = 0;
  bool pass = false;
 
  for (int iL = startLayer_; iL < endLayer_; ++iL) {
    layerSum += layerDigiE[iL];
  }
 
  pass = (layerSum <= layerESumCut);
  ldmx_log(debug) << "Got trigger energy sum " << layerSum
                  << "; and decision is pass = " << pass;
 
  ldmx::TriggerResult result;
  result.set(algoName_, pass, 4);
  result.setAlgoVar(0, layerSum);
  result.setAlgoVar(1, layerESumCut);
  result.setAlgoVar(2, endLayer_ - startLayer_);
  result.setAlgoVar(3, nElectrons);
 
  event.add(outputColl_, result);
 
  // mark the event
  if (pass)
    setStorageHint(framework::hint_shouldKeep);
  else
    setStorageHint(framework::hint_shouldDrop);
}

References algoName_, beamEnergy_, endLayer_, framework::hint_shouldDrop, framework::hint_shouldKeep, inputColl_, inputPass_, layerESumCuts_, mode_, outputColl_, ldmx::TriggerResult::set(), ldmx::TriggerResult::setAlgoVar(), framework::EventProcessor::setStorageHint(), and startLayer_.

Member Data Documentation

algoName_

TString recon::TriggerProcessor::algoName_

private

The name of the trigger algorithm used.

Definition at line 87 of file TriggerProcessor.h.

Referenced by configure(), and produce().

beamEnergy_

double recon::TriggerProcessor::beamEnergy_

private

The Beam energy [MeV].

Definition at line 66 of file TriggerProcessor.h.

Referenced by configure(), and produce().

endLayer_

int recon::TriggerProcessor::endLayer_ {0}

private

The endpoint layer of layer sum.

not inclusive - i.e. this is the last layer that is included in the layer sum.

Definition at line 84 of file TriggerProcessor.h.

{0};

Referenced by configure(), and produce().

inputColl_

std::string recon::TriggerProcessor::inputColl_

private

The name of the input collection (the Ecal hits).

Definition at line 90 of file TriggerProcessor.h.

Referenced by configure(), and produce().

inputPass_

std::string recon::TriggerProcessor::inputPass_

private

The pass name of the input (the Ecal hits).

Definition at line 93 of file TriggerProcessor.h.

Referenced by configure(), and produce().

layerESumCuts_

std::vector<double> recon::TriggerProcessor::layerESumCuts_

private

The energy sum to make cut on.

Definition at line 63 of file TriggerProcessor.h.

Referenced by configure(), and produce().

mode_

int recon::TriggerProcessor::mode_ {0}

private

The trigger mode to run in.

Mode zero sums over all cells in layer, while in mode 1 only cells in center module are summed over. (TODO)

Definition at line 73 of file TriggerProcessor.h.

{0};

Referenced by configure(), and produce().

outputColl_

std::string recon::TriggerProcessor::outputColl_

private

The name of the output collection (the trigger decision).

Definition at line 96 of file TriggerProcessor.h.

Referenced by configure(), and produce().

startLayer_

int recon::TriggerProcessor::startLayer_ {0}

private

The first layer of layer sum.

Definition at line 76 of file TriggerProcessor.h.

{0};

Referenced by configure(), and produce().

---

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

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