Inheritance diagram for ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer:

Public Member Functions
	HGCalSoARecHitsProducer (edm::ParameterSet const &config)

void	produce (device::Event &iEvent, device::EventSetup const &iSetup) override

	~HGCalSoARecHitsProducer () override=default

Public Member Functions inherited from ALPAKA_ACCELERATOR_NAMESPACE::stream::EDProducer<>
void	produce (edm::Event &iEvent, edm::EventSetup const &iSetup) final

Public Member Functions inherited from edm::ProducerBase< edm::stream::EDProducer, Args... >
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	computeThreshold ()

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_

std::vector< double >	dEdXweights_

int	deltasi_index_regemfac_

std::string	detector_

device::EDPutToken< HGCalSoARecHitsDeviceCollection > const	deviceToken_

double	ecut_

double	fcPerEle_

std::vector< double >	fcPerMip_

edm::EDGetTokenT< HGCRecHitCollection >	hits_token_

bool	initialized_

bool	isNose_

unsigned int	maxlayer_

unsigned	maxNumberOfThickIndices_

std::vector< double >	nonAgedNoises_

hgcal::RecHitTools	rhtools_

double	sciThicknessCorrection_

std::vector< double >	thicknessCorrection_

std::vector< std::vector< double > >	thresholds_

std::vector< std::vector< double > >	v_sigmaNoise_

Additional Inherited Members
Public Types inherited from edm::ProducerBase< edm::stream::EDProducer, Args... >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Protected Member Functions inherited from edm::ProducerBase< edm::stream::EDProducer, Args... >
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetterT< ProductType >	produces (std::string instanceName)

BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetterT< ProductType >	produces (std::string instanceName)

BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetterT< ProductType >	produces (std::string instanceName)

BranchAliasSetterT< ProductType >	produces ()

auto	produces () noexcept

ProducesCollector	producesCollector ()

Detailed Description

Definition at line 22 of file HGCalSoARecHitsProducer.cc.

Constructor & Destructor Documentation

◆ HGCalSoARecHitsProducer()

ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer::HGCalSoARecHitsProducer ( edm::ParameterSet const & config )

inline

Definition at line 24 of file HGCalSoARecHitsProducer.cc.

References edm::ProducerBase< edm::stream::EDProducer, Args... >::produces().

         : detector_(config.getParameter<std::string>("detector")),
           initialized_(false),
           isNose_(detector_ == "HFNose"),
           maxNumberOfThickIndices_(config.getParameter<unsigned>("maxNumberOfThickIndices")),
           fcPerEle_(config.getParameter<double>("fcPerEle")),
           ecut_(config.getParameter<double>("ecut")),
           fcPerMip_(config.getParameter<std::vector<double>>("fcPerMip")),
           nonAgedNoises_(config.getParameter<std::vector<double>>("noises")),
           dEdXweights_(config.getParameter<std::vector<double>>("dEdXweights")),
           thicknessCorrection_(config.getParameter<std::vector<double>>("thicknessCorrection")),
           caloGeomToken_(consumesCollector().esConsumes<CaloGeometry, CaloGeometryRecord>()),
           hits_token_(consumes<HGCRecHitCollection>(config.getParameter<edm::InputTag>("recHits"))),
           deviceToken_{produces()} {}

◆ ~HGCalSoARecHitsProducer()

ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer::~HGCalSoARecHitsProducer ( )

overridedefault

Member Function Documentation

◆ computeThreshold()

void ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer::computeThreshold ( )

inlineprivate

Definition at line 176 of file HGCalSoARecHitsProducer.cc.

References gather_cfg::cout, dEdXweights_, ecut_, fcPerEle_, fcPerMip_, initialized_, isNose_, maxlayer_, maxNumberOfThickIndices_, nonAgedNoises_, thicknessCorrection_, thresholds_, and v_sigmaNoise_.

Referenced by produce().

                             {
       // To support the TDR geometry and also the post-TDR one (v9 onwards), we
       // need to change the logic of the vectors containing signal to noise and
       // thresholds. The first 3 indices will keep on addressing the different
       // thicknesses of the Silicon detectors in CE_E , the next 3 indices will address
       // the thicknesses of the Silicon detectors in CE_H, while the last one, number 6 (the
       // seventh) will address the Scintillators. This change will support both
       // geometries at the same time.
 
       if (initialized_)
         return;  // only need to calculate thresholds once
 
       initialized_ = true;
 
       std::vector<double> dummy;
 
       dummy.resize(maxNumberOfThickIndices_ + !isNose_, 0);  // +1 to accomodate for the Scintillators
       thresholds_.resize(maxlayer_, dummy);
       v_sigmaNoise_.resize(maxlayer_, dummy);
 
       for (unsigned ilayer = 1; ilayer <= maxlayer_; ++ilayer) {
         for (unsigned ithick = 0; ithick < maxNumberOfThickIndices_; ++ithick) {
           float sigmaNoise = 0.001f * fcPerEle_ * nonAgedNoises_[ithick] * dEdXweights_[ilayer] /
                              (fcPerMip_[ithick] * thicknessCorrection_[ithick]);
           thresholds_[ilayer - 1][ithick] = sigmaNoise * ecut_;
           v_sigmaNoise_[ilayer - 1][ithick] = sigmaNoise;
 #if 0
             std::cout << "ilayer: " << ilayer << " nonAgedNoises: " << nonAgedNoises_[ithick]
               << " fcPerEle: " << fcPerEle_ << " fcPerMip: " << fcPerMip_[ithick]
               << " noiseMip: " << fcPerEle_ * nonAgedNoises_[ithick] / fcPerMip_[ithick]
               << " sigmaNoise: " << sigmaNoise << "\n";
 #endif
         }
       }
     }

◆ fillDescriptions()

static void ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

inlinestatic

Definition at line 140 of file HGCalSoARecHitsProducer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                              {
       edm::ParameterSetDescription desc;
       desc.add<std::string>("detector", "EE")->setComment("options EE, FH, BH,  HFNose; other value defaults to EE");
       desc.add<edm::InputTag>("recHits", edm::InputTag("HGCalRecHit", "HGCEERecHits"));
       desc.add<unsigned int>("maxNumberOfThickIndices", 6);
       desc.add<double>("fcPerEle", 0.00016020506);
       desc.add<std::vector<double>>("fcPerMip");
       desc.add<std::vector<double>>("thicknessCorrection");
       desc.add<std::vector<double>>("noises");
       desc.add<std::vector<double>>("dEdXweights");
       desc.add<double>("ecut", 3.);
       descriptions.addWithDefaultLabel(desc);
     }

◆ produce()

void ALPAKA_ACCELERATOR_NAMESPACE::HGCalSoARecHitsProducer::produce	(	device::Event &	iEvent,
		device::EventSetup const &	iSetup
	)

inlineoverridevirtual

Implements ALPAKA_ACCELERATOR_NAMESPACE::stream::EDProducer<>.

Definition at line 41 of file HGCalSoARecHitsProducer.cc.

References caloGeomToken_, hgcalTBTopologyTester_cfi::cells, computeThreshold(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), gather_cfg::cout, deltasi_index_regemfac_, DetId::det(), detector_, CaloRecHit::detid(), deviceToken_, CaloRecHit::energy(), relativeConstraints::geom, ALPAKA_ACCELERATOR_NAMESPACE::device::EventSetup::getHandle(), hgcal::RecHitTools::getLayerWithOffset(), hgcal::RecHitTools::getPosition(), hgcal::RecHitTools::getSiThickIndex(), DetId::HGCalHSi, HGCHEF, hfClusterShapes_cfi::hits, hits_token_, mps_fire::i, iEvent, isNose_, hgcal::RecHitTools::lastLayer(), maxlayer_, maxNumberOfThickIndices_, eostools::move(), hltrates_dqm_sourceclient-live_cfg::offset, position, edm::Handle< T >::product(), DetId::rawId(), rhtools_, hgcal::RecHitTools::setGeometry(), DetId::subdetId(), thresholds_, CaloRecHit::time(), HGCRecHit::timeError(), v_sigmaNoise_, and hgcal::RecHitTools::zside().

                                                                                {
       edm::Handle<HGCRecHitCollection> hits_h;
 
       edm::ESHandle<CaloGeometry> geom = iSetup.getHandle(caloGeomToken_);
       rhtools_.setGeometry(*geom);
       maxlayer_ = rhtools_.lastLayer(isNose_);
 
       hits_h = iEvent.getHandle(hits_token_);
       auto const& hits = *(hits_h.product());
       computeThreshold();
 
       // Count effective hits above threshold
       uint32_t index = 0;
       for (unsigned int i = 0; i < hits.size(); ++i) {
         const HGCRecHit& hgrh = hits[i];
         DetId detid = hgrh.detid();
         unsigned int layerOnSide = (rhtools_.getLayerWithOffset(detid) - 1);
 
         // set sigmaNoise default value 1 to use kappa value directly in case of
         // sensor-independent thresholds
         int thickness_index = rhtools_.getSiThickIndex(detid);
         if (thickness_index == -1) {
           thickness_index = maxNumberOfThickIndices_;
         }
         double storedThreshold = thresholds_[layerOnSide][thickness_index];
         if (hgrh.energy() < storedThreshold)
           continue;  // this sets the ZS threshold at ecut times the sigma noise
         index++;
       }
 
       // Allocate Host SoA will contain one entry for each RecHit above threshold
       HGCalSoARecHitsHostCollection cells(index, iEvent.queue());
       auto cellsView = cells.view();
 
       // loop over all hits and create the Hexel structure, skip energies below ecut
       // for each layer and wafer calculate the thresholds (sigmaNoise and energy)
       // once
       index = 0;
       for (unsigned int i = 0; i < hits.size(); ++i) {
         const HGCRecHit& hgrh = hits[i];
         DetId detid = hgrh.detid();
         unsigned int layerOnSide = (rhtools_.getLayerWithOffset(detid) - 1);
 
         // set sigmaNoise default value 1 to use kappa value directly in case of
         // sensor-independent thresholds
         float sigmaNoise = 1.f;
         int thickness_index = rhtools_.getSiThickIndex(detid);
         if (thickness_index == -1) {
           thickness_index = maxNumberOfThickIndices_;
         }
         double storedThreshold = thresholds_[layerOnSide][thickness_index];
         if (detid.det() == DetId::HGCalHSi || detid.subdetId() == HGCHEF) {
           storedThreshold = thresholds_.at(layerOnSide).at(thickness_index + deltasi_index_regemfac_);
         }
         sigmaNoise = v_sigmaNoise_.at(layerOnSide).at(thickness_index);
 
         if (hgrh.energy() < storedThreshold)
           continue;  // this sets the ZS threshold at ecut times the sigma noise
         // for the sensor
 
         const GlobalPoint position(rhtools_.getPosition(detid));
         int offset = ((rhtools_.zside(detid) + 1) >> 1) * maxlayer_;
         int layer = layerOnSide + offset;
         auto entryInSoA = cellsView[index];
         if (detector_ == "BH") {
           entryInSoA.dim1() = position.eta();
           entryInSoA.dim2() = position.phi();
         }  // else, isSilicon == true and eta phi values will not be used
         else {
           entryInSoA.dim1() = position.x();
           entryInSoA.dim2() = position.y();
         }
         entryInSoA.dim3() = position.z();
         entryInSoA.weight() = hgrh.energy();
         entryInSoA.sigmaNoise() = sigmaNoise;
         entryInSoA.layer() = layer;
         entryInSoA.recHitIndex() = i;
         entryInSoA.detid() = detid.rawId();
         entryInSoA.time() = hgrh.time();
         entryInSoA.timeError() = hgrh.timeError();
         index++;
       }
 #if 0
         std::cout << "Size: " << cells->metadata().size() << " count cells: " << index
           << " i.e. " << cells->metadata().size() << std::endl;
 #endif
 
       if constexpr (!std::is_same_v<Device, alpaka_common::DevHost>) {
         // Trigger copy async to GPU
         //std::cout << "GPU" << std::endl;
         HGCalSoARecHitsDeviceCollection deviceProduct{cells->metadata().size(), iEvent.queue()};
         alpaka::memcpy(iEvent.queue(), deviceProduct.buffer(), cells.const_buffer());
         iEvent.emplace(deviceToken_, std::move(deviceProduct));
       } else {
         //std::cout << "CPU" << std::endl;
         iEvent.emplace(deviceToken_, std::move(cells));
       }
     }