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	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::string	inProc_

std::string	quadCollName_

std::string	combinedQuadCollName_

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

Member Typedef Documentation

◆ e_t

typedef 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:252

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.

                                                               {
  inProc_ = ps.getParameter<std::string>("inputProc");
  quadCollName_ = ps.getParameter<std::string>("quadCollName");
  combinedQuadCollName_ = ps.getParameter<std::string>("combinedQuadCollName");
}

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

◆ 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(quadCollName_, inProc_)) {
    // std::cout << "missing collection! : " << quadCollName_ << " " << inProc_
    // << std::endl;
    return;
  }
 
  // auto
  // oneEndedQuads{event.getObject<ldmx::CaloTrigPrimCollection>(quadCollName_)};
  const std::vector<ldmx::CaloTrigPrim> oneEndedQuads =
      event.getCollection<ldmx::CaloTrigPrim>(quadCollName_, inProc_);
 
  //
  // sum bar ends to produce the combined quads
  std::map<int, ldmx::CaloTrigPrim> twoEndedQuadMap;
  for (const auto& oneEndedQuad : oneEndedQuads) {
    const ldmx::HcalTriggerID end_id(oneEndedQuad.getId());
    ldmx::HcalTriggerID combo_id(end_id.section(), end_id.layer(),
                                 end_id.superstrip(), 2);
    auto ptr = twoEndedQuadMap.find(combo_id.raw());
    if (ptr == twoEndedQuadMap.end()) {
      twoEndedQuadMap[combo_id.raw()] = oneEndedQuad;
    } else {
      ptr->second.setPrimitive(ptr->second.getPrimitive() +
                               oneEndedQuad.getPrimitive());
    }
  }
  ldmx::CaloTrigPrimCollection twoEndedQuads;
  for (auto p : twoEndedQuadMap) twoEndedQuads.push_back(p.second);
  event.add(combinedQuadCollName_, twoEndedQuads);
 
  //
  // Produce the layer-by-layer energy sums
  const unsigned int LayerMax = 50;
  const unsigned int SideLayerMax = 16;
  trigger::TrigEnergySumCollection backLayerSums;
  trigger::TrigEnergySumCollection sideLayerSums;
  backLayerSums.resize(LayerMax);
  for (int i = 0; i < LayerMax; i++) backLayerSums[i].setLayer(i);
  sideLayerSums.resize(SideLayerMax);
  for (int i = 0; i < SideLayerMax; i++) sideLayerSums[i].setLayer(i);
 
  int total_adc = 0;
  std::map<int, int> section_sum;
  for (auto p : twoEndedQuadMap) {
    auto tp = p.second;
    int adc = tp.getPrimitive();
    total_adc += adc;
    ldmx::HcalTriggerID combo_id(tp.getId());
    int ilayer = combo_id.layer();
    if (ilayer >= backLayerSums.size()) {
      std::cout << "[TrigHcalEnergySum.cxx] Warning(!), layer " << ilayer
                << " is out-of-bounds.\n";
      continue;
    }
    int isec = combo_id.section();
 
    if (isec == 0)
      backLayerSums[ilayer].setHwEnergy(adc + backLayerSums[ilayer].hwEnergy());
    else
      sideLayerSums[ilayer].setHwEnergy(adc + sideLayerSums[ilayer].hwEnergy());
 
    auto ptr = section_sum.find(isec);
    if (ptr == section_sum.end()) {
      section_sum[isec] = adc;
    } else {
      section_sum[isec] += adc;
    }
  }
  event.add(combinedQuadCollName_ + "BackLayerSums", backLayerSums);
  event.add(combinedQuadCollName_ + "SideLayerSums", sideLayerSums);
 
  trigger::TrigEnergySumCollection sectionSums;
  for (auto p : section_sum) {
    sectionSums.emplace_back(-1, p.first, p.second);
  }
  event.add(combinedQuadCollName_ + "SectionSums", sectionSums);
 
  // Also store total energy for now
  trigger::TrigEnergySum totalSum;
  totalSum.setLayer(-1);
  totalSum.setHwEnergy(total_adc);
  event.add(combinedQuadCollName_ + "Sum", totalSum);
}

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

Member Data Documentation

◆ combinedQuadCollName_

std::string trigger::TrigHcalEnergySum::combinedQuadCollName_

private

Definition at line 39 of file TrigHcalEnergySum.h.

◆ inProc_

std::string trigger::TrigHcalEnergySum::inProc_

private

Definition at line 37 of file TrigHcalEnergySum.h.

◆ quadCollName_

std::string trigger::TrigHcalEnergySum::quadCollName_

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

◆ combinedQuadCollName_

◆ inProc_

◆ quadCollName_