Public Types
typedef ap_ufixed< 16, 14 >	e_t

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

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

virtual void	produce (framework::Event &event)
	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	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	in_proc_

std::string	quad_coll_name_

std::string	combined_quad_coll_name_

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 24 of file TrigHcalEnergySum.h.

Member Typedef Documentation

◆ e_t

ap_ufixed<16, 14> trigger::TrigHcalEnergySum::e_t

Definition at line 33 of file TrigHcalEnergySum.h.

Constructor & Destructor Documentation

◆ TrigHcalEnergySum()

trigger::TrigHcalEnergySum::TrigHcalEnergySum	(	const std::string &	name,
		framework::Process &	process )

inline

Definition at line 26 of file TrigHcalEnergySum.h.

27 : 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 trigger::TrigHcalEnergySum::configure ( framework::config::Parameters & parameters )

virtual

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 11 of file TrigHcalEnergySum.cxx.

                                                               {
  in_proc_ = ps.get<std::string>("inputProc");
  quad_coll_name_ = ps.get<std::string>("quadCollName");
  combined_quad_coll_name_ = ps.get<std::string>("combinedQuadCollName");
}

References framework::config::Parameters::get().

◆ produce()

void trigger::TrigHcalEnergySum::produce ( framework::Event & event )

virtual

Process the event and put new data products into it.

Parameters

event The Event to process.

Implements framework::Producer.

Definition at line 16 of file TrigHcalEnergySum.cxx.

                                                     {
  // mV/ADC: 1.2
  // MeV/MIP: 4.66
  // PE/MIP: 68 (summed over BOTH ends, based on 1808.05219, p38)
  // mV/PE: 5
  // mV/MeV: 72.961 (= 5*68/4.66)
  //  const float mV_per_adc = 1.2;
  // adc gain
  // these are unused, should they be? FIXME
  // const float pe_per_adc = mV_per_adc / 5;
  // const float MeV_per_adc = mV_per_adc / 72.961;
  // const float samp_frac = 371/4e3; // ad-hoc, from a 4 GeV neutron sample
 
  // interaction length in Fe ('steel') = 16.77 cm (132.1 g/cm2)
  // polystyrene = 77.07 cm (81.7 g/cm2)
  // back hcal is 20mm bar, 25mm absorber
  // these are unused, should they be? FIXME
  // const float had_samp_frac =
  //    (20 / 77.07) / (20 / 77.07 + 25 / 16.77);  // 0.148266
  // const float em_samp_frac =
  //    (20 / 41.31) / (20 / 41.31 + 25 / 1.757);                    // 0.032906
  // const float samp_frac = (em_samp_frac + 2 * had_samp_frac) / 3;  //
  // 0.109813 const float attenuation = exp(-1 / 5.);  // 5m attenuation length,
  // 1m half-bar
 
  // for(auto t : event.searchProducts("","","")) std::cout << t.name() << " "
  // << t.passname() << " " << t.type() << std::endl;
 
  if (!event.exists(quad_coll_name_, in_proc_)) {
    // std::cout << "missing collection! : " << quad_coll_name_ << " " <<
    // in_proc_
    // << std::endl;
    return;
  }
 
  // auto
  // oneEndedQuads{event.getObject<ldmx::CaloTrigPrimCollection>(quad_coll_name_)};
  const std::vector<ldmx::CaloTrigPrim> one_ended_quads =
      event.getCollection<ldmx::CaloTrigPrim>(quad_coll_name_, in_proc_);
 
  //
  // sum bar ends to produce the combined quads
  std::map<int, ldmx::CaloTrigPrim> two_ended_quad_map;
  for (const auto& one_ended_quad : one_ended_quads) {
    const ldmx::HcalTriggerID end_id(one_ended_quad.getId());
    ldmx::HcalTriggerID combo_id(end_id.section(), end_id.layer(),
                                 end_id.superstrip(), 2);
    auto ptr = two_ended_quad_map.find(combo_id.raw());
    if (ptr == two_ended_quad_map.end()) {
      two_ended_quad_map[combo_id.raw()] = one_ended_quad;
    } else {
      ptr->second.setPrimitive(ptr->second.getPrimitive() +
                               one_ended_quad.getPrimitive());
    }
  }
  ldmx::CaloTrigPrimCollection two_ended_quads;
  for (auto p : two_ended_quad_map) two_ended_quads.push_back(p.second);
  event.add(combined_quad_coll_name_, two_ended_quads);
 
  //
  // Produce the layer-by-layer energy sums
  const unsigned int layer_max = 50;
  const unsigned int side_layer_max = 16;
  trigger::TrigEnergySumCollection back_layer_sums;
  trigger::TrigEnergySumCollection side_layer_sums;
  back_layer_sums.resize(layer_max);
  for (int i = 0; i < layer_max; i++) back_layer_sums[i].setLayer(i);
  side_layer_sums.resize(side_layer_max);
  for (int i = 0; i < side_layer_max; i++) side_layer_sums[i].setLayer(i);
 
  int total_adc = 0;
  std::map<int, int> section_sum;
  for (auto p : two_ended_quad_map) {
    auto tp = p.second;
    int adc = tp.getPrimitive();
    total_adc += adc;
    ldmx::HcalTriggerID combo_id(tp.getId());
    int ilayer = combo_id.layer();
    if (ilayer >= back_layer_sums.size()) {
      std::cout << "[TrigHcalEnergySum.cxx] Warning(!), layer " << ilayer
                << " is out-of-bounds.\n";
      continue;
    }
    int isec = combo_id.section();
 
    if (isec == 0)
      back_layer_sums[ilayer].setHwEnergy(adc +
                                          back_layer_sums[ilayer].hwEnergy());
    else
      side_layer_sums[ilayer].setHwEnergy(adc +
                                          side_layer_sums[ilayer].hwEnergy());
 
    auto ptr = section_sum.find(isec);
    if (ptr == section_sum.end()) {
      section_sum[isec] = adc;
    } else {
      section_sum[isec] += adc;
    }
  }
  event.add(combined_quad_coll_name_ + "BackLayerSums", back_layer_sums);
  event.add(combined_quad_coll_name_ + "SideLayerSums", side_layer_sums);
 
  trigger::TrigEnergySumCollection section_sums;
  for (auto p : section_sum) {
    section_sums.emplace_back(-1, p.first, p.second);
  }
  event.add(combined_quad_coll_name_ + "SectionSums", section_sums);
 
  // Also store total energy for now
  trigger::TrigEnergySum total_sum;
  total_sum.setLayer(-1);
  total_sum.setHwEnergy(total_adc);
  event.add(combined_quad_coll_name_ + "Sum", total_sum);
}

References framework::Event::exists(), ldmx::HcalTriggerID::layer(), ldmx::DetectorID::raw(), and ldmx::HcalTriggerID::superstrip().

Member Data Documentation

◆ combined_quad_coll_name_

std::string trigger::TrigHcalEnergySum::combined_quad_coll_name_

private

Definition at line 39 of file TrigHcalEnergySum.h.

◆ in_proc_

std::string trigger::TrigHcalEnergySum::in_proc_

private

Definition at line 37 of file TrigHcalEnergySum.h.

◆ quad_coll_name_

std::string trigger::TrigHcalEnergySum::quad_coll_name_

private

Definition at line 38 of file TrigHcalEnergySum.h.

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

Trigger/Algo/include/Trigger/TrigHcalEnergySum.h
Trigger/Algo/src/Trigger/TrigHcalEnergySum.cxx

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ e_t

Constructor & Destructor Documentation

◆ TrigHcalEnergySum()

Member Function Documentation

◆ configure()

◆ produce()

Member Data Documentation

◆ combined_quad_coll_name_

◆ in_proc_

◆ quad_coll_name_